GDAL相关小技巧

1)  在调用GDAL相关功能时,设置参数  CPL_DEBUG为 ON可以输出各种调试信息,示例如下:

       ogr2ogr -skipfailure -overwrite -f "GPKG" "1.gpkg" D:\1.gdb --config CPL_DEBUG ON

2) 调试GDAL的源代码中的方法

       ① 编译GDAL的debug版本;

       ② 将编译生成的调试文件,

         1)gdal18.dll

         2)gdal_i.exp

         3)gdal_i.lib

         4)gdal.lib

         5)gdal18.pdb

         6)gdal18.ilk

         7)gdal18.exp等文件

         ③ 将gdal18开头的文件拷贝到自己工程的生成目录中;

         ④ 调试自己的程序

         在执行到GDALOpen函数(或者其他GDAL的函数)时按F11键,就会进入到GDAL的源代码中进行调试GDAL代码。总结:进行代码调试需要具备下述条件。1. 具有下述文件(*.pdb ,*.lib,*.ilk,*.exp);2 这些文件被放到了中间目录下。

3) GDAL通过函数接口设置环境变量的接口:GDALSetConfigOption( "GDAL_DATA", "d:\\gdal\\data" );

4) 使用GDAL插件形式编译一些附加库时,如果目录中包含中文则不能正确连接到插件的dll文件

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
Python的GDAL库提供了裁剪栅格影像的功能,并且可以根据需要通过最小外接矩形来裁剪影像。下面是用Python GDAL进行最小外接矩形裁剪的步骤: 首先,导入需要的库,包括GDAL库和OGR库。GDAL库用于处理栅格数据,OGR库用于处理矢量数据。 ```python from osgeo import gdal, ogr ``` 然后,打开需要裁剪的栅格影像,并获取其投影信息和地理转换信息。 ```python # 打开栅格影像 ds = gdal.Open("input.tif") # 获取投影信息和地理转换信息 proj = ds.GetProjection() gt = ds.GetGeoTransform() ``` 接下来,定义最小外接矩形的位置和大小。可以根据需要手动指定,也可以通过矢量数据计算得到。 ```python # 手动指定最小外接矩形位置和大小 xmin, ymin, xmax, ymax = 10, 20, 30, 40 ``` 然后,创建最小外接矩形的矢量数据。 ```python # 创建最小外接矩形的矢量数据 ring = ogr.Geometry(ogr.wkbLinearRing) ring.AddPoint(xmin, ymin) ring.AddPoint(xmin, ymax) ring.AddPoint(xmax, ymax) ring.AddPoint(xmax, ymin) ring.AddPoint(xmin, ymin) poly = ogr.Geometry(ogr.wkbPolygon) poly.AddGeometry(ring) ``` 接下来,设置裁剪的参数,并执行裁剪操作。 ```python # 设置裁剪参数 options = ["ALL_TOUCHED=TRUE"] # 执行裁剪操作 gdal.Warp("output.tif", ds, format="GTiff", cutlineDSName=poly, dstSRS=proj, outputType=gdal.GDT_Float32, options=options) ``` 最后,关闭数据集并输出裁剪后的影像。 ```python # 关闭数据集 ds = None # 输出裁剪后的影像 output = gdal.Open("output.tif") ``` 以上是使用Python GDAL进行最小外接矩形裁剪的步骤。通过指定最小外接矩形的位置和大小,可以裁剪出指定区域的栅格影像。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值