SF, a inne biblioteki GIS
Wiecie już jak korzystać z nowej biblioteki sf. Jest to nowa biblioteka i jej obiekty nie współpracują np. z biblioteką raster. Pokażmy to na przykładzie. Wczytujemy np. czworobok zapisany w pliku SHP, którym będziemy chcieli wyciąć fragment rastra:
|
1 2 3 4 5 6 7 8 9 10 11 |
library(sf) shp <- st_read("d:/czworobok.shp") > shp Simple feature collection with 1 feature and 1 field geometry type: POLYGON dimension: XY bbox: xmin: 663753.4 ymin: 499031.9 xmax: 667507.3 ymax: 501790.8 epsg (SRID): NA proj4string: +proj=tmerc +lat_0=0 +lon_0=19 +k=0.9993 +x_0=500000 +y_0=-5300000 +ellps=GRS80 +units=m +no_defs id geometry 1 NA POLYGON ((663866.3 501747.8... |
Wczytaliśmy prosty poligon. Do działania na rastrach konieczna jest biblioteka raster. Wczytujemy bibliotekę oraz obraz pokrywający się z wczytanym wcześniej poligonem. W naszym przypadku będą to dane wysokościowe dla jednego z województw:
|
1 2 3 4 5 6 7 8 9 10 11 |
library(raster) r <- raster("d:/grid.tif") > r class : RasterLayer dimensions : 2747, 2638, 7246586 (nrow, ncol, ncell) resolution : 100, 100 (x, y) extent : 517650, 781450, 352550, 627250 (xmin, xmax, ymin, ymax) coord. ref. : NA data source : d:\grid.tif names : grid values : 51.64, 407.16 (min, max) |
Spróbujmy wyciąć fragment rastera r funkcją crop z biblioteki raster przy wykorzystaniu warstwy shp:
|
1 |
> fragment <- crop(r,shp) |
Otrzymujemy błąd, ponieważ funkcja crop nie obsługuje obiektów klasy sf:
|
1 |
Error in methods::validObject(y) : invalid class “Extent” object: TRUE |
Jak sobie z tym poradzić? Konwertując obiekt sf do znanego i wykorzystywanego przez nas sp. A jak? Zamieniając klasę sf na SpatialPolygonsDataFrame używając funkcji as rzutującej jeden obiekt w drugi:
|
1 2 3 4 5 6 7 8 9 |
shp_sp <- as(shp,"Spatial") > shp_sp class : SpatialPolygonsDataFrame features : 1 extent : 663753.4, 667507.3, 499031.9, 501790.8 (xmin, xmax, ymin, ymax) coord. ref. : +proj=tmerc +lat_0=0 +lon_0=19 +k=0.9993 +x_0=500000 +y_0=-5300000 +ellps=GRS80 +units=m +no_defs variables : 1 names : id value : Inf |
Sprawdzamy czy teraz działa wycinanie:
|
1 2 3 4 5 6 7 8 9 10 |
fragment <- crop(r,shp_sp) > fragment class : RasterLayer dimensions : 27, 38, 1026 (nrow, ncol, ncell) resolution : 100, 100 (x, y) extent : 663750, 667550, 499050, 501750 (xmin, xmax, ymin, ymax) coord. ref. : NA data source : in memory names : grid values : 105.57, 118.73 (min, max) |
Biblioteka sf jest jeszcze “młoda” i jej klasy obiektów nie zostały zaimplementowane w innych przestrzennych bibliotekach R. Z czasem na pewno się to zmieni, ale na razie musimy konwertować jej obiekty na inne bardziej standardowe.
