GDAL 地理坐标与图像行列坐标之间相互转换公式

最新推荐文章于 2024-04-11 22:45:32 发布
最新推荐文章于 2024-04-11 22:45:32 发布
阅读量9.1k 收藏 25
点赞数 2
CC 4.0 BY-SA版权
分类专栏： gdal 文章标签： gdal
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.youkuaiyun.com/zyf19930610/article/details/43487473
本文详细介绍了如何通过仿射变换获取地理坐标与行列坐标之间的转换方法，包括行列坐标转换为地理坐标以及地理坐标转换为行列坐标的步骤。

摘要生成于 C知道 ，由 DeepSeek-R1 满血版支持， 前往体验 >

通过GetGeoTransform(dGeoTrans)获取仿射信息

1,，行列坐标转换地理坐标

Xp = dGeoTrans [0] +Xpixel*dGeoTrans [1]+Yline*dGeoTrans [2];

Yp = dGeoTrans [3] + Xpixel*dGeoTrans [4] + YlineL*dGeoTrans[5];


2，地理坐标转换行列坐标

double dTemp = dGeoTrans[1] * dGeoTrans[5] - dGeoTrans[2] *dGeoTrans[4];

Xpixel= (dGeoTrans[5] * (Xp - dGeoTrans[0]) -dGeoTrans[2] * (Yp - dGeoTrans[3])) / dTemp + 0.5;

Yline = (dGeoTrans[1] * (Yp - dGeoTrans[3]) -dGeoTrans[4] * (Xp - dGeoTrans[0])) / dTemp + 0.5;


【遥感图像】使用GDAL进行影像投影坐标、地理坐标、图上坐标的转换
小苏打的学习博客
05-18 2637
使用GDAL进行影像投影坐标、地理坐标、图上坐标的转换 我使用GDAL库写了四个函数分别进行投影坐标与地理坐标（经纬度）之间的转换，投影坐标和图上坐标（行列号）之间的转换。有需要的朋友可以参考。 直接上代码吧，因为代码很简单（Python版本）。 # -*- encoding: utf-8 -*- from osgeo import gdal from osgeo import osr import numpy as np def getSRSPair(dataset): ''' 获得给
GDAL遥感影像像素坐标与地理（投影）坐标相互转换
小鱼的博客
01-26 2016
想直接用的代码，找不着了。。。只能翻出以前写的代码，记录一下 # 像素坐标转地理（投影）坐标 def pixel2geo(geotrans, px, py): ulX = geotrans[0] ulY = geotrans[3] pixelWidth = geotrans[1] pixelHeight = geotrans[5] rotate1 = geotrans[2] rotate2 = geotrans[4] Xgeo = ulX + px
gdal坐标变换（x,y变经纬度）
小松鼠的随感
06-30 1万+
gdal坐标变换最近因为一些原因需要读取地理数据，我使用的工具是gdal，在处理数据的时候遇到一些小的问题，自己折腾了会儿，现在分享出来，避免大家浪费时间。 tiff, geotif和gdal都可以提取tiff数据，对于地理数据gdal可能更方便一些 我使用了tiff,geotiff尝试提取过手头的tif数据，发现有一些tag是未定义的，会报warning，尽管也可以手工处理，但有些麻烦，gdal
GDAL 地理坐标与图像行列坐标之间相互转换公式
阿a楞的博客
06-04 3248
通过GetGeoTransform(dGeoTrans)获取仿射信息 1、行列坐标转换地理坐标 Xp = dGeoTrans [0] +Xpixel*dGeoTrans [1]+Yline*dGeoTrans [2]; Yp = dGeoTrans [3] + Xpixel*dGeoTrans [4] + Yline*dGeoTrans[5]; 2、地理坐标转换行列坐标 double dTem...
使用GDAL进行影像投影坐标、地理坐标、图上坐标的转换
weixin_44576131的博客
03-26 2987
使用GDAL库写了四个函数分别进行投影坐标与地理坐标（经纬度）之间的转换，投影坐标和图上坐标（行列号）之间的转换。 来源：https://cloud.tencent.com/developer/article/1386146 from osgeo import gdal from osgeo import osr import numpy as np def getSRSPair(dataset): ‘’’ 获得给定数据的投影参考系和地理参考系 :param dataset: GDAL地理数据 :retur
C++之GDAL读取tiff文件，获取行列坐标值
qq_35662333的博客
03-22 1293
原文链接：https://blog.youkuaiyun.com/weixin_43712770/article/details/122017317。tif中坐标计算的方法如下，其中Col表示该坐标点处图像的列号，ROW表示该坐标点处图像的行号。比如图像左上角Col为0，ROW为0，图像右下角Col为图像宽度，ROW为图像高度。
关于GDAL计算图像坐标的几个问题
qq_35662333的博客
04-11 1140
比如第一个，如何计算图像的四至范围，图像的四至范围从图1中可以看出，图像的四至其实就是图像左上角坐标和右下角的坐标为起来的矩形区域，那么就分别将左上角和右下角的行列号按照上面的公式进行转换即可得到四至范围；对于一个普通的标准图像来说（这里的标准图像是指GT(2)和GT(4)都为0），如图1所示，图像的行列号坐标为XOY，每个网格代表一个图像像素区域，i表示列号，j表示行号，淡蓝色右下角的行列坐标为（i，j），图中红色方块纵向长度为dy，横向长度为dx，分别为图像的分辨率；2：如何获取图像的四至范围？
gdal坐标转换总结(转换）
Orbita_wangtao
11-08 9128
gdal坐标转换总结(转换）   首先，在进行坐标转换之前，有必要先了解一下有关坐标系的几个基本概念。 地理坐标系（Geographic Coordinate Systems） 地理坐标系是一个球面的坐标系统，以经纬度为单位，它由椭球体和大地基准面两个部分组成。 椭球体（spheroid） 我们要将地理信息以球面坐标系的方式表达，首先需要找到一个可以量化计算的椭球体。一个椭球体的确定需...
gdal坐标转换总结
佳佳的博客
05-03 5506
首先，在进行坐标转换之前，有必要先了解一下有关坐标系的几个基本概念。 地理坐标系（Geographic Coordinate Systems） 地理坐标系是一个球面的坐标系统，以经纬度为单位，它由椭球体和大地基准面两个部分组成。 椭球体（spheroid） 我们要将地理信息以球面坐标系的方式表达，首先需要找到一个可以量化计算的椭球体。一个椭球体的确定需要以下参数：长半轴、短半轴、偏...
【Python栅格数据地理空间分析】：坐标转换与实时更新的秘密
[【Python栅格数据地理空间分析】：坐标转换与实时更新的秘密](http://deeplearning.lipingyang.org/wp-content/uploads/2019/03/img_5c9fa74c66f4e.png) # 1. Python栅格数据地理空间分析基础 ## 1.1 栅格数据概念...
【图像处理库介绍】GDAL库对栅格数据的支持与操作
GDAL（Geospatial Data Abstraction Library）是一个开源的翻译库，用于读取和写入栅格地理空间数据格式。它支持广泛的数据格式，并允许统一的访问，无需了解每个格式的内部结构。GDAL是GIS专业人员不可或缺的工具，...
GDAL中地理坐标与像素坐标之间的转换
qq_41595248的博客
07-04 244
遥感影像投影坐标系与地理坐标系之间的转换关系
Python遥感开发之使用GDAL进行影像投影坐标、地理坐标、图上坐标的转换
等待着冬天的风的博客
01-04 2203
前言：转载自《使用GDAL进行影像投影坐标、地理坐标、图上坐标的转换》 使用GDAL库写了四个函数分别进行投影坐标与地理坐标（经纬度）之间的转换，投影坐标和图上坐标（行列号）之间的转换。有需要的朋友可以参考。
影像经纬度与行列号互相转换
qq_35662333的博客
03-22 560
【代码】影像经纬度与行列号互相转换。
Gdal中GetGeoTransfrom的含义
schmiloo的博客
06-30 1万+
利用GetGeoTransfrom函数可以获得图像中任意一点的地理坐标。 示例代码如下： #include "gdal_priv.h" ... ... \\打开遥感影像文件 GDALDataset *ds = (GDALDataset*)GDALOpen(filename, GA_ReadOnly); \\geoTransform共包括6个参数 double geoTransform[6]; \\获取geoTransform参数 ds->GetGeoTransform(geoTransform);
使用GDAL进行坐标转换
alex1801
02-24 5057
我国的基本比例尺地形图（1：5千，1：1万，1：2.5万，1：10万，1：25万，1：50万，1：100万）中，大于或等于1：50万均采用高斯-克吕格投影（Gauss_Kruger），又叫横轴墨卡托投影（Transverse Mercator）；投影坐标系统是根据某种映射关系，将地理坐标系统中由经纬度确定的三维球面坐标投影到二维的平面上所使用的坐标系统。地图投影就是把地球表面的任意点，利用一定数学法则，转换到地图平面上的理论和方法，一般有两种坐标系来进行表示，分别是地理坐标系和投影坐标系。
GDAL综合整理--6：GDAL部分类说明
热门推荐
rrrrssss00的专栏
04-12 2万+
在网上看到这篇文章写得非常全面，而且是中文的，就转过来了，学习一下 其中有部分内容是重复的，不过不影响，非常好的资料 主要转自http://blog.youkuaiyun.com/liminlu0314?viewmode=contents GDALMajorObject 所有核心类的父类是GDALMajorObject，它定义了一些操作元数据的属性和方法供子类继承。 元数据：
GDAL获取tif的像元大小
最新发布
08-02
<think>我们正在使用GDAL来处理TIFF文件，用户想要获取TIFF图像文件的像元大小（即每个像素代表的地理范围）。 在遥感或地理信息系统中，像元大小通常指的是一个像素在现实世界中的大小，例如每个像素代表10米×10米。 在GDAL中，我们可以通过读取图像的地理变换（GeoTransform）来获取像元大小。 地理变换是一个包含6个元素的数组，其定义如下： adfGeoTransform[0] /* 左上角X坐标 */ adfGeoTransform[1] /* 东西方向像元大小（宽度） */ adfGeoTransform[2] /* 旋转系数，通常为0 */ adfGeoTransform[3] /* 左上角Y坐标 */ adfGeoTransform[4] /* 旋转系数，通常为0 */ adfGeoTransform[5] /* 南北方向像元大小（高度，通常为负值，因为图像从上到下扫描） */ 因此，像元在东西方向的大小（宽度）就是adfGeoTransform[1]的绝对值，而南北方向的大小（高度）是adfGeoTransform[5]的绝对值（因为高度通常为负值，表示Y坐标递减）。 步骤： 1. 使用GDAL打开TIFF文件。 2. 获取图像的地理变换（GeoTransform）。 3. 从地理变换数组中提取像元宽度和高度。 注意：如果地理变换不可用（例如返回的是[0,1,0,0,0,1]），则可能意味着该图像没有地理参考信息，此时像元大小可能为1（单位未知），但通常我们处理的是有地理参考的影像。 下面分别用Python和Java代码示例如何获取像元大小。 Python示例代码： ```python from osgeo import gdal # 打开TIFF文件 dataset = gdal.Open('path/to/your/image.tif', gdal.GA_ReadOnly) if dataset is None: print("打开文件失败") else: # 获取地理变换 geotransform = dataset.GetGeoTransform() if geotransform is None: print("未找到地理变换信息") else: # 像元宽度（东西方向）是geotransform[1]的绝对值 pixel_width = abs(geotransform[1]) # 像元高度（南北方向）是geotransform[5]的绝对值（因为通常是负值） pixel_height = abs(geotransform[5]) print("像元宽度: {0} 单位".format(pixel_width)) print("像元高度: {0} 单位".format(pixel_height)) dataset = None # 关闭数据集 ``` Java示例代码（参考GDAL Java绑定）： ```java import org.gdal.gdal.Dataset; import org.gdal.gdal.gdal; import org.gdal.gdalconst.gdalconstConstants; public class GetPixelSize { public static void main(String[] args) { gdal.AllRegister(); Dataset dataset = gdal.Open("path/to/your/image.tif", gdalconstConstants.GA_ReadOnly); if (dataset == null) { System.err.println("打开文件失败: " + gdal.GetLastErrorMsg()); } else { double[] geotransform = dataset.GetGeoTransform(); // 检查地理变换是否有效 if (geotransform == null || geotransform.length < 6) { System.err.println("未找到地理变换信息"); } else { double pixelWidth = Math.abs(geotransform[1]); double pixelHeight = Math.abs(geotransform[5]); System.out.println("像元宽度: " + pixelWidth + " 单位"); System.out.println("像元高度: " + pixelHeight + " 单位"); } dataset.delete(); // 释放资源 } } } ``` 注意：单位通常是地图单位（如米），具体取决于坐标参考系统（CRS）。如果图像是地理坐标（经纬度），则像元大小将是度。 另外，如果图像有旋转（即geotransform[2]和geotransform[4]不为0），那么像元大小就不是简单的矩形，但这种情况较少见，通常旋转系数为0。在旋转的情况下，像元大小仍然可以通过计算每个像素的实际地理范围得到，但比较复杂。不过，大多数地理参考图像都是北朝上的，所以旋转系数为0。 因此，在无旋转的情况下，我们直接取绝对值即可。 引用说明： 根据引用[3]中提到的函数`ImageRowCol2Projection`，它使用了地理变换来将行列号转换为投影坐标，这也印证了地理变换的重要性。 如果用户需要更精确的计算（例如考虑旋转），则需要使用更复杂的方法，但一般情况下，上述方法即可满足需求。 相关问题： 1. 如何获取TIFF文件的地理坐标信息？ 2. 如何使用GDAL读取TIFF文件中的某个特定位置的像元值？ 3. 如何判断TIFF文件是否具有地理参考信息？ 4. 当TIFF文件很大时（比如超过4GB），GDAL能否正确读取？如何解决大文件读取问题？（参考引用[2]） 5. 如何使用GDAL进行TIFF文件的裁剪和分块操作？（参考引用[1]）</think>### GDAL 获取 TIFF 文件像元大小的方法 在 GDAL 中获取 TIFF 文件的像元大小（即像素分辨率）主要通过地理变换矩阵（GeoTransform）实现。以下是具体方法和步骤： #### 核心方法 1. **获取地理变换矩阵** 使用 `GetGeoTransform()` 方法获取包含 6 个元素的数组： ```python from osgeo import gdal dataset = gdal.Open("image.tif") gt = dataset.GetGeoTransform() # 返回 [x_min, x_res, rot_x, y_max, rot_y, y_res] ``` - `gt[1]`：东西方向像元宽度（X 分辨率） - `gt[5]`：南北方向像元高度（Y 分辨率，通常为负值） 2. **计算实际像元大小** ```python pixel_width = abs(gt[1]) # 像元宽度（绝对值） pixel_height = abs(gt[5]) # 像元高度（绝对值） ``` #### 完整代码示例（Python） ```python from osgeo import gdal def get_pixel_size(tif_path): dataset = gdal.Open(tif_path) if dataset is None: raise ValueError("无法打开文件") # 获取地理变换矩阵 gt = dataset.GetGeoTransform() if gt is None: raise RuntimeError("无地理参考信息") # 计算像元大小 pixel_width = abs(gt[1]) pixel_height = abs(gt[5]) dataset = None # 关闭数据集 return (pixel_width, pixel_height) # 使用示例 try: width, height = get_pixel_size("input.tif") print(f"东西方向像元大小: {width} 单位") print(f"南北方向像元大小: {height} 单位") except Exception as e: print(f"错误: {str(e)}") ``` #### 关键说明 1. **单位含义** 返回值的单位取决于 TIFF 文件的坐标系： - 地理坐标系（经纬度）：单位为度 - 投影坐标系（如 UTM）：单位为米 - 可通过 `GetProjection()` 获取坐标系详情 2. **特殊情况处理** - 若图像有旋转（`gt[2]` 或 `gt[4]` 非零），需通过矩阵计算实际像元大小： $$ \text{实际宽度} = \sqrt{gt[1]^2 + gt[2]^2} $$ $$ \text{实际高度} = \sqrt{gt[4]^2 + gt[5]^2} $$ - 无地理参考时，`GetGeoTransform()` 返回 `None` 3. **Java 实现参考** ```java import org.gdal.gdal.Dataset; import org.gdal.gdal.gdal; double[] gt = dataset.GetGeoTransform(); double pixelWidth = Math.abs(gt[1]); double pixelHeight = Math.abs(gt[5]); ``` #### 注意事项 - 对于大型 TIFF 文件（>4GB），需确保 GDAL 版本 ≥ 2.0 且 libtiff ≥ 4.0[^2] - 旋转图像的像元大小需通过矩阵计算（见上文公式） - 无地理参考时，默认像元大小为 `(1, 1)`（无实际地理意义） --- ### 相关问题 1. 如何判断 TIFF 文件是否包含有效的地理参考信息？ 2. 当 TIFF 文件存在旋转时，如何准确计算像元的实际物理尺寸？ 3. GDAL 在处理超过 4GB 的大型 TIFF 文件时有哪些注意事项？ 4. 如何通过 GDAL 获取 TIFF 文件中特定地理坐标的像素值？ 5. 除了像元大小，GDAL 还能提取哪些重要的地理空间元数据？
评论 5
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值