#include "MagickCore/image-private.h"
#include "MagickCore/list.h"
#include "MagickCore/log.h"
+#include "MagickCore/matrix.h"
#include "MagickCore/memory_.h"
#include "MagickCore/memory-private.h"
#include "MagickCore/monitor.h"
% The format of the EdgeImage method is:
%
% Image *CannyEdgeImage(const Image *image,const double radius,
-% const double sigma,const double low_factor,const double high_factor,
-% const size_t threshold,ExceptionInfo *exception)
+% const double sigma,const double lower_precent,
+% const double upper_percent,ExceptionInfo *exception)
%
% A description of each parameter follows:
%
%
% o sigma: the sigma of the gaussian smoothing filter.
%
-% o low_factor: use this low factor in hysteresis.
+% o lower_precent: percentage of edge pixels in the lower threshold.
%
-% o hight_factor: use this high factor in hysteresis.
+% o upper_percent: percentage of edge pixels in the upper threshold.
%
% o exception: return any errors or warnings in this structure.
%
*/
+
+typedef struct _CannyInfo
+{
+ double
+ magnitude,
+ intensity;
+
+ int
+ orientation;
+
+ ssize_t
+ x,
+ y;
+} CannyInfo;
+
+static inline MagickBooleanType IsAuthenticPixel(const Image *image,
+ const ssize_t x,const ssize_t y)
+{
+ if ((x < 0) || (x >= (ssize_t) image->columns))
+ return(MagickFalse);
+ if ((y < 0) || (y >= (ssize_t) image->rows))
+ return(MagickFalse);
+ return(MagickTrue);
+}
+
+static MagickBooleanType TraceEdges(Image *edge_image,CacheView *trace_view,
+ MatrixInfo *pixel_cache,const ssize_t x,const ssize_t y,
+ const double lower_threshold,ExceptionInfo *exception)
+{
+ CannyInfo
+ pixel;
+
+ size_t
+ number_edges;
+
+ if (GetMatrixElement(pixel_cache,0,0,&pixel) == MagickFalse)
+ return(MagickFalse);
+ pixel.x=x;
+ pixel.y=y;
+ if (SetMatrixElement(pixel_cache,0,0,&pixel) == MagickFalse)
+ return(MagickFalse);
+ number_edges=1;
+ do
+ {
+ MagickBooleanType
+ status;
+
+ ssize_t
+ v,
+ x_offset,
+ y_offset;
+
+ number_edges--;
+ status=GetMatrixElement(pixel_cache,(ssize_t) number_edges,0,&pixel);
+ if (status == MagickFalse)
+ return(MagickFalse);
+ x_offset=pixel.x;
+ y_offset=pixel.y;
+ for (v=(-1); v <= 1; v++)
+ {
+ ssize_t
+ u;
+
+ for (u=(-1); u <= 1; u++)
+ {
+ Quantum
+ *q;
+
+ if ((u == 0) && (v == 0))
+ continue;
+ if (IsAuthenticPixel(edge_image,x_offset+u,y_offset+v) == MagickFalse)
+ continue;
+ /*
+ Not an edge if gradient value is below the lower threshold.
+ */
+ q=GetCacheViewAuthenticPixels(trace_view,x_offset+u,y_offset+v,1,1,
+ exception);
+ if (q == (Quantum *) NULL)
+ return(MagickFalse);
+ status=GetMatrixElement(pixel_cache,x_offset+u,y_offset+v,&pixel);
+ if (status == MagickFalse)
+ return(MagickFalse);
+ if ((pixel.intensity >= lower_threshold) &&
+ (GetPixelIntensity(edge_image,q) == 0))
+ {
+ *q=QuantumRange;
+ status=SyncCacheViewAuthenticPixels(trace_view,exception);
+ if (status == MagickFalse)
+ return(MagickFalse);
+ status=GetMatrixElement(pixel_cache,(ssize_t) number_edges,0,
+ &pixel);
+ if (status == MagickFalse)
+ return(MagickFalse);
+ pixel.x=x_offset+u;
+ pixel.y=y_offset+v;
+ status=SetMatrixElement(pixel_cache,(ssize_t) number_edges,0,
+ &pixel);
+ if (status == MagickFalse)
+ return(MagickFalse);
+ number_edges++;
+ }
+ }
+ }
+ } while (number_edges != 0);
+ return(MagickTrue);
+}
+
+
MagickExport Image *CannyEdgeImage(const Image *image,const double radius,
- const double sigma,const double low_factor,const double high_factor,
+ const double sigma,const double lower_precent,const double upper_percent,
ExceptionInfo *exception)
{
+ CacheView
+ *edge_view,
+ *trace_view;
+
+ CannyInfo
+ pixel;
+
+ char
+ geometry[MaxTextExtent];
+
+ double
+ lower_threshold,
+ upper_threshold;
+
Image
*edge_image;
+ KernelInfo
+ *kernel_info;
+
+ MagickBooleanType
+ status;
+
+ MatrixInfo
+ *pixel_cache;
+
+ register ssize_t
+ i;
+
+ size_t
+ *histogram,
+ number_pixels;
+
+ ssize_t
+ count,
+ y;
+
assert(image != (const Image *) NULL);
assert(image->signature == MagickSignature);
if (image->debug != MagickFalse)
(void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
assert(exception != (ExceptionInfo *) NULL);
assert(exception->signature == MagickSignature);
- edge_image=(Image *) NULL;
+ /*
+ Filter out noise.
+ */
+ (void) FormatLocaleString(geometry,MaxTextExtent,
+ "blur:%.20gx%.20g;blur:%.20gx%.20g+90",radius,sigma,radius,sigma);
+ kernel_info=AcquireKernelInfo(geometry);
+ if (kernel_info == (KernelInfo *) NULL)
+ ThrowImageException(ResourceLimitError,"MemoryAllocationFailed");
+ edge_image=MorphologyApply(image,ConvolveMorphology,1,kernel_info,
+ UndefinedCompositeOp,0.0,exception);
+ kernel_info=DestroyKernelInfo(kernel_info);
+ if (edge_image == (Image *) NULL)
+ return((Image *) NULL);
+ if (SetImageColorspace(edge_image,GRAYColorspace,exception) == MagickFalse)
+ {
+ edge_image=DestroyImage(edge_image);
+ return((Image *) NULL);
+ }
+ /*
+ Find the intensity gradient of the image.
+ */
+ pixel_cache=AcquireMatrixInfo(edge_image->columns,edge_image->rows,
+ sizeof(CannyInfo),exception);
+ if (pixel_cache == (MatrixInfo *) NULL)
+ {
+ edge_image=DestroyImage(edge_image);
+ return((Image *) NULL);
+ }
+ status=MagickTrue;
+ edge_view=AcquireVirtualCacheView(edge_image,exception);
+#if defined(MAGICKCORE_OPENMP_SUPPORT)
+ #pragma omp parallel for schedule(static,4) shared(status) \
+ magick_threads(edge_image,edge_image,edge_image->rows,1)
+#endif
+ for (y=0; y < (ssize_t) edge_image->rows; y++)
+ {
+ register const Quantum
+ *restrict p;
+
+ register ssize_t
+ x;
+
+ if (status == MagickFalse)
+ continue;
+ p=GetCacheViewVirtualPixels(edge_view,0,y,edge_image->columns+1,2,
+ exception);
+ if (p == (const Quantum *) NULL)
+ {
+ status=MagickFalse;
+ continue;
+ }
+ for (x=0; x < (ssize_t) edge_image->columns; x++)
+ {
+ CannyInfo
+ pixel;
+
+ double
+ dx,
+ dy;
+
+ int
+ orientation;
+
+ register const Quantum
+ *restrict kernel_pixels;
+
+ ssize_t
+ v;
+
+ static double
+ Gx[2][2] =
+ {
+ { -1.0, +1.0 },
+ { -1.0, +1.0 }
+ },
+ Gy[2][2] =
+ {
+ { +1.0, +1.0 },
+ { -1.0, -1.0 }
+ };
+
+ (void) ResetMagickMemory(&pixel,0,sizeof(pixel));
+ dx=0.0;
+ dy=0.0;
+ kernel_pixels=p;
+ for (v=0; v < 2; v++)
+ {
+ ssize_t
+ u;
+
+ for (u=0; u < 2; u++)
+ {
+ double
+ intensity;
+
+ intensity=GetPixelIntensity(edge_image,kernel_pixels+u);
+ dx+=0.5*Gx[v][u]*intensity;
+ dy+=0.5*Gy[v][u]*intensity;
+ }
+ kernel_pixels+=edge_image->columns+1;
+ }
+ pixel.magnitude=sqrt(dx*dx+dy*dy);
+ pixel.orientation=2;
+ if (fabs(dx) > MagickEpsilon)
+ {
+ double
+ theta;
+
+ theta=dy/dx;
+ if (theta < 0.0)
+ {
+ if (theta < -2.41421356237)
+ pixel.orientation=2;
+ else
+ if (theta < -0.414213562373)
+ pixel.orientation=1;
+ else
+ pixel.orientation=0;
+ }
+ else
+ {
+ if (theta > 2.41421356237)
+ pixel.orientation=2;
+ else
+ if (theta > 0.414213562373)
+ pixel.orientation=3;
+ else
+ pixel.orientation=0;
+ }
+ }
+ if (SetMatrixElement(pixel_cache,x,y,&pixel) == MagickFalse)
+ continue;
+ p+=GetPixelChannels(edge_image);
+ }
+ }
+ edge_view=DestroyCacheView(edge_view);
+ /*
+ Non-maxima suppression, remove pixels that are not considered to be part
+ of an edge.
+ */
+ (void) GetMatrixElement(pixel_cache,0,0,&pixel);
+ max=pixel.intensity;
+ min=pixel.intensity;
+ edge_view=AcquireAuthenticCacheView(edge_image,exception);
+#if defined(MAGICKCORE_OPENMP_SUPPORT)
+ #pragma omp parallel for schedule(static,4) shared(status) \
+ magick_threads(edge_image,edge_image,edge_image->rows,1)
+#endif
+ for (y=0; y < (ssize_t) edge_image->rows; y++)
+ {
+ register Quantum
+ *restrict q;
+
+ register ssize_t
+ x;
+
+ if (status == MagickFalse)
+ continue;
+ q=GetCacheViewAuthenticPixels(edge_view,0,y,edge_image->columns,1,
+ exception);
+ if (q == (Quantum *) NULL)
+ {
+ status=MagickFalse;
+ continue;
+ }
+ for (x=0; x < (ssize_t) edge_image->columns; x++)
+ {
+ CannyInfo
+ alpha_pixel,
+ beta_pixel,
+ pixel;
+
+ (void) GetMatrixElement(pixel_cache,x,y,&pixel);
+ switch (pixel.orientation)
+ {
+ case 0:
+ {
+ /*
+ 0 degrees.
+ */
+ (void) GetMatrixElement(pixel_cache,x-1,y,&alpha_pixel);
+ (void) GetMatrixElement(pixel_cache,x+1,y,&beta_pixel);
+ break;
+ }
+ case 1:
+ {
+ /*
+ 45 degrees.
+ */
+ (void) GetMatrixElement(pixel_cache,x-1,y-1,&alpha_pixel);
+ (void) GetMatrixElement(pixel_cache,x+1,y+1,&beta_pixel);
+ break;
+ }
+ case 2:
+ {
+ /*
+ 90 degrees.
+ */
+ (void) GetMatrixElement(pixel_cache,x,y-1,&alpha_pixel);
+ (void) GetMatrixElement(pixel_cache,x,y+1,&beta_pixel);
+ break;
+ }
+ case 3:
+ {
+ /*
+ 135 degrees.
+ */
+ (void) GetMatrixElement(pixel_cache,x-1,y+1,&alpha_pixel);
+ (void) GetMatrixElement(pixel_cache,x+1,y-1,&beta_pixel);
+ break;
+ }
+ }
+ pixel.intensity=pixel.magnitude;
+ if ((pixel.magnitude < alpha_pixel.magnitude) ||
+ (pixel.magnitude < beta_pixel.magnitude))
+ pixel.intensity=0;
+ (void) SetMatrixElement(pixel_cache,x,y,&pixel);
+#if defined(MAGICKCORE_OPENMP_SUPPORT)
+ #pragma omp critical (MagickCore_CannyEdgeImage)
+#endif
+ {
+ if (pixel.intensity < min)
+ min=pixel.intensity;
+ if (pixel.intensity > max)
+ max=pixel.intensity;
+ }
+ *q=0;
+ q+=GetPixelChannels(edge_image);
+ }
+ if (SyncCacheViewAuthenticPixels(edge_view,exception) == MagickFalse)
+ status=MagickFalse;
+ }
+ edge_view=DestroyCacheView(edge_view);
+ /*
+ Estimate hysteresis threshold.
+ */
+ lower_threshold=lower_percent*(max-min)+min;
+ upper_threshold=upper_percent*(max-min)+min;
+ /*
+ Hysteresis threshold.
+ */
+ edge_view=AcquireAuthenticCacheView(edge_image,exception);
+ trace_view=AcquireAuthenticCacheView(edge_image,exception);
+#if defined(MAGICKCORE_OPENMP_SUPPORT)
+ #pragma omp parallel for schedule(static,4) shared(status) \
+ magick_threads(edge_image,edge_image,edge_image->rows,1)
+#endif
+ for (y=0; y < (ssize_t) edge_image->rows; y++)
+ {
+ register ssize_t
+ x;
+
+ if (status == MagickFalse)
+ continue;
+ for (x=0; x < (ssize_t) edge_image->columns; x++)
+ {
+ CannyInfo
+ pixel;
+
+ register Quantum
+ *restrict q;
+
+ /*
+ Edge if pixel gradient higher than upper threshold.
+ */
+ status=GetMatrixElement(pixel_cache,x,y,&pixel);
+ if (status == MagickFalse)
+ break;
+ q=GetCacheViewAuthenticPixels(edge_view,x,y,1,1,exception);
+ if (q == (PixelPacket *) NULL)
+ {
+ status=MagickFalse;
+ continue;
+ }
+ if ((pixel.intensity >= upper_threshold) &&
+ (GetPixelIntensity(edge_image,q) == 0))
+ {
+ *q=QuantumRange;
+ status=SyncCacheViewAuthenticPixels(edge_view,exception);
+ if (status == MagickFalse)
+ continue;
+ status=TraceEdges(edge_image,trace_view,pixel_cache,x,y,
+ lower_threshold,exception);
+ if (status == MagickFalse)
+ continue;
+ }
+ }
+ if (SyncCacheViewAuthenticPixels(edge_view,exception) == MagickFalse)
+ status=MagickFalse;
+ }
+ trace_view=DestroyCacheView(trace_view);
+ edge_view=DestroyCacheView(edge_view);
+ pixel_cache=DestroyMatrixInfo(pixel_cache);
return(edge_image);
}
\f
projects / imagemagick / blobdiff