Déterminer les valeurs min et max dans un jeu de données raster ASCII à l'aide de Python?

J'ai un jeu de données raster au format ASCII. En utilisant Python, j'ai besoin de déterminer les valeurs minet maxà l'intérieur de l'ensemble de données. On m'a dit que les informations d'en-tête sont essentielles, qui contiennent des éléments comme le nombre de lignes / colonnes, la taille des cellules, etc.

Vous ne pouvez pas simplement ignorer les informations d'en-tête et lire l'ensemble de données complet pour déterminer les valeurs minet max?

Voilà ce que j'essaie de faire. J'ignore les deux premières lignes qui contiennent les informations d'en-tête et j'essaie de déterminer les valeurs à partir de là. Ce qui suit est en quelque sorte ce que j'ai, mais j'ai besoin de quelques conseils car je suis nouveau sur Python.

raster_file = open('data.asc', 'r') # Open the file
data = raster_file.readlines()[4:] # Read the lines in the file, and skip the first six lines

for lines in data:
    print max(data) # Find the max value in data
    print min(data) # Find hte min value in data

Aucune suggestion?

Vous pouvez utiliser numpy. Voir l'exemple ci-dessous. Un tableau masqué numpy peut être généré en tenant compte des valeurs sans données. Voir la rubrique d'aide numpy pour mafromtxt et genfromtxt

Below is a small ascii file with a nodata value of -999

ncols          3
nrows          3
xllcorner      0
yllcorner      0
cellsize       1
NODATA_value   -999
0 1 2
-999 4 5 
6 7 8

>>> import numpy as np
>>> ascii_file = "c:/temp/Ascii_3x3_1nodata.asc"
>>> an_array = np.mafromtxt(ascii_file, 'float', '#', None, 6, None, '-999')

>>> print an_array

[[0.0 1.0 2.0]
 [-- 4.0 5.0]  
 [6.0 7.0 8.0]]

>>>

à partir de là, il s'agit simplement de déterminer les statistiques que vous souhaitez

>>> print an_array.min()
0.0
>>> print an_array.max()
8.0
>>> print an_array.mean()
4.125
>>> 

Vous voulez des statistiques de données raster.
Voyez d'abord ce que vous faites dans l'interface graphique (pour les devoirs.)

Ensuite, vous pouvez utiliser une fenêtre python ou un script .

import arcpy
arcpy.CalculateStatistics_management("c:/data/image.tif", "4", "6", "0;255;21")

import sys

class Ascii_file(object):
    def __init__(self,file):
        self.raster_file = open(file, 'r') # Open the file
        self.max=sys.float_info.min
        self.min=sys.float_info.max
    def __minmax(self,value):
        if value>self.max:self.max=value
        if value<self.min:self.min=value
    def getMinMax(self):
        data = self.raster_file.readlines()
        data_values=data[6:]
        nodata=float(data[5].split()[1])
        for line in data_values:
            values=line.split(" ")
            for value in values:
                value=float(value)
                if value==nodata:continue
                else: self.__minmax(value)
        return self.min, self.max

if __name__=="__main__":
    myfile = Ascii_file('data.asc')
    print myfile.getMinMax()

Si vous ne voulez pas utiliser numpy (et vous devriez vraiment le faire, c'est parfait pour ce genre de chose), alors vous devrez:

• initialisez votre maximumvariable à un très grand nombre négatif et votre minimumvariable à un très grand nombre positif
• divisez chaque ligne pour obtenir une liste de chaînes et utilisez la compréhension de liste pour la convertir en une liste de flottants
• enfin utiliser quelque chose comme maximum = max(maximum, max(myfloatlist))et un équivalent pour la valeur minimale.

Je viens de le faire l'autre jour. J'ai utilisé arcpy.RasterToNumPyArray, converti le tableau numpy en liste, puis parcouru ma liste via une compréhension de liste pour trouver les valeurs min et max.

import arcpy
import numpy
myArray = arcpy.RasterToNumPyArray(r"D:\NED_93512417\NED_93512417_3DEM_RPRJ.TIF")
p = myArray.tolist()

max_elev = max([item for sublist in p for item in sublist])
min_elev = min([item for sublist in p for item in sublist])