Determinare i valori minimo e massimo in un set di dati raster ASCII usando Python?

Ho un set di dati raster in formato ASCII. Usando Python, devo determinare i valori mine maxall'interno del set di dati. Mi è stato detto che le informazioni di intestazione sono fondamentali, che contengono elementi come il numero di righe / colonne, dimensioni delle celle e così via.

Non puoi semplicemente saltare le informazioni dell'intestazione e leggere l'intero set di dati per determinare i valori mine max?

Questo è quello che sto cercando di fare. Sto saltando la prima coppia di righe che contengono le informazioni di intestazione e provo a determinare i valori da lì in poi. Quello che segue è una specie di quello che ho, ma ho bisogno di una guida dato che sono nuovo di 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

Eventuali suggerimenti?

Puoi usare numpy. Vedi l'esempio sotto. Un array mascherato intorpidito può essere generato tenendo conto dei valori di nessun dato. Vedere l'argomento della guida numpy per mafromtxt e 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]]

>>>

da lì è semplicemente una questione di determinare le statistiche desiderate

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

Vuoi statistiche sui dati raster.
Guarda prima cosa stai facendo nella gui (per i compiti.)

Quindi puoi usare una finestra di Python o uno 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()

Se non vuoi usare numpy (e dovresti davvero, è perfetto per questo genere di cose), allora dovrai:

• inizializza la tua maximumvariabile su un numero negativo molto grande e la tua minimumvariabile su un numero positivo molto grande
• dividere ogni riga per ottenere un elenco di stringhe e utilizzare la comprensione dell'elenco per convertirlo in un elenco di float
• infine usa qualcosa di simile maximum = max(maximum, max(myfloatlist))e un equivalente per il valore minimo.

L'ho fatto l'altro giorno. Ho usato arcpy.RasterToNumPyArray, convertito l'array numpy in un elenco, quindi ho ripetuto il mio elenco tramite una comprensione dell'elenco per trovare i valori minimo e massimo.

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])