Creating a 3D image from the 3D data points
显示 更早的评论
I have a dataset of 3D "grayscale" images (medical ultrasound dataset) in the form of cartesian coordinates of each voxel and according grayscale intesity values. I want to convert this data to the 3D image matrix (so that it can be used in matlab programming).
I have seen a similar function but for "plotting" named "pcolor" for "2D" images. I want that the output be a "3D" image "matrix" and I don't want only to visualize tha data, so this function is not useful.
-Do you know whether the matlab image processing toolbox (or file-exchange resources,...) has a function that can perform such conversion?
-If the suitable function does not exist in image processing toolbox, how I can program this conversion using matlab? I think that I must do some sort of interpolation, but I don't know exactly how to do it.
your help and comments are greatly appreciated I am looking forward to receiving responses
回答（3 个）
Bjorn Gustavsson
2011-9-4
If you have, or can put, your voxels in a simple plaid grid the matlab function to try first is slice. Then there are a bunch of tools for volume visualization/rendering on the file exchange:
For working with transparency I've had mixed success with matlabs alpha-functions, typically I've tried something like:
alphamap('rampdown');alphamap(0.3 .* alphamap);
8 个评论
leila
2011-9-4
Bjorn Gustavsson
2011-9-4
In addition to the tools Walter points to below you might have to use delaunay3 if you're using a slightly older version of matlab.
Or if your cartesian coordinates are on a plaid grid you could simply do something like this to build your voxel volume:
V = zeros([n1,n2,n3]); % Pre-allocate...
for i1 = 1:n1, % number of points along dim #1
for i2 = 1:n2,
for i3 = 1:n3
V(i1,i2,i3) = YourData(CorrespondingIndices);
end
end
end
Also, "want to work on this image with IP toolbox" is kind of vague - but for some operations you might have to work on 2-D slices, others support 3-D images.
Walter Roberson
2011-9-4
If your cartesian coordinates fall are multiples of a discrete interval apart, but are simply not complete (or not ordered), you can put them in to an gridded array (with no interpolation) by using accumarray() . The index you would use for any one coordinate would be 1 + round((actual coordinate minus minimum coordinate) / coordinate spacing)
leila
2011-9-4
Walter Roberson
2011-9-5
My earlier response about interp3 does not assume the data points are any particular distance apart. You did not respond to that, did not at that time provide the further information about the distribution. Then when Bjorn conditionally proposed a mechanism, making it clear under what conditions it would apply, I offered a refinement of that mechanism, also making it clear under what conditions that would apply.
When you provide all relevant information from the beginning, the risks of getting sidetracked are reduced.
You are not going to be able to use the image processing toolbox directly, as your array will be a 3D array that is not an RGB image. The image processing toolbox can only work directly on arrays that are either 2D (grayscale or binary), or 3D where the 3rd dimension is 3 (R, G, B channels.) (Though there are a couple of things that IPT can do with M x N x 4 arrays that are CMYK)
leila
2011-9-5
Bjorn Gustavsson
2011-9-5
Sure hell it is vague. You could analyse away on images in all sorts of directions with the IP toolbox and two persons setting off in different directions might never use the same special functions again. To correct myself I should have suggested griddata3 for you to use if your version of matlab does not have triscatteredinterp.
leila
2011-9-5
Walter Roberson
2011-9-4
0 个投票
You might be able to use interp3() to convert your data to a grid.
Depending exactly what you want to do with the data, you might find triscatteredinterp to be more useful.
2 个评论
leila
2011-9-5
Bjorn Gustavsson
2011-9-5
If your data points are randomly scattered in 3-D then your best start is to interpolate them to a plaid grid with griddata3. If I recall correctly interp3 expects the input data to be on a plaid grid already.
类别
在 帮助中心 和 File Exchange 中查找有关 Images 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
