Estoy desarrollando un software para calcular NDVI. Que es NDVI? es un calculo matemático con el cual se puede determinar el nivel de estrés de la vegetación NDVI = NIR-RED/NIR+RED Este resultado varia de -1 a 1 Estoy usando OpenCV como librería para visión artificial y Gdal para ayudar a los cálculos.
Mi idea es poder "colorear" los rangos que entregue NDVI, por ejemplo para valores < 0 colorear esos pixeles con un color rojo.
Mi problema es que no puedo recorrer la matriz de datos.
Les adjunto el codigo:
# Importamos Librerias a utilizar
from osgeo import gdal
import cv2
import numpy as np
# Abrimos la imagen a utilizar
img = cv2.imread('image01.JPG')
# Pasamos de RGB a HSV
hsv = cv2.cvtColor(img, cv2.COLOR_BGR2HSV)
# Claros
lw_range = np.array([4, 150, 0])
# Oscuros
up_range = np.array([34, 255, 255])
#Crear una mascara con solo los pixeles dentro del rango de verdes
mask = cv2.inRange(hsv, lw_range, up_range)
# resta para mostrar imagen original
res = cv2.bitwise_and(img,img, mask= mask)
cv2.imwrite('imglimpia.tif',res)
# Procesamiento con Gdal
imgGedal = gdal.Open('imglimpia.tif', gdal.GA_ReadOnly)
num_bands = imgGedal.RasterCount
rows = imgGedal.RasterYSize
cols = imgGedal.RasterXSize
#print('Numero de Bandas {n}\n'.format(n=num_bands))
# Separamos por banda
blue = imgGedal.GetRasterBand(3)
green = imgGedal.GetRasterBand(2)
red = imgGedal.GetRasterBand(1)
blue_data = blue.ReadAsArray()
green_data = green.ReadAsArray()
red_data = red.ReadAsArray()
cv2.imwrite('azul.tif',blue_data)
cv2.imwrite('gren.tif',green_data)
cv2.imwrite('red.tif',red_data)
# Obtener Valores Minimo y Maximo
band_min, band_max, band_mean, band_stddev = blue.GetStatistics(0, 1)
print('Rango Azul: Minimo: {minimo} - Medio: {medio} - Maximo: {maximo} - Tamano: {sz}- STDDEV: {stddev}'
.format(maximo=band_max,medio=blue_data.mean(),minimo=band_min,sz=blue_data.shape,stddev=band_stddev))
band_min, band_max, band_mean, band_stddev = green.GetStatistics(0, 1)
print('Rango Verde: Minimo: {minimo} - Medio: {medio} - Maximo: {maximo} - Tamano: {sz}- STDDEV: {stddev}'
.format(maximo=band_max,medio=blue_data.mean(),minimo=band_min,sz=blue_data.shape,stddev=band_stddev))
band_min, band_max, band_mean, band_stddev = red.GetStatistics(0, 1)
print('Rango Rojo: Minimo: {minimo} - Medio: {medio} - Maximo: {maximo} - Tamano: {sz}- STDDEV: {stddev}'
.format(maximo=band_max,medio=blue_data.mean(),minimo=band_min,sz=blue_data.shape,stddev=band_stddev))
# Calcular NDVI = ((Nir-Blue)/(Nir+Blue))
resta = cv2.subtract(red_data,blue_data)
suma = cv2.add(red_data,blue_data)
ndvi = cv2.divide(resta,suma)
print 'NDVI MIN : ',ndvi.min()
print 'NDVI MEDIO : ',ndvi.mean()
print 'NDVI MAX : ',ndvi.max()
Luego de esto, quiero recorrer el array de datos para colorear pixeles
# Creamos los colores
greenyellow = [4,255,173]
yellow = [0,255,255]
yellowred = [4,111,255]
green = [0,255,0]
greenyellow = [0,255,146]
blue = [255,0,0]
red = [0,0,255]
white = [255,255,255]
black = [0,0,0]
y = np.zeros((4000,3000,3),np.uint8) # img coloreada
cols,rows,band = y.shape
for i in range(rows):
for j in range(cols):
for k in ndvi[i,j]:
print (str([i,j,k]) + '\n')
if k < 0.0 :
y[i,j] = red
break
if 0.0 < k < 0.2 :
y[i,j] = yellowred
break
if 0.2 < k < 0.4 :
y[i,j] = yellow
break
if 0.4 < k < 0.6 :
y[i,j] = greenyellow
break
if k > 0.6 :
y[i,j] = green
break
else:
y[i,j] = black
break
cv2.imwrite('NDVI.jpg',y)
este es el error que me genera:
for k in ndvi[j,i]:
TypeError: 'numpy.uint8' object is not iterable
numpy
. Nigdal
niopenCV
te dan error.