X-Git-Url: http://git.euphorik.ch/?p=master-thesis.git;a=blobdiff_plain;f=Parasitemia%2FParasitemia%2FProgram.fs;h=09f43126106c852a7f65209546448e8db679f893;hp=ac17c82b0c9d00ab47b4752942338f7fde4ac9e3;hb=044b0ae69df3ac565432545b2fa934589016f9bd;hpb=e76da913cd58078ad2479357b2430ed62a6e0777

diff --git a/Parasitemia/Parasitemia/Program.fs b/Parasitemia/Parasitemia/Program.fs
index ac17c82..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
@@ -57,65 +60,79 @@ let parseArgs (args: string[]) : Arguments =
 let main args =
     match parseArgs args with
     | mode, debug ->
-        let scale = 1.
-        let config = {
-            debug = if debug then DebugOn "." else DebugOff
-
-            scale = scale
+        let config =
+            Config(
+              { scale = 1.
 
-            // RBC size range in px at scale = 1.0.
-            minRBCSize = 20.
-            maxRBCSize = 40.
+                initialAreaOpen = 2000
 
-            doGSigma1 = 1.5
-            doGSigma2 = 20.
-            doGLowFreqPercentageReduction = 0.75
+                minRbcRadius = -0.32f
+                maxRbcRadius = 0.32f
 
-            darkStainLevel = 0.5
+                preFilterSigma = 1.7 // 1.5
 
-            stainSigma = 10.
-            stainLevel = 0.9
-            stainSpreadRequired = 3.0
+                factorNbPick = 1.0
 
-            infectionSigma = 2.2
-            infectionLevel = 0.85
-            infectionPixelsRequired = 1
+                darkStainLevel = 0.22 // Lower -> more sensitive. 0.3. Careful about illumination on the borders.
+                maxDarkStainRatio = 0.1 // 10 %
 
-            percentageOfFgValidCell = 0.4
+                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 = 600. * scale ** 2. |> int
-            maxOffcenter = 0.5 }
+                standardDeviationMaxRatio = 0.5 // 0.55
+                minimumCellArea = 0.5f })
 
         match mode with
         | CmdLine (input, output) ->
-            let config = { config with debug = DebugOn output }
-            Utils.log <- (fun m -> Console.WriteLine m)
+            if debug
+            then
+                config.Debug <- DebugOn output
 
             Directory.CreateDirectory output |> ignore
 
+            use logFile = new StreamWriter(new FileStream(Path.Combine(output, "log.txt"), FileMode.Append, FileAccess.Write))
+            Utils.log <- (fun m -> logFile.WriteLine(m))
+            Utils.log (sprintf "=== New run : %A %A ===" DateTime.Now  (if debug then "[DEBUG]" else "[RELEASE]"))
+
             let files = match input with
                         | File file -> [ file ]
                         | Dir dir -> Directory.EnumerateFiles dir |> List.ofSeq
 
-            use resultFile = new StreamWriter(new FileStream(Path.Combine(output, "results.txt"), FileMode.Create, FileAccess.Write))
+            use resultFile = new StreamWriter(new FileStream(Path.Combine(output, "results.txt"), FileMode.Append, FileAccess.Write))
+
+            //try
+            let images = seq { for file in files -> Path.GetFileNameWithoutExtension(FileInfo(file).Name), new Image<Bgr, byte>(file) }
 
-            for file in files do
-                try
-                    let fileInfo = FileInfo(file)
-                    use img = new Image<Bgr, byte>(file)
-                    let cells = Utils.logTime "Whole analyze" (fun () -> ImageAnalysis.doAnalysis img fileInfo.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
-                | _ 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 ->
             let app = new Application()
             let mainWindow = Views.MainWindow()
 
+            if debug
+            then
+                config.Debug <- DebugOn "."
+
             Utils.log <- (fun m -> log mainWindow m)
 
             //display mainWindow img