ZBLOG

from osgeo import gdal gdal.UseExceptions() driver = gdal.GetDriverByName('GTiff') tifpath = "E:/空间数据统计分析/shuju" outpath = tifpath+"/"+'outtemp' band4 = outpath+'ASTGTMV003_N32E107_dem.tif' band3 = outpath+'ASTGTMV003_N32E106_dem.tif' band2 = ou...

根据你提供的代码,我看到了几个问题需要修正:

  1. 你指定的路径tifpath只是文件夹的路径,而不是包含具体文件名的完整路径。请确保在使用这些路径时加上文件名和扩展名。
band4 = outpath+'/ASTGTMV003_N32E107_dem.tif'
band3 = outpath+'/ASTGTMV003_N32E106_dem.tif'
band2 = outpath+'/ASTGTMV003_N32E105_dem.tif'
  1. 在创建输出数据集时,指定的宽度和高度应该来自于输入数据集(in_ds)的波段大小,而不是从某个波段(如in_band)获取。因此,可以修改为:
out_ds = gtiffDriver.Create(outpath+'/conband001.tif', in_ds.RasterXSize, in_ds.RasterYSize, 3, in_band.DataType)
  1. 将三个波段写入输出数据集时,请确保以正确的顺序进行。按照你当前代码中的顺序,修改如下:
in_ds = gdal.Open(band4) # band4
outband = out_ds.GetRasterBand(3)
outband.WriteArray(in_ds.ReadAsArray())

in_ds = gdal.Open(band3) # band3
outband = out_ds.GetRasterBand(2)
outband.WriteArray(in_ds.ReadAsArray())

in_ds = gdal.Open(band2) # band2
outband = out_ds.GetRasterBand(1)
outband.WriteArray(in_ds.ReadAsArray())

请注意,这只是针对你提供的代码修正了其中的一些问题。如果还存在其他问题,请提供完整的错误信息和相关数据,以便我们更好地帮助你进行调试。


内容由零声教学AI助手提供,问题来源于学员提问

本站部分文章来源于网络,版权归原作者所有,如有侵权请联系站长删除。
转载请注明出处:https://golang.0voice.com/?id=14989

分享:
扫描分享到社交APP
上一篇
下一篇
发表列表
游客 游客
此处应有掌声~
评论列表

还没有评论,快来说点什么吧~

联系我们

在线咨询: 点击这里给我发消息

微信号:3007537140

上班时间: 10:30-22:30

关注我们
x

注册

已经有帐号?