ฉันจะซ้อนแรสเตอร์ที่มีขอบเขตเกือบเท่ากันใน R ได้อย่างไร

ฉันกำลังพยายามซ้อนแรสเตอร์หลายแบนด์ สร้างขึ้นผ่านเครื่องมือเตรียมการประมวลผล ESA Sentinel 1 SNAP ไฟล์ tif ทุกไฟล์จะต้องมีเลเยอร์

ฉันโหลดแรสเตอร์สแต็กสองอันแล้วลองซ้อนมัน:

rs1 <- raster::stack("example/rs1.tif")
rs2 <- raster::stack("example/rs2.tif")
rsstack <- stack(rs1,rs2)

จากนั้นฉันได้รับข้อความแสดงข้อผิดพลาดต่อไปนี้:

ข้อผิดพลาดใน comparisonRaster(x) : ขอบเขตที่แตกต่างกัน

ขอบเขตของแรสเตอร์สแต็กเกือบจะเหมือนกัน:

> 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

ฉันจะทำอย่างไรเพื่อรวมทั้งสองเข้าด้วยกัน ฉันจะปรับขอบเขตไปที่อื่นได้อย่างไร?

สิ่งที่ฉันลอง: ฉันมีเรื่องที่สนใจด้วย ดังนั้นฉันจึงพยายามครอบตัดทั้งสองสแต็กไปยังพื้นที่ที่สนใจแล้วลองซ้อนอีกครั้ง:

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)

ข้อผิดพลาดใน comparisonRaster(x) : ขอบเขตที่แตกต่างกัน

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

ไม่มีใครมีวิธีแก้ปัญหา?


r gis r-raster satellite-image extent
person Muesgen    schedule 13.08.2020    source แหล่งที่มา

คำตอบ (1)


arrow_upward
1
arrow_downward

ฉันเห็นว่าแรสเตอร์ของคุณมีระบบพิกัดเดียวกัน (crs : +proj=longlat +datum=WGS84 +no_defs) ขอบเขตเท่านั้นที่แตกต่างกัน ขั้นแรก คุณสามารถครอบตัดมัน จากนั้นจึงใช้ฟังก์ชัน resample จากแพ็คเกจ raster ได้

library(raster)

#Crop the raster
rs1 <- raster::crop(rs1,shp)
rs2 <- raster::crop(rs2,shp)

#Conversion of rasters into same extent
rs2_resampled <- resample(rs2, rs1, method='bilinear')

#Stack the rasters
rsstack <- stack(rs1,rs2_resampled)
person Bappa Das    schedule 13.08.2020