I am writing a script for the processing of different hyperspectral datacubes. The datacubes contain flight lines that I've recorded with a UAV. In total there are 10 flight lines that need to be processed and mosaiced. The files are in .bsq format and contain 101 bands:
class : RasterBrick dimensions : 1293, 749, 968457, 101 (nrow, ncol, ncell, nlayers) resolution : 0.25, 0.25 (x, y) extent : 651041.9, 651229.1, 5687300, 5687623 (xmin, xmax, ymin, ymax) coord. ref. : +proj=utm +zone=31 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs data source : \HyperspectralDatacube_FlightLine02_PixelSize=0.25m_DSM=Photogrammetric.bsq names : X450.00000.Nanometers, X455.00000.Nanometers, X460.00000.Nanometers, X465.00000.Nanometers, X470.00000.Nanometers, ... I load all 10 hyperspectral rasters using lapply and brick to load all 101 layers:
# Load hyperspectral datacubes files
class : RasterBrick dimensions : 1293, 749, 968457, 101 (nrow, ncol, ncell, nlayers) resolution : 0.25, 0.25 (x, y) extent : 651041.9, 651229.1, 5687300, 5687623 (xmin, xmax, ymin, ymax) coord. ref. : +proj=utm +zone=31 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs data source : \HyperspectralDatacube_FlightLine02_PixelSize=0.25m_DSM=Photogrammetric.bsq names : X450.00000.Nanometers, X455.00000.Nanometers, X460.00000.Nanometers, X465.00000.Nanometers, X470.00000.Nanometers, ... I load all 10 hyperspectral rasters using lapply and brick to load all 101 layers:
# Load hyperspectral datacubes files