private unsafe void ProcessImage(UnmanagedImage image, byte* mask, int maskLineSize)
{
// get image dimension
int width = image.Width;
int height = image.Height;
pixels = pixelsWithoutBlack = 0;
int[] s = new int[256];
int[] l = new int[256];
int[] swb = new int[256];
int[] lwb = new int[256];
RGB rgb = new RGB();
HSL hsl = new HSL();
int pixelSize = (image.PixelFormat == PixelFormat.Format24bppRgb) ? 3 : 4;
int offset = image.Stride - width * pixelSize;
int maskOffset = maskLineSize - width;
// do the job
byte* p = (byte*)image.ImageData.ToPointer();
if (mask == null)
{
// for each line
for (int y = 0; y < height; y++)
{
// for each pixel
for (int x = 0; x < width; x++, p += pixelSize)
{
rgb.Red = p[RGB.R];
rgb.Green = p[RGB.G];
rgb.Blue = p[RGB.B];
// convert to HSL color space
Accord.Imaging.HSL.FromRGB(rgb, hsl);
s[(int)(hsl.Saturation * 255)]++;
l[(int)(hsl.Luminance * 255)]++;
pixels++;
if (hsl.Luminance != 0.0)
{
swb[(int)(hsl.Saturation * 255)]++;
lwb[(int)(hsl.Luminance * 255)]++;
pixelsWithoutBlack++;
}
}
p += offset;
}
}
else
{
// for each line
for (int y = 0; y < height; y++)
{
// for each pixel
for (int x = 0; x < width; x++, p += pixelSize, mask++)
{
if (*mask == 0)
continue;
rgb.Red = p[RGB.R];
rgb.Green = p[RGB.G];
rgb.Blue = p[RGB.B];
// convert to HSL color space
Accord.Imaging.HSL.FromRGB(rgb, hsl);
s[(int)(hsl.Saturation * 255)]++;
l[(int)(hsl.Luminance * 255)]++;
pixels++;
if (hsl.Luminance != 0.0)
{
swb[(int)(hsl.Saturation * 255)]++;
lwb[(int)(hsl.Luminance * 255)]++;
pixelsWithoutBlack++;
}
}
p += offset;
mask += maskOffset;
}
}
// create histograms
saturation = new ContinuousHistogram(s, new Range(0, 1));
luminance = new ContinuousHistogram(l, new Range(0, 1));
saturationWithoutBlack = new ContinuousHistogram(swb, new Range(0, 1));
luminanceWithoutBlack = new ContinuousHistogram(lwb, new Range(0, 1));
}