如何使用GDAL将操作的栅格值写入ASCII网格？

Question

我试图用GDAL操作网格（ASCII网格）中的栅格值。 但在继续执行此操作之前，我无法将新值写入文件。 调用slopeband.WriteArray时，我收到这些错误消息。

错误6：此数据集不支持slope.asc，band 1：WriteBlock（）。

错误1：slope.asc，band 1：写入脏块时发生错误

我很抱歉，如果这是非常基本的，但我仍然是python和GDAL的新手。 我在Mac OS X 10.6.8和Python 2.7上使用GDAL 1.9.0。 谢谢！

import numpy
import gdal
import gdalconst

dgm = gdal.Open("DGM_10_MR.asc", gdalconst.GA_ReadOnly)
driver = dgm.GetDriver()
geotransform = dgm.GetGeoTransform()
band = dgm.GetRasterBand(1)
data = band.ReadAsArray()

cols = dgm.RasterXSize
rows = dgm.RasterYSize
slope = driver.CreateCopy("slope.asc", dgm)
slope = None
dgm = None
slope = gdal.Open("slope.asc", gdalconst.GA_Update)
slope.SetGeoTransform(geotransform)
slopeband = slope.GetRasterBand(1)
s = slopeband.ReadAsArray()

for y in range(rows):
    for x in range(cols):
        s[y, x] = 0.0

slopeband.WriteArray(s)
slopeband.FlushCache()
del s

dgm = None
slope = None
print "done"

Answer 1

不幸的是，GDAL无法在所有文件类型中读取和写入相同的度数。 Arc ASCII网格恰好是GDAL无法写入的文件类型之一。 正如您的错误消息所示： WriteBlock() not supported for this dataset. ，因此您无法写入Arc ASCII网格。

作为替代方案，您可以将现有的ASCII数据集转换为不同的文件类型，这是GDAL更完全支持的文件类型，例如GeoTiff。 要转换文件类型，您可以使用gdal_translate命令行程序，如下所示：

gdal_translate -of GTiff DGM_10_R.asc DGM_10_R.tif

我能够在我的计算机上重现您的错误，只需更改文件类型即可修复错误。