写真から撮影方向を分析する－その２

前回の記事　写真から撮影方向を分析する－見えてきた問題点　の反省点を踏まえ、新たにデータ作成を行っています。

今回使ってみたのがヒッパルコス星表に代わるチコ２星表。
なんと約２５４万の星を含んでいます。

９等級までの星を抽出して約21万7千。
それらの距離が１．２度未満の組み合わせを1時間半かかって329万件作成。

インタープリタ言語では遅いのでC言語で書き直し、だいたい10秒程度で三角形を見つけられるようにはなったのですが、２つ問題が出てきました。

１つは大気差。
低空の星を取った場合、縦方向が縮んでしまい、正しい三角形が見つけられない場合があります。星の位置にプラスマイナス4秒の誤差を見込んでいますが、当然大気差のほうが大きいです。

2つ目はチコ２星表。
たったいまわかったのですが、１等星が含まれてません。
レグルスと8等星の作る三角形がどうしても見つけられず、違う場所を見つけてしまいます。
チコ１星表と排他関係なのかも知れないのでこれから調査します。

もし排他関係なら組み合わせ表の再作成からなのでめげます；；

対象データが大きいので sort や grep を駆使しているもののテストが大変です。