3

您好，我是 Gdal 的新手，我正在为我的代码苦苦挣扎。在我的代码中一切似乎都很顺利，但最后的输出波段是空的。当我指定 255 时，无数据值设置为 256，所以我真的不知道出了什么问题。感谢任何帮助将不胜感激！！！

这是我的代码

from osgeo import gdal
from osgeo import gdalconst
from osgeo import osr
from osgeo import ogr
import numpy

#graticule
src_ds = gdal.Open("E:\\NFI_photo_plot\\photoplotdownloadAllCanada\\provincial_merge\\Aggregate\\graticule1.tif")
band = src_ds.GetRasterBand(1)
band.SetNoDataValue(0)
graticule = band.ReadAsArray()
print('graticule done')
band="none"

#Biomass
dataset1 = gdal.Open("E:\\NFI_photo_plot\\photoplotdownloadAllCanada\provincial_merge\\Aggregate\\Biomass_NFI.tif")
band1 = dataset1.GetRasterBand(1)
band1.SetNoDataValue(-1)
Biomass = band1.ReadAsArray()
maskbiomass = numpy.greater(Biomass, -1).astype(int)
print("biomass done")
Biomass="none"
band1="none"
dataset1="none"

#Baseline
dataset2 = gdal.Open("E:\\NFI_photo_plot\\Baseline\\TOTBM_250.tif")
band2 = dataset2.GetRasterBand(1)
band2.SetNoDataValue(0)
baseline = band2.ReadAsArray()
maskbaseline = numpy.greater(baseline, 0).astype(int)
print('baseline done')
baseline="none"
band2="none"
dataset2="none"

#sommation
biosource=(graticule+maskbiomass+maskbaseline)
biosource1=numpy.uint8(biosource)
biosource="none"

#Écriture
dst_file="E:\\NFI_photo_plot\\photoplotdownloadAllCanada\\provincial_merge\\Aggregate\\Biosource.tif"
dst_driver = gdal.GetDriverByName('GTiff')
dst_ds = dst_driver.Create(dst_file, src_ds.RasterXSize,
src_ds.RasterYSize, 1, gdal.GDT_Byte)

#projection
dst_ds.SetProjection( src_ds.GetProjection() )
dst_ds.SetGeoTransform( src_ds.GetGeoTransform() )
outband=dst_ds.GetRasterBand(1)
outband.WriteArray(biosource1,0,0)
outband.SetNoDataValue(255)
biosource="none"
graticule="none"
4

1 回答 1

1

几点建议：

  • 在您拥有 的地方="none"，这些需要= None关闭/清理对象，否则您将对象设置为字符数组：n o n e，这不是您打算做的。
  • 为什么你有band1.SetNoDataValue(-1)，而其他 NoData 值为 0？此数据源是签名还是未签名？如果无符号，则 -1 不存在。
  • 当您在gdal.Open没有访问选项的情况下打开栅格时，它默认为gdal.GA_ReadOnly，这意味着您的后续SetNoDataValue调用什么也不做。如果要修改数据集，则需要gdal.GA_Update将gdal.Open.
  • 创建新栅格的另一种策略是使用driver.CreateCopy; 有关详细信息，请参阅教程。
于 2012-01-29T22:04:04.190 回答