Saya mencoba menumpuk raster multiband. Mereka dibuat melalui alat Prapemrosesan ESA Sentinel 1 SNAP. Setiap file tif harus dilapis.
Saya memuat dua rasterstack dan mencoba menumpuknya:
rs1 <- raster::stack("example/rs1.tif")
rs2 <- raster::stack("example/rs2.tif")
rsstack <- stack(rs1,rs2)
Lalu saya mendapat pesan kesalahan berikut:
Kesalahan dalam perbandinganRaster(x) : tingkat yang berbeda
Luas rasterstack hampir sama:
> rs1
class : RasterBrick
dimensions : 2273, 2100, 4773300, 2 (nrow, ncol, ncell, nlayers)
resolution : 8.983153e-05, 8.983153e-05 (x, y)
extent : 8.183134, 8.37178, 48.49076, 48.69495 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +no_defs
source : memory
names : S1A_IW_GRDH_1SDV_20180110T053421_20180110T053446_020088_0223E7_69E7_10.1, S1A_IW_GRDH_1SDV_20180110T053421_20180110T053446_020088_0223E7_69E7_10.2
min values : 1.729380e-07, 1.077101e-06
max values : 11.63158, 109.76797
> rs2
class : RasterBrick
dimensions : 2273, 2100, 4773300, 2 (nrow, ncol, ncell, nlayers)
resolution : 8.983153e-05, 8.983153e-05 (x, y)
extent : 8.183171, 8.371817, 48.49071, 48.6949 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +no_defs
source : memory
names : S1A_IW_GRDH_1SDV_20180106T171534_20180106T171559_020037_02223E_CE2A_10.1, S1A_IW_GRDH_1SDV_20180106T171534_20180106T171559_020037_02223E_CE2A_10.2
min values : 8.244981e-08, 5.691331e-06
max values : 6.012002, 64.965996
Apa yang dapat saya lakukan untuk menggabungkan keduanya? Bagaimana cara menyesuaikan sejauh mana yang lain?
Apa yang saya coba: Saya juga punya bidang minat. JADI saya mencoba memotong kedua tumpukan ke area yang diminati dan mencoba menumpuknya lagi:
shp <- readOGR(dsn=path.expand(example/area.shp)))
shp <- sp::spTransform(shp, CRS(proj4string(rs[[1]])))
rs1 <- raster::crop(rs1,shp)
rs2 <- raster::crop(rs2,shp)
rsstack <- stack(rsstack,r2)
Kesalahan dalam bandingkanRaster(x) : tingkat yang berbeda
shp
class : SpatialPolygonsDataFrame
features : 16
extent : 8.183144, 8.371817, 48.49075, 48.69491 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +no_defs
variables : 14
names : fid, Area, BoundLen, CentX, CentY, AreaIncI, BoundNotIn, CentXNotIn, CentYNotIn, PointInPol, PointInPo1, CompactRat, CompactRa1, rast
min values : 1, 1440.64063, 168.67161, 3442255.57469, 5377418.66407, 1440.64063, 168.67161, 3442255.57469, 5377418.66407, 3442259.72286, 5377409.4208, 1.14669, 1.14669, 1
max values : 16, 76089100.06641, 89693.52095, 3451427.74745, 5393682.39749, 76858585.26953, 76662.07157, 3451427.74745, 5393682.39749, 3451416.96957, 5393687.75063, 2.57399, 2.90065, 1
Apakah ada yang punya solusi?