Create a new image $[2-$D NumPy Array$]$ and fill it with data following this NDVI equation:

NDVI $=($NIR $$ RED$)/($NIR $+$ RED$)$

You will need to decide which of the bands to use for NIR and Red.

NOTE: Make sure you do not divide by ZERO, so you will need an exception handler $($logic in your code to skip division by ZERO$).$ It would be great if you could design a separate user defined function that you would call.

Plot this image and use the $$YlGn$$ or $$Greens$$ cmap colorbar$/$scale and add the color legend, labels, etc. to the image.

Scatterplot NDVI values for the objects, set the y$-$axis to NDVI and the x$-$axis to the index or a label of the point.