Hej. Jag tror som sagt att man kan fixa till detta med QGIS Raster
Calculator. Annars kan man skriva ett python-script som använder GDAL
för att fixa det manuellt, men jag tror att det blir mer jobb. Sen finns
det säkert fler sätt som jag inte känner till... :)
Jag håller med om att de pixlar som kommer att få samma taggar som
vektorer bör slås ihop, helst innan man konverterar från raster till
vektorer. Det borde ge mer enhetliga områden. Men frågan är bara ifall
det är så att taggarna borde vara olika, fast jag inte hittade det. Just
de olika typerna av lövskog hade jag svårt att hitta olika taggning för
i OSM.
Bra jobbat förresten! Jag ber om ursäkt för att jag inte varit så
delaktig i detta den senaste veckan, men jag har varit lite (för) upptagen.
Mvh Christian
Den 2019-04-06 kl. 00:30, skrev Grigory Rechistov via Talk-se:
Hej Peter,
Jo, jag har också märkt det och tycker att det är konstigt. Det
härstammar troligen från att olika rasterpixelsvärden motsvarar till
samma etikettsuppsättning. Till exempel, GeoTiffs värden 115, 116 och
117 blir till samma taggar:
## 115 Triviallövskog utanför våtmark
## 116 Ädellövskog utanför våtmark
## 117 Triviallövskog med ädellövinslag utanför våtmark
mapper['115']={"source": "NV NMD2018", "landuse": "forest",
"leaf_type": "broadleaved", "leaf_cycle": "deciduous"}
mapper['116']={"source": "NV NMD2018", "landuse": "forest",
"leaf_type": "broadleaved", "leaf_cycle": "deciduous"}
mapper['117']={"source": "NV NMD2018", "landuse": "forest",
"leaf_type": "broadleaved", "leaf_cycle": "deciduous"}
Det betyder att fler intilliggande rasterareor med olika DN-taggar kan
förvandlas till vektorvägar med identiska taggar. Det är nog korrekt
men redundant och leder till fler objekt än verkligen behövs.
Jag kan föreställa mig två möjliga lösningar:
1. Slå samman vägar i mitt skript. Men just nu har jag ingen aning hur
algoritmen för sådan transformation ska se ut.
2. Ändra värden i rasterbilden för att ha en-till-en överensstämmelse
mellan GeoTiff-värden och OSM-taggar. Det vill säga att ersätta 115 ->
117, 116 -> 117, och 117 -> 117. Hittills har jag inte lärt mig göra
det. Man tipsar att QGIS rasterkalkylator kan göra det, och just nyss
lärde jag att man kan ha ett slags "if-then-else"-yttryck i den:
https://gis.stackexchange.com/questions/141305/conditional-calculation-qgis-raster-calculator-or-grass-r-mapcalculator/141307
. Jag ska försöka tillämpa det imorgon.
Пятница, 5 апреля 2019, 23:18 +03:00 от Peter Svensson
<[email protected]>:
Jag har tittat på https://atakua.org/p/nmd/vinon-3.osm.gz och har
en fråga:
Om man tittar på områdena kring nod -1424890 (59.18809091653,
15.73068800312) så har området öster om (way -1434062), samt
väster om (way -1434100) denna nod exakt samma uppsättning tags:
"landuse"="forest"
"leaf_type"="broadleaved"
"source"="NV NMD2018"
"leaf_cycle"="deciduous"
Borde inte detta betyda att dessa områden med bör slås ihop till
ett enda större område?
mvh
On Fri, Apr 5, 2019 at 9:39 PM Grigory Rechistov via Talk-se
<[email protected]
<//e.mail.ru/compose/?mailto=mailto%3atalk%[email protected]>>
wrote:
Och samtidigt fortsätter jag förbättra mina skript. Skriptet
finns här
https://gist.github.com/grigory-rechistov/39c7e329cb1f9b42a97ca5960377173d
och det tar in en OSM fil som är direkt konvertering av en
GeoJSON-fil. Den sistnämnda filen kommer med de ursprungliga
"DN"-taggarna. Sedan förbättrar skriptet innehållet:
1. Ersätter "DN=nummer" taggar med "landuse=*" taggar enligt
schemat beskrivet här:
https://wiki.openstreetmap.org/wiki/WikiProject_Sweden/NMD_2018_Import_Plan#Tagging_Plans
.
2. Alla onödiga polygoner (vatten, byggnader osv) tappas bort
3. Alla duplicerade noder slås samman
4. Alla självkorsningar rättas till.
Den resulterande OSM-data har nu noll fel och betydligt mindre
varningar efter valideringen.
Här är min processen med dataimporteringen på Vinön
https://osm.org/go/0ZoM9_P?m= . Jämfört med Gränsön anser jag
att den har förbättrats, på stort sätt.
1. Den ursprungliga GeoJSON
https://atakua.org/p/nmd/vinon.geojson.gz och motsvarande OSM:
https://atakua.org/p/nmd/vinon-0.osm.gz har runt 10000 problem.
2. Efter att mitt skript körs finns det bara 137 varningar av
samma typ på ca 16000 nya noder. Se bilden:
https://atakua.org/p/nmd/vinon-1.png och datafilen:
https://atakua.org/p/nmd/vinon-1.osm.gz
3. Efter att jag manuellt söker alla mindre polygoner (färre
än 10 noder) och slänger dem kvarstår bara 18 varningar att
rätta till. Jag vill nämligen förbättra skriptet så att det
automatiskt raderar mindre polygoner. Bilden:
https://atakua.org/p/nmd/vinon-2.png och OSM-filen:
https://atakua.org/p/nmd/vinon-2.osm.gz
4. Samtliga kvarstående problem är att två polygoner
sammanfaller, varav den ena är inre väg (utan taggar) i en
multipolygon och den andra bär taggar. Jag anser att orsaken
är faktiskt en inkorrekt import från GeoJSON. Det går att
rätta till problemet automatiskt i mitt skript, men jag hann
inte göra det ännu.
5. Hur som helst, noll problem kvar efter 10 minuter manuellt
arbete! Bilden: https://atakua.org/p/nmd/vinon-3.png och
OSM-filen:
https://atakua.org/p/nmd/vinon-3.osm.gz
Nu minns jag att man har märkt att raka linjer längs vägar ser
ut som zig-zag eller sågtänder:
https://atakua.org/p/nmd/zigzagf%c3%a4lt.png . Det går enkelt
att rätta det till manuellt.
Det "förenkla yta (Skift-Y)" verktyget i JOSM med max-error=20
(https://josm.openstreetmap.de/wiki/Help/Action/SimplifyWay)
förvandlar den fula saken till en helt rak sträcka:
https://atakua.org/p/nmdett enkelt verktyg för
detta./f%c3%b6renkla-yta.png
<https://atakua.org/p/nmd/f%c3%b6renkla-yta.png>
Пятница, 5 апреля 2019, 17:19 +03:00 от Grigory Rechistov
<[email protected]
<//e.mail.ru/compose/?mailto=mailto%[email protected]>>:
Hej!
> Engligt
https://wiki.openstreetmap.org/wiki/Import/Guidelines ska
man även dokumentera sådana kommande importeringar och
tillkännage dem. Finns det någon som är intresserad att
fylla i Import/Catalogue-sidan och skriva en plan?
Okej, så här är mitt utkast på en importeringsplan:
https://wiki.openstreetmap.org/wiki/WikiProject_Sweden/NMD_2018_Import_Plan.
All feedback är välkommen!
Jag tänker att skriva till [email protected]
<http://e.mail.ru/compose/?mailto=mailto%[email protected]>
tillkännage importeringen i några dagar.
Med vänliga hälsningar,
Grigory Rechistov
With best regards,
Grigory Rechistov
_______________________________________________
Talk-se mailing list
[email protected]
<//e.mail.ru/compose/?mailto=mailto%3atalk%[email protected]>
https://lists.openstreetmap.org/listinfo/talk-se
С наилучшими пожеланиями,
Григорий Речистов.
Med vänliga hälsningar,
Grigory Rechistov
With best regards,
Grigory Rechistov
_______________________________________________
Talk-se mailing list
[email protected]
https://lists.openstreetmap.org/listinfo/talk-se
_______________________________________________
Talk-se mailing list
[email protected]
https://lists.openstreetmap.org/listinfo/talk-se
