How can I put Alaska map as a background behind my contour plots?

Question

Behrooz Daneshian 2023-2-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1914680-how-can-i-put-alaska-map-as-a-background-behind-my-contour-plots

回答： Star Strider 2023-3-15

POFDE.mat

Hi all,

Using the following code I could develope countour maps showing freezing depth across the Alaska states. How can I plot these contour maps in a way that the background would be the Alaska map? The x and y axis are longitude and latitude respectively.

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
    I = scatteredInterpolant(data(:,[1 2]), data(:,i+2));
    contourf(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))
end

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Ashfaq Ahmed 2023-3-15

You will find the map of Alaska and easily plot this using this toolbox -

Arctic Mapping Toolbox

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

the cyclist 2023-2-18

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1914680-how-can-i-put-alaska-map-as-a-background-behind-my-contour-plots#answer_1174420

If you have the Mapping Toolbox, you can use the usamap function to make an Alaska map, and superimpose your other data on it.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Behrooz Daneshian 2023-2-18

在 MATLAB Online 中打开

Thank you for your reply. However, when I use usamap function and hold on together, contours will not be displayed!

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