python - FITS 画像配列のピクセルが範囲外です

Question

以下のコードを使用して、FITS ファイル内のピクセルのしきい値を設定しようとしています。ただし、次のようなエラーが表示されます。

IndexError: Index 2620 is out of bounds for axis 0 with size 2620

これを修正する方法についてのアイデアはありますか?

これはコードです：

from astropy.io import fits
import numpy as np
hdulist = fits.open("12ratio.fits")
origImData = hdulist[0].data
newImData = origImData*0

for x, y in np.nditer(origImData.shape):
        curPixel = origImData[x, y]
        if curPixel > 0.28 or curPixel < 3.11:
                newImData[x, y] = curPixel
        else:
                newImData[x, y] = 0
newhdu = fits.PrimaryHDU(newImData)
newhdulist = fits.HDUList([newhdu])
newhdulist.writeto('modifiedratio12.fits')