JP2005331320A - 天空率および日照時間算出システム、および、算出プログラム - Google Patents
天空率および日照時間算出システム、および、算出プログラム Download PDFInfo
- Publication number
- JP2005331320A JP2005331320A JP2004148911A JP2004148911A JP2005331320A JP 2005331320 A JP2005331320 A JP 2005331320A JP 2004148911 A JP2004148911 A JP 2004148911A JP 2004148911 A JP2004148911 A JP 2004148911A JP 2005331320 A JP2005331320 A JP 2005331320A
- Authority
- JP
- Japan
- Prior art keywords
- sky
- image
- calculated
- calculating
- photographic 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
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
【課題】 熟練した技術無しに、簡易に、正確な天空率などを算出する。
【解決手段】 システム10は、少なくとも何れかの画像に、太陽が写し込まれ、少なくとも何れかの画像に垂直線が含まれ、かつ、それぞれの画像において、隣接する画像との間で2つの重複点が設定された複数の写真画像の画像データを記憶する記憶装置20と、写真画像中の太陽の位置、垂直線上の任意の2点の位置、および/または、重複点の位置に基づいて、写真画像の左右の向き、写真画像中、写真画像の重複点で囲まれた領域に相当する区分画像の任意の点の仰角および方位角を算出する高度/方位角算出部12と、区分画像中の画素値に基づいて、天空部分を抽出し、算出された左右方向の向き、仰角および方位角に基づいて、区分画像中の任意の微小領域の立体角投射率を算出し、天空部分に含まれる微小領域の立体角投射率を積算することにより天空率を算出する天空率算出部16とを備える。
【選択図】 図1
Description
本発明は、太陽を含んだ天空の撮影画像に基づいて、天空率および日照時間を算出するシステムおよびプログラムに関する。
建築基準法の改正により、建物の建築基準に天空率が導入されている。したがって、建築物を建築する際などにおいて、天空率を算出することが重要になっている。天空率とは、被照面の一点から天空を見るときに建物などの障害物を除いて実際に見える割合をいい、天空要素の水平面に対する立体角投射率に相当する。また、日照時間とは、日の出から日没までの間で、快晴時に実際に日照のある時間をいう。以下に述べるように、本実施の形態においては、一点(特に、屋内の日照時間を問題にする場合には、窓際)に直射日光が照射する時数が求められる。
「だれにもわかる日影規制とその対策」(武井正昭、オーム社、1978年) 第54頁〜第61頁、第111頁〜第129頁 「日本建築学会設計計画パンフレット 24 日照の測定と検討」(日本建築学会編、彰国社 昭和52年)
「だれにもわかる日影規制とその対策」(武井正昭、オーム社、1978年) 第54頁〜第61頁、第111頁〜第129頁 「日本建築学会設計計画パンフレット 24 日照の測定と検討」(日本建築学会編、彰国社 昭和52年)
たとえば、特許文献1や特許文献2には、魚眼レンズを装着したカメラを、光軸を垂直方向に向けて撮影し、撮影された写真画像から日照時間を算出する手法が紹介されている。しかしながら、方向を正確に測定し、或いは、光軸が垂直方向になるようにカメラの向きを設定するなど、撮影には熟練した技術が必要である。
また、得られた写真画像に基づいて、魚眼レンズの射影方式にしたがって、太陽軌道を記入する際にも熟練した技術が必要である。また、方向の測定誤差や光軸の誤差により、記入された太陽軌道から算出される日照時間も誤差が多いという問題点があった。
本発明は、熟練した技術無しに、簡易に、正確な天空率および日照時間を算出することができるシステムおよびプログラムを提供することを目的とする。
本発明の目的は、複数の写真画像であって、少なくとも何れかの画像に、太陽が写し込まれているとともに、少なくとも何れかの画像に垂直線が含まれ、かつ、それぞれの画像において、隣接する画像との間で2つの重複点が設定された写真画像の画像データを記憶する記憶装置と、前記写真画像中の太陽の位置、垂直線上の任意の2点の位置、および/または、重複点の位置に基づいて、写真画像の左右の向き、写真画像中、写真画像の重複点で囲まれた領域に相当する区分画像の任意の点の仰角および方位角を算出する位置情報算出手段と、前記区分画像中の画素値に基づいて、天空部分を抽出する天空部分抽出手段と、前記位置情報算出手段により算出された左右方向の向き、仰角および方位角に基づいて、区分画像中の任意の微小領域の立体角投射率を算出し、前記天空部分に含まれる微小領域の立体角投射率を積算することにより天空率を算出し、算出された天空率を記憶装置に記憶する天空率算出手段と、を備えたことを特徴とする天空率算出システムにより達成される。
好ましい実施態様においては、カメラの水平姿勢が保たれた状態で複数の写真画像の全てが撮影されている場合に、前記位置情報算出手段が、写真画像中の任意の点の仰角および方位角を算出するように構成されている。
別の好ましい実施態様においては、前記位置情報算出手段が、処理対象となる区分画像および隣接する区分画像の情報にしたがって、算出すべき未知数、および、当該未知数を得るための条件式を選択するように構成されている。
また、本発明の目的は、魚眼レンズにて撮影された写真画像であって、少なくとも太陽が写し込まれるとともに、垂直線が含まれる写真画像の画像データを記憶する記憶装置と、前記写真画像中の太陽の位置、垂直線上の任意の2点の位置に基づいて、前記魚眼レンズの左右方向の回転角、写真画像中の任意の点の迎角および方位角を算出する位置情報算出手段と、
前記写真画像の重複点で囲まれた領域に相当する区分画像において、前記区分画像中の画素値に基づいて、天空部分を抽出する天空部分抽出手段と、前記位置情報算出手段により算出された左右方向の向き、仰角および方位角に基づいて、区分画像中の任意の微小領域の立体角投射率を算出し、前記天空部分に含まれる微小領域の立体角投射率を積算することにより天空率を算出し、算出された天空率を記憶装置に記憶する天空率算出手段と、を備えたことを特徴とする天空率算出システムにより達成される。
前記写真画像の重複点で囲まれた領域に相当する区分画像において、前記区分画像中の画素値に基づいて、天空部分を抽出する天空部分抽出手段と、前記位置情報算出手段により算出された左右方向の向き、仰角および方位角に基づいて、区分画像中の任意の微小領域の立体角投射率を算出し、前記天空部分に含まれる微小領域の立体角投射率を積算することにより天空率を算出し、算出された天空率を記憶装置に記憶する天空率算出手段と、を備えたことを特徴とする天空率算出システムにより達成される。
本発明の別の実施態様においては、上述した天空率算出システムと、前記太陽の位置、撮影時刻、並びに、撮影場所の経度および緯度に基づいて、太陽軌道を算出する太陽軌道算出手段と、前記区分画像中の天空部分に相当する領域と、太陽軌道とを対応付けて、前記天空部分を通る太陽軌道の部分を特定して、日照時間を算出する日照時間を算出する日照時間算手段とを備えたことを特徴とする日照時間算出システムが提供される。
また、本発明の目的は、複数の写真画像であって、少なくとも何れかの画像に、太陽が写し込まれているとともに、少なくとも何れかの画像に垂直線が含まれ、かつ、それぞれの画像において、隣接する画像との間で2つの重複点が設定された写真画像の画像データを記憶する記憶装置を備えたコンピュータより読み出し可能なコンピュータプログラムであって、前記コンピュータを、前記写真画像中の太陽の位置、垂直線上の任意の2点の位置、および/または、重複点の位置に基づいて、写真画像の左右の向き、写真画像中、写真画像の重複点で囲まれた領域に相当する区分画像の任意の点の仰角および方位角を算出する位置情報算出手段、前記区分画像中の画素値に基づいて、天空部分を抽出する天空部分抽出手段、並びに、前記位置情報算出手段により算出された左右方向の向き、仰角および方位角に基づいて、区分画像中の任意の微小領域の立体角投射率を算出し、前記天空部分に含まれる微小領域の立体角投射率を積算することにより天空率を算出し、算出された天空率を記憶装置に記憶する天空率算出手段として機能させることを特徴とする天空率算出プログラムにより達成される。
好ましい実施態様においては、カメラの水平姿勢が保たれた状態で複数の写真画像の全てが撮影されている場合に、前記コンピュータを位置情報算出手段として機能させる際に、写真画像中の任意の点の仰角および方位角を算出するように、前記コンピュータを動作させる。
別の好ましい実施態様においては、前記コンピュータを位置情報算出手段として機能させる際に、処理対象となる区分画像および隣接する区分画像の情報にしたがって、算出すべき未知数、および、当該未知数を得るための条件式を選択するように、前記コンピュータを動作させる。
また、本発明の目的は、魚眼レンズにて撮影された写真画像であって、少なくとも太陽が写し込まれるとともに、垂直線が含まれる写真画像の画像データを記憶する記憶装置を備えたコンピュータにより読み出し可能なコンピュータプログラムであって、前記コンピュータを、前記写真画像中の太陽の位置、垂直線上の任意の2点の位置に基づいて、前記魚眼レンズの左右方向の回転角、写真画像中の任意の点の迎角および方位角を算出する位置情報算出手段、前記写真画像の重複点で囲まれた領域に相当する区分画像において、前記区分画像中の画素値に基づいて、天空部分を抽出する天空部分抽出手段、並びに、前記位置情報算出手段により算出された左右方向の向き、仰角および方位角に基づいて、区分画像中の任意の微小領域の立体角投射率を算出し、前記天空部分に含まれる微小領域の立体角投射率を積算することにより天空率を算出し、算出された天空率を記憶装置に記憶する天空率算出手段として機能させることを特徴とする天空率算出プログラムによっても達成される。
別の好ましい実施態様においては、前記コンピュータを、前述した天空率算出プログラムを構成する手段として機能させ、かつ、前記コンピュータを、前記太陽の位置、撮影時刻、並びに、撮影場所の経度および緯度に基づいて、太陽軌道を算出する太陽軌道算出手段、並びに、前記区分画像中の天空部分に相当する領域と、太陽軌道とを対応付けて、前記天空部分を通る太陽軌道の部分を特定して、日照時間を算出する日照時間を算出する日照時間算手段として機能させることを特徴とする日照時間算出プログラムが提供される。
本発明によれば、熟練した技術無しに、簡易に、正確な天空率および日照時間を算出することができるシステムおよびプログラムを提供することが可能となる。
以下、添付図面を参照して、本発明の実施の形態について説明する。図1は、本発明の第1の実施の形態にかかる天空率/日照時間算出システムの構成を示すブロックダイヤグラムである。図1に示すように、天空率/日照時間算出システム(以下、単に「システム」とも称する。)10は、天空率等の算出対象となった地点からの撮影画像の画像データ等を記憶する記憶装置12と、画像データを取り込んで、高度および方位角を算出する高度/方位角算出部14と、画像中の天空部分を分離する天空部分分離部16と、分離された天空部分の画像に基づいて、天空率を算出する天空率算出部18と、天空部分と太陽軌道とに基づいて日照時間を算出する日照時間算出部20とを備えている。また、天空率/日照時間算出システム10は、キーボードやマウスなどの入力装置22および表示装置24を備えている。
本実施の形態にかかる天空率/日照時間算出システムは、パーソナルコンピュータに、所定のプログラムをインストールすることにより実現される。高度/方位角算出部14、天空部分分離部16、天空率算出部18および日照時間算出部20の機能は、記憶装置12に記憶されたプログラムにより実現される。また、記憶装置12には、算出された天空率や日照時間が格納される。
まず、本実施の形態において用意される撮影画像について説明する。撮影者は、ディジタルカメラなどを領して以下の手順で画像を撮影する。
(1)まず、任意の位置で1枚目の画像を撮影する。
(2)次いで、1枚目の画像に隣接するように2枚目の画像を撮影する。ここでは、通常の連続写真と同様に、撮影者は同じ視点から撮影し、かつ、隣接する写真が、少なくとも一部で重なり合うようにしておく。
(3)2枚目の写真と同様に、3枚目以降も、隣接する写真と少なくとも一部が重なり合うようにして撮影を続ける。
(4)少なくとも1枚の写真には、太陽が写り込むようにする。なお、太陽が写り込む画像は、他の写真と同様にNDフィルタ未装着で撮影したもの、および、NDフィルタ装着で撮影したものの2種類を用意する。
(5)以下のいずれかのケースが生じた場合に、撮影は終了する。
i)写真に写り込む対象が天空のみとなった場合。
ii)撮影位置から見えるすべての天空の画像の撮影が終了した場合。
(1)まず、任意の位置で1枚目の画像を撮影する。
(2)次いで、1枚目の画像に隣接するように2枚目の画像を撮影する。ここでは、通常の連続写真と同様に、撮影者は同じ視点から撮影し、かつ、隣接する写真が、少なくとも一部で重なり合うようにしておく。
(3)2枚目の写真と同様に、3枚目以降も、隣接する写真と少なくとも一部が重なり合うようにして撮影を続ける。
(4)少なくとも1枚の写真には、太陽が写り込むようにする。なお、太陽が写り込む画像は、他の写真と同様にNDフィルタ未装着で撮影したもの、および、NDフィルタ装着で撮影したものの2種類を用意する。
(5)以下のいずれかのケースが生じた場合に、撮影は終了する。
i)写真に写り込む対象が天空のみとなった場合。
ii)撮影位置から見えるすべての天空の画像の撮影が終了した場合。
なお、本実施の形態においては、このようにして撮影者により撮影された個々の撮影画像から、後述する重複点で囲まれた領域を抽出することにより得られた画像を「区分画像」と称する。
太陽が写り込んだ区分画像の光軸点の方位角を求めるために、以下のいずれかの手法が採用される。
(a)上記(1)〜(5)において説明した撮影プロセスにおいて、少なくとも何れかの区分画像が、カメラの光軸或いはカメラの左右の傾きを水平に保持して撮影される。
(b)後述するシステムの処理において、ある写真画像において、垂直線(子午線)を手動で設定する。ここで、垂直線を設定する写真画像には、撮影の際に、建物の垂直線や下げ振りによる垂直線を写し込んでおくのが望ましい。また、複数の写真画像に垂直線を入れておくことにより、誤差を減じることが可能となる。なお、上記(a)の条件を満たして撮影された写真画像、或いは、(b)にあるように垂直線が含まれる写真画像を、本実施の形態において、「基準画像」と称する。
(a)上記(1)〜(5)において説明した撮影プロセスにおいて、少なくとも何れかの区分画像が、カメラの光軸或いはカメラの左右の傾きを水平に保持して撮影される。
(b)後述するシステムの処理において、ある写真画像において、垂直線(子午線)を手動で設定する。ここで、垂直線を設定する写真画像には、撮影の際に、建物の垂直線や下げ振りによる垂直線を写し込んでおくのが望ましい。また、複数の写真画像に垂直線を入れておくことにより、誤差を減じることが可能となる。なお、上記(a)の条件を満たして撮影された写真画像、或いは、(b)にあるように垂直線が含まれる写真画像を、本実施の形態において、「基準画像」と称する。
以下、本実施の形態にかかるシステムを利用した天空率および日照時間の算出原理を説明する。
(1)天球を、同一の視点から撮影した複数の平面の区分画像で表わす。
(2)隣接する区分画像は、双方の区分画像に写し込まれた共通の2点(これを「重複点」と称する)の高度および方位角の情報を用いて繋げられる。
(1)天球を、同一の視点から撮影した複数の平面の区分画像で表わす。
(2)隣接する区分画像は、双方の区分画像に写し込まれた共通の2点(これを「重複点」と称する)の高度および方位角の情報を用いて繋げられる。
この区分画像の接続において、最も効率がよいのは、ある写真画像のそれぞれの角付近の4点を設定し、隣接する右側の写真画像との間では、右列の2点を重複点とし、隣接する上側の写真画像との間では、上行の2点を重複点とする(さらに、左側および下側に隣接する写真画像との間でも同様に重複点を設定する)場合である。
つまり、写真画像中の4点を結んだ4つの線分(対角線を除く)で切り取った四角形を、区分画像として、高度方向および方位角方向に繋いでいくのが最も望ましい。なお、隣接する区分画像は、重複点における高度の情報および方位角の情報が一致するだけであり、実際の3次元空間上で繋がるものではないことに留意すべきである。
(3)それぞれの写真画像の任意の2点の、絶対空間における高度(水平0)、方位角(南0)が判れば、その写真画像内のすべての点の高度および方位角を決定することができる。たとえば、写真画像の中央点(光軸点に相当する)の高度、方位角を代表として利用することができる。
(4)一連の区分画像から、天空部分のみが抽出され、天空率が算出される。
(5)天空部分と太陽軌道の重なりから、日照時間が算出される。
(3)それぞれの写真画像の任意の2点の、絶対空間における高度(水平0)、方位角(南0)が判れば、その写真画像内のすべての点の高度および方位角を決定することができる。たとえば、写真画像の中央点(光軸点に相当する)の高度、方位角を代表として利用することができる。
(4)一連の区分画像から、天空部分のみが抽出され、天空率が算出される。
(5)天空部分と太陽軌道の重なりから、日照時間が算出される。
以下、より詳細に、本実施の形態にかかる天空率/日照時間算出システムにおいて実行される処理について説明する。
前述したように、操作者は、隣接する写真画像を選択し、表示装置の画面上に表示させて、双方の画像において2点の重複点を選び、これを登録する(図2のステップ201)。システム10は、2つの写真画像のそれぞれの重複点の座標(写真画像中の相対座標)を、写真画像と対応付けて、記憶装置12に記憶する(ステップ202)。この処理が、全ての隣接する画像について繰り返される。重複点の選択が終了すると(ステップ203でイエス(Yes))、操作者は、ある写真画像を基準画像と決定し、写真画像中、垂直線に相当する部分を指定する(ステップ204)。システム10は、決定された写真画像について、垂直線の座標を、当該写真画像と関連付けて、記憶装置12に記憶する(ステップ205)。なお、隣接する写真画像の間で全ての重複点が設定されることにより、写真画像中の4つの重複点で囲まれる領域である区分画像を繋げていくことが可能となる。
システム10は、操作者からの処理開始の指示を受信すると、図3に示すような処理を実行する。
図3に示すように、本実施の形態においては、システム10の高度/方位角算出部14は、記憶装置12から、1または複数の写真画像のデータ(画像データ)を取り込み(ステップ301)、画像データ中の未知数を特定し(ステップ302)、当該未知数を解くための条件式のパラメータの値を取得した上で(ステップ303)、条件式を解いて未知数を算出する(ステップ304)。算出された未知数は、記憶装置12に記憶される(ステップ305)。このような処理を、全ての写真画像の画像データについて実行する(ステップ306参照)。
区分画像中に太陽或いは垂直線が含まれているか否か、若しくは、撮影時にカメラの水平位置が確保されていたか否かなどに応じて、ステップ302の未知数の数が変動する。以下、具体的に、幾つかの場合について詳細に説明する。
[区分画像に太陽および垂直線が含まれる場合]
(1)この場合、未知数は、以下の3つのパラメータとなる。
[区分画像に太陽および垂直線が含まれる場合]
(1)この場合、未知数は、以下の3つのパラメータとなる。
θe:区分画像(カメラ)の左右方向の傾き
h0 :区分画像(カメラ)の仰角
w0:区分画像(カメラ)の方位角(真南0°)
(2)また、3つのパラメータを解くための3つの条件式は、以下の(式1)〜(式3)となる。
h0 :区分画像(カメラ)の仰角
w0:区分画像(カメラ)の方位角(真南0°)
(2)また、3つのパラメータを解くための3つの条件式は、以下の(式1)〜(式3)となる。
図4(a)および図4(b)は、それぞれ、写真画像における極座標、および、仰角と焦点距離を示す図である。式3は、写真画像に含まれる子午線(垂直線)の上端および下端の方位角が等しいことを示している。
また、図5は、写真画像に含まれる子午線(垂直線)および太陽の関係を示す図、図6は、各パラメータの関係を示す図である。図6に示すように、未知数がわかることにより、写真画像中の任意の1点の位置を表わすことが可能となる。
[隣接する区分画像(区分画像2)に、太陽或いは垂直線が含まれる場合]
(1)この場合には、未知数は、以下の6つのパラメータとなる。
[隣接する区分画像(区分画像2)に、太陽或いは垂直線が含まれる場合]
(1)この場合には、未知数は、以下の6つのパラメータとなる。
θe1 :区分画像1(カメラ)の左右方向の傾き
h01
:区分画像1(カメラ)の仰角
w01
:区分画像1(カメラ)の方位角(真南0°)
θe2 :区分画像2(カメラ)の左右方向の傾き
h02
:区分画像2(カメラ)の仰角
w02
:区分画像2(カメラ)の方位角(真南0°)
(2)上記6つのパラメータを解くための条件式は、以下の(式4)〜(式10)となる。
h01
:区分画像1(カメラ)の仰角
w01
:区分画像1(カメラ)の方位角(真南0°)
θe2 :区分画像2(カメラ)の左右方向の傾き
h02
:区分画像2(カメラ)の仰角
w02
:区分画像2(カメラ)の方位角(真南0°)
(2)上記6つのパラメータを解くための条件式は、以下の(式4)〜(式10)となる。
[カメラの水平姿勢が確保されていた場合]
(1)この場合には、区分画像の傾きは0°であるため、未知数は、以下の2つのパラメータとなる。
h0 :区分画像(カメラ)の仰角
w0:区分画像(カメラ)の方位角(真南0°)
(2)上記2つのパラメータを解くための条件式は、以下の(式11)および(式12)となる。
w0:区分画像(カメラ)の方位角(真南0°)
(2)上記2つのパラメータを解くための条件式は、以下の(式11)および(式12)となる。
なお、他の隣接する区分画像においても、重複点の条件(つまり、双方の方位角および高度が同一であること)を加えれば、同様に、全ての未知数(各区分画像における左右の傾き、方位角、仰角)を算出することが可能となる。
区分画像の数、未知数、条件式の数は、以下の関係を有する。
区分画像数=2:未知数:3×2=6
条件式の数:2×1(重複点の数)×2(高度および方位角)
+3(太陽の位置および垂直線)=7
区分画像数=3:未知数:3×3=9
条件式の数:2×2(重複点の数)×2(高度および方位角)
+3(太陽の位置および垂直線)=11
区分画像数=4:未知数:3×4=12
条件式の数:2×3(重複点の数)×2(高度および方位角)
+3(太陽の位置および垂直線)=15
区分画像数=5:未知数:3×5=15
条件式の数:2×4(重複点の数)×2(高度および方位角)
+3(太陽の位置および垂直線)=19
このように、未知数の数<条件式の数となるため、未知数を求めることが可能である。なお、区分画像が、基準画像と離間するにしたがって、上述するように未知数は増加するので、実際には、太陽が写し込まれた画像と、基準画像(垂直線が設定された画像)は隣接しているのが望ましい。
条件式の数:2×1(重複点の数)×2(高度および方位角)
+3(太陽の位置および垂直線)=7
区分画像数=3:未知数:3×3=9
条件式の数:2×2(重複点の数)×2(高度および方位角)
+3(太陽の位置および垂直線)=11
区分画像数=4:未知数:3×4=12
条件式の数:2×3(重複点の数)×2(高度および方位角)
+3(太陽の位置および垂直線)=15
区分画像数=5:未知数:3×5=15
条件式の数:2×4(重複点の数)×2(高度および方位角)
+3(太陽の位置および垂直線)=19
このように、未知数の数<条件式の数となるため、未知数を求めることが可能である。なお、区分画像が、基準画像と離間するにしたがって、上述するように未知数は増加するので、実際には、太陽が写し込まれた画像と、基準画像(垂直線が設定された画像)は隣接しているのが望ましい。
このようにして、区分画像の傾き、区分画像中の任意の位置の仰角、真南を0°とした場合の方位角を求めることが可能となる。
次に、システム10の天空部分分離部16の処理について説明する。図9(a)に示すように、天空部分分離部16は、写真画像の画像データを記憶装置12から取り出し(ステップ901)、画像データ中、区分画像を構成する画素のそれぞれの輝度および色度に基づいて、一定の輝度および色度を有するものを、天空部分として抽出する(ステップ902)。天空部分の画素の情報(たとえば、区分画像中の画素の位置)は、記憶装置12に記憶される(ステップ903)。このような処理を、全ての写真画像の画像データについて実行する(ステップ904参照)。
このような処理の後、天空率算出部18により天空率が算出される。より詳細には、図9(b)に示すように、区分画像の画像データのうち、天空部分を構成する画素に関する情報を取り出す(ステップ911)。次いで、天空率算出部18は、画素の面積dS、画素の高度h、および、当該画素が含まれる写真画像の画像中心の光軸と、画素と視点とを結ぶベクトルとのなす角αを取得して(ステップ912)、以下の(式13)を用いて、立体角投射率を算出する(ステップ913)。ここにいう画素の高度hは、画素の実際の高度である。
立体角投射率=dS*cosα*sinh (式13)
上記式13において、角αは、先に算出された画素の高度、方位角および左右方向の傾きを利用して得ることができる。また、「cosα*sinh」は、画素の面積に対する重みに相当する。
上記式13において、角αは、先に算出された画素の高度、方位角および左右方向の傾きを利用して得ることができる。また、「cosα*sinh」は、画素の面積に対する重みに相当する。
天空率算出部18は、取得された各画素の天空率照射率を記憶する(ステップ914)。このような処理を、全ての写真画像の画像データについて実行し(ステップ915参照)、最終的に、得られた各区分画像中の画素の天空率照射率を積算することにより、天空率を算出する(ステップ916)。なお、算出された天空率は、記憶装置12に記憶される(ステップ917)
天空率が算出されると、日照時間算出部20が、以下のような手順で太陽の方位角を求める。図10に示すように、日照時間算出部20は、撮影場所の緯度および経度を得る(ステップ1001)。緯度および経度を、操作者が入力しても良い。或いは、他の手段により、システム10の記憶装置12に、緯度および経度を予め記憶させておいても良い。次いで、撮影時刻の標準時Tが、真太陽時Tsに変換される(ステップ1002)。これは、以下の(式14)を用いて得ることができる。
天空率が算出されると、日照時間算出部20が、以下のような手順で太陽の方位角を求める。図10に示すように、日照時間算出部20は、撮影場所の緯度および経度を得る(ステップ1001)。緯度および経度を、操作者が入力しても良い。或いは、他の手段により、システム10の記憶装置12に、緯度および経度を予め記憶させておいても良い。次いで、撮影時刻の標準時Tが、真太陽時Tsに変換される(ステップ1002)。これは、以下の(式14)を用いて得ることができる。
Ts=T+e+(L−135)×4分 (式14)
ここに、eは均時差、Lは、撮影場所の軽度である。
ここに、eは均時差、Lは、撮影場所の軽度である。
次いで、日照時間算出部20は、時角tを算出する(ステップ1003)。時角tは、以下の(式15)に基づいて得ることができる。
Ts=12.0+t/15 (式15)
その後、太陽高度sinhsおよび太陽方位角sinAsが、それぞれ、(式16)および(式17)にしたがって算出される(ステップ1004)。
その後、太陽高度sinhsおよび太陽方位角sinAsが、それぞれ、(式16)および(式17)にしたがって算出される(ステップ1004)。
sinhs=sinφ・sinδ+cosφ・cosδ・cost (式16)
ここに、hsは太陽高度、θは緯度、δは太陽赤緯
sinAs=cosθ・sint/coshs (式17)
次いで、日照時間算出部20は、得られた太陽高度および方位角を、区分画像上の各画素の高度および方位角と対応付け、これにより区分画像上の太陽軌道を求める(ステップ1005)。次いで、区分画像上に現れる太陽軌道から、日照時間が算出される(ステップ1006)。算出された日照時間は、記憶装置12に記憶される。
ここに、hsは太陽高度、θは緯度、δは太陽赤緯
sinAs=cosθ・sint/coshs (式17)
次いで、日照時間算出部20は、得られた太陽高度および方位角を、区分画像上の各画素の高度および方位角と対応付け、これにより区分画像上の太陽軌道を求める(ステップ1005)。次いで、区分画像上に現れる太陽軌道から、日照時間が算出される(ステップ1006)。算出された日照時間は、記憶装置12に記憶される。
このように、本実施の形態によれば、通常のディジタルカメラなどにて撮影した写真画像を用いて、正確な天空率を算出し、かつ、天空部分を通る太陽軌道から日照時間を算出することが可能となる。
次に、本発明の第2の実施の形態について説明する。第1の実施の形態においては、通常のディジタルカメラなどを利用して撮影した画像に基づいて、区分画像をつなぎ合わせ、各区分画像中の位置を算出している。これに対して、第2の実施の形態では、カメラに魚眼レンズを装着し、レンズを上方に向けて撮影した画像を利用する。以下、魚眼レンズにより撮影した写真(魚眼写真)の座標系の変換について説明する。
絶対座標系において、X軸は東西軸、Y軸は南北軸、Z軸は垂直軸であると考える。魚眼レンズ中の座標系(魚眼座標系)のX軸、Y軸、Z軸のそれぞれと、絶対座標軸のX軸、Y軸、Z軸のそれぞれとの間のずれを見出すことにより、魚眼座標系の位置を、絶対座標系の位置に変換することが可能となる。
たとえば、図11に示すように、絶対系座標の魚眼写真は、当該絶対系座標のZ軸(Z’軸)まわりに方位角w0だけ回転することで、絶対座標系のY軸(Y’軸)と魚眼座標系のY軸(Y1軸)とを一致させ(この状態を「系1」と称する。)、X軸(X1軸)まわりに仰角h0だけ回転させることにより、絶対座標系のZ軸(Z’軸)と魚眼座標系のZ軸(Z2軸)とを一致させ(この状態を「系2」と称する。)、かつ、Y軸(Y2軸)まわりにθだけ回転させることにより、最終的に、絶対座標系の全ての軸と、魚眼座標系の全ての軸とを一致させることができる(この状態を、「系3」とも称する。)。
実際には、魚眼レンズを装着したカメラにより撮影された画像中の任意の点は、魚眼座標系で表わされる。したがって、
(1)X軸が水平方向になるように、Y軸まわりにθだけ回転
(2)Z軸が垂直方向になるように、X軸まわりに仰角h0だけ回転
(3)Y軸が南北方向となるように、Z軸まわりに方位角Woだけ回転
を実現することにより、任意の点を絶対座標系で表わすことが可能となる。なお、撮影時にカメラのX軸が水平になるように維持されている場合には、(2)のプロセスは不要となる。
(1)X軸が水平方向になるように、Y軸まわりにθだけ回転
(2)Z軸が垂直方向になるように、X軸まわりに仰角h0だけ回転
(3)Y軸が南北方向となるように、Z軸まわりに方位角Woだけ回転
を実現することにより、任意の点を絶対座標系で表わすことが可能となる。なお、撮影時にカメラのX軸が水平になるように維持されている場合には、(2)のプロセスは不要となる。
なお、本実施の形態においては、魚眼レンズの射影方式は等距離射影であり、図12の極座標によると、(式18)により、(式19)に示す「u」を求めることができる。
絶対座標系における任意の点pは、以下の(式20)〜(式22)で表わすことができる。
太陽が撮影された画像において、太陽方位角をWs、太陽高度をhsとする。なお、これらは、時刻、緯度および経度より算出され、絶対系座標で表わすことが出来る。また、太陽の位置の魚眼系座標を(Hsf、wsf)とする。
当該写真画像においては、第1の実施の形態と同様に、写真画像において垂直線が写し込まれるようにするのが望ましい。また、操作者は、入力装置を操作して、記憶装置12に記憶された写真画像のデータを読み出して、表示装置の画面上に写真画像を表示させ、垂直線を設定しておく。設定された垂直線の情報(垂直線の両端の魚眼系座標(h1f,w1f),(h2f,w2f))は、記憶装置12に記憶される。
第2の実施の形態にかかる高度/方位角算出部14も、第1の実施の形態と同様に、写真画像のデータ(画像データ)を取り込み、画像データ中の未知数を特定し、当該未知数を解くための条件式のパラメータの値を取得した上で、条件式を解いて未知数を算出する(ステップ304)。算出された未知数は、記憶装置12に記憶される。これらの処理は、図3のステップ301〜ステップ305に相当する。なお、第2の実施の形態においては、単一の写真画像のデータに対してのみ処理を施せばよい。
第2の実施の形態においては、未知数および条件式は、以下のように表わされる。
(1)未知数は、絶対座標系で、以下の6つとなる。
(1)未知数は、絶対座標系で、以下の6つとなる。
w0:カメラの方位角(真南0°)
h0:カメラの仰角
θ:カメラの左右方向の回転角
h1:垂直線上の第1の指示点(たとえば、一方の端点)
h2:垂直線上の第2の指示点(たとえば、他方の端点)
w12:垂直線上の双方の指示点の方位角
(2)上記6つのパラメータを解くための6つの条件式は、以下の(式35)〜(式40)となる。
h0:カメラの仰角
θ:カメラの左右方向の回転角
h1:垂直線上の第1の指示点(たとえば、一方の端点)
h2:垂直線上の第2の指示点(たとえば、他方の端点)
w12:垂直線上の双方の指示点の方位角
(2)上記6つのパラメータを解くための6つの条件式は、以下の(式35)〜(式40)となる。
第2の実施の形態において、魚眼レンズを利用することで、単一の写真画像(実際には、フィルタを使用して太陽位置が確認できるような写真画像と、通常の写真画像)のみを利用して、正確な天空率および日照時間を算出することが可能となる。
本発明は、以上の実施の形態に限定されることなく、特許請求の範囲に記載された発明の範囲内で、種々の変更が可能であり、それらも本発明の範囲内に包含されるものであることは言うまでもない。
10 天空率・日照時間算出システム
12 高度/方位角算出部
14 天空部分分離部
16 天空率算出部
18 日照時間算出部
20 記憶装置
12 高度/方位角算出部
14 天空部分分離部
16 天空率算出部
18 日照時間算出部
20 記憶装置
Claims (10)
- 複数の写真画像であって、少なくとも何れかの画像に、太陽が写し込まれているとともに、少なくとも何れかの画像に垂直線が含まれ、かつ、それぞれの画像において、隣接する画像との間で2つの重複点が設定された写真画像の画像データを記憶する記憶装置と、
前記写真画像中の太陽の位置、垂直線上の任意の2点の位置、および/または、重複点の位置に基づいて、写真画像の左右の向き、写真画像中、写真画像の重複点で囲まれた領域に相当する区分画像の任意の点の仰角および方位角を算出する位置情報算出手段と、
前記区分画像中の画素値に基づいて、天空部分を抽出する天空部分抽出手段と、
前記位置情報算出手段により算出された左右方向の向き、仰角および方位角に基づいて、区分画像中の任意の微小領域の立体角投射率を算出し、前記天空部分に含まれる微小領域の立体角投射率を積算することにより天空率を算出し、算出された天空率を記憶装置に記憶する天空率算出手段と、を備えたことを特徴とする天空率算出システム。 - カメラの水平姿勢が保たれた状態で複数の写真画像の全てが撮影されている場合に、前記位置情報算出手段が、写真画像中の任意の点の仰角および方位角を算出するように構成されたことを特徴とする請求項1に記載の天空率算出システム。
- 位置情報算出手段が、処理対象となる区分画像および隣接する区分画像の情報にしたがって、算出すべき未知数、および、当該未知数を得るための条件式を選択するように構成されたことを特徴とする請求項1に記載の天空率算出システム。
- 魚眼レンズにて撮影された写真画像であって、少なくとも太陽が写し込まれるとともに、垂直線が含まれる写真画像の画像データを記憶する記憶装置と、
前記写真画像中の太陽の位置、垂直線上の任意の2点の位置に基づいて、前記魚眼レンズの左右方向の回転角、写真画像中の任意の点の迎角および方位角を算出する位置情報算出手段と、
前記写真画像の重複点で囲まれた領域に相当する区分画像において、前記区分画像中の画素値に基づいて、天空部分を抽出する天空部分抽出手段と、
前記位置情報算出手段により算出された左右方向の向き、仰角および方位角に基づいて、区分画像中の任意の微小領域の立体角投射率を算出し、前記天空部分に含まれる微小領域の立体角投射率を積算することにより天空率を算出し、算出された天空率を記憶装置に記憶する天空率算出手段と、を備えたことを特徴とする天空率算出システム。 - 請求項1ないし請求項6の何れか一項に記載の天空率算出システムと、
前記太陽の位置、撮影時刻、並びに、撮影場所の経度および緯度に基づいて、太陽軌道を算出する太陽軌道算出手段と、
前記区分画像中の天空部分に相当する領域と、太陽軌道とを対応付けて、前記天空部分を通る太陽軌道の部分を特定して、日照時間を算出する日照時間を算出する日照時間算手段とを備えたことを特徴とする日照時間算出システム。 - 複数の写真画像であって、少なくとも何れかの画像に、太陽が写し込まれているとともに、少なくとも何れかの画像に垂直線が含まれ、かつ、それぞれの画像において、隣接する画像との間で2つの重複点が設定された写真画像の画像データを記憶する記憶装置を備えたコンピュータより読み出し可能なコンピュータプログラムであって、前記コンピュータを、
前記写真画像中の太陽の位置、垂直線上の任意の2点の位置、および/または、重複点の位置に基づいて、写真画像の左右の向き、写真画像中、写真画像の重複点で囲まれた領域に相当する区分画像の任意の点の仰角および方位角を算出する位置情報算出手段、
前記区分画像中の画素値に基づいて、天空部分を抽出する天空部分抽出手段、並びに、
前記位置情報算出手段により算出された左右方向の向き、仰角および方位角に基づいて、区分画像中の任意の微小領域の立体角投射率を算出し、前記天空部分に含まれる微小領域の立体角投射率を積算することにより天空率を算出し、算出された天空率を記憶装置に記憶する天空率算出手段として機能させることを特徴とする天空率算出プログラム。 - カメラの水平姿勢が保たれた状態で複数の写真画像の全てが撮影されている場合に、前記コンピュータを位置情報算出手段として機能させる際に、写真画像中の任意の点の仰角および方位角を算出するように、前記コンピュータを動作させることを特徴とする請求項6に記載の天空率算出プログラム。
- 前記コンピュータを位置情報算出手段として機能させる際に、処理対象となる区分画像および隣接する区分画像の情報にしたがって、算出すべき未知数、および、当該未知数を得るための条件式を選択するように、前記コンピュータを動作させることを特徴とする請求項6に記載の天空率算出プログラム。
- 魚眼レンズにて撮影された写真画像であって、少なくとも太陽が写し込まれるとともに、垂直線が含まれる写真画像の画像データを記憶する記憶装置を備えたコンピュータにより読み出し可能なコンピュータプログラムであって、前記コンピュータを、
前記写真画像中の太陽の位置、垂直線上の任意の2点の位置に基づいて、前記魚眼レンズの左右方向の回転角、写真画像中の任意の点の迎角および方位角を算出する位置情報算出手段、
前記写真画像の重複点で囲まれた領域に相当する区分画像において、前記区分画像中の画素値に基づいて、天空部分を抽出する天空部分抽出手段、並びに、
前記位置情報算出手段により算出された左右方向の向き、仰角および方位角に基づいて、区分画像中の任意の微小領域の立体角投射率を算出し、前記天空部分に含まれる微小領域の立体角投射率を積算することにより天空率を算出し、算出された天空率を記憶装置に記憶する天空率算出手段として機能させることを特徴とする天空率算出プログラム。 - 前記コンピュータを、請求項6ないし9の何れか一項に記載の手段として機能させ、かつ、前記コンピュータを、
前記太陽の位置、撮影時刻、並びに、撮影場所の経度および緯度に基づいて、太陽軌道を算出する太陽軌道算出手段、並びに、
前記区分画像中の天空部分に相当する領域と、太陽軌道とを対応付けて、前記天空部分を通る太陽軌道の部分を特定して、日照時間を算出する日照時間を算出する日照時間算手段として機能させることを特徴とする日照時間算出プログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2004148911A JP2005331320A (ja) | 2004-05-19 | 2004-05-19 | 天空率および日照時間算出システム、および、算出プログラム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2004148911A JP2005331320A (ja) | 2004-05-19 | 2004-05-19 | 天空率および日照時間算出システム、および、算出プログラム |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2005331320A true JP2005331320A (ja) | 2005-12-02 |
Family
ID=35486084
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2004148911A Pending JP2005331320A (ja) | 2004-05-19 | 2004-05-19 | 天空率および日照時間算出システム、および、算出プログラム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2005331320A (ja) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012111215A1 (ja) * | 2011-02-18 | 2012-08-23 | Jx日鉱日石エネルギー株式会社 | 日射評価装置、日射評価方法及び日射評価プログラム |
CN102721988A (zh) * | 2012-06-14 | 2012-10-10 | 中国气象科学研究院 | 基于天空可见光图像的日照时数测量方法 |
JP2015132539A (ja) * | 2014-01-14 | 2015-07-23 | 株式会社東芝 | 日照量計算装置、経路提案装置および日照量計算方法 |
CN104897130A (zh) * | 2015-06-18 | 2015-09-09 | 广西壮族自治区气象减灾研究所 | 天基遥感分块插值计算太阳高度角的方法 |
JP2019085056A (ja) * | 2017-11-10 | 2019-06-06 | 東日本旅客鉄道株式会社 | 軌道周辺の環境変化検出方法 |
CN112785678A (zh) * | 2019-11-09 | 2021-05-11 | 李本彦 | 基于三维模拟的日照分析方法和系统 |
CN113939710A (zh) * | 2019-05-29 | 2022-01-14 | 古野电气株式会社 | 信息处理系统、方法及程序 |
-
2004
- 2004-05-19 JP JP2004148911A patent/JP2005331320A/ja active Pending
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012111215A1 (ja) * | 2011-02-18 | 2012-08-23 | Jx日鉱日石エネルギー株式会社 | 日射評価装置、日射評価方法及び日射評価プログラム |
JP2012173057A (ja) * | 2011-02-18 | 2012-09-10 | Jx Nippon Oil & Energy Corp | 日射評価装置、日射評価方法及び日射評価プログラム |
CN102721988A (zh) * | 2012-06-14 | 2012-10-10 | 中国气象科学研究院 | 基于天空可见光图像的日照时数测量方法 |
CN102721988B (zh) * | 2012-06-14 | 2014-09-03 | 中国气象科学研究院 | 基于天空可见光图像的日照时数测量方法 |
JP2015132539A (ja) * | 2014-01-14 | 2015-07-23 | 株式会社東芝 | 日照量計算装置、経路提案装置および日照量計算方法 |
CN104897130A (zh) * | 2015-06-18 | 2015-09-09 | 广西壮族自治区气象减灾研究所 | 天基遥感分块插值计算太阳高度角的方法 |
CN104897130B (zh) * | 2015-06-18 | 2017-11-14 | 广西壮族自治区气象减灾研究所 | 天基遥感分块插值计算太阳高度角的方法 |
JP2019085056A (ja) * | 2017-11-10 | 2019-06-06 | 東日本旅客鉄道株式会社 | 軌道周辺の環境変化検出方法 |
CN113939710A (zh) * | 2019-05-29 | 2022-01-14 | 古野电气株式会社 | 信息处理系统、方法及程序 |
CN112785678A (zh) * | 2019-11-09 | 2021-05-11 | 李本彦 | 基于三维模拟的日照分析方法和系统 |
CN112785678B (zh) * | 2019-11-09 | 2023-12-12 | 李本彦 | 基于三维模拟的日照分析方法和系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107492069B (zh) | 基于多镜头传感器的图像融合方法 | |
CN110675450B (zh) | 基于slam技术的正射影像实时生成方法及系统 | |
JP5357902B2 (ja) | 撮像装置向き検出装置および当該装置を備える移動体 | |
JP2017175616A (ja) | 撮像装置、画像処理装置および方法 | |
WO2009093587A1 (ja) | オルソフォト画像の生成方法、および撮影装置 | |
JP6398472B2 (ja) | 画像表示システム、画像表示装置、画像表示方法、およびプログラム | |
US20040246333A1 (en) | Digital 3D/360 degree camera system | |
CN113329182A (zh) | 一种图像处理方法、无人机及系统 | |
CN106504192B (zh) | 一种输电线路走廊地质灾害勘探影像处理方法及系统 | |
CN103141078A (zh) | 图像显示装置及图像显示方法 | |
WO2022242395A1 (zh) | 图像处理方法、装置、电子设备和计算机可读存储介质 | |
WO2018198634A1 (ja) | 情報処理装置、情報処理方法、情報処理プログラム、画像処理装置および画像処理システム | |
CN109387186A (zh) | 测绘信息获取方法、装置、电子设备及存储介质 | |
WO2019100219A1 (zh) | 输出影像生成方法、设备及无人机 | |
CN113240615B (zh) | 图像处理方法、装置、电子设备和计算机可读存储介质 | |
JP5846172B2 (ja) | 画像処理装置、画像処理方法、プログラムおよび撮像システム | |
CN112857356B (zh) | 无人机水体环境调查和航线生成方法 | |
CN113450253B (zh) | 图像处理方法、装置、电子设备和计算机可读存储介质 | |
KR101574636B1 (ko) | 면형방식 디지털 항공카메라로 촬영한 시계열 항공사진을 입체시하여 좌표를 연동하고 변화지역을 판독하는 시스템 | |
WO2018214401A1 (zh) | 移动平台、飞行体、支持装置、便携式终端、摄像辅助方法、程序以及记录介质 | |
JP2005331320A (ja) | 天空率および日照時間算出システム、および、算出プログラム | |
CN108195359B (zh) | 空间数据的采集方法及系统 | |
WO2022040868A1 (zh) | 全景拍摄方法、电子设备及存储介质 | |
JP2016053978A (ja) | 画像処理装置 | |
JP2016223934A (ja) | 位置補正システム、位置補正方法、および位置補正プログラム |