* green/mln/io/plot/save_image_sh.hh: Add impl and internal namespace.
---
trunk/milena/sandbox/ChangeLog | 7 +
.../sandbox/green/mln/io/plot/save_image_sh.hh | 1083 ++++++++------------
2 files changed, 431 insertions(+), 659 deletions(-)
diff --git a/trunk/milena/sandbox/ChangeLog b/trunk/milena/sandbox/ChangeLog
index 634be0b..499af4c 100644
--- a/trunk/milena/sandbox/ChangeLog
+++ b/trunk/milena/sandbox/ChangeLog
@@ -1,3 +1,10 @@
+
+2009-10-02 Yann Jacquelet <jacquelet(a)lrde.epita.fr>
+
+ Add impl and internal namespace, correct static dispatch.
+
+ * green/mln/io/plot/save_image_sh.hh: Add impl and internal namespace.
+
2009-10-01 Fabien Freling <fabien.freling(a)lrde.epita.fr>
Update README.
diff --git a/trunk/milena/sandbox/green/mln/io/plot/save_image_sh.hh b/trunk/milena/sandbox/green/mln/io/plot/save_image_sh.hh
index 62f6652..0456d5b 100644
--- a/trunk/milena/sandbox/green/mln/io/plot/save_image_sh.hh
+++ b/trunk/milena/sandbox/green/mln/io/plot/save_image_sh.hh
@@ -33,19 +33,20 @@
/// \brief Define some functions to export to gnuplot format as a script shell.
///
/// Theses routines are dedicated to image visualization. The aim is to display
-/// image whatever the value used as pixel. The behaviour of gnuplot is more
+/// image whatever the value used as pixel. The behaviour of gnuplot is more
/// like xv than imageMagick. A Gnuplot script shell file is a text dump file
/// with a preambule to let gnuplot interpret data. As a script shell, you can
/// launch it (don't forget the permissions), and every thing is packed in it.
/// The script file call gnuplot in batch mode, the result window persists and
/// that's all.
-///
+///
#include <fstream>
#include <string>
#include <mln/trace/entering.hh>
#include <mln/trace/exiting.hh>
+
#include <mln/core/macros.hh>
#include <mln/core/contract.hh>
#include <mln/core/image/image1d.hh>
@@ -56,6 +57,9 @@
#include <mln/geom/max_row.hh>
#include <mln/geom/min_col.hh>
#include <mln/geom/max_col.hh>
+
+#include <mln/trait/value_.hh>
+
#include <mln/value/int_u.hh>
#include <mln/value/int_s.hh>
#include <mln/value/rgb.hh>
@@ -65,16 +69,16 @@
namespace mln
{
-
+
namespace io
{
-
+
namespace plot
{
/// \brief Save an image as a gnuplot script shell.
///
- /// Every thing is save. The image could be 1d, 2d or 3d. The value of
+ /// Every thing is save. The image could be 1d, 2d or 3d. The value of
/// the pixels could be int_u<n>, int_s<n>, float, double, hsl_f, hsl_d,
/// hsi_f, hsi_d and rgb<n>.
///
@@ -83,741 +87,510 @@ namespace mln
/// \return the status of the opening file operation.
///
/// The result depends on the permission to save the file with
- /// filename parameter as unix path. The script shell file must have the
- /// permission to execute (chmod 755). Launch the script shell to call
+ /// filename parameter as unix path. The script shell file must have the
+ /// permission to execute (chmod 755). Launch the script shell to call
/// gnuplot in batchmode with fine parameters.
-
+
template <typename I>
bool save_image_sh(const Image<I>& img, const std::string& filename);
- } // end of namespace mln::accu::stat
-
- } // end of namespace mln::accu
+ } // end of namespace mln::io::plot
+
+ } // end of namespace mln::io
namespace io
{
-
+
namespace plot
{
-
-#ifndef MLN_INCLUDE_ONLY
+#ifndef MLN_INCLUDE_ONLY
//------------------------------------------------------------------------
- // save_image_sh(const image2d<bool>&, const std::string&)
+ // Impl.
//------------------------------------------------------------------------
- inline
- bool save_image_sh_dispatch_overload(const image2d< bool >& img,
- const std::string& filename)
+ namespace impl
{
- trace::entering("mln::io::plot::save_image_sh<image2d<bool>>");
- mln_precondition(img.is_valid());
-
- std::ofstream out(filename.c_str());
- bool result = !out.fail();
- unsigned min_row = geom::min_row(img);
- unsigned max_row = geom::max_row(img);
- unsigned min_col = geom::min_col(img);
- unsigned max_col = geom::max_col(img);
- if (result)
- {
- mln_piter_(image2d< bool >) p(img.domain());
-
- // Output data prelude (terminal X11, image).
- out << "#!/bin/sh" << std::endl;
- out << "####################################" << std::endl;
- out << "# Columns = (x, y, bool) #" << std::endl;
- out << "####################################" << std::endl;
- out << std::endl;
- out << "gnuplot <<EOF" << std::endl;
- out << "set terminal x11 persist 1" << std::endl;
- out << "set palette gray" << std::endl;
- out << "set xrange [" << min_col;
- out << ":" << max_col;
- out << "]" << std::endl;
- out << "set yrange [-" << max_row;
- out << ":" << min_row;
- out << "]" << std::endl;
- out << "plot '-' using (\\$2):(-\\$1):(\\$3) with image" << std::endl;
-
- // Output data.
- for_all(p)
- {
- out << p.row() << " ";
- out << p.col() << " ";
- out << img(p) << std::endl;
- }
-
- // Close gnuplot data stream.
- out << "e" << std::endl;
- out << "EOF" << std::endl;
+ //----------------------------------------------------------------------
+ // save_image_sh_image2d_rgb(const image2d<rgb<n>>&, const string&)
+ //----------------------------------------------------------------------
- out.close();
- }
- else
+ template <unsigned n>
+ inline
+ bool save_image_sh_image2d_rgb(const image2d<value::rgb<n> >& img,
+ const std::string& filename)
{
- std::cerr << "ERROR[mln::io::plot::save_image_sh]:" << filename;
- std::cerr << " couldn't be opened !!" << std::endl;
- }
+ trace::entering("mln::io::plot::impl::save_image_sh_image2d_rgb");
+ mln_precondition(img.is_valid());
- trace::exiting("mln::io::plot::save_image_sh<image2d<bool>>");
- return result;
- }
+ std::ofstream out(filename.c_str());
+ bool result = !out.fail();
+ unsigned min_row = geom::min_row(img);
+ unsigned max_row = geom::max_row(img);
+ unsigned min_col = geom::min_col(img);
+ unsigned max_col = geom::max_col(img);
- //------------------------------------------------------------------------
- // save_image_sh(const image2d<float>&, const std::string&)
- //------------------------------------------------------------------------
-
- inline
- bool save_image_sh_dispatch_overload(const image2d<float>& img,
- const std::string& filename)
- {
- trace::entering("mln::io::plot::save_image_sh<image2d<float>>");
- mln_precondition(img.is_valid());
-
- std::ofstream out(filename.c_str());
- bool result = !out.fail();
- unsigned min_row = geom::min_row(img);
- unsigned max_row = geom::max_row(img);
- unsigned min_col = geom::min_col(img);
- unsigned max_col = geom::max_col(img);
-
- if (result)
- {
- mln_piter_(image2d< float >) p(img.domain());
-
- // Output data prelude (terminal X11, image).
- out << "#!/bin/sh" << std::endl;
- out << "####################################" << std::endl;
- out << "# Columns = (x, y, float) #" << std::endl;
- out << "####################################" << std::endl;
- out << std::endl;
- out << "gnuplot <<EOF" << std::endl;
- out << "set terminal x11 persist 1" << std::endl;
- out << "set palette gray" << std::endl;
- out << "set xrange [" << min_col;
- out << ":" << max_col;
- out << "]" << std::endl;
- out << "set yrange [-" << max_row;
- out << ":" << min_row;
- out << "]" << std::endl;
- out << "plot '-' using (\\$2):(-\\$1):(\\$3) with image" << std::endl;
-
- // Output data.
- for_all(p)
+ if (result)
{
- out << p.row() << " ";
- out << p.col() << " ";
- out << img(p) << std::endl;
+ mln_piter(image2d< value::rgb<n> >) p(img.domain());
+
+ // Output data prelude (terminal X11, image).
+ out << "#!/bin/sh" << std::endl;
+ out << "####################################" << std::endl;
+ out << "# Columns = (x, y, r, g, b) #" << std::endl;
+ out << "####################################" << std::endl;
+ out << std::endl;
+ out << "gnuplot <<EOF" << std::endl;
+ out << "set terminal x11 persist 1" << std::endl;
+ out << "set xrange [" << min_col;
+ out << ":" << max_col;
+ out << "]" << std::endl;
+ out << "set yrange [-" << max_row;
+ out << ":" << min_row;
+ out << "]" << std::endl;
+ out << "plot '-' using 2:(-\\$1):3:4:5 with rgbimage" << std::endl;
+
+ // Output data.
+ for_all(p)
+ {
+ out << p.row() << " ";
+ out << p.col() << " ";
+ out << img(p).red() << " ";
+ out << img(p).green() << " ";
+ out << img(p).blue() << std::endl;
+ }
+
+ // Close gnuplot data stream.
+ out << "e" << std::endl;
+ out << "EOF" << std::endl;
+
+ out.close();
}
-
- // Close gnuplot data stream.
- out << "e" << std::endl;
- out << "EOF" << std::endl;
-
- out.close();
- }
- else
- {
- std::cerr << "ERROR[mln::io::plot::save_image_sh]:" << filename;
- std::cerr << " couldn't be opened !!" << std::endl;
- }
-
- trace::exiting("mln::io::plot::save_image_sh<image2d<float>>");
- return result;
- }
-
- //------------------------------------------------------------------------
- // save_image_sh(const image2d<double>&, const std::string&)
- //------------------------------------------------------------------------
-
- inline
- bool save_image_sh_dispatch_overload(const image2d<double>& img,
- const std::string& filename)
- {
- trace::entering("mln::io::plot::save_image_sh<image2d<double>>");
- mln_precondition(img.is_valid());
-
- std::ofstream out(filename.c_str());
- bool result = !out.fail();
- unsigned min_row = geom::min_row(img);
- unsigned max_row = geom::max_row(img);
- unsigned min_col = geom::min_col(img);
- unsigned max_col = geom::max_col(img);
-
- if (result)
- {
- mln_piter_(image2d< double >) p(img.domain());
-
- // Output data prelude (terminal X11, image).
- out << "#!/bin/sh" << std::endl;
- out << "####################################" << std::endl;
- out << "# Columns = (x, y, double) #" << std::endl;
- out << "####################################" << std::endl;
- out << std::endl;
- out << "gnuplot <<EOF" << std::endl;
- out << "set terminal x11 persist 1" << std::endl;
- out << "set palette gray" << std::endl;
- out << "set xrange [" << min_col;
- out << ":" << max_col;
- out << "]" << std::endl;
- out << "set yrange [-" << max_row;
- out << ":" << min_row;
- out << "]" << std::endl;
- out << "plot '-' using (\\$2):(-\\$1):(\\$3) with image" << std::endl;
-
- // Output data.
- for_all(p)
+ else
{
- out << p.row() << " ";
- out << p.col() << " ";
- out << img(p) << std::endl;
+ std::cerr << "ERROR[mln::io::plot::save_image_sh]:" << filename;
+ std::cerr << " couldn't be opened !!" << std::endl;
}
-
- // Close gnuplot data stream.
- out << "e" << std::endl;
- out << "EOF" << std::endl;
- out.close();
- }
- else
- {
- std::cerr << "ERROR[mln::io::plot::save_image_sh]:" << filename;
- std::cerr << " couldn't be opened !!" << std::endl;
+ trace::exiting("mln::io::plot::impl::save_image_sh_image2d_rgb");
+ return result;
}
- trace::exiting("mln::io::plot::save_image_sh<image2d<double>>");
- return result;
- }
+ //----------------------------------------------------------------------
+ // save_image_sh_image2d_hsl(const image2d<hsl_<T,T,T>>&, const string&)
+ //----------------------------------------------------------------------
- //------------------------------------------------------------------------
- // save_image_sh(const image2d<int_u<n>>&, const std::string&)
- //------------------------------------------------------------------------
+ template <typename T>
+ inline
+ bool save_image_sh_image2d_hsl(const image2d<value::hsl_<T,T,T> >& img,
+ const std::string& filename)
+ {
+ trace::entering("mln::io::plot::impl::save_image_sh_image2d_hsl");
+ mln_precondition(img.is_valid());
- template <unsigned n>
- inline
- bool save_image_sh_dispatch_overload(const image2d<value::int_u<n> >& img,
- const std::string& filename)
- {
- trace::entering("mln::io::plot::save_image_sh<image2d<int_u<n>>");
- mln_precondition(img.is_valid());
+ std::ofstream out(filename.c_str());
+ bool result = !out.fail();
+ unsigned min_row = geom::min_row(img);
+ unsigned max_row = geom::max_row(img);
+ unsigned min_col = geom::min_col(img);
+ unsigned max_col = geom::max_col(img);
- std::ofstream out(filename.c_str());
- bool result = !out.fail();
- unsigned min_row = geom::min_row(img);
- unsigned max_row = geom::max_row(img);
- unsigned min_col = geom::min_col(img);
- unsigned max_col = geom::max_col(img);
+ typedef mln::value::hsl_<T,T,T> hsl;
- if (result)
- {
- mln_piter(image2d< value::int_u<n> >) p(img.domain());
-
- // Output data prelude (terminal X11, image).
- out << "#!/bin/sh" << std::endl;
- out << "####################################" << std::endl;
- out << "# Columns = (x, y, int_u<n>) #" << std::endl;
- out << "####################################" << std::endl;
- out << std::endl;
- out << "gnuplot <<EOF" << std::endl;
- out << "set terminal x11 persist 1" << std::endl;
- out << "set palette gray" << std::endl;
- out << "set xrange [" << min_col;
- out << ":" << max_col;
- out << "]" << std::endl;
- out << "set yrange [-" << max_row;
- out << ":" << min_row;
- out << "]" << std::endl;
- out << "plot '-' using (\\$2):(-\\$1):(\\$3) with image" << std::endl;
-
- // Output data.
- for_all(p)
+ if (result)
+ {
+ mln_piter(image2d< hsl >) p(img.domain());
+
+ // Output data prelude (terminal X11, image).
+ out << "#!/bin/sh" << std::endl;
+ out << "####################################" << std::endl;
+ out << "# Columns = (x, y, val) #" << std::endl;
+ out << "####################################" << std::endl;
+ out << std::endl;
+ out << "gnuplot <<EOF" << std::endl;
+ out << "q(l,s) = (l < 0.5)? (l*(1+s)):(l+s-(l*s))" << std::endl;
+ out << "p(l,s) = (2.0 * l - q(l,s))" << std::endl;
+ out << "n(x) = (x < 0)?(x+1):(x > 1)?(x-1) : (x)"<< std::endl;
+ out << std::endl;
+ out << "c(p,q,t) = (t <(1.0/6.0))?(p+(q-p)*6.0*t):\\"<< std::endl;
+ out << " (t <(1.0/2.0))?(q) :\\" << std::endl;
+ out << " (t<(2.0/3.0))?(p+(q-p)*6.0*((2.0/3.0)-t)):\\" << std::endl;
+ out << " (p)" << std::endl;
+ out << std::endl;
+ out << "r(h,s,l) = c(p(l,s),q(l,s),n(h/360.0+1.0/3.0))"<< std::endl;
+ out << "g(h,s,l) = c(p(l,s),q(l,s),n(h/360.0))" << std::endl;
+ out << "b(h,s,l) = c(p(l,s),q(l,s),n(h/360.0-1.0/3.0))"<< std::endl;
+ out << std::endl;
+ out << "set terminal x11 persist 1" << std::endl;
+ out << "set palette gray" << std::endl;
+ out << "set xrange [" << min_col;
+ out << ":" << max_col;
+ out << "]" << std::endl;
+ out << "set yrange [-" << max_row;
+ out << ":" << min_row;
+ out << "]" << std::endl;
+ out << "plot '-' using 2:(-\\$1):\\" << std::endl;
+ out << " (r(\\$3,\\$4,\\$5)):\\" << std::endl;
+ out << " (g(\\$3,\\$4,\\$5)):\\" << std::endl;
+ out << " (b(\\$3,\\$4,\\$5)) with rgbimage" << std::endl;
+
+ // Output data.
+ for_all(p)
+ {
+ out << p.row() << " ";
+ out << p.col() << " ";
+ out << img(p).hue() << " ";
+ out << img(p).sat() << " ";
+ out << img(p).lum() << std::endl;
+ }
+
+ // Close gnuplot data stream.
+ out << "e" << std::endl;
+ out << "EOF" << std::endl;
+
+ out.close();
+ }
+ else
{
- out << p.row() << " ";
- out << p.col() << " ";
- out << img(p) << std::endl;
+ std::cerr << "ERROR[mln::io::plot::save_image_sh]:" << filename;
+ std::cerr << " couldn't be opened !!" << std::endl;
}
-
- // Close gnuplot data stream.
- out << "e" << std::endl;
- out << "EOF" << std::endl;
- out.close();
- }
- else
- {
- std::cerr << "ERROR[mln::io::plot::save_image_sh]:" << filename;
- std::cerr << " couldn't be opened !!" << std::endl;
+ trace::exiting("mln::io::plot::impl::save_image_sh_image2d_hsl");
+ return result;
}
- trace::exiting("mln::io::plot::save_image_sh<image2d<int_u<n>>");
- return result;
- }
+ //----------------------------------------------------------------------
+ // save_image_sh_image2d_hsi(const image2d<hsi_<T,T,T>>&, const string&)
+ //----------------------------------------------------------------------
- //------------------------------------------------------------------------
- // save_image_sh(const image2d<int_s<n>>&, const std::string&)
- //------------------------------------------------------------------------
+ template <typename T>
+ inline
+ bool save_image_sh_image2d_hsi(const image2d<value::hsi_<T,T,T> >& img,
+ const std::string& filename)
+ {
+ trace::entering("mln::io::plot::impl::save_image_sh_image2d_hsi");
+ mln_precondition(img.is_valid());
- template <unsigned n>
- inline
- bool save_image_sh_dispatch_overload(const image2d<value::int_s<n> >& img,
- const std::string& filename)
- {
- trace::entering("mln::io::plot::save_image_sh<image2d<int_s<n>>");
- mln_precondition(img.is_valid());
+ std::ofstream out(filename.c_str());
+ bool result = !out.fail();
+ unsigned min_row = geom::min_row(img);
+ unsigned max_row = geom::max_row(img);
+ unsigned min_col = geom::min_col(img);
+ unsigned max_col = geom::max_col(img);
- std::ofstream out(filename.c_str());
- bool result = !out.fail();
- unsigned min_row = geom::min_row(img);
- unsigned max_row = geom::max_row(img);
- unsigned min_col = geom::min_col(img);
- unsigned max_col = geom::max_col(img);
+ typedef mln::value::hsi_<T,T,T> hsi;
- if (result)
- {
- mln_piter(image2d< value::int_s<n> >) p(img.domain());
-
- // Output data prelude (terminal X11, image).
- out << "#!/bin/sh" << std::endl;
- out << "####################################" << std::endl;
- out << "# Columns = (x, y, int_s<n>) #" << std::endl;
- out << "####################################" << std::endl;
- out << std::endl;
- out << "gnuplot <<EOF" << std::endl;
- out << "set terminal x11 persist 1" << std::endl;
- out << "set palette gray" << std::endl;
- out << "set xrange [" << min_col;
- out << ":" << max_col;
- out << "]" << std::endl;
- out << "set yrange [-" << max_row;
- out << ":" << min_row;
- out << "]" << std::endl;
- out << "plot '-' using 2:(-\\$1):3 with image" << std::endl;
-
- // Output data.
- for_all(p)
+ if (result)
+ {
+ mln_piter(image2d< hsi >) p(img.domain());
+
+ // Output data prelude (terminal X11, image).
+ out << "#!/bin/sh" << std::endl;
+ out << "####################################" << std::endl;
+ out << "# Columns = (x, y, val) #" << std::endl;
+ out << "####################################" << std::endl;
+ out << std::endl;
+ out << "gnuplot <<EOF" << std::endl;
+ out << "teta(h) = (h/180.0) * pi" << std::endl;
+ out << "alpha(s,h) = s*cos(teta(h))" << std::endl;
+ out << "beta(s,h) = s*sin(teta(h))" << std::endl;
+ out << "n(x) = (x < 0)?(x+1) : (x > 1)?(x-1):(x)"<< std::endl;
+ out << std::endl;
+ out << "c(p,q,t) = (t <(1.0/6.0))?(p+(q-p)*6.0*t):\\"<< std::endl;
+ out << " (t <(1.0/2.0))?(q) :\\" << std::endl;
+ out << " (t <(2.0/3.0))?(p+(q-p)*6.0*((2.0/3.0)-t)):\\"<< std::endl;
+ out << " (p)" << std::endl;
+ out << std::endl;
+ out << "r(h,s,i) = (sqrt(3.0)/3.0) * i \\" << std::endl;
+ out << " + (2.0/(sqrt(6.0))) * beta(s,h)" << std::endl;
+ out << "g(h,s,i) = (sqrt(3.0)/3.0) * i \\" << std::endl;
+ out << " + (2.0/(sqrt(2.0))) * alpha(s,h)\\" << std::endl;
+ out << " - (1.0/(sqrt(6.0))) * beta(s,h)" << std::endl;
+ out << "b(h,s,i) = (sqrt(3.0)/3.0) * i \\" << std::endl;
+ out << " - (2.0/(sqrt(2.0))) * alpha(s,h) \\"<< std::endl;
+ out << " - (1.0/(sqrt(6.0))) * beta(s,h)" << std::endl;
+ out << std::endl;
+ out << "set terminal x11 persist 1" << std::endl;
+ out << "set palette gray" << std::endl;
+ out << "set xrange [" << min_col;
+ out << ":" << max_col;
+ out << "]" << std::endl;
+ out << "set yrange [-" << max_row;
+ out << ":" << min_row;
+ out << "]" << std::endl;
+ out << "plot '-' using 2:(-\\$1):\\" << std::endl;
+ out << " (r(\\$3,\\$4,\\$5)):\\" << std::endl;
+ out << " (g(\\$3,\\$4,\\$5)):\\" << std::endl;
+ out << " (b(\\$3,\\$4,\\$5)) with rgbimage" << std::endl;
+
+ // Output data.
+ for_all(p)
+ {
+ out << p.row() << " ";
+ out << p.col() << " ";
+ out << img(p).hue() << " ";
+ out << img(p).sat() << " ";
+ out << img(p).inty() << std::endl;
+ }
+
+ // Close gnuplot data stream.
+ out << "e" << std::endl;
+ out << "EOF" << std::endl;
+
+ out.close();
+ }
+ else
{
- out << p.row() << " ";
- out << p.col() << " ";
- out << img(p) << std::endl;
+ std::cerr << "ERROR[mln::io::plot::save_image_sh]:" << filename;
+ std::cerr << " couldn't be opened !!" << std::endl;
}
-
- // Close gnuplot data stream.
- out << "e" << std::endl;
- out << "EOF" << std::endl;
- out.close();
- }
- else
- {
- std::cerr << "ERROR[mln::io::plot::save_image_sh]:" << filename;
- std::cerr << " couldn't be opened !!" << std::endl;
+ trace::exiting("mln::io::plot::impl::save_image_sh_image2d_hsi");
+ return result;
}
- trace::exiting("mln::io::plot::save_image_sh<image2d<int_s<n>>");
- return result;
- }
-
- //------------------------------------------------------------------------
- // save_image_sh(const image2d<rgb<n>>&, const std::string&)
- //------------------------------------------------------------------------
+ //----------------------------------------------------------------------
+ // save_image_sh_image2d(const image2d<I>&, const string&)
+ //----------------------------------------------------------------------
- template <unsigned n>
- inline
- bool save_image_sh_dispatch_overload(const image2d<value::rgb<n> >& img,
- const std::string& filename)
- {
- trace::entering("mln::io::plot::save_image_sh<image2d<rgb<n>>");
- mln_precondition(img.is_valid());
+ template <typename I>
+ inline
+ bool save_image_sh_image2d(const image2d<I>& img,
+ const std::string& filename)
+ {
+ trace::entering("mln::io::plot::impl::save_image_sh_image2d");
+ mln_precondition(img.is_valid());
- std::ofstream out(filename.c_str());
- bool result = !out.fail();
- unsigned min_row = geom::min_row(img);
- unsigned max_row = geom::max_row(img);
- unsigned min_col = geom::min_col(img);
- unsigned max_col = geom::max_col(img);
+ std::ofstream out(filename.c_str());
+ bool result = !out.fail();
+ unsigned min_row = geom::min_row(img);
+ unsigned max_row = geom::max_row(img);
+ unsigned min_col = geom::min_col(img);
+ unsigned max_col = geom::max_col(img);
- if (result)
- {
- mln_piter(image2d< value::int_u<n> >) p(img.domain());
-
- // Output data prelude (terminal X11, image).
- out << "#!/bin/sh" << std::endl;
- out << "####################################" << std::endl;
- out << "# Columns = (x, y, r, g, b) #" << std::endl;
- out << "####################################" << std::endl;
- out << std::endl;
- out << "gnuplot <<EOF" << std::endl;
- out << "set terminal x11 persist 1" << std::endl;
- out << "set xrange [" << min_col;
- out << ":" << max_col;
- out << "]" << std::endl;
- out << "set yrange [-" << max_row;
- out << ":" << min_row;
- out << "]" << std::endl;
- out << "plot '-' using 2:(-\\$1):3:4:5 with rgbimage" << std::endl;
-
- // Output data.
- for_all(p)
+ if (result)
{
- out << p.row() << " ";
- out << p.col() << " ";
- out << img(p).red() << " ";
- out << img(p).green() << " ";
- out << img(p).blue() << std::endl;
+ mln_piter(image2d<I>) p(img.domain());
+
+ // Output data prelude (terminal X11, image).
+ out << "#!/bin/sh" << std::endl;
+ out << "####################################" << std::endl;
+ out << "# Columns = (x, y, val) #" << std::endl;
+ out << "####################################" << std::endl;
+ out << std::endl;
+ out << "gnuplot <<EOF" << std::endl;
+ out << "set terminal x11 persist 1" << std::endl;
+ out << "set palette gray" << std::endl;
+ out << "set xrange [" << min_col;
+ out << ":" << max_col;
+ out << "]" << std::endl;
+ out << "set yrange [-" << max_row;
+ out << ":" << min_row;
+ out << "]" << std::endl;
+ out << "plot '-' using 2:(-\\$1):3 with image" <<std::endl;
+
+ // Output data.
+ for_all(p)
+ {
+ out << p.row() << " ";
+ out << p.col() << " ";
+ out << img(p) << std::endl;
+ }
+
+ // Close gnuplot data stream.
+ out << "e" << std::endl;
+ out << "EOF" << std::endl;
+
+ out.close();
+ }
+ else
+ {
+ std::cerr << "ERROR[mln::io::plot::save_image_sh]:"<< filename;
+ std::cerr << " couldn't be opened !!" << std::endl;
}
-
- // Close gnuplot data stream.
- out << "e" << std::endl;
- out << "EOF" << std::endl;
- out.close();
+ trace::exiting("mln::io::plot::impl::save_image_sh_image2d");
+ return result;
}
- else
- {
- std::cerr << "ERROR[mln::io::plot::save_image_sh]:" << filename;
- std::cerr << " couldn't be opened !!" << std::endl;
- }
-
- trace::exiting("mln::io::plot::save_image_sh<image2d<rgb<n>>");
- return result;
- }
-
- //------------------------------------------------------------------------
- // save_image_sh(const image2d<hsl_<T,T,T>>&, const std::string&)
- //------------------------------------------------------------------------
- template <typename T>
- inline
- bool save_image_sh(const image2d<mln::value::hsl_<T,T,T> >& img,
- const std::string& filename)
- {
- trace::entering("mln::io::plot::save_image_sh<image2d<hsl>>");
- mln_precondition(img.is_valid());
+ //----------------------------------------------------------------------
+ // save_image_sh_image1d(const image1d<I>&, const string&)
+ //----------------------------------------------------------------------
- std::ofstream out(filename.c_str());
- bool result = !out.fail();
- unsigned min_row = geom::min_row(img);
- unsigned max_row = geom::max_row(img);
- unsigned min_col = geom::min_col(img);
- unsigned max_col = geom::max_col(img);
+ template <typename I>
+ inline
+ bool save_image_sh_image1d(const image1d<I>& img,
+ const std::string& filename)
+ {
+ trace::entering("mln::io::plot::impl::save_image_sh_image1d");
+ mln_precondition(img.is_valid());
- typedef mln::value::hsl_<T,T,T> hsl;
+ std::ofstream out(filename.c_str());
+ bool result = !out.fail();
- if (result)
- {
- mln_piter(image2d< hsl >) p(img.domain());
-
- // Output data prelude (terminal X11, image).
- out << "#!/bin/sh" << std::endl;
- out << "####################################" << std::endl;
- out << "# Columns = (x, y, val) #" << std::endl;
- out << "####################################" << std::endl;
- out << std::endl;
- out << "gnuplot <<EOF" << std::endl;
- out << "q(l,s) = (l < 0.5)? (l*(1+s)) : (l+s-(l*s))" << std::endl;
- out << "p(l,s) = (2.0 * l - q(l,s))" << std::endl;
- out << "n(x) = (x < 0)?(x+1) : (x > 1)?(x-1) : (x)"<< std::endl;
- out << std::endl;
- out << "c(p,q,t) = (t < (1.0/6.0))?(p+(q-p)*6.0*t):\\" << std::endl;
- out << " (t < (1.0/2.0))?(q) :\\" << std::endl;
- out << " (t < (2.0/3.0))?(p+(q-p)*6.0*((2.0/3.0)-t)):\\" << std::endl;
- out << " (p)" << std::endl;
- out << std::endl;
- out << "r(h,s,l) = c(p(l,s),q(l,s),n(h/360.0+1.0/3.0))"<< std::endl;
- out << "g(h,s,l) = c(p(l,s),q(l,s),n(h/360.0))" << std::endl;
- out << "b(h,s,l) = c(p(l,s),q(l,s),n(h/360.0-1.0/3.0))"<< std::endl;
- out << std::endl;
- out << "set terminal x11 persist 1" << std::endl;
- out << "set palette gray" << std::endl;
- out << "set xrange [" << min_col;
- out << ":" << max_col;
- out << "]" << std::endl;
- out << "set yrange [-" << max_row;
- out << ":" << min_row;
- out << "]" << std::endl;
- out << "plot '-' using 2:(-\\$1):(r(\\$3,\\$4,\\$5)):\\" << std::endl;
- out << " (g(\\$3,\\$4,\\$5)):\\" << std::endl;
- out << " (b(\\$3,\\$4,\\$5)) with rgbimage"<< std::endl;
-
- // Output data.
- for_all(p)
+ if (result)
+ {
+ mln_piter(image1d<I>) p(img.domain());
+
+ // Output data prelude (terminal X11, impulse).
+ out << "#!/bin/sh" << std::endl;
+ out << "##########################" << std::endl;
+ out << "# Two columns = (x, val) #" << std::endl;
+ out << "##########################" << std::endl;
+ out << std::endl;
+ out << "gnuplot <<EOF" << std::endl;
+ out << "set terminal x11 persist 1" << std::endl;
+ out << "plot '-' with impulse" << std::endl;
+
+ // Output data.
+ for_all(p)
+ {
+ out << p.ind() << " ";
+ out << img(p) << std::endl;
+ }
+
+ // Close gnuplot data stream.
+ out << "e" << std::endl;
+ out << "EOF" << std::endl;
+
+ out.close();
+ }
+ else
{
- out << p.row() << " ";
- out << p.col() << " ";
- out << img(p).hue() << " ";
- out << img(p).sat() << " ";
- out << img(p).lum() << std::endl;
+ std::cerr << "ERROR[mln::io::plot::save_image_sh]:"<< filename;
+ std::cerr << " couldn't be opened !!" << std::endl;
}
-
- // Close gnuplot data stream.
- out << "e" << std::endl;
- out << "EOF" << std::endl;
- out.close();
- }
- else
- {
- std::cerr << "ERROR[mln::io::plot::save_image_sh]:" << filename;
- std::cerr << " couldn't be opened !!" << std::endl;
+ trace::exiting("mln::io::plot::impl::save_image_sh_image1d");
+ return result;
}
- trace::exiting("mln::io::plot::save_image_sh<image2d<hsl>>");
- return result;
- }
+ //----------------------------------------------------------------------
+ // save_image_sh_image3d(const image3d<I>&, const string&)
+ //----------------------------------------------------------------------
- //------------------------------------------------------------------------
- // save_image_sh(const image2d<hsi_<T,T,T>>&, const std::string&)
- //------------------------------------------------------------------------
-
- template <typename T>
- inline
- bool save_image_sh(const image2d<mln::value::hsi_<T,T,T> >& img,
- const std::string& filename)
- {
- trace::entering("mln::io::plot::save_image_sh<2d>vf");
- mln_precondition(img.is_valid());
- std::ofstream out(filename.c_str());
- bool result = !out.fail();
- unsigned min_row = geom::min_row(img);
- unsigned max_row = geom::max_row(img);
- unsigned min_col = geom::min_col(img);
- unsigned max_col = geom::max_col(img);
+ template <typename I>
+ inline
+ bool save_image_sh_image3d(const image3d<I>& img,
+ const std::string& filename)
+ {
+ trace::entering("mln::io::plot::impl::save_image_sh_image3d");
+ mln_precondition(img.is_valid());
- typedef mln::value::hsi_<T,T,T> hsi;
+ std::ofstream out(filename.c_str());
+ bool result = !out.fail();
- if (result)
- {
- mln_piter(image2d< hsi >) p(img.domain());
-
- // Output data prelude (terminal X11, image).
- out << "#!/bin/sh" << std::endl;
- out << "####################################" << std::endl;
- out << "# Columns = (x, y, val) #" << std::endl;
- out << "####################################" << std::endl;
- out << std::endl;
- out << "gnuplot <<EOF" << std::endl;
- out << "teta(h) = (h/180.0) * pi" << std::endl;
- out << "alpha(s,h) = s*cos(teta(h))" << std::endl;
- out << "beta(s,h) = s*sin(teta(h))" << std::endl;
- out << "n(x) = (x < 0)?(x+1) : (x > 1)?(x-1) : (x)"<< std::endl;
- out << std::endl;
- out << "c(p,q,t) = (t < (1.0/6.0))?(p+(q-p)*6.0*t):\\" << std::endl;
- out << " (t < (1.0/2.0))?(q) :\\" << std::endl;
- out << " (t < (2.0/3.0))?(p+(q-p)*6.0*((2.0/3.0)-t)):\\" << std::endl;
- out << " (p)" << std::endl;
- out << std::endl;
- out << "r(h,s,i) = (sqrt(3.0)/3.0) * i \\" << std::endl;
- out << " + (2.0/(sqrt(6.0))) * beta(s,h)" << std::endl;
- out << "g(h,s,i) = (sqrt(3.0)/3.0) * i \\" << std::endl;
- out << " + (2.0/(sqrt(2.0))) * alpha(s,h) \\"<< std::endl;
- out << " - (1.0/(sqrt(6.0))) * beta(s,h)" << std::endl;
- out << "b(h,s,i) = (sqrt(3.0)/3.0) * i \\" << std::endl;
- out << " - (2.0/(sqrt(2.0))) * alpha(s,h) \\"<< std::endl;
- out << " - (1.0/(sqrt(6.0))) * beta(s,h)" << std::endl;
- out << std::endl;
- out << "set terminal x11 persist 1" << std::endl;
- out << "set palette gray" << std::endl;
- out << "set xrange [" << min_col;
- out << ":" << max_col;
- out << "]" << std::endl;
- out << "set yrange [-" << max_row;
- out << ":" << min_row;
- out << "]" << std::endl;
- out << "plot '-' using 2:(-\\$1):(r(\\$3,\\$4,\\$5)):\\" << std::endl;
- out << " (g(\\$3,\\$4,\\$5)):\\" << std::endl;
- out << " (b(\\$3,\\$4,\\$5)) with rgbimage"<< std::endl;
-
- // Output data.
- for_all(p)
+ if (result)
{
- out << p.row() << " ";
- out << p.col() << " ";
- out << img(p).hue() << " ";
- out << img(p).sat() << " ";
- out << img(p).inty() << std::endl;
+ mln_piter(image3d<I>) p(img.domain());
+
+ // Output data prelude (terminal X11, pointtype 0).
+ out << "#!/bin/sh" << std::endl;
+ out << "####################################" << std::endl;
+ out << "# Columns = (x, y, z, val) #" << std::endl;
+ out << "####################################" << std::endl;
+ out << std::endl;
+ out << "gnuplot <<EOF" << std::endl;
+ out << "set terminal x11 persist 1" << std::endl;
+ out << "splot '-' with points palette pointtype 7" << std::endl;
+
+ // Output data.
+ for_all(p)
+ {
+ out << p.row() << " ";
+ out << p.col() << " ";
+ out << p.sli() << " ";
+ out << img(p) << std::endl;
+ }
+
+ // Close gnuplot data stream.
+ out << "e" << std::endl;
+ out << "EOF" << std::endl;
+
+ out.close();
+ }
+ else
+ {
+ std::cerr << "ERROR[mln::io::plot::save_image_sh]:"<< filename;
+ std::cerr << " couldn't be opened !!" << std::endl;
}
-
- // Close gnuplot data stream.
- out << "e" << std::endl;
- out << "EOF" << std::endl;
- out.close();
- }
- else
- {
- std::cerr << "ERROR[mln::io::plot::save_image_sh]:" << filename;
- std::cerr << " couldn't be opened !!" << std::endl;
+ trace::exiting("mln::io::plot::impl::save_image_sh_image3d");
+ return result;
}
- trace::exiting("mln::io::plot::save_image_sh<2d>vf");
- return result;
- }
+ } // end of namespace impl
+
//------------------------------------------------------------------------
- // save_image_sh(const image2d<I>&, const std::string&)
+ // Internal.
//------------------------------------------------------------------------
- template <typename I>
- inline
- bool save_image_sh_dispatch_overload(const image2d<I>& img,
- const std::string& filename)
+ namespace internal
{
- trace::entering("mln::io::plot::save_image_sh<image2d<I>>");
- mln_precondition(img.is_valid());
-
- std::ofstream out(filename.c_str());
- bool result = !out.fail();
-
- if (result)
+ template <unsigned n>
+ inline
+ bool save_image_sh_dispatch(const image2d<value::rgb<n> >& img,
+ const std::string& filename)
{
- mln_piter(image2d<I>) p(img.domain());
-
- // Output data prelude (terminal X11, image).
- out << "#!/bin/sh" << std::endl;
- out << "####################################" << std::endl;
- out << "# Columns = (x, y, builtin) #" << std::endl;
- out << "####################################" << std::endl;
- out << std::endl;
- out << "gnuplot <<EOF" << std::endl;
- out << "set terminal x11 persist 1" << std::endl;
- out << "set pm3d" << std::endl;
- out << "set dgrid3d 20,20,2" << std::endl;
- out << "splot '-' with lines" << std::endl;
-
- // Output data.
- for_all(p)
- {
- out << p.col() << " ";
- out << p.row() << " ";
- out << img(p) << std::endl;
- }
-
- // Close gnuplot data stream.
- out << "e" << std::endl;
- out << "EOF" << std::endl;
-
- out.close();
+ return impl::save_image_sh_image2d_rgb(img, filename);
}
- else
+
+ template <typename T>
+ inline
+ bool save_image_sh_dispatch(const image2d<value::hsl_<T,T,T> >& img,
+ const std::string& filename)
{
- std::cerr << "ERROR[mln::io::plot::save_image_sh]:"<< filename;
- std::cerr << " couldn't be opened !!" << std::endl;
+ return impl::save_image_sh_image2d_hsl(img, filename);
}
-
- trace::exiting("mln::io::plot::save_image_sh<image2d<I>>");
- return result;
- }
- //------------------------------------------------------------------------
- // save_image_sh(const image1d<I>&, const std::string&)
- //------------------------------------------------------------------------
+ template <typename T>
+ inline
+ bool save_image_sh_dispatch(const image2d<value::hsi_<T,T,T> >& img,
+ const std::string& filename)
+ {
+ return impl::save_image_sh_image2d_hsi(img, filename);
+ }
- template <typename I>
- inline
- bool save_image_sh_dispatch_overload(const image1d<I>& img,
- const std::string& filename)
- {
- trace::entering("mln::io::plot::save_image_sh<image1d<I>>");
- mln_precondition(img.is_valid());
- std::ofstream out(filename.c_str());
- bool result = !out.fail();
-
- if (result)
+ template <typename I>
+ inline
+ bool save_image_sh_dispatch(const image2d<I>& img,
+ const std::string& filename)
{
- mln_piter(image1d<I>) p(img.domain());
-
- // Output data prelude (terminal X11, impulse).
- out << "#!/bin/sh" << std::endl;
- out << "##########################" << std::endl;
- out << "# Two columns = (x, val) #" << std::endl;
- out << "##########################" << std::endl;
- out << std::endl;
- out << "gnuplot <<EOF" << std::endl;
- out << "set terminal x11 persist 1" << std::endl;
- out << "plot '-' with impulse" << std::endl;
-
- // Output data.
- for_all(p)
- {
- out << p.ind() << " ";
- out << img(p) << std::endl;
- }
-
- // Close gnuplot data stream.
- out << "e" << std::endl;
- out << "EOF" << std::endl;
-
- out.close();
+ return impl::save_image_sh_image2d(img, filename);
}
- else
+
+ template <typename I>
+ inline
+ bool save_image_sh_dispatch(const image1d<I>& img,
+ const std::string& filename)
{
- std::cerr << "ERROR[mln::io::plot::save_image_sh]:"<< filename;
- std::cerr << " couldn't be opened !!" << std::endl;
+ return impl::save_image_sh_image1d(img, filename);
}
-
- trace::exiting("mln::io::plot::save_image_sh<image1d<I>>");
- return result;
- }
-
- //------------------------------------------------------------------------
- // save_image_sh(const image3d<I>&, const std::string&)
- //------------------------------------------------------------------------
-
- template <typename I>
- inline
- bool save_image_sh(const image3d<I>& img, const std::string& filename)
- {
- trace::entering("mln::io::plot::save_image_sh<image3d<I>>");
- mln_precondition(img.is_valid());
-
- std::ofstream out(filename.c_str());
- bool result = !out.fail();
-
- if (result)
+ template <typename I>
+ inline
+ bool save_image_sh_dispatch(const image3d<I>& img,
+ const std::string& filename)
{
- mln_piter(image3d<I>) p(img.domain());
-
- // Output data prelude (terminal X11, pointtype 0).
- out << "#!/bin/sh" << std::endl;
- out << "####################################" << std::endl;
- out << "# Columns = (x, y, z, val) #" << std::endl;
- out << "####################################" << std::endl;
- out << std::endl;
- out << "gnuplot <<EOF" << std::endl;
- out << "set terminal x11 persist 1" << std::endl;
- out << "splot '-' with points palette pointtype 7" << std::endl;
-
- // Output data.
- for_all(p)
- {
- out << p.row() << " ";
- out << p.col() << " ";
- out << p.sli() << " ";
- out << img(p) << std::endl;
- }
-
- // Close gnuplot data stream.
- out << "e" << std::endl;
- out << "EOF" << std::endl;
-
- out.close();
+ return impl::save_image_sh_image3d(img, filename);
}
- else
+
+ template <typename I>
+ inline
+ bool save_image_sh_dispatch(const Image<I>& img,
+ const std::string& filename)
{
- std::cerr << "ERROR[mln::io::plot::save_image_sh]:"<< filename;
- std::cerr << " couldn't be opened !!" << std::endl;
+ return save_image_sh_dispatch(exact(img), filename);
}
-
- trace::exiting("mln::io::plot::save_image_sh<image3d<I>>");
- return result;
- }
+
+ } // end of namespace mln::io::plot::internal
//------------------------------------------------------------------------
// Facade.
@@ -825,19 +598,11 @@ namespace mln
template <typename I>
inline
- bool save_image_sh_dispatch(const Image<I>& img,
- const std::string& filename)
- {
- return save_image_sh_dispatch_overload(exact(img), filename);
- }
-
- template <typename I>
- inline
bool save_image_sh(const Image<I>& img, const std::string& filename)
{
trace::entering("mln::io::plot::save_image_sh");
- bool result = save_image_sh_dispatch(img, filename);
+ bool result = internal::save_image_sh_dispatch(img, filename);
trace::exiting("mln::io::plot::save_image_sh");
return result;
@@ -847,9 +612,9 @@ namespace mln
#endif // ! MLN_INCLUDE_ONLY
- } // end of namespace mln::accu::stat
+ } // end of namespace mln::io::plot
- } // end of namespace mln::accu
+ } // end of namespace mln::io
} // end of namespace mln
--
1.5.6.5
