GDALSetProjection使用的一个注意事项

最新推荐文章于 2022-06-10 19:31:51 发布
最新推荐文章于 2022-06-10 19:31:51 发布
阅读量1.6k 收藏 1
点赞数
文章标签： python java c/c++

GDALSetProjection 简述

GDALSetProjection是用来给GDALDataset设定投影信息(坐标系统)的接口，实际上是GDALDataset::SetProjection这个虚函数的转调而已。官网文档描述如下：

CPLErr GDALDataset::SetProjection (const char * pszProjection )
Set the projection reference string for this dataset.
The string should be in OGC WKT or PROJ.4 format. An error may occur because of incorrectly specified projection strings, because the dataset is not writable, or because the dataset does not support the indicated projection. Many formats do not support writing projections.
This method is the same as the C GDALSetProjection() function.
Parameters
pszProjection projection reference string.
Returns
CE_Failure if an error occurs, otherwise CE_None.

最低0.47元/天 解锁文章
weixin_34261739
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
python：GDAL库教程
养乐多的博客
04-14 3845
在本文将介绍如何使用Python中的GDAL库读取和保存遥感数据。对于读取遥感数据，我们通过gdal.Open()函数打开遥感数据集，然后使用GetRasterBand()函数获取数据集中的波段信息，并使用ReadAsArray()函数将数据读取为numpy数组。对于保存遥感数据，我们使用GDAL库中的gdal_array.SaveArray()函数，将numpy数组保存为遥感数据。
criteria.setProjection(Projectio…
随便起个不重复的名字
09-21 4501
criteria.setProjection(Projections.id()); 指criteria查询只返回表中id的数据， criteria.setProjection(Projections.id());和 criteria.setProjection(Projections.property("id"));效果一样！ 调用 setProjection()应用投影到一个查询，返回结果集合l
python+GDAL给图像设置地理参考（GeoTransform和Projection的设置）
KilllerQueen的博客
05-18 6050
使用GDAL给遥感影像设置地理参考和投影，分别设置两个参数：GeoTransform和Projection。对于没有附带地理信息的影像，需要从其它途径获取参数，如xml等元数据。 GDAL 中，GeoTransform是一个六个元素的元组。六个参数分别为： （左上角x坐标， 水平分辨率，旋转参数， 左上角y坐标，旋转参数，-垂直分辨率） 应该注意，最后一个参数是负值，旋转参数一般是0。而对于没有投影的图像（如WGS 84只有地理坐标系），它的分辨率不是以米位单位，而是经纬度，所以不能简单地指定分辨率这两个
Hibernate QBC查询
javacoffe的专栏
08-02 5446
 QBC查询：   QBC查询就是通过使用Hibernate提供的Query By Criteria API来查询对象，这种API封装了SQL语句的动态拼装，对查询提供了更加面向对象的功能接口。我们看下面的示例程序：Criteria criteria=session.createCriteria(User.class);criteria.add(Expression.eq(“nam
多个投影的设置
04-09 590
十五、投影、聚合和分组 投影运算实际上就是一个基于列的运算，通常用于投影到指定列（也就是过滤其他列，类似select子句的作用），还可以完成SQL语句中常用的分组、组筛选等功能。 Hibernate的条件过滤中使用Projection代表投影运算，Projection是一个接口，而Projections作为Projection的工厂，负责生成Projection对象。
python GDAL 影像处理系列之给遥感影像加地理信息（二）
向上生长的程序媛的博客
11-08 3139
在遥感影像中，地理坐标的信息异常重要，但是对于numpy array数据来说，没有地理位置信息的，因此如果需要地理信息的情况下，如果再将遥感影像作为普通多波段RGB影像来进行处理就有点不合适了，那怎么合适呢？我就来了。。。。（欢呼） 使用gdal的第一步，就需要安装gdal，安装教程在：gdal安装教程 为一张无地理信息的遥感影像添加地理信息，一个函数就可以呀！ from osg...
log4j 配置以及无法在控制台打印问题详解
TianXieZuoMaiKong的博客
03-15 5213
今天自己尝试着去配置了一下log4j，对于怎么配置，网上有一大把的详细说明，这里贴个链接也无妨（http://blog.csdn.net/dr_guo/article/details/50718063），所以，这篇博文，我就贴出我自己简单的配置文件代码，并讲解一下我遇到的问题。
【无标题】
hshxhdbdbhxh的博客
05-17 521
使用ArcGISPro2.5导出深度学习样本问题 最近需要在遥感影像上制作目标检测识别的VOC格式的深度学习样本数据集，在GF-2遥感影像上打好标签之后，利用export training data for deep learning工具导出样本，参数如下： 导出之后发现能导出，但是导出结果打不开，提示文件损坏， 设置导出格式为png或者tif时候可以正常打开，但是打开后是全黑的，python像元值都为0，但是在envi里面打开是正常显示的， 在ENVI里面显示效果 后来我发现是因为tiff数据的数据类
GDAL开发：地理坐标与投影坐标的转换
Tjmies的博客
06-30 3195
最近运用到QGis的开发，通过gdal库（操作各种栅格地理数据格式的库）解析shp文件得到想要的地理数据，开发过程中，想得到某段道路的折线长度，网上找很多资料都没有类似的代码示例（主要是讲解地理坐标和投影坐标的定义，很清晰、易理解，但没有找到自己想要的示例代码）。 开发环境：visual studio2013 （windows）gdal-2.3.1 例：需要得到图中这段道路的长度（gdal开发） 如图所示，gdal代码需要获取道路的长度921.130米 OGRLineString *lin...
Hibernate之criteria.setProjection(null)作用
qq_30908729的博客
09-30 1304
Hibernate之c.setProjection(null)作用: criteria.setProjection(null); 目的是为了获得行数，并设置投影为空，为的是返回List出来，详情地址：http://yayihouse.com/yayishuwu/chapter/1004
GDAL 几何矫正
fb_help的专栏
03-28 1443
一、概述 遥感影像和地理坐标进行关联的方式一般有好几种，一种是直接给出了仿射变换系数，即6个参数，左上角地理坐标，纵横方向上的分辨率，以及旋转系数。在这种情况下，求出某一像素点的地理坐标非常容易，直接用公式可以求出，具体代码如下： void CPL_STDCALL GDALApplyGeoTransform(double *padfGeoTransform, double dfPixel, double dfLine,
gdal使用常见问题
小胖七少爷的博客
11-19 2635
1. 编译proj高版本时，会发现proj.lib库是静态库，如果不配置环境变量，会出现以下错误： 1）报错 2）gdal库中的OGRCreateCoordinateTransformation方法，创建对象失败 解决方法： 找到proj.db文件，需要设置用户变量。 2. 使用gdal时候，需要设置环境变量，必须指定到data文件夹所在路径； 3. 从gdal2.0升级到gdal3.0在使用上的变化; ...
Hibernate Projections（投影、统计、不重复结果）
路。。。。。
07-29 1万+
Hibernate除了处理查询结果集中的对象之外，还可以将结果集中的结果当做行和列集来使用，这与通过JDBC执行select查询获得的数据的使用方式相似。因此，Hibernate也支持属性、统计函数和Group By等查询。要想使用Hibernate的投影统计功能，首先要从org
C++下GDAL的详细使用案例（含项目配置、tif读取为cv::Mat、Mat保存为tif）
a486259的博客
06-10 1万+
直接到下列链接下载即可，按照说明，将lib目录添加的系统环境变量中即可windows下GDAL322的库-深度学习文档类资源-CSDN下载包含目录中设置include目录 库目录中设置lib的路径 附加依赖项中设置gdal_i.lib GDAL读取数钱需要注册一下驱动（用于编码解码图像的驱动），同时可以设置一下支持中文路径。加载数据时需要注意，GA_Update和GA_ReadOnly两种模式。 2.1 获取图像的尺寸 这里所有的波段的size都是一样的 2.2 获取图像的通道数 2.3
GDAL影像投影转换
热门推荐
心如止水-GISer的成长之路
12-21 3万+
GDAL影像投影转换 又有一段时间没有写博客了，今天就把影像（DEM/遥感影像等）投影转换的经验和大家分享下。那么关于GDAL影像投影转换的文章在网上可以说是有很多，如果你还不清楚地图投影那么我建议你去看一看地图投影和坐标系统的基础知识，在这里就不详细讲解了。 一、影像投影转换的概念影像投影转换就是将一个地理坐标系统转换到另一个坐标系统，如果在同一个椭球基准面下的转换就是严密的转换，如果在同一个椭
GDAL获取投影坐标系注意问题
wanguangtong的专栏
02-14 2786
GDAL提供了获取投影坐标系的C函数GDALGetProjectionRef以及对应的C++函数GetProjectionRef，但在获取投影坐标系之前需要设置通过CPLSetConfigOption函数设置”GDAL_DATA”环境变量，例如：CPLSetConfigOption(“GDAL_DATA”, “D:/gdal/data”);其中“D:/gdal/data”为GDAL源码下的data；
使用C++和GDAL新建影像（地理坐标系）
Lianj.Qin
12-05 2645
使用GDAL创建含投影信息的影像，通过SetGeoTransform设置影像的左上角经纬度位置以及分辨率、OGRSpatialReference设置投影信息。
setProjectionMap(sContactsProjectMap)
我相信分享也是一种能力
11-28 3738
sContactsProjectMap是一个hashmap, 看到他的key和value是一样的, 顿生疑惑,搜索一下,发现如下文字,解惑之.   SQLiteQueryBuilder中的setProjectionMap, 这里的ProjectionMap必须是一个hashmap,更重要的是： If a projection map is set it must conta
gdal拼接tif影像的C++代码
最新发布
05-05
以下是使用GDAL C API拼接多个tif影像的示例代码： ```c #include "gdal.h" int main() { GDALAllRegister(); // 注册所有的驱动 // 打开第一个tif影像 GDALDatasetH hDataset1 = GDALOpen("path/to/first/tif", GA_ReadOnly); if (hDataset1 == NULL) { printf("Open dataset failed.\n"); return 1; } // 获取第一个tif影像的基本信息 GDALRasterBandH hBand1 = GDALGetRasterBand(hDataset1, 1); int nXSize1 = GDALGetRasterXSize(hDataset1); int nYSize1 = GDALGetRasterYSize(hDataset1); GDALDataType eDataType1 = GDALGetRasterDataType(hBand1); double adfGeoTransform1[6]; GDALGetGeoTransform(hDataset1, adfGeoTransform1); // 创建输出影像 GDALDriverH hDriver = GDALGetDriverByName("GTiff"); GDALDatasetH hOutDataset = GDALCreate(hDriver, "path/to/output/tif", nXSize1, nYSize1, 1, eDataType1, NULL); GDALSetGeoTransform(hOutDataset, adfGeoTransform1); GDALSetProjection(hOutDataset, GDALGetProjectionRef(hDataset1)); // 将第一个tif影像写入输出影像 GDALRasterBandH hOutBand = GDALGetRasterBand(hOutDataset, 1); GDALRasterIO(hBand1, GF_Read, 0, 0, nXSize1, nYSize1, GDALGetRasterBand(hOutDataset, 1), nXSize1, nYSize1, eDataType1, 0, 0); // 关闭第一个tif影像 GDALClose(hDataset1); // 循环打开并拼接其余的tif影像 for (int i = 2; i <= nCount; i++) { // 打开第i个tif影像 GDALDatasetH hDataset = GDALOpen("path/to/ith/tif", GA_ReadOnly); if (hDataset == NULL) { printf("Open dataset failed.\n"); return 1; } // 获取第i个tif影像的基本信息 GDALRasterBandH hBand = GDALGetRasterBand(hDataset, 1); int nXSize = GDALGetRasterXSize(hDataset); int nYSize = GDALGetRasterYSize(hDataset); GDALDataType eDataType = GDALGetRasterDataType(hBand); double adfGeoTransform[6]; GDALGetGeoTransform(hDataset, adfGeoTransform); // 判断第i个tif影像是否和第一个tif影像兼容 if (nXSize != nXSize1 || nYSize != nYSize1 || eDataType != eDataType1 || adfGeoTransform[0] != adfGeoTransform1[0] || adfGeoTransform[1] != adfGeoTransform1[1] || adfGeoTransform[2] != adfGeoTransform1[2] || adfGeoTransform[3] != adfGeoTransform1[3] || adfGeoTransform[4] != adfGeoTransform1[4] || adfGeoTransform[5] != adfGeoTransform1[5]) { printf("Incompatible raster dimensions or parameters.\n"); return 1; } // 将第i个tif影像写入输出影像 GDALRasterIO(hBand, GF_Read, 0, 0, nXSize, nYSize, GDALGetRasterBand(hOutDataset, 1), nXSize, nYSize, eDataType, 0, 0); // 关闭第i个tif影像 GDALClose(hDataset); } // 关闭输出影像 GDALClose(hOutDataset); return 0; } ``` 该示例代码使用GDAL C API打开、读取、创建、写入和关闭tif影像，拼接多个tif影像并生成输出影像。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值