topic How to create a point and see what it intersects with in Geoprocessing Questions https://community.esri.com/t5/geoprocessing-questions/how-to-create-a-point-and-see-what-it-intersects/m-p/350875#M12171 <HTML><HEAD></HEAD><BODY><SPAN>Can anyone see anything wrong with the following code?</SPAN><BR /><BR /><SPAN>I'm not getting any intersections even though I know that point is definitely within myFeatureLayer.</SPAN><BR /><BR /><SPAN>How does that deal with projection?</SPAN><BR /><BR /><BR /><PRE class="plain" name="code">pt = arcpy.Point(x,y) &nbsp;&nbsp;&nbsp; ptGeometry = arcpy.PointGeometry(pt) &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; cur = arcpy.SearchCursor(myFeaturelayer) &nbsp;&nbsp;&nbsp; for arow in cur: &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; geom = arow.shape &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if ptGeometry .touches(geom): &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.AddMessage("intersects") &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; break &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </PRE><DIV style="display:none;"> </DIV></BODY></HTML> Wed, 06 Oct 2010 02:30:23 GMT PeterTimmers 2010-10-06T02:30:23Z How to create a point and see what it intersects with https://community.esri.com/t5/geoprocessing-questions/how-to-create-a-point-and-see-what-it-intersects/m-p/350875#M12171 <HTML><HEAD></HEAD><BODY><SPAN>Can anyone see anything wrong with the following code?</SPAN><BR /><BR /><SPAN>I'm not getting any intersections even though I know that point is definitely within myFeatureLayer.</SPAN><BR /><BR /><SPAN>How does that deal with projection?</SPAN><BR /><BR /><BR /><PRE class="plain" name="code">pt = arcpy.Point(x,y) &nbsp;&nbsp;&nbsp; ptGeometry = arcpy.PointGeometry(pt) &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; cur = arcpy.SearchCursor(myFeaturelayer) &nbsp;&nbsp;&nbsp; for arow in cur: &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; geom = arow.shape &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if ptGeometry .touches(geom): &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.AddMessage("intersects") &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; break &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </PRE><DIV style="display:none;"> </DIV></BODY></HTML> Wed, 06 Oct 2010 02:30:23 GMT https://community.esri.com/t5/geoprocessing-questions/how-to-create-a-point-and-see-what-it-intersects/m-p/350875#M12171 PeterTimmers 2010-10-06T02:30:23Z Re: How to create a point and see what it intersects with https://community.esri.com/t5/geoprocessing-questions/how-to-create-a-point-and-see-what-it-intersects/m-p/350876#M12172 <HTML><HEAD></HEAD><BODY><SPAN>Could you try replacing ptGeometry.touches(geom) by ptGeometry.within(geom) or by ptGeometry.crosses(geom) because touches, according to the ArcGIS help, indicates if the boundaries of the geometries intersect. </SPAN><A href="http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//000v000000mv000000.htm">Within indicates if this geometry is contained within another geometry and crosses idicates if the two geometries intersect in a geometry of lesser dimension</A><SPAN>.</SPAN><BR /><BR /><SPAN>Succes,</SPAN><BR /><BR /><SPAN>Samuel</SPAN><BR /><A href="http://gissolved.blogspot.com">http://gissolved.blogspot.com</A></BODY></HTML> Wed, 06 Oct 2010 06:43:45 GMT https://community.esri.com/t5/geoprocessing-questions/how-to-create-a-point-and-see-what-it-intersects/m-p/350876#M12172 SamuelBosch 2010-10-06T06:43:45Z