このコードを適応させて、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/*.*'
誰かが何が起こっているのか、それを修正するために私が何をする必要があるのかを説明したり、代替案を提案したりできますか?ありがとう。