How can I superimpose my contour plot to a map using usamap function?

4 visualizaciones (últimos 30 días)
Using the following code I could develope countour maps showing freezing depth across the Alaska states. How can I superimpose these contours to the Alaska map?. When I use usamap function and hold on together, contours will not be displayed on the Alaska map.
clear
close all
clc
load('POFDE.mat');
data = cell2mat(POFDE);
[lat,lon] = meshgrid(unique(data(:,1)),unique(data(:,2)));
for i = 1:(size(data, 2)-2)
figure
usamap("Alaska")
states = readgeotable("usastatehi.shp");
row = states.Name == "Alaska";
alaska = states(row,:);
geoshow(alaska,"FaceColor",[0.3 1.0, 0.675])
textm(alaska.LabelLat,alaska.LabelLon,alaska.Name, ...
"HorizontalAlignment","center")
hold on;
I = scatteredInterpolant(data(:,[1 2]), data(:,i+2));
contourm(lon,lat,min(1,max(0,I(lat,lon))),'ShowText','on')
colorbar
title(['The probability of frost depth exceedance of \Omega_{\delta} = ' num2str(i) ' feet'])
exportgraphics(gca, sprintf('FrostPlot_%d_feet.png', i))
hold off
end

Respuesta aceptada

Walter Roberson
Walter Roberson el 19 de Feb. de 2023
contourm needs lat then long
  2 comentarios
Walter Roberson
Walter Roberson el 19 de Feb. de 2023
You cannot use contourf with an axesm based axes. contourf is only suitable for cartesian axes.
You can possibly use contourm with the 'fill' option.

Iniciar sesión para comentar.

Más respuestas (0)

Categorías

Más información sobre Contour Plots en Help Center y File Exchange.

Etiquetas

Community Treasure Hunt

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

Start Hunting!

Translated by