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

main.cpp 8.6KB

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