Parasitemia/ParasitemiaCore/Types.fs

   1 module ParasitemiaCore.Types
   2
   3 open System
   4 open System.Drawing
   5 open System.Collections.Generic
   6
   7 open Emgu.CV
   8 open Emgu.CV.Structure
   9
  10 open Const
  11 open UnitsOfMeasure
  12
  13 type Points = HashSet<Point>
  14
  15 type Ellipse (cx : float32, cy : float32, a : float32, b : float32, alpha : float32) =
  16     member this.Cx = cx
  17     member this.Cy = cy
  18     member this.A = a
  19     member this.B = b
  20     member this.Alpha = alpha
  21     member this.Area = a * b * PI
  22
  23     // Does the ellipse contain the point (x, y)?
  24     member this.Contains x y =
  25         ((x - cx) * cos alpha + (y - cy) * sin alpha) ** 2.f / a ** 2.f + ((x - cx) * sin alpha - (y - cy) * cos alpha) ** 2.f / b ** 2.f <= 1.f
  26
  27     member this.CutAVericalLine (y : float32) : bool =
  28         a ** 2.f + b ** 2.f - 2.f * y ** 2.f + 4.f * y * cx - 2.f * cx ** 2.f + a ** 2.f * cos (2.f * alpha) - b ** 2.f * cos (2.f * alpha) > 0.f
  29
  30     member this.CutAnHorizontalLine (x : float32) : bool =
  31         a ** 2.f + b ** 2.f - 2.f * x ** 2.f + 4.f * x * cy - 2.f * cy ** 2.f - a ** 2.f * cos (2.f * alpha) + b ** 2.f * cos (2.f * alpha) > 0.f
  32
  33     member this.IsOutside (width : float32) (height : float32) =
  34         this.Cx < 0.f || this.Cx >= width ||
  35         this.Cy < 0.f || this.Cy >= height ||
  36         this.CutAVericalLine 0.f || this.CutAVericalLine width ||
  37         this.CutAnHorizontalLine 0.f || this.CutAnHorizontalLine height
  38
  39     member this.Scale (factor : float32) : Ellipse =
  40         Ellipse (this.Cx, this.Cy, this.A * factor, this.B * factor, alpha)
  41
  42     // Approximation of Ramanujan.
  43     member this.Perimeter =
  44         PI * (3.f * (this.A + this.B) - sqrt ((3.f * this.A + this.B) * (this.A + 3.f * this.B)))
  45
  46     override this.ToString () =
  47         sprintf "{Ellipse: cx = %f, cy = %f, a = %f, b = %f, alpha = %f}" this.Cx this.Cy this.A this.B this.Alpha
  48
  49 [<Struct>]
  50 type CellClass = HealthyRBC | InfectedRBC | Peculiar
  51
  52 type Cell =
  53     {
  54         cellClass : CellClass
  55         center : Point
  56         nucleusArea : int
  57         parasiteArea : int
  58         elements : Matrix<byte>
  59     }
  60
  61 [<Struct>]
  62 type Line (a : float32, b : float32) =
  63     member this.A = a
  64     member this.B = b
  65
  66 type MaybeBuilder () =
  67     member this.Bind (x, f) =
  68         match x with
  69         | None -> None
  70         | Some a -> f a
  71
  72     member this.ReturnFrom x = x
  73
  74     member this.TryFinally (body, compensation) =
  75         try
  76             this.ReturnFrom (body ())
  77         finally
  78             compensation ()
  79
  80     member this.Using (disposable : 'a when 'a :> IDisposable, body) =
  81         let body' = fun () -> body disposable
  82         this.TryFinally (body', fun () ->
  83             match disposable with
  84             | null -> ()
  85             | disp -> disp.Dispose ())
  86
  87     member this.Zero () =
  88         None
  89
  90     member this.Return x =
  91         Some x
  92
  93 let maybe = MaybeBuilder ()
  94
  95 type Result<'a> =
  96     | Success of 'a
  97     | Fail of string // Error message.
  98
  99 type ResultBuilder () =
 100     member this.Bind (res, f) =
 101         match res with
 102         | Success value -> f value
 103         | fail -> fail
 104
 105     member this.ReturnFrom (x) = x
 106
 107 let result = ResultBuilder ()
 108
 109 type AnalysisResult =
 110     {
 111         Cells : Cell list
 112         RBCSize_μm : float<μm>
 113         RBCSize_px : float32
 114     }