複数のgeotiffを1つに集約したいができない

Question

現在複数のtiffファイルを1つに集約しようと下記のコードを書き、実行しているのですが、下記のようなエラーが発生し、どのように解決すればよいのかがわからず困っております。

#発生しているerror

0...10...20...30...40ERROR 5: temp\28363_13157_15.tif: GDALDataset::GetRasterBand(4) - Illegal band #

Traceback (most recent call last):
  File "C:\Users\lyon_\AppData\Local\Programs\Python\Python36\Scripts\gdal_merge.py", line 611, in <module>
    sys.exit(main())
  File "C:\Users\lyon_\AppData\Local\Programs\Python\Python36\Scripts\gdal_merge.py", line 596, in main
    fi.copy_into(t_fh, band, band, nodata)
  File "C:\Users\lyon_\AppData\Local\Programs\Python\Python36\Scripts\gdal_merge.py", line 343, in copy_into
    nodata_arg)
  File "C:\Users\lyon_\AppData\Local\Programs\Python\Python36\Scripts\gdal_merge.py", line 116, in raster_copy
    if s_band.GetMaskFlags() != gdal.GMF_ALL_VALID:
AttributeError: 'NoneType' object has no attribute 'GetMaskFlags'

#実行したコード

def merge_tiles(input_pattern, output_path):
    merge_command = ['gdal_merge.py', '-o', output_path]

    for name in glob.glob(input_pattern):
        merge_command.append(name)

    subprocess.call(merge_command , shell=True)


merge_tiles(temp_dir + '/*.tif', output_dir + '/merged.tif')

#考えたこと
調べても具体的にどのように修正すれば良いのかがわからず、具体的に修正したことはないです。
ただ、1点疑問に思ったことがあり、tiff画像が位置的に連続していないことが原因なのではないのかなと思いました。
というのは、現在使用しているtigg画像は、九州のある地域の一部だけです。下記の画像にあるように、九州の赤い部分で
囲まれたtiff画像を現在持っており、これらを1つに集約しようと考えているのですが、そもそもこの考え方がダメなのでしょうか。

<取得し、加工したpng画像>


<エラーが発生しているpng画像>

わからないことばかりで質問の意図がわかりくいかもしれませんが、何卒よろしくお願いたします。

Comments

[yuki23]

エラーメッセージは「28363_13157_15.tifがおかしい」といっていますが、地図でいうとどこの画像でしょうか？

[PINE1103]

お返事ありがとうございます。
エラー発生箇所について画像を添付させていただきましたのでご確認のほどよろしくお願いいたします。

https://maps.gsi.go.jp/#11/33.233519/131.458282/&base=std&ls=std%7Cafm&blend=0&disp=11&vs=c1j0h0k0l0u0t0z0r0s0m0f1

補足ですが、上記のサイトを利用しており、zoomする度合によって切り抜ける画像の大きさが変わってきます。

[yuki23]

画像ファイルがないと問題を再現できませんので、添付するか、作成方法を質問文に記載してください

[PINE1103]

エラーが生じている画像1枚と、正常に1つに集約することができたpngの画像を2枚添付させていただきました。
teratailの性質上tiff画像を挿入することができなかったのでpngの画像を添付させていただきました。

[yuki23]

geotiff は単なる画像ファイルではなく座標データなどが入った複雑なファイルですので、pngに変換してしまっては意味がないです

[PINE1103]

お返事ありがとうございます、
そうですよね...
tiff画像を添付することができないので、もう少し自分で問題を考え、かみ砕いてより詳細にご質問させていただきたいと思います。