public static Image <Bgr, float> Bgr2Hsv(Image <Gray, float> grayimage)
{
var numPoints = grayimage.Rows;
var b = new Image <Gray, float>(1, numPoints);
var g = new Image <Gray, float>(1, numPoints);
var r = new Image <Gray, float>(1, numPoints);
MatHelper.CopyTo(ref grayimage, 0, 0, 1, numPoints,
ref b, 0, 0, 1, numPoints);
MatHelper.CopyTo(ref grayimage, 1, 0, 1, numPoints,
ref g, 0, 0, 1, numPoints);
MatHelper.CopyTo(ref grayimage, 2, 0, 1, numPoints,
ref r, 0, 0, 1, numPoints);
return(Bgr2Hsv(b, g, r));
}
/// <summary>
/// function rgb=hsl2rgb(hsl_in)
/// </summary>
/// <param name="bgrimage"></param>
/// <returns></returns>
public static Image <Bgr, byte> Hsl2Bgr(Image <Bgr, byte> bgrimage)
{
var doubleBgr = bgrimage.Convert <Bgr, float>() / 255;
var numPoints = bgrimage.Rows * bgrimage.Cols;
var H = new Image <Gray, float>(1, numPoints);
var S = new Image <Gray, float>(1, numPoints);
var L = new Image <Gray, float>(1, numPoints);
var channels = doubleBgr.Split();
MatHelper.CopyTo(ref channels[2], 0, 0, channels[2].Cols, channels[0].Rows,
ref H, 0, 0, 1, numPoints);
MatHelper.CopyTo(ref channels[1], 0, 0, channels[1].Cols, channels[1].Rows,
ref S, 0, 0, 1, numPoints);
MatHelper.CopyTo(ref channels[0], 0, 0, channels[0].Cols, channels[2].Rows,
ref L, 0, 0, 1, numPoints);
Image <Gray, float> auxImage = new Image <Gray, float>(L.Size);
Image <Gray, Byte> auxImageByte = new Image <Gray, Byte>(L.Size);
auxImage.SetValue(new Gray(0.5));
Mat lowlidx = new Mat(L.Rows, L.Cols, Emgu.CV.CvEnum.DepthType.Cv8U, L.NumberOfChannels);
Mat nonlowlidx = new Mat(L.Rows, L.Cols, Emgu.CV.CvEnum.DepthType.Cv8U, L.NumberOfChannels);
CvInvoke.Compare(L, auxImage, lowlidx, Emgu.CV.CvEnum.CmpType.LessEqual);
CvInvoke.BitwiseNot(lowlidx, nonlowlidx);
//var lowLidx = L.ThresholdBinaryInv(new Gray(1.0 / 2.0), new Gray(1.0));
//Mat q = new Mat(L.Rows, L.Cols, , L.NumberOfChannels);
//Mat p = new Mat(L.Rows, L.Cols, Emgu.CV.CvEnum.DepthType.Cv8U, L.NumberOfChannels);
var partial1 = 1 + S;
partial1._Mul(L);
auxImageByte.SetValue(new Gray(255));
CvInvoke.Divide(lowlidx, auxImageByte, lowlidx);
Mat lowlidxFloat = new Mat();
lowlidx.ConvertTo(lowlidxFloat, Emgu.CV.CvEnum.DepthType.Cv32F);
CvInvoke.Multiply(partial1.Mat, lowlidxFloat, partial1.Mat);
CvInvoke.Divide(nonlowlidx, auxImageByte, nonlowlidx);
var partial2 = L.Mul(S);
Mat nonlowlidxFloat = new Mat();
nonlowlidx.ConvertTo(nonlowlidxFloat, Emgu.CV.CvEnum.DepthType.Cv32F);
CvInvoke.Multiply(partial2.Mat, nonlowlidxFloat, partial2.Mat);
partial2 = L + S - partial2;
var q = partial1 + partial2;
var p = 2 * L - q;
var hk = H;
var t1 = hk + 1 / 3.0;
var t2 = hk;
var t3 = hk - 1 / 3.0;
var underidx1 = t1.ThresholdBinaryInv(new Gray(0), new Gray(1));
var underidx2 = t2.ThresholdBinaryInv(new Gray(0), new Gray(1));
var underidx3 = t3.ThresholdBinaryInv(new Gray(0), new Gray(1));
var overidx1 = t1.ThresholdBinary(new Gray(1), new Gray(1));
var overidx2 = t2.ThresholdBinary(new Gray(1), new Gray(1));
var overidx3 = t3.ThresholdBinary(new Gray(1), new Gray(1));
var t = new Image <Bgr, float>(new Image <Gray, float>[] { t1, t2, t3 });
var underidx = new Image <Bgr, float>(new Image <Gray, float>[] { underidx2, underidx2, underidx3 });
var overidx = new Image <Bgr, float>(new Image <Gray, float>[] { overidx1, overidx2, overidx3 });
t = t + underidx - overidx;
t1.Dispose();
t2.Dispose();
t3.Dispose();
underidx1.Dispose();
underidx2.Dispose();
underidx3.Dispose();
overidx1.Dispose();
overidx2.Dispose();
overidx3.Dispose();
var tSplit = t.Split();
t1 = tSplit[0];
t2 = tSplit[1];
t3 = tSplit[2];
var range1_1 = t1.ThresholdBinaryInv(new Gray(1.0 / 6), new Gray(1.0));
var range1_2 = t2.ThresholdBinaryInv(new Gray(1.0 / 6), new Gray(1.0));
var range1_3 = t3.ThresholdBinaryInv(new Gray(1.0 / 6), new Gray(1.0));
var range2_1 = t1.InRange(new Gray(1.0 / 6), new Gray(1.0 / 2)).Convert <Gray, float>() / 255.0;
var range2_2 = t2.InRange(new Gray(1.0 / 6), new Gray(1.0 / 2)).Convert <Gray, float>() / 255.0;
var range2_3 = t3.InRange(new Gray(1.0 / 6), new Gray(1.0 / 2)).Convert <Gray, float>() / 255.0;
var range3_1 = t1.InRange(new Gray(1.0 / 2), new Gray(2.0 / 3)).Convert <Gray, float>() / 255.0;
var range3_2 = t2.InRange(new Gray(1.0 / 2), new Gray(2.0 / 3)).Convert <Gray, float>() / 255.0;
var range3_3 = t3.InRange(new Gray(1.0 / 2), new Gray(2.0 / 3)).Convert <Gray, float>() / 255.0;
var range4_1 = t1.ThresholdBinary(new Gray(2.0 / 3), new Gray(1.0));
var range4_2 = t2.ThresholdBinary(new Gray(2.0 / 3), new Gray(1.0));
var range4_3 = t3.ThresholdBinary(new Gray(2.0 / 3), new Gray(1.0));
var range1 = new Image <Bgr, float>(new Image <Gray, float>[] { range1_1, range1_2, range1_3 });
var range2 = new Image <Bgr, float>(new Image <Gray, float>[] { range2_1, range2_2, range2_3 });
var range3 = new Image <Bgr, float>(new Image <Gray, float>[] { range3_1, range3_2, range3_3 });
var range4 = new Image <Bgr, float>(new Image <Gray, float>[] { range4_1, range4_2, range4_3 });
underidx.Dispose();
underidx1.Dispose();
underidx2.Dispose();
underidx3.Dispose();
overidx.Dispose();
overidx1.Dispose();
overidx2.Dispose();
overidx3.Dispose();
t1.Dispose();
t2.Dispose();
t3.Dispose();
range1_1.Dispose();
range1_2.Dispose();
range1_3.Dispose();
range2_1.Dispose();
range2_2.Dispose();
range2_3.Dispose();
range3_1.Dispose();
range3_2.Dispose();
range3_3.Dispose();
range4_1.Dispose();
range4_2.Dispose();
range4_3.Dispose();
var P = new Image <Bgr, float>(new Image <Gray, float>[] { p, p, p });
var Q = new Image <Bgr, float>(new Image <Gray, float>[] { q, q, q });
var rgb_c = (P + ((6 * (Q - P)).Mul(t))).Mul(range1) +
Q.Mul(range2) +
(P + (((Q - P) * 6).Mul(2 / 3 - t))).Mul(range3) +
P.Mul(range4);
rgb_c *= 255.0;
P.Dispose();
Q.Dispose();
p.Dispose();
q.Dispose();
t.Dispose();
var rgb_cShape = new Image <Bgr, float>(bgrimage.Cols, bgrimage.Rows);
MatHelper.CopyTo(ref rgb_c, 0, 0, rgb_c.Cols, rgb_c.Rows, ref rgb_cShape, 0, 0, rgb_cShape.Cols, rgb_cShape.Rows);
//var H = hsv.Split()[0];
////MatHelper.CopyTo(ref hsplit0, 0, 0, hsplit0.Cols, hsplit0.Rows, ref H, 0, 0, H.Cols, H.Rows);
//var L_shape = new Image<Gray, float>(bgrimage.Cols, bgrimage.Rows);
//var S_shape = new Image<Gray, float>(bgrimage.Cols, bgrimage.Rows);
//MatHelper.CopyTo(ref S, 0, 0, S.Cols, S.Rows, ref S_shape, 0, 0, S_shape.Cols, S_shape.Rows);
//MatHelper.CopyTo(ref L, 0, 0, L.Cols, L.Rows, ref L_shape, 0, 0, L_shape.Cols, L_shape.Rows);
return(rgb_cShape.Convert <Bgr, byte>());
}
public static Image <Bgr, byte> Bgr2Hsl(Image <Bgr, byte> bgrimage)
{
var doubleBgr = bgrimage.Convert <Bgr, float>() / 255;
var numPoints = bgrimage.Rows * bgrimage.Cols;
var bgrDoubleReshaped = new Image <Gray, float>(3, numPoints);
var channels = doubleBgr.Split();
MatHelper.CopyTo(ref channels[0], 0, 0, channels[0].Cols, channels[0].Rows,
ref bgrDoubleReshaped, 0, 0, 1, numPoints);
MatHelper.CopyTo(ref channels[1], 0, 0, channels[1].Cols, channels[1].Rows,
ref bgrDoubleReshaped, 1, 0, 1, numPoints);
MatHelper.CopyTo(ref channels[2], 0, 0, channels[2].Cols, channels[2].Rows,
ref bgrDoubleReshaped, 2, 0, 1, numPoints);
var mx = MatHelper.maxInLines(bgrDoubleReshaped);
var mn = MatHelper.minInLines(bgrDoubleReshaped);
Mat nonzeroidx = new Mat(mx.Rows, mx.Cols, Emgu.CV.CvEnum.DepthType.Cv8U, mx.NumberOfChannels);
Mat zeroidx = new Mat(mx.Rows, mx.Cols, Emgu.CV.CvEnum.DepthType.Cv8U, mx.NumberOfChannels);
var L = (mx + mn) / 2;
var S = new Image <Gray, float>(L.Size);
CvInvoke.Compare(mx, mn, zeroidx, Emgu.CV.CvEnum.CmpType.Equal);
CvInvoke.BitwiseNot(zeroidx, nonzeroidx);
Mat lowlidx = new Mat(mx.Rows, mx.Cols, Emgu.CV.CvEnum.DepthType.Cv8U, mx.NumberOfChannels);
var calc = new Image <Gray, float>(mx.Size);
var mmsum = mx + mn;
var mmsub = mx - mn;
var auxImage = new Image <Gray, float>(L.Size);
auxImage.SetValue(new Gray(0.5));
////CvInvoke.Threshold(L, lowlidx, 0.5, 1.0, Emgu.CV.CvEnum.ThresholdType.BinaryInv);
CvInvoke.Compare(L, auxImage, lowlidx, Emgu.CV.CvEnum.CmpType.LessEqual);
//CvInvoke.Divide(mmsub, mmsum, calc);
calc = mmsub.Mul(1.0 / mmsum);
Mat idx = new Mat();
CvInvoke.BitwiseAnd(lowlidx, nonzeroidx, idx);
calc.Mat.CopyTo(S, idx);
Mat hilidx = new Mat(mx.Rows, mx.Cols, Emgu.CV.CvEnum.DepthType.Cv8U, mx.NumberOfChannels);
CvInvoke.Compare(L, auxImage, hilidx, Emgu.CV.CvEnum.CmpType.GreaterThan);
var mmsum2 = (2 - (mx + mn));
calc = mmsub.Mul(1.0 / mmsum2);
CvInvoke.BitwiseAnd(hilidx, nonzeroidx, idx);
calc.Mat.CopyTo(S, idx);
//var hsv = doubleBgr.Convert<Hsv, float>() / 255;
var hsv = Bgr2Hsv(bgrDoubleReshaped);
var H = hsv.Split()[2];
//MatHelper.CopyTo(ref hsplit0, 0, 0, hsplit0.Cols, hsplit0.Rows, ref H, 0, 0, H.Cols, H.Rows);
var L_shape = new Image <Gray, float>(bgrimage.Cols, bgrimage.Rows);
var S_shape = new Image <Gray, float>(bgrimage.Cols, bgrimage.Rows);
var H_shape = new Image <Gray, float>(bgrimage.Cols, bgrimage.Rows);
MatHelper.CopyTo(ref S, 0, 0, S.Cols, S.Rows, ref S_shape, 0, 0, S_shape.Cols, S_shape.Rows);
MatHelper.CopyTo(ref L, 0, 0, L.Cols, L.Rows, ref L_shape, 0, 0, L_shape.Cols, L_shape.Rows);
MatHelper.CopyTo(ref H, 0, 0, H.Cols, H.Rows, ref H_shape, 0, 0, H_shape.Cols, H_shape.Rows);
CvInvoke.Transpose(H_shape, H_shape);
CvInvoke.Transpose(S_shape, S_shape);
CvInvoke.Transpose(L_shape, L_shape);
var hsl = new Image <Bgr, float>(new Image <Gray, float>[] { L_shape, S_shape, H_shape });
hsl *= 255;
return(hsl.Convert <Bgr, byte>());
}
