脚本之家 服务器常用软件
快捷导航
您的位置：首页 → 脚本专栏 → python → Python GDAL坐标投影转换

使用Python和GDAL给图片加坐标系的实现思路(坐标投影转换)

 更新时间：2021年03月15日 10:28:52   作者：qq_41440099  
这篇文章主要介绍了使用Python和GDAL给图片加坐标系的实现思路(坐标投影转换),本文给大家介绍的非常详细，对大家的学习或工作具有一定的参考借鉴价值，需要的朋友可以参考下

**

使用Python和GDAL给图片加坐标系

**
假设你已经知道arcgis地理配准（如下图内容），懂一点python。

在这里插入图片描述

**

-目的和背景

1.从地图网站获得一张PNG格式的截图，已知坐标系为WGS84和左上角坐标。arcgis地理配准再定义投影即可给它加上原图的坐标系。
2.假设有上千张图片，可用Python和GDAL给图片加坐标系。

-实现思路

1.使用GDAL需要知道待投影图片的地理坐标信息、仿射矩阵参数。

仿射矩阵参数是干什么的？见：https://zhuanlan.zhihu.com/p/72184440
主要含义：
1）不同坐标系的转换，涉及到仿射变换，又称仿射映射，是指在几何中，一个向量空间进行一次线性变换并接上一个平移，变换为另一个向量空间。
2）仿射矩阵信息有六个参数，描述的是栅格行列号和地理坐标之间的关系：
‘''
0：图像左上角的X坐标；
1：图像东西方向分辨率；
2：旋转角度，如果图像北方朝上，该值为0；
3：图像左上角的Y坐标；
4：旋转角度，如果图像北方朝上，该值为0；
5：图像南北方向分辨率；
‘''

2.在arcgis使用一张图片和三个角点的坐标进行地理配准，再定义投影完成坐标转换。

使用下面的代码获取仿射矩阵和投影参数：

dataset = gdal.Open('a.png')
 print （dataset.GetGeoTransform()）#仿射矩阵
 print （dataset.GetProjection()）#地图投影信息
# 打印结果为：
# (116.33333, 8.321688443e-05, 0.0, 39.99999, 0.0, -6.223016769e-05)
# 'GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84", 6378137, 298.257223563, AUTHORITY["EPSG", "7030"]], AUTHORITY["EPSG", "6326"]], PRIMEM["Greenwich", 0], UNIT["degree", 0.0174532925199433], AUTHORITY["EPSG", "4326"]]'

3.批量获取图片的仿射矩阵

# coors是用来存储图片对应左上角坐标的字典。格式为{‘a.png‘'：[116.33333,39.6]，}
  image_list = os.listdir('D:\\dd')
  image_num = len(image_list)
  for k in range(image_num):
    if image_list[k].endswith('.png'):
      img_name = img_none_path + '/' + image_list[k]
      img_pos_transf = (float(coors[image_list[k]][0]), 8.321688443e-05,
               0.0, float(coors[image_list[k]][1]), 0.0, -6.223016769e-05)#根据第二步获得像元分辨率和投影
      print(img_pos_transf)
      img_pos_proj = 'GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84", 6378137, 298.257223563, AUTHORITY["EPSG", "7030"]], AUTHORITY["EPSG", "6326"]], PRIMEM["Greenwich", 0], UNIT["degree", 0.0174532925199433], AUTHORITY["EPSG", "4326"]]'
      def_geoCoordSys(img_name, img_pos_transf, img_pos_proj)#坐标转换的函数

4.给图片加坐标系的主要函数如下

来自文章 ：https://blog.csdn.net/nominior/article/details/102737294

def def_geoCoordSys(read_path, img_transf, img_proj):
    array_dataset = gdal.Open(read_path)
    img_array = array_dataset.ReadAsArray(
      0, 0, array_dataset.RasterXSize, array_dataset.RasterYSize)
    if 'int8' in img_array.dtype.name:
      datatype = gdal.GDT_Byte
    elif 'int16' in img_array.dtype.name:
      datatype = gdal.GDT_UInt16
    else:
      datatype = gdal.GDT_Float32

    if len(img_array.shape) == 3:
      img_bands, im_height, im_width = img_array.shape
    else:
      img_bands, (im_height, im_width) = 1, img_array.shape

    filename = read_path[:-4] + '_proj' + '.tif'
    driver = gdal.GetDriverByName("GTiff") # 创建文件驱动
    dataset = driver.Create(
      filename, im_width, im_height, img_bands, datatype)
    dataset.SetGeoTransform(img_transf) # 写入仿射变换参数
    dataset.SetProjection(img_proj) # 写入投影

    # 写入影像数据
    if img_bands == 1:
      dataset.GetRasterBand(1).WriteArray(img_array)
    else:
      for i in range(img_bands):
        dataset.GetRasterBand(i + 1).WriteArray(img_array[i])
    print(read_path, 'geoCoordSys get!')

到此这篇关于使用Python和GDAL给图片加坐标系的实现思路(坐标投影转换)的文章就介绍到这了,更多相关Python GDAL坐标投影转换内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家！

您可能感兴趣的文章:

相关文章

  • 如何在pycharm中快捷安装pip命令(如pygame)

    如何在pycharm中快捷安装pip命令(如pygame)

    这篇文章主要介绍了如何在pycharm中快捷安装pip命令(如pygame)，具有很好的参考价值，希望对大家有所帮助。
    2021-05-05
  • python3 如何解压缩.gz文件

    python3 如何解压缩.gz文件

    这篇文章主要介绍了python3 如何解压缩.gz文件，具有很好的参考价值，希望对大家有所帮助。如有错误或未考虑完全的地方，望不吝赐教
    2022-09-09
  • 详解利用上下文管理器扩展Python计时器

    详解利用上下文管理器扩展Python计时器

    本文将和大家一起了解什么是上下文管理器 和 Python 的 with 语句，以及如何完成自定义。然后扩展 Timer 以便它也可以用作上下文管理器，感兴趣的可以了解一下
    2022-06-06
  • Python Switch Case三种实现方法代码实例

    Python Switch Case三种实现方法代码实例

    这篇文章主要介绍了Python Switch Case2种实现方法代码实例,文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下
    2020-06-06
  • 如何基于python操作json文件获取内容

    如何基于python操作json文件获取内容

    这篇文章主要介绍了如何基于python操作json文件获取内容,文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下
    2019-12-12
  • python访问sqlserver示例

    python访问sqlserver示例

    这篇文章主要介绍了python访问sqlserver示例,需要的朋友可以参考下
    2014-02-02
  • python itchat给指定联系人发消息的方法

    python itchat给指定联系人发消息的方法

    这篇文章主要介绍了python itchat给指定联系人发消息的方法，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧
    2019-06-06
  • Python将多个list合并为1个list的方法

    Python将多个list合并为1个list的方法

    今天小编就为大家分享一篇Python将多个list合并为1个list的方法，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
    2018-06-06
  • pytest多线程与多设备并发appium

    pytest多线程与多设备并发appium

    这篇文章介绍了pytest多线程与多设备并发appium，文中通过示例代码介绍的非常详细。对大家的学习或工作具有一定的参考借鉴价值，需要的朋友可以参考下
    2022-06-06
  • 解决python super()调用多重继承函数的问题

    解决python super()调用多重继承函数的问题

    今天小编就为大家分享一篇解决python super()调用多重继承函数的问题，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
    2019-06-06

最新评论