仿真畸变测试

利用W2天区的229次曝光的06号探测器（最边缘），检查畸变是否随时间变化。

方法：利用0级数据hdu1中的线性wcs，把输入星表中的xy转化为RA，Dec，然后再与输入星表中的ra_orig, dec_orig比较。为了避免自行引起的偏差，只选择了obj_type=galaxy的源进行比较。

结果：无论是RA还是Dec，(x,y)线性转换后的结果都与(ra_orig, dec_orig)存在偏差。这个偏差在空间上的分布pattern是稳定的。但偏差的中值却随指向的不同发生变化，且与指向无关。见下面两图

图中每个方块表示一个探测器图像，总共229次曝光，包括了229个06号探测器方块。方块的颜色表示该次曝光的整个探测器ra/dec偏离线性度的中值。即median(ra_linearity - ra_orig)。按照理论，如果系统稳定，且指向无偏差，偏离线性的程度不应该有几十个角秒的变化。初步判断是线性关系参考点的计算有误，但需要仿真输出来检验。

此问题影响最大的不是主巡天位置定标，而是导星。一旦每次畸变都会有系统漂移，基本上就无法准确计算导星所在的准确位置（几十角秒显然是无法接受的）。

@weichengliang @zhangxin @liudezi @luochangqing 建议专门就WCS和畸变展开一次讨论