我有一个Fits文件,其中包含12种不同的图像信息。我可以通过使用该行找到我想要查看的数组的最大值MAG_AUTO。
a=pyfits.getdata(data1).MAG_AUTO
其中data1是来自fit文件的数据集。但是我在阵列中的位置之后发生了最大值。我怎么会发现这个?
答案 0 :(得分:1)
我想出了答案,以下命令为给定文件设置了最大值
for arg in sys.argv[1:]:
a=pyfits.getdata(arg).MAG_AUTO
arr=numpy.array(a)
indices = heapq.nlargest(10,xrange(len(arr)),key=arr.__getitem__)
print indices