python 栅格计算,python在栅格中运用函数
GDAL Python实现栅格影像处理之栅格矢量化及矢量栅格化相关概念使用方法效果展示
相关概念
栅格矢量化：将栅格数据生成矢量数据。
矢量栅格化：将矢量数据生成栅格数据。
使用方法栅格矢量化
使用gdal.Polygonize()方法参数说明（未列完)参数说明srcBand输入数据波段掩蔽带设置掩码波段（只对掩码区域进行处理)outLayer矢量化后的矢量图层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.RasterizeLayer()方法参数说明（未列完)参数说明资料组输出的栅格数据，注意该数据必须以更新模式打开法官指定要更新的波段个数（更新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的观点，版权归原作者所有，仅为传播更多信息之目的，如有侵权请联系，我们将第一时间修改或删除，多谢。