文章目录
- 1. ArcCatalog创建shp文件
- 2. ArcMap进行点线面的勾画
- 3. 矢量栅格化
1. ArcCatalog创建shp文件
- 先在目标文件夹中利用arccatalog创建矢量文件,如果是分割任务,就新建面矢量。
- 选择原始影像对应的坐标系
2. ArcMap进行点线面的勾画
- 将步骤一中生成的矢量文件利用arcmap打开,进行点线面的勾描
- 遇到需要在面状矢量打洞的情况
参考javascript:void(0)
3. 矢量栅格化
def vector2raster():
inputfilePath = '0000000001_V1/0000000001_V1_LINE.shp'
outputfile = 'raster_clip.tif'
templatefile = '0000000001_V1/0000000001.tif'
data = gdal.Open(templatefile, gdalconst.GA_ReadOnly)
x_res = data.RasterXSize
y_res = data.RasterYSize
vector = ogr.Open(inputfilePath)
layer = vector.GetLayer()
targetDataset = gdal.GetDriverByName('GTiff').Create(outputfile, x_res, y_res, 3, gdal.GDT_Byte)
targetDataset.SetGeoTransform(data.GetGeoTransform())
targetDataset.SetProjection(data.GetProjection())
band = targetDataset.GetRasterBand(1)
NoData_value = -999
band.SetNoDataValue(NoData_value)
band.FlushCache()
gdal.RasterizeLayer(targetDataset, [1, 2, 3], layer, )