印度洋涡旋统计

clear;clc;clf;close all;
%% 读取涡旋半径、时间范围、经纬度信息，以及涡旋点的经纬度；半径路径和涡旋极性；
file2=‘D:\matlab_work\accept_work\涡旋统计的平面分布\eddy_trajectory_2.0exp_19930101_20200307.nc’;
time=double(ncread(file2,‘time’));
lat=double(ncread(file2,‘latitude’));
lon=double(ncread(file2,‘longitude’));
cyclonic_type=double(ncread(file2,‘cyclonic_type’));
track=double(ncread(file2,‘track’));
days=time+datenum(1950,01,01);
datetimed=datevec(days);
clear time file2;
% 定义目标区域，用于寻找自己想找的涡旋区域
aim_eddy_area=[30 105;0 30 ];
% 年份开始日期和结束日期；
yea