Making a Map Image from a Matrix

Question

John Ziggs 2021-3-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/771733-making-a-map-image-from-a-matrix

评论： Chad Greene 2021-3-14

Hi all, I am trying to make a map based off data from a text documents containing daily surface air temperature. I attached 2 txt files as a sample. I am trying to create a map based off this.

This is my code:

inputdirTA = 'C:\Users\john\Desktop\TA2004daily'; %directory of folder containing data
S = dir(fullfile(inputdirTA,'*.txt')); 
for k = 1:numel(S) 
    fnm = fullfile(inputdirTA,S(k).name);
    vector = load(fnm);
    mtx = reshape(vector,72,144);
    if k == 1;
        DataTA = [mtx]; 
    else
        DataTA = [DataTA mtx];   
    end
end
figure(1);
rv = [0.4 90 0]; 
worldmap('world');
geoshow(DataTA, rv, 'displaytype', 'texturemap');
C = load('coast');
plotm(C.lat, C.long, 'k');
title('Map of Surface Air Temperature');

I am getting this map:

But I am trying to get this kind of map:

I understand the concept of having to create 3 planes of the matrix, but I am having troubles understanding the syntax.

Any help or advice is greatly appreciated.

Thanks.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Chad Greene 2021-3-13

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/771733-making-a-map-image-from-a-matrix#answer_647008

编辑：Chad Greene 2021-3-13

在 MATLAB Online 中打开

With the Climate Data Toolbox you could do it like this:

T = flipud(reshape(load('TS20040101.txt'),[72 144])); % loads, reshapes, and orients the data correctly
[lat,lon] = cdtgrid(2.5); % creates a 2.5 degree global grid
pcolor(lon,lat,T) % plots the grid
shading interp 
cmocean thermal % sets the colormap
borders('countries','color',rgb('gray'))
xlabel longitude
ylabel latitude 

You could even take it one step further and interpolate to a high-resolution grid, then plot with hillshaded topography:

% Create a 0.1 degree grid: 
[lati,loni] = cdtgrid(0.1); 
% Interpolate the temperature data to the high-res grid:
Ti = interp2(lon,lat,T,loni,lati); 
% Get elevation data on the high res grid: 
Zi = topo_interp(lati,loni); 
% Set ocean elevations to zero: 
Zi(Zi<0) = 0; 
figure
surf(loni,lati,Zi,Ti) 
shading interp 
view(2)
axis tight equal
cmocean thermal 
shadem(-10) % applies hillshade 