船舶AIS数据轨迹可视化python代码.rar

船舶AIS数据轨迹可视化python代码.rar

# -*- coding: utf-8 -*-
“““
Created on Tue Mar 10 22:50:05 2020

@author: TSON
“““
‘‘‘
-------------------------------------------------
  File_Name：    DouglasPeuker drawing
  Author :       L_ZONGYANG
  date：         2020/03/17
-------------------------------------------------
‘‘‘
import pylab as plt
import pandas as pd




data = pd.read_excel（‘C:/Users/TSON/Desktop/2019-01.xlsx‘）
print（data.index）#获取行的索引名称
print（data.columns）#获取列的索引名称

MMSIDATA=data[‘MMSI‘]
print（MMSIDATA）
SHIPSMMSI = list（set（MMSIDATA））
print（SHIPSMMSI）
lens=len（SHIPSMMSI）
print（lens）
for i in range（lens）:
    mmsi=SHIPSMMSI[i]
    ONESHIP=data.iloc[::][data[data.T.index[0]]==mmsi]
    ONESHIP=ONESHIP.sort_values（by=[“Unixtime“]axis=0）
    point_num=len（ONESHIP）
    start=0
    for j in range（point_num）:
        if point_num<2:
            break
        if j+2>point_num:
            break
        time=（ONESHIP.iloc[j+1]）[‘Unixtime‘]-（ONESHIP.iloc[j]）[‘Unixtime‘]
        if time >600:
    
            #print（‘oneshiptrack:time‘）
            start=j 
        else:
            track=ONESHIP[j:j+2]
            x_axis_data = track[‘Lon_d‘]
            y_axis_data = track[‘Lat_d‘]
            plt.plot（x_axis_data y_axis_data ‘-r‘ color=‘green‘ linewidth=1）
            
            #print（‘oneshiptrack:line‘）
            start=j
        

plt.title（‘AIS data‘fontsize=18）
plt.legend_.remove（）
plt.legend（loc=“upper right“）
plt.xlabel（‘longitude‘）
plt.ylabel（‘latitude‘）
plt.show（）   

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1622  2020-06-07 16:31  船舶AIS数据轨迹可视化python代码.py

----------- ---------  ---------- -----  ----

                 1622                    1