merge_DOM_ArcGIS.py
cf8c25f31e2eaw_solar

merge_DOM_ArcGIS.py
View Options

			import arcpy
			import os

			target_path = r'D:\Surface_data\DOM\DOM'
			path = r'D:\Surface_data\DOM\DOM master cells'

			# append all pixels to a list
			list_of_pixels = []
			for root, dirs, files in os.walk(path):
			for file in files:
			if file.endswith(".asc"): list_of_pixels.append(file)

			# create individual partial files for the merged DOM
			n_parts = 20

			target_files = [("DOM_merged_%02d" %n) for n in range(n_parts)]
			list_of_files = ["" for n in range(n_parts)]

			divisor = int( np.ceil( len(list_of_pixels)/n_parts ))
			for i, pixel in zip( range(len(list_of_pixels)), list_of_pixels):

			# find the right part
			n = int(np.floor( i/divisor ))
			if n == n_parts: n -= 1

			list_of_files[n] = ("%s;%s" %(list_of_files[n], os.path.join(path, pixel)))

			for n in range(n_parts):
			arcpy.MosaicToNewRaster_management(list_of_files[n], target_path, target_files[n], None, "32_BIT_FLOAT", 2, 1)

			###### MERGE 20 PIECES TO ONE LARGE ONE ###############

			total_file = "DOM_merged"

			list_of_aggregated_files = ""

			for file in target_files:
			list_of_aggregated_files = ("%s;%s" %(list_of_aggregated_files, os.path.join(target_path, file)))

			arcpy.MosaicToNewRaster_management(list_of_aggregated_files, target_path, total_file, None, "32_BIT_FLOAT", 2, 1)