Latest Contributions by stashavint

Is there a way I can create a flow map in ArcGis which shows movement from different point to one example if its an airport a flow map which shows all the traffic for Airport A from Airport B, Airport C, Airport D etc .. Thank you in advance. ... View more

I want to allow the user to choose several value from a list of values then later i want to use these values to iterate through a raster data set. But very time i keep getting errors saying the data set is not correct. My code is here: import arcpy,os inputraster = arcpy.GetParameterAsText(0) sea_level = arcpy.GetParameterAsText(1) fc = arcpy.GetParameterAsText(2) clipraster = arcpy.GetParameterAsText valuelist = [x.strip() for x in sea_level.split(";")] for value in valuelist:         arcpy.Clip_management(inputraster,"#",clipraster,fc,"#","ClippingGeometry","MAINTAIN_EXTENT")         clip_raster = arcpy.Raster(clipraster)         floodedarea = os.path.join(outputfolder + "floods" + str(value))         areasbelow = clip_raster <= value         areasbelow.save(floodedarea)         arcpy.AddMessage(areasbelow)‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

Hello, can anybody help with a python script that checks if cell A is connected to its neighboring cells and if it flooded then the cell connecting cells should be flooded too. I tried the code below but i cant seem to figure out how to in-cooperate my data into the script. My data include a DEM from which i will check inundated zone:   import heapq import numpy as np from scipy import ndimage def flood_fill(test_array, four_way=False): """""" input_array = np.copy(test_array) input_rows, input_cols = input_array.shape # Set h_max to a value larger than the array maximum to ensure #   that the while loop will terminate h_max = np.max(input_array*2.0) # Build mask of cells with data not on the edge of the image # Use 3x3 square structuring element inside_mask = ndimage.binary_erosion(     np.isfinite(input_array), structure=ndimage.generate_binary_structure(2, 2).astype(np.int)) # Initialize output array as max value test_array except edges output_array = np.copy(input_array) output_array[inside_mask] = h_max # Build priority queue and place edge pixels into priority queue # Last value is flag to indicate if cell is an edge cell put = heapq.heappush get = heapq.heappop fill_heap = [     (output_array[t_row, t_col], int(t_row), int(t_col), True)     for t_row, t_col in np.transpose(np.where(edge_mask))] heapq.heapify(fill_heap) def neighbors(row, col, four_way=False):     """Return indices of adjacent cells"""     if four_way:         return [             (row - 1, col), (row, col + 1),             (row + 1, col), (row, col - 1)]     else:         return [             (row - 1, col), (row - 1, col + 1),             (row, col + 1), (row + 1, col + 1),             (row + 1, col), (row + 1, col - 1),             (row, col - 1), (row - 1, col - 1)] # Iterate until priority queue is empty while True:     try:     h_crt, t_row, t_col, edge_flag = get(fill_heap) except IndexError:     break for n_row, n_col in neighbors(t_row, t_col, four_way):     # Skip cell if outside array edges     if edge_flag:         try:             if not inside_mask[n_row, n_col]:         continue except IndexError: ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

I have written a code to check if areas will be inundates when a certain value is added, the code doesn't give any errors but at the same time it does do anything. Just says running code complete, I am not even sure what the problem is. My code is here import arcpy import os, sys, traceback, shutil import numpy as np arcpy.env.overwriteOutput = True try:      outputfolder = arcpy.GetParameterAsText(0)      fc = arcpy.GetParameterAsText(1)      inputraster = arcpy.GetParameterAsText(2)      #check to see if the arcinfo license is on      if arcpy.CheckExtension("Spatial") == "Available":           arcpy.CheckOutExtension("Spatial")      else:           arcpy.AddError("Spatial analyst extension is not licensed.")           sys.exit()      #set the Progressor      #arcpy.SetProgessor("step","Checking input requirements....",0,9,1)      #Process :create an out put folder      arcpy.CreateFileGDB_management(outputfolder,"output.gdb")      outputgdb = outputfolder + "\\output.gdb"+ os.step      #create path to store raster data      clippedraster = outputfolder + "clippedraster"      flooded = outputfolder +"floodedarea"      #check to see if the inputs are in the same projection      desccs_ras = arcpy.Describe(inputraster).SpatialReference      desccs_fc = arcpy.Describe(fc).SpatialReference      if desccs_ras.name == "Unknown" and desccs_fc.name == "Unknown":           arcpy.AddMessage("Your data does not have a defined spatial reference")      else:           if desccs_ras.projectionName == desccs_fc.projectionName:                ""           else:                arcpy.AddError("The input Data should be in the same projection")                sys.exit()           del desccs_fc,desccs_ras      #remove the areas below sea level      currentRaster = arcpy.Raster(inputraster)      rastertoclip = Con(currentRaster>=0,currentRaster,1)      #clip the raster data to the interested zone extend      arcpy.Clip_management(rastertoclip,"#",clippedraster,fc,"#","ClippingGeometry","MAINTAIN_EXTENT")      clip_raster = arcpy.Raster("in_memory")      areasbelow = clip_raster <= 2      areasbelow.save(flooded)      md = arcpy.mapping.MapDocument("CURRENT")      df = arcpy.mapping.ListDataFrames(md)[0]      result = arcpy.MakeRasterLayer_management(flooded,"FloodedArea")      layer = result.getOutput(0)      arcpy.mapping.AddLayer(df, layer,"AUTO_ARRANGE") except:      tb = sys.exc_info()[2]      tbinfo = traceback.format_tb(tb)[0]      pymsg = "PYTHON ERRORS:\nTraceback Info:\n" + tbinfo + "\nError Info:\n    " + \             str(sys.exc_type)+ ": " + str(sys.exc_value) + "\n" ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

Hello everyone. is there a way to convert a polygon shape file into raster if it has values of precise decimal type. I used the polygon to raster tool and the produced raster does not have an attribute table and i cannot use the Build Raster attribute tool to build the table . ... View more