Xác định giá trị tối thiểu và tối đa trong bộ dữ liệu raster ASCII bằng Python?

Tôi có một bộ dữ liệu raster ở định dạng ASCII. Sử dụng Python, tôi cần xác định minvà maxcác giá trị bên trong tập dữ liệu. Tôi đã được thông báo rằng thông tin tiêu đề là khóa, chứa những thứ như số hàng / cột, kích thước ô, v.v.

Bạn có thể đơn giản bỏ qua thông tin tiêu đề và đọc toàn bộ dữ liệu để xác định minvà maxgiá trị không?

Đây là những gì tôi đang cố gắng làm. Tôi bỏ qua vài dòng đầu tiên chứa thông tin tiêu đề và cố gắng xác định các giá trị từ đó trở đi. Sau đây là những gì tôi có, nhưng cần một số hướng dẫn vì tôi chưa quen với 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

Bất kỳ đề xuất?

Bạn có thể sử dụng numpy. Xem ví dụ dưới đây. Một mảng mặt nạ numpy có thể được tạo ra chiếm không có giá trị dữ liệu. Xem chủ đề trợ giúp numpy cho mafromtxt và 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]]

>>>

từ đó đơn giản chỉ là vấn đề xác định số liệu thống kê bạn muốn

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

Bạn muốn thống kê dữ liệu raster.
Xem những gì bạn đang làm trong gui đầu tiên (cho bài tập về nhà.)

Sau đó, bạn có thể sử dụng một cửa sổ python hoặc một kịch bản .

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

Nếu bạn không muốn sử dụng numpy (và bạn thực sự nên, nó hoàn hảo cho loại điều này), thì bạn sẽ cần phải:

• khởi tạo maximumbiến của bạn thành một số âm rất lớn và minimumbiến của bạn thành một số dương rất lớn
• tách từng dòng để có được một danh sách các chuỗi và sử dụng sự hiểu biết danh sách để chuyển đổi nó thành một danh sách các float
• cuối cùng sử dụng một cái gì đó như maximum = max(maximum, max(myfloatlist))và tương đương với giá trị tối thiểu.

Tôi chỉ làm điều này vào ngày khác. Tôi đã sử dụng arcpy.RasterToNumPyArray, chuyển đổi mảng numpy thành một danh sách, sau đó lặp qua danh sách của mình thông qua việc hiểu danh sách để tìm các giá trị tối thiểu và tối đa.

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