最近在画一个CMAQ输出数据的填色图。数据,经度,纬度是171*171的数组,但数据的经纬度范围有点奇怪,是一个Lambert投影下的正方形,如图所示:
换成mercator投影是这样的
以下是代码
%% 读取数据
clc
clear
close all
load test.mat
%% 预览数据
fig = figure();
h = axesm('MapProjection','lambertstd'); % 创建一个Lambert投影的地图
framem on % 开启地图边框
axis off % 关闭axis
setm(gca,'MapLatLimit',[30,80],'MapLonLimit',[-40,70]) % 设置地图经纬度范围
geoshow(lat,lon,plotdata,'DisplayType','texturemap') % 绘制填色图
geoshow('landareas.shp', 'FaceColor', 'none','LineWidth',1) % 绘制海岸线
setm(gca,'MeridianLabel','on','MLabelLocation',10,'MLabelParallel&#