X-Git-Url: http://git.euphorik.ch/?p=master-thesis.git;a=blobdiff_plain;f=Parasitemia%2FParasitemia%2FProgram.fs;h=09f43126106c852a7f65209546448e8db679f893;hp=ac495b591819bcdb2ab086f441477cf59935d5ef;hb=044b0ae69df3ac565432545b2fa934589016f9bd;hpb=6d99ab0dadc57e91cd846844b620fc90b52a2f96 diff --git a/Parasitemia/Parasitemia/Program.fs b/Parasitemia/Parasitemia/Program.fs index ac495b5..09f4312 100644 --- a/Parasitemia/Parasitemia/Program.fs +++ b/Parasitemia/Parasitemia/Program.fs @@ -9,6 +9,9 @@ open System.Windows.Shapes open System.Windows.Controls open System.Drawing open System.Diagnostics +open System.Threading + +open FSharp.Collections.ParallelSeq open Emgu.CV open Emgu.CV.Structure @@ -61,26 +64,27 @@ let main args = Config( { scale = 1. - minRbcRadius = -0.3 - maxRbcRadius = 0.3 + initialAreaOpen = 2000 - preFilterSigma = 1.5 + minRbcRadius = -0.32f + maxRbcRadius = 0.32f - factorNbPick = 1.0 - factorWindowSize = 2.0 + preFilterSigma = 1.7 // 1.5 - darkStainLevel = 0.4 // Lower -> more sensitive. + factorNbPick = 1.0 - infectionArea = 0.015 - infectionLevel = 0.2 - parasitePixelsRequired = 1 + darkStainLevel = 0.22 // Lower -> more sensitive. 0.3. Careful about illumination on the borders. + maxDarkStainRatio = 0.1 // 10 % - stainArea = 0.06 - stainLevel = 0.2 + infectionArea = 0.012f // 1.2 % + infectionLevel = 1.12 // Lower -> more sensitive. - maxDarkStainRatio = 0.1 + stainArea = 0.08f // 8 % + stainLevel = 1.1 // Lower -> more sensitive. + maxStainRatio = 0.12 // 12 % - minimumCellArea = 0.3 }) + standardDeviationMaxRatio = 0.5 // 0.55 + minimumCellArea = 0.5f }) match mode with | CmdLine (input, output) -> @@ -100,16 +104,25 @@ let main args = use resultFile = new StreamWriter(new FileStream(Path.Combine(output, "results.txt"), FileMode.Append, FileAccess.Write)) - for file in files do - try - use img = new Image(file) - Utils.log (sprintf "== File: %A" file) + //try + let images = seq { for file in files -> Path.GetFileNameWithoutExtension(FileInfo(file).Name), new Image(file) } - let cells = Utils.logTime "Whole analyze" (fun () -> ImageAnalysis.doAnalysis img (Path.GetFileNameWithoutExtension(FileInfo(file).Name)) config) + let nbConcurrentTaskLimit = 4 + let n = Environment.ProcessorCount + + Utils.logTime "Whole analyze" (fun () -> + let results = + images + |> PSeq.map (fun (id, img) -> id, ImageAnalysis.doAnalysis img id (config.Copy())) + |> PSeq.withDegreeOfParallelism (if n > nbConcurrentTaskLimit then nbConcurrentTaskLimit else n) + + for id, cells in results do let total, infected = Utils.countCells cells - fprintf resultFile "File: %s %d %d %.2f\n" file total infected (100. * (float infected) / (float total)) - with - | :? IOException as ex -> Utils.log (sprintf "Unable to open the image '%A': %A" file ex) + fprintf resultFile "File: %s %d %d %.2f\n" id total infected (100. * (float infected) / (float total))) + + //Utils.log (sprintf "== File: %A" file) + //with + //| :? IOException as ex -> Utils.log (sprintf "Unable to open the image '%A': %A" file ex) 0 | Window ->