Questions on Isosurface() command
3 次查看(过去 30 天)
显示 更早的评论
Hello Everyone,
I have some questions regarding isosurface command. Image daset has 2D slices of CT images and can be represented by pixelsinXDir * pixelsinYDir * numberofSlices
1) Is there any specific research paper on which this command has been derived. I know that ultimately it uses interp3(). However, I want to learn theory behind it as well.
2) When we calculate fv=isosurface(imagename,isovalue); it generated structure with Faces and Vertices. Vertices matrix can easily be under stood with its X, Y and Z co-ordinates. What does matrix of face give?
3) Will the final volume be the Voxel data? if so then how the volume is to be calculated? (Using sum(imageset(:)) this?)
4) In order to calculate Centroid of the 3D volume how to proceed? I have written a code as well. Is it correct?
CC=bwconncomp(imgs, 26); % imgs is the image dataset 3D
Center=regionprops(CC,'Centroid');
Please find the web URL related to the above question here: https://www.mathworks.com/help/matlab/visualize/techniques-for-visualizing-scalar-volume-data.html
Thanking you in anticipation. I appreciate your time and efforts for this.
2 个评论
Walter Roberson
2018-4-14
I think you have given the wrong URL for the dataset. You just gave a link to my FAQ answer from several years ago about sites that supported file uploads back then. These days you would probably use dropbox or google drive for larger files, or would upload the data here if it is not so big.
采纳的回答
Walter Roberson
2018-4-14
2) The matrix of faces tells which vertices are connected to which.
3) What "final volume" ? isosurface() is for drawing, not for computing volume.
4) No, that code with bwconncomp calculates the centroid of each of the connected components independently, not the centroid of the 3D volume. For the centroid of the 3D volume you would regionprops Centroid of a binary version of the dataset, with the binary being whether the voxel was occupied or not.
7 个评论
Walter Roberson
2018-4-15
Yes, the unit would be in voxels.
delta x and delta y are typically the same for DICOM images, but delta z is not necessarily the same. I have not used much CT so I do not know the convention for those, but for MRI images it was not uncommon that delta z was twice as large as delta x or delta y.
You would get real-world volume by multiplying voxels by delta x * delta y * delta z.
更多回答(0 个)
另请参阅
产品
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!