【Matlab-M_Map】绘制海面温度分布图、等温线

【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(abs(lat-boundary(4)));
varname = input('输入变量：','s');
disp('时间层数为：')
disp(ticount);%如果是海浪之类的有多层时间的数据，还需要在下面加入一行代码，输入需要进行展示的时间层数
for t=1:ticount
    lon_scope = find(lon >= lon(lon1) & lon<=lon(lon2));%这两个经度之间的经度