D10, D50, D90 Partical Size Distribution calculation in Image-Pro

Hi Brijeshwar,

Statistical values of Median (D50), D90 and D10 will be available in the next version of Image/Pro. 

In the current version you can use the following macro to get percentile stats of the measurements on the active image (e.g. Diameter):

<pre>    
   Public Sub PrintPercentileStats
        Dim im As McImage=ThisApplication.ActiveImage
        If im Is Nothing Then Exit Sub
        Dim md As MediaCy.Addins.Measurements.McMMData=im.MeasurementsData
        If md.SubFeatures.Count=0 Then Exit Sub' no measurements
        Dim vals As New System.Collections.ArrayList
        'get diameter measurements (change for other measure if necessary)
        Dim measure As eMeasures=eMeasures.RgnMeanDiameter
        For Each sf As McMMSubFeature In md.SubFeatures
            vals.Add(sf.Value(measure))
        Next
        'measure 3 percentiles
        Dim percentiles As Double()={10,50,90}
        For Each p As Double In percentiles
            Debug.Print (String.Format("D{0} of {1} = {2}",p,measure.ToString,GetPercentileValue(vals.ToArray(),p)))
        Next
    End Sub

    Public Function GetPercentileValue(vals As Object, Percentile As Double) As Object
        If vals Is Nothing Then Return 0
        If Not IsArray(vals) Then
            Return vals
        End If
        Dim inLength As Integer = vals.GetUpperBound(0) + 1
        System.Array.Sort(vals) 'sort
        Dim ind As Integer = System.Math.Round(Percentile * inLength / 100, 0, System.MidpointRounding.AwayFromZero) - 1
        ind = System.Math.Max(0, ind)
        ind = System.Math.Min(inLength - 1, ind)
        Return vals(ind)
    End Function
</pre>

the result will look like this:

D10 of RgnMeanDiameter = 14,6595344023042
D50 of RgnMeanDiameter = 17,6075513462399
D90 of RgnMeanDiameter = 21,3247373718319

Yuri

Hi Brijeshwar,

I've created a project that includes these macros. You can just unzip the attached file and load  PercentileStats.ipx project to Image-Pro. Then you can run the PrintPercentileStats macro and the values (D10,D50,D90) will be printed to the Output window:



Let me know if you have questions,

Regards,

Yuri