Matlab绘制台风路径--数据来源：中国气象局热带气旋资料中心

% 读取台风数据
fid = fopen('CH2009BST.txt', 'r');
data = textscan(fid, '%s', 'Delimiter', '\n');
fclose(fid);
data = data{1};

% 提取台风Morakot数据
typhoon_data = [];
is_dora = false;
for i = 1:length(data)
    line = data{i};
    % 检查是否是Morakot台风的起始行
    if length(line) >= 66 && strcmp(line(1:5), '66666') && contains(line, 'Morakot')
        is_dora = true;
        continue;
    end
    % 检查是否到达下一个台风的起始行
    if length(line) >= 5 && strcmp(line(1:5), '66666') && is_dora
        break;
    end
    % 提取Morakot台风数据
    if is_dora && length(line) >= 33
        year = str2double(line(1:4));
        month = str2double(line(5:6));
        day = str2double(line(7:8));
        hour = str2double(line(9:10));
        intensity = str2double(line(12));
        lat = str2double(line(14:16)) / 10;