How can I superimpose my contour plot to a map using usamap function?
2 visualizzazioni (ultimi 30 giorni)
Mostra commenti meno recenti
Behrooz Daneshian
il 19 Feb 2023
Commentato: Walter Roberson
il 19 Feb 2023
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
0 Commenti
Risposta accettata
Walter Roberson
il 19 Feb 2023
contourm needs lat then long
2 Commenti
Walter Roberson
il 19 Feb 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.
Più risposte (0)
Vedere anche
Categorie
Scopri di più su Contour Plots in Help Center e File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!