欧美bbbwbbbw肥妇,免费乱码人妻系列日韩,一级黄片

Python實(shí)現(xiàn)常見(jiàn)的4種坐標(biāo)互相轉(zhuǎn)換

 更新時(shí)間:2023年02月09日 08:25:06   作者:費(fèi)弗里  
主流被使用的地理坐標(biāo)系并不統(tǒng)一,常用的有WGS84、GCJ02(火星坐標(biāo)系)、BD09(百度坐標(biāo)系)以及百度地圖中保存矢量信息的web墨卡托,本文利用Python編寫(xiě)相關(guān)類(lèi)以實(shí)現(xiàn)4種坐標(biāo)系統(tǒng)之間的互相轉(zhuǎn)換,希望對(duì)大家有所幫助

一、簡(jiǎn)介

主流被使用的地理坐標(biāo)系并不統(tǒng)一,常用的有WGS84、GCJ02(火星坐標(biāo)系)、BD09(百度坐標(biāo)系)以及百度地圖中保存矢量信息的web墨卡托,本文利用Python編寫(xiě)相關(guān)類(lèi)以實(shí)現(xiàn)4種坐標(biāo)系統(tǒng)之間的互相轉(zhuǎn)換。

二、代碼及說(shuō)明

import math

class LngLatTransfer():

    def __init__(self):
        self.x_pi = 3.14159265358979324 * 3000.0 / 180.0
        self.pi = math.pi  # π
        self.a = 6378245.0  # 長(zhǎng)半軸
        self.es = 0.00669342162296594323  # 偏心率平方
        pass

    def GCJ02_to_BD09(self, gcj_lng, gcj_lat):
        """
        實(shí)現(xiàn)GCJ02向BD09坐標(biāo)系的轉(zhuǎn)換
        :param lng: GCJ02坐標(biāo)系下的經(jīng)度
        :param lat: GCJ02坐標(biāo)系下的緯度
        :return: 轉(zhuǎn)換后的BD09下經(jīng)緯度
        """
        z = math.sqrt(gcj_lng * gcj_lng + gcj_lat * gcj_lat) + 0.00002 * math.sin(gcj_lat * self.x_pi)
        theta = math.atan2(gcj_lat, gcj_lng) + 0.000003 * math.cos(gcj_lng * self.x_pi)
        bd_lng = z * math.cos(theta) + 0.0065
        bd_lat = z * math.sin(theta) + 0.006
        return bd_lng, bd_lat


    def BD09_to_GCJ02(self, bd_lng, bd_lat):
        '''
        實(shí)現(xiàn)BD09坐標(biāo)系向GCJ02坐標(biāo)系的轉(zhuǎn)換
        :param bd_lng: BD09坐標(biāo)系下的經(jīng)度
        :param bd_lat: BD09坐標(biāo)系下的緯度
        :return: 轉(zhuǎn)換后的GCJ02下經(jīng)緯度
        '''
        x = bd_lng - 0.0065
        y = bd_lat - 0.006
        z = math.sqrt(x * x + y * y) - 0.00002 * math.sin(y * self.x_pi)
        theta = math.atan2(y, x) - 0.000003 * math.cos(x * self.x_pi)
        gcj_lng = z * math.cos(theta)
        gcj_lat = z * math.sin(theta)
        return gcj_lng, gcj_lat


    def WGS84_to_GCJ02(self, lng, lat):
        '''
        實(shí)現(xiàn)WGS84坐標(biāo)系向GCJ02坐標(biāo)系的轉(zhuǎn)換
        :param lng: WGS84坐標(biāo)系下的經(jīng)度
        :param lat: WGS84坐標(biāo)系下的緯度
        :return: 轉(zhuǎn)換后的GCJ02下經(jīng)緯度
        '''
        dlat = self._transformlat(lng - 105.0, lat - 35.0)
        dlng = self._transformlng(lng - 105.0, lat - 35.0)
        radlat = lat / 180.0 * self.pi
        magic = math.sin(radlat)
        magic = 1 - self.es * magic * magic
        sqrtmagic = math.sqrt(magic)
        dlat = (dlat * 180.0) / ((self.a * (1 - self.es)) / (magic * sqrtmagic) * self.pi)
        dlng = (dlng * 180.0) / (self.a / sqrtmagic * math.cos(radlat) * self.pi)
        gcj_lng = lat + dlat
        gcj_lat = lng + dlng
        return gcj_lng, gcj_lat


    def GCJ02_to_WGS84(self, gcj_lng, gcj_lat):
        '''
        實(shí)現(xiàn)GCJ02坐標(biāo)系向WGS84坐標(biāo)系的轉(zhuǎn)換
        :param gcj_lng: GCJ02坐標(biāo)系下的經(jīng)度
        :param gcj_lat: GCJ02坐標(biāo)系下的緯度
        :return: 轉(zhuǎn)換后的WGS84下經(jīng)緯度
        '''
        dlat = self._transformlat(gcj_lng - 105.0, gcj_lat - 35.0)
        dlng = self._transformlng(gcj_lng - 105.0, gcj_lat - 35.0)
        radlat = gcj_lat / 180.0 * self.pi
        magic = math.sin(radlat)
        magic = 1 - self.es * magic * magic
        sqrtmagic = math.sqrt(magic)
        dlat = (dlat * 180.0) / ((self.a * (1 - self.es)) / (magic * sqrtmagic) * self.pi)
        dlng = (dlng * 180.0) / (self.a / sqrtmagic * math.cos(radlat) * self.pi)
        mglat = gcj_lat + dlat
        mglng = gcj_lng + dlng
        lng = gcj_lng * 2 - mglng
        lat = gcj_lat * 2 - mglat
        return lng, lat


    def BD09_to_WGS84(self, bd_lng, bd_lat):
        '''
        實(shí)現(xiàn)BD09坐標(biāo)系向WGS84坐標(biāo)系的轉(zhuǎn)換
        :param bd_lng: BD09坐標(biāo)系下的經(jīng)度
        :param bd_lat: BD09坐標(biāo)系下的緯度
        :return: 轉(zhuǎn)換后的WGS84下經(jīng)緯度
        '''
        lng, lat = self.BD09_to_GCJ02(bd_lng, bd_lat)
        return self.GCJ02_to_WGS84(lng, lat)


    def WGS84_to_BD09(self, lng, lat):
        '''
        實(shí)現(xiàn)WGS84坐標(biāo)系向BD09坐標(biāo)系的轉(zhuǎn)換
        :param lng: WGS84坐標(biāo)系下的經(jīng)度
        :param lat: WGS84坐標(biāo)系下的緯度
        :return: 轉(zhuǎn)換后的BD09下經(jīng)緯度
        '''
        lng, lat = self.WGS84_to_GCJ02(lng, lat)
        return self.GCJ02_to_BD09(lng, lat)


    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 * self.pi) + 320 *
                math.sin(lat * self.pi / 30.0)) * 2.0 / 3.0
        return ret


    def _transformlng(self, lng, lat):
        ret = 300.0 + lng + 2.0 * lat + 0.1 * lng * lng + \
              0.1 * lng * lat + 0.1 * 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(lng * self.pi) + 40.0 *
                math.sin(lng / 3.0 * self.pi)) * 2.0 / 3.0
        ret += (150.0 * math.sin(lng / 12.0 * self.pi) + 300.0 *
                math.sin(lng / 30.0 * self.pi)) * 2.0 / 3.0
        return ret

    def WGS84_to_WebMercator(self, lng, lat):
        '''
        實(shí)現(xiàn)WGS84向web墨卡托的轉(zhuǎn)換
        :param lng: WGS84經(jīng)度
        :param lat: WGS84緯度
        :return: 轉(zhuǎn)換后的web墨卡托坐標(biāo)
        '''
        x = lng * 20037508.342789 / 180
        y = math.log(math.tan((90 + lat) * self.pi / 360)) / (self.pi / 180)
        y = y * 20037508.34789 / 180
        return x, y

    def WebMercator_to_WGS84(self, x, y):
        '''
        實(shí)現(xiàn)web墨卡托向WGS84的轉(zhuǎn)換
        :param x: web墨卡托x坐標(biāo)
        :param y: web墨卡托y坐標(biāo)
        :return: 轉(zhuǎn)換后的WGS84經(jīng)緯度
        '''
        lng = x / 20037508.34 * 180
        lat = y / 20037508.34 * 180
        lat = 180 / self.pi * (2 * math.atan(math.exp(lat * self.pi / 180)) - self.pi / 2)
        return lng, lat

整個(gè)模塊的使用方式可用下面的導(dǎo)圖概括,其中每個(gè)函數(shù)都只需要傳入經(jīng)緯度坐標(biāo)信息:

到此這篇關(guān)于Python實(shí)現(xiàn)常見(jiàn)的4種坐標(biāo)互相轉(zhuǎn)換的文章就介紹到這了,更多相關(guān)Python坐標(biāo)轉(zhuǎn)換內(nèi)容請(qǐng)搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家!

相關(guān)文章

  • python利用線程生成不同尺寸的縮略圖實(shí)例詳解

    python利用線程生成不同尺寸的縮略圖實(shí)例詳解

    這篇文章主要介紹了python利用線程生成不同尺寸的縮略圖,本文通過(guò)實(shí)例代碼給大家介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)或工作具有一定的參考借鑒價(jià)值,需要的朋友可以參考下
    2022-05-05
  • python2利用wxpython生成投影界面工具的圖文詳解

    python2利用wxpython生成投影界面工具的圖文詳解

    這篇文章主要介紹了python2利用wxpython生成投影界面工具的圖文詳解,本文通過(guò)實(shí)例圖文相結(jié)合給大家介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)或工作具有一定的參考借鑒價(jià)值,需要的朋友可以參考下
    2021-04-04
  • Python中Sorted()函數(shù)的key參數(shù)使用方法詳解

    Python中Sorted()函數(shù)的key參數(shù)使用方法詳解

    這篇文章主要介紹了關(guān)于Python中Sorted()函數(shù)的key參數(shù)使用方法 ,文中通過(guò)示例代碼介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下
    2019-06-06
  • Pycharm創(chuàng)建項(xiàng)目時(shí)如何自動(dòng)添加頭部信息

    Pycharm創(chuàng)建項(xiàng)目時(shí)如何自動(dòng)添加頭部信息

    這篇文章主要介紹了Pycharm創(chuàng)建項(xiàng)目時(shí) 自動(dòng)添加頭部信息,需要的朋友可以參考下
    2019-11-11
  • 使用實(shí)現(xiàn)pandas讀取csv文件指定的前幾行

    使用實(shí)現(xiàn)pandas讀取csv文件指定的前幾行

    下面小編就為大家分享一篇使用實(shí)現(xiàn)pandas讀取csv文件指定的前幾行,具有很好的參考價(jià)值,希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧
    2018-04-04
  • Pandas使用分隔符或正則表達(dá)式將字符串拆分為多列

    Pandas使用分隔符或正則表達(dá)式將字符串拆分為多列

    本文主要介紹了Pandas使用分隔符或正則表達(dá)式將字符串拆分為多列,文中通過(guò)示例代碼介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友們下面隨著小編來(lái)一起學(xué)習(xí)學(xué)習(xí)吧
    2023-02-02
  • python編程Flask框架簡(jiǎn)單使用教程

    python編程Flask框架簡(jiǎn)單使用教程

    這篇文章主要為大家介紹了python編程中Flask框架簡(jiǎn)單使用教程,有需要的朋友可以借鑒參考下希望能夠有所幫助,祝大家多多進(jìn)步早日升職加薪
    2021-11-11
  • 最新評(píng)論