くそー

■u'g'r'i'z' to BVRcIc
たとえば(prime 省略) Subaru system の R 等級を R、SDSS system をつかて作った R 等級を R'(r,i) と書く。
R = R'(r,i) ~ r - 0.2936*(r - i) - 0.1439
http://www.sdss.org/dr6/algorithms/sdssUBVRITransform.html
でもこれが違うから困る。

R'(r,i)~r-c1*(r-i)-c2
とする。
本来 R-R' = 0 のはずなので
<R-R'> = c2 -> 0 なるパラメータセットをとればよい。
ついでに言うと、c2 = 0 とすれば
<(R-R')^2>-> 0
でもあるはず。したがって、c1 を適当にふって一番 <(R-R')^2> が小さくなるときの値を選べばよい。

■そもそも
Rc だけ check してみるべし。Standard Star で。

■iraf
imcombine の sigclip ってどうやってるんだろう？