船舶AIS数据道格拉斯-普克DP压缩python代码

# -*- coding: utf-8 -*-

“““
  File Name：    DouglasPeuker
  Description :  道格拉斯-普克抽稀算法
  Author :        L_ZONGYANG
  date：          2020/03/17
“““

from __future__ import division
import pylab as plt
from math import sqrt pow
import pandas as pd 
 
THRESHOLD = 0.001  # 阈值


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

def point2LineDistance（point_a point_b point_c）:
    “““
    计算点a到点b c所在直线的距离
    :param point_a:
    :param point_b:
    :param point_c:
    :return:
    “““
    # 首先计算b c 所在直线的斜率和截距
    if point_b[0] == point_c[0]:
        return 9999999
    #斜率
    slope = （point_b[1] - point_c[1]） / （point_b[0] - point_c[0]）
    #截距
    intercept = point_b[1] - slope