X-Git-Url: http://git.euphorik.ch/?p=master-thesis.git;a=blobdiff_plain;f=Parasitemia%2FParasitemia%2FParasitesMarker2.fs;fp=Parasitemia%2FParasitemia%2FParasitesMarker2.fs;h=0000000000000000000000000000000000000000;hp=42645b4899f6a6ae55d4d7a02bb027d2dca1a65d;hb=f9ee9f2c0c1194a6a8f7475ee2ca0b7076b32261;hpb=785a7ea4885e87d76f0b1fb1c715230a08389558 diff --git a/Parasitemia/Parasitemia/ParasitesMarker2.fs b/Parasitemia/Parasitemia/ParasitesMarker2.fs deleted file mode 100644 index 42645b4..0000000 --- a/Parasitemia/Parasitemia/ParasitesMarker2.fs +++ /dev/null @@ -1,69 +0,0 @@ -module ParasitesMarker2 - -open System.Drawing -open System.Linq - -open Emgu.CV -open Emgu.CV.Structure - -open Utils - -type Result = { - darkStain: Image - infection: Image - stain: Image } - -// Create three binary markers : -// * 'Dark stain' corresponds to the colored pixel, it's independent of the size of the areas. -// * 'Stain' corresponds to the stain around the parasites. -// * 'Infection' corresponds to the parasite. It shouldn't contain thrombocytes. -let find (filteredGreen: Image) (filteredGreenFloat: Image) (config: Config.Config) : Result * Image * Image = - - let filteredGreenWithoutInfection = filteredGreen.Copy() - ImgTools.areaClose filteredGreenWithoutInfection (int config.InfectionArea) - - let filteredGreenWithoutStain = filteredGreenWithoutInfection.Copy() - ImgTools.areaClose filteredGreenWithoutStain (int config.StainArea) - - // We use the filtered image to find the dark stain. - let kmediansResults = logTime "Finding fg/bg (k-medians)" (fun () -> KMedians.kmedians (filteredGreenWithoutInfection.Convert()) 1.0) // FIXME: avoid converting again this in MainAnalysis - let { KMedians.fg = fg; KMedians.median_bg = median_bg; KMedians.median_fg = median_fg; KMedians.d_fg = d_fg } = kmediansResults - let darkStain = d_fg.Cmp(median_bg * config.Parameters.darkStainLevel, CvEnum.CmpType.GreaterThan) - darkStain._And(filteredGreenWithoutInfection.Cmp(median_fg, CvEnum.CmpType.LessThan)) - - let marker (img: Image) (closed: Image) (level: float) : Image = - let diff = closed - img - - let min = ref [| 0. |] - let minLocation = ref <| [| Point() |] - let max = ref [| 0. |] - let maxLocation = ref <| [| Point() |] - diff.MinMax(min, max, minLocation, maxLocation) - let max = (!max).[0] - - let threshold = 0.2 * max - let diff' = diff - threshold - - let m = MathNet.Numerics.Statistics.StreamingStatistics.Mean(seq { - for m in ImgTools.findMaxima diff' do - let p = m.First() - yield diff'.Data.[p.Y, p.X, 0] |> float }) + threshold - - let valueThreshold = m * level - - diff._ThresholdBinary(Gray(valueThreshold), Gray(255.)) - diff - - let infectionMarker = marker filteredGreen filteredGreenWithoutInfection config.Parameters.infectionLevel - let stainMarker = marker filteredGreenWithoutInfection filteredGreenWithoutStain config.Parameters.stainLevel - - { darkStain = darkStain - infection = infectionMarker - stain = stainMarker }, - filteredGreenWithoutInfection, - filteredGreenWithoutStain - - - - -