File Exchange

image thumbnail

loadLandSat8(metaFi​leName,bandList)

version 1.0 (4.16 KB) by

loading LandSat images using the metadata

9 Downloads

Updated

View License

This file loads landSat 8 images. Instead of individual TIF images, you need to point it to the meta data.
The usage is:

LandSatData=loadLandSat8('LC80400372014241LGN00_MTL.txt');

Note that the input parameter is the txt file. the output, i.e. LandSatData would be a structure:

Band: {11x1 cell}
BandInfo: {11x1 cell}
BQA: [7841x7691 uint16]
BQAInfo: [1x1 struct]
MetaData: [1x1 struct]

Above commands will load all bands including BQA band. If you want to load only certain bands then issue:

LandSatData=loadLandSat8('LC80400372014241LGN00_MTL.txt',[2,3,4,12]);

In above example, only band 2, 3 and 4 are loaded. The rest would be empty. Band 12 refers to Band Quality Data or BQA. If not provided it won't be loaded.

You can make a color composite as follows:

RGB=(cat(3, LandSatData.Band{4}, LandSatData.Band{3}, LandSatData.Band{2}));
imshow(imadjust(RGB,stretchlim(RGB)),[]);

Comments and Ratings (2)

Hi. This is a great function. How can I plot the data in geographic coordinates once it is loaded?

شكرا

MATLAB Release
MATLAB 8.4 (R2014b)

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video