Use float32 to reduce memory footprint.

[master-thesis.git] / Parasitemia / Parasitemia / Types.fs

diff --git a/Parasitemia/Parasitemia/Types.fs b/Parasitemia/Parasitemia/Types.fs
index 308fc6b..6240519 100644 (file)
--- a/Parasitemia/Parasitemia/Types.fs
+++ b/Parasitemia/Parasitemia/Types.fs
@@ -1,8 +1,63 @@
 module Types
 
 open System
+open System.Drawing
 
-type Ellipse = { cx: float; cy: float; a: float; b: float; alpha: float }
-    
+open Emgu.CV
+open Emgu.CV.Structure
+
+open Const
+
+type Ellipse (cx: float32, cy: float32, a: float32, b: float32, alpha: float32) =
+    member this.Cx = cx
+    member this.Cy = cy
+    member this.A = a
+    member this.B = b
+    member this.Alpha = alpha
+    member this.Area = a * b * PI
+
+    // Does the ellipse contain the point (x, y)?.
+    member this.Contains x y =
+        ((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
+
+    member this.CutAVericalLine (y: float32) : bool =
+        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
+
+    member this.CutAnHorizontalLine (x: float32) : bool =
+        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
+
+    member this.isOutside (width: float32) (height: float32) =
+        this.Cx < 0.f || this.Cx >= width ||
+        this.Cy < 0.f || this.Cy >= height ||
+        this.CutAVericalLine 0.f || this.CutAVericalLine width ||
+        this.CutAnHorizontalLine 0.f || this.CutAnHorizontalLine height
+
+    member this.Scale (factor: float32) =
+        Ellipse(this.Cx, this.Cy, this.A * factor, this.B * factor, alpha)
+
+    // Approximation of Ramanujan.
+    member this.Perimeter =
+        PI * (3.f * (this.A + this.B) - sqrt ((3.f * this.A + this.B) * (this.A + 3.f * this.B)))
+
+    override this.ToString () =
+        sprintf "(cx: %A, cy: %A, a: %A, b: %A, alpha: %A)" this.Cx this.Cy this.A this.B this.Alpha
+
+
+type CellClass = HealthyRBC | InfectedRBC | Peculiar
+
+type Cell = {
+    cellClass: CellClass
+    center: Point
+    elements: Matrix<byte> }
+
+[<Struct>]
+type Line (a: float32, b: float32) =
+    member this.A = a
+    member this.B = b
+    member this.Valid = not (Single.IsInfinity this.A)
+
+[<Struct>]
+type PointD (x: float32, y: float32) =
+    member this.X = x
+    member this.Y = y
 
-//type PointImg = { x: int; y: int }
\ No newline at end of file