python ,经纬度坐标转换 WGS84,BD09、GCJ02(火星坐标系)和UTM坐标系相互转换

"""
坐标转换工具类
"""
import math
from pyproj import Proj, transform, Transformer
import os
from warnings import simplefilter

simplefilter(action='ignore', category=FutureWarning)
simplefilter(action='ignore', category=DeprecationWarning)


# 坐标转换
class TransCoordinates:
    """
    WGS84,BD09、GCJ02(火星坐标系)和UTM坐标系相互转换
    """

    def __init__(self, longitude, latitude):
        """
        初始化
        :param longitude: 经度:float
        :param latitude: 纬度:float
        """
        self.longitude = longitude
        self.latitude = latitude
        self.pi = 3.1415926535897932384626  # π
        self.x_pi = 3.14159265358979324 * 3000.0 / 180.0
        self.a = 6378245.0  # 长半轴
        self.ee = 0.00669342162296594323  # 偏心率平方

    def out_of_china(self, longitude, latitude):
        """
        判断经纬度是否在国内
        :return: 在:True,不在:False
        """
        return not (73.66 < longitude < 135.05 and 3.86 < latitude < 53.55)

    def _transformlat(self, lng, lat):
        """
        纬度转换
        """
        ret = -100.0 + 2.0 * lng + 3.0 * lat + 0.2 * lat * lat + \
              0.1 * lng * lat + 0.2 * math.sqrt(math.fabs(lng))
        ret += (20.0 * math.sin(6.0 * lng * self.pi) + 20.0 *
                math.sin(2.0 * lng * self.pi)) * 2.0 / 3.0
        ret += (20.0 * math.sin(lat * self.pi) + 40.0 *
                math.sin(lat / 3.0 * self.pi)) * 2.0 / 3.0
        ret += (160.0 * math.sin(lat / 12.0
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值