for i in range(0, rows, yBlockSize):
if i + yBlockSize < rows:
numRows = yBlockSize
else:
numRows = rowsnumRows = rows –– ii
for j in range(0, cols, xBlockSize):
if j + xBlockSize < cols:
numCols = xBlockSize
else:
numCols = cols – j
data = band.ReadAsArray(j, i, numCols, numRows)
# do calculations here to create outData array
outBand.WriteArray(outData, j, i)
band对象可以设定NoData值
outBand.SetNoDataValue(-99)
还可以读取NoData值
ND = outBand.GetNoDataValue()
计算band的统计量
首先用FlushCache()把缓存数据写入磁盘
之后用GetStatistics(<approx_ok>, <force>)计算统计量。如果approx_ok=1那么计算是基于pyramid的，如果force=0那么当整幅图都要被重读一遍的时候就不计算统计量了。
outBand.FlushCache()
outBand.GetStatistics(0, 1)
设定新图的地理参考点
如果新图与另一张图的地理参考信息完全一致，那就很简单了
geoTransform = inDataset.GetGeoTransform()
outDataset.SetGeoTransform(geoTransform )
proj = inDataset.GetProjection()
outDataset.SetProjection(proj)
建立pyramids
设定Imagine风格的pyramids
gdal.SetConfigOption('HFA_USE_RRD', 'YES')
强制建立pyramids
outDataset.BuildOverviews(overviewlist=[2,4, 8,16,32,64,128])
图像的拼接
1. 对每张图：读取行数和列数，原点(minX,maxY)，像素长，像素宽，并计算坐标范围
maxX1 = minX1 + (cols1 * pixelWidth)
minY1 = maxY1 + (rows1 * pixelHeight)
2. 计算 输出图像的坐标范围：
minX = min(minX1, minX2, …) maxX = max(maxX1, maxX2, …)
minY = min(minY1, minY2, …) maxY = max(maxY1, maxY2, …)
3. 计算 输出图像的行数和列数：
cols = int((maxX – minX) / pixelWidth)
rows = int((maxY – minY) / abs(pixelHeight)
4. 建立并初始化 输出图像
5. 对每张待拼接的图：计算offset值
xOffset1 = int((minX1 - minX) / pixelWidth)
yOffset1 = int((maxY1 - maxY) / pixelHeight)
读入数据并按照上面计算的offset写入
6. 对 输出图像：计算统计量，设定geotransform ：[minX, pixelWidth, 0, maxY, 0, pixelHeight]，设定projection，建立pyramids
以上就是python gdal教程之：地图代数与栅格数据的写入的内容，更多相关内容请关注PHP中文网（HdhCmsTestgxlcms测试数据）！
查看更多关于pythongdal教程之：地图代数与栅格数据的写入的详细内容...