3

このコードを適応させて、LANDSAT衛星画像をGrass GISにインポートするPythonコードを作成しようとしています:http://grass.osgeo.org/wiki/LANDSAT

LANDSATタイルはフォルダーとしてダウンロードされ、各フォルダーには7つのtiff画像(バンド1〜7)が含まれています。したがって、いくつかのサブディレクトリ(LANDSATタイルごとに1つ)を含むディレクトリがあります。

現在の私のコードは次のとおりです。

#!/usr/bin/python

import os
import sys
import glob
import grass.script as grass

def import_tifs(dirpath):

    for dirpath, dirname, filenames in os.walk(dirpath):
        for dirname in dirpath:

        dirname = os.path.join(dirpath,dirname)

            for file in os.listdir(dirname):
                if os.path.splitext(file)[-1] != '.TIF':
                    continue
                ffile = os.path.join(dirname, file)
                name = os.path.splitext(file)[0].split(dirname)[-1]

                grass.message('Importing %s -> %s@%s...' % (file, name, dirpath))

                grass.run_command('r.in.gdal',
                                  flags = 'o',
                                  input = ffile,
                                  output = name,
                                  quiet = True,
                                  overwrite = True)

def main():                                 
    if len(sys.argv) == 1:
        for directory in filter(os.path.isdir, os.listdir(os.getcwd())):
            import_tifs(directory)
    else:
        import_tifs(sys.argv[1])

if __name__ == "__main__":
    main()

次のエラーが発生します。

Traceback (most recent call last):
  File "C:/Users/Simon/Documents/import_landsat2.py", line
40, in <module>
    main()
  File "C:/Users/Simon/Documents/import_landsat2.py", line
37, in main
    import_tifs(sys.argv[1])
  File "C:/Users/Simon/Documents/import_landsat2.py", line
17, in import_tifs
    for file in os.listdir(dirname):
WindowsError: [Error 3] The system cannot find the path
specified: 'dirpath\\C/*.*'

誰かが何が起こっているのか、それを修正するために私が何をする必要があるのか​​を説明したり、代替案を提案したりできますか?ありがとう。

4

2 に答える 2

2

私はあなたの主な問題は、がリスト(文字列ではない)dirnameos.walk()返すことであると信じています。そのため、後続の文字列(つまりdirname = os.path.join(dirpath,dirname))は少し奇形です。考えられる代替案の1つを次に示します。これをテストするために、ディレクトリへのフルパスをとして使用しましたsys.argv[1]が、ケースに合わせてより動的にすることができます。fileまた、 Pythonキーワードであるため、などの変数名は使用しないでください。私はあなたの関数をテストすることができませんでしgrass.*たが、うまくいけば、これはあなたが必要な方法を微調整できるように十分に明確な例になるでしょう。os.walk()多くの標準的な解析をネイティブに処理するため、ディレクトリ操作機能の一部を削除できます。

def import_tifs(dirpath):
  for dirpath, dirname, filenames in os.walk(dirpath):
    # Iterate through the files in the current dir returned by walk()
    for tif_file in filenames:
      # If the suffix is '.TIF', process
      if tif_file.upper().endswith('.tif'):
        # This will contain the full path to your file
        full_path = os.path.join(dirpath, tif_file)

        # tif_file will already contain the name, so you can call from here
        grass.message('Importing %s -> %s@%s...' % (full_path, tif_file, dirpath))

        grass.run_command('r.in.gdal',
                          flags = 'o',
                          input = full_path,
                          output = tif_file,
                          quiet = True,
                          overwrite = True)
于 2012-10-25T17:56:42.277 に答える
2

コードを書き直して、すべてのdirツリーを一覧表示し、ファイル拡張子(この場合は「.tif」)を見つけました。

#!/usr/bin/python

import os
import sys


def import_tifs(dirpath):

    for dirpath, dirname, filenames in os.walk(dirpath):
        for filename in filenames:
            name, extension = os.path.splitext(filename)
            if extension.lower() == ".tif":
                filepath = os.path.join(dirpath, filename)

                print(filepath, name, dirpath)


def main():
    if len(sys.argv) == 1:
            import_tifs(os.getcwd())
    else:
        import_tifs(sys.argv[1])


if __name__ == "__main__":
    main()

それがあなたが探しているものであるかどうかを確認してください...

于 2012-10-25T18:04:45.033 に答える