'PrimaryHDU' object is not subscriptable の対処法

Question

pythonでfits形式の画像とpng形式の画像の同座標のピクセル値を掛け合わせてリストに入れたいのですが
以下のようなエラーメッセージが表示されました。どのように対処すればよいでしょうか。
----エラーメッセージ--------------------------------------------------------------------------------

img3[i][j] = img1_arr[i][j] * img2_arr[i][j]
TypeError: 'PrimaryHDU' object is not subscriptable

---

画像のサイズはすでに合わせており、問題ないと思っていたのですが掛け合わせるところに問題があるのでしょうか。
ご教授よろしくお願いします。

import cv2
import glob
import numpy as np
import pyfits


cut_name_num=33    #eventxxxx/の部分が後の処理で邪魔なため文字列からカットするためのもの
img3 = []

def multi_pic(filename1,filename2):
    print(filename1[cut_name_num:])
    img1 = pyfits.open(filename1)
    img2 =  cv2.imread(filename2,cv2.IMREAD_GRAYSCALE)
    img1_arr=img1[0].data
    #print(img1[0].header) #画像情報
    img1.close()
    print(type(img1_arr[0][0]))
    img2_arr = np.asarray(img2).astype(np.float32) #型を合わせる
    #numpyで扱える配列をつくる
    print(type(img2_arr[0][0]))
    print(len(img1_arr),len(img1_arr))
    for i in range (len(img1_arr)):
      for j in range(len(img1_arr)):
        img3[i][j] = img1_arr[i][j] * img2_arr[i][j]  #### ここがエラー該当箇所
    print(img3)
    


    

if __name__ == "__main__":
    event_s = 0    #event0000から数えてevent_s+1枚目から
    event_e = 1    #event_e枚目まで読み込む
    path1= XXXXX #fits
    path2= YYYYY #png
    files1 = glob.glob(path1)[event_s:event_e]   
    files2 = glob.glob(path2)[event_s:event_e]   
    #files1 = glob.glob(path1)    #一括で行う場合はこれを使う
    #files2 = glob.glob(path2)    #一括で行う場合はこれを使う
    for name1,name2 in zip(files1,files2):
        multi_pic(name1,name2)   
    

Comments

[quickquip]

エラーとコードが合っていません。(エラーに載っている行がコードのどこにもない)
記載ミスなら編集すればいい話なのですが、そうでないなら実行しているファイルと見ているファイルに食い違いがないかを確認してみては?

[_sfgh3k]

書き間違えていましたので修正しました。

申し訳ありません。

[bsdfan]

img3 = []
で初期化されているので、img3[0][0]の以前に、img3[0]の段階でIndexErrorが発生するはず。

エラーがでたときはは、img3[0]になんかのPrimaryHDUクラスのオブジェクトが入っていそうです。
jupyterでやっているとかなら、一度kernelをリスタートしてエラーを確認しなおしてみては。

[_sfgh3k]

ありがとうございます。
試してみます。

Answer 1

Numpyってarr[y, x]って書きますよね。だから

img3[i, j] = img1_arr[i, j] * img2_arr[i, j]
```では？

Comments

[_sfgh3k]

試しましたがiもjも同じ数値までしかループしないので変わりませんでした。