# -*- 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()
