2

我使用 Sentinel2 图像,我正在尝试重新采样它们。

我尝试了以下代码:

import os, fnmatch

INPUT_FOLDER = "/d/afavro/Bureau/test_resampling/original"
OUTPUT_FOLDER = "/d/afavro/Bureau/test_resampling/resampling_10m"

    def findRasters (path, filter):
        for root, dirs, files in os.walk(path):
            for file in fnmatch.filter(files, filter):
                yield file

    for raster in findRasters(INPUT_FOLDER,'*.tif'):
        print(raster)
        inRaster = INPUT_FOLDER + '/' + raster
        print(inRaster)
        outRaster = OUTPUT_FOLDER + '/resample' + raster
        print (outRaster)
        cmd = "gdalwarp -tr 10 10 -r cubic " % (inRaster,outRaster)
        os.system(cmd)

但我仍然收到相同的错误消息:

def findRasters (path, filter): ^
IndentationError: unexpected indent

我已经尝试过使用相同类型的代码来制作一个子集并且它有效。我不明白我的错误来自哪里。

4

1 回答 1

1

错误类型IndentationError应该从字面上理解:您的缩进似乎是错误的。你的线

def findRasters (path, filter):

缩进太远,但需要与上一行的缩进级别相同

OUTPUT_FOLDER = "/d/afavro/Bureau/test_resampling/resampling_10m"

您提供的完整代码示例应如下所示:

import os, fnmatch

INPUT_FOLDER = "/d/afavro/Bureau/test_resampling/original"
OUTPUT_FOLDER = "/d/afavro/Bureau/test_resampling/resampling_10m"

def findRasters (path, filter):
    for root, dirs, files in os.walk(path):
        for file in fnmatch.filter(files, filter):
            yield file

for raster in findRasters(INPUT_FOLDER,'*.tif'):
    print(raster)
    inRaster = INPUT_FOLDER + '/' + raster
    print(inRaster)
    outRaster = OUTPUT_FOLDER + '/resample' + raster
    print (outRaster)
    cmd = "gdalwarp -tr 10 10 -r cubic " % (inRaster,outRaster)
    os.system(cmd)

此外,正如您在附加评论中所写,您的行

cmd = "gdalwarp -tr 10 10 -r cubic " % (inRaster,outRaster)

似乎是错误的,inRaster并且outRaster不会在字符串中使用。改用字符串格式

cmd = 'gdalwarp -tr 10 10 -r cubic "{}" "{}"'.format(inRaster, outRaster)
于 2019-04-24T11:40:55.910 回答