Volume of retinal lesion
1 view (last 30 days)
Show older comments
Mohammad Mehdi on 23 Mar 2023
Answered: John D'Errico on 23 Mar 2023
I am trying to measure the volume of a retinal lesions in an image obtained using OCT. Each section parameters (X and Z parameters) and the distance between each scans are available. I studied a few paper and it seems that Cavalieri principle is the main method that have been used to measure the volume using the area (by ImageJ). I just wanted to be sure that the proposed method is the accurate method so I can measure the volume of lesions.
My scans of retinal lesion are atttached.
Proposed method to measure the retinal lesions:
It would be great if you could give me your opinion or any suggestions.
John D'Errico on 23 Mar 2023
It seems clear. You have multiple slices. Compute the area of each slice. At that point, you could just sum the area of each slice, multiplied by the thickness of a slice. That would be the rectangle rule, an entirely valid estimator.
The paper suggests using trapezoidal rule, a little better estimator of the volume.
And the implementation is trivial. Compute the area of each slice, as a vector of areas, one element for each slice. Then call trapz, giving it the slice thickness as the X argument.
Find more on Medical Physics in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!Start Hunting!