その２の続きです。
チコ２星表は０～２等級と一部３等級の星が含まれていません。
とりあえずの手当てなのですが、ヒッパルコス星表にある３等級未満の明るい星全部を抽出して、チコ２星表に無いもの８５恒星を見比べながら追加しました。
そして２星間の離角表を作り直し（泣；；）
3000件ほど増えて3,296,000件程度になりました。
自宅で撮影する星空は高度１０～４０度と決まっていて１０～２０度は大気差の影響を受けることがわかってきました。縦画角１．１度で１．３秒/ピクセル、４秒の許容誤差だと大気差によって３ピクセル以上恒星がずれていることになります。
なので安直なやりかたですが、最初から高度２０度の写真を基準に、縦が縮んでいるという画角調整パラメータを導入しました。
許容誤差４秒を大きくしてしまうと見つかる三角形が増えすぎて目的を達成できないからです。
前回の２つの問題点を回避しました。
１．大気差 → 画角調整パラメータの導入（デフォルト高度２０度）
２．チコ２星表 → ヒッパルコスから明るい恒星を外挿（本来チコ１＋チコ２－カブリかも）
効果を写真でどうぞ
元画像で下記の3点の座標をプログラムに入力。
実行結果は２６件の三角形を発見。
そのうち赤緯-57（カノープス）以下に含まれる候補はカット（うちから見えませんから）。
４等級未満の明るい星を含まない候補はカット。
残ったのが2件。
1件目
2件目
1件目はおおいぬ座の後足の先にある Furud という３等星で近くの8等星と共に鏡像の三角形を作っています。
というわけで２件目のレグルスが正解となりました。
ソートやバイナリサーチなどを適宜使っているので実行時間は10秒未満です（２GHｚのPC）。
たぶん写真を90度回すと見つからないんだろうな・・・。
もう少しパラメータの調整や候補の絞りこみが自動化できないかやってみます。
チコ２星表は０～２等級と一部３等級の星が含まれていません。
とりあえずの手当てなのですが、ヒッパルコス星表にある３等級未満の明るい星全部を抽出して、チコ２星表に無いもの８５恒星を見比べながら追加しました。
そして２星間の離角表を作り直し（泣；；）
3000件ほど増えて3,296,000件程度になりました。
自宅で撮影する星空は高度１０～４０度と決まっていて１０～２０度は大気差の影響を受けることがわかってきました。縦画角１．１度で１．３秒/ピクセル、４秒の許容誤差だと大気差によって３ピクセル以上恒星がずれていることになります。
なので安直なやりかたですが、最初から高度２０度の写真を基準に、縦が縮んでいるという画角調整パラメータを導入しました。
許容誤差４秒を大きくしてしまうと見つかる三角形が増えすぎて目的を達成できないからです。
前回の２つの問題点を回避しました。
１．大気差 → 画角調整パラメータの導入（デフォルト高度２０度）
２．チコ２星表 → ヒッパルコスから明るい恒星を外挿（本来チコ１＋チコ２－カブリかも）
効果を写真でどうぞ
元画像で下記の3点の座標をプログラムに入力。
実行結果は２６件の三角形を発見。
そのうち赤緯-57（カノープス）以下に含まれる候補はカット（うちから見えませんから）。
４等級未満の明るい星を含まない候補はカット。
残ったのが2件。
1件目
TYC 7073-1071-1,,95.20901917,-30.63654028,9.7560000000000002,7.8099999999999996
TYC 7073-2425-1,,95.05722028,-30.01853833,9.2119999999999997,7.8109999999999999
TYC 7073-2426-1,30122,95.07828028,-30.06337694,2.7829999999999999,2.9609999999999999
2件目
TYC 833-134-1,,152.05388583,11.99694,9.2590000000000003,8.2430000000000003
TYC 833-847-1,,152.67269556,11.82495417,9.2769999999999992,8.4030000000000005
TYC 000,49669,152.093583,11.967194,1.360000,1.360000
1件目はおおいぬ座の後足の先にある Furud という３等星で近くの8等星と共に鏡像の三角形を作っています。
というわけで２件目のレグルスが正解となりました。
ソートやバイナリサーチなどを適宜使っているので実行時間は10秒未満です（２GHｚのPC）。
たぶん写真を90度回すと見つからないんだろうな・・・。
もう少しパラメータの調整や候補の絞りこみが自動化できないかやってみます。