static void
grain_stat(GwyContainer *data, GwyRunType run)
{
GtkWidget *dialog, *table, *hbox, *button;
GwyDataField *dfield, *mfield;
GwySIUnit *siunit, *siunit2;
GwySIValueFormat *vf, *vf2;
gint xres, yres, ngrains;
gdouble total_area, area, size, vol_0, vol_min, vol_laplace, bound_len, v;
gdouble *values = NULL;
gint *grains;
GString *str, *str2;
GPtrArray *report;
const guchar *title;
gchar *key, *value;
gint row, id;
guint i, maxlen;
g_return_if_fail(run & STAT_RUN_MODES);
gwy_app_data_browser_get_current(GWY_APP_DATA_FIELD, &dfield,
GWY_APP_MASK_FIELD, &mfield,
GWY_APP_DATA_FIELD_ID, &id,
0);
g_return_if_fail(dfield);
g_return_if_fail(mfield);
report = g_ptr_array_sized_new(20);
if (gwy_container_gis_string_by_name(data, "/filename", &title)) {
g_ptr_array_add(report, _("File:"));
g_ptr_array_add(report, g_strdup(title));
}
key = g_strdup_printf("/%d/data/title", id);
if (gwy_container_gis_string_by_name(data, key, &title)) {
g_ptr_array_add(report, _("Data channel:"));
g_ptr_array_add(report, g_strdup(title));
}
g_free(key);
/* Make empty line in the report */
g_ptr_array_add(report, NULL);
g_ptr_array_add(report, NULL);
xres = gwy_data_field_get_xres(mfield);
yres = gwy_data_field_get_yres(mfield);
total_area = gwy_data_field_get_xreal(dfield)
*gwy_data_field_get_yreal(dfield);
grains = g_new0(gint, xres*yres);
ngrains = gwy_data_field_number_grains(mfield, grains);
area = grains_get_total_value(dfield, ngrains, grains, &values,
GWY_GRAIN_VALUE_PROJECTED_AREA);
size = grains_get_total_value(dfield, ngrains, grains, &values,
GWY_GRAIN_VALUE_EQUIV_SQUARE_SIDE);
vol_0 = grains_get_total_value(dfield, ngrains, grains, &values,
GWY_GRAIN_VALUE_VOLUME_0);
vol_min = grains_get_total_value(dfield, ngrains, grains, &values,
GWY_GRAIN_VALUE_VOLUME_MIN);
vol_laplace = grains_get_total_value(dfield, ngrains, grains, &values,
GWY_GRAIN_VALUE_VOLUME_LAPLACE);
bound_len = grains_get_total_value(dfield, ngrains, grains, &values,
GWY_GRAIN_VALUE_FLAT_BOUNDARY_LENGTH);
g_free(values);
g_free(grains);
dialog = gtk_dialog_new_with_buttons(_("Grain Statistics"), NULL, 0,
GTK_STOCK_CLOSE, GTK_RESPONSE_CLOSE,
NULL);
gwy_help_add_to_proc_dialog(GTK_DIALOG(dialog), GWY_HELP_DEFAULT);
table = gtk_table_new(10, 2, FALSE);
gtk_container_add(GTK_CONTAINER(GTK_DIALOG(dialog)->vbox), table);
gtk_container_set_border_width(GTK_CONTAINER(table), 4);
row = 0;
str = g_string_new(NULL);
str2 = g_string_new(NULL);
g_string_printf(str, "%d", ngrains);
add_report_row(GTK_TABLE(table), &row, _("Number of grains:"),
str->str, str->str, report);
siunit = gwy_data_field_get_si_unit_xy(dfield);
siunit2 = gwy_si_unit_power(siunit, 2, NULL);
v = area;
vf = gwy_si_unit_get_format_with_digits(siunit2,
GWY_SI_UNIT_FORMAT_VFMARKUP,
v, 3, NULL);
vf2 = gwy_si_unit_get_format_with_digits(siunit2, GWY_SI_UNIT_FORMAT_PLAIN,
v, 3, NULL);
g_string_printf(str, "%.*f %s", vf->precision, v/vf->magnitude, vf->units);
g_string_printf(str2, "%.*f %s",
vf2->precision, v/vf2->magnitude, vf2->units);
add_report_row(GTK_TABLE(table), &row, _("Total projected area (abs.):"),
str->str, str2->str, report);
g_string_printf(str, "%.2f %%", 100.0*area/total_area);
add_report_row(GTK_TABLE(table), &row, _("Total projected area (rel.):"),
str->str, str->str, report);
v = area/ngrains;
gwy_si_unit_get_format_with_digits(siunit2, GWY_SI_UNIT_FORMAT_VFMARKUP,
v, 3, vf);
gwy_si_unit_get_format_with_digits(siunit2, GWY_SI_UNIT_FORMAT_PLAIN,
v, 3, vf2);
g_string_printf(str, "%.*f %s", vf->precision, v/vf->magnitude, vf->units);
g_string_printf(str2, "%.*f %s",
vf2->precision, v/vf2->magnitude, vf2->units);
add_report_row(GTK_TABLE(table), &row, _("Mean grain area:"),
str->str, str2->str, report);
v = size/ngrains;
gwy_si_unit_get_format_with_digits(siunit, GWY_SI_UNIT_FORMAT_VFMARKUP,
v, 3, vf);
gwy_si_unit_get_format_with_digits(siunit, GWY_SI_UNIT_FORMAT_PLAIN,
v, 3, vf2);
g_string_printf(str, "%.*f %s", vf->precision, v/vf->magnitude, vf->units);
g_string_printf(str2, "%.*f %s",
vf2->precision, v/vf2->magnitude, vf2->units);
add_report_row(GTK_TABLE(table), &row, _("Mean grain size:"),
str->str, str2->str, report);
siunit = gwy_data_field_get_si_unit_z(dfield);
gwy_si_unit_multiply(siunit2, siunit, siunit2);
v = vol_0;
gwy_si_unit_get_format_with_digits(siunit2, GWY_SI_UNIT_FORMAT_VFMARKUP,
v, 3, vf);
gwy_si_unit_get_format_with_digits(siunit2, GWY_SI_UNIT_FORMAT_PLAIN,
v, 3, vf2);
g_string_printf(str, "%.*f %s", vf->precision, v/vf->magnitude, vf->units);
g_string_printf(str2, "%.*f %s",
vf2->precision, v/vf2->magnitude, vf2->units);
add_report_row(GTK_TABLE(table), &row, _("Total grain volume (zero):"),
str->str, str2->str, report);
v = vol_min;
gwy_si_unit_get_format_with_digits(siunit2, GWY_SI_UNIT_FORMAT_VFMARKUP,
v, 3, vf);
gwy_si_unit_get_format_with_digits(siunit2, GWY_SI_UNIT_FORMAT_PLAIN,
v, 3, vf2);
g_string_printf(str, "%.*f %s", vf->precision, v/vf->magnitude, vf->units);
g_string_printf(str2, "%.*f %s",
vf2->precision, v/vf2->magnitude, vf2->units);
add_report_row(GTK_TABLE(table), &row, _("Total grain volume (minimum):"),
str->str, str2->str, report);
v = vol_laplace;
gwy_si_unit_get_format_with_digits(siunit2, GWY_SI_UNIT_FORMAT_VFMARKUP,
v, 3, vf);
gwy_si_unit_get_format_with_digits(siunit2, GWY_SI_UNIT_FORMAT_PLAIN,
v, 3, vf2);
g_string_printf(str, "%.*f %s", vf->precision, v/vf->magnitude, vf->units);
g_string_printf(str2, "%.*f %s",
vf2->precision, v/vf2->magnitude, vf2->units);
add_report_row(GTK_TABLE(table), &row, _("Total grain volume (laplacian):"),
str->str, str2->str, report);
v = bound_len;
gwy_si_unit_get_format_with_digits(siunit, GWY_SI_UNIT_FORMAT_VFMARKUP,
v, 3, vf);
gwy_si_unit_get_format_with_digits(siunit, GWY_SI_UNIT_FORMAT_PLAIN,
v, 3, vf2);
g_string_printf(str, "%.*f %s", vf->precision, v/vf->magnitude, vf->units);
g_string_printf(str2, "%.*f %s",
vf2->precision, v/vf2->magnitude, vf2->units);
add_report_row(GTK_TABLE(table), &row, _("Total projected boundary length:"),
str->str, str2->str, report);
gwy_si_unit_value_format_free(vf2);
gwy_si_unit_value_format_free(vf);
g_object_unref(siunit2);
maxlen = 0;
for (i = 0; i < report->len/2; i++) {
key = (gchar*)g_ptr_array_index(report, 2*i);
if (key)
maxlen = MAX(strlen(key), maxlen);
}
g_string_truncate(str, 0);
g_string_append(str, _("Grain Statistics"));
g_string_append(str, "\n\n");
for (i = 0; i < report->len/2; i++) {
key = (gchar*)g_ptr_array_index(report, 2*i);
if (key) {
value = (gchar*)g_ptr_array_index(report, 2*i + 1);
g_string_append_printf(str, "%-*s %s\n", maxlen+1, key, value);
g_free(value);
}
else
g_string_append_c(str, '\n');
}
g_ptr_array_free(report, TRUE);
g_object_set_data(G_OBJECT(dialog), "report", str->str);
hbox = gtk_hbox_new(FALSE, 0);
gtk_box_pack_start(GTK_BOX(GTK_DIALOG(dialog)->vbox), hbox,
FALSE, FALSE, 0);
button = grain_stats_add_aux_button(hbox, GTK_STOCK_SAVE,
_("Save table to a file"));
g_signal_connect_swapped(button, "clicked",
G_CALLBACK(grain_stat_save), dialog);
button = grain_stats_add_aux_button(hbox, GTK_STOCK_COPY,
_("Copy table to clipboard"));
g_signal_connect_swapped(button, "clicked",
G_CALLBACK(grain_stat_copy), dialog);
gtk_widget_show_all(dialog);
gtk_dialog_run(GTK_DIALOG(dialog));
gtk_widget_destroy(dialog);
g_string_free(str2, TRUE);
g_string_free(str, TRUE);
}
static void
grain_stat(G_GNUC_UNUSED GwyContainer *data, GwyRunType run)
{
GtkWidget *dialog, *table;
GwyDataField *dfield, *mfield;
GwySIUnit *siunit, *siunit2;
GwySIValueFormat *vf;
gint xres, yres, ngrains;
gdouble total_area, area, size, vol_0, vol_min, vol_laplace, v;
gdouble *values = NULL;
gint *grains;
GString *str;
gint row;
g_return_if_fail(run & STAT_RUN_MODES);
gwy_app_data_browser_get_current(GWY_APP_DATA_FIELD, &dfield,
GWY_APP_MASK_FIELD, &mfield,
0);
g_return_if_fail(dfield);
g_return_if_fail(mfield);
xres = gwy_data_field_get_xres(mfield);
yres = gwy_data_field_get_yres(mfield);
total_area = gwy_data_field_get_xreal(dfield)
*gwy_data_field_get_yreal(dfield);
grains = g_new0(gint, xres*yres);
ngrains = gwy_data_field_number_grains(mfield, grains);
area = grains_get_total_value(dfield, ngrains, grains, &values,
GWY_GRAIN_VALUE_PROJECTED_AREA);
size = grains_get_total_value(dfield, ngrains, grains, &values,
GWY_GRAIN_VALUE_EQUIV_SQUARE_SIDE);
vol_0 = grains_get_total_value(dfield, ngrains, grains, &values,
GWY_GRAIN_VALUE_VOLUME_0);
vol_min = grains_get_total_value(dfield, ngrains, grains, &values,
GWY_GRAIN_VALUE_VOLUME_MIN);
vol_laplace = grains_get_total_value(dfield, ngrains, grains, &values,
GWY_GRAIN_VALUE_VOLUME_LAPLACE);
g_free(values);
g_free(grains);
dialog = gtk_dialog_new_with_buttons(_("Grain Statistics"), NULL, 0,
GTK_STOCK_CLOSE, GTK_RESPONSE_CLOSE,
NULL);
gtk_dialog_set_has_separator(GTK_DIALOG(dialog), FALSE);
table = gtk_table_new(7, 2, FALSE);
gtk_container_add(GTK_CONTAINER(GTK_DIALOG(dialog)->vbox), table);
gtk_container_set_border_width(GTK_CONTAINER(table), 4);
row = 0;
str = g_string_new(NULL);
g_string_printf(str, "%d", ngrains);
add_report_row(GTK_TABLE(table), &row, _("Number of grains:"), str->str);
siunit = gwy_data_field_get_si_unit_xy(dfield);
siunit2 = gwy_si_unit_power(siunit, 2, NULL);
v = area;
vf = gwy_si_unit_get_format(siunit2, GWY_SI_UNIT_FORMAT_VFMARKUP, v, NULL);
g_string_printf(str, "%.*f %s", vf->precision, v/vf->magnitude, vf->units);
add_report_row(GTK_TABLE(table), &row, _("Total projected area (abs.):"),
str->str);
g_string_printf(str, "%.2f %%", 100.0*area/total_area);
add_report_row(GTK_TABLE(table), &row, _("Total projected area (rel.):"),
str->str);
v = area/ngrains;
gwy_si_unit_get_format(siunit2, GWY_SI_UNIT_FORMAT_VFMARKUP, v, vf);
g_string_printf(str, "%.*f %s", vf->precision, v/vf->magnitude, vf->units);
add_report_row(GTK_TABLE(table), &row, _("Mean grain area:"), str->str);
v = size/ngrains;
gwy_si_unit_get_format(siunit, GWY_SI_UNIT_FORMAT_VFMARKUP, v, vf);
g_string_printf(str, "%.*f %s", vf->precision, v/vf->magnitude, vf->units);
add_report_row(GTK_TABLE(table), &row, _("Mean grain size:"), str->str);
siunit = gwy_data_field_get_si_unit_z(dfield);
gwy_si_unit_multiply(siunit2, siunit, siunit2);
v = vol_0;
gwy_si_unit_get_format(siunit2, GWY_SI_UNIT_FORMAT_VFMARKUP, v, vf);
g_string_printf(str, "%.*f %s", vf->precision, v/vf->magnitude, vf->units);
add_report_row(GTK_TABLE(table), &row, _("Total grain volume (zero):"),
str->str);
v = vol_min;
gwy_si_unit_get_format(siunit2, GWY_SI_UNIT_FORMAT_VFMARKUP, v, vf);
g_string_printf(str, "%.*f %s", vf->precision, v/vf->magnitude, vf->units);
add_report_row(GTK_TABLE(table), &row, _("Total grain volume (minimum):"),
str->str);
v = vol_laplace;
gwy_si_unit_get_format(siunit2, GWY_SI_UNIT_FORMAT_VFMARKUP, v, vf);
g_string_printf(str, "%.*f %s", vf->precision, v/vf->magnitude, vf->units);
add_report_row(GTK_TABLE(table), &row, _("Total grain volume (laplacian):"),
str->str);
gwy_si_unit_value_format_free(vf);
g_string_free(str, TRUE);
g_object_unref(siunit2);
gtk_widget_show_all(dialog);
gtk_dialog_run(GTK_DIALOG(dialog));
gtk_widget_destroy(dialog);
}
