Landsat 8 Tools¶

The following tools can be used to easily convert individual Landsat 8 bands into multispectral composites for use in remote sensing.

Install Bassmap¶

Install Bassmap package to python notebook or conda environment. Make sure to install GDAL to your notebook or environment with:

notebook: !pip install bassmap gdal

conda env: pip install bassmap and conda install -c conda-forge gdal

In [ ]:

# !pip install bassmap
# !pip install gdal

# !pip install bassmap # !pip install gdal

Import bassmap module

In [1]:

import bassmap
from osgeo import gdal

import bassmap from osgeo import gdal

Define file path to Landsat 8 GeoTIFF bands

In [2]:

### If using local file directory

red_band = gdal.Open('C:/Users/bassn/temp_dir/LC09_L1TP_021035_20230419_20230420_02_T1_B4.TIF').ReadAsArray()
green_band = gdal.Open('C:/Users/bassn/temp_dir/LC09_L1TP_021035_20230419_20230420_02_T1_B3.TIF').ReadAsArray()
blue_band = gdal.Open('C:/Users/bassn/temp_dir/LC09_L1TP_021035_20230419_20230420_02_T1_B2.TIF').ReadAsArray()
nir_band = gdal.Open('C:/Users/bassn/temp_dir/LC09_L1TP_021035_20230419_20230420_02_T1_B5.TIF').ReadAsArray()
swir_band = gdal.Open('C:/Users/bassn/temp_dir/LC09_L1TP_021035_20230419_20230420_02_T1_B6.TIF').ReadAsArray()
swir2_band = gdal.Open('C:/Users/bassn/temp_dir/LC09_L1TP_021035_20230419_20230420_02_T1_B7.TIF').ReadAsArray()
georef = gdal.Open('C:/Users/bassn/temp_dir/LC09_L1TP_021035_20230419_20230420_02_T1_B4.TIF')

### If using local file directory red_band = gdal.Open('C:/Users/bassn/temp_dir/LC09_L1TP_021035_20230419_20230420_02_T1_B4.TIF').ReadAsArray() green_band = gdal.Open('C:/Users/bassn/temp_dir/LC09_L1TP_021035_20230419_20230420_02_T1_B3.TIF').ReadAsArray() blue_band = gdal.Open('C:/Users/bassn/temp_dir/LC09_L1TP_021035_20230419_20230420_02_T1_B2.TIF').ReadAsArray() nir_band = gdal.Open('C:/Users/bassn/temp_dir/LC09_L1TP_021035_20230419_20230420_02_T1_B5.TIF').ReadAsArray() swir_band = gdal.Open('C:/Users/bassn/temp_dir/LC09_L1TP_021035_20230419_20230420_02_T1_B6.TIF').ReadAsArray() swir2_band = gdal.Open('C:/Users/bassn/temp_dir/LC09_L1TP_021035_20230419_20230420_02_T1_B7.TIF').ReadAsArray() georef = gdal.Open('C:/Users/bassn/temp_dir/LC09_L1TP_021035_20230419_20230420_02_T1_B4.TIF')

If using Google Colab**

In [6]:

# Uncomment to use in Google Colab

# red_band = gdal.Open('/content/LC09_L1TP_021035_20230419_20230420_02_T1_B4.TIF').ReadAsArray()
# green_band = gdal.Open('/content/LC09_L1TP_021035_20230419_20230420_02_T1_B3.TIF').ReadAsArray()
# blue_band = gdal.Open('/content/LC09_L1TP_021035_20230419_20230420_02_T1_B2.TIF').ReadAsArray()
# nir_band = gdal.Open('/content/LC09_L1TP_021035_20230419_20230420_02_T1_B5.TIF').ReadAsArray()
# swir_band = gdal.Open('/content/LC09_L1TP_021035_20230419_20230420_02_T1_B6.TIF').ReadAsArray()
# swir2_band = gdal.Open('/content/LC09_L1TP_021035_20230419_20230420_02_T1_B7.TIF').ReadAsArray()
# georef = gdal.Open('/content/LC09_L1TP_021035_20230419_20230420_02_T1_B4.TIF')

# Uncomment to use in Google Colab # red_band = gdal.Open('/content/LC09_L1TP_021035_20230419_20230420_02_T1_B4.TIF').ReadAsArray() # green_band = gdal.Open('/content/LC09_L1TP_021035_20230419_20230420_02_T1_B3.TIF').ReadAsArray() # blue_band = gdal.Open('/content/LC09_L1TP_021035_20230419_20230420_02_T1_B2.TIF').ReadAsArray() # nir_band = gdal.Open('/content/LC09_L1TP_021035_20230419_20230420_02_T1_B5.TIF').ReadAsArray() # swir_band = gdal.Open('/content/LC09_L1TP_021035_20230419_20230420_02_T1_B6.TIF').ReadAsArray() # swir2_band = gdal.Open('/content/LC09_L1TP_021035_20230419_20230420_02_T1_B7.TIF').ReadAsArray() # georef = gdal.Open('/content/LC09_L1TP_021035_20230419_20230420_02_T1_B4.TIF')

Remote Sensing Applications¶

Using the following functions, you can create different spectral composites of Landsat 8 satellite data using local GeoTIFF files downloaded from USGS EarthExplorer.

Get True Color composite¶

In [3]:

from bassmap import get_truecolor
get_truecolor(red_band, green_band, blue_band, georef)

from bassmap import get_truecolor get_truecolor(red_band, green_band, blue_band, georef)

Out[3]:

'C:\\Users\\bassn\\AppData\\Local\\Temp\\tmp94oux5ua\\true_color_comp.tif'

Get Color Infrared (Buildings) composite¶

In [4]:

from bassmap import get_color_infrared
get_color_infrared(nir_band, red_band, green_band, georef)

from bassmap import get_color_infrared get_color_infrared(nir_band, red_band, green_band, georef)

Out[4]:

'C:\\Users\\bassn\\AppData\\Local\\Temp\\tmp94oux5ua\\color_infrared_comp.tif'

Get False Color (Vegetation) composite¶

In [5]:

from bassmap import get_false_color
get_false_color(swir2_band, swir_band, red_band, georef)

from bassmap import get_false_color get_false_color(swir2_band, swir_band, red_band, georef)

Out[5]:

'C:\\Users\\bassn\\AppData\\Local\\Temp\\tmp94oux5ua\\false_color_comp.tif'

Get Healthy Vegetation composite¶

In [6]:

from bassmap import get_health_veg
get_health_veg(nir_band, swir_band, blue_band, georef)

from bassmap import get_health_veg get_health_veg(nir_band, swir_band, blue_band, georef)

Out[6]:

'C:\\Users\\bassn\\AppData\\Local\\Temp\\tmp94oux5ua\\healthy_veg_comp.tif'

Get Normalized Difference Moisture Index (NDMI) composite¶

In [7]:

from bassmap import get_ndmi
get_ndmi(nir_band, swir_band, georef)

from bassmap import get_ndmi get_ndmi(nir_band, swir_band, georef)

Out[7]:

'C:\\Users\\bassn\\AppData\\Local\\Temp\\tmp94oux5ua\\ndmi_composite.tif'

Get Normalized Difference Vegetation Index (NDVI) composite¶

In [8]:

from bassmap import get_NDVI
get_NDVI(red_band, nir_band, georef)

from bassmap import get_NDVI get_NDVI(red_band, nir_band, georef)

Out[8]:

'C:\\Users\\bassn\\AppData\\Local\\Temp\\tmp94oux5ua\\ndvi_composite.tif'

---

Last update: 2023-05-11