How to map elevation in Matlab

elvis asong ZILEFAC

2013 4 5

0 个回答

18 次查看（30 天）

0 个投票

Hi,

I latitude and longitude as well as elevation. I want to show elevation using 'contourf'. Here is my script.

latlim = [47 62]; lonlim = [-124 -88]; figure('Color','white') axesm('miller','MapLatLimit',latlim,'MapLonLimit',lonlim, ... 'Frame','on','Grid','off','MeridianLabel','on', ... 'ParallelLabel','on') axis on geoshow('PrairiesBoundaries.shp', 'FaceColor','white') framem('FEdgeColor','blue','FLineWidth',0.5) hold on [LONG,LAT]=meshgrid(linspace(min(long),max(long)),linspace(min(lat),max(lat))); F = TriScatteredInterp(long,lat,elevation); qz = F(LONG,LAT); [C,h]=contourf(LONG,LAT,qz); colormap jet

My problems: - I need to transform

LONG,LAT]=meshgrid(linspace(min(long),max(long)),linspace(min(lat),max(lat))); F = TriScatteredInterp(long,lat,elevation); qz = F(LONG,LAT); [C,h]=contourf(LONG,LAT,qz);

so that it aligns with the projection

axesm('miller','MapLatLimit',latlim,'MapLonLimit',lonlim, ... 'Frame','on','Grid','off','MeridianLabel','on', ... 'ParallelLabel','on')

- secondly, I would like to improve my contour spacing and labelling.

Thanks so much for your help. Asong.