Map = shaperead('Tailand_n.shp') ;
geoshow(Map,'FaceColor',[0.5,1.0,0.5],'facealpha',0.0, 'EdgeColor','black');%把地图画出来
%geoshow函数:geoshow是用来显示地图数据的函数
geoshow的常用方法:
geoshow(lat,lon)直接输入经纬度,将以该经纬度为中心,生成一小片地图
geoshow(filename)直接输入地图文件名,生成该地图
geoshow(_,Name,Value,)输入相关属性并设定参数,生成所需要的地图。
geoshow的常用属性:
facecolor:定义了地图表面的颜色,需要输入三个参数且均在0至1之间。如果不写此属性,则默认是黄色。
'facealpha',0.0 透明度=0 空心
edgecolor:定义了边界的颜色,输入的三个参数值同上。也可以直接像上面代码中写的,直接写颜色:‘black’
% forest pattern
data=double(tc*100);
data(data==0)=nan; %-LY
[lon,lat]=meshgrid(rb(1,1)+(rb(2,1)-rb(1,1))/size(loss,2)/2:(rb(2,1)-rb(1,1))/size(loss,2):rb(2,1),...
rb(2,2)+(rb(1,2)-rb(2,2))/size(loss,1)/2:(rb(1,2)-rb(2,2))/size(loss,1):rb(1,2));
figure(1);
pcolor(lon,lat,data);
view([0 0 1]);
load liuyyy;%这里load时要添加的是mat格式的文件名称
colormap(liuyyy);
shading flat;
axis equal;
axis([97.2 101.9 14.9 20.6]);
caxis([0 100]);
set(gca,'box','on');
set(gca,'FontSize',13);
%title('Forest area of N. Thailand in 2000');
Map = shaperead('Tailand_n.shp') ;
geoshow(Map,'FaceColor',[0.5,1.0,0.5],'facealpha',0.0, 'EdgeColor','black');%把地图画出来
colorbar;
Result: