質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

ただいまの
回答率

90.45%

  • Python 2.7

    1473questions

    Python 2.7は2.xシリーズでは最後のメジャーバージョンです。Python3.1にある機能の多くが含まれています。

  • NumPy

    666questions

    NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

  • GIS

    13questions

    GISとは、地理情報システムの略称です。位置・空間などの様々なデータをコンピュータを使用して加工・管理することで、情報の分析や解析を行ったり、視覚的に表示します。行政や市民生活、ビジネスなどで利用されており、活用範囲が広がっています。

【Geotiff】フォルダ内の画像ファイルを全て読み込んで,合計を出す

解決済

回答 1

投稿

  • 評価
  • クリップ 0
  • VIEW 1,006

py-python

score 8

Geotiffを利用して解析を行ってるものです。
フォルダ内にあるGeotiffを読み込んで,全ての合計値を出す作業をしています。
avetest.tifを初期値として,残りの*.tifを読み込んで,積算するプログラムを作成したのですが,
2回目以降にエラーが起きてしまいます。原因を調べているのですが,分かりません。
どこの部分がおかしいのか教えてください。

import sys
import numpy
import numpy as np
import subprocess
from osgeo import gdal, gdalconst, gdal_array, osr
import matplotlib.pyplot as plt
import matplotlib.cm as cm
import glob
from PIL import Image
import PIL

path = "/media/student/data1"
path1 = "/media/student/data1/test/"

datalist = glob.glob(path1 + "*.tif")
datalist1 = glob.glob(path1 + "avetest.tif")
for data1 in datalist1:
    for data in datalist:
        name1 = data[len(path1)+0:-4]
        with open(data, 'r') as target_file:
                if target_file:
                                print name1+ "OK"
                                data = gdal.Open(data, gdalconst.GA_ReadOnly) # tif (read only)
                                data1= gdal.Open(data1, gdalconst.GA_ReadOnly) # tif (read only)
                                add = data.GetRasterBand(1).ReadAsArray() # Band1 numpy array
                           b1 = data1.GetRasterBand(1).ReadAsArray() #Band1numpyarray
                                sum_v = b1 + add #sum
                                pil_img = Image.fromarray(sum_v)
                                print "sum...ok"
                                outdata = path1+ "avetest.tif"
                                pil_img.save(outdata)
                                outdata = gdal.GetDriverByName('GTiff').Create(path1 + "avetest.tif" , 2500, 1500, 1, gdal.GDT_Float32) # empty_output path,xsize,ysize,band,dtype
                                xmin = 130 #latitude
                                ymax = 30 #longitude
                                resx = 0.02 #grid
                                resy = 0.02 #grid
                                outdata.SetGeoTransform([xmin, resx, 0, ymax, 0, -resy]) 
                                srs = osr.SpatialReference() 
                                srs.ImportFromEPSG(4326) # WGS84
                                outdata.SetProjection(srs.ExportToWkt()) # Combine spatial information
                                outdata.GetRasterBand(1).WriteArray(sum_v)   # b1_numpy_2matrix
                                outdata.FlushCache()                     # read file
                                outdata = None         
  • エラー部分
student@student:/media/student/data1$ python a_sum.py
20161001-0010OK
sum...ok
20161001-0020OK
Traceback (most recent call last):
  File "a_sum.py", line 27, in <module>
    data1= gdal.Open(data1, gdalconst.GA_ReadOnly) # tif (read only)
  File "/usr/lib/python2.7/dist-packages/osgeo/gdal.py", line 2943, in Open
    return _gdal.Open(*args)
RuntimeError: not a string


どうやら,2回目の27行目がエラーを起こしているのですが,原因が分かりません。

  • 気になる質問をクリップする

    クリップした質問は、後からいつでもマイページで確認できます。

    またクリップした質問に回答があった際、通知やメールを受け取ることができます。

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

回答 1

checkベストアンサー

0

for data1 in datalist1:ループ初回ではdata1はファイルパス文字列ですが
data1= gdal.Open(data1, gdalconst.GA_ReadOnly)で同じdata1GDALオブジェクト?が代入されています。
よってfor data in datalist:の2回目以降のループで提示エラーが発生していると思われます。

投稿

編集

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2017/11/14 16:26

    ありがとうございます。
    お手数ですが,どう書き換えればいいでしょうか?

    キャンセル

  • 2017/11/14 16:37

    data1= gdal.Open(data1, gdalconst.GA_ReadOnly) の中のdata1のうち、どちらかを別(の名前)の変数とすればよいです。
    また、たとえばdatalist,data,datalist1,data1などが何を表しているか変数名では分かりずらいので
    分かりやすい変数名をつけるように心がけた方がよいです。
    たとえば、data1 は file_path_base(初期値、基準となるデータファイルパス)など。

    キャンセル

  • 2017/11/15 11:39

    出来ました!ありがとうございます!

    キャンセル

15分調べてもわからないことは、teratailで質問しよう!

  • ただいまの回答率 90.45%
  • 質問をまとめることで、思考を整理して素早く解決
  • テンプレート機能で、簡単に質問をまとめられる

同じタグがついた質問を見る

  • Python 2.7

    1473questions

    Python 2.7は2.xシリーズでは最後のメジャーバージョンです。Python3.1にある機能の多くが含まれています。

  • NumPy

    666questions

    NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

  • GIS

    13questions

    GISとは、地理情報システムの略称です。位置・空間などの様々なデータをコンピュータを使用して加工・管理することで、情報の分析や解析を行ったり、視覚的に表示します。行政や市民生活、ビジネスなどで利用されており、活用範囲が広がっています。