4

我正在尝试将 TIFF 的波段读取为数组。问题是 GDAL 不会忽略 NoData 值。有没有办法告诉GDAL忽略它？

当我计算统计数据时，GDAL 会忽略 NoData 值。

import os, sys, gdal, ogr, numpy
from gdalconst import *

# register all of the drivers
gdal.AllRegister()

# open the image
ds = gdal.Open('test_slope.tif', GA_ReadOnly)

# get image size
rows = ds.RasterYSize
cols = ds.RasterXSize
bands = ds.RasterCount

# Set NoData Value
band = ds.GetRasterBand(1)
ndv = -3.4028230607371e+38
band.SetNoDataValue(ndv)

# Get Statistics
stats = band.ComputeStatistics(0)
print stats

# read in band as array
bandList = []
band.GetNoDataValue()
data = band.ReadAsArray(0, 0, cols, rows)
print data


>>> 
[0.0, 126.59918975830078, 25.757117870945123, 15.399812314100501]
[[ -3.40282306e+38  -3.40282306e+38  -3.40282306e+38 ...,  -3.40282306e+38
   -3.40282306e+38  -3.40282306e+38]
 [ -3.40282306e+38  -3.40282306e+38  -3.40282306e+38 ...,  -3.40282306e+38
   -3.40282306e+38  -3.40282306e+38]
 [ -3.40282306e+38  -3.40282306e+38  -3.40282306e+38 ...,  -3.40282306e+38
   -3.40282306e+38  -3.40282306e+38]
 ..., 
 [ -3.40282306e+38  -3.40282306e+38  -3.40282306e+38 ...,  -3.40282306e+38
   -3.40282306e+38  -3.40282306e+38]
 [ -3.40282306e+38  -3.40282306e+38  -3.40282306e+38 ...,  -3.40282306e+38
   -3.40282306e+38  -3.40282306e+38]
 [ -3.40282306e+38  -3.40282306e+38  -3.40282306e+38 ...,  -3.40282306e+38
   -3.40282306e+38  -3.40282306e+38]]
>>>
4

1 回答 1

4

我相信你可以numpy.MaskedArray从你的numpy.ndarray：

import numpy as np
import numpy.ma as ma

ndv = -3.40282306e+38

data = np.array([[0.0, 126.59918975830078, 25.757117870945123, 15.399812314100501],
                [-3.40282306e+38, -3.40282306e+38, -3.40282306e+38, -3.40282306e+38]])

masked_data = ma.masked_where(data == ndv, data)
print masked_data

结果：

[[0.0 126.599189758 25.7571178709 15.3998123141]
[-- -- -- --]]

在numpy中，anumpy.ndarray和掩码的组合用于处理丢失的数据。

于 2013-03-24T20:57:51.937 回答