基于AIS数据的船舶轨迹聚类方法

%2.DP压缩后的轨迹图 clear;close all;clc; % 清空工作区、关闭所有窗体、清空命令行 Fi = dir('E:\WHUT\matlab数据分析\ais\*.csv'); CF = []; % 设置一个空表，为总数据表 for i = 1:1:length(Fi) FN = strcat(Fi(i).folder,'\',Fi(i).name); % 每条船舶AIS数据的完全路径 name = strsplit(Fi(i).name,'.'); % 根据文件名分离出MMSI na = str2num(name{1}); % 获取船舶MMSI AIS = readtable(FN); % 根据完全路径读取AIS数据表 t1 = AIS.time; % 取出时间长字符串段 t= datetime(t1,'InputFormat','yyyy,MM,dd,HH,mm,ss'); %将时间string转为datetime [x,y] =GaussProjCal(AIS.lon,AIS.lat); % [~,idcf] = dp([x,y],50); s = repmat(na,[length(idcf) 1]); % 根据AIS轨迹个重复MMSI % cf = table(s,t(idcf),AIS.lon(idcf),AIS.lat(idcf)); % 构建一个当前船舶的AIS表 cf.Properties.VariableNames = {'MMSI','Time','Lon','Lat'}; % 创建表的列名 geoplot(cf.Lat, cf.Lon,"LineWidth",0.2,"Marker",".","Color","k"); % 绘制当前船舶的轨迹点连线 % geoscatter(cf.Lat, cf.Lon,'filled'); % 绘制当前船舶的轨迹点连线 hold on % 保留窗体 CF = ([CF;cf]); % 将当前船舶的AIS表格加入总数据表 end title('DP轨迹压缩图'); % 创建图标题 geobasemap topographic; % 加载地貌底图 % writetable(CF,'E:\WHUT\matlab数据分析\Demo_cf\data\AIS_D50.csv'); % 将总数据表保存