Latest Contributions by DavidPike

I would definitely project it into the metre coord system. The false eastings and northings let alone the metre foot difference mean the spatial location would be way off. What you mean by copy and paste I'm also not sure. ... View more

Hi Michelle, rampling can become quite confusing. I believe the function you're using is focal, ie it calculates from a moving window across all of your cells. I believe nlock statistics will suit your needs. How Block Statistics works—Help | ArcGIS for Desktop  ... View more

Then I just dont get why the field calculation isnt appropriate? ... View more

Are you saying you have multiple tables such as tracts and an ols table for each of them?  ... View more

Use arcpy.Raster? ... View more

So theres an OLS table for each row in the table in the second image? ... View more

Hi Jo,  Would you not want to 'chunk' your incidents rather than schools? ... View more

I'm sorry it's easily codeable but I still don't get the whole picture.  ... View more

If you send some pictures of the tables and a comprehensive explanation of what you want to achieve, I can write it out. ... View more

Is this your export web map service? If soRepublishe it with info level warnings and make sure its synchronous and you have added the get map template info result. Completely remove the widget, save, exit. Restart web editor and add the widget again. ... View more

Is this a default tool or one of your own? Perhaps resetting the normal.mxt might solve it. ... View more

Exactly, I would test this with another script tool with only a feature and the field map and see if you can grab a class attribute or type. Something like fieldmap = arcpy.GetParameter(1) property = fieldmap.inputFieldCount arcpy.AddMessage(property) ... View more

Would it be a field mapping object you would call as arcpy.GetParameter()? Good luck, field mapping objects give me nightmares. ... View more

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