[Talk-se] Naturvårdsverkets nya Nationella MarktäckeData

[Grigory Rechistov]

Hej, ikväll hade jag mer framgång med QGIS som hjälpverktyg och GeoJSON som vektorformat som QGIS kan skriva och JOSM kan läsa. Då lyckades jag se mutipolygoner med etiketter i JOSM. Med "ESRI shapefile" som mellanformatet hade jag samma problem som andra — etiketter överlevde inte exporteringsprocessen.

Så här gjorde jag.

1. Körde QGIS, öppnade nmd2018bas_ogeneraliserad_v1_0.tif som rasterlager.
2. Klippte en mindre rektangel från hela rasterlagret för experiment och sparade den som en ny fil. Man kan dra i QGIS kartfönstret för att välja klippgränsen, då körs gdal_translate som backend.
3. Öppnade klippet i QGIS och konverterade det till vektor. Då körs gdal_polygonize.py och man kan redigera kommandoraden innan det körs. Sen kunde jag se det nya vektorlagret i QGIS och alla linjer hade etiketter typ "DN=111".
4. Sedan kunde jag spara vektorlagret som en GeoJSON fil med WGS84 koordinatsystemet. Sedan öppnade jag den i JOSM (det krävdes att installera en insticksmodul i JOSM) och kunde se alla linjer och multipolygoner på rätt plats och med samma taggar "DN=111" i JOSM!

Eftersom GeoJSON är enkel textformat kan jag enkelt skriva ett program för att ersätta DN-taggarna till dem som OSM förstår ("landuse=forest" osv). Det går också att lägga till samtidigt lite brus till koordinater så att de importerade kurvorna inte ser ut som trappor. Och att filtrera linjer som man inte är intresserad i, till exempel vattenytor.

Mitt 140 kbyte TIFF-klipp producerade 2 Mbyte GeoJSON:en, så man bör jobba endast med lagom stora delområden. Nu är klockan mycket, så planerar jag fortsätta senare nästa vecka.


>Воскресенье, 31 марта 2019, 14:41 +03:00 от Grigory Rechistov via Talk-se <talk-se at openstreetmap.org>:
>
>>Det troliga är väl att det är något i JOSM (och pluginet "opendata") som är orsaken till problemen.
>
>Kanske är det .shp-filformatet själv som begränsar oss i det fallet. Jag vandrade genom wiki.openstreetmap.org och stötte på en anmärkning ( https://wiki.openstreetmap.org/wiki/Shapefiles ):
>
>NOTE: Shapefiles have many limitations, like cutting attribute names.
>
>Det ledde mig till  http://switchfromshapefile.org/ som tyder bland annat:
>
>10 Characters attribute names
>Attribute names are limited to 10 characters max. Longer names are usually automatically shortened. This leads to abbreviated and/or cryptic attribute names that are unintuitive to the recipient of the data.
>
>
>>Воскресенье, 31 марта 2019, 12:03 +03:00 от Christian Asker < christian.asker at gmail.com >:
>>
>>Hej. Jag kan som sagt se alla attributen i QGIS, så jag antar att dom 
>>har kommit med. Jag har använt liknande script för att skapa shapefiler 
>>som jag läser in i JOSM, tex för Terrängkartan och NVDB. Där har dock 
>>orginalfilerna varit vektorfiler, så det kanske är något i 
>>gdal_polygonize som är problemet.
>>
>>Öppnar jag shapefilen som skapas i första steget med JOSM så ser jag 
>>inga attribut där heller, trots att det finns attribut i shapefilen. Det 
>>är med andra ord redan där det blir problem. Det troliga är väl att det 
>>är något i JOSM (och pluginet "opendata") som är orsaken till problemen.
>>
>>
>>Jag ska prova att konvertera shapefilerna till OSM-format utanför JOSM 
>>för att se vad som händer. Återkommer ifall det löser problemen.
>>
>>Mvh Christian
>>
>>
>>
>>Den 2019-03-31 kl. 09:53, skrev Karl-Johan Karlsson:
>>> Jag använde Christian Asker skript för att skapa en shape-fil för ett 
>>> minimalt område efter det så konverterade jag shapefilen med 
>>> "ogr2osm.py -e 3006" till en osm-fil. Kollar man i osm-filen så ser 
>>> det ut som om alla attributen har tappats bort (vilket även ses när 
>>> man laddar in filen i JOSM). Jag kanske ska tillägga att jag inte har 
>>> verifierat att min genererade shape-fil har attributen, eftersom jag 
>>> inte har QGIS installerat.  Jag vet således inte i vilket steg 
>>> attributen har tappats bort (det kan ju vara handhavande fel från min 
>>> sida i de första stegen, då jag är nybörjare på GDAL).
>>
>>_______________________________________________
>>Talk-se mailing list
>>Talk-se at openstreetmap.org
>>https://lists.openstreetmap.org/listinfo/talk-se
>
>
>С наилучшими пожеланиями,
>Григорий Речистов.
>Med vänliga hälsningar,
>Grigory Rechistov
>With best regards,
>Grigory Rechistov
>_______________________________________________
>Talk-se mailing list
>Talk-se at openstreetmap.org
>https://lists.openstreetmap.org/listinfo/talk-se
 

С наилучшими пожеланиями,
Григорий Речистов.
Med vänliga hälsningar,
Grigory Rechistov
With best regards,
Grigory Rechistov
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.openstreetmap.org/pipermail/talk-se/attachments/20190401/29e341c6/attachment.html>