Combining hyperspectral image bands R, G, B and NIR

Good morning

I am working with hyperspectral images, I have 8 bands in TIF format, what happens is that I need to combine (concatenate) four of these bands, the band R, G, B and NIR to check the state of an atmospheric correction that I am performing but not I know how to do it, I know how to join the bands R, G and B but I do not know how to combine them with the NIR band.

Answer 1

[Subhadeep Koley]

% Read the TIF file
cube = imread('yourFile.tif'); % cube is your 8 band TIF file

% Extract R, G, B, and NIR band according to the band index. For example if
% R, G, B, and NIR bands corresponds to 1, 2, 3, and 4 band numbers
R = cube(:, :, 1);
G = cube(:, :, 1);
B = cube(:, :, 1);
NIR = cube(:, :, 1);

% Concatenate the bands
catImage = cat(3, R, G, B, NIR);