X-Git-Url: http://git.euphorik.ch/?p=master-thesis.git;a=blobdiff_plain;f=Parasitemia%2FParasitemia%2FKMedians.fs;fp=Parasitemia%2FParasitemia%2FKMedians.fs;h=82e09cb25e6192286fc8ad5609c84e46458cae8e;hp=09d6c8f70d525365721fe3390540e07c94c787c4;hb=e76da913cd58078ad2479357b2430ed62a6e0777;hpb=d9a6e072ecf299db691c05bb559a71265f812ba3 diff --git a/Parasitemia/Parasitemia/KMedians.fs b/Parasitemia/Parasitemia/KMedians.fs index 09d6c8f..82e09cb 100644 --- a/Parasitemia/Parasitemia/KMedians.fs +++ b/Parasitemia/Parasitemia/KMedians.fs @@ -22,28 +22,28 @@ let kmedians (img: Image) (fgFactor: float) : Result = let max = ref [| 0.0 |] let maxLocation = ref <| [| Point() |] img.MinMax(min, max, minLocation, maxLocation) - + let mutable median_bg = (!max).[0] - ((!max).[0] - (!min).[0]) / 4.0 let mutable median_fg = (!min).[0] + ((!max).[0] - (!min).[0]) / 4.0 let mutable d_bg = new Image(img.Size) let mutable d_fg = new Image(img.Size) let mutable fg = new Image(img.Size) - + for i in 1..nbIteration do CvInvoke.Pow(img - median_bg, 2.0, d_bg) CvInvoke.Pow(img - median_fg, 2.0, d_fg) fg <- (d_fg * fgFactor).Cmp(d_bg, CvEnum.CmpType.LessThan) - + median_fg <- MathNet.Numerics.Statistics.Statistics.Median(seq { for i in 0 .. h - 1 do for j in 0 .. w - 1 do if fg.Data.[i, j, 0] > 0uy then yield img.Data.[i, j, 0] |> float }) - + median_bg <- MathNet.Numerics.Statistics.Statistics.Median(seq { for i in 0 .. h - 1 do for j in 0 .. w - 1 do if fg.Data.[i, j, 0] = 0uy then yield img.Data.[i, j, 0] |> float }) - + CvInvoke.Sqrt(d_fg, d_fg) { fg = fg; median_bg = median_bg; median_fg = median_fg; d_fg = d_fg }