我在python中打开了一个目录,它有大约7万行不同对象的数据(ra、dec坐标和对象名称)。我还有另外一个列表,其中有大约15,000个感兴趣的对象,它们也出现在前面提到的目录中。对于这15,000个对象中的每一个,我想看看在70000个大列表中是否有任何其他对象在对象的10弧秒内具有ra,dec坐标。如果这是真的,我只想标记这个对象并转到下一个对象。obj_ind = list
我正在尝试将FITS文件的坐标系从其原始的赤道坐标系更改为银河坐标(以度为单位),以便使用这些坐标操作生成的FITS图像。为此,我需要提取一个包含每个像素赤道位置的数组,以便将它们转换为所需的银河系坐标。这就是我的知识有限的地方,似乎找不出如何提取该数组。最后,我想以如下方式根据纬度对图像进行切片:
f = pyfits.open("im1.fits"