利用python进行AIS轨迹绘制

python全部代码

import pylab as plt
import pandas as pd




data = pd.read_excel('/Users/dada/Desktop/2019-04-20.xls')
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[3]]==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()

模块导入

import pylab as plt
import pandas as pd

数据导入

data = pd.read_excel('/Users/dada/Desktop/2019-04-20.xls')
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[3]]==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()