python栅格转矢量,
栅格矢量化:将栅格数据生成矢量数据。
矢量栅格化:将矢量数据生成栅格数据。
使用方法
栅格矢量化
使用gdal .多边形化方法
参数说明(未列完)
参数说明
srcBand输入数据波段
掩蔽带设置掩码波段(只对掩码区域进行处理)
户外者矢量化后的矢量图层
iPixValField需要将(同向下)向下的值写入矢量字段的索引
选择算法选项
GDALProgressFunc进度条回调函数
回收进度条参数
代码实现def PolygonizeTheRaster():输入file=r g:/study projects/gdal/gdal study/Files/images/2016 cha。TIF的ds=gdal .打开(输入文件,gdal .GA_ReadOnly) srcband=ds .GetRasterBand(1)mask band=src band .GetMaskBand()dst _ filename= G:/study projects/gdal/gdal study/Files/plogon/test。shp drv=ogr .GetDriverByName(ESRI形状文件)dst_ds=drv .创建数据源(dst _ filename)SRS=None dst _ layer name= out dst _ layer=dst _ dsCreateLayer(dst _ layer name,SRS=SRS)dst _ field name= DN FD=ogr .FieldDefn(dst_fieldname,ogr .整数的)dst _ layer .create field(FD)dst _ field=0 Prog _ func=test()options=[]#参数输入栅格图像波段\掩码图像波段、矢量化后的矢量图层、需要将(同向下)向下的值写入矢量字段的索引、算法选项、进度条回调函数、进度条参数gdal .Polygonize(srcband,maskband,dst_layer,dst_field,options,callback=prog_func)
矢量栅格化
使用gdal .光栅层()方法
参数说明(未列完)
参数说明
资料组输出的栅格数据,注意该数据必须以更新模式打开
法官指定要更新的波段个数(更新123波段)
层指定的图层
股骨近端髓内针变压器几何图形坐标转换图像行列号函数
pTransformArg几何图形坐标转换图像行列号参数
刻录值输出图像的像素值
选择控制栅格化的参数
GDALProgressFunc进度条回调函数
回收进度条参数
代码实现def RasterzeTheVectorToRasterr():print( lets gogogo )#参数说明:输出的栅格数据,注意该数据必须以更新模式打开、指定要更新的波段个数(更新123波段)、指定的图层、几何图形坐标转换图像行列号函数、几何图形坐标转换图像行列号参数、以及图层中属性字段属性值输入文件路径= G:/study projects/gdal/gdal study/Files/plogon/test。shp 输出文件= G:/study projects/gdal/gdal study/Files/images/vector。TIF模板文件= G:/study projects/gdal/gdal study/Files/images/2016 cha。TIF的数据=gdal .打开(templetefile,gdalconst .GA_ReadOnly) x_res=数据RasterXSize y_res=data .RasterYSize vector=ogr .打开(输入文件路径)层=矢量.GetLayer() targetDataSet=gdal .GetDriverByName(GTiff ).创建(输出文件,x_res,y_res,3,gdal .GDT _字节)#目标数据集=gdal .GetDriverByName(GTiff ).创建副本(模板文件,数据)目标数据集.SetGeoTransform(数据GetGeoTransform())目标数据集.集合投影(数据get projection())band=目标数据集。GetRasterBand(1)no data _ value=-999波段SetNoDataValue(NoData_value)波段FlushCache() gdal .栅格化层(目标数据集,[1,2,3],图层,选项=[属性=DN])
郑重声明:本文由网友发布,不代表盛行IT的观点,版权归原作者所有,仅为传播更多信息之目的,如有侵权请联系,我们将第一时间修改或删除,多谢。