MATLAB绘制风速剖面图

Question

MATLAB编码
提取利奇马台风持续期间每一时刻最小风速的经纬度，作为台风位置，将各个位置作为点标出，并连接起来。数据采用era5的nc文件，文件名为“era5_liqima.nc”,绘制利奇马台风风速剖面图，按每一个台风所处位置的前后最近位置与本位置互相之间的连线和正北方向的夹角的平均值，在此位置上作平均夹角的角平分线，以这个角平分线为切面，绘制出以经度为x轴 最大风速为y轴的风速剖面图。在台风持续时间之内，每隔一小时绘制一张图。

Answer 1

% 读取era5_liqima.nc文件中的数据
ncfile = 'era5_liqima.nc';
lat = ncread(ncfile, 'latitude');
lon = ncread(ncfile, 'longitude');
time = ncread(ncfile, 'time');
wind_speed = ncread(ncfile, 'wind_speed');

% 提取利奇马台风持续期间每一时刻最小风速的经纬度
min_wind_speed = min(wind_speed, [], 3);
[min_speed, min_index] = min(min_wind_speed, [], 'all', 'linear');
[min_lat_index, min_lon_index] = ind2sub(size(min_wind_speed), min_index);
typhoon_lat = lat(min_lat_index);
typhoon_lon = lon(min_lon_index);

% 绘制台风路径图
figure;
plot(typhoon_lon, typhoon_lat, 'ro-');
xlabel('Longitude');
ylabel('Latitude');
title('Typhoon Path');

% 绘制风速剖面图
figure;
plot(time, squeeze(wind_speed(min_lat_index, min_lon_index, :)));
xlabel('Time');
ylabel('Wind Speed');
title('Wind Speed Profile');

% 计算角平分线的角度
north_angle = atan2d(diff(typhoon_lat), diff(typhoon_lon));
mean_angle = mean(north_angle);
angle_bisector = mean_angle + 90;

% 绘制风速剖面图的切面
figure;
max_wind_speed = squeeze(max(wind_speed, [], 3));
x = lon;
y = max_wind_speed(min_lat_index, :);
plot(x, y);
xlabel('Longitude');
ylabel('Max Wind Speed');
title('Wind Speed Profile along Angle Bisector');