Improve the thinning process performance.

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

diff --git a/Parasitemia/Parasitemia/Types.fs b/Parasitemia/Parasitemia/Types.fs
index f96573d..d7bacbd 100644 (file)
--- a/Parasitemia/Parasitemia/Types.fs
+++ b/Parasitemia/Parasitemia/Types.fs
@@ -1,6 +1,10 @@
 module Types
 
 open System
+open System.Drawing
+
+open Emgu.CV
+open Emgu.CV.Structure
 
 type Ellipse (cx: float, cy: float, a: float, b: float, alpha: float) =
     member this.Cx = cx
@@ -11,17 +15,36 @@ type Ellipse (cx: float, cy: float, a: float, b: float, alpha: float) =
     member this.Area = a * b * Math.PI
 
     // Does the ellipse contain the point (x, y)?.
-    member this.Contains x y = 
+    member this.Contains x y =
         ((x - cx) * cos alpha + (y - cy) * sin alpha) ** 2.0 / a ** 2.0 + ((x - cx) * sin alpha - (y - cy) * cos alpha) ** 2.0 / b ** 2.0 <= 1.0
 
-    // A line is defined as : y = mx + l
-    member this.CutALine (m: float) (l: float) : bool =
-        -2.0 * l ** 2.0 + a ** 2.0 + m ** 2.0 * a ** 2.0 + b ** 2.0 + m ** 2.0 * b ** 2.0 - 
-        4.0 * m * l * cx - 2.0 * m ** 2.0 * cx ** 2.0 + 4.0 * l * cy + 4.0 * m * cx * cy - 
-        2.0 * cy ** 2.0 + (-1.0 + m ** 2.0) * (a ** 2.0 - b ** 2.0) * cos (2.0 * alpha) - 2.0 * m * (a ** 2.0 - b ** 2.0) * sin (2.0 * alpha) > 0.0
-
-    member this.CutAVericalLine (y: float) : bool = 
+    member this.CutAVericalLine (y: float) : bool =
         a ** 2.0 + b ** 2.0 - 2.0 * y ** 2.0 + 4.0 * y * cx - 2.0 * cx ** 2.0 + a ** 2.0 * cos (2.0 * alpha) - b ** 2.0 * cos (2.0 * alpha) > 0.0
 
-    member this.CutAnHorizontalLine (x: float) : bool = 
-        a ** 2.0 + b ** 2.0 - 2.0 * x ** 2.0 + 4.0 * x * cy - 2.0 * cy ** 2.0 - a ** 2.0 * cos (2.0 * alpha) + b ** 2.0 * cos (2.0 * alpha) > 0.0
\ No newline at end of file
+    member this.CutAnHorizontalLine (x: float) : bool =
+        a ** 2.0 + b ** 2.0 - 2.0 * x ** 2.0 + 4.0 * x * cy - 2.0 * cy ** 2.0 - a ** 2.0 * cos (2.0 * alpha) + b ** 2.0 * cos (2.0 * alpha) > 0.0
+
+    member this.isOutside (width: float) (height: float) =
+        this.Cx < 0.0 || this.Cx >= width ||
+        this.Cy < 0.0 || this.Cy >= height ||
+        this.CutAVericalLine 0.0 || this.CutAVericalLine width ||
+        this.CutAnHorizontalLine 0.0 || this.CutAnHorizontalLine height
+
+
+type CellClass = HealthyRBC | InfectedRBC | Peculiar
+
+type Cell = {
+    cellClass: CellClass
+    center: Point
+    elements: Matrix<byte> }
+
+[<Struct>]
+type Line (a: float, b: float) =
+    member this.A = a
+    member this.B = b
+
+[<Struct>]
+type PointD (x: float, y: float) =
+    member this.X = x
+    member this.Y = y
+