Python 实现常见的坐标系之间的转换

最新推荐文章于 2024-07-23 14:36:35 发布

南山南上南山

最新推荐文章于 2024-07-23 14:36:35 发布

阅读量2.9k

点赞数 2

分类专栏： python 文章标签： python 几何学开发语言

本文链接：https://blog.csdn.net/sinat_35377574/article/details/125223196

版权

本文介绍了使用Python进行常见坐标系转换的方法，通过示例代码详细阐述了转换过程，旨在为GIS开发者提供实用帮助。

摘要由CSDN通过智能技术生成

常见的一些坐标系之间的转换
话不多说,直接上代码

import math

x_pi = 3.14159265358979324 * 3000.0 / 180.0
pi = 3.1415926535897932384626  # π
a = 6378245.0  # 长半轴
ee = 0.00669342162296594323  # 偏心率平方


class CoordinateTransform:

    @classmethod
    def gcj02_to_bd09(cls, lng, lat):
        """
        火星坐标系(GCJ-02)转百度坐标系(BD-09)
        谷歌、高德——>百度
        :param lng:火星坐标经度
        :param lat:火星坐标纬度
        :return:
        """
        z = math.sqrt(lng * lng + lat * lat) + 0.00002 * math.sin(lat * x_pi)
        theta = math.atan2(lat, lng) + 0.000003 * math.cos(lng * x_pi)
        bd_lng = z * math.cos(theta) + 0.0065
        bd_lat = z * math.sin(theta) + 0.006
        return [bd_lng, bd_lat]

    @classmethod
    def bd09_to_gcj02(cls, bd_lng, bd_lat):
        """
        百度坐标系(BD-09)转火星坐标系(GCJ-02)
        百度——>谷歌、高德
        :param bd_lat:百度坐标纬度
        :param bd_lon:百度坐标经度
        :return:转换后的坐标列表形式
        """
        x = bd_lng - 0.0065
        y = bd_lat - 0.006
        z = math.sqrt(x * x + y * y) - 0.00002 * math.sin(y * x_pi)
        theta = math.atan2(y, x) - 0.000003 * math.cos(x * x_pi)
        gcj_lng = z * math.cos(theta)
        gcj_lat = z * math.sin(theta)
        return [gcj_lng, gcj_lat]

    @classmethod
    def wgs84_to_gcj02(cls, lng, lat):
        """
        WGS84转GCJ02(火星坐标系)
        天地图> 高德谷歌
        :param lng:WGS84坐标系的经度
        :param lat:WGS84坐标系的纬度
        :return:
        """
        if cls.out_of_china(lng, lat):  # 判断是否在国内
            return [lng, lat]
        dlat = cls.transform_lat(lng - 105.0, lat - 35.0)
        dlng = cls.transform_lng(lng - 105.0