MATLAB 绘制点的地理空间分布，并用点的颜色或大小代表数值

article/2025/6/8 6:00:06

%% 修改日期 2021/12/8clc
clear
close all%% 测试数据，第一列代表纬度，第二列经度，第三列则是点的值
txt = [46.75296619	-69.1022775	0.31290975246.80735808	-69.06131914	0.39115348546.68438136	-68.40644856	-0.52491718547.68233472	-67.83951645	-0.36287248147.73347821	-67.48293127	0.37907463346.56927204	-66.63274775	0.62725311548.82562532	-89.90962194	-0.81023637449.0812938	-88.90238965	0.0249999349.35432864	-86.12624065	0.06561184748.68766235	-85.30347215	-0.90355251763.50718307	-106.5093842	0.04127716354.97197342	-94.08466339	0.03490436463.40093613	-104.4501114	0.07838092760.09253311	-99.99758911	0.48211839952.10863113	-88.39727783	0.11247205761.94615173	-100.0819702	0.02413884760.3741684	-102.3393326	0.35761289654.6315918	-111.2618713	-0.00014206555.32191086	-108.146904	0.07028974654.93228531	-109.2432938	0.75199462557.58321381	-96.2898941	0.0431657456.83683014	-104.2588959	-0.0073507954.33012772	-109.4531403	0.14377455454.71362686	-93.17708588	0.08516170852.77823257	-92.98579407	0.21854175651.54619598	-92.03131104	0.13476404552.16275024	-90.94496155	-0.76109150952.20027542	-89.87015533	0.2554320153.00246048	-89.05403137	0.53392649153.35476303	-89.2333374	0.007532292];%% 绘制点的地理空间分布，并用颜色和大小代表数值
%
figure(1)
worldmap([35 75],[-150 -50])               %加载指定区域的世界地图
geoshow('landareas.shp', 'FaceColor', [0.85 0.85 0.85], 'EdgeColor', [0.8 0.8 0.8]) %加载陆地hco = colorbar;                                   %显示colorbar
set(get(hco,'label'),'string','Trend of Q (mm/year)','FontSize',15);   %加载colorbar标签
load c_MPL_RdBu                               %加载指定colorbar
colormap(mycmap)
caxis([-1 1])                                        %colorbar上下限for i = 1:length(txt)lat = txt(i,1);                                   %点的纬度lon = txt(i,2);                                  %点的经度z =  txt(i,3);                                    %点的数值%% 关键一步，给点赋颜色，geoshow无法像scatter函数一样自动根据点据数据赋颜色% 计算点的颜色在colorbar中的位置，我采用的colorbar是128*3的rgbcolor_loc = round((z+1)/2*128);      color_loc(color_loc>128) = 128;      %说明点的值超过了color的上限（此处是1），赋以最蓝的颜色color_loc(color_loc<1) = 1;             %说明点的值超过了color的上限（此处是-1），赋以最洪的颜色geoshow(lat,lon,'displaytype','point','MarkerFaceColor',mycmap(color_loc,:),'MarkerEdgeColor',[0.5 0.5 0.5],'markersize',10,'Marker','o')
endh = figure(1);
set(h,'units','normalized','position',[0.1 0.1 0.4 0.35]);