Perhaps you’re\u00a0looking for this?<\/p>\n
ogr2ogr -f “CSV” “E:\\4_GIS\\NorthArkCartoData\\UnitedStates\\MO_wkt” “E:\\4_GIS\\NorthArkCartoData\\UnitedStates\\USStates.shp” -sql ” SELECT * FROM usstates WHERE STATE_NAME = ‘Missouri’ ” -lco “GEOMETRY=AS_WKT ” -lco “LINEFORMAT=CRLF” -lco “SEPARATOR=SEMICOLON”<\/p><\/blockquote>\n
My buddy at\u00a0work needed a way to get the WKT geometry definition for a single feature in a shapefile. I thought, “surely this is something we can do with OGR?” Lo’ and behold, yes it was. \ud83d\ude42<\/p>\n
The script above\u00a0uses OGR\u00a0SQL to interrogate\u00a0a shapefile for one lone feature and, when it’s found, exports the record to a comma separated values (CSV) file (or in this case, a semicolon delimited file).\u00a0Here’s a quick break down:<\/p>\n
-f “CSV” “E:\\4_GIS\\NorthArkCartoData\\UnitedStates\\MO_wkt”<\/em>\u00a0<\/strong> This means CREATE the directory MO_wkt<\/em> and export\u00a0the derived output to this directory. I wrapped the file path in double-quotes (“) so linebreaks could not be introduced by the terminal and confuse the argument parser.<\/p>\n
GOTCHA! DO NOT create the MO_wkt<\/em> directory before running the script. Let OGR create that directory. However, relative to my example, E:\\4_GIS\\NorthArkCartoData\\UnitedStates\\<\/em> could be on your file system already.<\/p><\/blockquote>\n
“E:\\4_GIS\\NorthArkCartoData\\UnitedStates\\USStates.shp”\u00a0<\/em><\/strong> This is the shapefile you want to interrogate. Yet again, the file path is wrapped in double-quotes (“).<\/p>\n
-sql ”\u00a0 SELECT * FROM usstates WHERE STATE_NAME = ‘Missouri’\u00a0 “<\/em><\/strong>\u00a0\u00a0 This is the OGR SQL expression I used to grab the entire record (*) for features where the state_name field was ‘Missouri’. Notice how the search term (‘Missouri’) is wrapped in single-quotes (‘). If\u00a0the field you’re searching on is a\u00a0VARCHAR or a TEXT field\u00a0(obviously a name requires non-numeric characters), you’ll need to wrap your search term in single-quotes.\u00a0If you were searching on a numeric field, you would not\u00a0<\/em>wrap your search term in\u00a0quotes.\u00a0Also, just like\u00a0the file paths shown above, I recommend always<\/em> wrapping long values and expressions in double-quotes to avoid confusing the argument parser.\u00a0Next, and this is a matter of user-preference, but notice how I added extra spaces between the double-quotes and my SQL expression (i.e. ”\u00a0 SELECT..\u00a0 “). This is\u00a0totally legit and it may help to\u00a0distinguish the SQL query from the rest of the OGR script.\u00a0<\/p>\n
I could have exported a subset of the fields using an expression like:<\/p>\n
”\u00a0 SELECT STATE_ABBR as abbr,\u00a0STATE_NAME as name FROM usstates WHERE STATE_NAME = ‘Missouri’\u00a0 ” <\/em><\/p>\n
This expression selects only the STATE_NAME<\/em> and the STATE_ABBR<\/em> fields and renames them name<\/em> and abbr<\/em>, respectively, by applying the AS operator. The fields for the subset\u00a0are\u00a0comma-separated, but the last\u00a0field is not<\/em> followed by a comma\u00a0(this is basic SQL syntax, but I’m explaining it in case the readership includes some first-timers).\u00a0If you don’t want\u00a0to rename the fields in your selection, just omit the “AS rename” from your expression. —when exporting to CSV, the geometry\u00a0will be\u00a0exported automatically according to the geometry flag used, which we set next. Unfortunately, though, this means there’s nothing you can do to omit the geometry from your selection.<\/p>\n
These last few parameters use -lco<\/strong> flags\u00a0to\u00a0signalize various\u00a0“creation options” used by the CSV driver. The creation options\u00a0are unique for each dataset OGR can export, and in this case I’m demonstrating a particular recipie. If you want to know more about the\u00a0CSV creation options, check out the OGR CSV driver page<\/a>. Alternatively, if you want to export to a different file format, visit the OGR formats page<\/a>, select the appropriate driver, and read-up on the various creation options.<\/p><\/blockquote>\n
-lco “GEOMETRY=AS_WKT “<\/em><\/strong>\u00a0 This creation option means “Export the Geometry value\u00a0as Well Known Text”. Other options that might interest you include GEOMETRY=AS_XYZ<\/strong>, GEOMETRY=AS_XY<\/strong> or GEOMETRY=AS_YX<\/strong>.<\/p>\n
-lco “LINEFORMAT=CRLF”\u00a0<\/em><\/strong> This creation option means “use a Windows-formatted linebreak”. If you’re running Linux you should use “LINEFORMAT=LF”. Honestly, I can’t remember if we really needed this flag\u00a0so you may want to try skipping this option.<\/p>\n
-lco “SEPARATOR=SEMICOLON”\u00a0<\/em><\/strong> This flag specifies the delimiter used to separates fields and values from one another in the CSV. We chose SEMICOLON because Polygon WKT already has a bunch of commas in it, so using\u00a0the SEMICOLON delimiter helped us to visually identify the WKT we were really looking for. Your other options include SEPARATOR=COMMA and SEPARATOR=TAB.<\/p>\n
Alternatively, do this\u00a0with ogrinfo..<\/h2>\n
You could also\u00a0use ogrinfo<\/a> to dump a feature’s WKT right into the console. This approach required\u00a0using one of the so-called special fields<\/em> recognized by OGR SQL. To\u00a0learn\u00a0more about OGR SQL, its particulars, and\u00a0these special fields, visit the OGR SQL help page<\/a>.\u00a0Without further adieu, here’s my ogrinfo script:<\/p>\n
ogrinfo “E:\\4_GIS\\NorthArkCartoData\\UnitedStates\\USStates.shp” -geom=yes -sql ” SELECT OGR_GEOM_WKT<\/strong> FROM usstates WHERE STATE_NAME = ‘Missouri’ “<\/p><\/blockquote>\n
In this case, the OGR SQL expression requests only the “OGR_GEOM_WKT” special field<\/em>. Because we’re using ogrinfo, our output will stay in the console. However, this\u00a0was ultimately undesirable for our task, as “real world” Polygon WKT typically spills beyond the console buffer. Moreover, copying the WKT geometry out of the console introduced hundreds of unwanted linebreaks into the WKT. So, for these reasons, we prefered to export our WKT results to a CSV because it allowed us to easily see, select, and copy the desired WKT data without any fuss using good ol’ Notepad with word wrap.<\/p>\n
I hope you like. \ud83d\ude42<\/p>\n
\/Elijah<\/p>\n