【Matlab-M_Map】绘制海面温度分布图、等温线
本人萌新,绘制海温图时发现可以查到的代码有很多地方比较晦涩不适合新手,所以在修改的基础上进行了注释,希望对其他同学有所帮助,不足和错误之处还请指正。
文件为nc格式,时间层数一层,海温数据为三维,经纬度和层数,在绘制海表温度分布时取第一层。
clc;
clear;
source1='C:\Users\wanjy\Desktop\BOA_Argo_2020_09.nc';
ncdisp(source1);
time=ncread(source1,'time');
ticount=length(time);
lon=double(ncread(source1,'lon'));
loncount=length(lon);
lat=double(ncread(source1,'lat'));
latcount=length(lat);
temp=double(ncread(source1,'temp'));
boundary = input('输入绘图范围,格式:[起始经度 终止经度 起始纬度 终止纬度]\n');
[~,lon1]=min(abs(lon-boundary(1)));%右边指的是取所有经度与boundary(1)之差的绝对值,取其最小值,即取最近的那个点
[~,lon2]=min(abs(lon-boundary(2)));%[Y,U]=min(A):返回行向量Y和U,Y向量记录A的每列的最小值,U向量记录每列最小值的行号。
[~,lat1]=min(abs(lat-boundary(3)));%当然了,“~”指的是忽略这个值
[~,lat2]=min(ab