Latest Contributions by RandyBurton

I've done some experimenting with the code you found for creating polygon fishnets.  It adds some steps to clip the fishnet to the shape of the original polygons.  I did not explore the group options in the code, so there may need to be some adjustments for these options to work correctly in the modified code. The basic changes involve tracking the object IDs of the starting polygons and clipping the fishnet boxes using the geometry intersect method. In my test, I had 3 polygons that I wanted to cover with a 20x20 grid and clip that to the shapes of the polygons.  The results looked like: Here's the code starting around line 63 (following the set-up and functions in the code linked to in your question): # this would be line 63 following the functions and other set-up code # ## read in_features geometry into dictionary for use with geometry intersect method later origPolys = { r[0] : r[1] for r in arcpy.da.SearchCursor(in_features,['OID@','SHAPE@']) } # ## list comprehension # Process the data scratchName = arcpy.CreateScratchName("MBG","ByWidth","FeatureClass") arcpy.MinimumBoundingGeometry_management(in_features,scratchName,"RECTANGLE_BY_WIDTH",group_option,group_fields,"MBG_FIELDS") arcpy.DeleteField_management(scratchName,["MBG_Width","MBG_Length"]) # ## ["MBG_Width","MBG_Length","ORIG_FID"]) remove ORIG_ID from delete field list arcpy.CreateFeatureclass_management(os.path.dirname(out_features),\ os.path.basename(out_features),"POLYGON",\ template=scratchName,\ spatial_reference=sR) cur = arcpy.da.InsertCursor(out_features,["SHAPE@","MBG_Orientation","ORIG_FID"]) # ## ["SHAPE@","MBG_Orientation"]) add ORIG_ID to field list if labels == "LABELS": arcpy.CreateFeatureclass_management(os.path.dirname(out_features),\ os.path.basename(out_features)+"_label","POINT",\ spatial_reference=sR) for row in arcpy.da.SearchCursor(scratchName,["OID@","SHAPE@","MBG_Orientation","ORIG_FID"]): # ## ["OID@","SHAPE@","MBG_Orientation"]): add ORIG_ID to field list oid = row[0] ofid = row[3] # ## original feature id ( in_features objectid ) arcpy.AddMessage("\nProcessing polygon " + str(oid)) mbgPoly = row[1] orientation = float(row[2]) pivotPoint = mbgPoly.firstPoint origin_coord = str(pivotPoint.X) + ' ' + str(pivotPoint.Y) orthogonalPoly = RotatePolygon(mbgPoly,pivotPoint,orientation) opposite_corner = str(orthogonalPoly.extent.XMax) + ' ' + str(orthogonalPoly.extent.YMax) yCoord = str(pivotPoint.X) + ' ' + str(pivotPoint.Y + 1.0) scratchName2 = arcpy.CreateScratchName("OrthogonalFishNet",str(oid),"FeatureClass") arcpy.CreateFishnet_management(scratchName2,origin_coord,yCoord,cell_width,cell_height,number_rows,number_columns,labels=labels,template=orthogonalPoly.extent,geometry_type="POLYGON") #arcpy.CreateFishnet_management(scratchName2,origin_coord,yCoord,cell_width,cell_height,number_rows,number_columns,opposite_corner,labels=labels,geometry_type="POLYGON") # ## create a buffer from the original polygon which is in the dictionary polyBuff = origPolys[ofid].buffer(0) # ## buffer distance either 0 or negative number in default spatial reference units # arcpy.AddMessage("\npolyBuff = {}".format(polyBuff.JSON)) # ## for debugging for cell in arcpy.da.SearchCursor(scratchName2,["OID@","SHAPE@"]): arcpy.AddMessage("Processing fishnet cell..." + str(cell[0])) outPoly = RotatePolygon(cell[1],pivotPoint,0-orientation) # arcpy.AddMessage("\noutPoly = {}".format(outPoly.JSON)) # ## for debugging (json shows spatial reference as null) # ## since outPoly has a null spatial reference, we need to specify it for the intersect method newPoly = arcpy.FromWKT(outPoly.WKT,sR) # ## use WKT and add default spatial reference http://desktop.arcgis.com/en/arcmap/latest/analyze/arcpy-functions/fromwkt.htm # arcpy.AddMessage("\nnewPoly = {}".format(newPoly.JSON)) # ## for debugging # ## clip outPoly as required using geometry's intersect method bufPoly = newPoly.intersect(polyBuff,4) # ## http://desktop.arcgis.com/en/arcmap/latest/analyze/arcpy-classes/geometry.htm (see intersect method) # arcpy.AddMessage("\nbufPoly = {}".format(bufPoly.WKT)) # ## for debugging # ## it is possible that the polygon falls outside the intersect clip # ## json will show: bufPoly = {"rings":[],"spatialReference":{"wkid":102100,"latestWkid":3857}} and WKT will indicate EMPTY if "EMPTY" not in bufPoly.WKT: # ## add only if geometry is not empty (null) cur.insertRow([bufPoly,orientation,ofid])# ## [outPoly,orientation]) added field to save original polygon's objectid, see insertcursor around line 75 arcpy.Delete_management(scratchName2) if labels == "LABELS": del cur #Otherwise two InsertCursors in the same workspace requires arcpy.da.Editor lcur = arcpy.da.InsertCursor(out_features+"_label",["SHAPE@"]) for label in arcpy.da.SearchCursor(scratchName2+"_label",["SHAPE@"]): lcur.insertRow([arcpy.PointGeometry(RotatePoint(label[0].firstPoint,pivotPoint,0-orientation),sR)]) arcpy.Delete_management(scratchName2+"_label") del lcur cur = arcpy.da.InsertCursor(out_features,["SHAPE@","MBG_Orientation"]) arcpy.Delete_management(scratchName) del cur arcpy.AddMessage("\nDone")‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ Hope this helps. ... View more

I believe for the add query or update (OBJECTID field required) query, the features section is all that is needed.  You shouldn't need "fields" or the other sections. { "features": [{ "attributes": { "RECORD_ID": 1, "DATE_TIME": 1546359300000, "SPEED": 1, "HEADING": 12 }, "geometry": { "paths": [ [ [-5871083.00245447, 6028331.87508284], [-5871713.60793782, 6027605.72331413] ] ] } }] }‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ I did notice that the spatial reference of the line layer returned is 3857 - web mercator (see below).  From your first json example, it looks like you are trying to use 4326 - lon/lat.  You may need to convert your coordinates before the update query with a projectAs. "spatialReference": { "wkid": 102100, "latestWkid": 3857 }‍‍‍‍‍‍‍‍ ... View more

I think since you are using a file geodatabase you might need to use something like (you'll want to make sure both sides are a string of 3 characters): SUBSTRING("OWNERSHIP_NUMBER" FROM 1 FOR 3) = SUBSTRING("Layer" FROM 1 FOR 3)‍ (This from an old thread: Definition query SQL won't accept Left() Function) And according to SQL reference for query expressions used in ArcGIS All SQL used by the file geodatabase is based on the SQL-92 standard. ... View more