def NDVI():
import time
import matplotlib.pyplot as plt
import matplotlib.image as mpimg
from astropy.nddata import NDDataArray
print("This programm will lead you through creating a NDVI map from images")
print("It works under assumption that you input two images into it: A Red band image and Near infrared image")
print("The images must be monochrome, cloudless and it would be best if they had a ocean mask applied.")
print("If they contain clouds/waterbodies, those will also be falsely ranged on scale.")
print("NIR is Band 4, R is Band 3 or broadband VIS")
# pathNIR=input("Please input the path to near infrared image: ")
# pathR=input("Please input the path to red image: ")
print("Please select a color pallete you want to use. Available palletes:")
print("Greens Greys hot seismic")
cmapN = input()
r = mpimg.imread('C:/Users/Karol/Desktop/NIR.png')
nir = mpimg.imread('C:/Users/Karol/Desktop/R.png')
numerator = NDDataArray.subtract(nir, r)
denominator = NDDataArray.add(nir, r)
ndvi = NDDataArray.divide(numerator, denominator)
plt.imshow(ndvi, cmap=cmapN)
plt.colorbar()
plt.show()
def NOAADUST(ch4,ch5):
numerator = NDDataArray.subtract(ch5, ch4)
denominator = NDDataArray.add(ch4, ch5)
dust = NDDataArray.divide(numerator, denominator)
plt.imshow(dust,cmap="hot")
plt.show()
return dust
