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