MATLAB Central Newsreader - Plot trajectories on a sphere?

Re: Plot trajectories on a sphere?

Charles Cuell — Wed, 09 Jan 2008 18:58:02 -0500

Re: Plot trajectories on a sphere?

Roger Stafford — Wed, 09 Jan 2008 18:27:01 -0500

roberson@ibd.nrc-cnrc.gc.ca (Walter Roberson) wrote in message <fm32q1
$rak$1@canopus.cc.umanitoba.ca>...
> In article <fm2ofj$evq$1@fred.mathworks.com>,
> Roger Stafford <ellieandrogerxyzzy@mindspring.com.invalid> wrote:
> >If you anticipate that x or y
> >will have numbers outside the range you describe, you would first
perform:
>
> > x = mod(x+200,400)-200;
> > y = mod(y+200,400)-200;
>
> >to do your "toroidal" thing.
>
> If you have a vector that crosses the boundary, then would not
> taking the mod like that would result in a large jump back rather
> than in smooth crossing of the boundary?
-------
Yes, that worried me too, Walter. It is to be hoped that in any given
"trajectory" all the points would be shifted by mod in the same way. I think
we need a lot more clarification from Anil.

Roger Stafford

Re: Plot trajectories on a sphere?

roberson@ibd.nrc-cnrc.gc.ca (Walter Roberson) — Wed, 09 Jan 2008 18:12:17 -0500

In article <fm2ofj$evq$1@fred.mathworks.com>,
Roger Stafford <ellieandrogerxyzzy@mindspring.com.invalid> wrote:
>If you anticipate that x or y
>will have numbers outside the range you describe, you would first perform:

> x = mod(x+200,400)-200;
> y = mod(y+200,400)-200;

>to do your "toroidal" thing.

If you have a vector that crosses the boundary, then would not
taking the mod like that would result in a large jump back rather
than in smooth crossing of the boundary?
--
"No one has the right to destroy another person's belief by
demanding empirical evidence." -- Ann Landers

Re: Plot trajectories on a sphere?

David — Wed, 09 Jan 2008 15:42:02 -0500

"Anil Seth" <a.k.seth.nospamplease@sussex.ac.uk> wrote in
message <fm2li8$ntr$1@fred.mathworks.com>...
> I would be VERY grateful for advice. Basically I want to
> plot an x,y trajectory on the surface of a sphere. The
data
> I have are in the following form:
>
> x: a vector (N) in [-200,200]
> y: a vector (N) in [-200,200]
>
> NB these vectors are toroidal in the sense that -201
maps to
> 199.
>
> What I want:
>
> plot these trajectories on a sphere of fixed radius (ie
non
> deformed).
>
> Simple question: If all this is too much, just tell me
how
> to plot a line on the surface of a sphere!
>
> Thanks!
>
> Thanks

since i can't figure out what kind of data you really have
from what you wrote above, here is something generic.
this plots a sphere with the earth topo on it then draws a
3 line trajectory with vectors at each point projected
onto the surface.

scale=100;

%plot a sphere with radius scale
load('topo.mat','topo','topomap1');
[X,Y,Z]=sphere(24);

X=X*scale;
Y=Y*scale;
Z=Z*scale;

h = surface(X,Y,Z,'FaceColor','texture','CData',topo);
rotate(h,[0 0 1],180)
colormap(topomap1)
hold on

%made up numbers for trajectory and vectors at 4 points
xt=[0;5;6;10];
vxt=[0;0;1;1];
yt=[5;10;15;16];
vyt=[1;1;0;0];
zt=[100;95;90;85];
vzt=[.1;.1;.1;.1];

%convert to spherical coords
[T,P,R]=cart2sph(xt,yt,zt);
%project onto surface of sphere by making radius constant
R(:)=scale;
%convert back to cartesian for plotting
[xt,yt,zt]=sph2cart(T,P,R);

%plot
plot3(xt,yt,zt,'r')
quiver3(xt,yt,zt,vxt,vyt,vzt,1,'y-');

hold off

Re: Plot trajectories on a sphere?

Roger Stafford — Wed, 09 Jan 2008 15:16:03 -0500

"Anil Seth" <a.k.seth.nospamplease@sussex.ac.uk> wrote in message
<fm2li8$ntr$1@fred.mathworks.com>...
> I would be VERY grateful for advice. Basically I want to
> plot an x,y trajectory on the surface of a sphere. The data
> I have are in the following form:
>
> x: a vector (N) in [-200,200]
> y: a vector (N) in [-200,200]
>
> NB these vectors are toroidal in the sense that -201 maps to
> 199.
>
> What I want:
>
> plot these trajectories on a sphere of fixed radius (ie non
> deformed).
>
> Simple question: If all this is too much, just tell me how
> to plot a line on the surface of a sphere!
>
> Thanks!
>
> Thanks
--------
I'm not sure what it is that puzzles you, Anil. If you anticipate that x or y
will have numbers outside the range you describe, you would first perform:

x = mod(x+200,400)-200;
y = mod(y+200,400)-200;

to do your "toroidal" thing.

You can find the corresponding z vector from:

z = z0 + sqrt(r^2-(x-x0).^2-(y-y0).^2);

for the upper surface (subtract the square root for the lower surface) of the
sphere with radius r and center at (x0,y0,z0). (Any z's that are imaginary
would be points that don't project onto the sphere.)

If your successive points in the (x,y,z) vectors are far apart and you wish to
connect them with a great circle arc on the sphere rather than a straight line,
that can be done in terms of cross products of three-dimensional vectors.
Let me know if this latter is your actual problem and I can give the necessary
formula.

Roger Stafford

Plot trajectories on a sphere?

Anil Seth — Wed, 09 Jan 2008 14:26:16 -0500

I would be VERY grateful for advice. Basically I want to
plot an x,y trajectory on the surface of a sphere. The data
I have are in the following form:

x: a vector (N) in [-200,200]
y: a vector (N) in [-200,200]

NB these vectors are toroidal in the sense that -201 maps to
199.

What I want:

plot these trajectories on a sphere of fixed radius (ie non
deformed).

Simple question: If all this is too much, just tell me how
to plot a line on the surface of a sphere!

Thanks!

Thanks