【Arcpy学习实践教程】wgs84坐标系和火星坐标系的转换中demo的对与错

最新推荐文章于 2024-05-06 18:18:34 发布
最新推荐文章于 2024-05-06 18:18:34 发布
阅读量470 收藏 1
点赞数
分类专栏： arcpy
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/ldl250058/article/details/104345662
版权

度娘和谷哥已经变成了我们学习工作生活中必不可少的工具。

更有甚者，甚至已经不用输入法来搜索，而直接通过语音识别来搜索。但是我们搜索的结果真的可靠？我们在找到我们想要的资源之后是否有认真检验一下，我们找到的代码，找到的资料是否正确。当没有思考力和鉴别力的搬运工进入大众视野时，我们就需要谨慎起来了。

最近，因为工作的原因需要对高德坐标（即火星坐标）和wgs84坐标系实现互转。一直记得有位大神曾经在网上发布过相关的代码，和大家一样，第一步做的也是搜索，信息铺天盖地，随便找几个阅读量大的原创文章，给大家展示一下代码吧。

1.在《博客园》找到了一篇文章《GCJ-02火星坐标系和WGS-84坐标系转换关系》（ps：一般喜欢在博客园里面找代码，因为这边的文章相对专业一些，附上链接，下面是上面贴的code，python哈）

# -*- coding: utf-8 -*-
import json
import math

x_pi = 3.14159265358979324 * 3000.0 / 180.0
pi = 3.1415926535897932384626  # π
a = 6378245.0  # 长半轴
ee = 0.00669342162296594323  # 扁率

def wgs84togcj02(lng, lat):
    """
    WGS84转GCJ02(火星坐标系)
    :param lng:WGS84坐标系的经度
    :param lat:WGS84坐标系的纬度
    :return:
    """
    if out_of_china(lng, lat):  # 判断是否在国内
        return lng, lat
    dlat = transformlat(lng - 105.0, lat - 35.0)
    dlng = transformlng(lng - 105.0, lat - 35.0)
    radlat = lat / 180.0 * pi
    magic = math.sin(radlat)
    magic = 1 - ee * magic * magic
    sqrtmagic = math.sqrt(magic)
    dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * pi)
    dlng = (dlng * 180.0) / (a / sqrtmagic * math.cos(radlat) * pi)
    mglat = lat + dlat
    mglng = lng + dlng
    return [mglng, mglat]


def gcj02towgs84(lng, lat):
    """
    GCJ02(火星坐标系)转GPS84
    :param lng:火星坐标系的经度
    :param lat:火星坐标系纬度
    :return:
    """
    if out_of_china(lng, lat):
        return lng, lat
    dlat = transformlat(lng - 105.0, lat - 35.0)
    dlng = transformlng(lng - 105.0, lat - 35.0)
    radlat = lat / 180.0 * pi
    magic = math.sin(radlat)
    magic = 1 - ee * magic * magic
    sqrtmagic = math.sqrt(magic)
    dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * pi)
    dlng = (dlng * 180.0) / (a / sqrtmagic * math.cos(radlat) * pi)
    mglat = lat + dlat
    mglng = lng + dlng
    return [lng * 2 - mglng, lat * 2 - mglat]


def transformlat(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 * pi) + 20.0 *
            math.sin(2.0 * lng * pi)) * 2.0 / 3.0
    ret += (20.0 * math.sin(lat * pi) + 40.0 *
            math.sin(lat / 3.0 * pi)) * 2.0 / 3.0
    ret += (160.0 * math.sin(lat / 12.0 * pi) + 320 *
            math.sin(lat * pi / 30.0)) * 2.0 / 3.0
    return ret


def transformlng(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 * pi) + 20.0 *
            math.sin(2.0 * lng * pi)) * 2.0 / 3.0
    ret += (20.0 * math.sin(lng * pi) + 40.0 *
            math.sin(lng / 3.0 * pi)) * 2.0 / 3.0
    ret += (150.0 * math.sin(lng / 12.0 * pi) + 300.0 *
            math.sin(lng / 30.0 * pi)) * 2.0 / 3.0
    return ret


def out_of_china(lng, lat):
    """
    判断是否在国内，不在国内不做偏移
    :param lng:
    :param lat:
    :return:
    """
    if lng < 72.004 or lng > 137.8347:
        return True
    if lat < 0.8293 or lat > 55.8271:
        return True
    return False


if __name__ == '__main__':
    [lng,lat]=[114.061202,22.529388]
    [dstlng, dstlat] = gcj02towgs84(lng, lat)
    print(dstlng, dstlat)

2.在csdn中找到一篇2017年的文章，阅读量大概过万了吧，反正应该还是比较权威的吧，或者是说比较早的---《火星坐标、百度坐标、WGS84坐标转换代码（JS、python版）》，连接不可少，哈哈。代码也是不能少的，这个是必须的

# -*- coding: utf-8 -*-
import json
import requests
import math
 
key = 'your key here'  # 这里填写你的百度开放平台的key
x_pi = 3.14159265358979324 * 3000.0 / 180.0
pi = 3.1415926535897932384626  # π
a = 6378245.0  # 长半轴
ee = 0.00669342162296594323  # 扁率
 
 
def geocode(address):
    """
    利用百度geocoding服务解析地址获取位置坐标
    :param address:需要解析的地址
    :return:
    """
    geocoding = {'s': 'rsv3',
                 'key': key,
                 'city': '全国',
                 'address': address}
    res = requests.get(
        "http://restapi.amap.com/v3/geocode/geo", params=geocoding)
    if res.status_code == 200:
        json = res.json()
        status = json.get('status')
        count = json.get('count')
        if status == '1' and int(count) >= 1:
            geocodes = json.get('geocodes')[0]
            lng = float(geocodes.get('location').split(',')[0])
            lat = float(geocodes.get('location').split(',')[1])
            return [lng, lat]
        else:
            return None
    else:
        return None
 
 
def gcj02tobd09(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]
 
 
def bd09togcj02(bd_lon, bd_lat):
    """
    百度坐标系(BD-09)转火星坐标系(GCJ-02)
    百度——>谷歌、高德
    :param bd_lat:百度坐标纬度
    :param bd_lon:百度坐标经度
    :return:转换后的坐标列表形式
    """
    x = bd_lon - 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)
    gg_lng = z * math.cos(theta)
    gg_lat = z * math.sin(theta)
    return [gg_lng, gg_lat]
 
 
def wgs84togcj02(lng, lat):
    """
    WGS84转GCJ02(火星坐标系)
    :param lng:WGS84坐标系的经度
    :param lat:WGS84坐标系的纬度
    :return:
    """
    if out_of_china(lng, lat):  # 判断是否在国内
        return lng, lat
    dlat = transformlat(lng - 105.0, lat - 35.0)
    dlng = transformlng(lng - 105.0, lat - 35.0)
    radlat = lat / 180.0 * pi
    magic = math.sin(radlat)
    magic = 1 - ee * magic * magic
    sqrtmagic = math.sqrt(magic)
    dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * pi)
    dlng = (dlng * 180.0) / (a / sqrtmagic * math.cos(radlat) * pi)
    mglat = lat + dlat
    mglng = lng + dlng
    return [mglng, mglat]
 
 
def gcj02towgs84(lng, lat):
    """
    GCJ02(火星坐标系)转GPS84
    :param lng:火星坐标系的经度
    :param lat:火星坐标系纬度
    :return:
    """
    if out_of_china(lng, lat):
        return lng, lat
    dlat = transformlat(lng - 105.0, lat - 35.0)
    dlng = transformlng(lng - 105.0, lat - 35.0)
    radlat = lat / 180.0 * pi
    magic = math.sin(radlat)
    magic = 1 - ee * magic * magic
    sqrtmagic = math.sqrt(magic)
    dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * pi)
    dlng = (dlng * 180.0) / (a / sqrtmagic * math.cos(radlat) * pi)
    mglat = lat + dlat
    mglng = lng + dlng
    return [lng * 2 - mglng, lat * 2 - mglat]
 
 
def transformlat(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 * pi) + 20.0 *
            math.sin(2.0 * lng * pi)) * 2.0 / 3.0
    ret += (20.0 * math.sin(lat * pi) + 40.0 *
            math.sin(lat / 3.0 * pi)) * 2.0 / 3.0
    ret += (160.0 * math.sin(lat / 12.0 * pi) + 320 *
            math.sin(lat * pi / 30.0)) * 2.0 / 3.0
    return ret
 
 
def transformlng(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 * pi) + 20.0 *
            math.sin(2.0 * lng * pi)) * 2.0 / 3.0
    ret += (20.0 * math.sin(lng * pi) + 40.0 *
            math.sin(lng / 3.0 * pi)) * 2.0 / 3.0
    ret += (150.0 * math.sin(lng / 12.0 * pi) + 300.0 *
            math.sin(lng / 30.0 * pi)) * 2.0 / 3.0
    return ret
 
 
def out_of_china(lng, lat):
    """
    判断是否在国内，不在国内不做偏移
    :param lng:
    :param lat:
    :return:
    """
    if lng < 72.004 or lng > 137.8347:
        return True
    if lat < 0.8293 or lat > 55.8271:
        return True
    return False
 
 
if __name__ == '__main__':
    lng = 128.543
    lat = 37.065
    result1 = gcj02tobd09(lng, lat)
    result2 = bd09togcj02(lng, lat)
    result3 = wgs84togcj02(lng, lat)
    result4 = gcj02towgs84(lng, lat)
    result5 = geocode('北京市朝阳区朝阳公园')
    print result1, result2, result3, result4, result5

对比了一下，虽然写法有些许不一样，但是算法核心都是一样的。看到这么多的代码都是一样的，我当然也一样照做了。和第一个方法一样，我也做了定位，在高德地图上面发现了一个小问题，其实转出来的结果和实际有非常小的偏差。根据转出的结果，又用代码转回到wgs84，做投影，神奇的事情发生了，相差2m多。

（转回去之后的wgs84投影差别2m）

 

这个时候意识到一个问题，一定是这个算法存在问题。网上继续检索查找，大神浮出水面，geohey大神。用了qgis上他开发的转化工具之后，逆转换没有任何问题，这个时候正常的人也只会做一个事了，趴代码。代码趴到手，这里我也做一下搬运工吧。

# -*- coding: utf-8 -*-
##########################################################################################
"""
/***************************************************************************
 OffsetWGS84Core
                                 A QGIS plugin
 Class with methods for geometry and attributes processing
                              -------------------
        begin                : 2016-10-11
        git sha              : $Format:%H$
        copyright            : (C) 2017 by sshuair
        email                : sshuair@gmail.com
 ***************************************************************************/

/***************************************************************************
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 ***************************************************************************/
"""
from __future__ import print_function
##########################################################################################
from builtins import zip
import math
from math import sin, cos, sqrt, fabs, atan2
from math import pi as PI
# from numba import jit


# =================================================sshuair=============================================================
# define ellipsoid
a = 6378245.0
f = 1 / 298.3
b = a * (1 - f)
ee = 1 - (b * b) / (a * a)

# check if the point in china
def outOfChina(lng, lat):
    return not (72.004 <= lng <= 137.8347 and 0.8293 <= lat <= 55.8271)

# @jit
def geohey_transformLat(x, y):
    ret = -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y + 0.1 * x * y + 0.2 * sqrt(fabs(x))
    ret = ret + (20.0 * sin(6.0 * x * PI) + 20.0 * sin(2.0 * x * PI)) * 2.0 / 3.0
    ret = ret + (20.0 * sin(y * PI) + 40.0 * sin(y / 3.0 * PI)) * 2.0 / 3.0
    ret = ret + (160.0 * sin(y / 12.0 * PI) + 320.0 * sin(y * PI / 30.0)) * 2.0 / 3.0
    return ret

# @jit
def geohey_transformLon(x, y):
    ret = 300.0 + x + 2.0 * y + 0.1 * x * x +  0.1 * x * y + 0.1 * sqrt(fabs(x))
    ret = ret + (20.0 * sin(6.0 * x * PI) + 20.0 * sin(2.0 * x * PI)) * 2.0 / 3.0
    ret = ret + (20.0 * sin(x * PI) + 40.0 * sin(x / 3.0 * PI)) * 2.0 / 3.0
    ret = ret + (150.0 * sin(x / 12.0 * PI) + 300.0 * sin(x * PI / 30.0)) * 2.0 / 3.0
    return ret

# @jit
def wgs2gcj(wgsLon, wgsLat):
    if outOfChina(wgsLon, wgsLat):
        return wgsLon, wgsLat
    dLat = geohey_transformLat(wgsLon - 105.0, wgsLat - 35.0)
    dLon = geohey_transformLon(wgsLon - 105.0, wgsLat - 35.0)
    radLat = wgsLat / 180.0 * PI
    magic = math.sin(radLat)
    magic = 1 - ee * magic * magic
    sqrtMagic = sqrt(magic)
    dLat = (dLat * 180.0) / ((a * (1 - ee)) / (magic * sqrtMagic) * PI)
    dLon = (dLon * 180.0) / (a / sqrtMagic * cos(radLat) * PI)
    gcjLat = wgsLat + dLat
    gcjLon = wgsLon + dLon
    return (gcjLon, gcjLat)


def gcj2wgs(gcjLon, gcjLat):
    g0 = (gcjLon, gcjLat)
    w0 = g0
    g1 = wgs2gcj(w0[0], w0[1])
    # w1 = w0 - (g1 - g0)
    w1 = tuple([x[0]-(x[1]-x[2]) for x in zip(w0,g1,g0)])  
    # delta = w1 - w0
    delta = tuple([x[0] - x[1] for x in zip(w1, w0)])
    while (abs(delta[0]) >= 1e-6 or abs(delta[1]) >= 1e-6):
        w0 = w1
        g1 = wgs2gcj(w0[0], w0[1])
        # w1 = w0 - (g1 - g0)
        w1 = tuple([x[0]-(x[1]-x[2]) for x in zip(w0,g1,g0)])
        # delta = w1 - w0
        delta = tuple([x[0] - x[1] for x in zip(w1, w0)])
    return w1


def gcj2bd(gcjLon, gcjLat):
    z = sqrt(gcjLon * gcjLon + gcjLat * gcjLat) + 0.00002 * sin(gcjLat * PI * 3000.0 / 180.0)
    theta = atan2(gcjLat, gcjLon) + 0.000003 * cos(gcjLon * PI * 3000.0 / 180.0)
    bdLon = z * cos(theta) + 0.0065
    bdLat = z * sin(theta) + 0.006
    return (bdLon, bdLat)


def bd2gcj(bdLon, bdLat):
    x = bdLon - 0.0065
    y = bdLat - 0.006
    z = sqrt(x * x + y * y) - 0.00002 * sin(y * PI * 3000.0 / 180.0)
    theta = atan2(y, x) - 0.000003 * cos(x * PI * 3000.0 / 180.0)
    gcjLon = z * cos(theta)
    gcjLat = z * sin(theta)
    return (gcjLon, gcjLat)


def wgs2bd(wgsLon, wgsLat):
    gcj = wgs2gcj(wgsLon, wgsLat)
    return gcj2bd(gcj[0], gcj[1])


def bd2wgs(bdLon, bdLat):
    gcj = bd2gcj(bdLon, bdLat)
    return gcj2wgs(gcj[0], gcj[1])


if __name__ == '__main__':
    # wgs2gcj
    # coord = (112, 40)
    # trans = WGS2GCJ()
    print(wgs2gcj(112, 40))
    print(gcj2wgs(112.00678230985764, 40.00112245823686))

    # gcj2wgs

很明显在火星坐标系转wgs84坐标系过程中的算法存在较大的差异，这是之前误差2m的绝对原因。

希望以后大家在搬运代码的时候回过去验证一下，做一个有判断力和有辨识力的搬运工，我们虽然可以不重复生产轮子，但是我们不能够把别人的方轮子拿来用吧！！！！

 

本文由元凿坊独创，如需转载，请获得作者的授权。如发现问题，请及时告知作者，

欢迎前来探讨交流！！！！

发布于 17:28

桃花浪子立
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
地球坐标系(WGS - 84)转火星坐标系(GCJ)
superbeyone
04-14 2835
地球坐标系(WGS - 84)转火星坐标系(GCJ) package com.tdt.convert.utils; import java.math.BigDecimal; /** * @author Mr.superbeyone * @project coordinate-convert-tools * @className LngLatConvertUtil * @descript...
如何实现wgs84(ep4326)坐标系与百度坐标系，火星坐标系的互转，java与JavaScript都有
weixin_42759398的博客
04-03 850
国家规定：国内出版的各种地图系统(包括电子形式)，必须至少采用GCJ-02对地理位置进行首次加密(高德地图、腾讯地图等)；工作中用到百度地图和高德地图服务，必须实现转换就成了难题了，现在它来了，java后台和纯js方法都有，直接上代码。众所周知国内的地址服务不让直接用wgs84，百度使用的是百度坐标系，高德地图使用的是火星坐标系。国际地图提供商使用的坐标系(谷歌地图国外、osm、mapbox)。百度标准，百度 SDK，百度地图，百度GeoCoding 使用；2.国测局坐标系(GCJ-02、火星坐标系)
Python ArcPy批量将大量栅格文件的投影坐标系转为地理坐标系
最新发布
疯狂学习GIS的博客
05-06 658
本文介绍基于Python语言中的ArcPy模块，批量将多个遥感影像由投影坐标系转为地理坐标系的方法~
火星坐标系，百度坐标系，wgs84坐标系转换算法
CUIT_IGIS 学而不精
11-21 9634
一 坐标介绍 1 wgs84坐标系：是国际上tongy
WGS84地球坐标系，GCJ02火星坐标系，BD09百度坐标系简介与转换(1)
2401_83915450的博客
04-18 679
System.out.println("火星坐标系 : " + wgs84ToGcj02);System.out.println("地球坐标系 : " + wgs84ToGcj02);System.out.println("百度坐标系 : " + wgs84ToBd09);System.out.println("百度坐标系 : " + wgs84ToBd09);System.out.println("地球坐标系 : " + point);// GCJ02=>WGS84 火星坐标系=>地球坐标系(粗略)
地球坐标系 (WGS-84) 到火星坐标系 (GCJ-02) 的转换算法
热门推荐
A Fork of Stuffs
03-18 11万+
原文 WGS-84 到 GCJ-02 的转换（即 GPS 加偏）算法是一个普通青年轻易无法接触到的“公开”的秘密。这个算法的代码在互联网上是公开的，详情请使用 Google 搜索 "wgtochina_lb" 。 整理后的算法代码请参考 https://on4wp7.codeplex.com/SourceControl/changeset/view/21483#353936 。知道了这个算
基于ArcPy Python工具箱 百度坐标、国测局坐标和WGS84坐标 空间图层相互转换
06-17
2. **GCJ-02坐标系（国测局坐标）**：中国国家测绘局采用的坐标系，也称为火星坐标，对WGS84进行了加密处理，以保护国家安全。 3. **百度坐标系**：百度地图使用的坐标系，是在GCJ-02基础上再次加密的坐标系统，...
基于arcgis影像设置坐标系和修改投影处理的工具箱
03-09
arcgis10.1、10.2工具箱，支持tif影像文件的坐标系批量设置，以及批量修改投影坐标系，例如加代号。
给img、tif影像添加坐标系
07-10
用GDAL给img影像添加坐标，并转化成tif
CGCS2000坐标系，地理坐标+投影坐标（3度+6度），包括：xml和prj两种格式
08-30
这个坐标系的建立基于地球椭球模型，采用了WGS84（世界大地坐标系）的几何参数，但在参考框架上考虑了中国的地壳运动情况，更加适应我国的地理特点。 CGCS2000坐标系的核心概念包括地理坐标和投影坐标。地理坐标...
坐标转换（支持百度坐标系，火星坐标系，WGS84之间的互转）
02-27
坐标转换小工具，支持百度坐标系，火星坐标系，WGS84之间的互转，批量转换，非常实用的小工具，对于GIS互联网行业地图很有帮助
（下载链接在博客）ArcGIS工具箱实现百度、火星与WGS84坐标互转（Excel版）
07-25
百度坐标、火星坐标、WGS84坐标转换工具全部升级为脱离ArcGIS环境，可在本人博客文末找到链接试用，链接呈上https://blog.csdn.net/yhiki/article/details/86476276
wgs84的geojson数据转火星坐标
04-16
双击exe文件运行，geojson源文件要wgs84的地理坐标数据，在exe目录下建立before文件夹存放源文件，生成的after文件夹是存放转换后的文件
js中gmt时间转换为普通时间格式
xvzhengyang的博客
12-01 784
renderTime(date) { var dateee = new Date(date).toJSON(); return new Date(+new Date(dateee) + 8 * 3600 * 1000) .toISOString() .replace(/T/g, " ") .replace(/\.[\d]{3}Z/, ""); } 时间格式转换 譬如:"2020-11-30T13:09:58.000+...
WGS84(GPS)、火星坐标系（GCJ02）、百度地图（BD09）坐标系转换案例教程（附转换工具下载）
「 刘一哥与GIS的故事」
11-13 3741
在做基于百度地图、高德地图等电子地图做为地图服务的二次开发时，通常需要将具有WGS84等坐标的矢量数据（如行政区划、地名、河流、道路等GIS地理空间数据）添加到地图上面，然而，在线地图大多使用的是火星坐标系，需要事先将矢量数据转为火星坐标系。本文以案例的形式，讲述WGS84(GPS)、火星坐标系（GCJ02）、百度地图（BD09）坐标系之间互相转换。 文章目录一、WGS转火星坐标系对比二、火星坐标系简介三、转换过程四、转换工具下载 一、WGS转火星坐标系对比 文中将行政区的WGS坐标转为火星坐标系，局.
Arcpy脚本——基于Arcgis实现shp/gdb等到火星/百度坐标系的转换
水梦雪幻
08-29 2241
  在生产中，会遇到将我们既有的坐标系转换到高德或者百度坐标系下，从而可以使用百度高德的基础地图。亦或者将获得的百度高德数据转换到普通的84坐标系下。本文基于Arcpy脚本实现了对shp或者gdb数据的坐标系之间的相互转换。 1、从常用坐标系到百度高德坐标系的转换 思路如下： 1）读取原始数据的坐标系； 2）利用Arcpy将原始坐标转换成WGS84经纬度； 3）对于线和面，利用Arcpy对数据进行...
WGS84坐标 火星坐标 BD09坐标转换
u013771380的博客
07-12 439
// 定义一些常量 const x_PI = (3.14159265358979324 * 3000.0) / 180.0 const PI = 3.1415926535897932384626 const a = 6378245.0 const ee = 0.00669342162296594323 /** * 百度坐标系 (BD-09) 与 火星坐标系 (GCJ-02)的转换 / 即百度转谷歌、高德 * @param { Number } bd_lon * @param { Number } .
WGS84坐标系转火星坐标系
BinChasing的博客
06-21 3355
from __future__ import division from math import pi,sqrt,sin,cos# a python binding of https://on4wp7.codeplex.com/SourceControl/changeset/view/21483#353936# Krasovsky 1940 # # a = 6378245.0, 1/f = 298.
使用arcpy转换成wgs84坐标系代码
06-07
下面是使用arcpy将数据...最后调用 `arcpy.management.Project()` 方法进行坐标系转换，将输入数据转换为WGS84坐标系下的数据并输出到指定路径。 需要注意的是，上述代码中的路径和坐标系需要根据实际情况进行修改。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值