JP2005242494A - 画像の数値的解析方法及びシステム並びにプログラム - Google Patents
画像の数値的解析方法及びシステム並びにプログラム Download PDFInfo
- Publication number
- JP2005242494A JP2005242494A JP2004048784A JP2004048784A JP2005242494A JP 2005242494 A JP2005242494 A JP 2005242494A JP 2004048784 A JP2004048784 A JP 2004048784A JP 2004048784 A JP2004048784 A JP 2004048784A JP 2005242494 A JP2005242494 A JP 2005242494A
- Authority
- JP
- Japan
- Prior art keywords
- color
- luminance value
- value
- image
- average luminance
- 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.)
- Granted
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 22
- 239000003086 colorant Substances 0.000 claims abstract description 43
- 230000007423 decrease Effects 0.000 claims abstract description 25
- 238000004891 communication Methods 0.000 claims abstract description 13
- 238000013480 data collection Methods 0.000 claims description 9
- 238000009434 installation Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 7
- 230000003247 decreasing effect Effects 0.000 claims description 6
- 238000000034 method Methods 0.000 claims description 4
- 238000007405 data analysis Methods 0.000 claims description 3
- 230000003203 everyday effect Effects 0.000 abstract description 5
- 238000012935 Averaging Methods 0.000 abstract 1
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 22
- 241000196324 Embryophyta Species 0.000 description 14
- 230000008859 change Effects 0.000 description 13
- 229910002092 carbon dioxide Inorganic materials 0.000 description 11
- 239000001569 carbon dioxide Substances 0.000 description 11
- 238000012544 monitoring process Methods 0.000 description 8
- 230000001932 seasonal effect Effects 0.000 description 8
- 241000973887 Takayama Species 0.000 description 5
- 238000010521 absorption reaction Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 239000002689 soil Substances 0.000 description 5
- 229920003266 Leaf® Polymers 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- 230000005855 radiation Effects 0.000 description 4
- 235000018185 Betula X alpestris Nutrition 0.000 description 3
- 235000018212 Betula X uliginosa Nutrition 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 240000007313 Tilia cordata Species 0.000 description 2
- 206010049155 Visual brightness Diseases 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 238000001454 recorded image Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 241000894007 species Species 0.000 description 2
- 230000007480 spreading Effects 0.000 description 2
- 230000036962 time dependent Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 241000208140 Acer Species 0.000 description 1
- 241000286209 Phasianidae Species 0.000 description 1
- 240000005499 Sasa Species 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 244000269722 Thea sinensis Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000005494 condensation Effects 0.000 description 1
- 238000009833 condensation Methods 0.000 description 1
- 230000002354 daily effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000035613 defoliation Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 235000008216 herbs Nutrition 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000003703 image analysis method Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000002844 melting Methods 0.000 description 1
- 230000008018 melting Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/13—Satellite images
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Remote Sensing (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Facsimile Image Signal Circuits (AREA)
- Color Image Communication Systems (AREA)
- Image Analysis (AREA)
Abstract
【解決手段】 毎日同時刻(正午前後)に撮影し、画素毎に該カラー画像を分解する複数色(RGB)のデジタル画像データとして保存する(ST1)。R、G、Bそれぞれの輝度値について、その画像データに含まれる全画素にわたり平均した値r,g,b(単色の平均輝度値)と全色の輝度値の平均値
I を求め(ST2)、r,g,bそれぞれをIで正規化したrav,gav,bavを求める(ST3)。次に、このrav,gav,bavの時系列を作成し、それぞれの値の1.0に対する大小、増減傾向の正負を求める(ST4)。そして、rav,gav,bavの1.0に対する大小および増減傾向の正負の組み合わせのパターンから、対応する森林の状態を判定する(ST5)。
【選択図】図1
Description
従来のリモートセンシング的手法による画像解析では、高価な放射量測定器を用いて画素単位で測定された特定波長の放射量(輝度)、もしくは放射量を別の指標に変換した値と、その画素を含む物体、例えば樹冠の葉なのか、積雪なのか,を,物体の放射特性から理論的に,もしくは画像の図柄あるいは実地における目視によって判別した結果とを関連付けた対応表(教師空間)を作成し、この関連付けにより画像の全画素を分類する。この分類により、被写体が未知である画素が何を写したものなのかを判別したり、分類ごとの画素数を画像全体にわたって計数して面積を求めたりすることができる。しかし、この手法は煩雑な処理を伴うものである。また、一般向けに市販されているビデオカメラなどの映像機器の出力は、本来は観賞を目的としたものであり、実際の光エネルギー量の絶対値に関する情報は含まれていないので、精密な解析には適さないが,簡易な解析により,有用な情報を取り出せる可能性がある。
そこで、本発明は、カラー画像を単純な数値で表現することによって、低速な通信手段でも、画像から得られるある程度の有用な情報を実時間で収集することを可能にすることを目的としている。
さらに,複雑な画像情報を単純な数値指標で表現することにより,即時に数値演算的な解析への適用を可能にすることを目的としている。
さらに,複雑な画像情報を単純な数値指標で表現することにより,即時に数値演算的な解析に適用することができる。
近年、デジタルスチルカメラなどのデジタル映像機器が安価に入手できるようになり、多点で同時に画像情報を記録することが可能になったことにより、画像情報を記録、蓄積する観測所が増えるものと予想される。精密な物理量測定ではないため、リモートセンシング的な定量解析は困難であるが、本発明のような簡易な解析は可能であり、有用な情報が得られるものと期待される。
植物季節:植生(樹木および草本)の状態の季節変化。たとえば樹木の葉が開く(展葉),紅葉,落葉など。
輝度値: 被写体から放射される光エネルギー量を指す「輝度」に対応して,映像機器が撮影後に記録媒体に整数値として記録した相対的な値。
単色の輝度値:カラーを分解した単色(RGBの場合は赤,緑,青それぞれ)の輝度値
全色の輝度値:カラーの輝度値
(全色の,単色の)平均輝度値:全色の輝度値あるいは単色の輝度値の単一の画像に含まれる全画素にわたる平均値
正規化平均輝度値:各色の単色の平均輝度値について,全色の平均輝度値に基づき正規化した値
まず、撮影記録を行う。図2は、カメラによる撮影を説明する図である。例示のカメラは、高さ約22mの鉄塔上にやや斜め下向きに設置した。より広い範囲を撮影するため、レンズには広角アタッチメントを装着している。画像に収められる範囲は地表で直径およそ40mの円内である。図3に示すように、ビデオカメラの撮像素子で捕らえられた画像は、エンコーダによりアナログのビデオ信号(NTSCコンポジット信号)に変換され、信号線を通じ、塔下に設置した記録装置に送信される。記録装置はパーソナルコンピュータとビデオキャプチャアダプタ(ビデオ信号用アナログ/デジタル変換器)からなり、カメラから受信したビデオ信号をデジタル画像データに変換し、記憶装置(パーソナルコンピュータ)に記録する仕組みになっている。デジタル化した画像データの解像度は横320画素、縦240画素である。なお、この景観カメラシステムは、記憶容量の軽量化を重視し、デジタル化した画像はJPEG形式に圧縮したデータファイルとして記録している。そのため、画像の曖昧化を伴う圧縮によりある程度の誤差(画質の劣化)を伴っている。精密な解析を行うには撮像素子が捕らえた画像を加工せずに記録することが望ましい。
分布を持つものを単一の数値で扱うもっとも簡単な方法は平均値を求めることである。そこで、各日の画像について、R、G、Bの単色それぞれについて輝度値の画像全体にわたる平均値(単色の平均輝度値)の経時変化を求める。R、G、Bの単色の平均輝度値
撮影した(高山サイト)画像では、年間を通じ、全色の平均輝度値が256階調表示で153〜204の範囲でおおむね一定であったが、日によって多少のばらつきがある。この全色の平均輝度値のばらつきを除去し、明るさを一定にした条件下での森林生態系の色の変化を扱うため、
rav、gav、bavの値の時系列により、画像の平均的な色の経時変化を求め,それぞれの値の所定値、例えば1.0に対する大小,増減傾向の正負を求める。画像周辺の遮蔽物の影響を除くため、平均値を求める範囲を図4のように画像中央部の横220、縦200の44000画素を対象として計算を行った。
解析結果を図5に示す。毎日14時に撮影、記録した画像を2001年1月から2002年12月までの2年間にわたって蓄積したものを解析した。各年を図中の上下2つの図に示している。各図中の折れ線は、rav、gav、bavの時系列である。rav、gav、bavを合成して作成した色は、各日の画像全体の平均の色、つまり図5の折れ線を日々について色で表したもの(以下、平均の色という)に相当する。大雑把に見て、冬季にはrav、gav、bavともにほぼ1.0、平均の色は灰色、4月(ユリウス日90日付近)以降に各色が1.0の線を離れ、平均の色はやや赤みを帯びる。6月初め(ユリウス日140日付近)から10月(ユリウス日270日付近)まではgavが最も大きく、平均の色も緑色を帯びている。その後、10月にravが約1ヶ月間卓越して大きく、11月(ユリウス日で300日付近)にはrav、gav、bavが再び1.0の線に収束しており、各年とも同じ季節変化がみられる。直感的には、これらの変化は、冬季の積雪の白色、雪解け後の地面と下草の茶ないし薄茶色、夏季の樹木の葉の緑色、紅葉および落葉の黄色ないし薄茶色と対応していると理解できる。ravは土や枯草、紅葉、落葉の茶系色、gavは下草や樹木の葉の緑色の主な色成分としてわかりやすく、それぞれの状況が発生する季節と正規化平均輝度値の増減が整合していると見ることができる、bav、つまり青色は森林の持つ色成分としては理解しにくい。しかし、bavの減少をその補色である黄色成分の増加として考えると、枯草や樹木の葉、紅葉と落葉に含まれる色として理解しやすい。なお、春から秋の間で、時折bavが0.9〜1.0まで急増することがあるが、これは霧の発生、もしくはカメラのレンズカバーに結露あるいは着霜が生じたことによって画像全体が白くなったためである。今後詳しい検討を行う必要があるが、これを霧の発生の検出に利用できる可能性もある。
I を求め、r,g,bそれぞれをIで正規化したrav,gav,bavを算出する。通信部は、この算出したrav,gav,bavを、送信する。
各日を含む前15日間のrr,rg,rbそれぞれの時間平均値を求める。改めてこれをrr,rg,rbとする。rr,rg,rbを移動平均値として扱うことにより、霧や降雪、撮影時の天気によって発生する1日単位の微細な変動(誤差)を除去し、季節的なトレンドを抽出する。
移動平均値rr,rg,rbそれぞれについて、15日前の移動平均値rr,rg,rbとの差を求める。これをdr,dg,dbとし、15日前におけるrr,rg,rbの増減傾向を表すものとする。
一定のしきい値を設定し、rr,rg,rb,dr,dg,dbの値が正か負か、あるいはゼロであるかを判定する。この例ではrr,rg,rbのしきい値を±0.03,dr,dg,dbのしきい値を±0.003とし、値がこの範囲内にある時はゼロ(0)、この範囲以上であれば正(+)、この範囲以下であれば負(-)とした。
以上6つの変数を3種類に分類した結果(+,-,0)をrr,rg,rb,dr,dg,dbの順に並べ、1つのパターンとする。たとえば、解析結果の図中にあるユリウス日90日に現れたパターンは(+0-+0-)である。なお、ある日に求められるrr,rg,rbは当日から7日前に関する値、dr,dg,dbは当日から15日前に関する値であるので、日付を合わせた上でパターンを作成する必要がある。
各パターンと当該日の画像を目視して判別される森林の状態の対応は図に示す通りである。同じパターンが繰り返し現れたパターン12からパターン19はいずれも樹木に夏の葉がついている状態であり、ひとまとまりのものとして扱って差しつかえない。
Claims (3)
- 撮影した複数画素からなるカラー画像を、画素毎にこのカラーを複数色に分解したデジタル画像データとして画素毎に保存し、
分解した色毎の単色の輝度値及び全色の輝度値について、その画像データに含まれる全画素にわたり平均した単色の平均輝度値及び全色の平均輝度値を求め、かつ、この全色の平均輝度値に基づき正規化した色毎の正規化平均輝度値を求め、
色毎に正規化平均輝度値の時系列を作成し、それぞれの値の所定値に対する大小及び増減傾向の正負を求め、
色毎の正規化平均輝度値の前記所定値に対する大小及び増減傾向の正負の組み合わせのパターンから、対応するカラー画像の状態を判定することから成る画像の数値的解析方法。 - 現場設置モジュールと、データ集積・解析モジュールとから成る画像の数値的解析システムにおいて、
前記現場設置モジュールは、複数画素からなるカラー画像を撮影して、画素毎にこのカラーを複数色に分解したデジタル画像データに変換する撮影部と、
分解した色毎の単色の輝度値及び全色の輝度値について、その画像データに含まれる全画素にわたり平均した単色の平均輝度値及び全色の平均輝度値を求め、かつ、この全色の輝度値の平均値に基づき正規化した色毎の正規化平均輝度値を算出するデータ収集及び演算部と、
この算出した色毎の正規化平均輝度値を送信する通信部とを有し、
前記データ集積・解析モジュールは、前記現場設置モジュールより送信された色毎の正規化平均輝度値を受信する通信部と、
色毎の正規化平均輝度値の時系列を作成し、かつ、それぞれの値の所定値に対する大小及び増減傾向の正負を求めて、色毎の正規化平均輝度値の前記所定値に対する大小及び増減傾向の正負の組み合わせのパターンから、対応するカラー画像の状態を判定するデータ収録及び演算部と、を有する、
ことから成る画像の数値的解析システム。 - 撮影した複数画素からなるカラー画像を、画素毎にこのカラーを複数色に分解したデジタル画像データとして画素毎に保存し、
その画像データ毎に含まれる全画素にわたり平均した単色の平均輝度値及び全色の平均輝度値を求め、かつ、この全色の平均輝度値に基づき正規化した色毎の正規化平均輝度値を求め、
色毎の正規化平均輝度値の時系列を作成し、それぞれの値の所定値に対する大小及び増減傾向の正負を求め、
色毎の正規化平均輝度値の前記所定値に対する大小及び増減傾向の正負の組み合わせのパターンから、対応するカラー画像の状態を判定する各手順を実行する画像の数値的解析プログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2004048784A JP4280823B2 (ja) | 2004-02-24 | 2004-02-24 | 画像の数値的解析方法及びシステム並びにプログラム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2004048784A JP4280823B2 (ja) | 2004-02-24 | 2004-02-24 | 画像の数値的解析方法及びシステム並びにプログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2005242494A true JP2005242494A (ja) | 2005-09-08 |
JP4280823B2 JP4280823B2 (ja) | 2009-06-17 |
Family
ID=35024199
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2004048784A Expired - Lifetime JP4280823B2 (ja) | 2004-02-24 | 2004-02-24 | 画像の数値的解析方法及びシステム並びにプログラム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4280823B2 (ja) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016518098A (ja) * | 2013-04-17 | 2016-06-20 | ビュルト エレクトロニク アイソス ゲゼルシャフト ミット ベシュレンクテル ハフツング ウント コンパニー コマンディトゲゼルシャフト | 通信装置 |
JP2016111939A (ja) * | 2014-12-11 | 2016-06-23 | 富士通株式会社 | 耕作放棄地判定方法、耕作放棄地判定プログラム及び耕作放棄地判定装置 |
JP2020161353A (ja) * | 2019-03-27 | 2020-10-01 | パナソニックIpマネジメント株式会社 | 照明システム、及び、照明方法 |
-
2004
- 2004-02-24 JP JP2004048784A patent/JP4280823B2/ja not_active Expired - Lifetime
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016518098A (ja) * | 2013-04-17 | 2016-06-20 | ビュルト エレクトロニク アイソス ゲゼルシャフト ミット ベシュレンクテル ハフツング ウント コンパニー コマンディトゲゼルシャフト | 通信装置 |
US9973046B2 (en) | 2013-04-17 | 2018-05-15 | Würth Elektronik eiSos Gmbh & Co. KG | Communication device |
JP2016111939A (ja) * | 2014-12-11 | 2016-06-23 | 富士通株式会社 | 耕作放棄地判定方法、耕作放棄地判定プログラム及び耕作放棄地判定装置 |
JP2020161353A (ja) * | 2019-03-27 | 2020-10-01 | パナソニックIpマネジメント株式会社 | 照明システム、及び、照明方法 |
JP7340788B2 (ja) | 2019-03-27 | 2023-09-08 | パナソニックIpマネジメント株式会社 | 照明システム、及び、照明方法 |
Also Published As
Publication number | Publication date |
---|---|
JP4280823B2 (ja) | 2009-06-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Winkler et al. | Assessing the effects of post‐pine beetle forest litter on snow albedo | |
US20190244022A1 (en) | Plant Stand Counter | |
US5130545A (en) | Video imaging plant management system | |
Toda et al. | Estimation of plant area index and phenological transition dates from digital repeat photography and radiometric approaches in a hardwood forest in the Northeastern United States | |
Inoue et al. | Utilization of ground-based digital photography for the evaluation of seasonal changes in the aboveground green biomass and foliage phenology in a grassland ecosystem | |
Abd Latif et al. | Determination of tree species using Worldview-2 data | |
Pitt et al. | Use of remote sensing for forest vegetation management: a problem analysis | |
Sakamoto et al. | Application of day and night digital photographs for estimating maize biophysical characteristics | |
Lang et al. | Digital photography for tracking the phenology of an evergreen conifer stand | |
Lim et al. | Ecological consideration for several methodologies to diagnose vegetation phenology | |
Soudani et al. | A survey of proximal methods for monitoring leaf phenology in temperate deciduous forests | |
AU2021204034B2 (en) | Information processing device, information processing method and program | |
Peltoniemi et al. | Webcam network and image database for studies of phenological changes of vegetation and snow cover in Finland, image time series from 2014 to 2016 | |
JP2010220569A (ja) | 植物体の生育変化の連続測定方法 | |
Zakaluk et al. | Predicting the leaf water potential of potato plants using RGB reflectance. | |
JP4280823B2 (ja) | 画像の数値的解析方法及びシステム並びにプログラム | |
Leblanc et al. | Recent advancements in optical field leaf area index, foliage heterogeneity, and foliage angular distribution measurements | |
Coops et al. | Procedures for predicting habitat and structural attributes in eucalypt forests using high spatial resolution remotely sensed imagery | |
US20230408889A1 (en) | Imaging apparatus and lens apparatus | |
Law | Estimation of leaf area index and light intercepted by shrubs from digital videography | |
Xiang et al. | Monitoring mangrove phenology using camera images | |
JP2004213627A (ja) | 植物活力変動の評価画像作成方法 | |
Everitt et al. | Light reflectance characteristics and video remote sensing of pricklypear. | |
Dewey et al. | A comparison of tools for remotely estimating leaf area index in loblolly pine plantations | |
Wilkinson et al. | Can upward-facing digital camera images be used for remote monitoring of forest phenology? |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20051004 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20090204 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20090210 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4280823 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
EXPY | Cancellation because of completion of term |