Friday, July 23, 2010

matlab 3d plot with a map background


load korea
load geoid
t=0:pi/100:10*pi;
xx=100+50*sin(t);
yy=60+50*cos(t);
pcolor(map);shading interp; colorbar
colormap(demcmap(map))
cmap=colormap;
for N=1:length(xx)
     C(N,1:3)=cmap(max(1,mod(N,64)),:);
end

hold on;
scatter3(xx(:),yy(:),t(:),10,C);
view(3)

2 comments:

Anonymous said...

Hey,
This plot is really good, I would like to make a similar plot but above the South Pole. I could not find a map similar to the one loaded by the command 'load korea', but for Antarctica.
Do you have any idea how to handle it?
Thank you

who said...

Hi,
I would suggest you extract the data from some global topo dataset, i.e., etopo

ShowCalendar