A tool to construct HDR-images from a series of exposures.

main.cpp 8.1KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264
  1. #include <math.h>
  2. #include <iostream>
  3. using namespace std;
  4. #include <gtk/gtk.h>
  5. #include <glib.h>
  6. #include "haader.h"
  7. #include "image.h"
  8. struct Gui {
  9. GtkWidget *window;
  10. GtkWidget *canvas;
  11. haader::HdrImageStack stack;
  12. haader::HdrImage hdr_image;
  13. haader::ResponseFunction response_function;
  14. haader::Image ldr_image;
  15. haader::Histogram hdr_histogram;
  16. haader::LdrHistogram ldr_histogram;
  17. cairo_surface_t *surface;
  18. double exposure_time;
  19. double compression;
  20. };
  21. static cairo_surface_t* image_to_cairo_surface(const haader::Image &image) {
  22. unsigned int width = image.get_width();
  23. unsigned int height = image.get_height();
  24. const unsigned char *image_data = image.get_const_image_data();
  25. cairo_surface_t *surf = cairo_image_surface_create(CAIRO_FORMAT_RGB24, width, height);
  26. int stride = cairo_image_surface_get_stride(surf);
  27. unsigned char *surf_data = cairo_image_surface_get_data(surf);
  28. unsigned int i = 0;
  29. unsigned int k = 0;
  30. for (unsigned int y = 0; y < height; y++) {
  31. for (unsigned int x = 0; x < width; x++) {
  32. surf_data[i + x * 4] = image_data[k + x * 3 + 2];
  33. surf_data[i + x * 4 + 1] = image_data[k + x * 3 + 1];
  34. surf_data[i + x * 4 + 2] = image_data[k + x * 3];
  35. }
  36. i += stride;
  37. k += width * 3;
  38. }
  39. return surf;
  40. }
  41. static void expose(Gui &gui) {
  42. gui.ldr_image = gui.hdr_image.expose(gui.exposure_time, gui.response_function, gui.compression);
  43. gui.ldr_histogram = gui.ldr_image.histogram();
  44. if (gui.surface) {
  45. cairo_surface_destroy(gui.surface);
  46. gui.surface = NULL;
  47. }
  48. gui.surface = image_to_cairo_surface(gui.ldr_image);
  49. }
  50. static gboolean canvas_draw(GtkWidget *widget, cairo_t *cr, gpointer user_data)
  51. {
  52. cairo_set_source_rgb(cr, 0.2, 0.2, 0.2);
  53. cairo_paint(cr);
  54. Gui *gui = (Gui*) user_data;
  55. if (!gui->surface) {
  56. return false;
  57. }
  58. if (gui->hdr_histogram.m_bins.size() != 256) {
  59. return false;
  60. }
  61. if (gui->ldr_histogram.m_bins.size() != 256) {
  62. return false;
  63. }
  64. cairo_set_source_surface(cr, gui->surface, 5, 5);
  65. cairo_paint(cr);
  66. // canvas rect
  67. double cx = gui->ldr_image.get_width() + 10.0;
  68. double cy = 5.0;
  69. double cw = 256;
  70. double ch = 256;
  71. cairo_set_source_rgb(cr, 0.0, 0.0, 0.0);
  72. cairo_rectangle(cr, cx, cy, cw, ch);
  73. cairo_fill(cr);
  74. {
  75. double hist_scale = (256.0 / (double)(gui->hdr_histogram.m_max_freq));
  76. for (unsigned int i = 0; i < 256; i++) {
  77. double y = gui->hdr_histogram.m_bins[i] * hist_scale;
  78. cairo_rectangle(cr, cx + i, cy + ch - y, 1, y);
  79. }
  80. cairo_set_source_rgba(cr, 1.0, 0.0, 0.2, 0.4);
  81. cairo_fill(cr);
  82. }
  83. {
  84. double hist_scale = (256.0 / (double)(gui->ldr_histogram.m_max_freq));
  85. for (unsigned int i = 0; i < 256; i++) {
  86. double y = gui->ldr_histogram.m_bins[i] * hist_scale;
  87. cairo_rectangle(cr, cx + i, cy + ch - y, 1, y);
  88. }
  89. cairo_set_source_rgba(cr, 1.0, 0.0, 0.2, 0.9);
  90. cairo_fill(cr);
  91. }
  92. double minl = gui->hdr_histogram.m_min_value;
  93. double maxl = gui->hdr_histogram.m_max_value;
  94. double log_exp = log(gui->exposure_time);
  95. double scale = (maxl - minl) / 256.0;
  96. for (int i = 0; i < 256; i++) {
  97. double lum = minl + i * scale;
  98. int ly = gui->response_function.lookup((lum + log_exp) * gui->compression);
  99. cairo_line_to(cr, cx + i, cy + ch - ly);
  100. }
  101. cairo_set_source_rgba(cr, 1.0, 1.0, 1.0, 0.8);
  102. cairo_set_line_width(cr, 1.0);
  103. cairo_stroke_preserve(cr);
  104. cairo_line_to(cr, cx + cw, cy + ch);
  105. cairo_line_to(cr, cx, cy + ch);
  106. cairo_close_path(cr);
  107. cairo_set_source_rgba(cr, 0.0, 0.5, 1.0, 0.2);
  108. cairo_fill(cr);
  109. return FALSE;
  110. }
  111. void scale_exposure_changed(GtkRange *range, gpointer user_data) {
  112. Gui *gui = (Gui*) user_data;
  113. gui->exposure_time = pow(2.0, gtk_range_get_value(range));
  114. expose(*gui);
  115. gtk_widget_queue_draw(gui->canvas);
  116. }
  117. void scale_compression_changed(GtkRange *range, gpointer user_data) {
  118. Gui *gui = (Gui*) user_data;
  119. gui->compression = pow(2.0, gtk_range_get_value(range));
  120. expose(*gui);
  121. gtk_widget_queue_draw(gui->canvas);
  122. }
  123. void file_select_clicked(GtkButton *button, gpointer user_data) {
  124. Gui *gui = (Gui*) user_data;
  125. GtkWidget *dialog = gtk_file_chooser_dialog_new(
  126. "Open File",
  127. GTK_WINDOW(gui->window),
  128. GTK_FILE_CHOOSER_ACTION_OPEN,
  129. "Cancel",
  130. GTK_RESPONSE_CANCEL,
  131. "Open",
  132. GTK_RESPONSE_ACCEPT,
  133. NULL);
  134. gtk_file_chooser_set_select_multiple (GTK_FILE_CHOOSER(dialog), true);
  135. bool ok = true;
  136. if (gtk_dialog_run(GTK_DIALOG(dialog)) == GTK_RESPONSE_ACCEPT) {
  137. GSList *list = gtk_file_chooser_get_filenames(GTK_FILE_CHOOSER(dialog));
  138. vector<string> new_list;
  139. GSList *element = list;
  140. while (element != NULL) {
  141. ok = ok && gui->stack.add_from_file_path((const char*)(element->data));
  142. g_free(element->data);
  143. element = g_slist_next(element);
  144. if (!ok) {
  145. break;
  146. }
  147. }
  148. g_slist_free(list);
  149. }
  150. if (ok) {
  151. // create HDR-image
  152. haader::InverseResponseFunction irf;
  153. if (gui->stack.get_inverse_response_function(irf, 4096)) {
  154. gui->hdr_image = gui->stack.get_hdr_image(irf);
  155. gui->ldr_image = gui->hdr_image.get_log_image();
  156. gui->hdr_histogram = gui->hdr_image.histogram(256);
  157. gui->ldr_histogram = gui->ldr_image.histogram();
  158. if (gui->surface) {
  159. cairo_surface_destroy(gui->surface);
  160. gui->surface = NULL;
  161. }
  162. gui->surface = image_to_cairo_surface(gui->ldr_image);
  163. gui->response_function = irf.to_response_function(1024);
  164. }
  165. gtk_widget_queue_draw(gui->canvas);
  166. }
  167. gtk_widget_destroy (dialog);
  168. }
  169. int main(int argc, char *argv[])
  170. {
  171. gtk_init(&argc, &argv);
  172. Gui gui;
  173. gui.surface = NULL;
  174. gui.exposure_time = 0.0;
  175. gui.compression = 1.0;
  176. gui.window = gtk_window_new(GTK_WINDOW_TOPLEVEL);
  177. gtk_window_set_title(GTK_WINDOW(gui.window), "haader");
  178. gtk_window_set_default_size(GTK_WINDOW(gui.window), 800, 600);
  179. GtkWidget *file_button = gtk_button_new_with_label("Select files");
  180. gui.canvas = gtk_drawing_area_new();
  181. GtkWidget *scale_exposure = gtk_scale_new_with_range(GTK_ORIENTATION_HORIZONTAL, -20.0, 4.0, 0.001);
  182. gtk_range_set_value(GTK_RANGE(scale_exposure), log(gui.exposure_time) / log(2));
  183. GtkWidget *scale_compression = gtk_scale_new_with_range(GTK_ORIENTATION_HORIZONTAL, -4, 4.0, 0.001);
  184. gtk_range_set_value(GTK_RANGE(scale_compression), log(gui.compression) / log(2));
  185. GtkWidget *vbox = gtk_box_new(GTK_ORIENTATION_VERTICAL, 5);
  186. gtk_container_add(GTK_CONTAINER(gui.window), vbox);
  187. gtk_box_pack_start(GTK_BOX(vbox), file_button, false, false, 0);
  188. gtk_box_pack_start(GTK_BOX(vbox), gui.canvas, true, true, 0);
  189. GtkWidget *grid = gtk_grid_new();
  190. gtk_box_pack_start(GTK_BOX(vbox), grid, false, false, 0);
  191. {
  192. GtkWidget *label = gtk_label_new("Exposure");
  193. gtk_grid_attach(GTK_GRID(grid), label, 0, 0, 1, 1);
  194. gtk_widget_set_size_request(scale_exposure, 600, -1);
  195. gtk_grid_attach(GTK_GRID(grid), scale_exposure, 1, 0, 1, 1);
  196. }
  197. {
  198. GtkWidget *label = gtk_label_new("Compression");
  199. gtk_grid_attach(GTK_GRID(grid), label, 0, 1, 1, 1);
  200. gtk_widget_set_size_request(scale_compression, 600, -1);
  201. gtk_grid_attach(GTK_GRID(grid), scale_compression, 1, 1, 1, 1);
  202. }
  203. g_signal_connect(gui.window, "destroy", G_CALLBACK(gtk_main_quit), NULL);
  204. g_signal_connect(gui.canvas, "draw", G_CALLBACK(canvas_draw), &gui);
  205. g_signal_connect(file_button, "clicked", G_CALLBACK(file_select_clicked), &gui);
  206. g_signal_connect(scale_exposure, "value-changed", G_CALLBACK(scale_exposure_changed), &gui);
  207. g_signal_connect(scale_compression, "value-changed", G_CALLBACK(scale_compression_changed), &gui);
  208. gtk_widget_show_all(gui.window);
  209. gtk_main ();
  210. return 0;
  211. }