Topography in a pcolor graph

53 views (last 30 days)
I created a pcolor graph with the magnitude of the velocity of the wind in a specific location. I want to know how to add the topography of the place, I tried two options: with geobasemap('topographic') and with contour(z) where z= height (msl) matrix, but it didn't work.
If I use geobasemap separated from pcolor, I can see the specific area with the topography but without the magnitude of velocity. Should I switch to use geobasemap('topographic') and try to display the magnitude of the wind in that graphic? How is it possible?
The magnitude of the velocity and the coordinates (latitude and longitude) are 3 matrixes that have the same dimensions.
Thanks in advance.

Accepted Answer

Adam Danz
Adam Danz on 28 Oct 2020
Edited: Adam Danz on 28 Oct 2020
The types of graphics objects that can be plotted to geographic axes are limited. For wind velocity, quiver arrows would be nice but neither quiver nor quiverm are accepted by geographic axes.
These options could visualize magnitudes at specific geo coordinates but they wouldn't specify the direction of a vector (ie, wind direction).
You're probably better off using map axes which allow for greater flexibility and support a wider range of graphics objects including quiver.
% Produce base map
load korea5c
% Show the *real* axes
% axis on
% Get axis limits
lonLim = xlim();
latLim = ylim();
% Produce (fake) wind data
[lon,lat] = meshgrid(linspace(lonLim(1),lonLim(2),10), linspace(latLim(1),latLim(2),10));
U = sin(lat);
V = cos(lon);
% Plot wind velocity
hold on
Adam Danz
Adam Danz on 31 Oct 2020
Indeed "topo60c" does not appear in the 2020a documentation. It must have been included from 20b, onward.
I don't know how you were using "egm96geoid" but it's been a function since r2006a.

Sign in to comment.

More Answers (0)

Community Treasure Hunt

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

Start Hunting!