(57)【要約】 【課題】光源からの距離が異なる被写体や背景が含まれるフリッカ成分が画像信号内で分布を持っている場合に対しても、適切なフリッカ検出枠を設定し、検出・補正に関わるスループットを下げることなく、精度良くフリッカ成分を検出することが可能な画像処理装置を提供すること。【解決手段】フィールド面内の垂直方向における、照明強度の変化情報を抽出するフリッカ検出手段と、検出した前記照明強度の変化情報に基づき、前記フリッカ検出手段で検出されたフリッカ成分の逆特性（逆数）を、フリッカ補正値として算出するフリッカ補正値算出手段と、前記フリッカ補正値算出手段で算出された補正値を画像信号に乗算し、フリッカ補正するフリッカ補正手段とを備え、距離画像を取得する距離画像取得手段と、その距離画像に基づいて前記フリッカ検出手段の検出枠と、前記フリッカ補正値生成手段の補正枠を設定する手段と、を有することを特徴とする。  【選択図】図１
2016122941.tif 000002
【特許請求の範囲】
【請求項１】
フィールド面内の垂直方向における、照明強度の変化情報を抽出するフリッカ検出手段と、検出した前記照明強度の変化情報に基づき、前記フリッカ検出手段で検出されたフリッカ成分の逆特性（逆数）を、フリッカ補正値として算出するフリッカ補正値算出手段と、前記フリッカ補正値算出手段で算出された補正値を画像信号に乗算し、フリッカ補正するフリッカ補正手段とを備え、距離画像を取得する距離画像取得手段と、その距離画像に基づいて前記フリッカ検出手段の検出枠と、前記フリッカ補正値生成手段の補正枠を設定する手段と、を有することを特徴とする画像処理装置。
【請求項２】
前記距離画像から同一距離における異なる被写体を検出する被写体検出手段と前記距離画像と前記被写体検出手段に基づいて前記フリッカ検出手段の検出枠と、前記フリッカ補正値生成手段の補正枠を設定する手段と前記フリッカ補正手段から各々の枠に対してフリッカ補正を行うことを特徴とする請求項１に記載の画像処理装置。
【請求項３】
フィールド面内の垂直方向における、照明強度の変化情報を抽出するフリッカ検出手段と、検出した前記照明強度の変化情報に基づき、前記フリッカ検出手段で検出されたフリッカ成分の逆特性（逆数）を、フリッカ補正値として算出するフリッカ補正値算出手段と、前記フリッカ補正値算出手段で算出された補正値を画像信号に乗算し、フリッカ補正するフリッカ補正手段とを備え、距離画像を取得する距離画像取得手段と、その距離画像に基づいて前記フリッカ検出手段の検出枠と、前記フリッカ補正値生成手段の補正枠を設定する手段と、を有することを特徴とする画像処理方法。
【請求項４】
前記距離画像から同一距離における異なる被写体を検出する被写体検出手段と前記距離画像と前記被写体検出手段に基づいて前記フリッカ検出手段の検出枠と、前記フリッカ補正値生成手段の補正枠を設定する手段と前記フリッカ補正手段から各々の枠に対してフリッカ補正を行うことを特徴とする請求項３に記載の画像処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、フリッカ検出装置及びフリッカ検出方法に関し、特に、ＸＹアドレス方式の撮像素子を用いて撮像された画像中のフリッカ成分を検出、補正する画像処理装置及び画像処理方法に関する。
【背景技術】
【０００２】
ＣＭＯＳイメージセンサに代表されるＸＹアドレス方式の個体撮像素子を用いた撮像装置において、蛍光灯照明下の被写体を動画撮影する場合、撮像画像に例えば図６に示すように垂直方向の周期的な明暗が現れることがある。この現象は蛍光灯フリッカと呼ばれ、商用電源周波数と撮像装置の垂直周期周波数との違いによるものである。
【０００３】
この蛍光灯フリッカの発生原理を図７を用いて説明する。図７は商用交流電源周波数が５０Ｈｚの地域における蛍光灯照明下で、ＣＭＯＳイメージセンサを利用した撮像装置で被写体を動画として撮影した場合の、撮像素子における電荷蓄積の様子を時間方向に模式的に表したものである。７０１は撮像装置の撮影方式で決定される垂直同期信号ＶＤ、７０２は世界各国および地域で各々決められている商用電源の周期（ここでは１／５０ｓ）である。また、７０３は商用電源の周期７０２に応じて蛍光発光している発光体の発行周期を示している。
【０００４】
ここで、ＣＭＯＳイメージセンサ上での電荷の蓄積時間が、垂直同期信号ＶＤの周期に従って垂直同期信号の周期と同じ時間（１／６０ｓ）間隔であるものとする。７０４に示したようにＣＭＯＳイメージセンサの場合、各画素の露光タイミングはライン毎に順次制御される為、単色の平面を撮影すると１ライン目の蓄積電荷量と２ライン目の蓄積電荷量とは異なる。このように同一フレーム内においても電荷蓄積時間内に受光する総光量がライン毎に異なる為、蛍光灯フリッカが発生する。
【０００５】
そこで、この問題を解決する為、フリッカの除去を行う方法が種々提案されており、特許文献１では、ＣＭＯＳイメージセンサの出力を水平方向に積分して垂直強度分布を求め、複数フレームの垂直強度分布から現フレームにおける垂直方向のフリッカ成分を算出する。そして、算出したフリッカ成分から補正係数を算出して該補正係数を現フレームの映像信号に乗じることにより、原画像（補正前の撮影画像）を補正している。
【０００６】
また、撮像素子の撮像領域の一部を、撮像以外の機能に使用しているものがある。特許文献２では、撮像素子の一部の画素の受光部を２分割することで瞳分割機能を付与している。そして、これらの画素を焦点検出用画素とし、撮像用画素群の間に所定の間隔で配置することで、位相差式焦点検出を可能としている。
【０００７】
また、特許文献３では、マイクロレンズとフォトダイオードの間に遮光層をもうけることで瞳分離を行う焦点検出用画素と通常の画素を併せ持つ撮像素子とする事で撮像装置でありながら焦点検出の機能を持つ構成としている。特許文献３では、左右または上下に瞳分離された焦点検出用画素からの信号のみを集めて瞳分割された２つの線像を形成し、その像のずれ量に基づいてデフォーカス量を算出する。一定のデフォーカス量において瞳分割された像の投影位置のずれ量は分割された２つの瞳間の基線長により決まる。基線長が長いほどずれ量は大きくなり焦点検出の精度が向上する。
【０００８】
また、特許文献４では、デフォーカス量からピント位置に対する被写体位置のずれの大きさ、すなわちピント位置を基準とした被写体までの距離に関する量について記載している。デフォーカス量に基づくデフォーカス情報を取得することにより、各焦点検出領域に対応する被写体までの距離に関する距離情報を取得し、その距離情報をメモリなどの記録媒体に記録する構成としている。
【０００９】
しかしながら、特許文献１記載の方法は、現在及び過去のフレーム（またはフィールド）の画像信号を水平方向に積分して生成した垂直方向の強度分布を用いてフリッカ成分を検出する方法であるため、例えば光源に対し奥に位置する被写体と、光源に対し手前に位置する被写体とで異なっている場合に、フリッカ成分が異なるため、正確にフリッカ成分を検出することができないという問題がある。この問題は、例えば図８に示す構図のような場合、図下にある撮像装置から撮影した時に、画面左側の人間は光源の手前に位置するため、逆光となりフリッカの影響を受けにくく、画面右側の人間は光源の奥に位置するため、フリッカの影響を受けるような場合に発生しうる。
【００１０】
また特許文献５では、複数の領域に分割してフリッカの検出枠を設定し、フリッカ補正を行うとしている。領域分割の最小単位として、１画面分の前記画像信号の水平方向において、各々が少なくとも１画素の画像信号を含む複数のフリッカ検出枠を設定すると記載している。
【００１１】
領域の分割に関して、分割数を増やせば増やすほどより精度の高いフリッカ検出が出来るとしているが、その分検出に使用する演算量も増えてしまい、更に分割数ごとに異なる補正ゲイン値を持つ為、使用するメモリの量が増えてしまい、結果スループットの低下と必要なメモリ量が増えるという課題がある。
【先行技術文献】
【特許文献】
【００１２】
【特許文献１】特開平１１−１２２５１３号公報
【特許文献２】特開２０００−２９２６８６号公報
【特許文献３】特開２０００−１５６８２３号公報
【特許文献４】特開２００９−５５０９３号公報
【特許文献５】特開２００９−１７２１４号公報
【発明の概要】
【発明が解決しようとする課題】
【００１３】
本発明は、このような従来技術の問題点に鑑みなされたものである。本発明の目的は、光源からの距離が異なる被写体や背景が含まれるフリッカ成分が画像信号内で分布を持っている場合に対しても、適切なフリッカ検出枠を設定し、検出・補正に関わるスループットを下げることなく、精度良くフリッカ成分を検出することが可能な画像処理装置及び画像処理方法を提供することにある。
【課題を解決するための手段】
【００１４】
上記の目的を達成する為に、本発明の画像処理装置は、撮像素子により得られた画像信号中において前述した焦点検出画素を用いて距離情報を取得し、その距離情報に基づいて複数のフリッカ検出枠を設定する設定手段と、画像信号に垂直方向の輝度レベルの変動として生じるフリッカ成分を複数のフリッカ検出枠の各々で検出するフリッカ検出手段と、フリッカ検出手段で算出された補正ゲイン値を画像信号に乗ずる事によってフリッカを補正するフリッカ補正手段とを有することを特徴とする。
【００１５】
また、上記の目的は、撮像素子により得られた画像信号中のフリッカ成分を検出するフリッカ検出方法であって、前述した焦点検出画素を用いて距離情報を取得し、その距離情報に基づいて複数のフリッカ検出枠を設定する設定工程と、画像信号に垂直方向の輝度レベルの変動として生じるフリッカ成分を複数のフリッカ検出枠の各々で検出するフリッカ検出工程と、フリッカ検出工程で算出された補正ゲイン値を画像信号に乗ずる事によってフリッカを補正するフリッカ補正工程とを有することを特徴とする画像処理方法によっても達成される。
【発明の効果】
【００１６】
本発明によれば、光源からの距離が異なる被写体や背景が含まれるフリッカ成分が画像信号内で分布を持っている場合や現フレームまたは現フィールドの垂直方向の強度分布に変動があった場合でも、精度良くフリッカ成分を検出することが可能な画像処理装置及び画像処理方法の提供を実現できる。
【図面の簡単な説明】
【００１７】
【図１】本発明の第１の実施形態における撮像装置の構成例を示した図である。
【図２】本発明の第１の実施形態におけるフリッカ検出、補正回路の構成例を示すブロック図である。
【図３】本発明の第１の実施形態におけるフリッカ検出部の構成例を示すブロック図である。
【図４】焦点検出用画素と通常の画素とを兼用している撮像素子の配置図の一例である。
【図５】縦に瞳分割を行った時の２つの像の図である。
【図６】ＸＹアドレス方式の固体撮像素子を用いた場合に発生するフリッカ現象を表す図である。
【図７】蛍光灯フリッカの発生原理を示した図である。
【図８】光源の手前と光源の奥に被写体が存在する画像を模した図である。
【図９】フリッカ検出、補正の動作の一例を示したフローチャートである。
【図１０】図８の画像に対して距離情報を算出し、その距離情報を画像に置き換えたものである。
【図１１】撮像装置からみて同一距離に異なる被写体が存在した時の図である。
【図１２】図１１の画像に対して距離情報を算出し、その距離情報を画像に置き換えたものである。
【図１３】同一距離で不連続な距離情報を含んだ１ラインのみの距離情報である。
【発明を実施するための形態】
【００１８】
以下、本発明の好適な実施形態について、図面を参照しながら詳細に説明する。
【００１９】
［実施例１］
図１は本発明の第１の実施形態における撮像装置の概略構成を示すブロック図である。被写体を撮影時、１０１の撮像手段は光電変換により光学情報を電気信号に換えて出力する。１０２のＣＤＳ（相関二重サンプリング）／ＡＧＣ（Ａｕｔｏ  Ｇａｉｎ  Ｃｏｎｔｒｏｌ）は、アナログの電気信号のＳ／Ｎを改善するための処理部であり、１０３のＡ／Ｄ変換手段によって、撮像電気信号は、デジタル画像信号に変換される。
【００２０】
１０９はシェーディング補正部であり、シェーディングによって発生する２像のレベル差を補正する。１０４はフリッカ検出、補正回路であり変換されたデジタル画像信号に対して距離情報を用いてフリッカの検出、補正を行う。また、１０４内の乗算器１０５は変換されたデジタル画像信号に対して算出された補正ゲイン値を乗算し、フリッカ補正を行う。
【００２１】
次に相関演算部１０６において相関演算を行い、デフォーカス量を算出し距離情報算出部１０８において距離情報を算出する。算出された距離情報は１０４フリッカ検出、補正回路内の制御部に送り、その距離情報に基づいてフリッカの検出枠を設定する。フリッカ補正を行った後、現像系データは画素加算１０７を行い、ベイヤ配列に変換し現像処理のブロックへデータを出力する。
【００２２】
図２は本発明の第１の実施形態に係るフリッカ検出装置の一例としての図１におけるフリッカ検出、補正回路１０４の構成例を示すブロック図である。図２において、フリッカ検出部２０１は、ＣＭＯＳイメージセンサのような、ＸＹアドレス方式の撮像素子によって撮像された画像信号からフリッカ成分を検出する。本実施形態において、画像信号は、撮像素子を一ライン毎に垂直方向に順次リセットすることによって順次蓄積が開始され、一ライン毎に順次水平走査することにより読みだされるものとする。フリッカ補正値生成部２０２はフリッカ検出部２０１で検出されたフリッカ成分の逆特性（逆数）を、フリッカ補正値として求め、画素単位のフリッカ補正値を生成する。
【００２３】
そして、乗算器２０３は、フリッカ補正値生成部２０２で得られたフリッカ補正値と画像信号を乗算し、補正後の画像信号として出力する。制御部２０４は、前述した距離情報を取得し、フリッカ検出・補正回路１０４全体の動作を制御する。また、制御部２０４は、距離情報からフリッカ検出部２０１で設定するべき検出枠数を算出し、画像信号内の画素がどの検出枠に該当するかを制御する。その制御情報をフリッカ検出部２０１に送り、フリッカ検出部２０１は前記制御情報に従い、各枠の検出値を算出する。
【００２４】
図３は図２におけるフリッカ検出部２０１の構成例を示すブロック図である。図３において、ブロック分割部３０１は、フリッカ成分の検出を行う領域（フリッカ検出枠）を、１画面分の画像信号に対して水平方向および垂直方向に設定する。加算平均部３０２は、ブロック分割部３０１が設定した検出領域の各々について、画素の輝度レベルの平均値を算出する。また、引算器３０３、乗算器３０４、加算器３０５およびメモリ３０６は、以下に示す演算を行うことで、いわゆる巡回型ローパスフィルタ３０８を形成している。
【００２５】
mem = ave × k + mout×(1-k)
ここで、aveは加算平均部３０１の出力を示し、moutはメモリ３０５からの出力を示す。memは加算器３０４からの出力を示し、メモリ３０５に新しく記憶される値を示す。また、kは巡回型ローパスフィルタ３０８のフィルタ係数である。除算器３０６は、加算平均部３０１の出力と、メモリ３０５からの出力を除算することによってフリッカ検出枠ごとのフリッカ成分を算出して出力する。つまり、水平方向の各々の検出領域（垂直方向の複数の検出領域からなる）では、垂直方向の輝度レベルの変動として生じるフリッカ成分が検出される。
【００２６】
次に、図４は焦点検出用画素と通常の画素とを兼用している撮像素子の配置図の一例である。４０１は４つの撮像素子とその色フィルタの前面に配置されるマイクロレンズである。このマイクロレンズ４０１は、明るい交換レンズの射出瞳（例えばF1.0）を通過する光速を全て受光するような形状に設計されている。この４つの撮像素子に各々図４のようにＡ、Ｂ、Ｃ、Ｄと割り当てていくと縦に瞳分割を行う際は、Ａ＋Ｂ、Ｃ＋Ｄと分ける事が出来る。
【００２７】
また、４０２の４つの撮像素子はＲ（赤）、４０３、４０４の各々４つの撮像素子はＧ（緑）、４０５の４つの撮像素子はＢ（青）の色フィルタが配置されている。
【００２８】
このマイクロレンズを配置した撮像素子において縦に瞳分割を行い、それを像にしたものが図５である。５０１がＡ＋Ｂ像、５０２がＣ＋Ｄ像であるとすると、ピントの合っている位置５０３においては像ずれが発生しない。
【００２９】
一方、ピントの合っていない位置５０４においては縦に瞳分割を行った為、縦方向に像ずれが発生している。
【００３０】
前述したずれ量から相関演算を行い、デフォーカス量を算出する事によって距離情報を取得する。
【００３１】
図１０は前述した図８に対して距離情報を取得し、それを画像情報に変換したものである。ここでは、一例として撮像装置からの距離手前の被写体の方が輝度が高く、奥の被写体の方が輝度が低いとしているが必ずしもこのような図になるとは限らない。
【００３２】
次に、本実施形態のフリッカ補正装置の動作の一例について、距離情報算出とフリッカ検出、補正に至るまでの手順をフローチャートに示したのが図９である。
【００３３】
Ｓ９０１によってシェーディングによって発生する２像のレベル差を補正する。シェーディング補正は光学系の特性による輝度ムラに対して、一様な明るさの画像になるように補正するものでセンサの特性によるフリッカよりも前に補正する。
【００３４】
Ｓ９０２で相関演算を行う。具体的には２つの像を重ね合わせて差の絶対値の和を求める処理を、重ね合わせ位置をずらしながら行い、最も差の絶対値が小さくなる位置を求める。相関演算結果を使用してＳ９０３において差の絶対値が小さくなる位置の前後の値からパラボラフィッティングという手法を用いてサブピクセルマッチングを行う。これにより像のずれ量の精度を向上させる。
【００３５】
最後に、Ｓ９０４で像ずれ量に基線長で決まる所定の係数を乗算する事でデフォーカス量が求まる。なお、デフォーカス量は、撮影レンズのピント位置に対する被写体位置のずれの大きさ、すなわちピント位置を基準とした被写体までの距離に関する量を表している。したがって、デフォーカス量に基づくデフォーカス情報を取得することにより、各焦点検出領域に対応する被写体までの距離に関する距離情報を取得するが出来る。すなわちこれが、Ｓ９０５であり、この距離情報をメモリにより一時的に記憶させる。
【００３６】
次に、Ｓ９０５で求められた距離情報を使用し、フリッカ検出枠を設定する。具体的には得られた距離情報の距離範囲（ｘｘ〜ｘｘ）から各画素がどの距離範囲に入っているかを見てその距離範囲毎に検出枠が存在し、従来のフリッカ検出と同様に積分・正規化を行う。
【００３７】
次に、Ｓ９０７では、前述したフリッカ検出と同様のフリッカ検出を行う。次に、Ｓ９０８では、得られたフリッカ検出値の逆数を取るようにし、フリッカ補正ゲイン値を算出する。ここで、得られる補正ゲイン値は、前述した距離範囲毎の検出枠別に算出されるものである。フリッカ検出、フリッカ補正ゲイン値算出は枠毎に行う。すなわち前述した距離情報から指定範囲により枠数が求まり、その枠数分だけのフリッカ検出、フリッカ補正ゲイン値が算出される。
【００３８】
Ｓ９０９では最後の検出枠かどうかの判断を行い、最後の検出枠で無ければ次の枠のフリッカ検出、補正ゲイン値算出を行う。最後の検出枠であれば全ての枠に対して検出、補正ゲイン値が算出された事になり、最後にＳ９１０では得られたフリッカ補正ゲイン値で画像に対してフリッカ補正を行う。また、上記フリッカ補正と並行して、Ｓ９１１ではデフォーカス量を用いてフォーカス処理を行う。
【００３９】
以上、本実施形態によれば、入力された画像信号の水平方向及び垂直方向の両方について、複数の領域でフリッカ成分を検出する。その為、フリッカ成分が画像信号内で分布を持っている場合や、現フレームまたは現フィールドの垂直方向の強度分布に変動があった場合でも、精度良くフリッカ成分を検出、補正することが出来る。
【００４０】
［実施例２］
以下、本発明を実施する為の最良の形態について、図面と共に説明する。図１１は撮像装置からの距離が同一な二人の被写体に対し、光源左に位置する人に対しては太陽光の光はあまり受けず、蛍光灯のみの光しか受けない。また、光源右に位置する人は蛍光灯と太陽光の両方の光を受ける。光源は蛍光灯と太陽光の二つ存在する構図であり、二人の被写体はフリッカ成分が異なっている事を示したものである。
【００４１】
また、図１２は前述した図１１に対して距離情報を取得し、それを画像情報に変換したものである。ここでは、一例として撮像装置からの距離手前の被写体の方が輝度が高く、奥の被写体の方が輝度が低いとしているが必ずしもこのような図になるとは限らない。
【００４２】
図１３は図１２の画像情報内の、被写体二人を含んだある１ラインを拡大して見た時の図である。ここで、１３０１は光源左側に位置する人の距離情報である。１３０２は光源の奥に位置する壁に対しての距離情報である。１３０３は光源右側に位置する人の距離情報である。
【００４３】
この距離情報から前述した制御部２０４で、ある一定の間隔を隔てて同一距離の情報が入力されたと判断し、フリッカ検出部２０１に異なる枠設定を行うよう指示する。
【００４４】
具体的には制御部２０４は、ライン単位で距離情報１が入力されてきたら距離フラグ１を立てておき（’1’の値）、次に距離情報２が入力されてきたら距離フラグ２を立て（’1’の値）、その情報を見て距離フラグ１に対する距離カウント１のカウントを行う。ここである値までカウントされ、再び距離情報１が入力された時、距離カウント１と距離フラグ１の値を見て、同一距離の異なる被写体が入力されてきたと判断する。
【００４５】
前述したのは制御部２０４において、距離情報から両者間が離れているか離れていないかを検出するための手段、方法である。
【００４６】
すなわち、図１２に示しているようにこの画像において、検出枠を３つ設定する事により、同一距離における異なる被写体それぞれに、フリッカ検出、補正を行う事が出来る。
【００４７】
以上のようにして同一距離であっても各被写体が離れていると判断された場合は、検出枠を増やす方向に設定を行い、異なる検出枠を設定する事によって、複数光源や光源からの距離が異なる被写体に対して精度の良いフリッカ検出、補正を行う事が出来る。また、適正なフリッカ検出枠を設定する事によりフリッカ検出、補正に関わる演算量と必要なメモリ量を減らすことが出来、画像処理に関わるスループットを低下させることなくフリッカ検出、補正を行うことが出来る。
【符号の説明】
【００４８】
１０１  撮像手段、１０２  ＣＤＳ／ＡＤＣ、１０３  Ａ／Ｄ変換、
１０４  フリッカ検出・補正回路、１０５，２０３  乗算器、１０６  相関演算部、
１０７  画素加算部、１０８  距離情報算出部、１０９  シェーディング補正部、
２０１  フリッカ検出部、２０２  フリッカ補正値生成部、２０４  制御部、
３０１  ブロック分割部、３０２  加算平均部、３０３  減算器、３０４  乗算器、
３０５  加算器、３０６  一時的なメモリ（バッファ）、３０７  除算器、
３０８  巡回型ローパスフィルタ、４０１  マイクロレンズ、
４０２  Ｒ（赤）のカラーフィルタ、４０３，４０４  Ｇ（緑）のカラーフィルタ、
４０５  Ｂ（青）のカラーフィルタ、
５０１  縦に瞳分割した時のマイクロレンズ上部の画素から生成した像（Ａ＋Ｂ像）、
５０２  縦に瞳分割した時のマイクロレンズ下部の画素から生成した像（Ｃ＋Ｄ像）、
５０３  ピントの合っている像、５０４  ピントの合っていない像、
７０１  垂直同期信号ＶＤ、７０２  商用電源の周期（ここでは１／５０ｓ）、
７０３  商用電源の周期７０２に応じて蛍光発光している発行体の周期、
７０４  各ラインの蓄積時間と蓄積電荷量、１３０１  光源左側の被写体の距離情報、
１３０２  光源奥の壁の距離情報、１３０３  光源右側の被写体の距離情報
【図１】
2016122941.tif 000003
【図２】
2016122941.tif 000004
【図３】
2016122941.tif 000005
【図４】
2016122941.tif 000006
【図５】
2016122941.tif 000007
【図６】
2016122941.tif 000008
【図７】
2016122941.tif 000009
【図８】
2016122941.tif 000010
【図９】
2016122941.tif 000011
【図１０】
2016122941.tif 000012
【図１１】
2016122941.tif 000013
【図１２】
2016122941.tif 000014
【図１３】
2016122941.tif 000015