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         $"{{{nameof Ellipse}: {nameof this.Cx} = %f{this.Cx}, {nameof this.Cy} = %f{this.Cy}, {nameof this.A} = %f{this.A}, {nameof this.B} = %f{this.B}, {nameof this.Alpha} = %f{this.Alpha}}}"
  48
  49     override this.Equals (other : obj) =
  50         match other with
  51         | :? Ellipse as otherEllipse ->
  52             otherEllipse.Cx = this.Cx &&
  53             otherEllipse.Cy = this.Cy &&
  54             otherEllipse.A = this.A &&
  55             otherEllipse.B = this.B &&
  56             otherEllipse.Alpha = this.Alpha
  57         | _ -> false
  58
  59     override this.GetHashCode () = HashCode.Combine (this.Cx, this.Cy, this.A, this.B, this.Alpha)
  60
  61 [<Struct>]
  62 type CellClass = HealthyRBC | InfectedRBC | Peculiar
  63
  64 type Cell =
  65     {
  66         cellClass : CellClass
  67         center : Point
  68         nucleusArea : int
  69         parasiteArea : int
  70         elements : Matrix<byte>
  71     }
  72
  73 [<Struct>]
  74 type Line (a : float32, b : float32) =
  75     member this.A = a
  76     member this.B = b
  77
  78 type MaybeBuilder () =
  79     member this.Bind (x, f) =
  80         match x with
  81         | None -> None
  82         | Some a -> f a
  83
  84     member this.ReturnFrom x = x
  85
  86     member this.TryFinally (body, compensation) =
  87         try
  88             this.ReturnFrom (body ())
  89         finally
  90             compensation ()
  91
  92     member this.Using (disposable : 'a when 'a :> IDisposable, body) =
  93         let body' = fun () -> body disposable
  94         this.TryFinally (body', fun () ->
  95             match disposable with
  96             | null -> ()
  97             | disp -> disp.Dispose ())
  98
  99     member this.Zero () =
 100         None
 101
 102     member this.Return x =
 103         Some x
 104
 105 let maybe = MaybeBuilder ()
 106
 107 type Result<'a> =
 108     | Success of 'a
 109     | Fail of string // Error message.
 110
 111 type ResultBuilder () =
 112     member this.Bind (res, f) =
 113         match res with
 114         | Success value -> f value
 115         | fail -> fail
 116
 117     member this.ReturnFrom (x) = x
 118
 119 let result = ResultBuilder ()
 120
 121 type AnalysisResult =
 122     {
 123         Cells : Cell list
 124         RBCSize_μm : float<μm>
 125         RBCSize_px : float32
 126     }