Geotorget produktstöd
Exemplet beskriver på hur du kan läsa in en ytmodellsfil (LAZ) och höjddata från Höjdmodell nedladdning och jämföra höjdvärdena. Skillnaden skrivs till en GeoTIFF.
Använd fmw-filer för att lösa uppgiften
För att lösa uppgiften finns det två fmw-filer du kan använda. Den första filen tar fram boundingboxen för det område som ska jämföras och sedan anropar den fil nummer två med max/min-koordinater som indata.
Så här ser innehållet i fil nummer ett ut:
Fil nummer ett: JamforDEMFranYtmodellMedWCSFranHojdmodellNedladdningWS1.fmw
Så här ser innehållet i fil nummer två ut:
Fil nummer två: JamforDEMFranYtmodellMedWCSFranHojdmodellNedladdningWS2.fmw
Anledningen till att det finns två filer är på grund av att vi inte kom på något bättre sätt att hämta ut WCS-höjder ur Höjdmodell nedladdning med boundingbox-koordinater från den ytmodellsfil som ska kontrolleras.
Det borde gå att fixa med till exempel en FeatureReader men det fungerade inte när vi testade det.
Beskrivning av innehållet i filerna
I fil nummer ett läses en ytmodellsfil (laz) och från den tas dess boundingbox ut.
Workspace-parametrarna i WorkspaceRunner-transformern i fil nummer ett ska vara de som behövs för att köra fil nummer två, se exempel i bilden nedan:
I fil nummer två skapas en raster DEM från ytmodellsfilen p.s.s. som har beskrivits i dokument KonverteraYtmodellsfilMedFME.docx. Dessutom läses höjddata från Höjddata nedladdning p.s.s. som du kan se exempel på i guiden Hämta höjddata med FME.
Höjdgriddet från WCS:en samplas sedan om för att få samma upplösning som DEM:en från ytmodellen. Sedan sätts NoData-värdet till -999 för båda dataströmmarna. Efter det beräknas skillnaden mellan höjddatavärdena genom att ta skillnaden mellan deras pixelvärden.
Exempel på utdata
På bilden nedan visas exempel på utdata. Ljusa delar = stora skillnader. I skillnadsbilden går det att se skogspartier och byggnader eftersom de skiljer sig mycket mot markhöjdmodellen från Höjddata nedladdning:
