| Commit message (Collapse) | Author | Age | Files |
|
|
|
|
|
|
|
|
|
|
| |
In some degenerate cases DelaunayTriangulation() returns a
GEOMETRYCOLLECTION where no triangle is fully contained in the input
geometry, which causes getRandomPoint() to fail as the list of
cumulative weights is an empty list.
We “overshoot” in that case and return the non-curated triangulation.
getRandomPoint(…, max_tries=1024) should be enough to eventually hit a
point in the input geometry.
|
| |
|
| |
|
|
|
|
|
| |
We always exit after a single WMS request if the layer has no SVG paths.
Don't print the request count in that case.
|
| |
|
| |
|
| |
|
|
|
|
| |
(Incl. failure.)
|
|
Unfortunately SGU/Bergsstaten doesn't offer layer files to download, but
it has an online webmap (WMS) at https://apps.sgu.se/kartvisare/kartvisare-mineralrattigheter.html
so we add a dedicated module to probe and fetch features from it.
Double checked that the resulting combination of GeoJSON files does not
result in data loss compared to the previous (private) script:
sort_features() { jq -S '.features |= sort_by(.properties.Name, .properties.Layer, .properties.Area)' }
diff -u --color=auto \
<(ogr2ogr -f GeoJSON -lco COORDINATE_PRECISION=2 /vsistdout/ $dir/Mineralrättigheter.gpkg \
SE.GOV.SGU.MRR.BEARBETNINGSKONCESSIONER_APPLIED_VY | jq '.name = "MRR:" + .name' | sort_features) \
<(ogr2ogr -f GeoJSON -lco COORDINATE_PRECISION=2 -nlt MULTIPOLYGON -nlt PROMOTE_TO_MULTI \
/vsistdout/ mrr/bearbetningskoncessioner_applied.geojson | sort_features)
(and similar for other layers).
|