module ParasitemiaCore.Granulometry
open System
-open System.IO
open System.Drawing
open Emgu.CV
/// <param name="img"></param>
/// <param name="range">Minimum radius * maximum radius</param>
/// <param name="scale">le 1.0, to speed up the process.</param>
+/// <param name="useOctagon"></param>
let findRadiusByClosing (img : Image<Gray, 'TDepth>) (range : int * int) (scale : float) (useOctagon : bool) : int =
- use scaledImg = if scale = 1. then img else img.Resize(scale, CvEnum.Inter.Area)
+ use scaledImg = if scale = 1. then img else img.Resize (scale, CvEnum.Inter.Area)
let r1, r2 = range
let r1', r2' = roundInt (float r1 * scale), roundInt (float r2 * scale)
// 's' must be odd.
let octagon (s : int) : Matrix<byte> =
if s % 2 = 0 then failwith "s must be odd"
- let m = new Matrix<byte>(Array2D.create s s 1uy)
+ let m = new Matrix<byte> (Array2D.create s s 1uy)
let r = (float s) / (Math.Sqrt 2. + 2.) |> roundInt
for i = 0 to r - 1 do
for j = 0 to r - 1 do
if useOctagon then
(octagon (2 * r - 1)).Mat // It doesn't speed up the process.
else
- CvInvoke.GetStructuringElement(CvEnum.ElementShape.Ellipse, Size(2 * r, 2 * r), Point(-1, -1))
+ CvInvoke.GetStructuringElement (CvEnum.ElementShape.Ellipse, Size (2 * r, 2 * r), Point (-1, -1))
- use closed = scaledImg.MorphologyEx(CvEnum.MorphOp.Close, se, Point(-1, -1), 1, CvEnum.BorderType.Replicate, MCvScalar(0.0))
+ use closed = scaledImg.MorphologyEx (CvEnum.MorphOp.Close, se, Point (-1, -1), 1, CvEnum.BorderType.Replicate, MCvScalar (0.0))
let n = closed.GetSum().Intensity
if r1 > r2 then
failwithf "'radiusRange' invalid : %O" radiusRange
- use imgCopy = img.Copy()
+ use imgCopy = img.Copy ()
let mutable maxDiff = 0.f
let mutable max_r = r1
- Morpho.areaCloseFWithFun imgCopy [ for r in r1 .. r2 -> Math.PI * float r ** 2. |> roundInt, r ] (fun r diff ->
- if r <> r1 && diff > maxDiff then
- maxDiff <- diff
- max_r <- r - 1)
+ Morpho.areaCloseFWithFun imgCopy [ for r in r1 .. r2 -> Math.PI * float r ** 2. |> roundInt, r ] (
+ fun r diff ->
+ if r <> r1 && diff > maxDiff then
+ maxDiff <- diff
+ max_r <- r - 1
+ )
+
max_r