1 files changed, 45 insertions, 0 deletions
diff --git a/src/boost/libs/gil/example/sobel_scharr.cpp b/src/boost/libs/gil/example/sobel_scharr.cpp
new file mode 100644
index 000000000..10e15c223
--- /dev/null
+++ b/src/boost/libs/gil/example/sobel_scharr.cpp
@@ -0,0 +1,45 @@
+#include <boost/gil/typedefs.hpp>
+#include <boost/gil/image_processing/numeric.hpp>
+#include <boost/gil/extension/io/png.hpp>
+#include <boost/gil/extension/numeric/convolve.hpp>
+#include <string>
+#include <iostream>
+
+namespace gil = boost::gil;
+
+int main(int argc, char* argv[])
+{
+    if (argc != 5)
+    {
+        std::cerr << "usage: " << argv[0] << ": <input.png> <sobel|scharr> <output-x.png> <output-y.png>\n";
+        return -1;
+    }
+
+    gil::gray8_image_t input_image;
+    gil::read_image(argv[1], input_image, gil::png_tag{});
+    auto input = gil::view(input_image);
+    auto filter_type = std::string(argv[2]);
+
+    gil::gray16_image_t dx_image(input_image.dimensions());
+    auto dx = gil::view(dx_image);
+    gil::gray16_image_t dy_image(input_image.dimensions());
+    auto dy = gil::view(dy_image);
+    if (filter_type == "sobel")
+    {
+        gil::detail::convolve_2d(input, gil::generate_dx_sobel(1), dx);
+        gil::detail::convolve_2d(input, gil::generate_dy_sobel(1), dy);
+    }
+    else if (filter_type == "scharr")
+    {
+        gil::detail::convolve_2d(input, gil::generate_dx_scharr(1), dx);
+        gil::detail::convolve_2d(input, gil::generate_dy_scharr(1), dy);
+    }
+    else
+    {
+        std::cerr << "unrecognized gradient filter type. Must be either sobel or scharr\n";
+        return -1;
+    }
+
+    gil::write_view(argv[3], dx, gil::png_tag{});
+    gil::write_view(argv[4], dy, gil::png_tag{});
+}