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

Question

0 个投票

POFDE.mat

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