Efficient display of 2D images in 3D space?

22 次查看（过去 30 天）
Prairie geophysics
回答： Mohamed Amine Henchir 2019-4-5
How can I plot a set of 2D ground penetrating radar sections in 3D, as rendered 2D images that have the appearance of 3D slices? [An analogy might be that I have a set of paintings, which I want to stack in a rack. The rack creates a 3D geometry for viewing the paintings, but the canvases themselves are still 2D.]
I have data from a series of parallel radar scans through a highway. Each scan generated a time series that can be arranged as xyzA where in a simple case: x is distance along the line, y is a constant for each line, z is depth, and A is the radar signal's amplitude. A typical scan comprises ~3 million data points.
A single scan can be displayed in 2D using imagesc (using the x,z, and A components). My task is to visualize the data in 3D. In theory, I could interpolate all of the scanned data to make a 3D volume, and slice it, but that will not differentiate between noise and signal, and is likely to create persistent artifacts. Instead, I'd like to just display a series of 2D slices in 3D space, where each slice is an individual scan.
This sounds like a trivial problem (after all, a 3 million point slice is equivalent to a low-resolution photo, and 3D perspective computer games handle this kind of problem in real time). So far, though I've had no success. The prospective solutions that I've found online can't handle the number of points involved, and produce memory errors. I've tried a couple of methods so far:
- 3D scatter plotting where each plotted point represents one data point from the file, and its colour is controlled by the signal amplitude; and - A slicing method that I found via Matlab code sharing. This used meshgrid to generate a distribution of data across a 3D plane, but could not cope with the number of points to be gridded.
I'm wondering whether the solution lies in the way the data are being managed. Would some approach involving projection of a jpeg onto a 3D surface be more efficient than trying to plot each value from a matrix of 3D data points?
An alternative is to go back to the start, and say "OK. If my data were in a commercial radar format, I could just use dedicated radar processing software. So then how can I flip the file from its current proprietary format to be readable by such a package?" I'll do that if this avenue turns out to be a dead end, but for now I'm looking for a Matlab visualization solution.

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

采纳的回答

Mike Garrity
Mike Garrity 2016-3-9
I wouldn't think that going to a scatter plot would be a good approach here.
What about something like this?
img = imread('ngc6543a.jpg');
for z=-100:10:100
g = hgtransform('Matrix',makehgtform('translate',[0 0 z]));
image(g,img)
end
view(3)
  7 个评论
显示 5更早的评论
Walter Roberson
Walter Roberson 2017-10-9
Guojin Feng comments to Prairie Geophysics:
Good solution
Con Ahern
Con Ahern 2018-2-9
I combined this solution with imagesc to scale the image I have to "paint" walls (using multiple rotations and translations) with a certain marble described by a JPEG file. Thanks very much!

请先登录，再进行评论。

更多回答（1 个）

Mohamed Amine Henchir
Hey Mike, Thanks for your solution! I was able to use it to show a stack of images updating every iteration. however I still have another challenge which is to plot detected circles in these images along with the original images on the same stack.
here is my code
[centersBright1, radiiBright1, metric1] = imfindcircles(I,[Rmin Rmid ],'ObjectPolarity',...
'bright','Method','PhaseCode','Sensitivity',slider2Value,'EdgeThreshold',slider1Value);
[centersBright2, radiiBright2, metric2] = imfindcircles(I,[Rmid Rmax],'ObjectPolarity',...
'bright','Method','PhaseCode','Sensitivity',slider2Value,'EdgeThreshold',slider1Value);
[x,y,r] = sort_cirlces(centersBright1,radiiBright1,metric1,centersBright2,radiiBright2,metric2,Rmin,Rmid);
center = [x y];
viscircles(center,r,'EdgeColor','g');
h = hgtransform('Matrix',makehgtform('translate',[0 0 p*10]));
imagesc(h,I)
colormap(gray(65536));
here is a picture of the current output, I would like to have the circles following its corresponding picture, but for some reason they keep being updated at the bottom.
Capture.PNG
any Idea how to fix this? thank you very much!

类别

在 Help Center 和 File Exchange 中查找有关 Vehicle Scenarios 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by