Topography in a pcolor graph
53 views (last 30 days)
Lydia Abigail Garza Coello on 27 Oct 2020
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.
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
% 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