Latest Contributions by MikeKnowles

I have the following for loop which subsets 420 DEM's into 20 DEM's.:   for i in range(0,420,21):     low = i     high = i + 21     interval = high / 21     AK_in = 'str_dems_AK' + str(interval)     AK_out = 'NED_AK' + str(interval)     subset_AK = AK_in + ' = str_dems_AK[' + str(low) + ':' + str(high) + ']'     exec subset_AK         arcpy.MosaicToNewRaster_management(AK_in, r"H:\dmf\pds\ned\grid\raster_mosaics", AK_out, "", "", "", "1", "", "")     arcpy.BuildPyramids_management(AK_out, "", "", "", "", "", "SKIP_EXISTING")   which results in the following error:   Traceback (most recent call last):   File "H:\NED_BuildPyramids_Mosaic_test.py", line 85, in <module> arcpy.MosaicToNewRaster_management(AK_in, r"H:\dmf\pds\ned\grid\raster_mosaics", AK_out, "", "", "", "1", "", "")   File "C:\Program Files (x86)\ArcGIS\Desktop10.1\arcpy\arcpy\management.py", line 11619, in MosaicToNewRaster     raise e ExecuteError: ERROR 000582: Error occurred during execution.   If I manually run the following command it works so I assume the problem is passing the value of “AK_out” into the mosaic command: arcpy.MosaicToNewRaster_management(AK_in, r"H:\dmf\pds\ned\grid\raster_mosaics",”NED_AK1”, "", "", "", "1", "", "")   If have tried a variety of alternatives but I always get an error so I'm hoping that someone can tell me how to (a) pass the value of AK_out into the mosiac command or (b) if that's not the problem, identify what is causing the eror.   Thanks, Mike. ... View more

I'm working with a NETCDF file and am running the following code and the nodata option of the CopyRaster command doesn't appear to be working as I expect.  The CopyRaster successfully creates the raster but in the properties the maximum value is still displayed as the value in the CopyRaster command.  I've also tried entering the 36-digit number instead of the scientific notation and I get the same result.  I'm not sure if it's the large value that's the problem or if I'm overlooking something. The nc variable I'm working with has the large value for water bodies and the remaining cells between 0-1 (floating point) so if the nodata option worked I'd expect a max <= 1. import arcpy from arcpy import env arcpy.CheckOutExtension("spatial") from arcpy.sa import * arcpy.env.overwriteOutput = True arcpy.env.workspace = r"G:\data\NETCDF" NetCDFFile = r"G:\data\NETCDF\year.nc" GDB = r"G:\data\NETCDF\GCM\Hist_Output_Part_Burn4.gdb" NetCDFFileProp = arcpy.CreateObject("NetCDFFileProperties", NetCDFFile) ##for i in range(0,int(NetCDFFileProp.getDimensionSize("year"))): for i in range(0,1): ## # convert month from a floating point to an integer (prevents error in naming raster) DIMENSION_VALUE = int(NetCDFFileProp.getDimensionValue("year",i)) print 'Working on Dimension Value: ', DIMENSION_VALUE, '(index ', i, ')' # delete existing raster if it exists if arcpy.Exists( GDB + "\\" + "year" + str(DIMENSION_VALUE )):   print '--- deleting ', "year" + str(DIMENSION_VALUE)   arcpy.Delete_management( GDB + "\\" + "year" + str(DIMENSION_VALUE )) # create raster layer from NetCDF file arcpy.MakeNetCDFRasterLayer_md(NetCDFFile, "part_burn", "lon", "lat", "year_"+str(i), "", "year " + str(DIMENSION_VALUE)) # create raster arcpy.CopyRaster_management("year_"+str(i), GDB + "\\" + "Hist_PB_"+str(DIMENSION_VALUE+1895), "", "", "9.96921e+036") #reclass raster #Notify when finished print "Finished!" Thanks, Mike ... View more

I have a co-worker that is running a Python script that continually quits responding (i.e., Task Manager shows that Python is Not Responding).  The script takes a NetCDF file and creates rasters for future analysis and when it quits responding it is at random rasters within the NetCDF file  (occurs in both PythonWin and IDLE).  We have run this same script/data on another "identical" PC (i.e., same make/model, RAM, and processing speed) as well as my PC (same make/different model, RAM, processing speed) and it successfully runs.  This implies, along with the random locations where the script quits responding, that the problem is not with the script or the data.  Does anyone have any thoughts on what might be causing Python to quit responding?  We have gone so far as to remove the PC from our network, disable the anti-virus (from workstation aml days), no other programs running. Running ArcGIS 10.0, Python 2.6.5 on Windows XP. Thanks, Mike ... View more

I have a script which I've added to a toolbox in ArcGIS 10.0.  The script asks the user for a geodatabase (I'm creating and storing rasters from a NetCDF file) and either uses an existing geodatabase or creates a new one if it doesn't exist.  In adding the script to a toolbox the problem I'm having is what parameter data type to select and/or what parameter properties need to be set to allow me to put in a geodatabase that doesn't exist. The script runs if an existing geodatabase is selected but not if a new one is selected. Thanks, Mike ... View more

I'm new to geoprocessing with Python (but with a basic knowledge of Python) and have been trying to loop through a shapefile to perform analysis on each polygon in the shapefile.  I have seen references to looping through directories. etc. but haven't been able to translate them to working with polygons within a shapefile. Thanks, Mike Using ArcGIS 9.3 on XP. ... View more

We have polygons for protected areas which are entirely contained within a study area polygon.  We would like to randomly move, and rotate, each protected areas polygon such that they do not overlap, and are fully contained, within the study area polygon.  We would like to repeat this randomization process x number of times. We currently have a corridor network based on the the existing configuration of protected areas and would like to compare this corridor network with the random protected areas configurations. Is this possible?  This is currently outside my expertise (but willing to learn) and am looking for direction/terms/concepts/etc. I should be aware of to get me started. Thanks, Mike ... View more

I have a point shapefile of nest sites from which I created a TIN and then converted the TIN to a TIN edge (3D Analyst -- Conversion -- From TIN -- TIN Edge).  This creates a line shapefile with no node identifiers and I would like to be able to identify the distance between the nodes (nest sites) that comprise each of the edge lines (I've added a length field).  Is there a procedure or set of steps that could be recommended to accomplish this? I have done this by exporting the tables for the point shapefile, TIN edge, and Generate Near table and analyzing in SAS.  Being new to 3D Analyst, I'm wondering if there is an alternative way(s) to do it within ArcMap? I've attached an image showing the point ID and the edge lines for a specific nest site. Running ArcMap 9.3.1 on Windows XP. Thanks, Mike ... View more