Plotting data on curvilinear coord projection
4 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
mashtine
am 1 Jul. 2016
Kommentiert: Dushantha Sandaruwan WIJENDRA NAIDHELAGE
am 3 Aug. 2022
Hello,
I have two lat and lon matrices, each 622x810 and they are of a curvilinear projection. I previously just used vectors of lon and lat to plot my data with surfacem (see below) but I am not sure how to do this with lat and lon as a matrix of curvilinear coordinates.
Any suggestions?
ax = worldmap(latlim, lonlim);
S = shaperead('landareas','UseGeoCoords',true);
surfacem(lat, lon, inpdata)
shading interp
geoshow([S.Lat], [S.Lon],'Color','black');
0 Kommentare
Antworten (3)
KSSV
am 1 Jul. 2016
As you have matrices in hand, it shall be very easy to plot what you want. If you want to get vectors of lon, lat from 622x810 matrices of each, try using unique(); this will give you the vectors.
2 Kommentare
Dushantha Sandaruwan WIJENDRA NAIDHELAGE
am 3 Aug. 2022
Is ther anyway to convert the covarian projection (meshgird) into noramal Longitude and latitudes?
José-Luis
am 1 Jul. 2016
It is not clear to me how you want those arrays displayed. If you only want to show the points:
geoshow(S.Lat(:), S.Lon(:),'Color','black');
2 Kommentare
Jonathan Eliashiv
am 18 Okt. 2016
Bearbeitet: Jonathan Eliashiv
am 18 Okt. 2016
Super simple actually.
Go ahead and reshape your coordinates and data into 1-d vectors:
X_curvi = reshape(lon,[],1);
Y_curvi = reshape(lat,[],1);
data_curvi = reshape(inpdata,[],1);
Make a meshgrid that you want to interpolate into:
[lon_grid,lat_grid] = meshgrid(minlon:dx:maxlon,minlat:dy:maxlat)
and then you can use the griddata function:
[~,~,data_rectilinear] = griddata(X_curvi,Y_curvi,data_curvi,...
lon_grid,lat_grid)
and plot away
pcolor(lon_grid,lat_grid,data_rectilinear); shading flat
0 Kommentare
Siehe auch
Kategorien
Mehr zu Map Display finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!