MATLAB Answers

How can I project a 3-D sphere onto a 2-D surface?

113 views (last 30 days)
luc
luc on 11 Dec 2014
Commented: luc on 16 Dec 2014
I got a matrix M with 3 col, n rows. The points created by the matrix are all on a sphere with radius r. In a seperate matrix, let's call it "U", same size. Has either 0's or the exact same row. (indicating parts of the the sphere that have been selected).
So a sphere where nothing has been selected is a 3xn matrix filled with 0's.
I want to project the sphere onto a 2D surface (lets ignore the z, and plot x and y), which gives me a "double" sphere, spanned by points. I used the scatter command to get this result. I shifted half of the circle down in y, so u would get two circles next to eachother.
Now... here's my question:
Can I create a memory efficient 2-D plot, in which the points are translated to "surface parts", in order to better indicate which parts of the sphere has been seen.
The reason for the memory efficient part is due to a constant refresh of the data.
See attached picture for an more visual explanation.
Added question: say I don't cut the circle in half in the middle, but say at 1/3. What would be the most easy way to "fold" the data points outward that would otherwise be obscured by the datapoints that have the same x&y, but different z points.
(see explanation in the left down corner)
Thanks in advance,
  5 Comments
luc
luc on 16 Dec 2014
update.
The problem now lies with some flaw in my program.
The circle should be fully colored in, but it aint.
What is the problem here?
clc;clear all; close all;
N=1; %counter
Point=1; %secondary counter
Q=25; %amount of data points the sphere has
r=20; %radius of the sphere
[X,Y,Z] = sphere(Q); %generating coordinates
%count=Q+1;
count=length(X);
M=zeros(((Q+1)^2),3); %assigning a matrix M to speed up code.
for N=1:(count)
M(((N-1)*count)+1:N*count,1)=r*X(N,:); %%rewriting to one matrix
M(((N-1)*count)+1:N*count,2)=r*Y(N,:);
M(((N-1)*count)+1:N*count,3)=r*Z(N,:);
N=N+1;
end
%%end of code of sphere
count=sqrt(length(M));
count2=length(M);
seen=nan(count2,3); %filling a "seen" matrix in which the points are located that have been "seen".
tic
dt=delaunayTriangulation(M(:,1),M(:,2));
toc
UU=1;
counterrrr=1
for TT=1:count^2
Point=M(TT,:)+0.1*[rand rand rand];
p=rand
if p<=3 %selecting random points -> currently selected EVERY point instead of a random one.
seen=Point;
vertexId = nearestNeighbor(dt, seen(1),seen(2))
triangleId=pointLocation(dt, seen(1),seen(2))
try %this loop needs to run because for some strange reason the triangleId sometimes gives a NaN... While this should not be possible...
tri11 = dt(triangleId, [1:end]);
catch
triangleId=1
counterrrr=counterrrr+1 %the amount of errors RANDOMLY CREATED?!
end
hold on
patch(M(tri11,1), M(tri11,2), 'r', 'LineWidth',1, 'FaceColor','g')
triplot(dt.ConnectivityList(triangleId,:),M(:,1),M(:,2),'c')
tri1(UU,:)=dt(vertexId,:);
UU=UU+1;
end
end

Sign in to comment.

Accepted Answer

Mohammad Abouali
Mohammad Abouali on 11 Dec 2014
Edited: Mohammad Abouali on 11 Dec 2014
This is called map projection. There are multiple methods. check any map projection resource.
The best resource is perhaps this:
MATLAB has a built-in command [x,y] = mfwdtran(lat,lon) There are some other command that might be helpful too, such as projfwd() .
These commands are part of "Mapping Toolbox".
  1 Comment
luc
luc on 16 Dec 2014
I found something called "Delaunay triangulation".
This "sort of" solves my problems, combined with the (very usefull) "patch" code.
But I can't get it to work properly, something is wrong with the indices that it selects to plot.
The current code generates a sphere, ignores the zcoordinates, creates a delaunay triangulation matrix, enters a loop that selects random points from the dataset and then plots the new points. however, when selecting the points something goes wrong.
clc;clear all; close all;
N=1; %counter
Point=1; %secondary counter
Q=25; %amount of data points the sphere has
r=20; %radius of the sphere
[X,Y,Z] = sphere(Q); %generating coordinates
%count=Q+1;
count=length(X);
M=zeros(((Q+1)^2),3); %assigning a matrix M to speed up code.
for N=1:(count)
M(((N-1)*count)+1:N*count,1)=r*X(N,:); %%rewriting to one matrix
M(((N-1)*count)+1:N*count,2)=r*Y(N,:);
M(((N-1)*count)+1:N*count,3)=r*Z(N,:);
N=N+1;
end
%%end of code of sphere
count=sqrt(length(M));
count2=length(M);
seen=nan(count2,3); %filling a "seen" matrix in which the points are located that have been "seen".
dt=delaunayTriangulation(M(:,1),M(:,2));
UU=1;
for TT=1:count^2
Point=M(TT,:)+[0.001 0.001 0.001];
p=rand
if p<=0.5 %selecting random points
seen=Point;
vertexId = nearestNeighbor(dt, seen(1),seen(2));
tri11 = dt(vertexId, [1:end 1]);
hold on
patch(M(tri11,1), M(tri11,2), 'r', 'LineWidth',1, 'FaceColor','g')
tri1(UU,:)=dt(vertexId,:);
UU=UU+1;
end
end
triplot(dt);

Sign in to comment.

More Answers (1)

matt dash
matt dash on 11 Dec 2014
Edited: matt dash on 11 Dec 2014
Ok... based on your code above, if your question is just how you can quickly update the data to reflect new information about which points are on/off, the answer is that you can use NaNs to turn points off (looks like you already know this) and then set the x/y/z data of the trisurf or surf object to quickly update the sphere. Here is an example... on my computer this gives me about 400 frames per second in 2014b, or 60 frames per second in earlier versions... so this should be fast enough to keep up with the max possible speed of your monitor:
[x,y,z]=sphere(50);
data=[x(:) y(:) z(:)];
figure('renderer','opengl')
axes('xlim',[-1 1],'ylim',[-1 1],'zlim',[-1 1])
L = surface(x,y,z);
view(3)
axis square
n=size(data,1);
tic
for i = 1:1000
npts = randi(n);
idx=false(n,1);
idx(1:npts)=true;
idx=idx(randperm(n));
thisx=x;
thisy=y;
thisz=z;
thisx(idx)=nan;
thisy(idx)=nan;
thisz(idx)=nan;
set(L,'xdata',thisx,'ydata',thisy,'zdata',thisz);
drawnow;
end
t=toc
disp(['Frames per second = ',num2str(1000/t)])
  2 Comments
matt dash
matt dash on 11 Dec 2014
I'm not sure about the 2d vs 3d thing. All matlab axes are 3d, they just appear 2d when you view them from above. So a "2D" plot is identical to a 3D plot with all the z values set to 0. I'm not sure if this really provides any performance benefit.
As far as plotting both halves of the sphere separately, i would just define a plane, and use it in a condition that checks which side of the plane each point of the sphere is. Then have one plot show the points on one side, one plot show the points on the other side. Same idea applies if you had a more complex condition (splitting the sphere into 3 regions etc). I'm still not clear if you want to apply some additional transformation, like a map projection, but in any case I can't help you with that. If you happen to have the mapping toolbox it has many transformations built in.

Sign in to comment.

Community Treasure Hunt

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

Start Hunting!