muttenz's blog

スイス星空だより

観測時刻の修正

ある観測の後でコンピュータの時計をチェックしたら合わせたはずが合っておらず、 ds秒進んでいるとする。AIP4WINのログで出てくる時刻はdate(YYYY-MM-DD)とtime(hh:mm:ss.sss)のフォーマットで、ds秒ずらすのは単純に引き算ではできない。

Rでこういう風にh m sを変換するようにしてみた。

(ただし、秒の単位を四捨五入したら、32時が33時になってしまうという場合はこれではだめですが。まあ滅多にないでしょう。スイスの場合は時差の関係でこれも考えに入れないといけないけど。)

 

  S <- round(h * 3600 + m * 60 + s - ds)

  s <- S %% 60

  m <- (S %/% 60) %% 60

  h <- S %/% 3600