GDAL C++ API 学习之路 （2）Dataset篇 代码示例 翻译 自学

已于 2023-04-19 14:18:49 修改
阅读量533 收藏 1
点赞数
分类专栏： GDAL 文章标签： 学习 c++
于 2023-04-01 22:33:31 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_69574549/article/details/129900938
版权

开头

紧接着第2章未完成的，继续GDALDataset 的API学习

GetGeoTransform

virtual CPLErr GetGeoTransform（double *padfTransform)         此方法与 C GDALGetGeoTransform（） 函数执行相同的操作

参数:  padfTransform – 一个现有的六个双缓冲区，转换将被放置到其中。

返回: CE_None成功时，或者CE_Failure无法获取转换             默认变换为 （0，1，0，0，0，1），即使返回CE_Failure错误（例如不支持转换为投影坐标的格式），也应返回默认变换。

double adfGeoTransform[6];         //定义了一个大小为 6 的双精度浮点型数组
dataset->GetGeoTransform(adfGeoTransform);       //函数将会将获取到的地理变换参数值存储到这个数组中

无论使用一维数组还是二维数组，都需要保证数组大小至少为6，以存储6个变换系数

地理变换参数的含义如下

  1. adfGeoTransform[0]：左上角的 X 坐标，即图像左上角像素的地理坐标的 X 坐标值。
  2. adfGeoTransform[1]：图像的水平分辨率，即每个像素在 X 方向上的实际宽度。
  3. adfGeoTransform[2]：X 方向上的旋转角度，通常为 0，即图像平行于 X 轴。
  4. adfGeoTransform[3]：左上角的 Y 坐标，即图像左上角像素的地理坐标的 Y 坐标值。
  5. adfGeoTransform[4]：Y 方向上的旋转角度，通常为 0，即图像平行于 Y 轴。
  6. adfGeoTransform[5]：图像的垂直分辨率，即每个像素在 Y 方向上的实际高度，通常为负数，表示 Y 轴坐标值的减小方向是从上到下。

其中，前两个参数（左上角的 X 坐标和图像的水平分辨率）用于计算像素的地理坐标的 X 坐标值，后两个参数（左上角的 Y 坐标和图像的垂直分辨率）用于计算像素的地理坐标的 Y 坐标值，中间两个参数（X 方向上的旋转角度和 Y 方向上的旋转角度）通常为 0，表示图像没有旋转。 

具体操作：http://t.csdn.cn/ADkbN

SetGeoTransform

virtual CPLErr SetGeoTransform（double *padfTransform)             此方法与 C GDALSetGeoTransform（） 函数执行相同的操作

设置仿射变换系数

参数:padfTransform – 一个六双精度缓冲区，其中包含要与数据集一起写入的转换系数。

返回:CE_None成功，如果无法写入此转换，则CE_Failure。

double adfGeoTransform[2][3] = {{pixelWidth, 0, xOrigin}, {0, pixelHeight, yOrigin}};
poDataset->SetGeoTransform(adfGeoTransform[0]);

AddBand

virtual CPLErr AddBand（GDALDataType eType， char **papszOptions = nullptr)

参数:eType – 新波段中像素的数据类型。

        papszOptions – NAME=VALUE 选项字符串的列表。支持的选项特定于格式。默认情况下可能传递 NULL。

返回:CE_None成功或CE_Failure失败

在成功完成后并不会直接返回新的GDALRasterBand对象指针，而是需要使用GDALDataset：：GetRasterBand（GDALDataset：：GetRasterCount（））方法获取新添加的波段。因为新添加的波段始终是数据集的最后一个波段，而GetRasterCount（）方法可以返回当前数据集的波段数量。所以，通过这种方法可以获取到新添加的波段的指针，然后再进行后续的读写操作等

 GDALDataset* poDS = static_cast<GDALDataset*>(GDALOpen("example.tif", GA_Update));
 poDS->AddBand(GDT_Byte, NULL);
 GDALRasterBand* poBand = poDS->GetRasterBand(poDS->GetRasterCount() + 1);

papszOptions

  1. “PIXELTYPE={BYTE/UINT16/INT16/UINT32/INT32/FLOAT32/FLOAT64/CPL_STRING}”：设置像素类型。默认为 BYTE。

  2. “INTERLEAVE={ BAND / PIXEL }”：设置波段排列方式。默认为 BAND，表示按波段排列。

  3. “COMPRESS={ NONE / LZW / JPEG / DEFLATE / ZSTD / PACKBITS }”：设置压缩方式。默认为 NONE。

  4. “PHOTOMETRIC={ MINISBLACK / MINISWHITE / RGB / PALETTE / CIELAB / ICCLAB / ITULAB / YCBCR }”：设置波段颜色编码方式。默认为 MINISBLACK。

  5. “NBITS=n”：设置像素的位数。n 可以为 1 到 64。默认为 8

  6. “BLOCKXSIZE=n” 和 “BLOCKYSIZE=n”：设置块的大小。默认为 256。

  7. “TILED=YES”：启用瓦片模式。默认为 BAND 模式。

  8. “BIGTIFF=YES”：启用 BigTIFF 格式。默认为普通 TIFF 格式。

GetInternalHandle

virtual void *GetInternalHandle（const char *pszHandleName)            此方法与 C GDALGetInternalHandle（） 方法相同

提取特定于格式的内部有意义的句柄

参数: pszHandleName – 所需的句柄名称。有意义的名称将特定于文件格式。

返回: 所需的句柄值，如果无法识别/支持，则为 NULL

Handle：

  • “METADATA”： 光栅波段的元数据句柄，用于存储和访问数据集的描述性信息，例如名称、作者、时间戳、坐标参考系统、地理变换信息等。

  • “IMAGE_STRUCTURE”： 光栅波段的影像结构句柄，用于存储和访问影像的附加结构信息，例如压缩类型、块大小、缩略图等。

  • “RPC”： 光栅波段的RPC（Rational Polynomial Coefficients）模型句柄，用于存储和访问RPC模型信息，该模型可用于对图像进行几何校正和重投影

  • “GCP”： 光栅波段的GCP（Ground Control Points）模型句柄，用于存储和访问GCP模型信息，该模型可用于对图像进行几何校正和重投影。

  • “DIMENSION”： 光栅波段的尺寸信息句柄，用于存储和访问波段的像素宽度、高度和深度等信息。

  • “NBITS”： 光栅波段的位深度句柄，用于存储和访问波段的位深度信息，例如8位、16位、32位等。

  • “SCALE”： 光栅波段的比例尺句柄，用于存储和访问波段的比例尺信息，例如像素宽度和高度的比例尺。

  • “OFFSET”： 光栅波段的偏移量句柄，用于存储和访问波段的偏移量信息，例如像素坐标的偏移量。

GetDriver

virtual GDALDriver *GetDriver(void)       此方法与 C GDALGetDatasetDriver（） 函数相同

提取与此数据集相关的驱动程序

返回：GDALDriver类的指针

const char* filepath = "example.tif";

// 打开数据集
dataset = (GDALDataset*) GDALOpen(filepath, GA_ReadOnly);

// 获取数据集所使用的驱动程序的指针
driver = dataset->GetDriver();

// 输出驱动程序的名称
printf("Driver name: %s\n", driver->GetDescription());

GetFileList

virtual char **GetFileList(void)            此方法与 C GDALGetFileList（） 函数相同

获取形成数据集的文件  

返回:NULL 或以 NULL 结尾的文件名数组

获取当前数据集所包含的文件列表，是指获取一个数据集所有文件的路径列表，这些文件通常包含该数据集的各个组成部分或属性。在地理空间数据处理中，数据集通常不是一个单独的文件，而是由多个文件组成的，每个文件存储数据集的不同部分或属性。例如，多波段遥感影像数据通常由多个单波段影像文件组成，矢量数据则可能由多个.shp、.dbf、.shx等文件组成。使用GetFileList函数可以获取数据集的所有文件路径，以便在需要时对其进行读取、写入或复制等操作

 GDALDataset* poDataset = (GDALDataset*)GDALOpen("sample.tif", GA_ReadOnly);
 char** papszFileList = poDataset->GetFileList();
 CSLDestroy(papszFileList);

GetDriverName

virtual const char *GetDriverName() 

返回驱动程序名称

返回:驱动程序名称

 GDALDataset* poDataset = (GDALDataset*)GDALOpen("sample.tif", GA_ReadOnly);
 // 获取驱动程序名称
 const char* pszDriverName = poDataset->GetDriverName();
 cout << "Dataset is using driver: " << pszDriverName << endl;

持续更新，下次继续GCP

场主不吃鍋巴
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
专栏目录
GDAL 自学入门使用文档
11-06
GDAL是一个用于栅格和矢量地理空间数据格式的转换程序库，在X/MIT风格的开放源代码下发布 许可 由 Open Source Geospatial Foundation . 作为一个库，它为所有支持的格式的调用应用程序提供单个栅格抽象数据模型和...
Python | 图像坐标、地理坐标和投影坐标的相互转化
Prince的博客
04-29 2878
from osgeo import gdal from osgeo import osr import numpy as np def getSRSPair(dataset): ''' 获得给定数据的投影参考系和地理参考系 :param dataset: GDAL地理数据 :return: 投影参考系和地理参考系 ''' prosrs = os...
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的GDALDataset::GetGeoTransform()函数
热门推荐
时日随心的专栏
08-03 2万+
GDALDataset中有函数GetGeoTransform()，函数声明为： CPLErr GDALDataset::GetGeoTransform  ( double *  padfTransform  )。         GDAL官网对该函数说明为：         Fetch the affine transformation coefficients.         Fet
【python函数】GetGeoTransform()
weixin_43711893的博客
01-16 1000
GeoTransform[0]，左上角横坐标（应该是投影坐标） GeoTransform[2]，行旋转 GeoTransform[1]，像元宽度(影像在水平空间的分辨率) GeoTransform[3]，左上角纵坐标（应该是投影坐标） GeoTransform[4]，列旋转 GeoTransform[5]，像元高度(影像在垂直空间的分辨率)， 如果影像是指北的，GeoTransform[2]和GeoTransform[4]这两个参数的值为0，GeoTransform[5]为负； 如果图像不含地
gdal 图层内相交要素合并
qq_42022928的博客
04-22 1487
一、整体思路 将图层内所有要素合并要一个集合类型中，集合类型进行内部合并，遍历集合类型找原图层内查询属性信息，进行属性信息补充。 二、代码实现 # pgurl pg数据库连接地址 # lyrname 合并图层名 from osgeo import gdal, ogr, osr pg_driver = ogr.GetDriverByName("PostgreSQL") ds = pg_driver.Open(pgurl,gdal.GA_Update) sr = osr.SpatialReference
（源代码）Qt GDAL C++ 遥感图像显示 处理
07-26
曾经的程序源代码，运行环境：Qt5.12，GDAL2.3.1，纯C++编写
GDAL-API手册
04-01
GDAL开发手册API详细说明，内含详细示例代码，适合初学者使用，非常棒的入门书！In pythe Import go from gdalconst import s datasct-gdal. Open( filename, GA ReadOnly if dataset is none 如果 GDALOpen(函数...
GDALtest_gdaldem_GDAL高程_C++_开源dem编辑_GDAL_
10-04
使用开源图像GDAL实现数字高程DEM文件的大规模快速读取，添加了数据格式匹配
Qt GDAL C++ 遥感图像显示 处理
05-05
软件包括了显示遥感图像，及鼠标拖动，随鼠标位置缩放，及部分栅格图像处理功能(图像锐化，平滑等)，和显示矢量图像的功能。
gdal--矢量求交
secyb的博客
05-08 6638
前言 最近在帮同学写几个小程序，包括矢量数据求交（Intersection）。网上相关文章很少，所以我参考GDAL官网写了一个完整的程序（裁剪、合并、擦除、更新等功能只要修改代码的求交函数部分，下面会有说明）。 代码 #include "gdal_priv.h" #include "ogrsf_frmts.h" //for ogr bool VectorIntersection(...
web项目上传不了文件Linux,C#实现Web文件上传的两种方法
weixin_29491259的博客
05-14 310
在Web编程中,我们常需要把一些本地文件上传到Web服务器上,上传后,用户可以通过浏览器方便地浏览这些文件，应用十分广泛。那么使用C#如何实现文件上传的功能呢?下面笔者简要介绍一下。首先,在你的Visual C# web project 中增加一个上传用的Web Form,为了要上传文件,需要在ToolBox中选择HTML类的File Field控件,将此控件加入到Web Form中,然而此时该控...
建立矢量缓冲区+合并矢量要素
XinemaChen的博客
03-03 389
程序包含把单个shapefile所有的要素都合并在一起，基于shapefile建立缓冲区两个功能
多个矢量文件合并为一个矢量文件
John的博客
08-10 2503
见如下代码。 这个代码是根据我自己业务需要写的，整体是可以根据自己的需要进行适当调整和修改。许多坑已经都跨过去了。 def mutil_shapefiles_merge(input_folder_path, output_shpfile_abspath, delete_area_threshold=0.03): '''本函数用于多个矢量文件的合并 第一个参数为需要合并矢量文件的文件夹路径 第二个是合并后矢量文件的完整路径''' ptdriver = ogr.GetDri
c++ GDAL 完成PCA
smallcubelo的博客
09-02 736
c++ GDAL 完成PCA int main() { const char* imagePath = "D:/data/before.img"; const char* txtPath = "D:/data/result.txt"; UpdateRasterFile(imagePath, txtPath); return 0; } 通过UpdateRasterFile(imagePath, txtPath);完成PCA的计算和输出 void UpdateRasterFile(const cha
c++ gdal 矢量转栅格_python gdal 矢量转栅格
weixin_28761983的博客
01-26 136
data = gdal.Open(templateTifFileName, gdalconst.GA_ReadOnly)geo_transform = data.GetGeoTransform()x_min = geo_transform[0]y_min = geo_transform[3]x_res = data.RasterXSizey_res = data.RasterYSizemb_v =...
gdal栅格数据处理(python)
陨星落云的博客
06-30 2016
栅格数据处理 关闭栅格数据集 这里显示如何关闭栅格数据集。在脚本的中间非常有用，它可以恢复访问数据集所拥有的资源，删除文件等。在脚本的末尾没有必要，因为Python垃圾收集器会在脚本运行时自动执行相同的操作。脚本退出。 from osgeo import gdal # 打开数据集 ds = gdal.Open("jianda.tif") # 关闭数据集 ds = None 获取栅格数据集的元数据 获取栅格元数据以进行快速的分辨率检查 from osgeo import gdal gtif = gda
python 读取TIFF信息
weixin_43711893的博客
01-16 5251
用python读取TIFF的行列数、投影转换信息、地图投影信息以及栅格像元值。
【C++】【GADL】读取栅格数据获取信息
胡侃有料的博客
08-11 1271
every blog every motto: Light tomorrow with today. 0. 前言 简单记录下用c++读取栅格数据 1. 正文 #include<iostream> #include "gdal_priv.h" using namespace std; // 读取栅格数据 void loadRasterData() { const char * fileName = "D:\\Programme\\Viual Studio\\C++\\learning\\
gdal加载影像数据的代码示例
最新发布
05-24
以下是使用GDAL加载影像数据的Python代码示例： ```python import gdal # 打开影像文件 ds = gdal.Open('path/to/image.tif') # 获取影像的宽度和高度 width = ds.RasterXSize height = ds.RasterYSize # 获取...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

场主不吃鍋巴

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值