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

Question

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行目がエラーを起こしているのですが，原因が分かりません。

Answer 1

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

Comments

[py-python]

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

[can110]

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

[py-python]

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