Latest Contributions by DavidPike

Raster calculator the difference between the true DEM and your interpolated one, square it, the run something like zonal stats as table to grab the mean, then take the root. here's some code. I've not tested it. may have to save the raster objects as interim data. let me know if it doesn't work and i'll amend the code. import arcpy #path to your original accurate DEM - make sure path is enclosed by r' ' original_dem = r'c:\\myproject\mygeodatabase\originalDEM' #change the below paths to your interpolated DEM paths # - make sure path is enclosed by r' ' #interpolated DEM 1 interp_dem1 = r'c:\\myproject\mygeodatabase\interpDEM1' #interpolated DEM 2 interp_dem2 = r'c:\\myproject\mygeodatabase\interpDEM2' #interpolated DEM 3 interp_dem3 = r'c:\\myproject\mygeodatabase\interpDEM3' #interpolated DEM 4 interp_dem4 = r'c:\\myproject\mygeodatabase\interpDEM4' #interpolated DEM 5 interp_dem5 = r'c:\\myproject\mygeodatabase\interpDEM5' #make a list of the interpolated DEM paths interp_dem_list = [interp_dem1, interp_dem2, interp_dem3, interp_dem4, interp_dem5] #turn the original DEM into a raster object, #iterate through the list, turn paths into raster objects to perform #raster algebra on them, output raster cell values will be #(interp_dem - original_dem) ^ 2 original_dem_rasobj = arcpy.Raster(original_dem) #set a counter to label result with iteration number iteration_counter = 1 for interp_dem_path in interp_dem_list: ras_object = arcpy.Raster(interp_dem_path) ras_diff = (ras_object - original_dem_rasobj) ** 2 #get the mean of this raster ras_stats = arcpy.CalculateStatistics_management(ras_diff) ras_prop = arcpy.GetRasterProperties_management(ras_stats, "MEAN") ras_mean = ras_prop.getOutput(0) rmse = float(ras_mean) ** 0.5 #print out the rmse in the python window print("RMSE of Interpolated DEM Raster " + str(iteration_counter) + " = " + str(rmse)) #add 1 to the iteration counter iteration_counter += 1 print("Processing Complete") ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

You likely just need a return statement. Myscript() Return outputpath ... View more

You could add the xls as a gdb table then use search and update cursors to do the calculations. ... View more

Hi Stephen, There is often alot of confusion with tables which hold lat long figures. A table can exist but hold no geospatial location. Changing these values in the table will not create  or update the location of points. You need to add the table as xy data, then save it as a feature class or shapefile. Search for adding xy data in the help files or Google and you will find a trove of info. ... View more

I think the mxd needs to be in the datastore as you're setting the workspace then using list layers I guess? ... View more

Is the mxd in the data store and the data referenced in the mxd through a full network path? ... View more

Is this 1 mxd? Is it registered to the data store? ... View more

Maybe reclassify the 0 to NoData ... View more

What's the source of your rasters? ... View more

Other than setting the display to nearest neighbour I dont know. Would be good to find out if and when you solve it. Good luck! ... View more

And unique values set in the symbology? ... View more

Could be a DRA issue? Are they all set to the same symbology? ... View more

Hi Kathy. Date fields are python datetime objects which can be handy. Field 1 - Field 2     produces a TimeDelta object. TimeDelta.total_seconds()  will produce the number of seconds of the timedelta object. I believe the timedelta will take into account leap years so the result may be off several days from what you may expect. Total_seconds = (TimeField1 - TimeField2).total_seconds() SecInYear = (60*60*24*365.25) Years = Total_seconds/SecInYear This is very basic and prone to error. Hopefully someone can give a better answer ... View more