Latest Contributions by KONPETROV

Hi I have severous viewsheds (view1, view2, view3 etc) and I want to calculate the percent of visible areas per viewshed. I have the following code but I always get zero "0" at the second for loop and i have been trying that for hours, any help? Thanks in advance import arcpy import os from arcpy import env from arcpy.sa import * #workspace = arcpy.GetParameterAsText(0) arcpy.env.workspace = "C:Visibility/Calc_View/View" ras = arcpy.ListRasters("view*", "GRID") #cellsize = arcpy.GetParameterAsText(1) fields = ["VALUE", "COUNT"] for r in ras:     with arcpy.da.UpdateCursor(r, fields) as cursor:         total = 0         for row in cursor:             total = total + row[1]             print r, total         C = arcpy.da.UpdateCursor(r, fields, """"VALUE" = 1""")         for row in C:             percent = row[1] / total             print r, percent ... View more

Hello, I have the following code that I am trying to make it work as a tool. If I pass the paths to variables it works perfectly but when I pass them as Tool Parameters when I don't have any result. I don't know what to correct, but there must be some problem with the combination of tool parameter name + "the file a create". I tried also to import it to variable outside tool and then import the variable for use in the tool but nothing changed. This is the code: import arcpy import os from arcpy import env from arcpy.sa import * workspace = arcpy.GetParameterAsText(0) values = arcpy.GetParameterAsText(1) num_pnts = int(values) demodel = arcpy.GetParameterAsText(2) outReclass = Reclassify(demodel, "Value", RemapRange([[100,300,1],[301,500,0]]), "NODATA") outReclass.save(workspace + "rcls") n = 0 for i in range(num_pnts):     n += 1     # Spatially Balanced Points     arcpy.CreateSpatiallyBalancedPoints_ga(workspace + "rcls", 1, workspace + "point" + str(n))     inRaster = workspace + "rcls"     inObserverFeatures = "point" + str(n) + ".shp"     outViewshed = Viewshed(inRaster, inObserverFeatures)     outViewshed.save(workspace + "view" + str(n)) ... View more

Hello I have the following code (its a bit big but very simple), and the problem I am facing is that I cannot pass the raster rcls to the Create Spatially Balanced tool, unless I add it to Table of Contents (lines 32- 42). I am building a new tool, so I don't want to interrupt the proccess by adding something to the Talbe of Contents, cracking my tool to 2 tools. The "values" and "num_pnts" variable (line 13-14) show the number of points that will be created. Thanks in advance. import arcpy import os from arcpy import env from arcpy.sa import * #workspace = arcpy.GetParameterAsText(0) workspace = "C:/Results/ArcMap/Visibility/Calc_View/" # arithmos simion kai viewsheds #value = arcpy.GetParameterAsText(1) values = 2 num_pnts = int(values) #height_val = arcpy.GetParameterAsText(2) height_val = 100 height = int(height_val) #demodel = arcpy.GetParameterAsText(3) #demMAXResult = arcpy.GetRasterProperties_management(demodel, "MAXIMUM") #demMAXResult = arcpy.GetRasterProperties_management(demodel, "MINIMUM") demMAXResult = arcpy.GetRasterProperties_management("C:/Results/ArcMap/Visibility/Calc_View/tin2ras3g_1", "MAXIMUM") demMINResult = arcpy.GetRasterProperties_management("C:/Results/ArcMap/Visibility/Calc_View/tin2ras3g_1", "MINIMUM") demMAX = demMAXResult.getOutput(0) demMIN = demMINResult.getOutput(0) LIMIT = float(demMAX) - height # Reclassify outReclass = Reclassify("tin2ras3g_1", "Value", RemapRange([[float(demMIN),LIMIT,1],[LIMIT,demMAX,0]]), "NODATA") outReclass.save("C:/Results/ArcMap/Visibility/Calc_View/rcls") n = 0 for i in range(num_pnts):     n += 1     # Spatially Balanced Points     inRaster = "rcls"     outPoints = "C:/Results/ArcMap/Visibility/Calc_View/point" + str(n)     arcpy.CreateSpatiallyBalancedPoints_ga(inRaster, 1, outPoints) ... View more