diff --git a/notebooks/05_cmw.Rmd b/notebooks/05_cmw.Rmd index 9058df2..681751b 100644 --- a/notebooks/05_cmw.Rmd +++ b/notebooks/05_cmw.Rmd @@ -318,7 +318,8 @@ plot(nDSM, main = "Wysokość drzew [m]", uwzględnisz: * wizualizacje NMT i NMPT -* statystyki opisowe wysokości terenu +* statystyki opisowe wysokości terenu uwzględniając wartość minimalną, +maksymalną, średnią oraz odchylenie standardowe * profil wysokościowy * wizualizacje wysokości obiektów (znormalizowany NMPT) diff --git a/notebooks/05_cmw.html b/notebooks/05_cmw.html index 1748021..15ce0ce 100644 --- a/notebooks/05_cmw.html +++ b/notebooks/05_cmw.html @@ -1591,17 +1591,17 @@

Pozyskanie danych

View().

# wyświetl 10 pierwszych wierszy i 6 pierwszych kolumn
 dane[1:10, 1:6]
-
##             sheetID year product              format resolution avgElevErr
-## 1    N-33-130-D-b-1 2007     DTM      Intergraph TTN       <NA>       1.50
-## 2    N-33-130-D-b-1 2007     DTM            ESRI TIN       <NA>       1.50
-## 3       6.179.11.14 2021     DTM ARC/INFO ASCII GRID      0.5 m       0.10
-## 4       6.179.11.14 2021     DSM ARC/INFO ASCII GRID      0.5 m       0.10
-## 5    N-33-130-D-b-1 2011     DTM           ASCII TBD     25.0 m       0.90
-## 6    N-33-130-D-b-1 2007     DTM           ASCII TBD     40.0 m       1.50
-## 7    N-33-130-D-b-1 2014     DTM           ASCII TBD     25.0 m       0.90
-## 8  N-33-130-D-b-1-1 2021     DTM ARC/INFO ASCII GRID      5.0 m       0.50
-## 9  N-33-130-D-b-1-1 2021     DTM ARC/INFO ASCII GRID      1.0 m       0.15
-## 10 N-33-130-D-b-1-1 2020     DTM ARC/INFO ASCII GRID      5.0 m       0.50
+
##             sheetID year    product              format resolution avgElevErr
+## 1  N-33-130-D-b-1-1 2019        DTM      ASCII XYZ GRID      5.0 m       0.60
+## 2       6.179.11.14 2021        DSM ARC/INFO ASCII GRID      0.5 m       0.10
+## 3     6.179.11.14.1 2021 PointCloud                 LAZ    24 p/m2       0.10
+## 4  N-33-130-D-b-1-1 2012        DSM      ASCII XYZ GRID      0.5 m       0.10
+## 5  N-33-130-D-b-1-1 2022        DTM ARC/INFO ASCII GRID      1.0 m       0.20
+## 6  N-33-130-D-b-1-1 2020        DTM ARC/INFO ASCII GRID      5.0 m       0.50
+## 7  N-33-130-D-b-1-1 2021        DTM ARC/INFO ASCII GRID      1.0 m       0.15
+## 8  N-33-130-D-b-1-1 2019        DTM ARC/INFO ASCII GRID      1.0 m       0.10
+## 9  N-33-130-D-b-1-1 2023        DTM ARC/INFO ASCII GRID     1.00 m       0.50
+## 10 N-33-130-D-b-1-1 2012        DTM      ASCII XYZ GRID      1.0 m       0.10

Jak możemy zauważyć powyższe metadane opisują produkty o różnych formatach, aktualności, rozdzielczości oraz dokładności. Do naszej analizy potrzebujemy numerycznego modelu terenu (DTM) i numerycznego @@ -1619,8 +1619,8 @@

Pozyskanie danych

dane_sel = rbind(DTM_sel, DSM_sel) dane_sel[, 1:6] 
##             sheetID year product              format resolution avgElevErr
-## 13 N-33-130-D-b-1-1 2019     DTM ARC/INFO ASCII GRID      1.0 m        0.1
-## 11 N-33-130-D-b-1-1 2019     DSM ARC/INFO ASCII GRID      0.5 m        0.1
+## 8 N-33-130-D-b-1-1 2019 DTM ARC/INFO ASCII GRID 1.0 m 0.1 +## 15 N-33-130-D-b-1-1 2019 DSM ARC/INFO ASCII GRID 0.5 m 0.1

Wykorzystajmy funkcję tile_download() do pobrania tych dwóch produktów.

options(timeout = 600)
@@ -1760,7 +1760,7 @@ 
Przetworzenie danych

 lapp().

 
plot(nDSM, main = "Wysokość drzew [m]",
      col = hcl.colors(9, palette = "Greens", rev = TRUE))

-

+

 
 
 

@@ -1769,7 +1769,8 @@ 
Zadanie

 obszaru, w której uwzględnisz:

 
    
 
  • wizualizacje NMT i NMPT
    • 
-
  • statystyki opisowe wysokości terenu
    • 
+
  • statystyki opisowe wysokości terenu uwzględniając wartość minimalną,
+maksymalną, średnią oraz odchylenie standardowe
    • 
 
  • profil wysokościowy
    • 
 
  • wizualizacje wysokości obiektów (znormalizowany NMPT)
    •