分类：gdal

2020
09-24

在Python中用GDAL实现矢量对栅格的切割实例

概述：本文讲述如何在Python中用GDAL实现根据输入矢量边界对栅格数据的裁剪。效果：裁剪前矢量边界裁剪后实现代码：#-*-coding:utf-8-*-"""@authorlzugis@date2017-06-02@brief利用shp裁剪影像"""fromosgeoimportgdal,gdalnumeric,ogrfromPILimportImage,ImageDrawimportosimportoperatorgdal.UseExceptions()#Thisfunctionwillconverttherasterizedclippershapefile#toamaskforusewithinGDAL... 继续阅读 >

2020
09-24

GDAL 矢量属性数据修改方式(python)

Case：需要给一个现有的shp数据创建一个字段，并将属性表中原有的一个文本类型的属性转换为整型后填入新创建的字段。Problem：新字段创建成功，但是赋值操作无效，即无法成功给字段写入值。solution：对字段进行赋值后需要，重新写入Feature，否则赋值无效，即layer0.SetFeature(feature)。特别注意：在对数据进行读写操作，一定要以读写的方式打开，即Open(filePath,1),该方法的原型为Open(pszName,intbUpdate=false),并且具有... 继续阅读 >

2020
09-24

python使用gdal对shp读取,新建和更新的实例

昨天要处理一个shp文件，读取里面的信息，做个计算然后写到后面新建的field里面。先写个外面网上都能找到的新建和读取吧。1.读取shp文件#-*-coding:cp936-*-try:fromosgeoimportgdalfromosgeoimportogrexceptImportError:importgdalimportogrdefReadVectorFile():#为了支持中文路径，请添加下面这句代码gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8","NO")#为了使属性表字段支持中文，请添加... 继续阅读 >

2020
09-24

python+gdal+遥感图像拼接(mosaic)的实例

作为摄影测量与遥感的从业者，笔者最近开始深入研究gdal，为工作打基础！个人觉得gdal也是没有什么技术含量，调用别人的api。但是想想这也是算法应用的一个技能，多学无害！关于遥感图像的镶嵌，主要分为6大步骤：step1：1）对于每一幅图像，计算其行与列；2）获取左上角X,Y3）获取像素宽和像素高4）计算maxX和minY，切记像素高是负值maxX1=minX1+(cols1*pixelWidth)minY1=maxY1+(rows1*pixelHeight）step2：计... 继续阅读 >