他のフォーマットへの変換 WKT、WKB、sp
よく知られたテキスト (wkt) との変換
シンプルフィーチャが表示されているのは、通常、Well Known Text (wkt) です。
st_point(c(0,1))
## POINT (0 1)
st_linestring(matrix(0:9,ncol=2,byrow=TRUE))
## LINESTRING (0 1, 2 3, 4 5, 6 7, 8 9)これらの Well Known Text 文字列は、st_as_text を使用して明示的に作成することができます。
x = st_linestring(matrix(0:9,ncol=2,byrow=TRUE))
str = st_as_text(x)
x
## LINESTRING (0 1, 2 3, 4 5, 6 7, 8 9)WKT からは st_as_sfc を使って戻すことができます。
st_as_sfc(str)
## Geometry set for 1 feature
## Geometry type: LINESTRING
## Dimension: XY
## Bounding box: xmin: 0 ymin: 1 xmax: 8 ymax: 9
## CRS: NA
## LINESTRING (0 1, 2 3, 4 5, 6 7, 8 9)well-known binary への変換と well-known binary からの変換
Well-known binary は st_as_binary によって、シンプルフィーチャから生成されます。
x = st_linestring(matrix(0:9,ncol=2,byrow=TRUE))
(x = st_as_binary(x))
## [1] 01 02 00 00 00 05 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 f0 3f
## [26] 00 00 00 00 00 00 00 40 00 00 00 00 00 00 08 40 00 00 00 00 00 00 10 40 00
## [51] 00 00 00 00 00 14 40 00 00 00 00 00 00 18 40 00 00 00 00 00 00 1c 40 00 00
## [76] 00 00 00 00 20 40 00 00 00 00 00 00 22 40
class(x)
## [1] "raw"st_as_binary が返すオブジェクトはクラス WKB で、生のベクトルを含むリストか、単一の生のベクトルのどちらかです。これらは rawToHex を用いて、16進数の文字ベクトルに変換することができます。
rawToHex(x)
## [1] "0102000000050000000000000000000000000000000000f03f000000000000004000000000000008400000000000001040000000000000144000000000000018400000000000001c4000000000000020400000000000002240"sf に戻すには、 st_as_sfc を使用します。
x = st_as_binary(st_sfc(st_point(0:1), st_point(5:6)))
st_as_sfc(x)
## Geometry set for 2 features
## Geometry type: POINT
## Dimension: XY
## Bounding box: xmin: 0 ymin: 1 xmax: 5 ymax: 6
## CRS: NA
## POINT (0 1)
## POINT (5 6)sp への変換と sp からの変換
パッケージ sp が保持する空間オブジェクトは、 st_as_sf と st_as_sfc によってそれぞれシンプルフィーチャオブジェクトとジオメトリに変換することができます。
methods(st_as_sf)
## [1] st_as_sf.data.frame* st_as_sf.lpp* st_as_sf.map*
## [4] st_as_sf.owin* st_as_sf.ppp* st_as_sf.ppplist*
## [7] st_as_sf.psp* st_as_sf.s2_geography* st_as_sf.sf*
## [10] st_as_sf.sfc* st_as_sf.Spatial* st_as_sf.SpatVector*
## see '?methods' for accessing help and source code
methods(st_as_sfc)
## [1] st_as_sfc.bbox* st_as_sfc.blob*
## [3] st_as_sfc.character* st_as_sfc.dimensions*
## [5] st_as_sfc.factor* st_as_sfc.list*
## [7] st_as_sfc.map* st_as_sfc.owin*
## [9] st_as_sfc.pq_geometry* st_as_sfc.psp*
## [11] st_as_sfc.raw* st_as_sfc.s2_geography*
## [13] st_as_sfc.sf* st_as_sfc.SpatialLines*
## [15] st_as_sfc.SpatialMultiPoints* st_as_sfc.SpatialPixels*
## [17] st_as_sfc.SpatialPoints* st_as_sfc.SpatialPolygons*
## [19] st_as_sfc.tess* st_as_sfc.WKB*
## see '?methods' for accessing help and source code例を挙げると、次のようになります。
library(sp)
data(meuse)
coordinates(meuse) = ~x+y
m.sf = st_as_sf(meuse)
opar = par(mar=rep(0,4))
plot(m.sf)
## Warning: plotting the first 10 out of 12 attributes; use max.plot = 12 to plot
## all
クラス sf または sfc のシンプルフィーチャオブジェクトから、対応する Spatial* オブジェクトへの変換は、 as メソッドを使用して、 Spatial に強制的に変換されます。
x = st_sfc(st_point(c(5,5)), st_point(c(6,9)), crs = 4326)
as(x, "Spatial")
## SpatialPoints:
## coords.x1 coords.x2
## [1,] 5 5
## [2,] 6 9
## Coordinate Reference System (CRS) arguments: +proj=longlat +datum=WGS84
## +no_defs