NDVI calculation from Landsat8 images with Python 3 and Rasterio - Tutorial

Satellite images are georasters, these images are a regular array of columns and rows (a matrix per band) with a georeferenciation. Python is a programming and data analysis language very versatile for the matrix algebra with the Numpy library, however there was no efective and simple way to process a georaster until the development of the Rasterio package.
Rasterio is a library to open, write, explore and analyze georasters in Python. The library uses GeoTIFF images along with other formats and is capable to work with satellite images, digital elevation models, and drone generated imagery.
This tutorial show the complete procedure to analyse the NDVI from a Landsat 8 image with Python 3 and Rasterio. The scripting and representation was performed on a interactive enviroment called Jupyter Notebook, finally the result georaster was opened in QGIS and compared with some background images.
Input data
You can download the input data and scripts used for the tutorial on this link:
www.hatarilabs.com/ih-en/ndvi...

[@FRANKFIFORM]

This is exactly what I needed, clear and simple! Thank you so much!

[@raulborda5056]

Excelente video y muchas gracias.

[@alosplanetas]

Super helpful video, thank you!

[@bi4lwt5]

Love your video!

[@olumideakinlabi1519]

This is apparently helpful. Thank you so much

[@smurfk7678]

It's a nice video

[@islandmonkey87]

Are the bands radiometrically corrected prior to ndvi calculation?

[@raulborda5056]

Hola Saul. Es posible con rasterio obtener el corte de un poligono que muestre el NDVI?

[@smurfk7678]

Can you please make videos for landuse landcover classification from multiple Landsat 8 images for a specific study area including agriculture, built-up, water bodies, forest , vegetation cover, snow and glacier etc?

[@hatarilabs]

We will think about ti!

[@c.ravindra3626]

hello sir, can we stack panchromatic band(b8.tif) with (b2,b3,b4). if you have a code or any information related to these please reply.

[@ayebmalek2954]

Hello, can someone helps me to interprate the last plot (yellow and green) please , how can we read the graph ?

[@smurfk7678]

How can we calculate ndvi from multiple images? Can you provide a script?

[@shreyashchoubey4530]

can you save the ndviImage tiff file with cmap ? if yes how?

[@yashpurbey567]

hello sir, everytime i try to open such bands using rasterio and print out my arrays using open and read they all come out as zeros. could you help me out with why thats happening?

[@TheYashO]

Actually I want to ask can we apply ndvi to single plant and get comments if it's healthy or stressed ??

[@luigic6544]

Can some here help me to download rasterio lib?

[@jacobokomo1880]

It is pretty easy to download and install depending on your OS and programming language. in windows, i would recommend you run these commands: pip install geopandas this will install fiona, rTree, rasterio, e.t.c in Linux open the terminal and run these commands : sudo apt-get update sudo apt-get install geoapandas. this will install fiona, rTree, rasterio, e.t.c

[@bhupendrakumar2277]

Hatari, very good job. May God/Allah Bless you. Can you help me solving a problem: I have images with huge volume with total size being around 25 to 190 GB i have to merge them together whiling taking care of no data value. Secondly, from same images i want to create TMS(tiled Map Service ) a lighter and faster version of raster for uploading on any platform- that u might be knowing much better than me. So, could you help in solving these issue.

[@arnabpaul4073]

what kind of merging you want? Kindly clear the problem

[@islandmonkey87]

Are the bands radiometrically corrected prior to ndvi calculation?

[@hatarilabs]

No