Python坐标转换并保存文件

最新推荐文章于 2024-03-23 16:33:48 发布

不秃头的LT

最新推荐文章于 2024-03-23 16:33:48 发布

阅读量727

点赞数 1

分类专栏：笔记文章标签： python excel

本文链接：https://blog.csdn.net/firstlt0217/article/details/108355212

版权

笔记专栏收录该内容

15 篇文章 2 订阅

订阅专栏

运行前需先导入：

import math
import xlrd #excle模块
import xlwt
安装方式：pip install *

import math
import xlrd    #excle模块
import xlwt
class LngLatTransfer():

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

    def GCJ02_to_BD09(self, gcj_lng, gcj_lat):
        # """
        # 实现GCJ02向BD09坐标系的转换
        # :param lng: GCJ02坐标系下的经度
        # :param lat: GCJ02坐标系下的纬度
        # :return: 转换后的BD09下经纬度
        # """
        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):
        # """
        # 实现BD09坐标系向GCJ02坐标系的转换
        # :param bd_lng: BD09坐标系下的经度
        # :param bd_lat: BD09坐标系下的纬度
        # :return: 转换后的GCJ02下经纬度
        # """
        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):
        # """
        # 实现WGS84坐标系向GCJ02坐标系的转换
        # :param lng: WGS84坐标系下的经度
        # :param lat: WGS84坐标系下的纬度
        # :return: 转换后的GCJ02下经纬度
        # """
        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):
        # """
        # 实现GCJ02坐标系向WGS84坐标系的转换
        # :param gcj_lng: GCJ02坐标系下的经度
        # :param gcj_lat: GCJ02坐标系下的纬度
        # :return: 转换后的WGS84下经纬度
        # """
        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):
        # """
        # 实现BD09坐标系向WGS84坐标系的转换
        # :param bd_lng: BD09坐标系下的经度
        # :param bd_lat: BD09坐标系下的纬度
        # :return: 转换后的WGS84下经纬度
        # """
        lng, lat = self.BD09_to_GCJ02(bd_lng, bd_lat)
        return self.GCJ02_to_WGS84(lng, lat)


    def WGS84_to_BD09(self, lng, lat):
        # """
        # 实现WGS84坐标系向BD09坐标系的转换
        # :param lng: WGS84坐标系下的经度
        # :param lat: WGS84坐标系下的纬度
        # :return: 转换后的BD09下经纬度
        # """
        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):
        # """
        # 实现WGS84向web墨卡托的转换
        # :param lng: WGS84经度
        # :param lat: WGS84纬度
        # :return: 转换后的web墨卡托坐标
        # """
        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):
        # """
        # 实现web墨卡托向WGS84的转换
        # :param x: web墨卡托x坐标
        # :param y: web墨卡托y坐标
        # :return: 转换后的WGS84经纬度
        # """
        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

    def open_excel(self, excel_file):
        # '''
        #     读取excel函数
        #     args：excel_file（excel文件，目录在py文件同目录）
        #     returns：book
        # '''
        try:
            book = xlrd.open_workbook(excel_file)  # 文件名，把文件与py文件放在同一目录下
        except:
            print("open excel file failed!")
        try:
            sheet = book.sheet_by_name("Sheet1")  # execl里面的worksheet1
            return sheet
        except:
            print("locate worksheet in excel failed!")

    def set_style(self, name, height, bold=False):

        style = xlwt.XFStyle()  # 初始化样式

        font = xlwt.Font()  # 创建字体, Font定义字体的大小、颜色
        font.name = name  # 字体名称
        font.bold = bold  # 粗体
        # font.italic = True  #斜体
        font.height = height  # 字体大小
        # font.colour_index = 4  # 字体颜色

        borders = xlwt.Borders()  # 单元格边框线
        borders.left = xlwt.Borders.THIN  # 设置边框线粗细
        borders.right = xlwt.Borders.THIN
        borders.top = xlwt.Borders.THIN
        borders.bottom = xlwt.Borders.THIN
        # borders.right_colour = 1  #设置边框线颜色
        # borders.left_colour = 1
        # borders.top_colour = 1
        # borders.bottom_colour = 1

        alignment = xlwt.Alignment()  # 设置字体在单元格中的位置
        alignment.horz = xlwt.Alignment.HORZ_CENTER  # 水平居左
        alignment.vert = xlwt.Alignment.VERT_CENTER  # 垂直居中

        pat = xlwt.Pattern()  # 设置单元格背景颜色
        pat.pattern = xlwt.Pattern.SOLID_PATTERN
        pat.pattern_fore_colour = 5  # 黄色

        style.font = font
        style.borders = borders
        style.alignment = alignment
        style.pattern = pat

        return style

    def write_file(self, data):
        book = xlwt.Workbook(encoding='utf-8') #创建Workbook，相当于创建Excel

        # 创建sheet，Sheet1为表的名字，cell_overwrite_ok为是否覆盖单元格
        sheet1 = book.add_sheet(u'Sheet1', cell_overwrite_ok=True)

        # 向表中添加数据
        # sheet1.write(0, 0, 'Englishname')  # 第0行第0列
        # sheet1.write(1, 0, 'Hellen')  # 第一行第0列
        # sheet1.write(0, 1, '中文名字')
        # sheet1.write(1, 1, '海伦')
        r = 0
        for i, j in data.items():  # i表示data中的key，j表示data中的value
            le = len(j)  # values返回的列表长度
            if r == 0:
                sheet1.write(r, 0, i, self.set_style('Arial', 220, True))  # 添加第 0 行 0 列数据单元格背景设为黄色
            else:
                sheet1.write(r, 0, i)  # 添加第 1 列的数据

            for c in range(1, le + 1):  # values列表中索引
                if r == 0:
                    sheet1.write(r, c, j[c - 1], self.set_style('Arial', 220, True))  # 添加第 0 行，2 列到第 5 列的数据单元格背景设为黄色
                else:
                    sheet1.write(r, c, j[c - 1])

            r += 1  # 行数
        book.save('write.xls')


try:
    lt = LngLatTransfer()
    result = lt.BD09_to_GCJ02(119, 40)
    sheet = lt.open_excel('坐标.xlsx')
    print(result)
    row_num = sheet.nrows
    list=[]
    data={}
    data[0] = ["序号", "bdlng", "bdlat", "wgs84lng", "wgs84lat"]
    for i in range(1, row_num):  # 第一行是标题名，对应表中的字段名所以应该从第二行开始，计算机以0开始计数，所以值是1
        row_data = sheet.row_values(i)  # 按行获取excel的值
        value = (row_data[0], row_data[1])
        result = lt.BD09_to_GCJ02(row_data[0], row_data[1])
        # print(i, result[0])
        data[i]=[i, row_data[0], row_data[0], result[0], result[0]]
        # print(i, data)
    print(data)
    lt.write_file(data)
except:
    print("失败")

不秃头的LT

关注

1
点赞
踩
4

收藏

觉得还不错? 一键收藏
打赏
0
评论
Python坐标转换并保存文件

var GPS = { PI : 3.14159265358979324, x_pi : 3.14159265358979324 * 3000.0 / 180.0, delta : function (lat, lon) { // Krasovsky 1940 // // a = 6378245.0, 1/f = 298.3 // b = a * (1 - f) // ee = (a^2 - b^2) /
复制链接

扫一扫