Making a new Colored "TIF" and "RGB" image out of Landsat-8 Geotiff image(11 separate bands)

Hi There,

I've downloaded a Landsat-8 image from USGS website, and now I have all the 11 bands(12 images and Metadata file), But all of them are gray-scale images, I want to:

1. Create an RGB Image and, 2. Create a new colored TIF image

out of proper bands,

To the best of my knowledge, color channels/Bands are as follows: Blue: Band_2 Green: Band_3 Redd: Band_4.

In order to make an RGB file I tried to combine them using the following code:

File_Name = File_Name_B2.TIF
[B2, ~] = geotiffread(File_Name);
File_Name = File_Name_B3.TIF
[B3, ~] = geotiffread(File_Name);
File_Name = File_Name_B4.TIF
[B4, ~] = geotiffread(File_Name);
RGB_Image = cat(3, B4, B3, B2);

but result image is still in gray-scale format.

Any help is appreciated.

[Image Analyst]

Then B2, B3, and B4 are all identical.

[Cyrus]

Yes, that's right, They are separate images.

Answer 1

[Chad Greene]

I wrote a function to do this for Landsat 8, but it's basically the same method you use. Check out L8read.

[Cyrus]

Hi Chad,

Thank you so much for the code, I used it in my code, as follows:

RGB_Image = L8read(Folder_Name);
figure(), imshow(RGB_Image, []), title('RGB Image')

But the result is still in gray-scale format. Do you have any suggestions about it?

Thanks

[Chad Greene]

Without an example it's hard to know exactly what's going on. Are you sure it's not just low light that makes it appear grayscale? Try

RGB_Image = L8read(Folder_name,'stretch');