のんびり天体趣味

今回もあぷらなーとさんの処理ロジックをSirilで実行してみます。 試すのは ハイパーイーブンオッド 法 です！ ロジックの詳細はあぷらなーとさんの記事で。 人生初の星ぐる写真 今年の３月に天体撮影を始めたので、星ぐる写真も撮ったことがありません。 という事でまずは試しに撮影してみました。 機材：ZV-E10、SELP1650 焦点距離16mm F6.5 まぁなんというか、とりあえずロジックを試せれば良いやという感じでテキトーに撮りました。 光害は甚大ですが、案外映るもんですね笑 ただ、Siril 1.4.0-beta3で比較明合成が出来ない。 memory allocation error が発生して比較明合成が動きません。 メモリの割当量を空きメモリの半分くらいに設定してもダメ。 こんな基本的な処理がなんで出来ないんだ！と不思議に思っていたのですが… すっかり忘れていましいたが、 今使っているのはbeta版 だった！ という事でSiril 1.2.6 Portable版をダウンロードして上の処理を行いました。 いやー、何の問題もなく色々処理してくれるのでbeta3だという事を完全に忘れていました。 そりゃ不具合もあるよね。 まずは手動でイーブンオッドを試す 今回の処理はSirilの機能で実現可能だと思っているので、ひとまず手動で処理してみます。 順序は下の通り。 ①撮影データをイーブン群とオッド群に分ける ②それぞれの群を比較明合成する ③比較明合成したものを加算合成する ④最大輝度値データと比較暗合成する ①以外は全部Sirilの標準機能で出来ると思うので、ひとまず手動でやってみます。 左が普通の比較明合成、左がイーブンオッド法です（≠ハイパー） 輝度の低い星が途切れる現象が見事に改善できています。 ひとまずイーブンオッドはクリア。 が、ここから「ハイパー」にするために比較暗合成で処理が必要ですが、できません。 1.4.0-beta3でも、1.2.6でも、実行した瞬間に Sirilがクラッシュ してしまいます。 これはダメだ… pythonまで使わなくても単純なスクリプトで実行可能なハズだったのに、結局pythonスクリプトで処理しなければならないという事になりました。 スクリプト作成 先日、chatGPTが5にアップデートしたという記事を見たので...

 前回、あぷらなーとさんのノイズ補正ロジック３手法をsiril環境で実行する挑戦に成功しました。 ノイズの多い画像で絶大な効果が確認出来ています。 この最後の方に「私のカメラZV-E10はピクセルマッピングが丁寧に行われていて、ノイズ補正スクリプトの出番が少ない」と書きました。 もうひとつ、追尾エラーが大きすぎて天然ディザリングになっている可能性もあり、余計にノイズが目立たないと思っています。 が、よく考えれば ノイズが少ない上に 散っているから 気付かないだけ であって、存在している はず です。 ということは、その後の デノイズやデコンボリューションへの影響が少なからず出るはずでは？ と思いました。 人の目で見て気にならないレベルでも、こういった計算処理において 無視できない可能性 があります。 ということで、ノイズがあまり目立っていない画像に対してどこまで効果があるのかを調べてみます。 今までの処理方法とピクセルマッピング&クールファイル補正法を追加したものを比べてみたいと思います。 なお、いつも使っているDrizzleはMADクリッピングしないと色がおかしくなるので、コスミカットは使っていません。 M8 こないだSEL55210で撮影しましたが、SE102でも撮影しました。 機材：SE102, ZV-E10未改造, QBP, AZ-GTi経緯台モード 条件：ISO1250, 30s露光×270枚 処理：トリミング、SPCC、AI デノイズ&デコンボリューション (object) あり そしてこっちがピクセルマッピング&クールファイル補正版です。 若干の色味に差があるものの、大きな差はありません。 周辺の淡い部分を拡大してみると… 左側今までの処理、右側がピクセルマッピング&クールファイル補正版です。 若干ノイズ感が改善している様に見えます。 差は大きくないものの、想定通り多少は改善が見られるということかな？ M20 続いてM20です。 これも、青色が出ると良いなと期待して再度撮影しました。 機材：SE102, ZV-E10未改造, QBP, AZ-GTi経緯台モード 条件：ISO1250, 20s露光×426枚 処理：トリミング、SPCC、AI デノイズ&デコンボリューション (object)あり これをピ...

ここ最近、あぷらなーとさん考案のソフトウェアピクセルマッピング、クールファイル補正法、コスミカット法をSiril環境で実行する事に挑戦していました。 先日、最後のコスミカット法が完成し 挑戦に成功 しました！ これで私のSirilがちょっと邪悪になりました！笑 ３手法全てにおいてSiril標準の機能では実現不能で、pythonを利用して初めて実現可能になります。 なぜPythonかというと、SirilでpythonスクリプトがサポートされたのでSirilに無い機能を自由に作ることが出来るようになったからです。 ちなみに、tomはプログラミングがほぼ出来ません。 本業で必要に駆られてVBAを軽くかじっただけです。 （ループ処理がやっと書ける程度の万年初心者です） なのでコードは全てAIで書いています。 主にGeminiで、たまにCopilotを使っています。 プログラミングが出来ない人でも、AIを使う事でやりたい処理のスクリプトを作り上げることが出来ます！ では、それぞれどのようにして実現したかの話に移ります。 ソフトウェアピクセルマッピング ロジックの詳細はあぷらなーとさんの記事を参照してください。 邪崇帝主とステライメージ10にしか無い機能だと思うので、当然Sirilで出来るものではありません。 ピクセルマッピングを行うアプリを作る必要があります。 なお、出来上がったマップデータをSirilの処理の時に使用するだけので、Sirilのpythonスクリプトを使う必要はありません。 単独のアプリとして作成する方が楽です。 処理の流れは次の通りです。 ①ダークフレームを取り込んで各ピクセル毎にメジアンと標準偏差を算出し、散布図にする。 ②散布図からホットピクセルの所を選択して、この座標データをSirilで読み込む形に変換する。 Sirilで読み込む形式は documentation の通りです。 ①は指示次第ですが意外とすぐ作れました。 処理に必要なライブラリはGeminiが勝手に選定してくれます。 ②は指示の難易度が上がるのと、処理が重くなりやすいので出来るだけ簡単なGUIを検討した方が良いです。 なんせ、この散布図の中には画素数と同じ数のドットがプロットされているので下手に組むとフリーズしまくります。 私はメジアンの上限下限のパーセンタイル指定と、2点のXY座標...

 またしても綺麗に晴れたので、撮影が捗ります。 M8とM20 少し前に工作してAZ-GTiにカメラを載せられるようにしたので、SEL55210で撮影してみました。 早速リザルトです。 機材：SEL55210(焦点距離210mm), ZV-E10未改造, QBP, AZ-GTi経緯台モード 条件：ISO1250, 30s露光×245枚 処理：トリミング、SPCC、AI デノイズ&デコンボリューション ダブルズームキットのレンズとカメラ（未改造）、かつ光害地のベランダ撮影 でここまで写れば文句はないです。 望遠鏡を買わなくてもここまで写せる のは少々驚きました。 エントリークラスの機種でもダブルズームキットを持ってる人なら経緯台とQBPフィルターを買うだけで星雲が写せる！ 周囲に光害が少なければQBPすら不要ですし。 天体撮影って思ってたより ハードルが低かった 。 望遠鏡が必須という思い込みがありましたが、むしろ経緯台さえあれば何とかなる事が分かりました。 ちなみに、QBPを使っているので三裂星雲の青色が苦しいと思っていましたが、ギリギリ何とか写りました。ぼんやりしていて青色と言えるレベルでは無いですが笑 また、エントリークラスの未改造カメラでこのレベルですから、天体用カメラを買うとどこまで綺麗になるのか楽しみです。 ただ、三裂星雲の右上に違和感があるので調べてみたところ、円形のゴースト(？)が出てしまいました。 スタックしただけの画像を思いっきり強調しています。 レンズの先にQBPフィルターを取り付けたので、レンズとフィルターが近すぎたせいかな？と妄想しています。 解決策が無いか、のんびり考えたいと思います。 これを撮影していて、200mm前後の鏡筒がたくさん出ている理由がちょっと分かりました。 複数の星雲が１枚に写るのってすごく良い。 やっぱりいろんな焦点距離の鏡筒が欲しくなってしまうんですね… 天文が 沼 と呼ばれる理由もホントによく分かりました。 pythonスクリプト 撮影ではない方のネタに移ります。 ここ最近は、あぷらなーとさん考案の処理法をSirilに実装しようと頑張っています。 前回悩んでいたseqcosme（cosmeのループ処理）を使用する方法からは一旦離れて、別の方法をとりました。 これでピクセルマッピングとクールファイル補正...

  土曜の夜、風呂に入る前に外を見ると完全に曇っていて、撮影を諦めました。 が、風呂から上がってしばらくして、ふと空を見ると綺麗に晴れてる！ という事で急遽撮影しました。 狙いは前回ゴミで酷いことになったM16です。 早速リザルトです。 機材：SE102, ZV-E10未改造, QBP, AZ-GTi経緯台モード 条件：ISO1250, 30s露光×179枚 処理：トリミング、SPCC、AI デノイズ&デコンボリューションあり 処理ですが、GraXpert の deconvolution（stellar）は使っていません。 というのも、stellar で処理すると画像の外周ばかりが強烈に補正されてしまいます。 M17の時も悩みましたが、今回は避けようがないくらいに酷かったので採用出来ませんでした。 どうすれば外周だけでなく全体的に補正がかかるのか… 何か良い条件が無いか探ってみたいと思います。 GraXpertの不具合なら逆にありがたいくらい。 また、以前効果が無いと言った deconvolution の object ですが、めっちゃ効果ありました。 創造の柱がクッキリしました。 左がdeconolution(object)なし、右がありです。 意味無いと言ってスイマセン…これからはobjectも使います。 ただ、以前の検討だとこんなに分かりやすい変化は無かったので、何か違いがあるんでしょうね。 ちなみに、objectでも strength を上げすぎるとアーティファクトが出ます。 何事もやりすぎは良くないです。 下は strength = 0.8まで上げたら出てきました。 続いて、前回の最後にちょろっと書いたクールファイル補正法のpythonスクリプトです。 引き続き検討を進めています。 この記事を書いた後、あぷらなーとさんのXの投稿で「星ナビ８月号に ノイズまみれの画像が付いてくる 」 と知る事ができたので、急いで購入して検証を始めてみました。 ノイズまみれの画像が手に入ると分かった瞬間に 猛烈に喜んで しまいました。 で、画像をダウンロードして圧倒的なノイズを見たら、 更に喜んでしまいました ｗ ノイズに喜んでしまうあたり、いい感じにイカれてきた気がしています。 が、やはり実装は思うようにいきません。 ピクセルマッピングを適用した画像にクール...

 天体撮影をする時に、天気予報が信用しにくいです。 普通の天気予報の「晴れ」は全く当てにならなくて、いつもwindyと睨めっこしています。 これを何とか指標化出来ないかなともがいてみます。 しかも、PCではなくandroidアプリで。 結論を最初に言うとアプリは出来上がりましたが使い物になっていません。 windy は API を提供しているので、ひとまずアカウント登録します。 point forcastのAPI キーを取得しました。 取得できるデータを 確認 し、計算式を考えて、計算アプリを作ってみます（もちろんコーディングは全部Gemini任せ）。 言語はpythonです。 早速、現段階の状況です。 見た目はともかく、私の端末（pixel6a）では狙い通り動いてます。 算出ロジックは下の通り。 総合スコア 総合スコア = (雲スコア * 0.75) + (シーイングスコア * 0.1) + (湿度スコア * 0.05) + (風スコア * 0.1) この計算式は、tomが何となく決めた式です。 何か根拠があるわけではありません。 結果を見ながら修正を加えようと思っています。 雲スコア (配点75%) 下層雲・中層雲・上層雲の合計雲量から算出します。 合計雲量 = min(下層雲[%] + 中層雲[%] + 上層雲 [%], 100)  ※足して100を超える場合、100を上限とする 雲スコア = 100 - 合計雲量 シーイングスコア (配点10%) シーイングはwindy APIで取得できるデータから算出できそうだったので頑張ってみました。 ここでは、大気屈折率構造定数(C n 2 )を、複数の大気層(850hPa〜200hPa)にわたって積分することでシーイングを推定します。 ...

  tomは軍資金が限られた環境で天体撮影をしています。 SE102の１本だけで撮影をしていましたが、他の焦点距離も試したいという気持ちが出てきています。 そこで、ZV-E10 ダブルズームキットに付属のSEL55210を使って撮影しようと思います。 望遠側はF6.3と暗いですが、安いからしょーがない。 まずは猛烈な光害をカットする為にQBPフィルターをレンズに取り付けます。 49-48mm ステップダウンリングで問題なく取り付けられました。 あとはカメラをAZ-GTiに取り付けるだけ。 ただし、L型アリミゾプレートって安いモノでも4400円なんですよね… という事で 安く作成してみました 。 木材 (348円) L型ブラケット (798円) 低頭ねじ (151円) プレート（25円） 木材が若干高くつきましたが… まぁ1300円程度で出来上がれば上等です。 木材はアリガタになる面を鉋で削ってテーパーをつけます。 本来は75°らしいですが、そんな細かい所は気にしません。 テキトーです。目分量で削りました。 もう片方にはプレートを付けて、ネジ締めした時に木材が凹まないようにしました。 木材は購入した長さの1/4くらいしか使わなかったので、今後DIYする際の木製アリ型プレートの予備にします。 これで55-210mmの望遠ズームも土俵に加わりました。 戦力になるかどうか、テストが楽しみです。 自作L字アリミゾブラケットの、もうひとつの使いみちも考えていました。 国立天文台望遠鏡キットを持ってたのですが、SE102を買って出番がなくなっていました。 接続は1/4ユニファイネジなので、今回作ったモノがそのまま流用出来ます。 これにチップスターの筒を投影版として望遠鏡に取り付けて黒点観測をしてみました。 あまりにやっつけ仕事で見た目がアレなので、スマホで撮影した結果だけ。 写真だと微妙ですが、目視だと黒点がそこそこ見えて満足です。 話は変わりますが、L字ブラケットと同時に延長筒も届きました（ 過去投稿参照 ）。 これで望遠鏡設置時のズレも落ち着けば嬉しいです。 ソフト面でもハード面でも着々と準備が整っています。 テスト&撮影が本当に待ち遠しい。