Latest Contributions by NeilAyres

I think you are expecting "interpolation" to do something it is not designed to do. It will interpolate between low and high depending on the decay function. So, if you do not think this is "real", do as Dan_Patterson says, make a clip polygon and clip out the bit that you "like". ... View more

But what is the problem with what you have got? The raster interp seems to be doing exactly what is required. I take it that the white dots have very low values of whatever it is you are measuring and the blue dots are higher values. ... View more

"scallop shaped". ?? Can you show us a picture of what you mean? ... View more

The mxd is one thing. The data layers within it are another. These layers point to specific datasets somewhere in your syste. Are these in an enterprise db? SDE or something. If you want to "export" certain layers but only for a certain extent your will need to either use Clip as suggested above, or a graphical select. Using Select By Graphics—Help | ArcGIS Desktop  Clip—Help | ArcGIS Desktop  ... View more

If this is a mosaic dataset, have you built overviews. But generally that results in the crazy checkerboard view. ... View more

IOError: "Input_Union" This line you are supposed to edit to reflect the name of your Union feature class. Ditto the path for your environment. ... View more

Well this doesn't quite do what you want, but makes a condensed version of the union output with each poly now recording the count of overlaps and a list of the original polygons that went into them. import sys, os, arcpy from arcpy import env HomeDir = r"path to directory" scratchFgdb = os.path.join(HomeDir, "scratch.gdb") # input fgdb fc = "Input_Union" # input polygons (result of a union) env.workspace = scratchFgdb env.overwriteOutput = True # check the sr of the input desc = arcpy.Describe(fc) sr = desc.spatialReference poly_dict = {} # union has "ORIG_FID" as part of its output with arcpy.da.SearchCursor(fc, ["ORIG_FID", "SHAPE@", "SHAPE@XY"]) as Cur:     for row in Cur:         # read in the data         origFid = row[0]         geom = row[01]         cent = row[2]         # making a Key variable, just round off the centroid         x = round(cent[0], 1)         y = round(cent[1], 1)         Key = "{:.1f}/{:.1f}".format(x, y)         if Key in poly_dict:             poly_dict[Key][1] += 1             poly_dict[Key][2].append(origFid)         else:             poly_dict[Key] = [geom, 1, [origFid]] # output polygon outFc = "tmpOverlapPoly" arcpy.CreateFeatureclass_management("in_memory", outFc, "POLYGON", "", "", "", sr) arcpy.AddField_management(os.path.join("in_memory", outFc), "Key", "TEXT", "", "", 50) arcpy.AddField_management(os.path.join("in_memory", outFc), "COUNT", "SHORT") arcpy.AddField_management(os.path.join("in_memory", outFc), "FidList", "TEXT", "", "", 255) inCur = arcpy.da.InsertCursor(os.path.join("in_memory", outFc), ["Key", "COUNT", "FidList", "SHAPE@"]) for k, v in poly_dict.items():     Key = k     geom = v[0]     cnt = v[1]     fids = v[2]     fidStr = ",".join(str(f) for f in fids)     data = [Key, cnt, fidStr, geom]     inCur.insertRow(data) del inCur arcpy.CopyFeatures_management(os.path.join("in_memory", outFc), outFc)‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ Looks like this : ... View more

But, then again.... Would you want each overlap poly referring back to others. ie in your example... Poly #1 overlaps with Polys 2 & 3. Therefore Poly # 2 also overlaps 1.... Might be trickier than I think to get what you want. I have used in the past a fairly simple script just to output the polygon, with a count of the overlaps for that piece. ... View more

Well the start would be a Union in the Analysis / Overlay tool box. That gives you multiple polygon pieces overlying each other where there was once an overlap. After that, to get the output you want, I think you would have to do some arcpy / python scripting. Not too difficult. ... View more

Why is the extent all "NaN", or Nulls? Is the data getting published at all? ... View more

Where is the data - fgdb or SDE? Normally for publishing purposes (and not for editable features), we copy over from an enterprise system to a fgdb every night. That uses truncate and append. Doesn't affect the services based on these features. Not an expert on how all this works myself, just aware of the process. ... View more

Its exactly as the error says, there is no attribute or method with the GIS object for Tools. arcgis.gis module — arcgis 1.0.1 documentation  You should look at the geocoding module. arcgis.geocoding module — arcgis 1.0.1 documentation  ... View more

Generally when updating data which is used by a service we truncate then append the new data. That doesn't seem to upset things. Deleting is a no-no. ... View more

Is it possible that you are zoomed way out when looking at that service. Note that in ArcMap you are at a scale of 1:3000, the service has 1:650k+ ... View more

Why not just use Create Fishnet to do this? Create Fishnet—Help | ArcGIS Desktop  ... View more