//Saving the actual image
G_MODULE_EXPORT void cb_save_clicked(GtkButton *button)
{
//Check if there is an image
if(camera_params.image_number==0)
{
gtk_widget_show( camera_params.objects->no_image_dialog );
return;
}
camera_params.wand_data.saving_wand=NULL;
if(strcmp("processed_save_button",gtk_buildable_get_name(GTK_BUILDABLE(button)))==0)
{
pthread_mutex_lock(&camera_params.wand_data.processed_img_mutex);
camera_params.wand_data.saving_wand=CloneMagickWand(camera_params.wand_data.processed_magick_wand);
pthread_mutex_unlock(&camera_params.wand_data.processed_img_mutex);
}
else if(strcmp("raw_save_button",gtk_buildable_get_name(GTK_BUILDABLE(button)))==0)
{
pthread_mutex_lock(&camera_params.wand_data.raw_img_mutex);
camera_params.wand_data.saving_wand=CloneMagickWand(camera_params.wand_data.raw_magick_wand);
pthread_mutex_unlock(&camera_params.wand_data.raw_img_mutex);
//we scale the image on 16 bits before saving
imagemagick_levelimage(camera_params.wand_data.saving_wand, 0, 1<<((int)camera_params.sensorbits));
}
if(camera_params.wand_data.saving_wand)
{
gtk_widget_show( camera_params.objects->imagesavedialog );
g_print("save button %s clicked\n",gtk_buildable_get_name(GTK_BUILDABLE(button)));
}
}
static PyObject *
im_CloneImage(PyObject *self, PyObject *args) {
void *magick_pointer;
MagickWand *magick_wand;
MagickWand *wand_clone;
if (!PyArg_ParseTuple(args, "O", &magick_pointer)){
Py_INCREF(Py_None);
return Py_None;
}
magick_wand = (MagickWand *) PyCObject_AsVoidPtr(magick_pointer);
wand_clone = CloneMagickWand(magick_wand);
return Py_BuildValue("O", PyCObject_FromVoidPtr((void *)wand_clone,
(void *)DestroyMagickWand));
}
void test_wand(void)
{
MagickWand *magick_wand = NULL;
MagickWand *c_wand = NULL;
DrawingWand *d_wand = NULL;
PixelWand *p_wand = NULL;
// Used for text effect #3
double dargs[1] = {120.};
// Used for text effect #5
double d_args[8] = {
-0.02,0.0,
0.0,1.02,
0.0,0.0,
-0.5,1.9
};
MagickWandGenesis();
// Text effect 1 - shadow effect using MagickShadowImage
// This is derived from Anthony's Soft Shadow effect
// convert -size 300x100 xc:none -font Candice -pointsize 72 \
// -fill white -stroke black -annotate +25+65 'Anthony' \
// \( +clone -background navy -shadow 70x4+5+5 \) +swap \
// -background lightblue -flatten -trim +repage font_shadow_soft.jpg
//NOTE - if an image has a transparent background, adding a border of any colour other
// than "none" will remove all the transparency and replace it with the border's colour
magick_wand = NewMagickWand();
d_wand = NewDrawingWand();
p_wand = NewPixelWand();
PixelSetColor(p_wand,"none");
// Create a new transparent image
MagickNewImage(magick_wand,350,100,p_wand);
// Set up a 72 point white font
PixelSetColor(p_wand,"white");
DrawSetFillColor(d_wand,p_wand);
DrawSetFont (d_wand, "Verdana-Bold-Italic" ) ;
DrawSetFontSize(d_wand,72);
// Add a black outline to the text
PixelSetColor(p_wand,"black");
DrawSetStrokeColor(d_wand,p_wand);
// Turn antialias on - not sure this makes a difference
DrawSetTextAntialias(d_wand,MagickTrue);
// Now draw the text
DrawAnnotation(d_wand,25,65,(const unsigned char *)"Magick");
// Draw the image on to the magick_wand
MagickDrawImage(magick_wand,d_wand);
// Trim the image down to include only the text
MagickTrimImage(magick_wand,0);
// equivalent to the command line +repage
MagickResetImagePage(magick_wand,"");
// Make a copy of the text image
c_wand = CloneMagickWand(magick_wand);
// Set the background colour to blue for the shadow
PixelSetColor(p_wand,"blue");
MagickSetImageBackgroundColor(magick_wand,p_wand);
// Opacity is a real number indicating (apparently) percentage
MagickShadowImage(magick_wand,70,4,5,5);
// Composite the text on top of the shadow
MagickCompositeImage(magick_wand,c_wand,OverCompositeOp,5,5);
if(c_wand)c_wand = DestroyMagickWand(c_wand);
c_wand = NewMagickWand();
// Create a new image the same size as the text image and put a solid colour
// as its background
PixelSetColor(p_wand,"rgb(125,215,255)");
MagickNewImage(c_wand,MagickGetImageWidth(magick_wand),MagickGetImageHeight(magick_wand),p_wand);
// Now composite the shadowed text over the plain background
MagickCompositeImage(c_wand,magick_wand,OverCompositeOp,0,0);
// and write the result
MagickWriteImage(c_wand,"text_shadow.png");
/* Clean up */
if(magick_wand)magick_wand = DestroyMagickWand(magick_wand);
if(c_wand)c_wand = DestroyMagickWand(c_wand);
if(d_wand)d_wand = DestroyDrawingWand(d_wand);
if(p_wand)p_wand = DestroyPixelWand(p_wand);
// Text effect 2 - tiled text using the builtin checkerboard pattern
// Anthony's Tiled Font effect
// convert -size 320x100 xc:lightblue -font Candice -pointsize 72 \
// -tile pattern:checkerboard -annotate +28+68 'Anthony' \
// font_tile.jpg
magick_wand = NewMagickWand();
d_wand = NewDrawingWand();
p_wand = NewPixelWand();
set_tile_pattern(d_wand,"#check","pattern:checkerboard");
PixelSetColor(p_wand,"lightblue");
// Create a new transparent image
MagickNewImage(magick_wand,320,100,p_wand);
// Set up a 72 point font
DrawSetFont (d_wand, "Verdana-Bold-Italic" ) ;
DrawSetFontSize(d_wand,72);
// Now draw the text
DrawAnnotation(d_wand,28,68,(const unsigned char *)"Magick");
// Draw the image on to the magick_wand
MagickDrawImage(magick_wand,d_wand);
// Trim the image
MagickTrimImage(magick_wand,0);
// Add a transparent border
PixelSetColor(p_wand,"lightblue");
MagickBorderImage(magick_wand,p_wand,5,5);
// and write it
MagickWriteImage(magick_wand,"text_pattern.png");
/* Clean up */
if(magick_wand)magick_wand = DestroyMagickWand(magick_wand);
if(d_wand)d_wand = DestroyDrawingWand(d_wand);
if(p_wand)p_wand = DestroyPixelWand(p_wand);
// Text effect 3 - arc font (similar to http://www.imagemagick.org/Usage/fonts/#arc)
//convert -size 320x100 xc:lightblue -font Candice -pointsize 72 \
// -annotate +25+65 'Anthony' -distort Arc 120 \
// -trim +repage -bordercolor lightblue -border 10 font_arc.jpg
magick_wand = NewMagickWand();
d_wand = NewDrawingWand();
p_wand = NewPixelWand();
// Create a 320x100 lightblue canvas
PixelSetColor(p_wand,"lightblue");
MagickNewImage(magick_wand,320,100,p_wand);
// Set up a 72 point font
DrawSetFont (d_wand, "Verdana-Bold-Italic" ) ;
DrawSetFontSize(d_wand,72);
// Now draw the text
DrawAnnotation(d_wand,25,65,(const unsigned char *)"Magick");
// Draw the image on to the magick_wand
MagickDrawImage(magick_wand,d_wand);
MagickDistortImage(magick_wand,ArcDistortion,1,dargs,MagickFalse);
// Trim the image
MagickTrimImage(magick_wand,0);
// Add the border
PixelSetColor(p_wand,"lightblue");
MagickBorderImage(magick_wand,p_wand,10,10);
// and write it
MagickWriteImage(magick_wand,"text_arc.png");
/* Clean up */
if(magick_wand)magick_wand = DestroyMagickWand(magick_wand);
if(d_wand)d_wand = DestroyDrawingWand(d_wand);
if(p_wand)p_wand = DestroyPixelWand(p_wand);
// Text effect 4 - bevelled font http://www.imagemagick.org/Usage/fonts/#bevel
// convert -size 320x100 xc:black -font Candice -pointsize 72 \
// -fill white -annotate +25+65 'Anthony' \
// -shade 140x60 font_beveled.jpg
magick_wand = NewMagickWand();
d_wand = NewDrawingWand();
p_wand = NewPixelWand();
// Create a 320x100 canvas
PixelSetColor(p_wand,"gray");
MagickNewImage(magick_wand,320,100,p_wand);
// Set up a 72 point font
DrawSetFont (d_wand, "Verdana-Bold-Italic" ) ;
DrawSetFontSize(d_wand,72);
// Set up a 72 point white font
PixelSetColor(p_wand,"white");
DrawSetFillColor(d_wand,p_wand);
// Now draw the text
DrawAnnotation(d_wand,25,65,(const unsigned char *)"Magick");
// Draw the image on to the magick_wand
MagickDrawImage(magick_wand,d_wand);
// the "gray" parameter must be true to get the effect shown on Anthony's page
MagickShadeImage(magick_wand,MagickTrue,140,60);
#ifdef COLORIZE
PixelSetColor(p_wand,"yellow");
DrawSetFillColor(d_wand,p_wand);
cp_wand = NewPixelWand();
PixelSetColor(cp_wand,"gold");
MagickColorizeImage(magick_wand,p_wand,cp_wand);
#endif
// and write it
MagickWriteImage(magick_wand,"text_bevel.png");
/* Clean up */
if(magick_wand)magick_wand = DestroyMagickWand(magick_wand);
if(d_wand)d_wand = DestroyDrawingWand(d_wand);
if(p_wand)p_wand = DestroyPixelWand(p_wand);
#ifdef COLORIZE
if(cp_wand)cp_wand = DestroyPixelWand(cp_wand);
#endif
// Text effect 5 and 6 - Plain text and then Barrel distortion
// This one uses d_args
magick_wand = NewMagickWand();
d_wand = NewDrawingWand();
p_wand = NewPixelWand();
// Create a 320x100 transparent canvas
PixelSetColor(p_wand,"none");
MagickNewImage(magick_wand,320,100,p_wand);
// Set up a 72 point font
DrawSetFont (d_wand, "Verdana-Bold-Italic" ) ;
DrawSetFontSize(d_wand,72);
// Now draw the text
DrawAnnotation(d_wand,25,65,(const unsigned char *)"Magick");
// Draw the image on to the magick_wand
MagickDrawImage(magick_wand,d_wand);
MagickWriteImage(magick_wand,"text_plain.png");
// Trim the image
MagickTrimImage(magick_wand,0);
// Add the border
PixelSetColor(p_wand,"none");
MagickBorderImage(magick_wand,p_wand,10,10);
// MagickSetImageMatte(magick_wand,MagickTrue);
// MagickSetImageVirtualPixelMethod(magick_wand,TransparentVirtualPixelMethod);
// d_args[0] = 0.1;d_args[1] = -0.25;d_args[2] = -0.25; [3] += .1
// The first value should be positive. If it is negative the image is *really* distorted
d_args[0] = 0.0;
d_args[1] = 0.0;
d_args[2] = 0.5;
// d_args[3] should normally be chosen such the sum of all 4 values is 1
// so that the result is the same size as the original
// You can override the sum with a different value
// If the sum is greater than 1 the resulting image will be smaller than the original
d_args[3] = 1 - (d_args[0] + d_args[1] + d_args[2]);
// Make the result image smaller so that it isn't as likely
// to overflow the edges
// d_args[3] += 0.1;
// 0.0,0.0,0.5,0.5,0.0,0.0,-0.5,1.9
d_args[3] = 0.5;
d_args[4] = 0.0;
d_args[5] = 0.0;
d_args[6] = -0.5;
d_args[7] = 1.9;
// DON'T FORGET to set the correct number of arguments here
MagickDistortImage(magick_wand,BarrelDistortion,8,d_args,MagickTrue);
// MagickResetImagePage(magick_wand,"");
// Trim the image again
MagickTrimImage(magick_wand,0);
// Add the border
PixelSetColor(p_wand,"none");
MagickBorderImage(magick_wand,p_wand,10,10);
// and write it
MagickWriteImage(magick_wand,"text_barrel.png");
/* Clean up */
if(magick_wand)magick_wand = DestroyMagickWand(magick_wand);
if(d_wand)d_wand = DestroyDrawingWand(d_wand);
if(p_wand)p_wand = DestroyPixelWand(p_wand);
// Text effect 7 - Polar distortion
// This one uses d_args[0]
magick_wand = NewMagickWand();
d_wand = NewDrawingWand();
p_wand = NewPixelWand();
// Create a 320x200 transparent canvas
PixelSetColor(p_wand,"none");
MagickNewImage(magick_wand,320,200,p_wand);
// Set up a 72 point font
DrawSetFont (d_wand, "Verdana-Bold-Italic" ) ;
DrawSetFontSize(d_wand,72);
// Now draw the text
DrawAnnotation(d_wand,25,65,(const unsigned char *)"Magick");
// Draw the image on to the magick_wand
MagickDrawImage(magick_wand,d_wand);
d_args[0] = 0.0;
// DON'T FORGET to set the correct number of arguments here
MagickDistortImage(magick_wand,PolarDistortion,1,d_args,MagickTrue);
// MagickResetImagePage(magick_wand,"");
// Trim the image again
MagickTrimImage(magick_wand,0);
// Add the border
PixelSetColor(p_wand,"none");
MagickBorderImage(magick_wand,p_wand,10,10);
// and write it
MagickWriteImage(magick_wand,"text_polar.png");
/* Clean up */
if(magick_wand)magick_wand = DestroyMagickWand(magick_wand);
if(d_wand)d_wand = DestroyDrawingWand(d_wand);
if(p_wand)p_wand = DestroyPixelWand(p_wand);
// Text effect 8 - Shepard's distortion
// This one uses d_args[0]
magick_wand = NewMagickWand();
d_wand = NewDrawingWand();
p_wand = NewPixelWand();
// Create a 320x200 transparent canvas
PixelSetColor(p_wand,"none");
MagickNewImage(magick_wand,640,480,p_wand);
// Set up a 72 point font
DrawSetFont (d_wand, "Verdana-Bold-Italic" ) ;
DrawSetFontSize(d_wand,72);
// Now draw the text
DrawAnnotation(d_wand,50,240,(const unsigned char *)"Magick Rocks");
// Draw the image on to the magick_wand
MagickDrawImage(magick_wand,d_wand);
d_args[0] = 150.0;
d_args[1] = 190.0;
d_args[2] = 100.0;
d_args[3] = 290.0;
d_args[4] = 500.0;
d_args[5] = 200.0;
d_args[6] = 430.0;
d_args[7] = 130.0;
// DON'T FORGET to set the correct number of arguments here
MagickDistortImage(magick_wand,ShepardsDistortion,8,d_args,MagickTrue);
// Trim the image
MagickTrimImage(magick_wand,0);
// Add the border
PixelSetColor(p_wand,"none");
MagickBorderImage(magick_wand,p_wand,10,10);
// and write it
MagickWriteImage(magick_wand,"text_shepards.png");
/* Clean up */
if(magick_wand)magick_wand = DestroyMagickWand(magick_wand);
if(d_wand)d_wand = DestroyDrawingWand(d_wand);
if(p_wand)p_wand = DestroyPixelWand(p_wand);
MagickWandTerminus();
}
/*
* given a source image, a destination filepath and a set of image metadata thresholds,
* search for the lowest-quality version of the source image whose properties fall within our
* thresholds.
* this will produce an image file that looks the same to the casual observer, but which
* contains much less information and results in a smaller file.
* typical savings on unoptimized images vary widely from 10-80%, with 25-50% being most common.
*/
MagickWand * search_quality(MagickWand *mw, const char *dst,
const struct imgmin_options *opt)
{
MagickWand *tmp = NULL;
char tmpfile[MAX_PATH] = "/tmp/imgminXXXXXX";
if (0 == strcmp("-", dst))
{
if (-1 == mkstemp(tmpfile))
{
perror("mkstemp");
return CloneMagickWand(mw);
}
} else {
strcpy(tmpfile, dst);
}
/*
* The overwhelming majority of JPEGs are TrueColorType; it is those types, with a low
* unique color count, that we must avoid.
*/
if (!enough_colors(mw, opt))
{
fprintf(stdout, " Color count is too low, skipping...\n");
return CloneMagickWand(mw);
}
if (quality(mw) < opt->quality_in_min)
{
fprintf(stdout, " Quality < %u, won't second-guess...\n", opt->quality_in_min);
return CloneMagickWand(mw);
}
size_t width = MagickGetImageWidth(mw);
size_t height = MagickGetImageHeight(mw);
dssim_info *dssim = dssim_init(1);
void *convert_data = convert_row_start(mw);
dssim_set_original_float_callback(dssim, width, height, convert_row_callback, convert_data);
convert_row_finish(convert_data);
{
ExceptionInfo *exception = AcquireExceptionInfo();
const double original_density = color_density(mw);
unsigned qmax = min(quality(mw), opt->quality_out_max);
unsigned qmin = opt->quality_out_min;
unsigned steps = 0;
/*
* binary search of quality space for optimally lowest quality that
* produces an acceptable level of distortion
*/
while (qmax > qmin + 1 && steps < opt->max_steps)
{
double density_ratio;
unsigned q;
steps++;
q = (qmax + qmin) / 2;
/* change quality */
tmp = CloneMagickWand(mw);
MagickSetImageCompressionQuality(tmp, q);
/* apply quality change */
MagickWriteImages(tmp, tmpfile, MagickTrue);
DestroyMagickWand(tmp);
tmp = NewMagickWand();
MagickReadImage(tmp, tmpfile);
void *convert_data = convert_row_start(tmp);
dssim_set_modified_float_callback(dssim, width, height, convert_row_callback, convert_data);
convert_row_finish(convert_data);
double error = 20.0 * dssim_compare(dssim, NULL); // scaled to threshold of previous implementation
density_ratio = fabs(color_density(tmp) - original_density) / original_density;
/* color density ratio threshold is an alternative quality measure.
If it's exceeded, pretend MSE was higher to increase quality */
if (density_ratio > opt->color_density_ratio) {
error *= 1.25 + density_ratio; // fudge factor
}
/* eliminate half search space based on whether distortion within thresholds */
if (error > opt->error_threshold)
{
qmin = q;
} else {
qmax = q;
}
if (opt->show_progress)
{
fprintf(stdout, "%.2f/%.2f@%u ", error, density_ratio, q);
}
/* Stop searching if close enough to the target */
if (fabs(error - opt->error_threshold) < opt->error_threshold * ERROR_THRESHOLD_INACCURACY) {
qmax = q;
break;
}
}
if (opt->show_progress)
{
putc('\n', stdout);
}
MagickSetImageCompressionQuality(mw, qmax);
/* "Chroma sub-sampling works because human vision is relatively insensitive to
* small areas of colour. It gives a significant reduction in file sizes, with
* little loss of perceived quality." [3]
*/
#if MagickLibVersion >= 0x630 /* FIXME: available in 0x660, not available in 0x628, not sure which version it was introduced in */
(void) MagickSetImageProperty(mw, "jpeg:sampling-factor", "2x2");
#endif
/* strip an image of all profiles and comments */
(void) MagickStripImage(mw);
MagickWriteImages(mw, tmpfile, MagickTrue);
(void) DestroyMagickWand(tmp);
tmp = NewMagickWand();
MagickReadImage(tmp, tmpfile);
exception = DestroyExceptionInfo(exception);
}
return tmp;
}
/*
* given a source image, a destination filepath and a set of image metadata thresholds,
* search for the lowest-quality version of the source image whose properties fall within our
* thresholds.
* this will produce an image file that looks the same to the casual observer, but which
* contains much less information and results in a smaller file.
* typical savings on unoptimized images vary widely from 10-80%, with 25-50% being most common.
*/
MagickWand * search_quality(MagickWand *mw, const char *dst,
const struct imgmin_options *opt)
{
MagickWand *tmp = NULL;
char tmpfile[MAX_PATH] = "/tmp/imgminXXXXXX";
if (0 == strcmp("-", dst))
{
if (-1 == mkstemp(tmpfile))
{
perror("mkstemp");
return CloneMagickWand(mw);
}
} else {
strcpy(tmpfile, dst);
}
/*
* The overwhelming majority of JPEGs are TrueColorType; it is those types, with a low
* unique color count, that we must avoid.
*/
if (unique_colors(mw) < opt->min_unique_colors && MagickGetType(mw) != GrayscaleType)
{
fprintf(stderr, " Color count is too low, skipping...\n");
return CloneMagickWand(mw);
}
if (quality(mw) < opt->quality_in_min)
{
fprintf(stderr, " Quality < %u, won't second-guess...\n", opt->quality_in_min);
return CloneMagickWand(mw);
}
{
ExceptionInfo *exception = AcquireExceptionInfo();
const double original_density = color_density(mw);
unsigned qmax = min(quality(mw), opt->quality_out_max);
unsigned qmin = opt->quality_out_min;
unsigned steps = 0;
/*
* binary search of quality space for optimally lowest quality that
* produces an acceptable level of distortion
*/
while (qmax > qmin + 1 && steps < opt->max_steps)
{
double distortion[CompositeChannels+1];
double error;
double density_ratio;
unsigned q;
steps++;
q = (qmax + qmin) / 2;
/* change quality */
tmp = CloneMagickWand(mw);
MagickSetImageCompressionQuality(tmp, q);
/* apply quality change */
MagickWriteImages(tmp, tmpfile, MagickTrue);
DestroyMagickWand(tmp);
tmp = NewMagickWand();
MagickReadImage(tmp, tmpfile);
/* quantify distortion produced by quality change
* NOTE: writes to distortion[], which we don't care about
* and to image(tmp)->error, which we do care about
*/
(void) GetImageDistortion(GetImageFromMagickWand(tmp),
GetImageFromMagickWand(mw),
MeanErrorPerPixelMetric,
distortion,
exception);
/* FIXME: in perlmagick i was getting back a number [0,255.0],
* here something else is happening... the hardcoded divisor
* is an imperfect attempt to scale the number back to the
* perlmagick results. I really need to look into this.
*/
error = GetImageFromMagickWand(tmp)->error.mean_error_per_pixel / 380.;
density_ratio = fabs(color_density(tmp) - original_density) / original_density;
/* eliminate half search space based on whether distortion within thresholds */
if (error > opt->error_threshold || density_ratio > opt->color_density_ratio)
{
qmin = q;
} else {
qmax = q;
}
if (opt->show_progress)
{
fprintf(stderr, "%.2f/%.2f@%u ", error, density_ratio, q);
}
}
if (opt->show_progress)
{
putc('\n', stderr);
}
MagickSetImageCompressionQuality(mw, qmax);
/* "Chroma sub-sampling works because human vision is relatively insensitive to
* small areas of colour. It gives a significant reduction in file sizes, with
* little loss of perceived quality." [3]
*/
(void) MagickSetImageProperty(mw, "jpeg:sampling-factor", "2x2");
/* strip an image of all profiles and comments */
(void) MagickStripImage(mw);
MagickWriteImages(mw, tmpfile, MagickTrue);
(void) DestroyMagickWand(tmp);
tmp = NewMagickWand();
MagickReadImage(tmp, tmpfile);
exception = DestroyExceptionInfo(exception);
}
return tmp;
}
int main(int argc, char **argv)
{
printf("Hello World\n");
#define QuantumScale ((MagickRealType) 1.0/(MagickRealType) QuantumRange)
#define SigmoidalContrast(x) (QuantumRange*(1.0/(1+exp(10.0*(0.5-QuantumScale*x)))-0.0066928509)*1.0092503)
#define ThrowWandException(wand)\
{\
char *description;\
ExceptionType\
severity;\
\
description=MagickGetException(wand,&severity);\
(void)fprintf(stderr,"%s %s %lu %s\n",GetMagickModule(),description);\
exit(-1);\
}
long y;
MagickBooleanType status;
MagickPixelPacket pixel;
MagickWand *contrast_wand, *image_wand;
PixelIterator *contrast_iterator,*iterator;
PixelWand **contrast_pixels, **pixels;
register long x;
unsigned long width;
if (argc != 3)
{
(void)fprintf(stdout,"Usage: %s image sigmoidal-image\n",argv[0]);
exit(0);
}
/*
Read an Image
*/
MagickWandGenesis();
image_wand = NewMagickWand();
status = MagickReadImage(image_wand,argv[1]);
if(status == MagickFalse)
ThrowWandException(image_wand);
contrast_wand = CloneMagickWand(image_wand);
/*
Sigmoidal non-linearity contrast control.
*/
iterator = NewPixelIterator(image_wand);
contrast_iterator=NewPixelIterator(contrast_wand);
if((iterator == (PixelIterator*)NULL) || (contrast_iterator == (PixelIterator*)NULL))
ThrowWandException(image_wand);
for(y=0; y<(long)MagickGetImageHeight(image_wand); y++)
{
pixels = PixelGetNextIteratorRow(iterator,&width);
contrast_pixels = PixelGetNextIteratorRow(contrast_iterator,&width);
if((pixels==(PixelWand**)NULL) || (contrast_pixels == (PixelWand**)NULL))
break;
for(x = 0; x<(long)width; x++)
{
PixelGetMagickColor(pixels[x],&pixel);
pixel.red = SigmoidalContrast(pixel.red);
pixel.green = SigmoidalContrast(pixel.green);
pixel.blue = SigmoidalContrast(pixel.blue);
pixel.index = SigmoidalContrast(pixel.index);
PixelSetMagickColor(contrast_pixels[x],&pixel);
}
(void)PixelSyncIterator(contrast_iterator);
}
if(y < (long)MagickGetImageHeight(image_wand))
ThrowWandException(image_wand);
contrast_iterator = DestroyPixelIterator(contrast_iterator);
iterator = DestroyPixelIterator(iterator);
image_wand=DestroyMagickWand(image_wand);
status = MagickWriteImages(contrast_wand,argv[2],MagickTrue);
if(status == MagickFalse)
ThrowWandException(image_wand);
contrast_wand = DestroyMagickWand(contrast_wand);
MagickWandTerminus();
return 0;
}
