protected override Image Load(string filename)
{
if (File.Exists(filename))
{
ecw_wrapper_init();
var view = FileView.Open(filename);
view.Set();
var info = view.Info;
int width = (int) (info.SizeX - 1);
int height = (int) (info.SizeY - 1);
var image = NewImage(width, height, ImageBaseType.Rgb,
ImageType.Rgb, filename);
var rgn = new PixelRgn(image.Layers[0], true, false);
var line = new byte[width * 3];
for (int y = 0; y < height; y++)
{
view.ReadLineRGB(line);
rgn.SetRow(line, 0, y, width);
}
view.Close();
return image;
}
return null;
}
override protected Image Load(string filename)
{
if (File.Exists(filename))
{
ecw_wrapper_init();
var view = FileView.Open(filename);
view.Set();
var info = view.Info;
int width = (int)(info.SizeX - 1);
int height = (int)(info.SizeY - 1);
var image = NewImage(width, height, ImageBaseType.Rgb,
ImageType.Rgb, filename);
var rgn = new PixelRgn(image.Layers[0], true, false);
var line = new byte[width * 3];
for (int y = 0; y < height; y++)
{
view.ReadLineRGB(line);
rgn.SetRow(line, 0, y, width);
}
view.Close();
return(image);
}
return(null);
}
override protected Image Load()
{
if (ReadHeader())
{
var colormap = new RGB[] {
new RGB(255, 255, 255),
new RGB(255, 0, 0),
new RGB(0, 255, 0),
new RGB(255, 255, 0),
new RGB(0, 0, 255),
new RGB(255, 0, 255),
new RGB(0, 255, 255),
new RGB(181, 181, 181),
new RGB(84, 84, 84),
new RGB(127, 0, 0),
new RGB(0, 127, 0),
new RGB(127, 127, 0),
new RGB(0, 0, 127),
new RGB(127, 0, 127),
new RGB(0, 127, 127),
new RGB(0, 0, 0)
};
var image = NewImage(_imageWidth, _imageHeight,
ImageBaseType.Indexed,
ImageType.Indexed, Filename);
image.Colormap = colormap;
var rgn = new PixelRgn(image.Layers[0], true, false);
int bparrow = (_imageWidth + 7) / 8;
for (int y = 0; y < _imageHeight;)
{
// byte[] line = new byte[bparrow * 8];
var line = new byte[_imageWidth];
int count = ReadLine(line);
do
{
rgn.SetRow(line, 0, y++);
}while (--count > 0);
}
return(image);
}
return(null);
}
override protected void Render(Image image, Drawable drawable)
{
var progress = new Progress(_("Colorizing..."));
var rectangle = drawable.MaskIntersect;
// Fix me: replace with x1, y1, x2, y2
int i = rectangle.X1;
int j = rectangle.Y1;
int ii = rectangle.X2;
int jj = rectangle.Y2;
bool hasSel = (rectangle != null);
if (!hasSel || _useEntireImage)
{
j = i = 0;
jj = image.Width;
ii = image.Height;
}
Drawable sel = null;
int threshGuc = (int)(_thresh * 255);
PixelRgn selRgn = null;
if (hasSel)
{
sel = image.Selection;
selRgn = new PixelRgn(sel, rectangle, false, false);
}
var srcRgn = new PixelRgn(drawable, rectangle, false, false);
var dstRgn = new PixelRgn(drawable, rectangle, true, true);
var markRgn = new PixelRgn(_marked, rectangle, false, false);
int h = srcRgn.H;
int w = srcRgn.W;
var A = new double[WindowPixels, h *w];
var AI = new int[WindowPixels * h * w];
var AJ = new int[WindowPixels * h * w];
var Y = new double[h, w];
var I = new double[h, w];
var Q = new double[h, w];
double[,] inI = null;
double[,] inQ = null;
if (_useChroma)
{
inI = new double[h, w];
inQ = new double[h, w];
}
bool[,] mask = new bool[h, w];
Tile.CacheDefault(drawable);
if (sel != null)
{
// Retarded check for selections, because gimp doesn't
// _REALLY_ return FALSE when there's no selection.
if (j == 0 && i == 0 && jj == drawable.Width &&
ii == drawable.Height)
{
bool goodSelection = false;
foreach (var pixel in
new ReadPixelIterator(sel, RunMode.Noninteractive))
{
if (pixel[0] != 0)
{
goodSelection = true;
break;
}
}
if (!goodSelection)
{
sel.Detach();
sel = null;
}
}
}
Pixel[] selRow = null;
var whitePixel = new Pixel(255, 255, 255);
for (i = 0; i < h; i++)
{
var imgRow = srcRgn.GetRow(srcRgn.X, srcRgn.Y + i, w);
var markRow = markRgn.GetRow(markRgn.X, markRgn.Y + i, w);
if (sel != null)
{
selRow = selRgn.GetRow(selRgn.X, selRgn.Y + i, w);
}
for (j = 0; j < w; j++)
{
var imgPixel = imgRow[j];
var markPixel = markRow[j];
int selIdx = (sel != null) ? j : 0;
double iY, iI, iQ;
double mY;
rgb2yiq(imgPixel, out iY, out iI, out iQ);
if (_useChroma)
{
inI[i, j] = iI;
inQ[i, j] = iQ;
}
if (_includeOriginal)
{
var diff = imgPixel - markPixel;;
int delta = Math.Abs(diff.Red) + Math.Abs(diff.Green) +
Math.Abs(diff.Blue);
}
// big dirty if statement
if (_pureWhite && markPixel.IsSameColor(whitePixel))
{
mask[i, j] = true;
}
else if ((_includeOriginal && !imgPixel.IsSameColor(markPixel)) ||
(!_includeOriginal && markPixel.Alpha >= threshGuc))
{
mask[i, j] = true;
rgb2yiq(markPixel, out mY, out iI, out iQ);
}
else if (sel != null && selRow[selIdx].Red < threshGuc)
{
mask[i, j] = true;
}
else
{
mask[i, j] = false;
iI = iQ = 0;
}
Y[i, j] = iY;
I[i, j] = iI;
Q[i, j] = iQ;
}
}
if (sel != null)
{
sel.Detach();
}
progress.Update(0.1);
int n = 0;
for (i = 0; i < h; i++)
{
for (j = 0; j < w; j++)
{
if (!mask[i, j])
{
int min_ii = Math.Max(0, i - WindowRadius);
int max_ii = Math.Min(h - 1, i + WindowRadius);
int min_jj = Math.Max(0, j - WindowRadius);
int max_jj = Math.Min(w - 1, j + WindowRadius);
var vary = new int[WindowPixels];
var varx = new int[WindowPixels];
var var = new double[WindowPixels];
int count = 0;
double sum_sq = 0;
double sum = 0;
double min_variance = 1.0;
for (ii = min_ii; ii <= max_ii; ii++)
{
for (jj = min_jj; jj <= max_jj; jj++)
{
double val = Y[ii, jj];
sum += val;
sum_sq += val * val;
if (ii == i && jj == j)
{
continue;
}
vary[count] = i * w + j;
varx[count] = ii * w + jj;
var[count] = val - Y[i, j];
var[count] *= var[count];
if (_useChroma)
{
val = inI[ii, jj] - inI[i, j];
var[count] += val * val;
val = inQ[ii, jj] - inQ[i, j];
var[count] += val * val;
}
if (var[count] < min_variance)
{
min_variance = var[count];
}
++count;
}
}
double sigma =
(sum_sq - (sum * sum) / (double)(count + 1)) / (double)count;
if (sigma < 0.000002)
{
sigma = 0.000002;
}
else if (sigma < (min_variance / LN_100))
{
sigma = min_variance / LN_100;
}
sum = 0;
for (ii = 0; ii < count; ii++)
{
var[ii] = Math.Exp(-var[ii] / sigma);
sum += var[ii];
}
for (ii = 0; ii < count; ii++)
{
AI[n] = vary[ii];
AJ[n] = varx[ii];
// Fix me: just A[i, j]?
A[n / (h * w), n % (h * w)] = -var[ii] / sum;
++n;
}
}
AI[n] = AJ[n] = i * w + j;
// Fix me: just A[i, j]?
A[n / (h * w), n % (h * w)] = 1.0;
++n;
}
}
UmfPack umf = new UmfPack();
umf.Defaults();
var Ax = new double[WindowPixels, h *w];
var Ap = new int[h * w + 1];
var Ai = new int[WindowPixels * h * w];
var Map = new int[WindowPixels * h * w];
umf.TripletToCol(h * w, h * w, n, AI, AJ, A, Ap, Ai, Ax, Map);
umf.Symbolic(h * w, h * w, Ap, Ai, Ax);
umf.Numeric(Ap, Ai, Ax);
umf.FreeSymbolic();
progress.Update(0.3);
var outI = new double[h, w];
var outQ = new double[h, w];
umf.Solve(Ap, Ai, Ax, outI, I);
progress.Update(0.6);
umf.Solve(Ap, Ai, Ax, outQ, Q);
umf.FreeNumeric();
progress.Update(0.9);
for (i = 0; i < h; i++)
{
// FIXME: This is only for the alpha channel..
var imgRow = srcRgn.GetRow(srcRgn.X, srcRgn.Y + i, w);
for (j = 0; j < w; j++)
{
yiq2rgb(Y[i, j], outI[i, j], outQ[i, j], imgRow[j]);
}
dstRgn.SetRow(imgRow, dstRgn.X, dstRgn.Y + i);
}
drawable.Flush();
drawable.MergeShadow(true);
drawable.Update(dstRgn.X, dstRgn.Y, dstRgn.W, dstRgn.H);
progress.Update(1.0);
}
protected override void Render(Image image, Drawable drawable)
{
var progress = new Progress(_("Colorizing..."));
var rectangle = drawable.MaskIntersect;
// Fix me: replace with x1, y1, x2, y2
int i = rectangle.X1;
int j = rectangle.Y1;
int ii = rectangle.X2;
int jj = rectangle.Y2;
bool hasSel = (rectangle != null);
if (!hasSel || _useEntireImage)
{
j = i = 0;
jj = image.Width;
ii = image.Height;
}
Drawable sel = null;
int threshGuc = (int) (_thresh * 255);
PixelRgn selRgn = null;
if (hasSel)
{
sel = image.Selection;
selRgn = new PixelRgn(sel, rectangle, false, false);
}
var srcRgn = new PixelRgn(drawable, rectangle, false, false);
var dstRgn = new PixelRgn(drawable, rectangle, true, true);
var markRgn = new PixelRgn(_marked, rectangle, false, false);
int h = srcRgn.H;
int w = srcRgn.W;
var A = new double[WindowPixels, h * w];
var AI = new int[WindowPixels * h * w];
var AJ = new int[WindowPixels * h * w];
var Y = new double[h, w];
var I = new double[h, w];
var Q = new double[h, w];
double[,] inI = null;
double[,] inQ = null;
if (_useChroma)
{
inI = new double[h, w];
inQ = new double[h, w];
}
bool[,] mask = new bool[h, w];
Tile.CacheDefault(drawable);
if (sel != null)
{
// Retarded check for selections, because gimp doesn't
// _REALLY_ return FALSE when there's no selection.
if (j == 0 && i == 0 && jj == drawable.Width &&
ii == drawable.Height)
{
bool goodSelection = false;
foreach (var pixel in
new ReadPixelIterator(sel, RunMode.Noninteractive))
{
if (pixel[0] != 0)
{
goodSelection = true;
break;
}
}
if (!goodSelection)
{
sel.Detach();
sel = null;
}
}
}
Pixel[] selRow = null;
var whitePixel = new Pixel(255, 255, 255);
for (i = 0; i < h; i++)
{
var imgRow = srcRgn.GetRow(srcRgn.X, srcRgn.Y + i, w);
var markRow = markRgn.GetRow(markRgn.X, markRgn.Y + i, w);
if (sel != null)
{
selRow = selRgn.GetRow(selRgn.X, selRgn.Y + i, w);
}
for (j = 0; j < w; j++)
{
var imgPixel = imgRow[j];
var markPixel = markRow[j];
int selIdx = (sel != null) ? j : 0;
double iY, iI, iQ;
double mY;
rgb2yiq(imgPixel, out iY, out iI, out iQ);
if (_useChroma)
{
inI[i, j] = iI;
inQ[i, j] = iQ;
}
if (_includeOriginal)
{
var diff = imgPixel - markPixel;;
int delta = Math.Abs(diff.Red) + Math.Abs(diff.Green) +
Math.Abs(diff.Blue);
}
// big dirty if statement
if (_pureWhite && markPixel.IsSameColor(whitePixel))
{
mask[i, j] = true;
}
else if ((_includeOriginal && !imgPixel.IsSameColor(markPixel))
|| (!_includeOriginal && markPixel.Alpha >= threshGuc))
{
mask[i, j] = true;
rgb2yiq(markPixel, out mY, out iI, out iQ);
}
else if (sel != null && selRow[selIdx].Red < threshGuc)
{
mask[i, j] = true;
}
else {
mask[i, j] = false;
iI = iQ = 0;
}
Y[i, j] = iY;
I[i, j] = iI;
Q[i, j] = iQ;
}
}
if (sel != null)
{
sel.Detach();
}
progress.Update(0.1);
int n = 0;
for (i = 0; i < h; i++)
{
for (j = 0; j < w; j++)
{
if (!mask[i, j])
{
int min_ii = Math.Max(0, i - WindowRadius);
int max_ii = Math.Min(h - 1, i + WindowRadius);
int min_jj = Math.Max(0, j - WindowRadius);
int max_jj = Math.Min(w - 1, j + WindowRadius);
var vary = new int[WindowPixels];
var varx = new int[WindowPixels];
var var = new double[WindowPixels];
int count = 0;
double sum_sq = 0;
double sum = 0;
double min_variance = 1.0;
for (ii = min_ii; ii <= max_ii; ii++)
{
for (jj = min_jj; jj <= max_jj; jj++)
{
double val = Y[ii, jj];
sum += val;
sum_sq += val * val;
if (ii == i && jj == j)
continue;
vary[count] = i * w + j;
varx[count] = ii * w + jj;
var[count] = val - Y[i, j];
var[count] *= var[count];
if (_useChroma)
{
val = inI[ii, jj] - inI[i, j];
var[count] += val * val;
val = inQ[ii, jj] - inQ[i, j];
var[count] += val * val;
}
if (var[count] < min_variance)
min_variance = var[count];
++count;
}
}
double sigma =
(sum_sq - (sum * sum)/(double)(count + 1)) / (double)count;
if (sigma < 0.000002)
sigma = 0.000002;
else if (sigma < (min_variance / LN_100))
sigma = min_variance / LN_100;
sum = 0;
for (ii = 0; ii < count; ii++)
{
var[ii] = Math.Exp(-var[ii] / sigma);
sum += var[ii];
}
for (ii = 0; ii < count; ii++)
{
AI[n] = vary[ii];
AJ[n] = varx[ii];
// Fix me: just A[i, j]?
A[n / (h * w) , n % (h * w)] = -var[ii] / sum;
++n;
}
}
AI[n] = AJ[n] = i * w + j;
// Fix me: just A[i, j]?
A[n / (h * w), n % (h * w)] = 1.0;
++n;
}
}
UmfPack umf = new UmfPack();
umf.Defaults();
var Ax = new double[WindowPixels, h * w];
var Ap = new int[h * w + 1];
var Ai = new int[WindowPixels * h * w];
var Map = new int[WindowPixels * h * w];
umf.TripletToCol(h * w, h * w, n, AI, AJ, A, Ap, Ai, Ax, Map);
umf.Symbolic(h * w, h * w, Ap, Ai, Ax);
umf.Numeric(Ap, Ai, Ax);
umf.FreeSymbolic();
progress.Update(0.3);
var outI = new double[h, w];
var outQ = new double[h, w];
umf.Solve(Ap, Ai, Ax, outI, I);
progress.Update(0.6);
umf.Solve(Ap, Ai, Ax, outQ, Q);
umf.FreeNumeric();
progress.Update(0.9);
for (i = 0; i < h; i++)
{
// FIXME: This is only for the alpha channel..
var imgRow = srcRgn.GetRow(srcRgn.X, srcRgn.Y + i, w);
for (j = 0; j < w; j++)
{
yiq2rgb(Y[i, j], outI[i, j], outQ[i, j], imgRow[j]);
}
dstRgn.SetRow(imgRow, dstRgn.X, dstRgn.Y + i);
}
drawable.Flush();
drawable.MergeShadow(true);
drawable.Update(dstRgn.X, dstRgn.Y, dstRgn.W, dstRgn.H);
progress.Update(1.0);
}
protected override Image Load()
{
if (ReadHeader())
{
var colormap = new RGB[] {
new RGB(255,255,255),
new RGB(255,0,0),
new RGB(0,255,0),
new RGB(255,255,0),
new RGB(0,0,255),
new RGB(255,0,255),
new RGB(0,255,255),
new RGB(181,181,181),
new RGB(84,84,84),
new RGB(127,0,0),
new RGB(0,127,0),
new RGB(127,127,0),
new RGB(0,0,127),
new RGB(127,0,127),
new RGB(0,127,127),
new RGB(0,0,0)
};
var image = NewImage(_imageWidth, _imageHeight,
ImageBaseType.Indexed,
ImageType.Indexed, Filename);
image.Colormap = colormap;
var rgn = new PixelRgn(image.Layers[0], true, false);
int bparrow = (_imageWidth + 7) / 8;
for (int y = 0; y < _imageHeight; )
{
// byte[] line = new byte[bparrow * 8];
var line = new byte[_imageWidth];
int count = ReadLine(line);
do
{
rgn.SetRow(line, 0, y++);
}
while (--count > 0);
}
return image;
}
return null;
}
