尝试了一下在Matlab下,在地图上绘制数据,但是效果没有python的basemap这个好。
这里以卫星重力数据画图为例,记录一下使用过程。
数据来源 :采用Eigen6C4全球重力位模型,用开源软件GrafLab,计算g r g_rgr, g λ g_\lambdagλ, g ϕ g_\phigϕ, 并画出其中的g r g_rgr。该软件计算时自动加上了向心力。
将计算的结果输出到total.txt文件里,画图流程如下
clear all
close all
data=load('total.txt');
% data=load('total_magnitude.txt');
latitude=data(:,1);
longitude=data(:,2);
gx_topo_normal=data(:,3);
gy_topo_normal=data(:,4);
gz_topo_normal=data(:,5);
LATI=reshape(latitude,181,360);
LONGI=reshape(longitude,181,360);
GZ=reshape(gz_topo_normal,181,360);
GZ=-GZ*1e-5;%把方向改为朝下
%画图
clear all
close all
data=load('total.txt');
% data=load('total_magnitude.txt');
latitude&#