MontageError:无效的表文件

Question

我在我的Debian服务器上使用Python 2.7.3。以下是输入代码：

import aplpy
import pyfits
from numpy import *
import matplotlib as plt
import montage as montage_wrapper
import PIL

r = pyfits.open('3c324IR.fits')
b = pyfits.open('3c324UVIS.fits')
g = pyfits.open('3c5GHZ.fits')

r.info()
g.info()
b.info()

print(r[1].header)
print(g[1].header)
print(b[1].header)

r1 = r[1].data
g1 = g[0].data
b1 = b[1].data

hdu = pyfits.PrimaryHDU(r1)
hdulist = pyfits.HDUList([hdu])
hdulist.writeto('r.fits')

hdu = pyfits.PrimaryHDU(g1)
hdulist = pyfits.HDUList([hdu])
hdulist.writeto('g.fits')

hdu = pyfits.PrimaryHDU(b1)
hdulist = pyfits.HDUList([hdu])
hdulist.writeto('b.fits')

aplpy.make_rgb_cube(['r.fits','g.fits','b.fits'], '3c_324_rgb.fits')

aplpy.make_rgb_image('3c_324_rgb.fits','3c_324_rgb.png')

f = aplpy.FITSFigure('3c_324_rgb.fits')

f.show_rgb()
f.save('3c_324_rgb2.png')

print 'END'

在我开始代码的aplpy.make_rgb_cube()部分之前，一切都很正常。我之所以知道这一点，是因为当我启动它并将其剪切出来时，它可以正常工作(直到它意识到这段代码没有输出可供继续)。以下是终端中的混乱情况：

Traceback (most recent call last):
  File "test9.py", line 47, in <module>
    aplpy.make_rgb_cube(['r.fits','g.fits','b.fits'], '3c_324_rgb.fits')
  File "/usr/local/lib/python2.7/dist-packages/aplpy/rgb.py", line 309, in make_rgb_cube
    montage.mMakeHdr(images_raw_tbl, header_hdr, north_aligned=north, system=system, equinox=equinox)
  File "/usr/local/lib/python2.7/dist-packages/montage_wrapper/commands.py", line 1468, in mMakeHdr
    return status.parse_struct("mMakeHdr", p.stdout.read().strip())
  File "/usr/local/lib/python2.7/dist-packages/montage_wrapper/status.py", line 33, in parse_struct
    result = Struct(command, string)
  File "/usr/local/lib/python2.7/dist-packages/montage_wrapper/status.py", line 70, in __init__
    raise MontageError("%s: %s" % (command, self.msg))
montage_wrapper.status.MontageError: mMakeHdr: Invalid table file: /tmp/tmpmyYyN7/images_raw.tbl

我不知道上面的是什么意思。我在谷歌中到处寻找它，但一无所获。这个错误是什么?我如何修复它？

Answer 1

在生成立方体的过程中，aplpy.make_rgb_cube会创建几个由蒙太奇mMakeHdr使用的文件。

aplpy.make_rgb_cube跳过了mMakeHdr的所有步骤。

mMakeHdr尝试读取其中一个创建的名为/tmp/tmpmyYyN7/images_raw.tbl的文件，读取该文件失败，因为该文件无效。无效的文件可能包含错误的内容，也可能意味着该文件不存在。

1. 是否真的创建了/tmp/tmpmyYyN7 /

_raw.tbl？

1. 如果创建了，它是否包含图像文件名列表？
2. 检查创建日期，查看这是您创建的文件还是不能对/tmp和/tmp/tmpmyYyN7
3. 拥有读写权限的旧文件。如果文件不在那里，但您对临时文件夹有写权限。如果您的系统上有多个images_raw.tbl安装，请检查并查看是否在您的harddrive.
4. Check上的任意位置写入了一个名为mProjectPP的文件:从shell命令行运行“mProjectPP

在Python中，do: import os os.system(‘mProjectPP')应该会提供相同的答案。如果不是这样，您就有两个蒙太奇安装，而且它们正在混淆。

Answer 2

1. 我转到/tmp/tmpmyYyN7/，images_raw.tbl在那里。
2. 如果我输入“images_raw.tbl images_raw.tbl”，它包含nothing...does，这意味着它不包含任何图像文件名？
3. 我用"chmod 777“更改了读/写权限，但当我运行它时，它仍然不起作用，只是产生了相同的输出。
4. 输入”nothing...does images_raw.tbl“，结果如下：

根@jarvis:/tmp/tmpmyYyN7# stat images_raw.tbl文件：`images_raw.tbl‘大小:0块:0 IO块: 4096常规空文件设备: 801h/2049d Inode: 5772965链接:1访问：(0644/-rw-r--r--) Uid：( 0/根) Gid：( 0/根)访问: 2014-08-15 04:40:25.433245271 -0400Modify: 2014-08-15 04:40:25.429245244 -0400 Change: 2014-08-15 04:40:25.429245244 -0400出生：-

• 这是输出：

/home/python_lib/Montage_v3.3/Montage//mProjectPP @jarvis:/tmp/tmpmyYyN7#哪个mProjectPP root

Answer 3

以下是使用调试器的输出：

> /home/python_lib/test10.py(59)<module>()
-> aplpy.make_rgb_cube(['/tmp/tmpmyYyN7/r.fits','/tmp/tmpmyYyN7/g.fits'], '3c_324_rgb.fits')
(Pdb) n
MontageError: MontageE...aw.tbl',)
> /home/python_lib/test10.py(59)<module>()
-> aplpy.make_rgb_cube(['/tmp/tmpmyYyN7/r.fits','/tmp/tmpmyYyN7/g.fits'], '3c_324_rgb.fits')
(Pdb) n
--Return--
> /home/python_lib/test10.py(59)<module>()->None
-> aplpy.make_rgb_cube(['/tmp/tmpmyYyN7/r.fits','/tmp/tmpmyYyN7/g.fits'], '3c_324_rgb.fits')

然后，当我再次点击"n“"enter”时，它会弹出错误：

(Pdb) n
Traceback (most recent call last):
File "test10.py", line 59, in <module>
aplpy.make_rgb_cube(['/tmp/tmpmyYyN7/r.fits','/tmp/tmpmyYyN7/g.fits'], '3c_324_rgb.fits')
File "/usr/local/lib/python2.7/dist-packages/aplpy/rgb.py", line 309, in make_rgb_cube
montage.mMakeHdr(images_raw_tbl, header_hdr, north_aligned=north, system=system, equinox=equinox)
File "/usr/local/lib/python2.7/dist-packages/montage_wrapper/commands.py", line 1468, in mMakeHdr
return status.parse_struct("mMakeHdr", p.stdout.read().strip())
File "/usr/local/lib/python2.7/dist-packages/montage_wrapper/status.py", line 33, in parse_struct
result = Struct(command, string)
File "/usr/local/lib/python2.7/dist-packages/montage_wrapper/status.py", line 70, in __init__
raise MontageError("%s: %s" % (command, self.msg))
montage_wrapper.status.MontageError: mMakeHdr: Invalid table file: /tmp/tmpQILi76/images_raw.tbl
Exception AttributeError: "'NoneType' object has no attribute 'path'" in <function _remove at 0x7f7ad03c2668> ignored

我会检查它，看看我是否能理解它……我插入了

pdb.set_trace()

就在前面的一行

aplpy.make_rgb_cube

果然，这就是错误所在。有没有可能是因为我在前面的代码中对图像进行了“切片”？我遇到的.fits镜像的问题是，每个镜像都有多个“镜像”文件：

Filename: 3c324IR.fits
No.    Name         Type      Cards   Dimensions   Format
0    PRIMARY     PrimaryHDU     748   ()
1    SCI         ImageHDU       144   (1091, 966)   float32
2    WHT         ImageHDU        48   (1091, 966)   float32
3    CTX         ImageHDU        43   (1091, 966)   int32
4    HDRTAB      BinTableHDU    557   2R x 274C    [9A, 3A, J, D, D, D, D, D, D, D, D, D, D, D, D, J, 6A, 9A, 7A, 18A, D, D, D, D, D, D, D, D, D, D, D, D, D, D, D, D, D, D, J, J, D, D, D, D, D, D, J, J, 8A, 23A, 11A, 18A, 4A, D, D, J, J, D, D, D, D, 23A, D, D, D, D, J, J, D, 23A, 8A, J, D, D, D, 23A, 1A, D, D, D, D, D, D, 12A, 12A, 8A, 23A, D, D, 10A, 10A, D, D, D, 2A, 3A, 3A, 4A, 8A, 7A, D, J, D, 6A, 9A, D, D, D, 4A, 18A, 3A, J, 5A, D, D, D, 8A, D, D, D, D, D, 23A, 1A, D, 23A, D, D, D, 3A, J, 1A, 4A, D, 3A, 6A, D, D, D, D, D, 23A, D, D, D, D, D, 1A, J, J, J, J, 8A, 23A, J, J, 10A, 7A, D, D, D, D, D, D, D, D, D, D, D, D, D, D, D, D, D, 12A, D, 24A, 23A, D, 2A, 8A, D, J, D, D, 6A, 4A, D, D, D, D, D, D, D, D, D, D, D, 19A, D, J, 23A, D, 8A, D, D, 13A, D, D, J, D, D, 5A, 9A, 1A, 4A, D, 1A, J, 1A, 1A, 1A, 1A, D, D, D, D, D, D, 4A, D, D, 4A, 3A, 1A, J, 3A, 6A, J, D, D, D, 4A, J, D, D, D, J, J, J, 1A, D, D, D, D, J, J, 6A, J, 8A, D, D, 1A, 5A, J, 3A, 8A, 8A, D, J, D, 8A, 8A]
5    WCSCORR     BinTableHDU     59   7R x 24C     [40A, I, 1A, 24A, 24A, 24A, 24A, D, D, D, D, D, D, D, D, 24A, 24A, D, D, D, D, J, 40A, 128A]

上面是我打印哈勃红外.fits头时的输出，它显示有几个HDU，但我唯一想要的是#1，这是场的图像。UVIS和VLA图像也是如此。我能想到的选择其中之一的唯一方法是用PyFITS打开它们，然后用如下所示的行选择每一个

r1 = r[1].data

这是正确的方法吗？我猜大多数.fits镜像都没有将多镜像设置到一个文件中。当我尝试以另一种方式运行它时，打开原始文件，尽管它肯定不喜欢它。

在使用上面的行进行选择后，我只需将其单独保存为.fits图像。

布兰登·多伊尔