%给定地图投影
m_proj('miller','longitudes',[-180 180],'latitudes',[-90 90]);
m_coast('linewidth',1,'color',[0.6 0.6 0.6]);
%利用ginput函数手动选择想要的经纬度坐标点,可以得到在图像上的像素值
m_grid('xtick',[-80 -35],'ytick',[-40 10]); %给定所需区域经纬度坐标
[Amazon_x,Amazon_y] = ginput(4); %利用ginput函数手动选择所需的经纬度坐标
%利用已得到的图上像素值绘制矩形框
clf;
m_proj('miller','longitudes',[-180 180],'latitudes',[-90 90]);
m_coast('linewidth',1,'color',[0.6 0.6 0.6]);
m_grid('xtick',-360:60:360,'ytick',-90:20:90);
%Amazon
%利用line函数,依照点点的顺序绘制矩形框
line([Amazon_x(1),Amazon_x(2)],[Amazon_y(1),Amazon_y(2)],'linewidth',2);
line([Amazon_x(2),Amazon_x(3)],[Amazon_y(2),Amazon_y(3)],'linewidth',2);
line([Amazon_x(3),Amazon_x(4)],[Amazon_y(3),Amazon_y(4)],'linewidth',2);
line([Amazon_x(4),Amazon_x(1)],[Amazon_y(4),Amazon_y(1)],'linewidth',2);
text(Amazon_x(4)+0.1,Amazon_y(4)+0.5,'Amazon','Fontsize',15,'color','b');