Volume of retinal lesion
4 views (last 30 days)
Show older comments
Greetings everyone,
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:
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1334520/image.jpeg)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1334525/image.jpeg)
It would be great if you could give me your opinion or any suggestions.
Thank you
0 Comments
Answers (1)
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.
help trapz
0 Comments
See Also
Categories
Find more on Read, Write, and Modify Image 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!