I want to pan sharpen my Landsat8 scene (according to FIHS algorithm, see here).
According to the paper the following equation should do the pansharpening:
MS_high <- MS_low + (pan - I)
"MS_high" is the pansharpened MS band and "MS_low" the original MS band.
The I image is intended to be a simulated image of the Pan band at the spatial resolution of the MS bands. Spatial details from Pan are obtained by subtracting the lower resolution I band from the higher resolution Pan band.
When doing the substraction (Pan - I) the following output occurs:
Error in compareRaster(e1, e2, extent = FALSE, rowcol = FALSE, crs = TRUE, : different resolution
Does anyone used the FIHS pansharpening algorithm in R so far?
Any suggestions are welcome.