Skip to content

Instantly share code, notes, and snippets.

@palmerj

palmerj/get_centroid.bash

Last active December 20, 2017 04:47
Show Gist options
  • Save palmerj/07df8f09f389f01cb6b39dd66f93ef3d to your computer and use it in GitHub Desktop.
Save palmerj/07df8f09f389f01cb6b39dd66f93ef3d to your computer and use it in GitHub Desktop.
Download ZIP
OGR Create JSON centroid for labelling
Raw
get_centroid.bash
ogr2ogr -f GeoJSON labels.json WFS:"https://data.linz.govt.nz/services;key=$LDS_API_KEY/wfs/layer-51572" -dialect sqlite -sql 'SELECT id, lease_name AS name, ST_Centroid(geometry) from "data.linz.govt.nz:layer-51572"' -lco COORDINATE_PRECISION=4 -t_srs "EPSG:4326"
@palmerj
Copy link
Author

palmerj commented Dec 20, 2017 •
edited
Loading

Even better way to remove GeoJSON schema with jq:

ogr2ogr -f GeoJSON /vsistdout WFS:"https://data.linz.govt.nz/services;key=$LDS_API_KEY/wfs/layer-51572" -dialect sqlite -sql "SELECT id, lease_name AS name, ST_Y(ST_Transform(ST_Centroid(geometry), 4326)) || ', ' || ST_X(ST_Transform(ST_Centroid(geometry), 4326)) AS latlon FROM \"data.linz.govt.nz:layer-51572\"" | jq "[.features[].properties]" > 51572.json

ogr2ogr -f GeoJSON /vsistdout WFS:"https://data.linz.govt.nz/services;key=$LDS_API_KEY/wfs/layer-53564" -dialect sqlite -sql "SELECT napalis_id AS id, name, ST_Y(ST_Transform(ST_Centroid(geometry), 4326)) || ', ' || ST_X(ST_Transform(ST_Centroid(geometry), 4326)) AS latlon FROM \"data.linz.govt.nz:layer-53564\"" | jq "[.features[].properties]" > 53564.json

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment