{"id":58060,"date":"2024-06-21T12:22:39","date_gmt":"2024-06-21T16:22:39","guid":{"rendered":"https:\/\/www.aaroads.com\/?p=58060"},"modified":"2024-11-11T10:21:00","modified_gmt":"2024-11-11T15:21:00","slug":"working-with-geometries-in-python","status":"publish","type":"post","link":"https:\/\/www.aaroads.com\/blog\/working-with-geometries-in-python\/","title":{"rendered":"Working with Geometries in Python"},"content":{"rendered":"

Module 6 introduces both working with geometries with Python and external text files or tabular data. Using geometry objects in memory can be used for input and output in geoprocessing. This also circumvents the need to make temporary feature classes. Utilizing the\u00a0open()<\/span>\u00a0function,\u00a0fileinput<\/span>\u00a0module, and\u00a0write<\/span>\u00a0module, text from .txt, .csv and other files can be imported into Python or exported to a file outside of ArcGIS Pro. The\u00a0readlines()<\/span>\u00a0and\u00a0writelines()<\/span>\u00a0methods also work with external files, but on a line by line basis, sometimes with iteration through a loop. These aspects further add to the utility of automating processes with Python.<\/p>\n

The Lab assignment for Module 6 is to write a Python script that creates a new text file outside of ArcGIS, and to populate it with geometry data from a feature class of polylines for rivers in Maui, Hawaii.<\/p>\n<\/div>

\"Maui,<\/a><\/span>

The Maui, Hawaii rivers feature class for Module 6<\/h4><\/div><\/div><\/div>

Following is a graphic I made in an attempt to better explain how to work with geometries in Python. The\u00a0getPart() method<\/span>, which returns an array of point objects for each geometry part, is the crux of the assignment.\u00a0<\/p>\n<\/div>

\"Diagram<\/a><\/span><\/div>

The SearchCursor function is used to read feature class geometries. Using the SearchCursor function requires several iterations of a feature class to return the points or vertices that make up the geometry of a feature. All polylines and polygons in a feature class are comprised of point objects.<\/p>\n

Every feature class is comprised of individual features, things we can see and interact with in ArcGIS Pro. Each feature has a record, which is tabulated in an attribute table. The rows of an attribute table contain the\u00a0records<\/b>\u00a0of every feature found within the feature class. The columns contain the attributes based upon a field.<\/p>\n

When using a SearchCursor function to work with geometries, the initial loop iterates through the records and returns a tuple with two objects for each feature. The first variable is the row index number and the second is the\u00a0polyline object<\/b>. A polyline object is a shape defined by one or more paths.<\/p>\n

A second iteration within the first iteration is required to ultimately return point objects from each polyline object. This is what is known as a\u00a0nested loop<\/i>. The\u00a0for loop<\/i>\u00a0that is the nested loop iterates through each polyline object and returns an\u00a0array object<\/b>. Part of the ArcPy package, arrays are groupings of points that comprise the geometry of a feature.\u00a0An array in other programming languages\u00a0 is a list.<\/i><\/p>\n

A second nested loop iterates through the array object to extract the\u00a0point object<\/b>, which contains the coordinate information of the vertices that make up each feature. The\u00a0getPart()<\/span>\u00a0method however allows a script to bypass this third iteration by coupling the operation with the previous\u00a0for loop<\/i>\u00a0that returns the array object of a feature. That is because the polyline object has this method, which receives an input parameter and index, and returns an array of point objects.<\/p>\n

The\u00a0getPart()<\/span>\u00a0method iterates through the tuple output by the\u00a0SearchCursor<\/span>\u00a0function. The index inputs the first element of the tuple, which for the SearchCursor is based upon the\u00a0field names<\/span>\u00a0argument.\u00a0Part<\/i>\u00a0in the function name references the return of an array object of point objects for a particular part of the geometry, which for Module 6 is index 0, the FID number. The method iterates off the polyline object returned from the tuple.<\/p>\n

The flowchart for the Module 6 script outlines the SearchCursor function and one nested loop used to obtain point objects for each part of the array object of a feature:<\/p>\n<\/div>

\"Nested<\/a><\/span><\/div>

This becomes a little more complicated with multipart features, where a feature is comprised of multiple arrays of point objects making up one set of attributes. Scenarios of this are where a feature contains polygons that are not physically connected, or where empty (interior) polygons fall within a larger (exterior) polygon. Think of St. Martin Parish, Louisiana, which consists of two areas separated by Iberia Parish in between, or Shively, Kentucky, a city encircled on all sides by Louisville.<\/p>\n

With multipart features, an iteration over all distinct parts comprising the overall feature is required to return its geometry. In summary, the polyline object returned from the SearchCursor contains array objects for each part. Within those array objects is another array which contains the point objects.<\/p>\n

The output of the script compiled for Module 6 produced a formatted text file for each river feature in the Maui rivers feature class:<\/p>\n<\/div>

\"Module<\/a><\/span><\/div>

Data output from point objects in the\u00a0writelines()<\/span>\u00a0method produced the ID number and the X-coordinate and Y-coordinate for each vertex of all records in the rivers feature class. This in turn was written to the new text file created at the beginning of the script with the\u00a0open()<\/span>\u00a0function.<\/p>\n

The vertex number for the vertices were assigned as part of the nested loop on the array object, where a variable\u00a0vertexnum<\/i>\u00a0increased by one for each pass. The\u00a0getPart()<\/span>\u00a0method also returned the name object of each vertex, which was derived from the SearchCursor iteration.<\/p>\n<\/div>

\"\u2018Ohe\u2019o<\/a><\/span>

\u2018Ohe\u2019o Gulch, downriver from Palikea Stream, which is one of the features from the Module 6 dataset<\/h4><\/div><\/div><\/div>

My wife and I visited Maui in January 2011 and drove the iconic Hana Highway to Haleakal\u0101 National Park. The above photo is one of \u2018Ohe\u2019o Gulch, which is connects Palikea Stream to the Pacific Ocean at the southeast end of the island.<\/p>\n<\/div><\/div><\/div><\/div><\/div>\n","protected":false},"excerpt":{"rendered":"","protected":false},"author":1,"featured_media":58061,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[1582],"tags":[1581,1602,1603,1601],"class_list":["post-58060","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-gis","tag-geographic-information-systems","tag-gis4102","tag-programming","tag-python"],"_links":{"self":[{"href":"https:\/\/www.aaroads.com\/blog\/wp-json\/wp\/v2\/posts\/58060","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.aaroads.com\/blog\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.aaroads.com\/blog\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.aaroads.com\/blog\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/www.aaroads.com\/blog\/wp-json\/wp\/v2\/comments?post=58060"}],"version-history":[{"count":3,"href":"https:\/\/www.aaroads.com\/blog\/wp-json\/wp\/v2\/posts\/58060\/revisions"}],"predecessor-version":[{"id":58068,"href":"https:\/\/www.aaroads.com\/blog\/wp-json\/wp\/v2\/posts\/58060\/revisions\/58068"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.aaroads.com\/blog\/wp-json\/wp\/v2\/media\/58061"}],"wp:attachment":[{"href":"https:\/\/www.aaroads.com\/blog\/wp-json\/wp\/v2\/media?parent=58060"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.aaroads.com\/blog\/wp-json\/wp\/v2\/categories?post=58060"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.aaroads.com\/blog\/wp-json\/wp\/v2\/tags?post=58060"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}