JP2022032358A - 天体写真撮影システム及びプログラム - Google Patents
天体写真撮影システム及びプログラム Download PDFInfo
- Publication number
- JP2022032358A JP2022032358A JP2020136045A JP2020136045A JP2022032358A JP 2022032358 A JP2022032358 A JP 2022032358A JP 2020136045 A JP2020136045 A JP 2020136045A JP 2020136045 A JP2020136045 A JP 2020136045A JP 2022032358 A JP2022032358 A JP 2022032358A
- Authority
- JP
- Japan
- Prior art keywords
- astrophotography
- image
- light
- light pollution
- intermediate image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000012545 processing Methods 0.000 claims abstract description 68
- 238000000034 method Methods 0.000 claims description 31
- 230000008569 process Effects 0.000 claims description 29
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 230000002194 synthesizing effect Effects 0.000 claims 2
- 238000012937 correction Methods 0.000 abstract description 38
- 230000009466 transformation Effects 0.000 abstract description 9
- 238000012935 Averaging Methods 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 27
- 238000001514 detection method Methods 0.000 description 8
- 230000008859 change Effects 0.000 description 6
- 238000013519 translation Methods 0.000 description 6
- 239000002131 composite material Substances 0.000 description 5
- 239000000203 mixture Substances 0.000 description 5
- 230000008832 photodamage Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 238000009499 grossing Methods 0.000 description 3
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 2
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- FFBHFFJDDLITSX-UHFFFAOYSA-N benzyl N-[2-hydroxy-4-(3-oxomorpholin-4-yl)phenyl]carbamate Chemical compound OC1=C(NC(=O)OCC2=CC=CC=C2)C=CC(=C1)N1CCOCC1=O FFBHFFJDDLITSX-UHFFFAOYSA-N 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 235000012308 Tagetes Nutrition 0.000 description 1
- 241000736851 Tagetes Species 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03B—APPARATUS OR ARRANGEMENTS FOR TAKING PHOTOGRAPHS OR FOR PROJECTING OR VIEWING THEM; APPARATUS OR ARRANGEMENTS EMPLOYING ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ACCESSORIES THEREFOR
- G03B15/00—Special procedures for taking photographs; Apparatus therefor
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N23/00—Cameras or camera modules comprising electronic image sensors; Control thereof
- H04N23/60—Control of cameras or camera modules
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Image Processing (AREA)
- Studio Devices (AREA)
Abstract
【課題】 光害に基づく光成分とノイズが除去されて綺麗な星空の天体写真を生成する天体写真撮影システムを提供すること【解決手段】 カメラの連写機能を用いて同一領域の星空を100枚撮影した写真データを記憶する記憶手段と、その写真データから点状の星が存在する星空の天体写真を作成する処理手段を備える。処理手段は、複数の写真データを加算平均して光害に基づく光成分(光害パターン)を求め(S10)、全ての写真データに対し光害パターンを減算して第一中間画像を生成し(S11)、第一中間画像に対し同じ星に対応するピークの位置を合わせるように画像全体をシフトする平面位置補正をして第二中間画像を生成し(S12,S13)、第二中間画像における同じ星に対応するピークの位置を合わせるように射影変換してゆがみを解消した第三中間画像を形成し(S14,S15)、第三中間画像をコンポジット合成して天体写真を生成する機能を備える。【選択図】 図6
Description
本発明は、天体写真撮影システム及びプログラムに関する。
デジタルカメラの高性能化、スマートフォン等の携帯端末に実装されるカメラ機能の高性能化にともない、例えば、夜間等においてもストロボやフラッシュを使用することなく夜景等を綺麗に撮影することが可能となっている。
一方、このように高性能化された一般のデジタルカメラや携帯端末のカメラ機能等を用い、そのまま星空等の天体写真を撮ることは難しい。すなわち、星空は夜景よりもずっと暗いため、例えば、数分~数10分、さらには数時間という長時間の露出時間が必要になる場合がある。そのため、例えば、星空を点像に写すには、シャッターを開けている間、デジタルカメラ等を星の動いている方向に合わせてカメラの向きを少しずつスムーズに変えていくといった追尾撮影を行う必要がある。係る追尾撮影を精度良く行うためには、例えば、モーターが内蔵され自動で星を追いかける機器である赤道儀を用いる必要がある。
上述した赤道儀を用いた天体写真撮影システムは、高価かつ大型化する。そして、そのように大きいことから、当該システムは、持ち運びが煩雑となる。また、例えば予め天体を撮る目的で所望の地域に出かける場合には、荷物が重く大きくなったとしても当該赤道儀を含めた天体写真撮影ステムを持ち運ぶことで天体撮影は可能となる。しかし、係る天体写真撮影システムを常時携帯することは難しい。そのため、例えば外出先でたまたま撮影する環境・状況が良いなど天体写真を撮りたい状況になった場合に、デジタルカメラなどの撮影機器はあるものの赤道儀などのシステム全体を用意できずにシャッターチャンスを逃すおそれもある。特に天体写真の撮影は天候に左右されるため、例えば常時携帯することができるコンパクトで軽量な天体写真撮影システムの提供が望まれる。
また、都市圏で天体写真を撮るためには、光害の影響を考慮する必要がある。すなわち、夜間も経済活動が活発で、街中に人工光が溢れている都市圏の場合、夜空が明るくなり、星が見えにくくなってしまう。そのため、例えば2等星以下の星の明るさは、光害による光の明るさ以下となり写真に写らないことがある。係る光害に対する対策としては、例えば、街灯など主な光害の波長だけをカットする「光害カットフィルタ」をデジタルカメラに装着する方法がある。しかし、ユーザが既に所有しているデジタルカメラに装着可能な光害フィルターが存在しない場合、光害カットフィルタを装着可能なデジタルカメラを改めに購入することになる。また、仮に既に所有しているデジタルカメラに装着可能な光害カットフィルタが存在していたとしても、当該光害カットフィルタを購入する必要がある。そして、天体写真の撮影時には、デジタルカメラに光害カットフィルタを装着し、一般撮影時には光害カットフィルタをデジタルカメラから取り外すとともに保管することになる。このように、光害カットフィルタの着脱や保管をする作業が煩雑となる。さらに、光害カットフィルタは、全ての光害の波長をカットするものではない。よって、例えばカットできない波長の光成分が多い地域では、光害の影響を充分に抑制できないおそれがある。
上述した課題はそれぞれ独立したものとして記載しているものであり、本発明は、必ずしも記載した課題の全てを解決できる必要はなく、少なくとも一つの課題が解決できればよい。またこの課題を解決するための構成についても単独で分割出願・補正等により権利取得する意思を有する。
(1)本発明に係る天体写真撮影システムは、同一領域の星空を撮影した複数の写真データを記憶する記憶手段と、前記複数の写真データから点状の星が存在する星空の天体写真を作成する処理手段を備える。そして前記処理手段は、前記記憶手段に記憶された前記複数の写真データの一部または全部に基づいて光害に基づく光成分の光害パターンを求める機能と、前記記憶手段に記憶された前記複数の写真データに対し、それぞれ前記光害パターンを減算して複数の第一中間画像を生成する機能と、複数の前記第一中間画像間で、同じ星に対応するピークの位置を合わせる処理を行う機能と、位置合わせをして生成した複数の画像をコンポジット合成して前記天体写真を作成する機能を備える。
(2)前記処理手段は、前記複数の写真データの一部または全部をコンポジット合成して前記光害パターンを作成する機能を備える生成するように構成するとよい。
(3)前記位置を合わせる処理は、処理対象の前記第一中間画像を、前記ピークの位置が合うように相対的に画像全体を平面内で移動させることで行うようにするとよい。
(4)前記位置を合わせる処理は、処理対象の前記第一中間画像を、前記ピークの位置が合うように相対的に画像全体を平面内で移動させて第二中間画像を生成し、その生成した複数の前記第二中間画像に対し、前記ピークの位置が合うように射影変換を行い第三中間画像を生成することで行うようにするとよい。
(5)前記天体写真を作成する機能は、前記位置合わせをして生成した複数の画像をコンポジット合成した画像に対し、一定輝度以下の画素を0に変換する処理を行うようにするとよい。
(6)前記天体写真を作成する機能は、前記一定輝度以下でない画素に対し、輝度を強調する処理を行うとよい。
(7)本発明のプログラムは、上記のいずれかに記載の天体写真撮影システムの機能をコンピュータに実現させるためのプログラムである。
本発明は、光害に基づく光成分とノイズが除去或いは低減されるため、綺麗な星空の天体写真を生成することができる。また、係る天体写真は、同一領域を撮影した複数の写真データに基づいて作成されるため、デジタルカメラ等の撮影機器等で済み、星の移動を追いかけるための特別な装置等も不要となる。本発明は、安価・汎用的な機材を用いて、都市近郊等の光害の影響が大きい地域で星空の天体写真をとることができる。
以下、本発明の好適な実施形態について図面に基づき、詳細に説明する。なお、本発明は、これに限定されて解釈されるものではなく、本発明の範囲を逸脱しない限りにおいて、当業者の知識に基づいて、種々の変更、修正、改良を加え得るものである。
図1は、本発明に係る天体写真撮影システムの好適な一実施形態を示している。同図に示すように、本実施形態の天体写真撮影システム10は、デジタルカメラ11を用いて撮影した天体写真(星景写真)の生データに基づいて、所定の画像処理等を実行し、星空を点像に写した天体写真(星景写真)を生成する処理装置13と、その処理装置13に接続される入力装置14、出力装置15並びに記憶装置16等を備える。
処理装置13は、パーソナルコンピュータその他の演算処理能力を有するコンピュータ等であり、実装されたアプリケーションプログラムを実行する機能等を備える。入力装置14は、例えばキーボード、マウス、タッチパネルなどの処理装置13に対して情報,命令を入力する装置である。出力装置15は、例えば、モニタ装置であり、例えば取得したデータ(例えば、天体の生写真データ等)や、取得した生写真データに対して処理装置13が所定の処理をした結果(例えば、星空を点像に写した天体写真等)を表示する装置である。記憶装置16は、処理装置13内のハードディスクその他の内部記憶装置でもよいし、処理装置13に接続される外付けの外部記憶装置でもよい。
デジタルカメラ11で撮影した写真データは、例えばデジタルカメラ11と処理装置13を通信ケーブルで接続したり、無線通信を用いて接続したりした状態で処理装置13にデータ転送することで、処理装置13に取り込む。また、例えば、デジタルカメラ11に着脱自在に装着されたメモリカード等の記録媒体を用いて処理装置13に写真データを入力するようにしてもよい。処理装置13は、このようにして取得した生写真データを、記憶装置16に記憶する。連写した100枚以上の生写真データは、例えば同一のホルダに保存するなど、別に撮影した生写真データと区別がつくように記憶するとよい。
[天体写真の撮影]
本実施形態では、ユーザは、デジタルカメラ11を用いて以下のようにして天体写真(星景写真)を撮影する。撮影地域は、例えば、都市圏等の光害の影響が大きい地域とするとよい。本実施形態では、後述するように光害の影響が大きい地域で撮影した星景写真であっても、星空を点像に写した写真を生成することができる。但し、天体写真撮影システム10は、光害の影響が小さい地域で撮影した生写真データに対する処理を排除するものではない。光害の影響が小さい地域で撮影した生写真データであっても、当該小さい光害の影響やその他の影響をなくし、より綺麗な天体写真を生成することができる。
本実施形態では、ユーザは、デジタルカメラ11を用いて以下のようにして天体写真(星景写真)を撮影する。撮影地域は、例えば、都市圏等の光害の影響が大きい地域とするとよい。本実施形態では、後述するように光害の影響が大きい地域で撮影した星景写真であっても、星空を点像に写した写真を生成することができる。但し、天体写真撮影システム10は、光害の影響が小さい地域で撮影した生写真データに対する処理を排除するものではない。光害の影響が小さい地域で撮影した生写真データであっても、当該小さい光害の影響やその他の影響をなくし、より綺麗な天体写真を生成することができる。
撮影は、光をできるだけ多く取り込むようにする。具体的には、高感度撮影(例えばISO感度は最大)を行い、F値は最小で、シャッター速度は遅めに設定するとよい。シャッター速度が遅めは、例えば、撮影した天体写真の全体が白っぽくなる程度とするとよい。また、手振れ補正はなしとする。そして、固定撮影により、同一の視野内の天体写真を多数枚取得する。具体的には、デジタルカメラ11を三脚12に固定し、連写機能を用いて例えば100枚以上の天体写真を撮影する。また、色温度も固定して撮影を行う。また本実施形態では、連写により100枚程度撮影しており、この撮影に要する時間は2から3分程度となる。
また、撮影は、月のない日に行うのが望ましい。さらにカメラの撮影範囲内に大型の光点(月・街頭)がなく、夜空のみで、地上の風景は画角内に入れないように設定する。さらに、撮影時には、飛行機・人工衛星などの高速に移動する光点が画角内にないようにすると良い。
[光の種類]
図2は、都市圏で撮影した天体写真を示す。同図に示すよう、光害による影響から、点在するはずの星が、ほとんど見える状態で撮影できていない。図2に示す写真に写されている光の種類は、例えば、光害と、星と、ノイズなどがある。写真の各位置の光の強度は、それら3つの光の内の一つまたは複数の光の成分が載ったものとなる。
図2は、都市圏で撮影した天体写真を示す。同図に示すよう、光害による影響から、点在するはずの星が、ほとんど見える状態で撮影できていない。図2に示す写真に写されている光の種類は、例えば、光害と、星と、ノイズなどがある。写真の各位置の光の強度は、それら3つの光の内の一つまたは複数の光の成分が載ったものとなる。
図3は、図2中の矢印で示す各位置における光の強度を示すグラフである。すなわち、写真のX軸上のある点におけるY軸を横軸にとり、各位置の光の強度(明るさ)を縦軸に採った模式図である。また、図4は、図2中の矢印で示す各位置における光害の光成分、星の光成分と、ノイズの光成分のそれぞれの光の強度を示す模式図である。図3に示す模式図は、図4に示す光害と星とノイズの光成分を合わせたものである。
光害は、地上の光が大気散乱したものであり、図4に示すように、連続的に変化する。この光害は、巨視的量依存のため、2~3分程度の短時間では変化しない。すなわち、連写した100枚の写真における光害の光のパターンは、ほとんど変化しない。そして、この光害に基づく光成分は、除外対象である。
ノイズは、センサーに発生するノイズである。このノイズは、熱が原因のため、図4に示すように発生位置はランダムである。ランダム変化のため、写真毎にノイズに基づく光成分の発生位置,光のパターンは異なる。そして、このノイズに基づく光成分は、除外対象である。
星が存在する部分には、その星からの光が存在する。図4に示すように、星が存在する位置にピーク状に光の強度が強くなり、その大きさは星により異なる。1等星以上は十分に明るく肉眼でも視認できる。1等星以上の星の光の強さは、ノイズと比べて十分大きい(図4中、破線枠参照)。一方、2等星以下の星は暗く、ノイズと同程度かそれ以下の強度となる(図4中、一点鎖線枠参照)。よって、2等星以下の星は、熱的ノイズと区別が困難である。また、星は日周運動により常に位置が変わるが相対位置関係はほぼ変わらない。2~3分の時間経過によっても、星の位置が異動するため、撮影した100枚の各写真における同じ星の存在位置が変わる(図5参照)。この星に基づく光成分は、残す対象である。なお、図5は、星の位置を示すピークの移動を模式的に表したものであり、天体における実際の移動は円状に動く。
上述したように、生写真データは、3種類の光の成分を合わせたものであり、図3に示すように、ランダムの振幅で振動しながら全体的にはなだらか変化した波形が現れ、ところどころ(星のある部分)にピークが現れる。図3からも明らかなように、肉眼で確認できる一等星に対応する部分は、明確なピークとして現れるが、その他の星はノイズに隠れた状態となる。
[画像処理]
上述したように、3種類の光の成分が存在すること、及び各光の成分の特徴に基づき、処理装置13は、以下に示す処理アルゴリズムを実行することで、不要な光成分を除去しつつ、2等級以下の星も点像として見える星空の天体写真を生成する。具体的には、処理装置13は、図6に示すフローチャートを実行する。
上述したように、3種類の光の成分が存在すること、及び各光の成分の特徴に基づき、処理装置13は、以下に示す処理アルゴリズムを実行することで、不要な光成分を除去しつつ、2等級以下の星も点像として見える星空の天体写真を生成する。具体的には、処理装置13は、図6に示すフローチャートを実行する。
すなわち、処理装置13は、処理対象の100枚の生写真データにおける光害パターンを求める(S10)。図4に模式図を示したように、光害のパターンはなだらかに変化し、処理対象となる全ての生写真データにおいて共通する。そこで、処理装置13は、100枚の生写真データから複数枚例えば20枚程度をランダムにサンプリングし、抽出した生写真データを位置補正なしでコンポジット合成する。例えば、図7に模式図として示すように、抽出した20枚の生写真データにおける図2中の矢印で示す各位置における光の強度を示すグラフを加算平均すると、上から2番目に示すようなグラフとなる。すなわち、ノイズに伴う光成分は、平均化するとほぼフラットになり、振動したとしてもその振幅は非常に小さくなる。また、星に伴うピークの位置は、時間経過に伴い変化するため、平均化すると影響が少なくなる。画像としてみた場合には、見えにくくなる。一方、光害は、時間経過によって位置変化しない定常パターンのため、結果として、複数の生写真データを平均化すると、主成分として現れる。単純に加算平均処理をしただけだと、ノイズの影響を完全に消すことができず、わずかに振動が残る。さらに例えば星の大きいピークなどの影響が残り、わずかに膨らんだ部分も生じる。そこで、処理装置13は、滑らかになる処理を行う。この滑らかにする処理は、例えば、Medianフィルターを用いたり、関数近似を行ったりするとよい。これにより、図7に示すように、上から2番目に示すグラフは、上から3番目に示すようになだらかな直線からなる光成分が生成される。これにより、光害パターンが抽出される。なお、抽出された光害パターンには、正確にはノイズの平均値も含まれるか、その値はわずかであることから、光害パターンと称する。
係る処理における実際の写真の一例を示すと、図8は、生写真データの一部を拡大して示す図である。同図に示すように、まだらな光害と熱的ノイズに加え、点状の星が存在する。この写真データを含む複数(例えば20枚)の生写真データを加算平均すると、図9に示すようになめらかな光害に線上に星が残った画像となる。そして、係る画像に対し、例えばMedianフィルターで星を消すと、図10に示すようになめらかな光害からなる画像が生成される。
次に処理装置13は、光害減算処理を実行する(S11)。すなわち、処理装置13は、処理対象の100枚の生画像データから、光害に基づく光成分を除去する。具体的には、求めた光害パターンを用いて、処理対象の全ての生写真データから光害パターンを減算し、第一中間画像を生成する。処理対象の生写真データが100枚存在する場合、第一中間画像は100枚生成される。
例えば図2中の矢印で示す各位置における光強度が、図11の上側に示すようになっているとすると、光害減算処理は、その図から、同一箇所における光害パターン(図7の上から3番目の図等)を減算する。すなわち、横軸の同じ座標位置の光強度をそれぞれ減算する。これにより、図11の下側に示すように、ベースラインが横軸と平行なグラフに変換される。
この光害減算処理により生成された第一中間画像の一例は、図12の左側の図となる。この図は、第一中間画像の一部を拡大して示す図であり、図8に示す生写真データが光害減算処理をした結果である。図12の右側の図は、第一中間画像に対して明るさ調整を行った結果である。光害による光成分を除去したため、全体的に暗くなるが、星の部分が明るくなる。但し、一枚の画像を見ただけだと、星の明るさが弱く、左側の第一中間画像では、星空に存在する点状の星の存在が見にくいが、明るさ調整を行った左側の図から明らかなように、光害パターンを除去することで、星が点状に存在した星空が再現されていることが確認できる。
次に処理装置13は、平面位置補正用ピーク検出処理を実行する(S12)。ピークは、星の移動に基づく各画像間の位置補正用のマーカーである。例えば第一中間画像内における1等星に対応する位置をピークとする。具体的には、処理装置13は、処理対象の第一中間画像の各画素のRGBの合計値をそれぞれ求め、第一閾値以上の画素をピーク候補として抽出する。第一閾値は、例えば1等級の星の明るさに対応する値とするとよい。
図13に示すように、上に示す第一中間画像には、同図の下に示すように四角の枠で囲む位置に第一閾値以上のピーク候補が存在する。そして最終的に検出されるピークは、第一閾値以上の全てを抽出するのではなく、20~30程度とするとよい。本実施形態では、20個のピークを検出するようにした。この20個のピークは、例えば、RGBの合計値が大きいものから順に20個をピックアップしてもよいが、検出するピークは、画像全体に散らばっているようにするとよい。そこで、係る20個を検出するためには、例えば、第一閾値以上のピーク候補の中から、20個のピークをランダムに抽出したり、N個のピークはRGBの合計値が大きいものから順に抽出し、残りは抽出されなかった領域からピックアップしたりする等、各種の手法が採れる。処理装置13は、このピーク検出処理を全ての第一中間画像に対して行う。
次に処理装置13は、平面位置補正を行う(S13)。時間経過に伴い夜空に見える星の位置は移動するので、同じ領域を撮影した100枚の写真データ中に存在する同じ星の位置は、全て異なる。そこで、この平面位置補正は、画像間で同じ星が同一位置に来るように位置補正を行う。具体的には、処理装置13は、100枚の第一中間画像の一枚(例えば1枚目に撮影した画像)を読み出して基準画像(Base)にし、別の1枚の画像を読み出して位置補正対象画像(Target)とする(図14中、左側参照)。
次いで、処理装置13は、基準画像Bと位置補正対象画像Tのそれぞれの画像から、ピークを二つ選ぶ。図14の右側の図は、選んだ2つのピークを結ぶ線分を便宜上示している。その線分は、基準画像では破線で示し、位置補正対象画像では一点鎖線で示している。選択した2つのピークが同じ星に基づくものの場合、ピーク間の距離すなわち2つの線分の長さが等しいか、近似し、さらにその存在位置も比較的近い位置にある。そこで、各画像において、ピックアップされた20個中の2つのピークの全ての組み合わせについてピーク間距離を算出し、2つの画像間で距離が類似しているピークの組み合わせを抽出し、以下に示す平面位置補正を行う。なお、ピーク間距離が短すぎる場合には、正確に位置補正できない場合があるため、ピーク間距離が所定の閾値以上のものに対してのみ位置補正を行うようにしてもよい。
平面位置補正は、以下に示す3つの処理を繰り返し実行し、最終的に適正な補正量に基づいて位置補正を行う。図15に示すように、まず、平行移動補正を行い、片方のピークを合わせる(処理1)。すなわち、位置補正対象画像Tを縦方向や横方向に移動し、2つの画像間で片方のピークの位置を重ならせる。このときの平行移動量を記録する。
次いで、回転移動量補正を行い、残りのピークを合わせる(処理2)。すなわち、処理1であわせたピークの位置を回転中心として、位置補正対象画像Tを回転させ、2つの画像間で残りのピークの位置を重ならせる。これにより、2つのピークが重なる。このときの回転移動量を記録する。夜空における星は、北極星を中心に回転するため、単純な平行移動ではなく、上述したように平行移動補正と回転移動補正を行う。
このように平行移動補正と回転移動補正を行った後の補正対象画像Tにおける処理ステップS11で検出した各ピークの位置が、基準画像Bの星の座標上に存在するか確認する(処理3)。そして、係る存在する個数、すなわち、両画像でピークの位置が一致した個数を計数し、記録する。
処理装置13は、全てのピークの組み合わせについて、上記の処理1~3を繰り返し実行し、ピークの位置が一致した個数が最大のものに関連付けられた平行移動量と回転移動量を、平面位置補正移動量に決定する。そして、処理装置13は、位置補正対象画像Tに対し、決定した平面位置補正移動量で補正処理を行い、第二中間画像を生成し、記録する。また、処理装置13は、基準画像以外の99枚の全ての第一中間画像に対して個別に平面位置補正を行い、対応する第二中間画像を生成する。
図11に模式図で示したように、光害減算処理後のグラフは、ベースラインはフラットになるが、ノイズ成分が含まれている。よって、例えば2等星以下の暗い星のピークは、ノイズ成分に隠れてしまい、検出しにくい。そこで、例えば上述した平面位置補正等を行い、ピーク位置をあわせた画像を加算平均すると、ノイズ成分は小さくなるが、星に伴うピークは変わらない。これにより、2等星以下の星であっても、星が存在する座標の明るさが目立つ。
図11に模式図に従って概念を説明すると、光害減算処理後のグラフにおいて一等星の位置を合わせるように画像全体をシフトさせると、図16に示すようになる。そして、係るグラフを加算平均処理すると、図17に示すようになる。すなわち、星のピークは、基準線Lより上にくる。これには、補正後の星の位置は固定されて一致しているため平均化するとピークが明確になる。図示の例では、確認できるピーク個数は、図3に示す補正前では5個であったが、補正後の図17では9個に増え、等級の低い星も検出できる。一方、ノイズによる光成分は基準線Lより下に位置する。これは、ノイズの光強度はランダムのため、平均化すると低い値でフラット化するためである。
ところで、図18に示すように、星は3次元の天球上を移動する。それをカメラで撮影し2次元の写真画像にした場合、3次元→2次元の変換により時間変化により相対位置関係が微妙に変化する。そのため、例えば第二中間画像をそのまま単純にコンポジット合成を行うと、図19に示すように特定領域R1では星は正確にフィッティングされるが、別の領域R2では位置ずれのため流れて合成される。星は連続的に変異することから、正確にフィッティングされない領域R2では、星は線状に延びた広がりをもった状態となる。
そこで本実施形態では、以下に示す射影変換を行い、3次元から2次元に変換する際に生じるゆがみを除去するようにした。具体的には、まず処理装置13は、変換補正用ピーク検出処理を行う(S14)。この変換補正用ピーク検出処理は、まず、前処理で生成した100枚の全ての第二中間画像を用いて比較明合成を行い、得られた画像から変換補正に使用する4つのピークを選定する。
この4つのピークの選定は、処理装置13が、比較明合成した画像に対し、第二閾値を用いたピーク検出処理を行い、検出された多数のピークの中から、下記の条件を満たすピークを検出する。比較明合成をした画像では、日周運動の影響から一つの星が点状ではなく、線状など一定の広がりを持った状態で存在することがある。そして、星は連続的に移動することから、同じ星に基づく光は、比較明合成をした画像においてつながった状態にある。そこで、上記のピークを検出する際に、星は連続的に移動することから、同じ星に基づく光は、比較明合成をした画像においてつながった状態にある。そこで、上記のピークを検出する際に、係る繋がりのある部位は、一つの星に対応するピークとするとよい。
条件1:画像の端にない。日周運動により星が画像外に移動する場合があるため、排除する。
条件2:他のピークから一定以上の距離がある。誤選定防止のためである。
条件3:画像の四隅に最も近いもの。4点で描く四角形の面積がなるべく大きくすると、精度向上するためである。
条件1:画像の端にない。日周運動により星が画像外に移動する場合があるため、排除する。
条件2:他のピークから一定以上の距離がある。誤選定防止のためである。
条件3:画像の四隅に最も近いもの。4点で描く四角形の面積がなるべく大きくすると、精度向上するためである。
図20は、比較明合成をした画像に対し、上記のピークの検出処理をした結果である。図示するように、多数のピーク候補が検出され、最終的に4つの射影変換補正用ピークPが決定される。また、この比較明合成をした画像の右隅を拡大して示す図から明らかなように、星に対応する部分が連続した広がりを持って明るくなっている。
さらに、処理装置13は、射影点の範囲、すなわち、決定した4つのピーク毎に星が動きうる範囲を決定する。具体的には、例えばピークの周りの第三閾値以上の領域を抽出する処理を行うと良い。
次に処理装置13は、射影変換を行う(S15)。すなわち処理装置13は、100枚の第二中間画像の一枚(例えば1枚目に撮影した画像)を読み出して基準画像(Base)にし、別の1枚の画像を読み出して変換補正対象画像(Target)とする。そして、それぞれの画像における4つのピークを選定する(図21参照)。この選定により、図22に模式図で示すように、基準画像(Base)と変換補正対象画像(Target)では、4つのピークを頂点とする四角形が形成される。
次いで処理装置13は、選定した4つの点を結ぶ仮想四角形の各辺の長さ、すなわち、隣接するピーク間距離を求め、基準画像(Base)のものと変換補正対象画像(Target)のものとで、差が一定値以下であることを確認する。このとき、差が一定値よりも大きいと、ピークを誤選定した可能性があるので、処理装置13は、射影変換を行わない。そして、差が一定値以下の場合には、処理装置13は、4つの点を合わせるように画像全体を射影変換する。このようにすることで、1枚の基準画像と、それぞれに射影変換をしてゆがみを直した99枚の射影変換後の画像が生成され、それら全て(100枚)を第三中間画像として記録する。この100枚の第三中間画像における同一の星に対応する光の位置は、3次元から2次元に変換する際のゆがみが解消或いは低減され、一致する。
次に処理装置13は、最終合成処理を行う(S16)。すなわち、処理装置13は、100枚の全ての第三中間画像をコンポジット合成する。ここでは、コンポジット合成は、加算平均を行う。そして、一定輝度以下の画素のRGBの各値を0に変換する。一定輝度以下は、例えば、RGBの各値が設定値未満のものとするとよい。これにより、ノイズを除去する。そして処理装置13は、それ以外の画素の輝度を強調する。この強調処理は、例えば、強調係数を掛けることで行うとよい。このようにすることで、星以外の領域は暗くなり、星の部分の明るさが強調される。特に射影変換をしてゆがみが低減されるため、画像の全体にわたり、星の光を炙り出することができる。これにより、星の部分が点状に現れ、星空の画像が形成される(図23等参照)。処理装置13は、このようにして生成された最終画像を、記憶装置16に記録し、表示装置に出力するとよい。
上述した実施形態では、処理装置13は、スタンドアローンタイプのパーソナルコンピュータを用い、カメラで撮影した画像データを記録し、処理するようにしたが、本発明はこれに限ることはなく、例えば、インターネット等のネットワークに接続されたサーバ等で構成してもよい。また、カメラが例えばスマートフォン等の携帯端末で構成し、その携帯端末に上述した処理アルゴリズムを実行するアプリを実装する。このようにすると、携帯端末を用いて所定枚数の夜景を撮影し、当該アプリで処理することで、綺麗な星空を生成し、表示等することができるので良い。
以上、本発明の様々な側面を実施形態並びに変形例を用いて説明してきたが、これらの実施形態や説明は、本発明の範囲を制限する目的でなされたものではなく、本発明の理解に資するために提供されたものであることを付言しておく。本発明の範囲は、明細書に明示的に説明された構成や製法に限定されるものではなく、本明細書に開示される本発明の様々な側面の組み合わせをも、その範囲に含むものである。本発明のうち、特許を受けようとする構成を、添付の特許請求の範囲に特定したが、現在の処は特許請求の範囲に特定されていない構成であっても、本明細書に開示される構成を、将来的に特許請求する可能性があることを、念のために申し述べる。
10 :天体写真撮影システム
11 :デジタルカメラ
12 :三脚
13 :処理装置
14 :入力装置
15 :出力装置
16 :記憶装置
11 :デジタルカメラ
12 :三脚
13 :処理装置
14 :入力装置
15 :出力装置
16 :記憶装置
Claims (7)
- 同一領域の星空を撮影した複数の写真データを記憶する記憶手段と、
前記複数の写真データから点状の星が存在する星空の天体写真を作成する処理手段を備え、
前記処理手段は、
前記記憶手段に記憶された前記複数の写真データの一部または全部に基づいて光害に基づく光成分の光害パターンを求める機能と、
前記記憶手段に記憶された前記複数の写真データに対し、それぞれ前記光害パターンを減算して複数の第一中間画像を生成する機能と、
複数の前記第一中間画像間で、同じ星に対応するピークの位置を合わせる処理を行う機能と、
位置合わせをして生成した複数の画像をコンポジット合成して前記天体写真を作成する機能を備える天体写真撮影システム。 - 前記処理手段は、前記複数の写真データの一部または全部をコンポジット合成して前記光害パターンを作成する機能を備える生成する請求項1に記載の天体写真撮影システム。
- 前記位置を合わせる処理は、処理対象の前記第一中間画像を、前記ピークの位置が合うように相対的に画像全体を平面内で移動させることで行う請求項1または2に記載の天体写真撮影システム。
- 前記位置を合わせる処理は、処理対象の前記第一中間画像を、前記ピークの位置が合うように相対的に画像全体を平面内で移動させて第二中間画像を生成し、その生成した複数の前記第二中間画像に対し、前記ピークの位置が合うように射影変換を行い第三中間画像を生成することで行う請求項1または2に記載の天体写真撮影システム。
- 前記天体写真を作成する機能は、前記位置合わせをして生成した複数の画像をコンポジット合成した画像に対し、一定輝度以下の画素を0に変換する処理を行うようにした請求項1~4のいずれか1項に記載の天体写真撮影システム。
- 前記天体写真を作成する機能は、前記一定輝度以下でない画素に対し、輝度を強調する処理を行う請求項5に記載の天体写真撮影システム。
- 請求項1から6のいずれかに記載の天体写真撮影システムの機能をコンピュータに実現させるためのプログラム。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020136045A JP2022032358A (ja) | 2020-08-11 | 2020-08-11 | 天体写真撮影システム及びプログラム |
PCT/JP2020/031519 WO2022034697A1 (ja) | 2020-08-11 | 2020-08-20 | 天体写真撮影システム及び記録媒体 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020136045A JP2022032358A (ja) | 2020-08-11 | 2020-08-11 | 天体写真撮影システム及びプログラム |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2022032358A true JP2022032358A (ja) | 2022-02-25 |
Family
ID=80247047
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020136045A Pending JP2022032358A (ja) | 2020-08-11 | 2020-08-11 | 天体写真撮影システム及びプログラム |
Country Status (2)
Country | Link |
---|---|
JP (1) | JP2022032358A (ja) |
WO (1) | WO2022034697A1 (ja) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013186879A (ja) * | 2012-03-12 | 2013-09-19 | Casio Comput Co Ltd | 撮像装置、その画像合成方法、及びプログラム |
WO2014038629A1 (ja) * | 2012-09-07 | 2014-03-13 | 株式会社Ihi | 移動体検出方法 |
JP2017208016A (ja) * | 2016-05-20 | 2017-11-24 | 三菱電機株式会社 | 移動体検出装置、観測システム及び移動体検出方法 |
-
2020
- 2020-08-11 JP JP2020136045A patent/JP2022032358A/ja active Pending
- 2020-08-20 WO PCT/JP2020/031519 patent/WO2022034697A1/ja active Application Filing
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013186879A (ja) * | 2012-03-12 | 2013-09-19 | Casio Comput Co Ltd | 撮像装置、その画像合成方法、及びプログラム |
WO2014038629A1 (ja) * | 2012-09-07 | 2014-03-13 | 株式会社Ihi | 移動体検出方法 |
JP2017208016A (ja) * | 2016-05-20 | 2017-11-24 | 三菱電機株式会社 | 移動体検出装置、観測システム及び移動体検出方法 |
Also Published As
Publication number | Publication date |
---|---|
WO2022034697A1 (ja) | 2022-02-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5871862B2 (ja) | 3d深度情報に基づいたイメージぼかし | |
KR101643607B1 (ko) | 영상 데이터 생성 방법 및 장치 | |
CA2662355C (en) | Mosaic oblique images and methods of making and using same | |
US10447999B2 (en) | Alignment of images of a three-dimensional object | |
US20110222757A1 (en) | Systems and methods for 2D image and spatial data capture for 3D stereo imaging | |
CN113191954B (zh) | 一种基于双目相机的全景图像拼接方法 | |
CN109299662B (zh) | 深度数据计算设备与方法及人脸识别设备 | |
WO2015164025A1 (en) | Generation and use of a 3d radon image | |
CN109544620A (zh) | 图像处理方法和装置、计算机可读存储介质和电子设备 | |
EP3494693A1 (en) | Combining images aligned to reference frame | |
WO2008114264A2 (en) | A method and apparatus for video image stabilization | |
Přibyl et al. | Feature point detection under extreme lighting conditions | |
KR100834157B1 (ko) | 영상 합성을 위한 조명환경 재구성 방법 및 프로그램이기록된 기록매체 | |
CA2870480A1 (en) | Hybrid precision tracking | |
Okumura et al. | Augmented reality based on estimation of defocusing and motion blurring from captured images | |
JPWO2016175044A1 (ja) | 画像処理装置及び画像処理方法 | |
JP2014051212A (ja) | 移動体検出方法 | |
CN111899345A (zh) | 一种基于2d视觉图像的三维重建方法 | |
JP2012075088A (ja) | 画像処理システムおよび画像処理方法 | |
Song et al. | Advanced underwater image restoration in complex illumination conditions | |
JP4843787B2 (ja) | 被写体追尾方法及び装置 | |
JP6292968B2 (ja) | 疑似hdr画像推定装置及び方法 | |
WO2022034697A1 (ja) | 天体写真撮影システム及び記録媒体 | |
JP2016110312A (ja) | 画像処理方法、画像処理装置及びプログラム | |
Li et al. | Improving high dynamic range image based light measurement |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20200828 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20210929 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20220330 |