calculate volume under of 3D isosurface plot

5 次查看（过去 30 天）
Abhishek Saini
Abhishek Saini 2019-11-27
编辑： Constantino Carlos Reyes-Aldasoro 2019-11-29
I have a 3D matrix data V(201,101,101). I want to calculate the isosurface volume for a particular isovalue.
How to calculate the volume.
I used following code to generate isosurface.
figure;
p = patch(isosurface(x,y,z,V,-2));
isonormals(x,y,z,V,p);
cdata = smooth3(rand(size(V)),'box',7);
isocolors(x,y,z,cdata,p)
p.Facecolor = 'interp';
p.Edgecolor = 'none';
  4 个评论
显示 2更早的评论
Abhishek Saini
Abhishek Saini 2019-11-28
编辑：Abhishek Saini 2019-11-28
How to get the values of x,y, and z.
I used
p = issurface(x,y,z,V,isovalue).
It returns
p.vertices and
p.faces.
What should I use to calculate the colum and how does it will give volume.
I have attached isosurface image.
Capture.PNG

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Constantino Carlos Reyes-Aldasoro
编辑：Constantino Carlos Reyes-Aldasoro 2019-11-29
What isosurface and isovalue are doing are generating a surface around those points that are above a certain value. All code you showed is good for visualising, but what you really want here is not a visual display but a number of the volume.
If you want the volume that is the same as how many points satisfy your condition of voxel intensity > value, so what you have to do is the following:
volumeOfInterest = (V >-2 );
That will generate a matrix of same dimensions as V with 1 for those points above the isovalue and 0 below, then simply sum those values
volume = sum (volumeOfInterest(:));
And then, if needed you will need to callibrate that value to the physical quantity of your interest (say mm^3) with a conversion between voxels and that quantity.
Hope that helps
  3 个评论
显示 1更早的评论
Walter Roberson
Walter Roberson 2019-11-29
Abhishek Saini: please only use Flag to notify the moderators of problems with posts.
Constantino Carlos Reyes-Aldasoro
Constantino Carlos Reyes-Aldasoro 2019-11-29
True, the code I wrote was only calculating the number of voxels, but if you want to callibrate that number into a certain physical quantity, like m^3 or mm^3 you will need to do that multiplication. Will edit my previous response so that it has the complete answer in case someone new arrives to this post.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Scalar Volume Data 的更多信息

标签

产品


版本

R2018b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by