主要用的函数
- 拼接用到的函数是镶嵌至新栅格的函数arcpy.MosaicToNewRaster_management()
官方api:https://resources.arcgis.com/zh-cn/help/main/10.1/index.html#/na/001700000098000000/ - 裁剪用到的函数是按掩膜提取arcpy.gp.ExtractByMask_sa()
官方api:https://resources.arcgis.com/zh-cn/help/main/10.1/index.html#/na/009z0000002n000000/
代码
import arcpy
import os
arcpy.env.overwriteOutput = True
arcpy.env.workspace = 'D:\\testData\\ganziTest'
mosaic_tif = "D:\\testData\\ganzipytest\\mosaic.tif"
output_location = "D:\\testData\\ganzipytest"
testclip_shp = "D:\\testData\\testclip.shp"
clip_mosaict_tif = "D:\\testData\\ganzipytest\\clip_mosaict.tif"
try:
coordinate_system_for_the_raster = ""
rasters = arcpy.ListRasters("*", 'tif')
tifName = "mosaic"
strList = ''
for i in range(0, len(rasters)):
strList += rasters[i]
strList += ";"
arcpy.MosaicToNewRaster_management(strList + tifName + '.tif', output_location, "mosaic.tif",
coordinate_system_for_the_raster, "16_BIT_UNSIGNED", "", "1", "LAST", "FIRST")
print "finished"
arcpy.gp.ExtractByMask_sa(mosaic_tif, testclip_shp, clip_mosaict_tif)
except Exception as e:
print(e)