GDAL
文章平均质量分 60
Hexk_
这个作者很懒,什么都没留下…
展开
-
GDAL Python 过滤Shape Polygon中的面积小于某个阈值的小图斑
【代码】GDAL Python 过滤Shape Polygon中的面积小于某个阈值的小图斑。原创 2023-09-06 11:10:31 · 246 阅读 · 0 评论 -
GDAL循环读取矢量文件中的geometry
GetArea和Area在geom尺度上是同一个函数。按照feature数量进行循环。3就代表polygon。原创 2023-09-06 10:12:32 · 228 阅读 · 0 评论 -
关于gdal numpy计算NDVI等指数的问题
在利用python gdal numpy计算NDVI等运算栅格数组的指数时出现nan和inf的问题。会出现nan和inf,在QGIS中加载出现问题。去除nan和inf,并进行二值化。以上得到np多维数组。原创 2023-09-05 10:24:33 · 182 阅读 · 0 评论 -
12_OSR模块和空间参考系
定义在WGS84基准面下的北半球UTM17带。设置一个定义用户名字的投影坐标系统并确定系统被投影过,应该是参数是用户随意输入,只表示用于定义坐标系的名称;设定一个坐标系统,EPSG中已有的各种参数;设置投影转换参数细节,包含两个参数,zone表示投影带号,north表示是否为北半球,默认是True。原创 2022-11-30 15:41:17 · 1186 阅读 · 0 评论 -
11_空间计算
包含一些空间分析的方法和地理数据处理函数。原创 2022-11-30 15:40:44 · 464 阅读 · 0 评论 -
10_空间过滤器Spatial filters
使用来实现根据空间位置的筛选。先设置两个shp文件,其中的要素需要相交。首先需要获取重叠两个shp文件的图层。获得叠加图层的空间坐标使用函数对原有图层进行筛选,此时的是已经被筛选过后的。直接CopyLayer这个图层形成新的shp文件。这个函数可以设定周围四个边界点形成矩形框。效果和shp筛选一致。注意:经过后的layer并不是完整的。若需要对layer执行其他操作,那么需要将其进行筛选清空。原创 2022-11-30 15:38:11 · 354 阅读 · 0 评论 -
9_OGR创建Shapefile
使用ogr创建一个shp文件,并且生成一个点坐标,添加两个字段并填入内容。整个逻辑过程如下:首先,生成SHP的Driver;判断当前路径是否已经有同名的SHP文件存在,若存在则通过os模块删除;通过Driver创建数据集,即创建矢量文件;为矢量文件新建图层,同时在参数中设定图层名称、空间参考系、几何类型;设定字段,并在矢量文件的图层中添加字段;为图层创建feature;创建几何图形,并设定几何图形的坐标,将其添加到feature中;为feature设定字段属性;将feature写入。原创 2022-11-30 15:36:57 · 1036 阅读 · 0 评论 -
8_读取矢量数据
读取矢量数据的方式就是上一章所讲的两种方法。原创 2022-11-09 14:56:25 · 470 阅读 · 0 评论 -
7_使用OGR模块操作矢量数据
其实参数都是一样的,第一个参数表示传入的路径,第二个表示读取方式。在实际的编程过程中,最好是先声明打开文件类型的。文件而言,这个是缺省的打开方式。函数之间有点小小的不同。一样的操作,但是对于。原创 2022-11-09 14:56:00 · 241 阅读 · 0 评论 -
6_GDAL和Pillow互操作
GDAL和PIL很相像,它们处理和操作的对象都是栅格图像。区别在于GDAL主要重点放在地理或遥感数据的读写和数据建模以及地理定位和转换;PIL的重点是放在图像本身处理上的。在底层逻辑上,两个库都是对numpy的数据进行处理,因此理论上可以相互共享和交换数据。GDAL的核心在波段band, 一切操作的基础和核心都在波段。波段可以单独拿出来操作,至于波段在数据集中的顺序无关紧要。遥感图像大多比RGB图像的波段要多,而每个波段单独都是一个完整的整体, 每个波段单独拿出来都是一个数据集。Pillow。原创 2022-11-06 20:54:30 · 255 阅读 · 0 评论 -
5_GDAL创建影像
该函数用于创建影像。指保存的路径和文件名称;指生成影像的长和宽;int bands指生成影像的波段数目,默认是1;指影像的DataType,一般与原始影像的DT相同。生成影像的基本流程是:先声明一个指定图像格式的Driver;通过函数生成初始化的影像,同时指明路径、大小、波段和DT;为Create后的影像设置Projection、;写入具有值的矩阵到新影像中。gdal。原创 2022-11-06 20:53:28 · 4077 阅读 · 0 评论 -
4_获取影像像元数据
本文中的方法可以对第三章进行补充,能够访问像元灰度值,并对其进行处理。GDAL提供了下面两个函数来访问影像的数值。读取图像数据(以二进制的形式)读取图像数据(以数组的形式)xoff, yoff指定想要读取的部分原点位置在整张图像中距离全图原点的位置;指定要读取部分图像的矩形的长和宽;~~放缩比例,最后读取影像的长宽会根据这个比例因子进行缩放;~~测试发现这两个参数并不是缩放比例因子,并不是缩放比例因子,而是新的指定区域,比如指定为(20, 20),那么。原创 2022-11-05 21:05:43 · 866 阅读 · 0 评论 -
3_读取遥感影像信息
本文中的所有方法可以读取遥感影像信息,但是仅限于概括地获取元数据信息,无法具体到对像元进行处理。原创 2022-11-05 21:04:25 · 2034 阅读 · 0 评论 -
2_开始使用GDAL
GDAL和OGR都归于osgeo除了GDAL包,还有一个gdalconst包最好也导入。gdalconst也是osgeo的一个包,它只是在代码中对GDAL中用到的一些常量进行了绑定。其中gdalconst中的常量都加了前缀,力图与其他模块冲突最小。原创 2022-11-05 21:03:47 · 641 阅读 · 0 评论 -
1_GDAL简介
包含GDAL的许多基本概念,感觉是机翻,但都是GIS和RS的部分基础衍生。GDAL——处理栅格空间数据的类库OGR——处理矢量数据的类库。原创 2022-11-05 21:02:11 · 432 阅读 · 0 评论