/* [<][>][^][v][top][bottom][index][help] */
DEFINITIONS
This source file includes following definitions.
- PlasmaPixel
- PlasmaImage
/*
% Copyright (C) 2003 - 2009 GraphicsMagick Group
% Copyright (C) 2003 ImageMagick Studio
% Copyright 1991-1999 E. I. du Pont de Nemours and Company
%
% This program is covered by multiple licenses, which are described in
% Copyright.txt. You should have received a copy of Copyright.txt with this
% package; otherwise see http://www.graphicsmagick.org/www/Copyright.html.
%
% GraphicsMagick Plasma Methods.
%
*/
/*
Include declarations.
*/
#include "magick/studio.h"
#include "magick/monitor.h"
#include "magick/pixel_cache.h"
#include "magick/random.h"
#include "magick/utility.h"
/*
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %
% %
% P l a s m a I m a g e %
% %
% %
% %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% PlasmaImage() initializes an image with plasma fractal values. The image
% must be initialized with a base color and the random number generator
% seeded before this method is called.
%
% The format of the PlasmaImage method is:
%
% MagickPassFail PlasmaImage(Image *image,const SegmentInfo *segment,
% unsigned long attenuate,unsigned long depth)
%
% A description of each parameter follows:
%
% o image: The image.
%
% o segment: Define the region to apply plasma fractals values.
%
% o attenuate: Define the plasma attenuation factor.
%
% o depth: Limit the plasma recursion depth.
%
%
*/
static Quantum PlasmaPixel(const double pixel,const double noise)
{
double
value;
value=pixel+noise*MagickRandomReal()-noise/2;
if (value <= 0.0)
return(0);
if (value >= MaxRGB)
return(MaxRGB);
return((Quantum) (value+0.5));
}
MagickExport MagickPassFail PlasmaImage(Image *image,const SegmentInfo *segment,
unsigned long attenuate,unsigned long depth)
{
double
plasma;
long
x,
x_mid,
y,
y_mid;
PixelPacket
u,
v;
register PixelPacket
*q;
assert(image != (Image *) NULL);
assert(image->signature == MagickSignature);
assert(segment != (SegmentInfo *) NULL);
if (depth != 0)
{
SegmentInfo
local_info;
/*
Divide the area into quadrants and recurse.
*/
depth--;
attenuate++;
x_mid=(long) (segment->x1+segment->x2+0.5)/2;
y_mid=(long) (segment->y1+segment->y2+0.5)/2;
local_info=(*segment);
local_info.x2=x_mid;
local_info.y2=y_mid;
(void) PlasmaImage(image,&local_info,attenuate,depth);
local_info=(*segment);
local_info.y1=y_mid;
local_info.x2=x_mid;
(void) PlasmaImage(image,&local_info,attenuate,depth);
local_info=(*segment);
local_info.x1=x_mid;
local_info.y2=y_mid;
(void) PlasmaImage(image,&local_info,attenuate,depth);
local_info=(*segment);
local_info.x1=x_mid;
local_info.y1=y_mid;
return(PlasmaImage(image,&local_info,attenuate,depth));
}
if ((image->colorspace != RGBColorspace) ||
(image->storage_class != DirectClass))
{
(void) SetImageType(image,TrueColorType);
}
x_mid=(long) (segment->x1+segment->x2+0.5)/2;
y_mid=(long) (segment->y1+segment->y2+0.5)/2;
if ((segment->x1 == x_mid) && (segment->x2 == x_mid) &&
(segment->y1 == y_mid) && (segment->y2 == y_mid))
return(False);
/*
Average pixels and apply plasma.
*/
plasma=MaxRGB/(2.0*attenuate);
if ((segment->x1 != x_mid) || (segment->x2 != x_mid))
{
/*
Left pixel.
*/
x=(long) (segment->x1+0.5);
(void) AcquireOnePixelByReference(image,&u,x,(long) (segment->y1+0.5),&image->exception);
(void) AcquireOnePixelByReference(image,&v,x,(long) (segment->y2+0.5),&image->exception);
q=SetImagePixels(image,x,y_mid,1,1);
if (q == (PixelPacket *) NULL)
return(True);
q->red=PlasmaPixel(((double) u.red+v.red)/2,plasma);
q->green=PlasmaPixel(((double) u.green+v.green)/2,plasma);
q->blue=PlasmaPixel(((double) u.blue+v.blue)/2,plasma);
(void) SyncImagePixels(image);
if (segment->x1 != segment->x2)
{
/*
Right pixel.
*/
x=(long) (segment->x2+0.5);
(void) AcquireOnePixelByReference(image,&u,x,(long) (segment->y1+0.5),&image->exception);
(void) AcquireOnePixelByReference(image,&v,x,(long) (segment->y2+0.5),&image->exception);
q=SetImagePixels(image,x,y_mid,1,1);
if (q == (PixelPacket *) NULL)
return(True);
q->red=PlasmaPixel(((double) u.red+v.red)/2,plasma);
q->green=PlasmaPixel(((double) u.green+v.green)/2,plasma);
q->blue=PlasmaPixel(((double) u.blue+v.blue)/2,plasma);
(void) SyncImagePixels(image);
}
}
if ((segment->y1 != y_mid) || (segment->y2 != y_mid))
{
if ((segment->x1 != x_mid) || (segment->y2 != y_mid))
{
/*
Bottom pixel.
*/
y=(long) (segment->y2+0.5);
(void) AcquireOnePixelByReference(image,&u,(long) (segment->x1+0.5),y,&image->exception);
(void) AcquireOnePixelByReference(image,&v,(long) (segment->x2+0.5),y,&image->exception);
q=SetImagePixels(image,x_mid,y,1,1);
if (q == (PixelPacket *) NULL)
return(True);
q->red=PlasmaPixel(((double) u.red+v.red)/2,plasma);
q->green=PlasmaPixel(((double) u.green+v.green)/2,plasma);
q->blue=PlasmaPixel(((double) u.blue+v.blue)/2,plasma);
(void) SyncImagePixels(image);
}
if (segment->y1 != segment->y2)
{
/*
Top pixel.
*/
y=(long) (segment->y1+0.5);
(void) AcquireOnePixelByReference(image,&u,(long) (segment->x1+0.5),y,&image->exception);
(void) AcquireOnePixelByReference(image,&v,(long) (segment->x2+0.5),y,&image->exception);
q=SetImagePixels(image,x_mid,y,1,1);
if (q == (PixelPacket *) NULL)
return(True);
q->red=PlasmaPixel(((double) u.red+v.red)/2,plasma);
q->green=PlasmaPixel(((double) u.green+v.green)/2,plasma);
q->blue=PlasmaPixel(((double) u.blue+v.blue)/2,plasma);
(void) SyncImagePixels(image);
}
}
if ((segment->x1 != segment->x2) || (segment->y1 != segment->y2))
{
/*
Middle pixel.
*/
x=(long) (segment->x1+0.5);
y=(long) (segment->y1+0.5);
(void) AcquireOnePixelByReference(image,&u,x,y,&image->exception);
x=(long) (segment->x2+0.5);
y=(long) (segment->y2+0.5);
(void) AcquireOnePixelByReference(image,&v,x,y,&image->exception);
q=SetImagePixels(image,x_mid,y_mid,1,1);
if (q == (PixelPacket *) NULL)
return(True);
q->red=PlasmaPixel(((double) u.red+v.red)/2,plasma);
q->green=PlasmaPixel(((double) u.green+v.green)/2,plasma);
q->blue=PlasmaPixel(((double) u.blue+v.blue)/2,plasma);
(void) SyncImagePixels(image);
}
if (((segment->x2-segment->x1) < 3.0) && ((segment->y2-segment->y1) < 3.0))
return(True);
return(False);
}