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

main.cpp 10KB

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