異なる鏡筒･カメラの画像をコンポジット_その3（関数電卓からエクセルへ）

わが越後の梅雨明けは8月に入ってからになりそうなので、引き続き同じテーマでの記事になります。今回「異なる鏡筒・カメラの画像をコンポジット」の題材に選んだのは、カシオペア座の端にある「バブル（しゃぼんだま）星雲 付近」の過去の撮影画像です。
以下の３つが今回選んだ題材で、以前の処理画像になります。撮影時期が新しい順に並べると
２０２３年８月１２日　撮影鏡筒：R200SS（PH）fl=760mm　カメラ：ASI ５３３ＭＣ　↑
２０１７年9月２６日　撮影鏡筒：R200SS（PH）fl=760mm　カメラ：Cooled ６０Ｄ　↑
２０１５年9月１０日　撮影鏡筒：R200SS（PH）fl=760mm　カメラ：Cooled ６０Ｄ　↑

作業手順は　前回記事　に記載したように①．各画像に共通する2つの星を選んで･････････････････････････････････････････････②．それぞれのピクセル座標をメモする････････････････････････････････････････････➂．「ピタゴラスの定理」を使って関数電卓で２星間のピクセル数を算出･･････････････④．２星間のピクセル数が同じになるよう、ステライメージで「画像解像度」を変更する⑥．画像解像度をそろえた画像をコンポジットする･･････････････････････････････････⑦．コンポジットできた画像に対して、通常の画像処理を行う････････････････････････なのですが関数電卓を使って２星間のピクセル数を求める作業がとても疲れました。まして対象画像が３枚、４枚にもなれば、コンポジットが終わった時点でぐったり。改めて画像処理をやりなおすエネルギーが残っていません。そこで、”電卓”ではなくエクセルを使うことにしました。最近は一覧表の作成などで使うことが多かったのですが、エクセル本来の表計算として久しぶりに使ってみました。
以下が作成したエクセルのシートで、今回のコンポジットまでの作業にさっそく使ってみました。入力が必要なのは、使用する各元画像のﾋﾟｸｾﾙサイズと、各画像の２つの星のピクセル座標だけ（あとは必要によりメモとして記入）以上を入力するだけで①．各画像の２星間のピクセル数↓②．各画像の画像解像度の比率↓➂．各画像のステライメージでの「画像解像度」設定値（どの画像を基準にするかで選択）
今回は２星間のピクセル数の多い画像Ａ（画像解像度の高い画像）を基準画像としたので、画像Ａはそのまま変更せずに、画像Ｂ・画像Ｃの画像解像度を表示値に変更しました。
（参考）：[２星間のﾋﾟｸｾﾙ数] で 使用しているエクセル関数は　=SQRT((H7-L7)^2+(I7-M7)^2)　･･･ 画像Aの場合
今回の３つの画像のうち、２０１７年・２０１５年の画像は同じ鏡筒。同じカメラで撮影しているのになぜかこれまでコンポジット処理された形跡がありません。そこで今回 画像B・画像Cとして別画像として処理したのですが、上のエクセル計算表を見ると、２星間のピクセル数が微妙に異なっています。比率でみればわずか0.2％の差ですが、ﾋﾟｸｾﾙ数では７ピクセルになり重ね合わせると星のズレがわかります。撮影時の気温の差などで鏡筒の焦点距離が微妙に変化したのでしょうか？
３年分の撮影画像を画像解像度変更後コンポジット、画像処理しなおした結果です。↓（　各年の画像を加算平均　完成画像を約６０％に縮小したもの　）
今回の画像統合の効果を検証するため、バブル星雲付近の等倍切り出し（縮小なし）画像で比較します。
評価のポイントとしては、2023年のASI533MCで撮った画像と比較してどうか？　というところでしょうか。

= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =今回で３回目となるこの連載記事、ふとそのむかし同じような事を書いた気がして過去の記事を調べてみたら、10年前の2014年の1月に同じような連載記事を書いていました。当時のカメラは初期のデジ一眼canon kiss DXで、2代目のCooled６０Dに代わってまもないころ。10年も経ってカメラはフルサイズのEOS６Ｄ、
更に昨年からはCMOSカメラのASI５３３MCに代わりましたが、10年前は雪雲で撮れない冬休み、今回は同じく雲の多い夏休みで、ブログネタがなくて書いた記事というのは共通しています。
先日は孫と一緒に花火を見に地元の夏祭りへ、今回は夜店のくじは はずれ だったみたい

にほんブログ村
↑ いつもポチッとありがとうございます。
= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =

異なる鏡筒･カメラの画像をコンポジット その2（作業手順）

前回記事 に続き、今回は過去に撮った「鏡筒・カメラの異なる画像のコンポジット」の作業手順について。使用するのは ステライメージ９ と 関数電卓。
コンポジット（重ね合わせ）というのは撮影時に何枚も撮っておいて、それを画像処理時に重ね合わせる事によりノイズを低減するという基本的な方法ですが、 ”同じ鏡筒（含む補正レンズ）と同じカメラ” で撮影したものであれば画像毎の多少のずれがあってもステライメージなどを使って重ね合わせることができます。しかしながら同じ撮影対象であっても、鏡筒の焦点距離やカメラが変わると重ね合わせはできません。（わたしが知らないだけで、世の中にはこれができるソフトもあるのかも･･）
今回、この通常はできない ”異なる鏡筒と異なるカメラ” で撮った撮影画像のコンポジット（重ね合わせ）に挑戦した理由は、前回記事にも書いたように●この十余年で同じ対象を撮った画像は増えたが、撮影鏡筒（補正レンズを含む）およびカメラが異なっているため過去画像の有効活用ができない･･･････････････････●年間を通じて夜間晴天日数の少ない日本海側では、総露光時間で10時間も撮るには
何年もかかってしまう　（その間ずっと同じ光学系・カメラというのはつらい）･････････････
そして、もっとも大きな理由が■現在 「星のふるさと館 」で展示中のわたしの写真展で、９月に予定している差し替え画像の準備で、過去画像を加えた画像の再処理が可能なら活用したい、というもの。
前置きが長くなりましたが、それでは 今回行った作業手順について説明します。□　過去の撮影画像の選択9月からの「秋・冬 編」の対象でもある、NGC7331 & ｽﾃﾌｧﾝの五つ子 から 2017年・2020年に撮影したものをピックアップ。詳細な撮影データは画像拡大で見てもらうとして、「異なる鏡筒と異なるカメラ」に関するデータについては◇2017/ 9･･VC200L(旧ﾚﾃﾞﾕｰｻ) fl=1,278㎜　Cooled６０D 4.3μ/pixel　総露光時間 １３０分◇2020/ 8･･VC200L(ﾚﾃﾞﾕｰｻHD) fl=1,386㎜　EOS６D(HKIR) 6.5μ/pixel　総露光時間 ９２分
鏡筒の焦点距離が長いと画像解像度は大きくなり、カメラの１ﾋﾟｸｾﾙのサイズが大きくなると画像解像度は小さくなる。＊画像解像度･･･ある写野角をどれだけのピクセル数で描いているか ということで、画像解像度が合っていない画像ではコンポジットでズレが発生する。
上の二つの画像を比較すると、焦点距離では2020年の方が 画像解像度 が大きく他方、カメラの1ﾋﾟｸｾﾙのサイズでは2017年の方が 画像解像度 は大きくなります。結局のところ、この二つの画像をコンポジットするにはそれぞれトータルでの画像解像度を正確に求める事が必要となります。
STEP -Ⅰ ．各画像の画像解像度を求める
以下の作業はステライメージ９と 関数電卓 を使用しました（＊ 左側が2017年、右側が2020年の画像で、それぞれダーク・フラット処理およびコンポジットまで終えて保存されていた画像を利用しました）
手順 ①　ステライメージの基準点指定を使って、両画像に共通な２つの星をマーキングします
手順 ②　画像を拡大してそれぞれ（４つの星）の（X、Y）のピクセル座標を読み取りメモします
手順 ➂　次に画像解像度 算出のもとになる2つの星の間隔（ﾋﾟｸｾﾙ数）を、直角三角形の斜辺算出でおなじみの「ピタゴラスの定理」を使って求めます。（お忘れの方はこちら）→　2星間のピクセル数（斜辺）＝√ ΔX（底辺）² ＋ ΔY（高さ）²上の2017年の例では　√ (2942-395) ²　＋(4422-568) ²　＝4619.58　となります
手順 ④　算出した２星間のピクセル数から、2017年の画像の方がﾄータルで画像解像度が高い事がわかりました

STEP -Ⅱ ．二つの画像の画像解像度を合わせる
ステライメージの 画像 ﾒﾆｭｰの「画像解像度」を使って両画像の画像解像度を合わせるのですが、そろえる方法が二通りあります。（方法Ａ）　画像解像度の低いほうを上げる（低い方の画像が拡大される）（方法Ｂ）　画像解像度の高いほうを下げる（高い方の画像が縮小される）わたしが選んだのは（方法Ａ）でした。（理由. 縮小するのは後でも可能だから）
（　ステライメージ９の画像解像度変更の画面　）
手順 ①　手順Ⅰで求められた２星間のピクセル数の値から、解像度の低い2020年の画像解像度を高くします。では現在のピクセル数をどれだけ増やせば画像解像度が同じになるのか比率を求めます４６１９．５８ ÷ ３２２７．１６＝ １．４３１４６・・・現在の「幅」または「高さ」のピクセル値を ×１．４３１４６･･倍 すれば2017年の画像の解像度と同じになるはずです。＊「縦横比」が”固定”になっていれば、縦・または横のピクセル数を変えるだけで もう一方も自動的に変わります
手順 ②　画像解像度変更後の両画像を比較します。
表示拡大率をそろえて２星間の間隔を比較すると大きな計算間違いはないようです

STEP -Ⅲ ．二つの画像をコンポジット（重ね合わせ）する
関数電卓まで使ってここまでやってきたのは、鏡筒もカメラも異なる画像をコンポジットするためでした。はたしてピッタリ重なってくれるのか？
実はコンポジットを行う その前に、もう一つ選択しなければならないことが。それはコンポジットする際に「加算平均」か「加算」か？いつもなら「加算平均」なのですが、鏡筒が異なれば画像の明るさ（Ｆ値）も大きく違います。おまけに設定感度も様々で、極端に撮影画像の明るさに差がある事も考えられます。それを撮影データから推測することはかなりやっかいなので、画像の明るさが大きく異なっても無難そうな「加算」で行くことにしました。
実際にコンポジットで加算をおこなった画像です。 ↓これまで５、６タイトルほど実施したのですが、不思議とズレずにうまく重なってくれました。ズレる事もありましたが、その原因は星のピクセル位置の読み間違いでした。（老眼なもので）
ここまでやれれば、あとはコンポジットができた範囲をトリミングしてステライメージなどで画像処理を行うだけ。画像処理を終えた統合画像になります。↓総露光時間が２２０分に増えたことにより、従来よりきつめの画像処理が可能となりました
肝心なのは、これだけ手間をかけて元の画像より改善されたかどうかですが ･･画像の一部を切り取って比較してみました。ＮＧＣ ７３３１ 付近左から　２０１７年・２０２０年の完成画像　そして右側が今回 両画像を統合して処理した画像です
ステファンの五つ子 画像の順番はＮＧＣ７３３１と同じ
実は画像解像度を上げる際の画像拡大でダメージが出るのではと思ったのですが、意外に健闘しています。ただ、鏡筒・カメラの組み合わせによっては画像解像度に差がありすぎて破綻することも考えられます。
そこで 鏡筒・カメラの組み合わせによる画像解像度を推定できる一覧表を作ってみました。これまでに撮影した画像が多そうなVC200L（HD）とEOS ６D(HKIR改造)の組み合わせを基準として作成してあります
数値が近い組み合わせ同士なら、画像解像度の変換も控えめに済むのではないかと考えています。
引き続き差し替え用画像に使える画像が得られるか検証を進めていきます。

= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =余り期待していなかったのですが、もしかしたら･･という結果が。すべての画像とはいかないでしょうが、過去の画像もコンポジットに使えるとなれば少なくともプリントが必要な写真展示には活用できそうです。
現在使用している真四角センサーのカメラ ASI５３３ＭＣはプリントが苦手です。

にほんブログ村
↑ いつもポチッとありがとうございます。
= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =

過去の異なる鏡筒･カメラの画像をコンポジット（PixInsightじゃないけど･･）

関東甲信は梅雨明けしたようですが、わが「越後」は置いてけぼり。もっとも梅雨が明けても夜に雲の多い天気は毎年経験済み。となれば9月から予定している「星のふるさと館」の展示写真の入れ替えの準備でもと考えて始めました。
現在展示中の「春・夏」編の展示レイアウト （展示してかなり経ったので今回はボカシなしで）現在展示されている写真は「春・夏」編なので、次は「秋・冬」編ということでこれまで10年余りで撮りだめた中から見栄えのする対象を選んで、撮影画像の再処理後にプリント。となるのですが、”再処理”といっても同じステライメージを使用するわけで、マジックのような画像の改善は期待できません。そこで新たな試みとして同じ対象を撮った 鏡筒・カメラ の異なる画像をコンポジットして総露光時間を増やしてからより強い画像処理を実施してみてどう変わるか？
退職後に「天体写真」を始めてじきにわかったことは、何枚も撮って「総露光時間」を増やす事の効果は絶大であるという事。天候に恵まれない日本海側でも、何年もかけて撮り続ければ総露光時間10時間越えも可能である。と思っていたのですが、しかし、現実は･･何年もかけて同じ対象を撮影したものの、その間に撮影鏡筒（含む補正レンズ）および カメラが変わりそのままではコンポジットができない画像が増えていきました。
主な 撮影鏡筒（含む補正レンズ）および 撮影カメラ の変遷･･････････････････････････････････････････････････････････････････■2008. 4 EOS kissDX (ﾌｲﾙﾀｰ改造) APS-C ﾋﾟｸｾﾙｻｲｽﾞ=5.7μm･･････●2009.5_ Vixen R200SS反射鏡筒・ｺﾏｺﾚｸﾀｰ3(fl= 800㎜）･･･････■2012. 7 SEO Cooled60D APS-C ﾋﾟｸｾﾙｻｲｽﾞ=4.3μm･････････････ 〇 2013. 2_ MPCC-MARKⅢｺﾚｸﾀｰ (R200SS) (fl=800㎜）●2014. 9_ Vixen VC200L反射鏡筒(fl=1,800㎜）中古････････････ 〇 2014. 9_ ﾚﾃﾞﾕｰｻｰ (VC200L) (fl=1,278㎜）･･････････ 〇 2015. 7_ ｺﾚｸﾀｰPH (R200SS) (fl=760㎜）･･･････････ 〇 2018. 2_ ﾚﾃﾞﾕｰｻｰHD (VC200L) (fl=1.386㎜）･･･････■2018. 6 EOS 6D (HKIR改造) ﾌﾙｻｲｽﾞﾋﾟｸｾﾙｻｲｽﾞ=6.5μm　･････････■2023. 4 ASI 533MC Pro CMOSｶﾒﾗ ﾋﾟｸｾﾙｻｲｽﾞ=3.76μm････････････････････････････････････････････････････････････････････････････約１５年の歳月とは言え、ずいぶんポチッたものですが一番お金がかかったのはここにはないタカハシの赤道儀だと思います
今回はいかにポチッてきたかという話ではないので、話を戻すと上の履歴から撮影鏡筒が増え、補正レンズが変わるたびに鏡筒の焦点距離が変わり撮影画像の拡大率も変化。（焦点距離が長くなると狭い範囲を拡大して撮影）それに加えてカメラが変わると１ﾋﾟｸｾﾙの画素ｻｲｽﾞも違って、同じ範囲の描画画素数も変化。せめてカメラがずっと同じなら、鏡筒が変わっても写る範囲が変わるだけで基本的にコンポジット（重ね合わせ）は容易だと思われるのですが、最も技術の進化の恩恵を受けるデジタルカメラでそれは無理な話です。
結局、撮影対象は同じでも撮った時点での鏡筒（補正レンズを含む）と撮影カメラの組み合わせによって、そのままでは重ねることのできない画像が増えてしまいました。これを重ねるためには一方の画像を他方の画像に合わせて画像の解像度を変更する必要があります。その後に重ね合わせる事により、総露光時間が増えた効果が得られるのではないか？
どうせ当分撮影できないのだし、成果が無くても差し替え画像の選択には役立つだろうと考えました。「秋・冬」編の中心タイトルになるであろう「アンドロメダ銀河」「オリオン大星雲」は後回しにして他の小粒でも個性的な星雲・銀河の数タイトルを終えてみた感想では、組み合わせる画像同士が良好なら、総露光時間が増えるメリットが得られるという結果が。（あたりまえ？）
次回はその作業手順と、すでに検証が済んだタイトルの結果報告の予定です。
= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =実は当初、写真展の差し替え画像の作成にPixInsightのフリーライセンスを使おうかと思いました。定期購読している「星ナビ」の読者の入選者はほとんど使っています。フリーライセンスは45日間無料で使用できることからどんなものか「習うより、慣れろ」と考えたのですが、今月の「上越天体写真友の会(J-APA)」の月例会でPixInsightを1年使用している人に聞いたところ未だ機能の全貌把握を行っている段階で、45日間では成果は無理、とのこと。お盆は東京の孫たちも来るし、と今回は断念しました。いつかは使う日がくるのでしょうか？

にほんブログ村
↑ いつもポチッとありがとうございます。
= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =

14年後の皆既日食は見れるか？

（お詫び）記事タイトルには”14年後”とありますが、記事作成時点では”11年後”が正解でした
前回撮影したのが5月17日の上弦過ぎの月ですから、早2ヶ月近く撮影できていません。ずいぶん長い期間 趣味の「天体写真」から遠ざかっているのですが、梅雨時につき毎年の事ではあります。これまでは せめてブログだけでも･･と企画を立案して過去に撮った画像を使ったシリーズものを掲載したりして凌（しの）いでいたのですが･･今年はモチベーション（動機づけ・やる気 ）が一向に生まれてきません。（ 人はそれを”老化”という ? ）
そんなある日小二になった孫娘が学校から「星のふるさと館」の夏休み行事のポスターをもらってきてそこに 「じぃじ の名前がのってるよ」 というのです。上越市内の小学校が夏休みに入る前に「星のふるさと館」が配布したようです。「星のふるさと館」では今年4月から11月末まで私の写真展が開催中になっています。写真展「雲の上には宇宙(そら)」　裏庭で撮った大宇宙　春・夏　編
さすがに一年間同じ写真じゃ見飽きちゃうだろうと、９月には作品の総入れ替えをやる予定とは言ったものの、展示する写真についてはまだ何も検討していません。過去の画像再処理やレイアウト検討、そしてプリントアウトまで考えるとそろそろ始めないと。
そこで退職により単身赴任生活の終わった２０１２年ごろから現在に至る撮影画像を再チェック。次は「秋・冬」シーズンに撮ったものが中心になるのですが、太陽や月や惑星の写真も。そこで目にとまったのが、『星ナビ 7月号』の「星空人生１００年計画」という特集記事。その中でこれまでに起こった「Ｓ級天文現象リスト(1950年以降)」の一覧がありその中で、退職した２０１２年に起こった天文現象の２つを自宅で撮影していました。
その一つ目が２０１２年５月２１日の日本の太平洋側で見られた金環日食。東京をはじめ太平洋側の大都市では金環日食となったものの、天候は曇りで雲の合間から見えたようです一方で新潟県上越市では部分日食となったものの、まあまあの天気で日食経過の大部分を撮影することができました撮影ＤＡＴＡ　：　２０１２／ ５／２１　タカハシ セミアポ９０㎜　（ｆl＝１２００㎜）ISO ２００　　NDﾌｨﾙﾀｰ　EOS kiss Ｘ２　タカハシ ＥＭ－２００ ＵＳＤ　　

二つ目 は　上記の日食の翌月２０１２年６月６日に起きた金星の日面通過です。撮影鏡筒、カメラは　上の日食と同じ　撮影時刻の記録が残っていなくて、もっとも深く入り込んだ頃の撮影と思われます
当時の撮影ノートのメモ書きを見ると、”見た目も地味で特に準備もせずに撮影した” とありました。実は見た目の地味さとは裏腹に金星の太陽面通過は非常にまれな現象で、この撮影以後は 地球上では２１１７年12月まで見ることができないそうです。その時わたしは１７０歳を数えており、とうてい現世で見ることはできません。
一方で先に紹介した日食の方はというと、あと１１年待てば本当の皆既日食が上越市の自宅で見れる可能性があります。その時わたしは男性の平均余命を越えているのですが、なぜか見れそうな気がしています。
2012年の金環食ほか、この先起こる皆既・金環日食帯の予報です。↓富山天文台のＨＰより
ちなみに今回掲載した２つの写真の撮影画像は残っていないため、再処理せずに差し替え後の写真展に使うかもしれません。

= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =7月の新月期もまったく撮影できずに終わろうとしています。また新しい企画でも考えないとブログの維持ができそうにありません。我が家の『かたてま野菜』はこのところの雨模様と しつこい蚊で草取りもできず。今では草に埋もれてしまったスイカは15回以上も人工授粉したのに、育ったのは2個だけ。このスイカは6月27日に７番目に受粉したスイカ、もう一つはこの二日前に受粉した６番目になります
カラスにやられなければ 今月下旬から来月上旬にかけてが収穫時期。

にほんブログ村
↑ いつもポチッとありがとうございます。
= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =