今日记录

1. https://reproject.readthedocs.io/en/stable/

一个常见的例子是，您有两个fits图像，并且希望将其中一个重新投影到与另一个相同的坐标中。

这可以很轻松的用reproject得到，我们来用下面的例子证明。我们先下载两个图像，

在 http://data.astropy.org下载。

reprodect总共有三个算法:

from reproject import reproject_interp
from reproject import reproject_adaptive
from reproject import reproject_exact

重投影函数有两个主要参数。第一个参数是要重投影的图像，以及有关该图像的WCS信息。这可以是：

• The name of a FITS file

• An HDUList object

• An image HDU object such as a PrimaryHDU, ImageHDU, or CompImageHDU instance

• A tuple where the first element is a ndarray and the second element is either a WCS or a Header object

 

2.http://montage.ipac.caltech.edu/docs/algorithms.html  

用montage做reproject

图像在做reproject的时候发生了什么？

图像重投影涉及从一组输入像素到一组输出像素的信息再分配。对于天文数据，输入像素代表从天空某个区域接收到的总能量，在将这些信息重新分配到输出像素时，保存这些信息至关重要。同样在天文学中，保持能量分布的位置（天体测量）精度也很重要，因此，将输入像素的所有能量添加到“最近的”输出像素等常用技术是不够的。

相反，我们必须根据这些像素的精确重叠重新分配输入像素能量到输出，甚至可能使用基于原始仪器的点扩散函数的像素加权函数。其目标是创建一个输出图像，该图像尽可能接近如果使用具有输出图像像素模式的仪器观察天空时所创建的图像。我们也致力于建立一个同样能处理所有天文投影和坐标系统的系统。

montage是怎样做图像reproject的？

 

3.https://danmoser.github.io/notes/gai_fits-imgs.html

上面的网址和fits文件的头文件和坐标有关。

fits n维图像，每一个像素的物理坐标或wcs。

地图投影的几种方式，也就是立体转为平面的方式：

 

fits头文件中的wcs坐标系统，定义了像素和天空坐标之间的关系。

投影关系总共有25个，主要有

TAN: Gnomonic = Tangent Plane  切面

CAR: Cartesian 笛卡尔

CTYPE1 and CTYPE2 indicate the coordinate type and projection

CTYEP1 = xxxxyyyy / 'Comment' CTYEP2 = xxxxyyyy / 'Comment'

The first four characters are RA-- and DEC-, GLON and GLAT, or ELON and ELAT, for equatorial, galactic, and ecliptic coordinates, respectively.

赤道，银河，黄道。

第二个四个字符yyyy包含投影的四个字符代码

CRPIX1和CRPIX2是投影和旋转所参考的参考点的像素坐标

crval1 crval2,以十进制度数表示的赤经和赤纬或经纬度的中心坐标。

CDELT1和CDELT2用来表示每像素度数的图版比例

CROTA2用于表示水平轴和垂直轴的旋转度

 4.https://www.codenong.com/21415059/

关于python：通过PyFITS / AstroPy在FITS图像中的笛卡尔投影问题