def hillshade_generator(): arcpy.env.workspace = "C:\\Users\\michal\\Desktop\\_BP" arcpy.env.overwriteOutput = 1 print("Creating list of DEM files...") dem_pasy2 = [] for root, dirs, files in os.walk( r'C:\Users\michal\Desktop\_BP\gis\las_pasy_dem\arcgis\full'): for file in files: if file.endswith('_dem.tif'): pas = ("gis\\las_pasy_dem\\arcgis\\full\\" + file, "gis\\las_pasy_dem\\arcgis\\full\\hillshade\\" + file[:-8] + "_hillshade.tif", "gis\\las_pasy_dem\\arcgis\\area_n81\\hillshade\\" + file[:-8] + "_area81_hillshade.tif") dem_pasy2.append(pas) print("List of DEM files created!") for dem in dem_pasy2: print("Generating hillshade: ", dem[0]) arcpy.env.snapRaster = "Krkonose\\2019_06_14_Lucni_2_DSM.tif" arcpy.env.cellSize = "Krkonose\\2019_06_14_Lucni_2_DSM.tif" arcpy.env.extent = "MAXOF" arcpy.HillShade_3d(dem[0], dem[1]) arcpy.env.extent = "Krkonose\\area_n81.shp" arcpy.HillShade_3d(dem[0], dem[2]) print("Done!")
def make_topo_layers(topo_folder): """ Writes the layers :param topo_folder: We want to make layers for the stuff in this folder :return: """ source_code_folder = os.path.dirname(os.path.abspath(__file__)) symbology_folder = os.path.join(source_code_folder, 'BRATSymbology') dem_symbology = os.path.join(symbology_folder, "DEM.lyr") slope_symbology = os.path.join(symbology_folder, "Slope.lyr") hillshade_symbology = os.path.join(symbology_folder, "Hillshade.lyr") for folder in os.listdir(topo_folder): dem_folder_path = os.path.join(topo_folder, folder) dem_file = None for file_name in os.listdir(dem_folder_path): if file_name.endswith(".tif"): dem_file = os.path.join(dem_folder_path, file_name) make_layer(dem_folder_path, dem_file, "DEM", dem_symbology, is_raster=True) hillshade_folder = make_folder(dem_folder_path, "Hillshade") hillshade_file = os.path.join(hillshade_folder, "Hillshade.tif") try: arcpy.HillShade_3d(dem_file, hillshade_file) make_layer(hillshade_folder, hillshade_file, "Hillshade", hillshade_symbology, is_raster=True) except arcpy.ExecuteError as err: if get_execute_error_code(err) == "000859": arcpy.AddWarning( "Warning: Unable to create hillshade layer. Consider modifying your DEM input if you need a hillshade." ) else: raise arcpy.ExecuteError(err) slope_folder = make_folder(dem_folder_path, "Slope") slope_file = os.path.join(slope_folder, "Slope.tif") try: out_slope = arcpy.sa.Slope(dem_file) out_slope.save(slope_file) make_layer(slope_folder, slope_file, "Slope", slope_symbology, is_raster=True) except arcpy.ExecuteError as err: if get_execute_error_code(err) == "000859": arcpy.AddWarning( "Warning: Unable to create hillshade layer. Consider modifying your DEM input if you need a hillshade." ) else: raise arcpy.ExecuteError(err)
def calculateShadow(self): """ Program description: INPUT_PARAMETERS: inputValue_ - COMMENTS: """ #arcpy.CheckOutExtension("3D") print "Start creating shadow image..." inFileName = tkFileDialog.askopenfilename( defaultextension='TIFF', filetypes=[('ERDAS IMAGINE', '*.img'), ('TIFF', '*.tif')], initialdir=self.workspace, initialfile=INITFILE_DEM, multiple=False, parent=tkRoot, title= 'Chose input elevation raster file for processing shadow (ArcPy)') inFile = sa.Raster(inFileName) #Determine Shadow Area: Spatial Analyst Tools --> Solar Radiation: Area Solar Radiation #--> Not applicable here #3D Anaylist Tools --> Raster Surface: HillShade azimuth = float(143.1156527) #SUN_AZIMUTH from Landsat altitude = float(51.2588465) #SUN_ELEVATION from Landsat outFileName = "out_shadow.tif" # Problem length of filename! arcpy.HillShade_3d(inFile, outFileName, azimuth, altitude, "SHADOWS", 1) print "Hillshade shadow image '" + str( outFileName) + "' with azimut '" + str( azimuth) + "' and altitude '" + str(altitude) + "' created." return
def makeTopoLayers(topoFolder): """ Writes the layers :param topoFolder: We want to make layers for the stuff in this folder :return: """ sourceCodeFolder = os.path.dirname(os.path.abspath(__file__)) symbologyFolder = os.path.join(sourceCodeFolder, 'BRATSymbology') demSymbology = os.path.join(symbologyFolder, "DEM.lyr") slopeSymbology = os.path.join(symbologyFolder, "Slope.lyr") layers = [] for folder in os.listdir(topoFolder): demFolderPath = os.path.join(topoFolder, folder) demFile = None for fileName in os.listdir(demFolderPath): if fileName.endswith(".tif"): demFile = os.path.join(demFolderPath, fileName) layers.append( makeLayer(demFolderPath, demFile, os.path.basename(demFile), demSymbology, isRaster=True)) hillshadeFolder = makeFolder(demFolderPath, "Hillshade") hillshadeFile = os.path.join(hillshadeFolder, "Hillshade.tif") arcpy.HillShade_3d(demFile, hillshadeFile) makeLayer(hillshadeFolder, hillshadeFile, "Hillshade", isRaster=True) slopeFolder = makeFolder(demFolderPath, "Slope") slopeFile = os.path.join(slopeFolder, "Slope.tif") outSlope = arcpy.sa.Slope(demFile) outSlope.save(slopeFile) makeLayer(slopeFolder, slopeFile, "Slope", slopeSymbology, isRaster=True)
#PolyMask = str(workspaceGDB + "\\" + "Seg_" + inSegment + "_" + inYear + "_" + CellSizeName + "_mask_poly") PolyMask = str(workspaceGDB + "\\" + "Seg_" + inSegment + "_" + inYear + "_" + CellSizeName + "_mask_poly") arcpy.AddMessage("output polygon mask: " + PolyMask) arcpy.AddMessage(PolyMask) arcpy.RasterToPolygon_conversion(RasterMaskInt, PolyMask, "NO_SIMPLIFY", "Value") # ############################################### # Create hillshade # azimuth = 315 altitude = 45 modelShadows = "NO_SHADOWS" zFactor = 1 outHillshade = str(outFolder + "\\" + "Seg_" + inSegment + "_" + inYear + "_" + CellSizeName + "_hillshade.tif") arcpy.HillShade_3d(outRaster, outHillshade, azimuth, altitude, modelShadows, zFactor) ############################################### # ## Create data source polygon file and source raster file # # copy mb boundary to workspace gdb and assign short name inMBbasename = os.path.basename(inMBboundary) MBboundaryCopy = str(workspaceGDB + "\\" + inMBbasename) arcpy.CopyFeatures_management(inMBboundary, MBboundaryCopy) inMBbdy = os.path.basename(MBboundaryCopy) # assign short name to poly mask inPolyMask = os.path.basename(PolyMask) # copy topo boundary and assign short name TopoBoundaryFC = str(workspaceGDB + "\\" + "TopoBND_Seg_" + inSegment + "_" + inYear) arcpy.CopyFeatures_management(inTopoBoundaryShape, TopoBoundaryFC)
import arcpy from arcpy import env try: if arcpy.CheckExtension("3D") == "Available": arcpy.CheckOutExtension("3D") else: #Raise a custom exception raise LicenseError env.workspace = "D:/GrosMorne" arcpy.HillShade_3d("WesternBrook", "westbrook_hill", 300) arcpy.Aspect_3d("WesternBrook", "westbrook_aspect") except LicenseError: print "3D Analyst license is unavailable" except: print arcpy.GetMessages(2) finally: # Check in the 3D Analyst extension arcpy.CheckInExtension("3D")
#Check out extensions arcpy.CheckOutExtension("Spatial") arcpy.CheckOutExtension("3D") arcpy.env.workspace = gdbworkspace #convert img to TIFFs with elevation raster in feet and hillshade elev = arcpy.ListRasters("img*") for raster in elev: rasterObject = arcpy.Describe(raster) name = rasterObject.basename elevFeet = published_data + park_name + "_elev_feet.tif" hillshade = published_data + park_name + "_hillshade.tif" r1 = arcpy.sa.Raster(raster) * 3.281 r1.save(elevFeet) arcpy.HillShade_3d(elevFeet, hillshade, 270, 55) ################# EXTRA: Create NDVI and canopy layer ################# print """############# Creating NDVI and canopy rasters as TIFFs ############# """ #Create NDVI (extra stuff) naip = arcpy.ListRasters("NAIP*") for raster in naip: rasterObject = arcpy.Describe(raster) name = rasterObject.basename print name RedBand = published_data + rasterObject.basename + "_1.tif" NIRBand = published_data + rasterObject.basename + "_4.tif" compRed = rasterObject.catalogPath + "/Band_1" print "Extracting Near Infrared Band to {0}".format(RedBand)
# This sets the workspace and shows the raster files in the database try: arcpy.env.workspace = geodatabase rasterList = arcpy.ListRasters() print ("Raster: ", rasterList) except: print("Not a valid database") arcpy.env.overwriteOutput = True # This asks the user to choose what layer they want to use to create the hillshade in_raster = input("Please choose what Raster Layer you want to use: ") # This defines the variables that will be used to create the hillshade out_dir = geodatabase i = 90 azimuth = 90 altitude = 55 # This creates a while loop that will create the 3 hillshades # The loop will run as long as azimuth is less that 271 while azimuth < 271: # This determines what the hillshade output will be called out_raster = "hillshade" + str(i) print("Creating hillshade from " + in_raster + " ...") # This runs the arcpy hillshade function arcpy.HillShade_3d(in_raster, out_raster, azimuth, altitude, "SHADOWS", 1) # This adds 90 to the variabke azimuth azimuth += 90 # This adds 90 to the variable i i += 90
polygon = arcpy.Polygon(arcpy.Array(boundaryPoints)) arcpy.CreateFeatureclass_management(arcpy.env.workspace, name + '_Boundary.shp', 'POLYGON', polygon) arcpy.CopyFeatures_management(polygon, name + '_Boundary.shp') arcpy.AddMessage(name + '_Boundary.shp generated.') # Create TIN based on topo data and boundary arcpy.CreateTin_3d(name + '_TIN') arcpy.EditTin_3d( name + '_TIN', [[name + '_Points.shp', 'Z', '<None>', 'Mass_Points', False], [name + '_Boundary.shp', '<None>', '<None>', 'Soft_Clip', False]]) arcpy.AddMessage(name + '_TIN generated.') # Create raster out of TIN arcpy.TinRaster_3d(name + '_TIN', name + '_Raster.tif', 'FLOAT', 'LINEAR', sample_distance='CELLSIZE {}'.format(cellsize), z_factor=1) arcpy.AddMessage(name + '_Raster.tif generated.') # Create hillshade out of raster arcpy.HillShade_3d(name + '_Raster.tif', name + '_Hillshade.tif', azimuth, altitude, z_factor=1) arcpy.AddMessage(name + '_Hillshade.tif generated.')
Flow_Velocity = r"C:\Users\NReid\Documents\ArcGIS\Projects\TrailProjectModel\modelData.gdb\MaxFlow_Vel" TC_Reclas_in = r"C:\Users\NReid\Documents\ArcGIS\Projects\TrailProjectModel\modelData.gdb\TC_Reclassify" TrailPT_OVel = r"C:\Users\NReid\Documents\ArcGIS\Projects\TrailProjectModel\modelData.gdb\Excess_Vel" Point_Vel = r"C:\Users\NReid\Documents\ArcGIS\Projects\TrailProjectModel\modelData.gdb\Velocity_pt" # Process: Fill - works fill = arcpy.sa.Fill(DEM) fill.save(Fill_DEM) print ("Fill complete.") # Process: Flow Direction - works flow = arcpy.gp.FlowDirection_sa(Fill_DEM, DEM_FlowDir, "Normal", "D8") print ("Flow Direction complete.") # Process: Hillshade - works hshade = arcpy.HillShade_3d(Fill_DEM, Hillshade, "315", "45", "NO_SHADOWS", "1") print ("Hillshade complete.") # Process: Slope Raster slope = arcpy.Slope_3d(Fill_DEM, Slope_Percent, "PERCENT_RISE", 0.3048, "PLANAR", "") print ("Slope Raster complete.") # Process: Aspect Raster aspect = arcpy.sa.Aspect(Fill_DEM, "PLANAR", "") aspect.save (DEM_Aspect) print ("Aspect Raster complete.") # Process: Raster Domain - works arcpy.RasterDomain_3d(Fill_DEM, DEM_Domain, "POLYGON") print ("Raster Domain complete.")
def terreno(self, path): # Paso 1: # TERRENO: DEM, DIRECCION DE FLUJO, ACUMULACION, LONGITUD t_inputs = self.inputs() print('TERRENO: DEM, DIRECCION DE FLUJO, ACUMULACION, LONGITUD') # Crear carpeta de salidas print('Crear carpeta de salidas') os.system('md modeloFONAG\\%s' % t_inputs[16]) os.system('md modeloFONAG\\%s\\tables' % t_inputs[16]) env.workspace = '%s/modeloFONAG/%s' % (path, t_inputs[16]) # Lectura de DEM y relleno de vacios r_DEM = arcpy.Raster('%s/modeloFONAG/inputs/%s' % (path, t_inputs[0])) print('Lectura DEM') print('%s/modeloFONAG/inputs/%s' % (path, t_inputs[0])) # Leer Mascara inMaskData = '%s/modeloFONAG/inputs/%s/%s' % (path, t_inputs[1], t_inputs[2]) print('Lectura Mascara') print(inMaskData) # Extraer por mascara r_DEM_mask = ExtractByMask(r_DEM, inMaskData) r_DEM_mask.save('%s/modeloFONAG/%s/r_DEM_mask.tif' % (path, t_inputs[16])) print('Extraer por mascara') # Relleno de vacios r_SurfaceRaster = Fill(r_DEM_mask) r_SurfaceRaster.save('%s/modeloFONAG/%s/r_DEM_fill.tif' % (path, t_inputs[16])) print('Relleno de vacios') # Variables de entorno env.cellSize = '%s/modeloFONAG/%s/r_DEM_fill.tif' % (path, t_inputs[16]) env.extent = '%s/modeloFONAG/%s/r_DEM_fill.tif' % (path, t_inputs[16]) print('Variables de entorno') # Generar raster de ceros y de unos para ciertos procesos r_DEM_ceros = r_SurfaceRaster * 0 r_DEM_unos = r_SurfaceRaster / r_SurfaceRaster # r_DEM_ceros.save('%s/modeloFONAG/%s/r_DEM_ceros.tif' % (path, t_inputs[16])) # r_DEM_unos.save('%s/modeloFONAG/%s/r_DEM_unos.tif' % (path, t_inputs[16])) print('Extraer por mascara') # Genera hillshade r_Hillshade = '%s/modeloFONAG/%s/r_DEM_shd.tif' % (path, t_inputs[16]) arcpy.HillShade_3d(r_SurfaceRaster, r_Hillshade, z_factor=3) print('Generar hillshade') # Direccion de flujo r_FlowDirection = FlowDirection(r_SurfaceRaster, 'NORMAL') r_FlowDirection.save('%s/modeloFONAG/%s/r_DEM_dir.tif' % (path, t_inputs[16])) print('Direccion de flujo') # Acumulacion de flujo inWeightRaster = '' r_FlowAccumulation = FlowAccumulation(r_FlowDirection, inWeightRaster) + 1 r_FlowAccumulation.save('%s/modeloFONAG/%s/r_DEM_acc.tif' % (path, t_inputs[16])) print('Acumulacion de flujo') # Area de aporte cellXResult = arcpy.GetRasterProperties_management(r_DEM_mask, 'CELLSIZEX') cellYResult = arcpy.GetRasterProperties_management(r_DEM_mask, 'CELLSIZEY') cellX = float(cellXResult.getOutput(0)) cellY = float(cellYResult.getOutput(0)) cellArea = cellX * cellY r_area = r_FlowAccumulation * cellArea r_area.save('%s/modeloFONAG/%s/r_DEM_area.tif' % (path, t_inputs[16])) print('Area de aporte') # Longitud de flujo r_FlowLength = FlowLength(r_FlowDirection, 'UPSTREAM', inWeightRaster) r_FlowLength.save('%s/modeloFONAG/%s/r_DEM_length.tif' % (path, t_inputs[16])) print('Longitud de flujo') r_LengthAccumulation = FlowAccumulation(r_FlowDirection, r_FlowLength) r_Length_mean = r_LengthAccumulation / r_FlowAccumulation r_Length_norm = (r_FlowLength + 1) / (r_Length_mean + 1) c_factordist = Exp(-1 * r_Length_norm) c_factordist.save('%s/modeloFONAG/%s/r_factordist.tif' % (path, t_inputs[16])) print('Factor de ponderacion por distancia') print('FIN DE MODULO DE TERRENO')
arcpy.CreateTin_3d(output_folder + "\\tin") #edit tin arcpy.EditTin_3d(output_folder + "\\tin", spot + ";" + contours + ";" + breaklines) #tin to raster conversion data_type = "FLOAT" interp_method = "LINEAR" samp_method = "CELLSIZE 10" factor = 1 arcpy.TinRaster_3d(output_folder + "\\tin", output_folder + "\\grid", data_type, interp_method, samp_method, factor) #slope computation arcpy.Slope_3d(output_folder + "\\grid", output_folder + "\\slope", "DEGREE", 1) #hillshading arcpy.HillShade_3d() arcpy.CheckInExtension('3D') else: print('Licence not available') #------------------------------------------------------------------------------- # PRINT MAPS #------------------------------------------------------------------------------- import arcpy import arcpy.mapping as map path_mxd = r"C:\DAS\SESSION4\template2.mxd" #mxd reference mxd = map.MapDocument(path_mxd)
def Model1(): # Model 1 # To allow overwriting outputs change overwriteOutput option to True. arcpy.env.overwriteOutput = False # Check out any necessary licenses. arcpy.CheckOutExtension("3D") arcpy.ImportToolbox(r"F:\software\LAStools\ArcGIS_toolbox\LAStools.tbx") # Model Environment settings with arcpy.EnvManager(cellSize="0,05", snapRaster=r"KRNAP\2019_06_14_Lucni_2_DSM.tif"): jtsk_pas3_area81_las_2_ = "C:\\Users\\michal\\Desktop\\_BP\\Krkonose\\LiDAR\\area_81\\jtsk_pas3_area81.las" grid_size = 0, 05 keep = "lowest" jtsk_pas4_area81_las_2_ = "C:\\Users\\michal\\Desktop\\_BP\\Krkonose\\LiDAR\\area_81\\jtsk_pas4_area81.las" jtsk_pas5_area81_las = "C:\\Users\\michal\\Desktop\\_BP\\Krkonose\\LiDAR\\area_81\\jtsk_pas5_area81.las" jtsk_pas6_area81_las = "C:\\Users\\michal\\Desktop\\_BP\\Krkonose\\LiDAR\\area_81\\jtsk_pas6_area81.las" jtsk_pas7_area81_las = "C:\\Users\\michal\\Desktop\\_BP\\Krkonose\\LiDAR\\area_81\\jtsk_pas7_area81.las" jtsk_pas8_area81_las = "C:\\Users\\michal\\Desktop\\_BP\\Krkonose\\LiDAR\\area_81\\jtsk_pas8_area81.las" # Process: lasthin (lasthin) lasthin(input_file=jtsk_pas3_area81_las_2_, grid_size=grid_size, keep=keep, output_format="", output_file=jtsk_pas3_area81_thin_min_las, output_directory="", output_appendix="", additional_command - line_parameters="", verbose=True) # Process: las2dem (las2dem) las2dem(input_file=jtsk_pas3_area81_thin_min_las, step=0, 05, kill_triangles_with_edges_longer_than=100, item="elevation", output="actual values", sun_direction_for_hillside_shading="north east", sun_position_for_hillside_shading="1 pm", set_min__for_gray_ramp_and_false_colors_=None, set_max__for_gray_ramp_and_false_colors_=None, triangulate="all points", use_tile_bounding_box=False, closed_breaklines="", open_breaklines="", output_format="tif", output_file=jtsk_pas3_area81_dem_tif, output_directory="", output_appendix="", additional_command - line_parameters="", verbose=True) # Process: HillShade (HillShade) jtsk_pas3_area81_hillshade_tif = "C:\\Users\\michal\\Desktop\\_BP\\gis\\las_pasy_dem\\lastools\\area_n81\\hillshade\\jtsk_pas3_area81_hillshade.tif" arcpy.HillShade_3d(in_raster=jtsk_pas3_area81_dem_tif, out_raster=jtsk_pas3_area81_hillshade_tif, azimuth=315, altitude=45, model_shadows="NO_SHADOWS", z_factor=1) # Process: lasthin (2) (lasthin) lasthin(input_file=jtsk_pas4_area81_las_2_, grid_size=grid_size, keep=keep, output_format="", output_file=jtsk_pas4_area81_thin_min_las, output_directory="", output_appendix="", additional_command - line_parameters="", verbose=True) # Process: las2dem (2) (las2dem) las2dem(input_file=jtsk_pas4_area81_thin_min_las, step=0, 05, kill_triangles_with_edges_longer_than=100, item="elevation", output="actual values", sun_direction_for_hillside_shading="north east", sun_position_for_hillside_shading="1 pm", set_min__for_gray_ramp_and_false_colors_=None, set_max__for_gray_ramp_and_false_colors_=None, triangulate="all points", use_tile_bounding_box=False, closed_breaklines="", open_breaklines="", output_format="tif", output_file=jtsk_pas4_area81_dem_tif, output_directory="", output_appendix="", additional_command - line_parameters="", verbose=True) # Process: HillShade (2) (HillShade) jtsk_pas4_area81_hillshade_tif = "C:\\Users\\michal\\Desktop\\_BP\\gis\\las_pasy_dem\\lastools\\area_n81\\hillshade\\jtsk_pas4_area81_hillshade.tif" arcpy.HillShade_3d(in_raster=jtsk_pas4_area81_dem_tif, out_raster=jtsk_pas4_area81_hillshade_tif, azimuth=315, altitude=45, model_shadows="NO_SHADOWS", z_factor=1) # Process: lasthin (3) (lasthin) lasthin(input_file=jtsk_pas5_area81_las, grid_size=grid_size, keep=keep, output_format="", output_file=jtsk_pas5_area81_thin_min_las, output_directory="", output_appendix="", additional_command - line_parameters="", verbose=True) # Process: las2dem (3) (las2dem) las2dem(input_file=jtsk_pas5_area81_thin_min_las, step=0, 05, kill_triangles_with_edges_longer_than=100, item="elevation", output="actual values", sun_direction_for_hillside_shading="north east", sun_position_for_hillside_shading="1 pm", set_min__for_gray_ramp_and_false_colors_=None, set_max__for_gray_ramp_and_false_colors_=None, triangulate="all points", use_tile_bounding_box=False, closed_breaklines="", open_breaklines="", output_format="tif", output_file=jtsk_pas5_area81_dem_tif, output_directory="", output_appendix="_dem", additional_command - line_parameters="", verbose=True) # Process: HillShade (3) (HillShade) jtsk_pas5_area81_hillshade_tif = "C:\\Users\\michal\\Desktop\\_BP\\gis\\las_pasy_dem\\lastools\\area_n81\\hillshade\\jtsk_pas5_area81_hillshade.tif" arcpy.HillShade_3d(in_raster=jtsk_pas5_area81_dem_tif, out_raster=jtsk_pas5_area81_hillshade_tif, azimuth=315, altitude=45, model_shadows="NO_SHADOWS", z_factor=1) # Process: lasthin (4) (lasthin) lasthin(input_file=jtsk_pas6_area81_las, grid_size=grid_size, keep=keep, output_format="", output_file=jtsk_pas6_area81_thin_min_las, output_directory="", output_appendix="", additional_command - line_parameters="", verbose=True) # Process: las2dem (4) (las2dem) las2dem(input_file=jtsk_pas6_area81_thin_min_las, step=0, 05, kill_triangles_with_edges_longer_than=100, item="elevation", output="actual values", sun_direction_for_hillside_shading="north east", sun_position_for_hillside_shading="1 pm", set_min__for_gray_ramp_and_false_colors_=None, set_max__for_gray_ramp_and_false_colors_=None, triangulate="all points", use_tile_bounding_box=False, closed_breaklines="", open_breaklines="", output_format="tif", output_file=jtsk_pas6_area81_dem_tif, output_directory="", output_appendix="_dem", additional_command - line_parameters="", verbose=True) # Process: HillShade (4) (HillShade) jtsk_pas6_area81_hillshade_tif = "C:\\Users\\michal\\Desktop\\_BP\\gis\\las_pasy_dem\\lastools\\area_n81\\hillshade\\jtsk_pas6_area81_hillshade.tif" arcpy.HillShade_3d(in_raster=jtsk_pas6_area81_dem_tif, out_raster=jtsk_pas6_area81_hillshade_tif, azimuth=315, altitude=45, model_shadows="NO_SHADOWS", z_factor=1) # Process: lasthin (5) (lasthin) lasthin(input_file=jtsk_pas7_area81_las, grid_size=grid_size, keep=keep, output_format="", output_file=jtsk_pas7_area81_thin_min_las, output_directory="", output_appendix="", additional_command - line_parameters="", verbose=True) # Process: las2dem (5) (las2dem) las2dem(input_file=jtsk_pas7_area81_thin_min_las, step=0, 05, kill_triangles_with_edges_longer_than=100, item="elevation", output="actual values", sun_direction_for_hillside_shading="north east", sun_position_for_hillside_shading="1 pm", set_min__for_gray_ramp_and_false_colors_=None, set_max__for_gray_ramp_and_false_colors_=None, triangulate="all points", use_tile_bounding_box=False, closed_breaklines="", open_breaklines="", output_format="tif", output_file=jtsk_pas7_area81_dem_tif, output_directory="", output_appendix="_dem", additional_command - line_parameters="", verbose=True) # Process: HillShade (5) (HillShade) jtsk_pas7_area81_hillshade_tif = "C:\\Users\\michal\\Desktop\\_BP\\gis\\las_pasy_dem\\lastools\\area_n81\\hillshade\\jtsk_pas7_area81_hillshade.tif" arcpy.HillShade_3d(in_raster=jtsk_pas7_area81_dem_tif, out_raster=jtsk_pas7_area81_hillshade_tif, azimuth=315, altitude=45, model_shadows="NO_SHADOWS", z_factor=1) # Process: lasthin (6) (lasthin) lasthin(input_file=jtsk_pas8_area81_las, grid_size=grid_size, keep=keep, output_format="", output_file=jtsk_pas8_area81_thin_min_las, output_directory="", output_appendix="", additional_command - line_parameters="", verbose=True) # Process: las2dem (6) (las2dem) las2dem(input_file=jtsk_pas8_area81_thin_min_las, step=0, 05, kill_triangles_with_edges_longer_than=100, item="elevation", output="actual values", sun_direction_for_hillside_shading="north east", sun_position_for_hillside_shading="1 pm", set_min__for_gray_ramp_and_false_colors_=None, set_max__for_gray_ramp_and_false_colors_=None, triangulate="all points", use_tile_bounding_box=False, closed_breaklines="", open_breaklines="", output_format="tif", output_file=jtsk_pas8_area81_dem_tif, output_directory="", output_appendix="_dem", additional_command - line_parameters="", verbose=True) # Process: HillShade (6) (HillShade) jtsk_pas8_area81_hillshade_tif = "C:\\Users\\michal\\Desktop\\_BP\\gis\\las_pasy_dem\\lastools\\area_n81\\hillshade\\jtsk_pas8_area81_hillshade.tif" arcpy.HillShade_3d(in_raster=jtsk_pas8_area81_dem_tif, out_raster=jtsk_pas8_area81_hillshade_tif, azimuth=315, altitude=45, model_shadows="NO_SHADOWS", z_factor=1) # Process: lasmerge (lasmerge) jtsk_pas3_8_area81_las = "C:\\Users\\michal\\Desktop\\_BP\\Krkonose\\LiDAR\\area_81\\jtsk_pas3-8_area81.las" lasmerge(input_file_1=jtsk_pas3_area81_las_2_, input_file_2=jtsk_pas4_area81_las_2_, input_file_3=jtsk_pas5_area81_las, input_file_4=jtsk_pas6_area81_las, input_file_5=jtsk_pas7_area81_las, input_file_6=jtsk_pas8_area81_las, input_file_7="", input_file_8="", input_file_9="", output_file=jtsk_pas3_8_area81_las, additional_command - line_parameters="", verbose=True)