static void
dwt(GwyContainer *data, GwyRunType run)
{
GtkWidget *dialog;
GwyDataField *dfield;
GwyDataLine *wtcoefs;
DWTArgs args;
gboolean ok;
gint xsize, ysize, newsize, newid, oldid, i;
g_return_if_fail(run & DWT_RUN_MODES);
gwy_app_data_browser_get_current(GWY_APP_DATA_FIELD, &dfield,
GWY_APP_DATA_FIELD_ID, &oldid,
0);
g_return_if_fail(dfield);
xsize = gwy_data_field_get_xres(dfield);
ysize = gwy_data_field_get_yres(dfield);
if (xsize != ysize) {
dialog = gtk_message_dialog_new
(gwy_app_find_window_for_channel(data, oldid),
GTK_DIALOG_DESTROY_WITH_PARENT,
GTK_MESSAGE_ERROR,
GTK_BUTTONS_OK,
_("%s: Data must be square."), "DWT");
gtk_dialog_run(GTK_DIALOG(dialog));
gtk_widget_destroy(dialog);
return;
}
for (newsize = 1, i = xsize-1; i; i >>= 1, newsize <<= 1)
;
dwt_load_args(gwy_app_settings_get(), &args);
if (run == GWY_RUN_INTERACTIVE) {
ok = dwt_dialog(&args, xsize, newsize);
dwt_save_args(gwy_app_settings_get(), &args);
if (!ok)
return;
}
dfield = gwy_data_field_new_resampled(dfield, newsize, newsize,
args.interp);
gwy_data_field_add(dfield, -gwy_data_field_get_avg(dfield));
wtcoefs = GWY_DATA_LINE(gwy_data_line_new(10, 10, TRUE));
wtcoefs = gwy_dwt_set_coefficients(wtcoefs, args.wavelet);
gwy_data_field_dwt(dfield, wtcoefs, 1, 4);
newid = gwy_app_data_browser_add_data_field(dfield, data, TRUE);
g_object_unref(dfield);
gwy_app_set_data_field_title(data, newid, _("DWT"));
gwy_app_sync_data_items(data, data, oldid, newid, FALSE,
GWY_DATA_ITEM_PALETTE, 0);
g_object_unref(wtcoefs);
}
/**
* gwy_data_line_psdf:
* @data_line: data line
* @target_line: result data line
* @windowing: windowing method
* @interpolation: interpolation method
*
* Copmutes power spectral density function and stores the values in
* @target_line.
**/
void
gwy_data_line_psdf(GwyDataLine *data_line, GwyDataLine *target_line, gint windowing, gint interpolation)
{
GwyDataLine *iin, *rout, *iout;
gint i, order, newres, oldres;
g_return_if_fail(GWY_IS_DATA_LINE(data_line));
oldres = data_line->res;
order = (gint) floor(log ((gdouble)data_line->res)/log (2.0)+0.5);
newres = (gint) pow(2,order);
iin = (GwyDataLine *)gwy_data_line_new(newres, data_line->real, TRUE);
rout = (GwyDataLine *)gwy_data_line_new(newres, data_line->real, TRUE);
iout = (GwyDataLine *)gwy_data_line_new(newres, data_line->real, TRUE);
/*resample to 2^N (this could be done within FFT, but with loss of precision)*/
gwy_data_line_resample(data_line, newres, interpolation);
gwy_data_line_fft(data_line, iin, rout, iout, gwy_data_line_fft_hum,
windowing, 1, interpolation,
TRUE, TRUE);
gwy_data_line_resample(target_line, newres/2.0, GWY_INTERPOLATION_NONE);
/*compute module*/
for (i = 0; i < (newres/2); i++) {
target_line->data[i] = (rout->data[i]*rout->data[i] + iout->data[i]*iout->data[i])
*data_line->real/(newres*newres*2*G_PI);
}
target_line->real = 2*G_PI*target_line->res/data_line->real;
/*resample to given output size*/
gwy_data_line_resample(target_line, oldres/2.0, interpolation);
g_object_unref(rout);
g_object_unref(iin);
g_object_unref(iout);
}
static void
line_correct_median(GwyContainer *data, GwyRunType run)
{
GwyDataField *dfield;
GwyDataLine *line, *modi;
gint xres, yres, i;
GQuark dquark;
gdouble median;
g_return_if_fail(run & LINECORR_RUN_MODES);
gwy_app_data_browser_get_current(GWY_APP_DATA_FIELD, &dfield,
GWY_APP_DATA_FIELD_KEY, &dquark,
0);
g_return_if_fail(dfield && dquark);
gwy_app_undo_qcheckpointv(data, 1, &dquark);
xres = gwy_data_field_get_xres(dfield);
line = gwy_data_line_new(xres, 1.0, FALSE);
yres = gwy_data_field_get_yres(dfield);
modi = gwy_data_line_new(yres, 1.0, FALSE);
for (i = 0; i < yres; i++) {
gwy_data_field_get_row(dfield, line, i);
median = gwy_math_median(xres, gwy_data_line_get_data(line));
gwy_data_line_set_val(modi, i, median);
}
median = gwy_data_line_get_median(modi);
for (i = 0; i < yres; i++) {
gwy_data_field_area_add(dfield, 0, i, xres, 1,
median - gwy_data_line_get_val(modi, i));
}
g_object_unref(modi);
g_object_unref(line);
gwy_data_field_data_changed(dfield);
}
static gboolean
height_dist(GwyContainer *data, GwyRunType run)
{
GString *lab;
GtkWidget *graph;
GwyDataWindow *data_window;
GwyGraphAutoProperties prop;
GwyDataLine *dataline;
GwyDataField *dfield;
GwyDataField *mfield;
GwySIValueFormat *units;
g_return_val_if_fail(run & DIST_RUN_MODES, FALSE);
g_return_val_if_fail(gwy_container_contains_by_name(data, "/0/mask"),
FALSE);
graph = gwy_graph_new();
gwy_graph_get_autoproperties(GWY_GRAPH(graph), &prop);
prop.is_point = 0;
prop.is_line = 1;
gwy_graph_set_autoproperties(GWY_GRAPH(graph), &prop);
dfield = GWY_DATA_FIELD(gwy_container_get_object_by_name(data, "/0/data"));
mfield = GWY_DATA_FIELD(gwy_container_get_object_by_name(data, "/0/mask"));
dataline = (GwyDataLine*)gwy_data_line_new(10, 10, TRUE);
gwy_data_field_grains_get_height_distribution(dfield, mfield, dataline);
lab = g_string_new(_("Grain height histogram"));
units = gwy_si_unit_get_format(dfield->si_unit_xy, dataline->real, NULL);
gwy_graph_add_dataline_with_units(GWY_GRAPH(graph), dataline, 0, lab, NULL,
units->magnitude, 1, units->units, " ");
data_window = gwy_app_data_window_get_for_data(data);
gwy_app_graph_window_create_for_window(GWY_GRAPH(graph), data_window,
_("Grain height distribution"));
g_string_free(lab, TRUE);
g_object_unref(dataline);
g_free(units);
return FALSE;
}
static void
line_correct_match(GwyContainer *data,
GwyRunType run)
{
GwyDataField *dfield;
GwyDataLine *shifts;
gint xres, yres, i;
gdouble *d, *s;
GQuark dquark;
g_return_if_fail(run & LINECORR_RUN_MODES);
gwy_app_data_browser_get_current(GWY_APP_DATA_FIELD, &dfield,
GWY_APP_DATA_FIELD_KEY, &dquark,
0);
g_return_if_fail(dfield && dquark);
gwy_app_undo_qcheckpointv(data, 1, &dquark);
yres = gwy_data_field_get_yres(dfield);
xres = gwy_data_field_get_xres(dfield);
d = gwy_data_field_get_data(dfield);
shifts = gwy_data_line_new(yres, 1.0, TRUE);
s = gwy_data_line_get_data(shifts);
s[0] = 0.0;
for (i = 1; i < yres; i++) {
s[i] = find_shift(xres, d + i*xres, d + (i - 1)*xres);
g_printerr("%d %g\n", i, s[i]);
}
gwy_data_line_cumulate(shifts);
for (i = 1; i < yres; i++)
gwy_data_field_area_add(dfield, 0, i, xres, 1, s[i]);
gwy_data_field_add(dfield, -s[yres-1]/(xres*yres));
g_object_unref(shifts);
gwy_data_field_data_changed(dfield);
}
int main(int argc, char *argv[])
{
GdkPixbuf *pxb;
GError *error=NULL;
guchar *buffer;
gsize size, pos;
FILE *fh;
gint i, j;
GwyDataField *a, *b, *c, *d, *p;
GwyDataLine *k, *m, n;
GwyPalette *pal;
g_type_init();
g_message("preparing data_field...");
a = (GwyDataField *) gwy_data_field_new(512, 512, 512, 512, 1);
b = (GwyDataField *) gwy_data_field_new(512, 512, 512, 512, 1);
c = (GwyDataField *) gwy_data_field_new(512, 512, 512, 512, 1);
d = (GwyDataField *) gwy_data_field_new(512, 512, 512, 512, 1);
p = (GwyDataField *) gwy_data_field_new(512, 512, 512, 512, 1);
k = (GwyDataLine *) gwy_data_line_new(20, 20, 1);
m = (GwyDataLine *) gwy_data_line_new(500, 20, 1);
gwy_data_line_initialize(&n, 500, 20, 1);
make_test_image(a);
make_test_image(p);
make_test_line(m);
// gwy_data_line_hhcf(m, k);
gwy_data_field_get_line_stat_function(a, &n,
0, 0, 300, 300, /*GwySFOutputType*/ GWY_SF_OUTPUT_PSDF, /*GtkOrientation*/ 0,
/*interpolation*/ 1, /*GwyWindowingType*/ 1, /*nstats*/ 200);
printf("Dataline:\n");
for (i=0; i<10; i++) printf("%d %f\n", i, n.data[i]);
/*
gwy_data_field_2dfft(a, b, c, d, gwy_data_line_fft_hum, GWY_WINDOWING_RECT, 1, GWY_INTERPOLATION_BILINEAR, 0, 0);
gwy_data_field_copy(c, a);
gwy_data_field_copy(d, b);
gwy_data_field_fill(c, 0);
gwy_data_field_fill(d, 0);
gwy_data_field_2dfft(a, b, c, d, gwy_data_line_fft_hum, GWY_WINDOWING_RECT, -1, GWY_INTERPOLATION_BILINEAR, 0, 0);
*/
/* for (i=0; i<512; i++)
{
printf("%f %f\n", a->data[20 + 512*i], //sqrt(a->data[20 + 512*i]*a->data[20 + 512*i]+ b->data[20 + 512*i]*b->data[20 + 512*i]),
p->data[20 + 512*i]);
}
*/
/*
gwy_palette_def_setup_presets();
pal = (GwyPalette*) gwy_palette_new(NULL);
gwy_palette_set_by_name(pal, GWY_PALETTE_OLIVE);
pxb = gdk_pixbuf_new(GDK_COLORSPACE_RGB, FALSE, 8, a->xres, a->yres);
gwy_pixfield_do(pxb, a, pal);
gdk_pixbuf_save(pxb, "a.jpg", "jpeg", &error, "quality", "100", NULL);
gwy_pixfield_do(pxb, b, pal);
gdk_pixbuf_save(pxb, "b.jpg", "jpeg", &error, "quality", "100", NULL);
gwy_pixfield_do(pxb, c, pal);
gdk_pixbuf_save(pxb, "c.jpg", "jpeg", &error, "quality", "100", NULL);
gwy_pixfield_do(pxb, d, pal);
gdk_pixbuf_save(pxb, "d.jpg", "jpeg", &error, "quality", "100", NULL);
*/
return 0;
}
static void
prof_dialog(ProfArgs *args,
GwyContainer *data,
GwyDataField *dfield,
gint id)
{
GtkWidget *hbox, *hbox2, *vbox, *label;
GwyDataChooser *chooser;
GtkTable *table;
GtkDialog *dialog;
GwyGraph *graph;
GwyGraphArea *area;
ProfControls controls;
gint response;
gint row;
controls.args = args;
controls.psdffield = prof_psdf(dfield, args);
controls.modfield = prof_sqrt(controls.psdffield);
controls.hx = gwy_data_field_get_xmeasure(dfield)/(2*G_PI);
controls.hy = gwy_data_field_get_xmeasure(dfield)/(2*G_PI);
controls.dialog = gtk_dialog_new_with_buttons(_("Radial PSDF Section"),
NULL, 0,
GTK_STOCK_CLEAR,
RESPONSE_CLEAR,
GTK_STOCK_CANCEL,
GTK_RESPONSE_CANCEL,
GTK_STOCK_OK,
GTK_RESPONSE_OK,
NULL);
dialog = GTK_DIALOG(controls.dialog);
gtk_dialog_set_default_response(dialog, GTK_RESPONSE_OK);
gtk_dialog_set_response_sensitive(dialog, GTK_RESPONSE_OK, FALSE);
gwy_help_add_to_proc_dialog(GTK_DIALOG(dialog), GWY_HELP_DEFAULT);
hbox = gtk_hbox_new(FALSE, 2);
gtk_box_pack_start(GTK_BOX(GTK_DIALOG(dialog)->vbox), hbox,
FALSE, FALSE, 4);
controls.mydata = gwy_container_new();
gwy_container_set_object_by_name(controls.mydata, "/0/data",
controls.modfield);
g_object_unref(controls.modfield);
gwy_container_set_string_by_name(controls.mydata, "/0/base/palette",
g_strdup("DFit"));
gwy_container_set_enum_by_name(controls.mydata, "/0/base/range-type",
GWY_LAYER_BASIC_RANGE_AUTO);
gwy_app_sync_data_items(data, controls.mydata, id, 0, FALSE,
GWY_DATA_ITEM_REAL_SQUARE,
0);
controls.view = create_preview(controls.mydata, 0, PREVIEW_SIZE, FALSE);
controls.selection = create_vector_layer(GWY_DATA_VIEW(controls.view), 0,
"Point", TRUE);
gwy_selection_set_max_objects(controls.selection, NLINES);
g_object_set(gwy_data_view_get_top_layer(GWY_DATA_VIEW(controls.view)),
"draw-as-vector", TRUE,
NULL);
g_signal_connect_swapped(controls.selection, "changed",
G_CALLBACK(prof_selection_changed), &controls);
gtk_box_pack_start(GTK_BOX(hbox), controls.view, FALSE, FALSE, 4);
vbox = gtk_vbox_new(FALSE, 0);
gtk_box_pack_start(GTK_BOX(hbox), vbox, TRUE, TRUE, 4);
controls.gmodel = gwy_graph_model_new();
g_object_set(controls.gmodel,
"title", _("PSDF Section"),
"axis-label-bottom", "k",
"axis-label-left", "W",
NULL);
init_graph_model_units(controls.gmodel, controls.psdffield);
controls.line = gwy_data_line_new(1, 1.0, FALSE);
controls.graph = gwy_graph_new(controls.gmodel);
graph = GWY_GRAPH(controls.graph);
gtk_widget_set_size_request(controls.graph, GRAPH_WIDTH, -1);
gwy_graph_set_axis_visible(graph, GTK_POS_LEFT, FALSE);
gwy_graph_set_axis_visible(graph, GTK_POS_RIGHT, FALSE);
gwy_graph_set_axis_visible(graph, GTK_POS_TOP, FALSE);
gwy_graph_set_axis_visible(graph, GTK_POS_BOTTOM, FALSE);
gwy_graph_enable_user_input(graph, FALSE);
area = GWY_GRAPH_AREA(gwy_graph_get_area(graph));
gwy_graph_area_enable_user_input(area, FALSE);
gtk_box_pack_start(GTK_BOX(vbox), controls.graph, TRUE, TRUE, 0);
table = GTK_TABLE(gtk_table_new(4, 4, FALSE));
gtk_table_set_row_spacings(table, 2);
gtk_table_set_col_spacings(table, 6);
gtk_container_set_border_width(GTK_CONTAINER(table), 4);
gtk_box_pack_start(GTK_BOX(vbox), GTK_WIDGET(table), FALSE, FALSE, 0);
row = 0;
controls.resolution = gtk_adjustment_new(controls.args->resolution,
MIN_RESOLUTION, MAX_RESOLUTION,
1, 10, 0);
gwy_table_attach_hscale(GTK_WIDGET(table), row, _("_Fix res.:"), NULL,
controls.resolution,
GWY_HSCALE_CHECK | GWY_HSCALE_SQRT);
g_signal_connect_swapped(controls.resolution, "value-changed",
G_CALLBACK(prof_resolution_changed), &controls);
controls.fixres = gwy_table_hscale_get_check(controls.resolution);
gtk_toggle_button_set_active(GTK_TOGGLE_BUTTON(controls.fixres),
controls.args->fixres);
g_signal_connect_swapped(controls.fixres, "toggled",
G_CALLBACK(prof_fixres_changed), &controls);
gwy_table_hscale_set_sensitive(controls.resolution, controls.args->fixres);
row++;
controls.separate = gtk_check_button_new_with_mnemonic(_("_Separate curves"));
gtk_toggle_button_set_active(GTK_TOGGLE_BUTTON(controls.separate),
args->separate);
gtk_table_attach(table, controls.separate,
0, 3, row, row+1, GTK_EXPAND | GTK_FILL, 0, 0, 0);
g_signal_connect_swapped(controls.separate, "toggled",
G_CALLBACK(prof_separate_changed), &controls);
row++;
hbox2 = gtk_hbox_new(FALSE, 6);
gtk_table_attach(table, hbox2,
0, 3, row, row+1, GTK_EXPAND | GTK_FILL, 0, 0, 0);
label = gtk_label_new_with_mnemonic(_("_Interpolation type:"));
gtk_misc_set_alignment(GTK_MISC(label), 0.0, 0.5);
gtk_box_pack_start(GTK_BOX(hbox2), label, FALSE, FALSE, 0);
controls.interpolation
= gwy_enum_combo_box_new(gwy_interpolation_type_get_enum(), -1,
G_CALLBACK(prof_interpolation_changed),
&controls,
controls.args->interpolation, TRUE);
gtk_label_set_mnemonic_widget(GTK_LABEL(label), controls.interpolation);
gtk_box_pack_end(GTK_BOX(hbox2), controls.interpolation, FALSE, FALSE, 0);
row++;
controls.target_hbox = hbox2 = gtk_hbox_new(FALSE, 6);
gtk_widget_set_sensitive(hbox2, !args->separate);
gtk_table_attach(GTK_TABLE(table), hbox2,
0, 4, row, row+1, GTK_EXPAND | GTK_FILL, 0, 0, 0);
label = gtk_label_new_with_mnemonic(_("Target _graph:"));
gtk_misc_set_alignment(GTK_MISC(label), 0.0, 0.5);
gtk_box_pack_start(GTK_BOX(hbox2), label, FALSE, FALSE, 0);
controls.target_graph = gwy_data_chooser_new_graphs();
chooser = GWY_DATA_CHOOSER(controls.target_graph);
gwy_data_chooser_set_none(chooser, _("New graph"));
gwy_data_chooser_set_active(chooser, NULL, -1);
gwy_data_chooser_set_filter(chooser, filter_target_graphs, &controls, NULL);
gwy_data_chooser_set_active_id(chooser, &args->target_graph);
gwy_data_chooser_get_active_id(chooser, &args->target_graph);
gtk_label_set_mnemonic_widget(GTK_LABEL(label), controls.target_graph);
gtk_box_pack_end(GTK_BOX(hbox2), controls.target_graph, FALSE, FALSE, 0);
g_signal_connect_swapped(controls.target_graph, "changed",
G_CALLBACK(target_graph_changed), &controls);
row++;
gtk_widget_show_all(controls.dialog);
do {
response = gtk_dialog_run(dialog);
switch (response) {
case GTK_RESPONSE_CANCEL:
case GTK_RESPONSE_DELETE_EVENT:
gtk_widget_destroy(controls.dialog);
case GTK_RESPONSE_NONE:
goto finalize;
break;
case GTK_RESPONSE_OK:
break;
case RESPONSE_CLEAR:
gwy_selection_clear(controls.selection);
break;
default:
g_assert_not_reached();
break;
}
} while (response != GTK_RESPONSE_OK);
gtk_widget_destroy(controls.dialog);
prof_execute(&controls, data);
finalize:
g_object_unref(controls.psdffield);
g_object_unref(controls.mydata);
g_object_unref(controls.gmodel);
g_object_unref(controls.line);
}
static void
dwt_anisotropy(GwyContainer *data, GwyRunType run)
{
GtkWidget *dialog;
GwyDataField *dfield, *mask;
GQuark dquark, mquark;
GwyDataLine *wtcoefs;
DWTAnisotropyArgs args;
gboolean ok;
gint xsize, ysize, newsize, limit, id, i;
g_return_if_fail(run & DWT_ANISOTROPY_RUN_MODES);
gwy_app_data_browser_get_current(GWY_APP_DATA_FIELD_KEY, &dquark,
GWY_APP_DATA_FIELD, &dfield,
GWY_APP_DATA_FIELD_ID, &id,
GWY_APP_MASK_FIELD_KEY, &mquark,
GWY_APP_MASK_FIELD, &mask,
0);
g_return_if_fail(dfield && dquark);
xsize = gwy_data_field_get_xres(dfield);
ysize = gwy_data_field_get_yres(dfield);
if (xsize != ysize) {
dialog = gtk_message_dialog_new
(gwy_app_find_window_for_channel(data, id),
GTK_DIALOG_DESTROY_WITH_PARENT,
GTK_MESSAGE_ERROR,
GTK_BUTTONS_OK,
_("%s: Data must be square."), _("DWT Anisotropy"));
gtk_dialog_run(GTK_DIALOG(dialog));
gtk_widget_destroy(dialog);
return;
}
dwt_anisotropy_load_args(gwy_app_settings_get(), &args);
if (run == GWY_RUN_INTERACTIVE) {
ok = dwt_anisotropy_dialog(&args);
dwt_anisotropy_save_args(gwy_app_settings_get(), &args);
if (!ok)
return;
}
for (newsize = 1, i = xsize-1; i; i >>= 1, newsize <<= 1)
;
dfield = gwy_data_field_new_resampled(dfield, newsize, newsize,
args.interp);
gwy_data_field_add(dfield, -gwy_data_field_get_avg(dfield));
gwy_app_undo_qcheckpoint(data, dquark, mquark, 0);
if (!mask) {
mask = gwy_data_field_new_alike(dfield, FALSE);
gwy_container_set_object(data, mquark, mask);
g_object_unref(mask);
}
gwy_data_field_resample(mask, newsize, newsize, GWY_INTERPOLATION_NONE);
wtcoefs = gwy_data_line_new(10, 10, TRUE);
wtcoefs = gwy_dwt_set_coefficients(wtcoefs, args.wavelet);
/*justo for sure clamp the lowlimit again*/
limit = pow(2, CLAMP(args.lowlimit, 1, 20));
gwy_data_field_dwt_mark_anisotropy(dfield, mask, wtcoefs, args.ratio,
limit);
gwy_data_field_resample(mask, xsize, ysize, GWY_INTERPOLATION_ROUND);
g_object_unref(wtcoefs);
g_object_unref(dfield);
gwy_data_field_data_changed(mask);
gwy_app_channel_log_add_proc(data, id, id);
}
static void
preview(EntropyControls *controls)
{
EntropyArgs *args = controls->args;
GwyDataField *dfield = controls->dfield;
GwyDataField *mfield = controls->mfield;
GwyGraphModel *gmodel;
GwyGraphCurveModel *gcmodel;
GwyDataLine *ecurve;
gchar buf[24];
gdouble S, s, Smax = 0.0;
ecurve = gwy_data_line_new(1, 1.0, FALSE);
if (args->mode == ENTROPY_VALUES) {
S = gwy_data_field_area_get_entropy_at_scales(dfield, ecurve,
mfield, args->masking,
0, 0,
dfield->xres,
dfield->yres,
0);
s = gwy_data_field_area_get_rms_mask(dfield,
mfield, args->masking,
0, 0, dfield->xres, dfield->yres);
Smax = ENTROPY_NORMAL + log(s);
}
else {
GwyDataField *xder = gwy_data_field_new_alike(dfield, FALSE);
GwyDataField *yder = gwy_data_field_new_alike(dfield, FALSE);
compute_slopes(controls->dfield,
args->fit_plane ? args->kernel_size : 0, xder, yder);
xder = fake_mask(xder, mfield, args->masking);
yder = fake_mask(yder, mfield, args->masking);
if (args->mode == ENTROPY_ANGLES)
transform_to_sphere(xder, yder);
S = gwy_data_field_get_entropy_2d_at_scales(xder, yder, ecurve, 0);
if (args->mode == ENTROPY_SLOPES) {
s = calculate_sigma2_2d(xder, yder);
Smax = ENTROPY_NORMAL_2D + log(s);
}
g_object_unref(xder);
g_object_unref(yder);
}
g_snprintf(buf, sizeof(buf), "%g", S);
gtk_label_set_text(GTK_LABEL(controls->entropy), buf);
if (args->mode != ENTROPY_ANGLES) {
g_snprintf(buf, sizeof(buf), "%g", Smax - S);
gtk_label_set_text(GTK_LABEL(controls->entropydef), buf);
}
else
gtk_label_set_text(GTK_LABEL(controls->entropydef), _("N.A."));
gmodel = gwy_graph_get_model(GWY_GRAPH(controls->graph));
gwy_graph_model_remove_all_curves(gmodel);
g_object_set(gmodel,
"axis-label-bottom", "log h",
"axis-label-left", "S",
"label-position", GWY_GRAPH_LABEL_NORTHWEST,
NULL);
if (gwy_data_line_get_min(ecurve) > -0.5*G_MAXDOUBLE) {
gcmodel = gwy_graph_curve_model_new();
g_object_set(gcmodel,
"description", _("Entropy at scales"),
"mode", GWY_GRAPH_CURVE_LINE_POINTS,
"color", gwy_graph_get_preset_color(0),
NULL);
gwy_graph_curve_model_set_data_from_dataline(gcmodel, ecurve, 0, 0);
gwy_graph_model_add_curve(gmodel, gcmodel);
g_object_unref(gcmodel);
}
if (S > -0.5*G_MAXDOUBLE) {
GwyDataLine *best = gwy_data_line_duplicate(ecurve);
gdouble *ydata = gwy_data_line_get_data(best);
guint i, res = gwy_data_line_get_res(best);
for (i = 0; i < res; i++)
ydata[i] = S;
gcmodel = gwy_graph_curve_model_new();
g_object_set(gcmodel,
"description", _("Best estimate"),
"mode", GWY_GRAPH_CURVE_LINE,
"color", gwy_graph_get_preset_color(1),
NULL);
gwy_graph_curve_model_set_data_from_dataline(gcmodel, best, 0, 0);
gwy_graph_model_add_curve(gmodel, gcmodel);
g_object_unref(gcmodel);
g_object_unref(best);
}
g_object_unref(ecurve);
zoom_in_changed(controls, GTK_TOGGLE_BUTTON(controls->zoom_in));
}
static GwySpectra*
rhkspm32_read_spectra(RHKPage *rhkpage)
{
guint i, j;
gdouble *data;
GwySIUnit *siunit = NULL;
GwyDataLine *dline;
GwySpectra *spectra = NULL;
GPtrArray *spectrum = NULL;
// i'm leaving this alone, though it probably doesn't make sense,
// and i should just create graphs straight away - but in case of
// future use, i'll just convert the data later to graphs
// xres stores number of data points per spectra,
// yres stores the number of spectra
// reading data
gwy_debug("rhk-spm32: %d spectra in this page\n", rhkpage->yres);
for (i = 0; i < rhkpage->yres; i++) {
dline = gwy_data_line_new(rhkpage->xres, rhkpage->x.scale, FALSE);
gwy_data_line_set_offset(dline, (rhkpage->x.offset));
data = gwy_data_line_get_data(dline);
// store line data in physical units - which are the z values, not y
if ((rhkpage->data_type) == RHK_DATA_INT16) {
const guint16 *p = (const guint16*)(rhkpage->buffer
+ rhkpage->data_offset);
for (j = 0; j < rhkpage->xres; j++) {
data[j] = GINT16_FROM_LE(p[i*(rhkpage->xres) + j])
*(rhkpage->z.scale)+(rhkpage->z.offset);
}
}
else if ((rhkpage->data_type) == RHK_DATA_SINGLE) {
const guchar *p = (const guchar*)(rhkpage->buffer
+ rhkpage->data_offset);
for (j = 0; j < rhkpage->xres; j++) {
data[j] = gwy_get_gfloat_le(&p)*rhkpage->z.scale
+ rhkpage->z.offset;
}
}
siunit = gwy_si_unit_new(rhkpage->x.units);
gwy_data_line_set_si_unit_x(dline, siunit);
g_object_unref(siunit);
// the y units (and data) for a 1D graph are stored in Z in the rhk
// spm32 format!
/* Fix "/\xfbHz" to "/Hz".
* XXX: It might be still wrong as the strange character might mean
* sqrt. */
if (g_str_has_suffix(rhkpage->z.units, "/\xfbHz")) {
gchar *s = gwy_strkill(g_strdup(rhkpage->z.units), "\xfb");
siunit = gwy_si_unit_new(s);
g_free(s);
}
else
siunit = gwy_si_unit_new(rhkpage->z.units);
gwy_data_line_set_si_unit_y(dline, siunit);
g_object_unref(siunit);
if (!spectrum)
spectrum = g_ptr_array_sized_new(rhkpage->yres);
g_ptr_array_add(spectrum, dline);
}
gwy_debug("rhk-spm32: finished parsing sps data\n");
spectra = gwy_spectra_new();
for (i = 0; i < rhkpage->yres; i++) {
dline = g_ptr_array_index(spectrum, i);
// since RHK spm32 does not record where it took the spectra,
// i'm setting these to zero
gwy_spectra_add_spectrum(spectra, dline, 0, 0);
g_object_unref(dline);
}
gwy_spectra_set_title(spectra, rhkpage->label);
if (spectrum)
g_ptr_array_free(spectrum, TRUE);
return spectra;
}
static GwySpectra*
omicron_read_cs_data(OmicronFile *ofile,
OmicronSpectroChannel *channel,
GError **error)
{
GError *err = NULL;
GwyDataLine *dline;
GwySIUnit *siunit = NULL, *coord_unit = NULL;
GwySpectra *spectra = NULL;
GPtrArray *spectrum = NULL;
gchar *filename;
gdouble *data, x, y;
gdouble *coords = NULL;
gchar *buffer;
gdouble scale;
guint i, j;
gint power10 = 0;
gint ncurves = 0;
gchar* line;
filename = omicron_fix_file_name(ofile->filename, channel->filename, error);
if (!filename)
return NULL;
gwy_debug("Succeeded with <%s>", filename);
if (!g_file_get_contents(filename, &buffer, NULL , &err)) {
g_free(filename);
err_GET_FILE_CONTENTS(error, &err);
return NULL;
}
g_free(filename);
scale = channel->resolution; /* can also be extracted from min&max
raw and phys settings */
while ((line = gwy_str_next_line(&buffer))) {
if (strstr(line, ";n_curves")) {
/* Find number of curves this should appear first in file */
ncurves = g_ascii_strtod(strchr(line, ':')+1, NULL);
}
if (strstr(line, "BEGIN COORD")) {
/* Read in cordinates Spectroscopy Curves */
i = 0;
coord_unit = gwy_si_unit_new_parse("nm", &power10);
while ((line = gwy_str_next_line(&buffer))) {
gchar *val2;
if (strstr(line, "END")) {
if (i != ncurves) {
gwy_debug("Less coords than ncurves");
}
break;
}
if (i == ncurves) {
g_critical("More coords than ncurves.");
break;
}
if (!coords) {
if (!(coords = g_new0(gdouble, ncurves*2))) {
gwy_debug("Failed to allocate mem: coords");
return NULL;
}
}
val2 = line+16;
x = g_ascii_strtod(line, &val2) * pow10(power10);
y = g_ascii_strtod(val2, NULL) * pow10(power10);
gwy_debug("Coord %i: x:%g y:%g", i, x, y);
coords[2*i] = x;
coords[2*i+1] = y;
i++;
}
/* i is set to 0 and used as a counter for the dline */
i = 0;
}
if (strstr(line, "BEGIN") && !strstr(line, "COORD")) {
/* Read spectroscopy points */
dline = gwy_data_line_new(channel->npoints,
channel->end - channel->start,
FALSE);
gwy_data_line_set_offset(dline, (channel->start));
data = gwy_data_line_get_data(dline);
j = 0;
while ((line = gwy_str_next_line(&buffer))) {
gchar *val2;
if (strstr(line, "END") || j >= channel->npoints)
break;
val2 = line+13;
x = g_ascii_strtod(line, &val2);
y = g_ascii_strtod(val2, NULL)*scale;
data[j] = y;
j++;
}
/* Set Units for the parameter (x) axis */
if ((channel->param[0] == 'V') || (channel->param[0] == 'E')) {
siunit = gwy_si_unit_new("V");
power10 = 0;
}
else if (channel->param[0] == 'I')
siunit = gwy_si_unit_new_parse("nA", &power10);
else if (channel->param[0] == 'Z')
siunit = gwy_si_unit_new_parse("nm", &power10);
else {
gwy_debug("Parameter unit not recognised");
}
if (siunit) {
gwy_data_line_set_si_unit_x(dline, siunit);
g_object_unref(siunit);
}
if (power10) {
gdouble offset = 0;
gdouble realsize = 0;
offset = gwy_data_line_get_offset(dline)*pow10(power10);
realsize = gwy_data_line_get_real(dline)*pow10(power10);
gwy_data_line_set_offset(dline, offset);
gwy_data_line_set_real(dline, realsize);
}
/* Set Units for the Value (y) Axis */
siunit = gwy_si_unit_new_parse(channel->units, &power10);
gwy_data_line_set_si_unit_y(dline, siunit);
g_object_unref(siunit);
if (power10)
gwy_data_line_multiply(dline, pow10(power10));
if (!spectrum)
spectrum = g_ptr_array_sized_new(ncurves);
g_ptr_array_add(spectrum, dline);
}
}
if (!spectrum)
spectrum = g_ptr_array_new();
if (spectrum->len < ncurves) {
gwy_debug("Less actual spectra than ncurves");
ncurves = spectrum->len;
}
if (spectrum->len > ncurves) {
gwy_debug("More actual spectra than ncurves, "
"remaining pos will be set at (0.0,0.0)");
coords = g_renew(gdouble, coords, spectrum->len*2);
if (!coords) {
g_critical("Could not reallocate mem for coords.");
return NULL;
}
while (spectrum->len > ncurves) {
coords[ncurves*2] = 0.0;
coords[ncurves*2+1] = 0.0;
ncurves++;
}
}
spectra = gwy_spectra_new();
if (coord_unit) {
gwy_spectra_set_si_unit_xy(spectra, coord_unit);
g_object_unref(coord_unit);
}
for (i = 0; i < ncurves; i++) {
dline = g_ptr_array_index(spectrum, i);
gwy_spectra_add_spectrum(spectra, dline,
coords[i*2], ofile->yreal - coords[i*2+1]);
g_object_unref(dline);
}
g_ptr_array_free(spectrum, TRUE);
g_free(coords);
g_free(buffer);
return spectra;
}
static void
line_correct_match(GwyContainer *data,
GwyRunType run)
{
GwyDataField *dfield;
GwyDataLine *shifts;
gint xres, yres, i, j;
gdouble m, wsum, lambda, x;
gdouble *d, *s, *w;
const gdouble *a, *b;
GQuark dquark;
g_return_if_fail(run & LINECORR_RUN_MODES);
gwy_app_data_browser_get_current(GWY_APP_DATA_FIELD, &dfield,
GWY_APP_DATA_FIELD_KEY, &dquark,
0);
g_return_if_fail(dfield && dquark);
gwy_app_undo_qcheckpointv(data, 1, &dquark);
yres = gwy_data_field_get_yres(dfield);
xres = gwy_data_field_get_xres(dfield);
d = gwy_data_field_get_data(dfield);
shifts = gwy_data_line_new(yres, 1.0, TRUE);
s = gwy_data_line_get_data(shifts);
w = g_new(gdouble, xres-1);
for (i = 1; i < yres; i++) {
a = d + xres*(i - 1);
b = d + xres*i;
/* Diffnorm */
wsum = 0.0;
for (j = 0; j < xres-1; j++) {
x = a[j+1] - a[j] - b[j+1] + b[j];
wsum += fabs(x);
}
if (wsum == 0)
continue;
m = wsum/(xres-1);
/* Weights */
wsum = 0.0;
for (j = 0; j < xres-1; j++) {
x = a[j+1] - a[j] - b[j+1] + b[j];
w[j] = exp(-(x*x/(2.0*m)));
wsum += w[j];
}
/* Correction */
lambda = (a[0] - b[0])*w[0];
for (j = 1; j < xres-1; j++)
lambda += (a[j] - b[j])*(w[j-1] + w[j]);
lambda += (a[xres-1] - b[xres-1])*w[xres-2];
lambda /= 2.0*wsum;
gwy_debug("%g %g %g", m, wsum, lambda);
s[i] = lambda;
}
gwy_data_line_cumulate(shifts);
for (i = 1; i < yres; i++)
gwy_data_field_area_add(dfield, 0, i, xres, 1, s[i]);
gwy_data_field_add(dfield, -s[yres-1]/(xres*yres));
g_object_unref(shifts);
g_free(w);
gwy_data_field_data_changed(dfield);
}
static void
dwt_denoise(GwyContainer *data, GwyRunType run)
{
GtkWidget *dialog;
GwyDataField *dfield;
GwyDataLine *wtcoefs;
DWTDenoiseArgs args;
gboolean ok;
gint xsize, ysize, newsize;
gint oldid, newid;
g_return_if_fail(run & DWT_DENOISE_RUN_MODES);
gwy_app_data_browser_get_current(GWY_APP_DATA_FIELD, &dfield,
GWY_APP_DATA_FIELD_ID, &oldid,
0);
g_return_if_fail(dfield);
xsize = gwy_data_field_get_xres(dfield);
ysize = gwy_data_field_get_yres(dfield);
if (xsize != ysize) {
dialog = gtk_message_dialog_new
(GTK_WINDOW(gwy_app_data_window_get_current()),
GTK_DIALOG_DESTROY_WITH_PARENT,
GTK_MESSAGE_ERROR,
GTK_BUTTONS_OK,
_("%s: Data must be square."), _("DWT Denoise"));
gtk_dialog_run(GTK_DIALOG(dialog));
gtk_widget_destroy(dialog);
return;
}
dwt_denoise_load_args(gwy_app_settings_get(), &args);
if (run == GWY_RUN_INTERACTIVE) {
ok = dwt_denoise_dialog(&args);
dwt_denoise_save_args(gwy_app_settings_get(), &args);
if (!ok)
return;
}
dfield = gwy_data_field_duplicate(dfield);
newsize = gwy_fft_find_nice_size(xsize);
gwy_data_field_add(dfield, -gwy_data_field_get_avg(dfield));
gwy_data_field_resample(dfield, newsize, newsize,
GWY_INTERPOLATION_BILINEAR);
wtcoefs = gwy_data_line_new(10, 10, TRUE);
wtcoefs = gwy_dwt_set_coefficients(wtcoefs, args.wavelet);
gwy_data_field_dwt_denoise(dfield, wtcoefs, TRUE, 20, args.method);
if (args.preserve)
gwy_data_field_resample(dfield, xsize, ysize, args.interp);
newid = gwy_app_data_browser_add_data_field(dfield, data, TRUE);
gwy_app_copy_data_items(data, data, oldid, newid,
GWY_DATA_ITEM_GRADIENT,
GWY_DATA_ITEM_MASK_COLOR,
0);
g_object_unref(dfield);
gwy_app_set_data_field_title(data, newid, _("DWT denoised"));
g_object_unref(wtcoefs);
}
static void
unrotate(GwyContainer *data, GwyRunType run)
{
enum { nder = 4800 };
GwyDataField *dfield, *mfield, *sfield;
UnrotateArgs args;
gdouble correction[GWY_SYMMETRY_LAST];
GwyPlaneSymmetry symm;
GwyDataLine *derdist;
GQuark dquark, mquark, squark;
gdouble phi;
gboolean ok = TRUE;
gint id;
g_return_if_fail(run & UNROTATE_RUN_MODES);
gwy_app_data_browser_get_current(GWY_APP_DATA_FIELD_KEY, &dquark,
GWY_APP_DATA_FIELD, &dfield,
GWY_APP_DATA_FIELD_ID, &id,
GWY_APP_MASK_FIELD_KEY, &mquark,
GWY_APP_MASK_FIELD, &mfield,
GWY_APP_SHOW_FIELD_KEY, &squark,
GWY_APP_SHOW_FIELD, &sfield,
0);
g_return_if_fail(dfield && dquark && mquark && squark);
load_args(gwy_app_settings_get(), &args);
derdist = GWY_DATA_LINE(gwy_data_line_new(nder, 2*G_PI, FALSE));
gwy_data_field_slope_distribution(dfield, derdist, 5);
symm = gwy_data_field_unrotate_find_corrections(derdist, correction);
g_object_unref(derdist);
if (run == GWY_RUN_INTERACTIVE) {
ok = unrotate_dialog(&args, data, dfield, id, correction, symm);
save_args(gwy_app_settings_get(), &args);
if (!ok)
return;
}
if (args.symmetry)
symm = args.symmetry;
phi = correction[symm];
if (!mfield)
mquark = 0;
if (!sfield)
squark = 0;
if (!mfield && sfield)
GWY_SWAP(GQuark, mquark, squark);
gwy_app_undo_qcheckpoint(data, dquark, mquark, squark, 0);
gwy_data_field_rotate(dfield, phi, args.interp);
gwy_data_field_data_changed(dfield);
if (mfield) {
gwy_data_field_rotate(mfield, phi, GWY_INTERPOLATION_ROUND);
gwy_data_field_data_changed(mfield);
}
if (sfield) {
gwy_data_field_rotate(sfield, phi, args.interp);
gwy_data_field_data_changed(sfield);
}
gwy_app_channel_log_add_proc(data, id, id);
}
static gboolean
curvature_plot_graph(GwyDataField *dfield,
const Intersection *i1,
const Intersection *i2,
GwyGraphModel *gmodel)
{
GwyGraphCurveModel *gcmodel;
GwyDataLine *dline;
gint xres, yres;
guint i;
if (!gwy_graph_model_get_n_curves(gmodel)) {
GwySIUnit *siunitxy, *siunitz;
gchar *s;
siunitxy = gwy_si_unit_duplicate(gwy_data_field_get_si_unit_xy(dfield));
siunitz = gwy_si_unit_duplicate(gwy_data_field_get_si_unit_z(dfield));
g_object_set(gmodel,
"title", _("Curvature Sections"),
"si-unit-x", siunitxy,
"si-unit-y", siunitz,
NULL);
g_object_unref(siunitxy);
g_object_unref(siunitz);
for (i = 0; i < 2; i++) {
gcmodel = gwy_graph_curve_model_new();
s = g_strdup_printf(_("Profile %d"), (gint)i+1);
g_object_set(gcmodel,
"description", s,
"mode", GWY_GRAPH_CURVE_LINE,
"color", gwy_graph_get_preset_color(i),
NULL);
g_free(s);
gwy_graph_model_add_curve(gmodel, gcmodel);
g_object_unref(gcmodel);
}
}
else {
g_assert(gwy_graph_model_get_n_curves(gmodel) == 2);
}
dline = gwy_data_line_new(1, 1.0, FALSE);
xres = gwy_data_field_get_xres(dfield);
yres = gwy_data_field_get_yres(dfield);
for (i = 0; i < 2; i++) {
gint col1 = gwy_data_field_rtoj(dfield, i1[i].x);
gint row1 = gwy_data_field_rtoi(dfield, i1[i].y);
gint col2 = gwy_data_field_rtoj(dfield, i2[i].x);
gint row2 = gwy_data_field_rtoi(dfield, i2[i].y);
gwy_data_field_get_profile(dfield, dline,
CLAMP(col1, 0, xres-1),
CLAMP(row1, 0, yres-1),
CLAMP(col2, 0, xres-1),
CLAMP(row2, 0, yres-1),
-1, 1, GWY_INTERPOLATION_BILINEAR);
gwy_data_line_set_offset(dline,
i1[i].t/(i2[i].t - i1[i].t)
* gwy_data_line_get_real(dline));
gcmodel = gwy_graph_model_get_curve(gmodel, i);
gwy_graph_curve_model_set_data_from_dataline(gcmodel, dline, 0, 0);
}
g_object_unref(dline);
return TRUE;
}
