1,888 reputation
620
bio website n/a
location Tampa, FL
age 34
visits member for 2 years, 2 months
seen 7 hours ago
stats profile views 147

Currently working with ArcGIS, MODFLOW, and Python to develop local and regional groundwater flow models.


7h
comment While looping with arcpy: ERROR 010240: Could not save raster dataset
Are still getting exactly the same error message?
7h
revised While looping with arcpy: ERROR 010240: Could not save raster dataset
deleted 278 characters in body
7h
comment While looping with arcpy: ERROR 010240: Could not save raster dataset
I didn't read the documentation on the slope tool well enough...looks like the output_measurement variable you had was good, needs to be one of "DEGREE" or "PERCENT RISE". Sorry to confuse you there. Not sure where the problem is - can you update your question with your latest code? I will edit my answer to remove the part where it changes the output_measurement variable to a non-acceptable input value for the slope tool.
8h
comment While looping with arcpy: ERROR 010240: Could not save raster dataset
The slope tool may be trying to save the output to the out_measurement variable, even though you are working with it in terms of an object, out_slope. You can try setting the out_measurement variable to use memory space explicitly like this: r"in_memory\PERCENT_RISE_{0}".format(actual_timestep)
10h
revised While looping with arcpy: ERROR 010240: Could not save raster dataset
added 261 characters in body
10h
revised While looping with arcpy: ERROR 010240: Could not save raster dataset
added 277 characters in body
10h
comment While looping with arcpy: ERROR 010240: Could not save raster dataset
Ok, the last bit of the error where is says Could not save raster dataset to C:\SedMod\SedMod\Workspace.gdb\numpy_ras6 with output format FGDBR looks to be signalling the fact that you are trying to save multiple rasters to the same GDB with the same name ("PERCENT_RISE"). See help.arcgis.com/EN/arcgisdesktop/10.0/help/index.html#//…. You must either devise a way to come up with unique names (e.g. add your counter-value to the name using string formatting) or set acrpy.overwriteOutput = True.
10h
answered While looping with arcpy: ERROR 010240: Could not save raster dataset
11h
suggested suggested edit on While looping with arcpy: ERROR 010240: Could not save raster dataset
14h
comment Copying netCDF coordinate system
Hi Rich, thanks for the code. The only problem with going straight to ESRI grid is that my x,y locations the point of origin is unknown (e.g., 0,0; 4,4; 20,92, etc.). I have no way of placing these in a GCS without the projection info. I'll try the netcdf4 package from Gohlke's site to see if there are additional bits of information that I'm not getting with the scipy module.
Jun
14
comment Iterating through an numpy array and index to a value in another numpy array
I see you've used the shape property of the Use the sediment_transport_np array to create another array of the same dimensions, flow_direction_np. Use the shape property of the flow_direction_np array to make sure it has the same dimensions.
Jun
14
comment Iterating through an numpy array and index to a value in another numpy array
It looks like your flow direction array has different dimensions than one or more of the other two arrays which you are indexing with i and/or j.
Jun
14
comment How to travers fields of records without introducing them all in the function?
With my above comment, you also replace for field in fieldlist: with for i in range(len(fieldlist)):
Jun
14
comment How to travers fields of records without introducing them all in the function?
You could iterate over the list of fields using the index number and call something like row.getValue(fields[i].name)
Jun
13
revised Copying netCDF coordinate system
edited title
Jun
13
comment Copying netCDF coordinate system
I'm not entirely sure how to get the projection information from the netCDF files to define the projection for the resulting rasters. The FlatEarth variable value appears to be " " - this is the value returned when I call ncfile.variables['FlatEarth'][()]. The NOAA Weather and Climate Toolkit appears to do a bit of magic when you convert the raw NEXRAD files using that GUI - it converts 0-255 values to hourly rates of rainfall (for the digital precip product I have). Maybe it also applies the projection information as well?
Jun
12
comment Copying netCDF coordinate system
Thanks for the links.
Jun
12
revised Copying netCDF coordinate system
added 1727 characters in body
Jun
12
asked Copying netCDF coordinate system
Jun
10
comment 'Pausing' a Python script tool to allow the user to provide input before continuing
You should be able to achive this by strategically placing the raw_input function in your script. docs.python.org/2/library/functions.html#raw_input