Rutina de Python para obtener los valores de Percentil...

Rutina de Python para obtener los valores de Percentil...

de Luis Alberto Calderón Ruiz -
Número de respuestas: 0

# Importar biblioteca de Python.

import numpy as np


# Valor de Percentil buscado.

vValorPercentil = 5.0


# Raster del analisis en la Tabla de Contenidos

vNombreRaster = "rasSolucion"


# Que valor le vamos a dar al NoData ?

vValorNoData = -1000


# Convierto mi raster en una matriz de dobles.

vArray = arcpy.RasterToNumPyArray(vNombreRaster, nodata_to_value = vValorNoData)


# Filtro la matriz: excluyo los valores que son iguales a NoData.

vArray = vArray [vArray != vValorNoData]


# Ordeno la Matriz.

vArray = np.sort(vArray)


# Calculo el percentil.

vPercentil = vArray[int(np.size(vArray) * vValorPercentil / 100.0)]


# Presento el Resultado.

print('Valor de Percentil Buscado = ', vPercentil)