Hi ALL,
I want to mosaic a set of rasters in separate folders but I am stuck with my code. Basically, I have two sets of modis tiles located in two separate folders, and I want to create a mosaic dataset.
case:
Folder 1 = mod_h26v06_2012001.tif, mod_h26v06_2012009.tif, mod_h26v06_2012017.tif, mod_h26v06_2012025.tif,......., mod_h26v06_2012361.tif
Folder 2 = mod_h27v06_2012001.tif, mod_h27v06_2012009.tif, mod_h27v06_2012017.tif, mod_h27v06_2012025.tif,......., mod_h27v06_2012361.tif
What I want to do:
mod_h26v06_2012001.tif merge with mod_h27v06_2012001.tif, ........, mod_h26v06_2012361.tif merge with mod_h27v06_2012361.tif
This is my initial code:
# import necessary modules of ArcGIS
import arcpy, os, sys
from arcpy import env
from arcpy.sa import *
# this will overwrite output. Important in testing codes
arcpy.env.overwriteOutput = True
# activate the spatial analyst extention of ArcGIS
arcpy.CheckOutExtension("Spatial")
# specify your workspace. This is where raster is stored
raster1 = 'X:/raster_ws'
arcpy.env.workspace = raster1
h26v06 = [os.path.join(raster1, l) for l in arcpy.ListRasters("*2012*", "TIF")]
raster2 = 'X:/raster_ws'
arcpy.env.workspace = raster2
h27v06 = [os.path.join(raster2, r) for r in arcpy.ListRasters("*2012*", "TIF")]
out_ws = 'X:/output_ws/'
# (input rasters,outputloc,name w/ ext, coor, pixeltype, cellsize, #of bands
# mosaic method, mosaic color map)
for h in range(0,46):
parse = l[18:25]
filename = "modis_evi_" + parse + ".tif" #'modis_8day_2012_{:03d}.tif'.format(h+1)
mosaic = arcpy.Mosaic_management(Raster(h26v06);Raster(h27v06), out_ws, filename, "", "16_BIT_SIGNED", "", "1", "LAST","FIRST")
print mosaic
Any help is much appreciated.
Thanks,
-Leo