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');
Thanking you in anticipation. I appreciate your time and efforts for this.
  2 个评论
Walter Roberson
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.
Mahi
Mahi 2018-4-14
Hi Walter, Thanks for letting me know. I have changed that link, it was my mistake. Sorry for that...

请先登录,再进行评论。

采纳的回答

Walter Roberson
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
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.
Mahi
Mahi 2018-4-15
Thanks very much, Walter...!

请先登录,再进行评论。

更多回答(0 个)

产品

Community Treasure Hunt

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

Start Hunting!

Translated by