命令行记录-矢量向栅格转换
1、 (1)重点学习如何应用 gdal.RasterizeLayer 函数 gdal.RasterizeLayer( target_ds, [3, 2, 1], source_lyr, burn_values = [10,10, 55], options = ["ALL_TOUCHED=TRUE"])注意:栅格化时,将波段 1,2,3 上,分别赋值 55,10,10,source_lyr为打开的文件,burn_values和options这两个属性是非必须的。 (2)关键是选项的选择问题,有4种赋值方式,如: options=["ATTRIBUTE=LUCODE"],表示字段" LUCODE "的属性值将为栅格值,如果没有设定字段名,则赋值为0; options=["BURN_VALUE_FROM=Z"],表示 3D 图形的栅格值为其 Z 值(高程值) ; options=["ALL_TOUCHED=TRUE"],表示图形接触到的像素均将输出; options=["MERGE_ALG=ADD","MERGE_ALG=REPLACE"], 表示多图形覆盖同一像素的取值方式 2、 from osgeo import gdal, ogr, osr #定义投影 sr = osr.SpatialReference('LOCAL_CS["arbitrary"]') #在内存中