MATLAB文件读取与地图可视化：地理信息展示

gmt软件学习专业性强，对于一般人很难迅速上手。而m_map工具箱提供了较好的制图风格。

配置过程可参考： m_map配置

m_map官网截图

下面将分享一下如何在matlab中读取带经纬度的文件，并绘制在地图中。

本例以台风轨迹数据为例。

for j =1 :6 %6个台风轨迹文件

A = textread([add,'s_',num2str(j),'.txt'])

subplot(2,3,j) %循环绘制子图

plot(A(:,1),A(:,2))

m_proj('Miller Cylindrical','longtitude',[60 180],'latitude',[-20 60]); %设置投影方式与显示范围

m_coast; %绘制底图

% m_grid('box','plain','tickdir','out');

m_grid('box','fancy','tickdir','in');

% A=textread('yangtze.txt'); 读取长江边界

for i=1:length(A)

    [X,Y]=m_ll2xy(A(i,1),A(i,2));%convert cordinates  %经纬度转xy

    line(X,Y,'marker','.','markersize',10,'color','b');

end

[X,Y]=m_ll2xy(106,30);

text(X,Y,'Yangtze Basin','vertical','top'); %添加文字

end

得到以下结果：

6个台风轨迹

欢迎与我联系交流！1565903211