buffer distance

buffer distance example

library(GSIF)
## GSIF version 0.5-4 (2017-04-25)
## URL: http://gsif.r-forge.r-project.org/
library(sp)
library(ranger)
library(raster)
demo(meuse)
## 
## 
##  demo(meuse)
##  ---- ~~~~~
## 
## > require(sp)
## 
## > crs = CRS("+init=epsg:28992")
## 
## > data("meuse")
## 
## > coordinates(meuse) <- ~x+y
## 
## > proj4string(meuse) <- crs
## 
## > data("meuse.grid")
## 
## > coordinates(meuse.grid) <- ~x+y
## 
## > gridded(meuse.grid) <- TRUE
## 
## > proj4string(meuse.grid) <- crs
## 
## > data("meuse.riv")
## 
## > meuse.riv <- SpatialPolygons(list(Polygons(list(Polygon(meuse.riv)),"meuse.riv")))
## 
## > proj4string(meuse.riv) <- crs
## 
## > data("meuse.area")
## 
## > meuse.area = SpatialPolygons(list(Polygons(list(Polygon(meuse.area)), "area")))
## 
## > proj4string(meuse.area) <- crs
test=buffer.dist(meuse[1,'om'],meuse.grid[1,],as.factor(1))
test@data
##       layer
## 69 169.7056
test=buffer.dist(meuse["om"], meuse.grid[1],
                 classes=as.factor(1:nrow(meuse)))
head(test@data)
##       layer.1  layer.2  layer.3  layer.4  layer.5  layer.6  layer.7
## 69  169.70563 256.1250 200.0000 268.3282 417.6123 520.0000 360.0000
## 146 113.13708 200.0000 164.9242 256.1250 393.9543 501.1986 322.4903
## 147 144.22205 226.2742 160.0000 233.2381 379.4733 483.3218 320.0000
## 148 178.88544 256.1250 164.9242 215.4066 368.7818 468.1880 322.4903
## 223  56.56854 144.2221 144.2221 256.1250 377.3592 488.2622 291.2044
## 224  89.44272 169.7056 126.4911 226.2742 357.7709 466.4762 282.8427
##      layer.8  layer.9 layer.10 layer.11 layer.12 layer.13 layer.14
## 69  393.9543 533.6666 561.4268 640.0000 737.5636 512.2499 520.0000
## 146 341.7601 486.6210 526.1179 601.3319 690.5071 456.0702 468.1880
## 147 357.7709 494.7727 521.5362 600.0000 698.5700 481.6638 483.3218
## 148 377.3592 505.9644 520.0000 601.3319 708.8018 509.1169 501.1986
## 223 291.2044 441.8144 494.7727 565.6854 644.9806 400.0000 417.6123
## 224 304.6309 447.2136 486.6210 561.4268 651.1528 425.2058 430.8132
##     layer.15 layer.16 layer.17 layer.18  layer.19 layer.20 layer.21
## 69  582.4088 600.0000 754.7185 898.8882 1043.0724 1221.147 1162.755
## 146 533.6666 544.0588 699.7142 843.8009  987.9271 1166.190 1109.234
## 147 544.0588 568.5068 721.1103 865.3323 1009.5544 1187.266 1127.120
## 148 557.1355 594.6427 744.0430 888.1441 1032.2790 1209.297 1146.124
## 223 486.6210 488.2622 644.9806 788.9233  932.9523 1111.396 1056.030
## 224 494.7727 512.2499 665.7327 809.9383  954.1488 1132.078 1073.313
##     layer.22 layer.23 layer.24 layer.25 layer.26 layer.27  layer.28
## 69  1145.426 1181.863 800.9994 920.8692 960.0000 973.2420 1059.0562
## 146 1092.520 1129.248 760.0000 880.0000 920.8692 927.7931 1012.7191
## 147 1109.234 1145.426 761.0519 880.9086 920.0000 933.8094 1019.8039
## 148 1127.120 1162.755 764.1989 883.6289 920.8692 941.4882 1028.3968
## 223 1040.000 1077.033 721.1103 840.9518 883.6289 883.6289  967.4709
## 224 1056.030 1092.520 720.0000 840.0000 880.9086 888.1441  973.2420
##      layer.29 layer.30 layer.31 layer.32 layer.33 layer.34 layer.35
## 69  1038.4604 815.8431 1160.689 1280.625 1329.662 1381.014 1459.863
## 146  989.5454 785.8753 1120.000 1232.234 1280.625 1335.066 1414.214
## 147 1000.0000 776.6595 1120.714 1241.934 1291.201 1341.641 1420.422
## 148 1011.9289 769.4154 1122.854 1252.837 1302.920 1349.370 1427.725
## 223  941.4882 758.9466 1080.740 1184.567 1232.234 1289.961 1369.379
## 224  950.7891 747.2617 1080.000 1193.315 1241.934 1295.531 1374.627
##     layer.36 layer.37 layer.38 layer.39 layer.40 layer.41 layer.42
## 69  1480.000 1397.140 1360.000 1365.284 1509.437 1556.406 1671.407
## 146 1429.405 1344.619 1305.986 1310.267 1454.373 1503.064 1619.877
## 147 1442.221 1360.588 1324.840 1331.465 1475.669 1520.526 1634.136
## 148 1456.022 1377.534 1344.619 1353.514 1497.732 1538.831 1649.242
## 223 1379.275 1292.440 1252.198 1255.388 1399.428 1449.966 1568.694
## 224 1391.402 1307.823 1270.590 1276.244 1420.422 1466.970 1582.403
##     layer.43 layer.44 layer.45 layer.46 layer.47 layer.48 layer.49
## 69  1718.139 1880.425 1667.093 1826.034 1932.046 2075.765 2237.856
## 146 1665.173 1827.786 1612.452 1771.779 1878.297 2022.672 2185.040
## 147 1681.904 1843.909 1632.667 1791.089 1896.523 2039.608 2201.454
## 148 1699.412 1860.753 1653.602 1811.077 1915.411 2057.183 2218.468
## 223 1612.452 1775.387 1557.947 1717.673 1824.719 1969.772 2132.416
## 224 1628.742 1791.089 1577.847 1736.663 1842.607 1986.353 2148.488
##     layer.50 layer.51 layer.52 layer.53 layer.54 layer.55 layer.56
## 69  2380.252 2433.434 1697.528 1800.000 1905.361 2039.608 2150.349
## 146 2327.230 2380.252 1641.950 1744.018 1849.108 1983.129 2093.800
## 147 2344.014 2397.332 1664.692 1768.163 1874.460 2009.975 2121.509
## 148 2361.356 2414.953 1688.076 1792.875 1900.316 2037.253 2149.605
## 223 2274.379 2327.230 1586.443 1688.076 1792.875 1926.655 2037.253
## 224 2290.851 2344.014 1608.975 1712.075 1818.131 1953.458 2064.946
##     layer.57 layer.58 layer.59 layer.60 layer.61 layer.62 layer.63
## 69  2217.747 2080.385 1864.618 2292.597 2435.077 2581.472 2704.958
## 146 2161.481 2024.253 1808.867 2236.068 2378.571 2525.074 2648.773
## 147 2186.870 2049.000 1832.157 2263.449 2405.660 2551.235 2673.724
## 148 2212.691 2074.223 1856.017 2291.201 2433.105 2577.751 2699.037
## 223 2105.232 1968.146 1753.169 2179.541 2322.068 2468.684 2592.605
## 224 2130.540 1992.787 1776.288 2206.898 2349.127 2494.795 2617.480
##     layer.64 layer.65 layer.66 layer.67 layer.68 layer.69 layer.70
## 69  2848.719 2992.524 3136.367 3117.435 3167.081 3204.497 3464.333
## 146 2792.562 2936.392 3080.260 3062.025 3113.069 3151.634 3408.460
## 147 2817.375 2961.081 3104.835 3084.153 3131.773 3168.091 3432.084
## 148 2842.534 2986.101 3129.728 3106.638 3150.873 3184.965 3456.009
## 223 2736.421 2880.278 3024.169 3006.659 3059.150 3098.903 3352.611
## 224 2761.159 2904.892 3048.672 3028.663 3077.661 3115.124 3376.152
##     layer.71 layer.72 layer.73 layer.74 layer.75 layer.76 layer.77
## 69  3529.306 3576.143 3696.701 3650.534 3761.914 3749.773 3871.950
## 146 3473.557 3520.227 3640.879 3594.885 3706.211 3694.320 3816.595
## 147 3496.741 3544.009 3664.314 3617.734 3729.236 3716.558 3838.541
## 148 3520.227 3568.081 3688.035 3640.879 3752.546 3739.091 3860.777
## 223 3417.836 3464.333 3585.080 3539.265 3650.534 3638.901 3761.276
## 224 3440.930 3488.037 3608.435 3562.022 3673.472 3661.038 3783.120
##     layer.78 layer.79 layer.80 layer.81 layer.82 layer.83 layer.84
## 69  3880.206 3736.308 3560.899 3448.594 3555.503 3224.159 526.1179
## 146 3824.134 3680.217 3504.740 3392.462 3505.653 3168.091 481.6638
## 147 3848.532 3704.700 3529.533 3417.133 3517.272 3192.491 486.6210
## 148 3873.190 3729.450 3554.603 3442.092 3529.306 3217.204 494.7727
## 223 3768.076 3624.141 3448.594 3336.345 3456.009 3112.041 440.0000
## 224 3792.413 3648.561 3473.327 3360.952 3467.333 3136.367 441.8144
##     layer.85 layer.86 layer.87 layer.88 layer.89 layer.90 layer.91
## 69  2985.297 2841.127 3600.000 3744.009 3984.972 4138.792 3994.396
## 146 2929.573 2785.390 3544.009 3688.035 3929.173 4083.136 3939.137
## 147 2952.694 2808.558 3568.081 3712.034 3952.518 4105.996 3960.808
## 148 2976.172 2832.102 3592.436 3736.308 3976.129 4129.116 3982.763
## 223 2873.882 2729.689 3488.037 3632.079 3873.396 4027.505 3883.916
## 224 2896.895 2752.744 3512.036 3656.009 3896.665 4050.284 3905.483
##     layer.92 layer.93 layer.94 layer.95 layer.96 layer.97 layer.98
## 69  4192.374 4048.506 3847.077 4025.916 3617.071 3423.449 3506.565
## 146 4136.279 3992.393 3793.890 3971.901 3562.696 3368.798 3451.377
## 147 4160.769 4016.964 3810.932 3990.589 3582.178 3388.923 3472.866
## 148 4185.499 4041.782 3828.316 4009.589 3601.999 3409.399 3494.682
## 223 4080.196 3936.293 3740.802 3917.959 3508.390 3314.212 3396.233
## 224 4104.632 3960.808 3757.659 3936.496 3527.719 3334.187 3417.601
##     layer.99 layer.100 layer.101 layer.102 layer.103 layer.104 layer.105
## 69  3638.901  3889.884  3632.079  3492.850  3193.994  2794.280  2563.123
## 146 3583.518  3835.205  3578.827  3440.233  3141.974  2741.678  2510.777
## 147 3605.551  3855.386  3595.998  3456.241  3156.961  2757.680  2526.341
## 148 3627.892  3875.874  3613.530  3472.636  3172.381  2774.167  2542.440
## 223 3528.172  3780.582  3525.677  3387.743  3090.113  2689.238  2458.617
## 224 3550.099  3800.632  3542.654  3403.528  3104.835  2704.958  2473.863
##     layer.106 layer.107 layer.108 layer.109 layer.110 layer.111 layer.112
## 69   2390.983  2732.325  2909.295  3113.326  3115.124  3130.495  3011.710
## 146  2340.940  2682.089  2858.811  3062.025  3062.287  3076.882  2957.026
## 147  2352.871  2694.290  2871.376  3075.841  3078.701  3094.770  2977.247
## 148  2365.418  2707.028  2884.441  3090.113  3095.545  3113.069  2997.866
## 223  2291.201  2632.109  2808.558  3010.913  3009.585  3023.376  2902.413
## 224  2302.694  2643.936  2820.780  3024.434  3025.756  3041.052  2922.465
##     layer.113 layer.114 layer.115 layer.116 layer.117 layer.118 layer.119
## 69   2867.752  2718.235  2888.044  2712.047  2512.051  3463.409  2164.809
## 146  2813.112  2662.931  2832.102  2656.012  2456.013  3415.026  2110.924
## 147  2833.231  2684.772  2856.011  2680.298  2480.323  3424.617  2129.413
## 148  2853.769  2707.028  2880.278  2704.958  2505.035  3434.647  2148.488
## 223  2758.550  2607.681  2776.184  2600.000  2400.000  3366.898  2057.183
## 224  2778.489  2629.373  2800.000  2624.195  2424.211  3376.152  2075.379
##     layer.120 layer.121 layer.122 layer.123 layer.124 layer.125 layer.126
## 69   2346.742  2312.055  2355.929  2074.223  1708.801  1963.670  1938.659
## 146  2292.248  2256.014  2299.565  2017.920  1653.602  1908.298  1886.796
## 147  2312.055  2280.351  2325.511  2043.526  1675.231  1930.389  1901.578
## 148  2332.381  2305.125  2351.850  2069.589  1697.528  1953.049  1917.081
## 223  2237.856  2200.000  2243.212  1961.632  1598.499  1852.998  1835.211
## 224  2257.432  2224.230  2269.097  1987.159  1619.877  1874.887  1849.540
##     layer.127 layer.128 layer.129 layer.130 layer.131 layer.132 layer.133
## 69   2050.561  2206.173  715.5418  967.4709  2771.570  2059.126  912.1403
## 146  1999.200  2154.066  662.1178  913.0170  2720.294  2004.794  863.4813
## 147  2013.157  2169.239  680.0000  932.9523  2734.081  2024.253  873.6132
## 148  2027.807  2185.040  699.7142  954.1488  2748.381  2044.309  885.4377
## 223  1948.127  2102.189  609.2618  858.8364  2669.232  1950.590  815.8431
## 224  1961.632  2116.979  626.0990  878.1799  2682.685  1969.772  824.6211
##     layer.134 layer.135 layer.136 layer.137 layer.138 layer.139 layer.140
## 69   2163.331  2720.000  2828.427  2936.392  3761.276  3827.689  3917.754
## 146  2112.439  2665.933  2774.167  2883.609  3705.995  3772.108  3862.642
## 147  2125.653  2684.772  2793.421  2899.931  3727.734  3794.733  3883.916
## 148  2139.533  2704.071  2813.112  2916.710  3749.773  3817.643  3905.483
## 223  2061.844  2611.972  2720.000  2830.972  3650.753  3716.558  3807.571
## 224  2074.608  2630.589  2739.051  2847.034  3672.383  3739.091  3828.733
##     layer.141 layer.142 layer.143 layer.144 layer.145 layer.146 layer.147
## 69   4169.412  4062.118  4138.792  4280.748  4294.182  4473.030  4545.943
## 146  4115.337  4008.990  4086.270  4228.096  4240.943  4419.774  4491.770
## 147  4134.150  4025.916  4102.097  4244.149  4258.075  4436.936  4510.787
## 148  4153.264  4043.167  4118.252  4260.516  4275.512  4454.391  4530.077
## 223  4061.330  3955.958  4033.857  4175.548  4187.792  4366.601  4437.657
## 224  4080.000  3972.707  4049.494  4191.420  4204.759  4383.606  4456.546
##     layer.148 layer.149 layer.150 layer.151 layer.152 layer.153 layer.154
## 69   4596.346  4348.241  4232.635  4083.332  4019.950  4149.217  3774.652
## 146  4541.365  4294.182  4177.846  4028.896  3965.047  4093.702  3720.860
## 147  4562.280  4312.957  4198.285  4048.506  3985.775  4116.115  3739.091
## 148  4583.448  4332.020  4219.005  4068.415  4006.794  4138.792  3757.659
## 223  4486.424  4240.189  4123.106  3974.519  3910.192  4038.217  3667.151
## 224  4507.239  4258.826  4143.429  3993.995  3930.801  4060.542  3685.214
##     layer.155
## 69   3603.776
## 146  3558.202
## 147  3564.267
## 148  3570.770
## 223  3512.947
## 224  3518.636
plot(test)

classes 는 나중에 column의 개수와 관련된다. width는 buffer 계산에서 반경에 해당된다. 아래의 그림과 동일한 방식은 아니나 직관적으로 볼 수 있다.

classes=as.factor(1:nrow(meuse))

i=1
(width=areaSpatialGrid(meuse.grid))
## [1] 4964800
plot(buffer(meuse,width=5))
## Loading required namespace: rgeos

plot(buffer(meuse,width=50))

plot(buffer(meuse,width=500))

rasterize에 function의 기본값은 last이므로 아래는 meuse의 coords에 대한 인덱스를 value로 가지고있고 coordinates가 meuse.grid와 동일한 레스터를 생성한다.

즉, meuse.grid에서 meuse의 범주가 i번째에 해당하는 관측소에 값을 raster로 반환하고 나머지 meuse.grid에 해당하는 관측소에 값은 NA를 반환

rasterize(meuse[which(classes==levels(classes)[i]),]@coords,
raster(meuse.grid))
## class       : RasterLayer 
## dimensions  : 104, 78, 8112  (nrow, ncol, ncell)
## resolution  : 40, 40  (x, y)
## extent      : 178440, 181560, 329600, 333760  (xmin, xmax, ymin, ymax)
## coord. ref. : +init=epsg:28992 +proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +towgs84=565.4171,50.3319,465.5524,-0.398957,0.343988,-1.87740,4.0725 +units=m +no_defs 
## data source : in memory
## names       : layer 
## values      : 1, 1  (min, max)
meuse[which(classes==levels(classes)[i]),]@coords
##        x      y
## 1 181072 333611
coordinates(raster(meuse.grid))[300,]
##      x      y 
## 181060 333620

boundaries 는 외곽부분을 찾아주지만 del이 1개 관측소만 가지므로 del과 동일한 output을 나타냄.

rasterToPoints는 function조건을 만족하는 위경도와 값을 제공

del=rasterize(meuse[which(classes==levels(classes)[i]),]@coords,
                raster(meuse.grid))
(r=boundaries(del, classes=FALSE, type='inner') )
## class       : RasterLayer 
## dimensions  : 104, 78, 8112  (nrow, ncol, ncell)
## resolution  : 40, 40  (x, y)
## extent      : 178440, 181560, 329600, 333760  (xmin, xmax, ymin, ymax)
## coord. ref. : +init=epsg:28992 +proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +towgs84=565.4171,50.3319,465.5524,-0.398957,0.343988,-1.87740,4.0725 +units=m +no_defs 
## data source : in memory
## names       : layer 
## values      : 1, 1  (min, max)
(pts=rasterToPoints(r, fun=function(z){ z>0 } )[,1:2, drop=FALSE])
##           x      y
## [1,] 181060 333620

xy는 결측값에 해당하는 셀의 중심좌표를 추출한 값.

out <- raster(del)
x <- values(del)
  index <- which(is.na(x))
  x[] <- 0
  #결측에 해당하는 셀의 중심 좌표
  xy <- xyFromCell(out, i)

고도를 고려하기위해 유클리드 거리식 부분을 수정한 곳 원래는 Cpp로 구성되어 있음. 같은 결과 나오는지 확인하기 위해 아직 고도 반영을 하지 않음.

col=0;dis=NULL
while(col<ncol(xy)){
  col=col+1
  dis=cbind(dis,(xy[,col]-pts[,col])^2)
  
}
classes=as.factor(1:nrow(meuse))
s=list();s2=list()
i=1
for(i in 1:length(levels(classes))){
#반경
  width=areaSpatialGrid(meuse.grid)
#레벨에 해당하는 지역 추출
  del=rasterize(meuse[which(classes==levels(classes)[i]),]@coords,
                raster(meuse.grid))
  
  r <- boundaries(del, classes=FALSE, type='inner') 
  #해당 지역 위경도
  pts <- try(  rasterToPoints(r, fun=function(z){ z>0 } )[,1:2, drop=FALSE] )
  out <- raster(del)
  
  canProcessInMemory(out, 6)
  x <- values(del)
  index <- which(is.na(x))
  x[] <- 0
  #결측에 해당하는 셀의 중심 좌표
  xy <- xyFromCell(out, index)
col=0
dis=NULL
while(col<ncol(xy)){
  col=col+1
  dis=cbind(dis,(xy[,col]-pts[,col])^2)
  
}
x[index]=apply(dis,1,function(x){sqrt(sum(x))})
out <- setValues(out, x)
s2[[i]]=out
s[[i]]=raster::distance(del,width=width)
# range(x)
# range(values(raster::distance(del,width=width)))

# plot(out)
# plot(raster::distance(del,width=width))

# .Call('_raster_distanceToNearestPoint', xy[1,,drop=F], pts, longlat, 6378137.0, 1/298.257223563, PACKAGE='raster')
  
}


s <- s[sapply(s, function(x){!is.null(x)})]
s <- brick(s)
s <- as(s, "SpatialPixelsDataFrame")
s <- s[meuse.grid@grid.index,]

s2 <- s2[sapply(s2, function(x){!is.null(x)})]
s2 <- brick(s2)
s2 <- as(s2, "SpatialPixelsDataFrame")
s2 <- s2[meuse.grid@grid.index,]

head(s,1)
##     layer.1 layer.2 layer.3  layer.4  layer.5 layer.6 layer.7  layer.8
## 69 169.7056 256.125     200 268.3282 417.6123     520     360 393.9543
##     layer.9 layer.10 layer.11 layer.12 layer.13 layer.14 layer.15 layer.16
## 69 533.6666 561.4268      640 737.5636 512.2499      520 582.4088      600
##    layer.17 layer.18 layer.19 layer.20 layer.21 layer.22 layer.23 layer.24
## 69 754.7185 898.8882 1043.072 1221.147 1162.755 1145.426 1181.863 800.9994
##    layer.25 layer.26 layer.27 layer.28 layer.29 layer.30 layer.31 layer.32
## 69 920.8692      960  973.242 1059.056  1038.46 815.8431 1160.689 1280.625
##    layer.33 layer.34 layer.35 layer.36 layer.37 layer.38 layer.39 layer.40
## 69 1329.662 1381.014 1459.863     1480  1397.14     1360 1365.284 1509.437
##    layer.41 layer.42 layer.43 layer.44 layer.45 layer.46 layer.47 layer.48
## 69 1556.406 1671.407 1718.139 1880.425 1667.093 1826.034 1932.046 2075.765
##    layer.49 layer.50 layer.51 layer.52 layer.53 layer.54 layer.55 layer.56
## 69 2237.856 2380.252 2433.434 1697.528     1800 1905.361 2039.608 2150.349
##    layer.57 layer.58 layer.59 layer.60 layer.61 layer.62 layer.63 layer.64
## 69 2217.747 2080.385 1864.618 2292.597 2435.077 2581.472 2704.958 2848.719
##    layer.65 layer.66 layer.67 layer.68 layer.69 layer.70 layer.71 layer.72
## 69 2992.524 3136.367 3117.435 3167.081 3204.497 3464.333 3529.306 3576.143
##    layer.73 layer.74 layer.75 layer.76 layer.77 layer.78 layer.79 layer.80
## 69 3696.701 3650.534 3761.914 3749.773  3871.95 3880.206 3736.308 3560.899
##    layer.81 layer.82 layer.83 layer.84 layer.85 layer.86 layer.87 layer.88
## 69 3448.594 3555.503 3224.159 526.1179 2985.297 2841.127     3600 3744.009
##    layer.89 layer.90 layer.91 layer.92 layer.93 layer.94 layer.95 layer.96
## 69 3984.972 4138.792 3994.396 4192.374 4048.506 3847.077 4025.916 3617.071
##    layer.97 layer.98 layer.99 layer.100 layer.101 layer.102 layer.103
## 69 3423.449 3506.565 3638.901  3889.884  3632.079   3492.85  3193.994
##    layer.104 layer.105 layer.106 layer.107 layer.108 layer.109 layer.110
## 69   2794.28  2563.123  2390.983  2732.325  2909.295  3113.326  3115.124
##    layer.111 layer.112 layer.113 layer.114 layer.115 layer.116 layer.117
## 69  3130.495   3011.71  2867.752  2718.235  2888.044  2712.047  2512.051
##    layer.118 layer.119 layer.120 layer.121 layer.122 layer.123 layer.124
## 69  3463.409  2164.809  2346.742  2312.055  2355.929  2074.223  1708.801
##    layer.125 layer.126 layer.127 layer.128 layer.129 layer.130 layer.131
## 69   1963.67  1938.659  2050.561  2206.173  715.5418  967.4709   2771.57
##    layer.132 layer.133 layer.134 layer.135 layer.136 layer.137 layer.138
## 69  2059.126  912.1403  2163.331      2720  2828.427  2936.392  3761.276
##    layer.139 layer.140 layer.141 layer.142 layer.143 layer.144 layer.145
## 69  3827.689  3917.754  4169.412  4062.118  4138.792  4280.748  4294.182
##    layer.146 layer.147 layer.148 layer.149 layer.150 layer.151 layer.152
## 69   4473.03  4545.943  4596.346  4348.241  4232.635  4083.332   4019.95
##    layer.153 layer.154 layer.155
## 69  4149.217  3774.652  3603.776
head(s2,1)
##     layer.1 layer.2 layer.3  layer.4  layer.5 layer.6 layer.7  layer.8
## 69 169.7056 256.125     200 268.3282 417.6123     520     360 393.9543
##     layer.9 layer.10 layer.11 layer.12 layer.13 layer.14 layer.15 layer.16
## 69 533.6666 561.4268      640 737.5636 512.2499      520 582.4088      600
##    layer.17 layer.18 layer.19 layer.20 layer.21 layer.22 layer.23 layer.24
## 69 754.7185 898.8882 1043.072 1221.147 1162.755 1145.426 1181.863 800.9994
##    layer.25 layer.26 layer.27 layer.28 layer.29 layer.30 layer.31 layer.32
## 69 920.8692      960  973.242 1059.056  1038.46 815.8431 1160.689 1280.625
##    layer.33 layer.34 layer.35 layer.36 layer.37 layer.38 layer.39 layer.40
## 69 1329.662 1381.014 1459.863     1480  1397.14     1360 1365.284 1509.437
##    layer.41 layer.42 layer.43 layer.44 layer.45 layer.46 layer.47 layer.48
## 69 1556.406 1671.407 1718.139 1880.425 1667.093 1826.034 1932.046 2075.765
##    layer.49 layer.50 layer.51 layer.52 layer.53 layer.54 layer.55 layer.56
## 69 2237.856 2380.252 2433.434 1697.528     1800 1905.361 2039.608 2150.349
##    layer.57 layer.58 layer.59 layer.60 layer.61 layer.62 layer.63 layer.64
## 69 2217.747 2080.385 1864.618 2292.597 2435.077 2581.472 2704.958 2848.719
##    layer.65 layer.66 layer.67 layer.68 layer.69 layer.70 layer.71 layer.72
## 69 2992.524 3136.367 3117.435 3167.081 3204.497 3464.333 3529.306 3576.143
##    layer.73 layer.74 layer.75 layer.76 layer.77 layer.78 layer.79 layer.80
## 69 3696.701 3650.534 3761.914 3749.773  3871.95 3880.206 3736.308 3560.899
##    layer.81 layer.82 layer.83 layer.84 layer.85 layer.86 layer.87 layer.88
## 69 3448.594 3555.503 3224.159 526.1179 2985.297 2841.127     3600 3744.009
##    layer.89 layer.90 layer.91 layer.92 layer.93 layer.94 layer.95 layer.96
## 69 3984.972 4138.792 3994.396 4192.374 4048.506 3847.077 4025.916 3617.071
##    layer.97 layer.98 layer.99 layer.100 layer.101 layer.102 layer.103
## 69 3423.449 3506.565 3638.901  3889.884  3632.079   3492.85  3193.994
##    layer.104 layer.105 layer.106 layer.107 layer.108 layer.109 layer.110
## 69   2794.28  2563.123  2390.983  2732.325  2909.295  3113.326  3115.124
##    layer.111 layer.112 layer.113 layer.114 layer.115 layer.116 layer.117
## 69  3130.495   3011.71  2867.752  2718.235  2888.044  2712.047  2512.051
##    layer.118 layer.119 layer.120 layer.121 layer.122 layer.123 layer.124
## 69  3463.409  2164.809  2346.742  2312.055  2355.929  2074.223  1708.801
##    layer.125 layer.126 layer.127 layer.128 layer.129 layer.130 layer.131
## 69   1963.67  1938.659  2050.561  2206.173  715.5418  967.4709   2771.57
##    layer.132 layer.133 layer.134 layer.135 layer.136 layer.137 layer.138
## 69  2059.126  912.1403  2163.331      2720  2828.427  2936.392  3761.276
##    layer.139 layer.140 layer.141 layer.142 layer.143 layer.144 layer.145
## 69  3827.689  3917.754  4169.412  4062.118  4138.792  4280.748  4294.182
##    layer.146 layer.147 layer.148 layer.149 layer.150 layer.151 layer.152
## 69   4473.03  4545.943  4596.346  4348.241  4232.635  4083.332   4019.95
##    layer.153 layer.154 layer.155
## 69  4149.217  3774.652  3603.776
head(coordinates(s))
##           x      y
## [1,] 181180 333740
## [2,] 181140 333700
## [3,] 181180 333700
## [4,] 181220 333700
## [5,] 181100 333660
## [6,] 181140 333660
head(coordinates(s2))
##           x      y
## [1,] 181180 333740
## [2,] 181140 333700
## [3,] 181180 333700
## [4,] 181220 333700
## [5,] 181100 333660
## [6,] 181140 333660
plot(s)

plot(s2)