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

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

Answer 1

[Walter Roberson]

contourm needs lat then long

[Behrooz Daneshian]

sorry that must be contourf

[Walter Roberson]

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.

See also the discussion at https://www.mathworks.com/matlabcentral/answers/1771865-error-in-the-simultaneous-use-of-contourf-and-geolimits#comment_2294530