Python 如何在循环中链接堆栈
我正在尝试创建一个代码,允许我为不同文件夹中包含的不同图像堆叠同一图像的10个波段。 我正在尝试创建一个代码,允许我为不同文件夹中包含的不同图像堆叠同一图像的10个波段。 我在论坛上找到了一段代码,我试图从中获得灵感,但我无法将其应用到我的案例中。。。 有关守则如下:Python 如何在循环中链接堆栈,python,stack,sentinel2,Python,Stack,Sentinel2,我正在尝试创建一个代码,允许我为不同文件夹中包含的不同图像堆叠同一图像的10个波段。 我正在尝试创建一个代码,允许我为不同文件夹中包含的不同图像堆叠同一图像的10个波段。 我在论坛上找到了一段代码,我试图从中获得灵感,但我无法将其应用到我的案例中。。。 有关守则如下: from osgeo import gdal, gdal_array import numpy as np b1 = gdal.Open("LT50250232011160PAC01_sr_band1.tif") b2 = gda
from osgeo import gdal, gdal_array
import numpy as np
b1 = gdal.Open("LT50250232011160PAC01_sr_band1.tif")
b2 = gdal.Open("LT50250232011160PAC01_sr_band2.tif")
array1 = b1.ReadAsArray()
array2 = b2.ReadAsArray()
stacked = np.array([array1,array2])
gdal_array.SaveArray(stacked.astype("int"), "b12_stacked.tif", "GTiff", gdal.Open("LT50250232011160PAC01_sr_band1.tif"))
我的代码是:
import os
from osgeo import gdal, gdal_array
import numpy as np
src_directory = "/d/afavro/Bureau/test_subset/stack_resample_subset/"
dossier = os.listdir (src_directory)
print(dossier)
for fichier in dossier:
print (fichier)
ssrc_directory = "/d/afavro/Bureau/test_subset/stack_resample_subset/" + fichier + "/"
rasters = os.listdir (ssrc_directory)
print(rasters)
for img in rasters:
print (img)
INPUT_FOLDER = ssrc_directory
print(INPUT_FOLDER)
OUTPUT_FOLDER = "/d/afavro/Bureau/test_stack/" + fichier
print(OUTPUT_FOLDER)
if not os.path.exists(OUTPUT_FOLDER):
os.makedirs(OUTPUT_FOLDER)
# Define output file name
if img[-4:-3] == '.' :
file_img = OUTPUT_FOLDER+'/'+os.path.basename(img[:-4])
else:
file_img = OUTPUT_FOLDER+'/'+os.path.basename(img)
# Stack
band = gdal.Open(img)
array = band.ReadAsArray()
stacked = np.array(array)
gdal_array.SaveArray(stacked.astype("int"),"stack_" + fichier, "GTiff", prototype = None )
问题就在这里,我收到了错误消息:
band = gdal.Open(img)
^
IndentationError: unindent does not match any outer indentation level
我认为我们需要创建一个空矩阵或类似的东西,但我真的不知道如何创建,你有什么想法吗?错误通知很清楚,是关于最后4行代码的缩进。如果它们位于光栅中img的
循环中,请为它们再添加一个空间:
循环。哦,是的,对不起!我看不见,但有个错误。。。现在我有以下消息:AttributeError:'NoneType'对象没有属性'ReadAsArray'。错误通知很清楚,这是关于最后4行代码的缩进。如果它们位于光栅中img的循环中,请为它们再添加一个空间:
循环。哦,是的,对不起!我看不见,但有个错误。。。现在我有以下消息:AttributeError:'NoneType'对象没有属性'ReadAsArray'