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

Christian Asker christian.asker at gmail.com
Mon Apr 1 04:31:18 UTC 2019


Det låter lovande, bra jobbat!

Mvh Christian

Grigory Rechistov via Talk-se <talk-se at openstreetmap.org> skrev: (1 april 2019 00:17:11 CEST)
>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

-- 
Skickat från min Android-enhet med K-9 Mail. Ursäkta min fåordighet.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.openstreetmap.org/pipermail/talk-se/attachments/20190401/cd50c573/attachment-0001.html>


More information about the Talk-se mailing list