のんびり星あそび

  ３部作最後の加工編になります。 １部手動前処理編 ２部スクリプト編 スクリプト編の終盤にコメントしましたが、加工に関しては勉強中です。 特にストレッチが難しくて、ここ最近はGHSを練習していますが何が正解なのか正直よく分かりません… まだまだ試行錯誤中なので、 tomが使えている（ = 素人でも比較的簡単に使える）ツールに重点を絞って 進めていきます。 1.4.0の新機能を使うため、SPCCとGraXpertを使用する形で進めます。 追記：Siril 1.4.0 リリースに伴って内容を一部修正しています。 主にUIの変更や仕様変更で生じた操作法の変更を反映し ています。 GraXpertのインストール siril 1.4.0の新機能のため、GraXpertをインストールしておきます。 記事作成時点ではGraXpert 3.1.0 rc2です。 rc版なので、ダウンロードは github になります。 このバージョンからAIデコンボリューションが使えます。 ただし、リリース候補版なのでどんな不具合があるか分かりません。 安定版を使いたい場合は、 GraXpert のHPからダウンロードします。 Sirilでの設定は ゴンカネさんのブログ に詳しく載っています。 リニア処理 スタッキングが終わった画像を読み込みます。 画面はほとんど真っ黒のままです。 右下にある「線形」を「自動イコライゼーション」に変えると撮影したものが見えます。 画像の内、必要な所をトリミングします。 （もちろん、画面全体を使う場合は必要ありません） ここでGraXpert-AI.pyを起動させます。 Background Extractionを選択します。 Smoothingは0.2程度で良いと思います。 applyを押せば補正されます。 イマイチであれば戻るボタンで処理を取り消して数値を変更しやり直します。 続いてSPCCを実行します。 siril 1.4.0から導入されたSPCCは、画像の座標データが必須になります。 まず ツールから「天体測量」をクリックし、その先の「アストロメトリ」をクリックします。 まずは望遠鏡の焦点距離をFocal lengthに、カメラの１ピクセルのサイズをピクセルサイズに入力します。 画像パラメータに撮影した天体の名称を入れ、findをクリックします。 す...

  前の投稿 でスケアリングが問題になったので、どれくらい傾いているのか調査から始めます。 何も処理していないraw画像をディベイヤーしてfitsファイルに変換します。 この画像を読み込んで Tool → image analysis → show tilt で下の画像が出てきます。 コンソールに表示される数値の意味は次の通りです。 stars：検出された星の数 Truncated mean：センサー全体で検出された星のFWHMの平均値です。外れ値を取り除くため、トランケート平均が使われています。 Sensor tilt：センサーの傾き。画像四隅のfwhmの最良値と最悪値の差です。括弧内には差の割合が表示されており、これが 10%以上だと調整が奨励 されます。 Off-axis aberration：中心部と周辺部のfwhmの差。コマ収差や像面歪曲を表す指標です。 分かっていたことですが、猛烈に傾いていますね。 調整が奨励される 閾値の ３倍 ですか・・・ 像面歪曲は仕方がありません。F5 格安アクロマートですから。 では、 30%を超えている傾きを10%未満になるまで頑張ります 。 と言いたいところですが、星が出てないと直せないのは厳しいです。 星が出てたら撮影したいのに、素人がスケアリング調整してたら撮影時間無くなっちゃうよ… という事で、星が出てなくてもスケアリングする方法を探りたかったのですが… tomの家の中ではピントを合わせられる直線距離がありません。 集合住宅だと、撮影時だけでなくこういう所も制限が多くなっちゃいますね。 どーすっぺ。 悩んでいる間に、薄曇りながら少しだけ星が見えたのでスケアリング調整にトライしました。 狙いとしては、Tリングマウントアダプターにアルミホイルを挟んでみようと思っています。 さて撮影。 お！ 15%まで改善した。 わりと良い感じ？と思ったのも束の間。 ここから迷走が始まります・・・ 調整したアクションと、チルトの結果が全く一致しないんです。 四苦八苦した結果、アルミホイルは思うような効果が出てない事が分かりました。 というのも、 他に支配的な因子が判明 したからです。 ネジを緩めて締め直せば傾きが変わる という、ごく初歩的な問題でした。 そりゃアルミホイルの厚みの約12μmよりも、ブレ幅は大きいでしょうな。 つまり...

タイトルの通り、スタック時の重み付けオプションが気になったので調べてみます。 もうひとつ、 前回の投稿 の後半で -equlize_cfa の有無では差異無しという結果でしたが、等高線グラフにして僅かな差も見えるようにしてみました。 （ちょっとは差があって欲しいという願望で追加調査です） スタック時の重み付け tomはスタック時の重み付けについて -weight=fwhm を使っています（ 備忘録 Siril 1.4.0 スクリプト編 参照）。 これはこれで良いとは思っているのですが、他のオプションとの違いは何だろう？と疑問に思いました。 対象によって向き不向きがあるかもしれません。 という事で、確認していきます。 オプションは次の３種類で、ざっくりした意味です。 noise：背景ノイズの値によって重み付けを行う nbstack：露光時間によって重み付けを行う nbstar or fwhm：fwhmに基づいて重み付けを行う。fwhmが同じ場合は星の数が多い方が良いと判断される。 tomは露光時間に差をつけて撮影していないので、nbstackは意味がなさそうです。 いつもはfwhmで処理しているので、今回はnoiseと比較してみます。 スタック後にノイズ評価と動的SPFで比較します。 fwhmで重みづけ noiseで重みづけ 確かに、若干ですがnoiseで重み付けした方はノイズが少なく、fwhmで重み付けした方がfwhmが小さいです。 ですが、これがどれくらいの差なのか素人には分かりません。 数値上の話だけかもしれないので、実際に画像処理してみます。 左がfwhm、右がnoiseです。 画像だと何にも変わらないです。 んー、オメガ星雲は全体に何かしら写っているので、バックグラウンドノイズの判定が難しいかも？ という事で、３月撮影のオリオン大星雲でも処理してみます。 左がfwhm、右がnoiseです。 変わらない… 数字上は僅かながら変化はありますが、私の撮影した画像では差が出ませんでした。 星団とかを写すときはfwhm、星雲や銀河の時はnoiseにするべきなのかなぁと思います。 が、差は大きくないと言うことですね。 -equlize_cfaの追加調査問題 前回の投稿で、-equlize_cfaの有無ではパっと見た感じでは変化がありませんでした。 有料ソフトでは等高線図...

 ようやく梅雨が明け、待ちに待った撮影が6/28に出来ました。 結果から先にいきます。 機材：SE102, ZV-E10未改造, QBP, AZ-GTi経緯台モード 条件：ISO1250, 25s露光×298枚 処理：トリミング、SPCC、AI デノイズ&デコンボリューションあり 撮影を始めてから７天体目。 初めて露光時間が２時間を超えました。 撮影よりも処理が難しかったです・・・ tomとしては満足のいく仕上がりになりました。 今回の撮影で試そうと思っていた事が２つあります。 ①撮影時にフラットを撮る 今まで撮影後に ダーク を撮って、 フラット は別日に撮っていました。 今回は撮影後に フラット を撮って、 ダーク を別日に撮りました。 ダークの温度合わせは出来ませんが、光学系は完全に同一なのでフラットは精度が上がるはずです。 どちらの方が良いのか調べていきます。 ②フラット補正時に、-equalize_cfa の有無で差異の確認 tomはQBPフィルターを用いているので、撮影した映像もフラットもRGBが偏っています。 このRGBが偏ったマスターフラットで補正をする場合、 -equlize_cfa でマスターフラットのRGB強度を強制的に均一化させる事で何かしらの影響があるかを確認します。 では①からです。 これはスタッキングが終わった瞬間に結果が分かりました。 上は過去撮影したM42、下は今回撮影したM17です。 M42は露光時間１時間、M17は２時間なので写野回転の影響はM17の方が大きいです。 が、M42にあった円形のカブリがM17ではなくなっています。 明らかにフラットを当日撮影した方が良くなっています。 フラットを当日撮影したデメリットとして、ダークの温度が合わなくなります。 これに加え、夏場の高温でダーク減算が更に合いにくくなるはずです。 …と思っていましたが、過補正や補正不足となっているピクセルが分かりません。 なんでだ…？？ 私はAZ-GTiを使っているので、写野回転が発生します。 過補正や補正不足によるノイズは写野回転に伴って同心円状に発生すると考えています。 また、家が鉄筋コンクリートのマンションなのでスマホのコンパスが狂ってしまい、真北が正確に分かりません。 追尾エラーはそこそこ大きいです。 今回の撮影でも、途中で鏡筒の向き...

それでは Siril 操作手順のスクリプト編になります。 前回 は「自動前処理編」と言いましたが、「スクリプト編」とした方が良い内容になってきたのでタイトルを変えました。 Siril スクリプトを使用方法についてのブログはそこそこありますが、スクリプトの中身を扱うブログは少ないと思います。 スクリプトは自分向けにカスタマイズすれば利便性が格段に上がりますが、GUI操作に比べてハードルは高めです。 と言うことで、簡単ですがスクリプトの解説を含めた内容にしました。 備忘録と銘打っていますが、二の足を踏んでいた人の背中を押すことが出来たら嬉しいです。 なお、本ページで記載しているスクリプトはSiril公式のスクリプトを元に、tomが非公式に改良したものです。 元スクリプトの配布元は下記のリンク先にあります。 元のスクリプトは GNU General Public License v3（GPL v3） に基づいて公開されており、自由に利用・改変・再配布が可能です。 ただし、再配布の際には同じライセンスを適用する必要がありますので、私の改良版もまたGPL v3に基づいて公開しています。 著作権は元スクリプト作者である Cyril Richard氏に帰属します。 改良後のスクリプトになります。 あくまでも「tomが使っているスクリプト」の備忘録なので、もっと良い方法があるかもしれません。 こういうやり方もあるんだ程度の軽いノリで読んで頂ければと思います。 Siril のバージョンは1.4.0 beta2です。 追記：Siril 1.4.0 までのUI変更等に対応して画像などを差し替えています。 前提は前回から太字の所を加えます。 ・ベイヤー構造のカラーカメラで撮影した画像 ・ダークフレーム、 フラットフレーム 、ライトフレームがある ・ディベイヤーではなく、bayer drizzleを使う バイアスについては、バイアスフレームではなく合成バイアスを使っています。 撮影条件によってバイアスフレームを準備するのも大変なので、条件を変えて撮影したバイアスフレームから合成バイアスを算出します。 算出についても説明します。 その後、マスターフラット、マスターダークを作成して、最後に天体写真の前処理を行います。 合成バイアスの算出 この合成バイアスを使うメリットとしては、１度算出してしまえ...

最近は撮影が全く出来ないので、Sirilについての記事ばかりになってきました。 いろいろ触って使い方にも慣れてきたので、Sirilの操作手順を備忘録的な感じでまとめたいと思います。 あくまでも「tomが行っている処理方法」の記録で、Sirilのバージョンは1.4.0 beta2です。 追記：Siril 1.4.0 でも問題ない事を確認済みです。 今のところは手動前処理編、自動前処理編、加工編の３部作になる予定です。 インストールの前に まずはインストール時のちょっとしたポイントです。 edgeでダウンロードすると下のアラートが出ますが、Sirilに関しては問題ありません。 もっと見るから保持を選びます。 単純に使っている人が少ないので「このソフト大丈夫？」と聞かれているだけです。 使用者が少ないという事は安全性の確認が出来ているか怪しいので、このアラートで一旦立ち止まらせる機能です。 びっくりしますし面倒ですが、安全には変えられません。 私は良いセキュリティ機能だと思います。 こういうアラートが出た場合は、インストールする前にセキュリティソフトでスキャンすれば、ウイルスに感染するリスクを下げられると思います。 では無事にインストールが終わったとして前処理を始めたいと思います。 今回は下記３点の前提で進めていきます。 特に1.4.0の新機能であるbayer drizzleを使うところがポイントです。 ・ベイヤー構造のカラーカメラで撮影した画像 ・ダークフレームとライトフレームがある ・ディベイヤーではなく、bayer drizzleを使う フラットとバイアスは次回の自動処理編で話をしようと思います。 マスターダーク作成 画面左上のホームボタンを押し、撮影したファイルがあるフォルダを指定します。 まず、 撮影した画像ファイルを、sirilで処理するfitsファイルに変換します 。 同時にシーケンスファイルというモノも作成されるのですが、これがSirilの前処理にとって必要です。 最初からfitsファイルの場合もシーケンスファイルを生成する為に行ってください。 画面右上の変換タブをクリックします。 「+」ボタンを押して、撮影したダークフレームを全て選択します。 シンボリックとディベイヤーのチェックを外し、シーケンス名を入力します。 今回はダークフレームなので、dark に...