寻找某一mmsi船号的所有轨迹点


import os
import glob
from pandas.core.frame import DataFrame
import csv
import pandas as pd
import numpy as np
import datetime
mmsi = '219670000'
files = './data/dataset/ais'
outfiles = './data/dataset/3.1'
def findmmsi(mmsi,files,outfiles):
# 按照mmsi号查找报文信息
 filelist = os.listdir(files)
 speed = []
 lon = []
 lat = []
 time=[]
# 遍历文件列表中的每个文件
# 遍历文件列表中的每个文件
 for file in filelist:
    # 构建当前文件的完整路径
    input_file_path = os.path.join(files, file)

    # 打开当前文件
    with open(input_file_path, 'r') as table:
# 逐行读取文件内容,分割成列表并筛选出符合条件的记录
        for line in table:
            data = line.strip().split(',')
            if len(data) > 5 and data[2] == mmsi:
              time.append(datetime.datetime.strptime(data[0], '%d/%m/%Y %H:%M:%S'))
              if data[8].isnumeric():
                      speed.append(float(data[8]))
              else:
                      speed.append(' ')
              lon.append(float(data[4]))
              lat.append(float(data[3]))
    df = pd.DataFrame({'time': time, 'lon': lon, 'lat': lat, 'cog': lon})
    df = df.sort_values(by='time')
    print(df)

    input_file_with_name = file.replace(".csv", f"_{mmsi}.csv")

# 拼接完整的文件路径
    full_path = os.path.join(outfiles, input_file_with_name)

# 保存文件
    df.to_csv(full_path, index=False, encoding='utf-8')
 print("finish")
findmmsi(mmsi,files,outfiles)

数据样例:

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值