双流机场到天府机场ADS-B数据导入MATLAB

Evand J2024-02-29 13:37

MATLAB导入数据

导入的数据Excel部分截图：

一些处理

matlab 复制代码

% 导入外部轨迹数据并转成标准形式
clear;clc;
%% 导入&预处理
[NUM,TXT,RAW]=xlsread('2021年10月31日CTU-TFU');
time_cell = RAW(3:end,1); %拉取时间数据（cell）
time_char = char(time_cell);
time = zeros(length(time_cell),1);
for i2 = 1:length(time_cell) %时分秒转换为secend
    if time_cell{1}(5) == '1' && time_cell{length(time_cell)}(5) == '0'
        if time_cell{i2}(5) == '1'
            time(i2) = (time_char(i2,5)*10+time_char(i2,6))*3600+(time_char(i2,8)*10+time_char(i2,9))*60+time_char(i2,11)*10+time_char(i2,12);
        else
            time(i2) = 12*3600+(time_char(i2,5)*10+time_char(i2,6))*3600+(time_char(i2,8)*10+time_char(i2,9))*60+time_char(i2,11)*10+time_char(i2,12);
        end
    else
        time(i2) = (time_char(i2,5)*10+time_char(i2,6))*3600+(time_char(i2,8)*10+time_char(i2,9))*60+time_char(i2,11)*10+time_char(i2,12);
    end
end
time = time-time(1);

longitude = NUM(:,1); %拉取维度数据
latitude = NUM(:,2); %拉取经度数据

direction_cell = RAW(3:end,4); %拉取方向角数据（cell）
direction = zeros(length(direction_cell),1);
for i2 = 1:length(direction_cell)
    if(direction_cell{i2}(5) == '°')
        direction(i2) = str2num(direction_cell{i2}(3:4));
    else
        direction(i2) = str2num(direction_cell{i2}(3:5));
    end
end

kts = NUM(:,3); %速度，节（海里每小时）
kmph = NUM(:,4); %速度，千米每小时
height = NUM(:,5); %高度（米）
%% 转换
% 将时间间隔不一致的数据转换为间隔一致的

结果

经度和维度

高度：

源代码

https://download.csdn.net/download/callmeup/88887561