(57)【要約】【課題】農作物の内部の状態を非破壊、非接触で計測すること。【解決手段】中心波長１７００ｎｍ、半値幅２４２ｎｍ、帯域１４００ｎｍ〜１９５０ｎｍのスーパーコンティニューム光（ＳＣ光）を用いて、米のＯＣＴ計測を行った。図７（ａ）は、米の外観を示した光学顕微鏡写真であり、図７（ｂ）は、ＯＣＴ画像（深さ方向の断面画像であり、（ａ）におけるｘｚ面での断面画像）である。このＯＣＴ画像の濃淡は、デンプンの濃度や結晶化度を反映していると考えられる。このように、ＯＣＴ画像の濃淡によって、米の内部の状態を非破壊、非接触で観測することができる。【選択図】図７
2017181300.tif 000002
【特許請求の範囲】
【請求項１】
農作物の内部の状態を非破壊で観測する農作物観測方法であって、
中心波長が０．７〜２．２μｍのスーパーコンティニューム光を用いたＯＣＴ計測を行うことにより、前記農作物の内部の状態を観測する、
ことを特徴とする農作物観測方法。
【請求項２】
前記農作物内部の糖類の濃度分布を計測する、ことを特徴とする請求項１に記載の農作物観測方法。
【請求項３】
前記農作物内部のデンプンの濃度分布を計測する、ことを特徴とする請求項２に記載の農作物観測方法。
【請求項４】
前記農作物は米であり、デンプンの濃度分布の計測により登熟度を評価する、ことを特徴とする請求項３に記載の農作物観測方法。
【請求項５】
スーパーコンティニューム光として、中心波長が１．３〜２．０μｍのものを用いて、デンプンの濃度分布を計測する、ことを特徴とする請求項４または請求項５に記載の農作物観測方法。
【請求項６】
前記農作物内部の糖類の濃度分布は、多変量解析により推定する、ことを特徴とする請求項２ないし請求項５のいずれか１項に記載の農作物観測方法。
【請求項７】
前記農作物の内部の細胞壁の状態を計測する、ことを特徴とする請求項１に記載の農作物観測方法。
【請求項８】
前記農作物の葉の葉脈の状態を計測する、ことを特徴とする請求項１に記載の農作物観測方法。
【請求項９】
前記ＯＣＴ計測は、周波数領域で行うものであり、合わせて前記農作物の分光情報を取得する、ことを特徴とする請求項１ないし請求項８のいずれか１項に記載の農作物観測方法。
【請求項１０】
請求項１ないし請求項９のいずれか１項に記載の農作物観測方法を用いて、農作物の内部の状態を観測し、その観測結果を元に、前記農作物の育成方法を調整して品質評価を行いながら農作物を育成する農作物の生産方法。
【請求項１１】
中心波長が０．７〜２．２μｍのスーパーコンティニューム光を生成する光源と、
前記光源からのスーパーコンティニューム光を分岐して、農作物に照射して反射光させる測定光を得るための測定光系と、
その反射光との干渉光を得るための参照光を得る参照光系と、
前記測定光と前記参照光との干渉光を得て、農作物の内部の状態を反映した干渉画像を得る画像生成装置と、
を有することを特徴とする農作物観測装置。
【請求項１２】
前記参照光系と前記測定光系とを、手持ちが可能な筐体に一体に納めたハンディプローブを有する、ことを特徴とする請求項１０に記載の農作物観測装置。
【請求項１３】
前記農作物観測装置は、周波数領域で前記干渉光の波形を取得する装置である、ことを特徴とする請求項１２に記載の農作物観測装置。
【請求項１４】
前記農作物の分光情報を得る手段をさらに有する、ことを特徴とする請求項１３に記載の農作物観測装置。
【請求項１５】
中心波長が０．７〜２．２μｍのスーパーコンティニューム光を生成する光源であって、ＯＣＴ計測によって農作物の内部の状態を非破壊で観測するためのスーパーコンティニューム光源。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、農作物の内部の状態を非破壊で計測することができる農作物観測方法および農作物観測装置である。また、その農作物観測方法を用いた農作物の生産方法である。
【背景技術】
【０００２】
農作物の糖度などを非破壊で計測する方法として、農作物に近赤外線を照射し、その反射スペクトルや透過スペクトルを測定して分析する方法が広く知られている。特に、果物の糖度を分析する方法として広く用いられている。また、果物を一定の厚さにスライスして、そのスライスしたものに近赤外線を照射してスペクトル測定することで、果物の３次元的な糖度分布を計測することも行われている。
【０００３】
また、光の干渉を利用して対象物内部の表面近傍の断層画像を得る光干渉断層法（Optical Coherence Tomography；ＯＣＴ）が知られている。ＯＣＴは、医療分野において広く活用されており、たとえば眼底検査において網膜の３次元画像を得ることに利用されている。特許文献１、２には、スーパーコンティニューム光を用いたＯＣＴが記載されている。スーパーコンティニューム光は、超短パルス光を高非線形光ファイバーに導入することにより得られる波長帯域が数百ｎｍの広帯域光である。特許文献１には、スーパーコンティニーム光を用いて歯髄の断層画像を取得する歯髄撮像装置が記載されている。特許文献２には、スーパーコンティニューム光を用いて、生体の皮膚や皮下組織の断層画像を取得する光断層計測装置が記載されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１４−６３０４２号公報
【特許文献２】特開２０１５−１６７８２３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかし、非破壊で農作物の反射スペクトルや透過スペクトルを解析する方法では、農作物内部のより詳細な状態、たとえば農作物の内部構造や糖度分布、成熟度などを計測することはできない。また、農作物を切断して農作物の内部構造や糖度分布を計測するのは非常に手間や時間がかかり、商品となる農産物への適用は困難である。
【０００６】
また、特許文献１、２に記載のスーパーコンティニューム光によるＯＣＴ計測は、農作物の内部状態を観測することには従来利用されていなかった。
【０００７】
そこで本発明の目的は、農作物の内部の状態を非破壊でより詳細に計測することができる農作物観測方法およびその農作物観測装置を提供することである。また、その農作物観測方法を用いた農作物の生産方法を提供する。
【課題を解決するための手段】
【０００８】
本発明は、農作物の内部の状態を非破壊で観測する農作物観測方法であって、中心波長が０．７〜２．２μｍのスーパーコンティニューム光を用いたＯＣＴ計測を行うことにより、前記農作物の内部の状態を計測する、ことを特徴とする農作物観測方法である。
【０００９】
農作物の内部の状態は、たとえば、農作物内部の糖類の濃度分布や、農作物の内部の細胞壁の状態、農作物の葉の葉脈の状態である。特に、農作物内部のデンプンの濃度分布の計測が可能であり、米に対してこれを計測することで、登熟度の評価を非破壊、非接触で行うことが可能となる。デンプンの濃度分布の計測には、スーパーコンティニューム光の中心波長が１．３〜２．０μｍであることがより好ましい。糖類の濃度分布は、たとえば多変量解析などの統計的な推定方法を用いて推定することができる。
【００１０】
ＯＣＴ計測は、時間領域でもよいし周波数領域でもよいが、測定精度向上や測定時間短縮の点では周波数領域ＯＣＴが望ましい。また、周波数領域ＯＣＴにおいて、合わせて農作物の分光情報を取得することで、農作物の内部の状態をより詳細に観測することができる。
【００１１】
また本発明は、中心波長が０．７〜２．２μｍのスーパーコンティニューム光を生成する光源と、光源からのスーパーコンティニューム光を分岐して、農作物に照射して反射光させる測定光を得るための測定光系と、その反射光との干渉光を得るための参照光を得る参照光系と、測定光と参照光との干渉光を得て、農作物の内部の状態を反映した干渉画像を得る画像生成装置と、を有することを特徴とする農作物観測装置である。
【００１２】
本発明の農作物観測装置は、時間領域で干渉波形を取得する装置であってもよいし、周波数領域で干渉波形を取得する装置であってもよい。測定精度向上や測定時間短縮の点では周波数領域が望ましい。
【００１３】
参照光系と測定光系とを、手持ちが可能な筐体に一体に納めたハンディプローブを有する構成としてもよい。ハンディプローブとすることでその場計測が容易となる。また、手持ちとすることによる手振れによる影響が、参照光系と測定光系で等しく発生するため、キャンセルさせて手振れの影響を軽減することができる。特に、周波数領域の農作物観測装置において、このようなハンディプローブを構成するとよい。手振れの影響がより低減される。
【００１４】
また本発明は、中心波長が０．７〜２．２μｍのスーパーコンティニューム光を生成する光源であって、ＯＣＴ計測によって農作物の内部の状態を非破壊で観測するためのスーパーコンティニューム光源である。
【発明の効果】
【００１５】
本発明によれば、農作物の内部の状態（登熟度、細胞壁、葉脈など）を非破壊、非接触で計測することができる。そのため、農作物の状態を随時計測して、その計測結果を農作物の育成方法にフィードバックすることにより、農作物の品質を管理しながら農作物を育成することが可能となる。
【図面の簡単な説明】
【００１６】
【図１】実施例１の農作物観測装置の構成を示した図。
【図２】スーパーコンティニューム光源１００の構成を示した図。
【図３】スーパーコンティニューム光のスペクトルを示した図。
【図４】干渉光の波形を示した図。
【図５】変形例の農作物観測装置の構成を示した図。
【図６】変形例の農作物観測装置の構成を示した図。
【図７】コシヒカリ（玄米、正常米）の光学顕微鏡写真およびＯＣＴ画像を示した図。
【図８】コシヒカリ変異体（玄米、正常米）の光学顕微鏡写真およびＯＣＴ画像を示した図。
【図９】コシヒカリ変異体（玄米、腹白米）の光学顕微鏡写真およびＯＣＴ画像を示した図。
【図１０】コシヒカリ変異体（白米、未熟米）のＳＥＭ画像およびＯＣＴ画像を示した図。
【図１１】コシヒカリ変異体（白米、正常米）のＳＥＭ画像およびＯＣＴ画像を示した図。
【図１２】コシヒカリ（白米、正常米）のＳＥＭ画像およびＯＣＴ画像を示した図。
【図１３】イチゴの葉のＯＣＴ画像を示した図。
【図１４】トマトの葉のＯＣＴ画像を示した図。
【図１５】トマトの葉の表側のＯＣＴ画像を示した図。
【図１６】トマトの葉の表側のＯＣＴ画像を示した図。
【図１７】トマトの実のＯＣＴ画像を示した図。
【図１８】トマトの実のＯＣＴ画像を示した図。
【発明を実施するための形態】
【００１７】
以下、本発明の具体的な実施例について図を参照に説明するが、本発明は実施例に限定されるものではない。
【実施例１】
【００１８】
図１は、実施例１の農作物観測装置の構成を示した図である。実施例１の農作物観測装置は、主として、スーパーコンティニューム光源１００と、光ファイバーカプラ６１ａ、６１ｂ、６２と、ミラー６５と、駆動装置６６と、Ｘ−Ｙガルバノスキャナ７０と、受光装置であるフォトディテクタ（以下、「ＰＤ」と記す）７４ａ、７４ｂと、差動増幅器７５と、干渉波形解析手段であるコンピュータ７６と、表示装置７７と、によって構成されている。
【００１９】
（スーパーコンティニューム光源１００の構成）
まず、スーパーコンティニューム光源１００の構成について、図２を参照に詳しく説明する。
【００２０】
図２のように、図１に示すスーパーコンティニューム光源１００は、主として、ファイバーレーザ１０、ファイバー増幅器２０、波長シフトファイバー３０、スーパーコンティニューム光を生成する高非線形光ファイバー５０、シングルモードファイバー５１とから成る。
【００２１】
スーパーコンティニューム光源１００のファイバーレーザ１０は、エルビウムドープ光ファイバー（ＥＤＦ）１１を用いたリング型レーザ発振器である。エルビウムドープ光ファイバー（ＥＤＦ）１１の発振するレーザ光の中心波長は、１５５８ｎｍである。光ファイバー１１の長さは、３．５ｍである。光ファイバー１１の光路には、カーボンナノ構造体が塗布された樹脂フィルム１２が挿入されている。カーボンナノ構造体は、グラフェンが筒状になったカーボンナノチューブ、グラフェンが壁状に単層、多層に形成されたカーボンナノウォール、フラーレン、グラフェンなどである。本実施例では、カーボンナノチューブ粉末をポリイミドフィルムに塗布したものを用いた。励起光源１４には、波長９８０ｎｍのレーザダイオードを用いた。励起光源１４の出力する連続光は、波長分割多重カプラ（ＷＤＭ）１５を介して、光ファイバー１１に入力される。連続光は、アイソレータ１３により、光ファイバー１１の中を一方向にのみ循環する。
【００２２】
カーボンナノチューブは、１５６０ｎｍ帯の光を吸収し、光強度が大きくなると透明となる可飽和吸収特性を有している。発振するレーザ光の振幅は雑音成分により時間的に高周波数で変動しているために、瞬時的に透明となる閾値を越える瞬間の光だけがカーボンナノチューブで吸収されずに、通過することになる。このパルス光が光ファイバー１１を巡回する連続光に重畳されて、誘導放出が促進されてそのパルス光の強度が大きくなり、さらに、カーボンナノチューブを透過し易くなる。このようにしてパルス光が正帰還により成長し、光ファイバー１１を巡回する間に、安定したソリトンパルスとなる。
【００２３】
樹脂フィルム１２およびそれに接続する光ファイバーのコネクタ部分は、窒素や希ガスなどの不活性ガスの雰囲気とするか、あるいは真空とするのがよい。樹脂フィルム１２中のカーボンナノチューブが酸化などによって劣化し、その結果ファイバーレーザ１０からのソリトンパルスの出力が経時的に変動してしまうのを抑制することができる。
【００２４】
本実施例では、リング型レーザ発振器の全分散値の総和が負になるように設定されており、ソリトンパルスが発生する。なお、リング型レーザ発振器の全分散値の総和をゼロとした場合は、ストレッチパルスとなり、正とした場合には、散逸性ソリトンパルスとなる。ファイバーレーザ１０は、必ずしもソリトンパルスを発生させるものでなくともよいが、ストレッチパルスや散逸性ソリトンパルスと比較して、ソリトンパルスを用いた方が高い効率でスーパーコンティニューム光を得ることができる。
【００２５】
ファイバーレーザ１０の光路は、樹脂フィルム１２を除いて全て偏波保持型のものを用いるとよい。ソリトンパルスの偏光方向を時間的に安定させることができる。
【００２６】
このソリトンパルスは、光ファイバーカプラ１６により分岐され、パルス拡幅ファイバー２２に入力される。パルス拡幅ファイバー２２は、ソリトンパルスのパルス幅をやや拡大して、波形を生成するためのシングルモードファイバーである。パルス拡幅ファイバー２２から出力されるソリトンパルスは、偏波制御素子２７を通過して、Ｓ波、又は、Ｐ波の直線偏光となる。
【００２７】
パルス拡幅ファイバー２２の出力端はアイソレータ２８と波長分割多重カプラ（ＷＤＭ）２３とを介して、長さ２ｍのエルビウムドープ光ファイバー（ＥＤＦ）２１に接続されている。ファイバー増幅器２０は、波長分割多重カプラ（ＷＤＭ）２３、２４、光ファイバー（ＥＤＦ）２１、発振波長１５６０ｎｍのエルビウムを励起させるための波長１４８０ｎｍのレーザ光を発振するレーザダイオード（ＬＤ）２５、２６とから成る。レーザダイオード２５、２６からの励起光は、それぞれ、波長分割多重カプラ（ＷＤＭ）２３、２４を介して、光ファイバー（ＥＤＦ）２１に、相互に対向する方向から入射する。したがって、レーザダイオード（ＬＤ）２５からの励起光は、光ファイバー（ＥＤＦ）２１で生成されるソリトンパルスの進行方向と同一方向に伝搬し、レーザダイオード（ＬＤ）２６からの励起光は、そのパルス光の進行方向とは反対方向に伝搬する。
【００２８】
光ファイバー（ＥＤＦ）２１の出力端は、波長分割多重カプラ（ＷＤＭ）２４とアイソレータ３２とを介して、波長シフトファイバー３０に接続されている。波長シフトファイバー３０はシングルモードファイバーであり、ソリトンパルスの波長帯で異常分散の特性を持つ。ファイバー中で誘起される非線形効果によって、強度が高いほどより長波長側にソリトンパルスの中心波長がシフトする。また、ファイバー長が長い方が波長がより長波長側にシフトするが、パルス幅が広がるため、適切な長さに調節されている。アイソレータ３２の両端は光ファイバーを直結し、波長シフトファイバー３０とアイソレータ３２は光ファイバーで直結する構造となっている。また、励起光源１４から波長シフトファイバー３０の出力端までは、全光伝送路が光ファイバーで直結されており、レンズなどの空間伝搬部を排除した。これにより、出力されるスーパーコンティニューム光を安定させることができた。なお、スーパーコンティニューム光の中心波長をシフトさせる必要がないのであれば、波長シフトファイバー３０は省略することができる。
【００２９】
波長シフトファイバー３０の出力端は、レンズ３３、半波長板（λ／２）３４、λ／４板３５、ローパスフィルタ（ＬＰＦ）３６、レンズ３７を介して、ソリトンパルスをスーパーコンティニューム光に変換するための高非線形光ファバー５０に接続されている。
【００３０】
高非線形光ファイバー５０には、入力するソリトンパルスの中心波長を基準に±２００ｎｍの範囲内の波長において、分散が零である正常分散高非線形ファイバー（ND-HNLF）を用いた。正常分散高非線形ファイバー（ND-HNLF）としては、非線形性γが２３Ｗ-1ｋｍ-1、長さ５ｍのファイバーを用いた。この高非線形光ファイバー５０の出力端は、光源の使用性を向上させるために、通常のシングルモードファイバー５１に直結されている。シングルモードファイバー５１の出力端からスーパーコンティニューム光が出力される。
【００３１】
尚、光ファイバー１１、パルス拡幅ファイバー２２、エルビウムドープ光ファイバー（ＥＤＦ）２１、波長シフトファイバー３０、高非線形光ファイバー５０、シングルモードファイバー５１、波長分割多重カプラ（ＷＤＭ）１５、２３、２４などの全光伝送路に、偏波保持型のファイバーや素子を用いても良い。偏波保持型を用いることで、スーパーコンティニューム光の出力を安定化させることができる。また、全光路を偏波保持型とすれば、レンズ３３とローパスフィルタ３６の間に設けている半波長板（λ／２）３４、λ／４板３５を省略してより構成を簡素にすることができる。
【００３２】
特に、光ファイバー（ＥＤＦ）２１は、偏波保持型であって、ソリトンパルスの波長帯（すなわち中心波長をλ０、半値幅をΔλとして、λ０−Δλからλ０＋Δλの波長帯域）で正常分散特性を有しているものが望ましい。偏波面およびパルス波形を保持したまま十分な強度まで増幅することができる。
【００３３】
ファイバーレーザ１０の出力するソリトンパルスの中心波長における分散値β２は、０ｐｓ2 ／ｋｍより大きければ任意の値でよい。ただし、分散値β２は０ｐｓ2 ／ｋｍより大きく＋６５ｐｓ2 ／ｋｍ以下とすることが望ましい。分散値β２がこの範囲外であると、パルス波形を保持したまま十分な強度まで増幅することが難しくなるためである。さらに望ましくは＋５ｐｓ2 ／ｋｍ以上＋６０ｐｓ2 ／ｋｍ以下である。
【００３４】
また、実施例１では光ファイバー（ＥＤＦ）２１の励起光源としてレーザダイオード２５、２６の２つを用いているが、上記のような正常分散特性の光ファイバー（ＥＤＦ）２１を用いれば、エルビウムのドープ量を高濃度とすることができ、励起光源を１つ設けるのみで十分に増幅することができ、構成の簡素化、低コスト化を図ることができる。励起光の吸収の大きさは、４０ｄＢ／ｍ〜１４０ｄＢ／ｍとするのが望ましい。この範囲であれば、励起光源２５が１つでも十分な増幅が可能となる。より望ましい励起光の吸収の大きさは７０ｄＢ／ｍ〜１２０ｄＢ／ｍである。
【００３５】
スーパーコンティニューム光源１００から出力されるスーパーコンティニューム光（ＳＣ光）の特性は次の通りである。ＳＣ光のスペクトルは、図３に示すように、中心波長１７００ｎｍ、半値幅２４２ｎｍ、帯域１４００ｎｍ〜１９５０ｎｍである。また、パルス半値幅は、５ｐｓ、繰り返し周波数は１１０ＭＨｚ、最大出力強度は６０ｍＷである。
【００３６】
ＳＣ光のコヒーレンス長は、０．５〜１０μｍが望ましい。この範囲であると、干渉ピークの幅が広がったり、干渉ピーク同士が重なったりするなどの問題を防止することができ、また干渉ピークも良好に測定することができるため、測定精度を向上させることができる。より望ましいＳＣ光のコヒーレンス長は、０．５〜２μｍである。
【００３７】
ＳＣ光の波長は、上記に限るものではないが、農作物７２の内部の状態を高精度に観測のためには０．７〜２．２μｍとすることが望ましい。農作物７２の内部の観測においては、農作物７２における糖類（デンプン、セルロース、果糖、ショ糖、などの単糖類や多糖類）の分布を計測することが重要となる。これら糖類は近赤外線領域（０．８〜２．２μｍ）に吸収がある。そこで、ＳＣ光の中心波長は、０．７〜２．２μｍとすることが望ましい。この範囲とすることで、農作物における糖類の分布を高精度に計測することができる。特に、農作物７２内部のデンプンの濃度分布を計測するためには、１．３〜２．０μｍとするのがよい。また、分解能向上のためにＳＣ光のスペクトルの帯域幅は広いことが望ましく、１００ｎｍ以上が望ましい。より望ましくは２００〜２０００ｎｍである。
【００３８】
次に、実施例１の農作物観測装置におけるスーパーコンティニューム光源１００以外の各構成について、順に説明する。
【００３９】
光ファイバーカプラ６１ａ、６１ｂは、２×２ポートの構造であり、光ファイバーカプラ６１ａの入力ポートの一方にはスーパーコンティニューム光源１００からのＳＣ光が入力され、光ファイバーカプラ６１ａの出力端は光ファイバーカプラ６２に接続されている。光ファイバーカプラ６２に入力した光は、測定光と参照光に分配されて２つの出力ポートから出力される。２つの出力ポートには、偏波面の方位を一定方向に安定させる偏波コントローラ６３ａ、６３ｂが接続されており、それらの出力には、それぞれコリメータ６４ａ、６４ｂが接続されている。
【００４０】
コリメータ６４ａを透過した参照光は、コーナーキューブプリズム６７で曲げられてミラー６５に入射する。また、コーナーキューブプリズム６７は駆動装置６６により得るべき画像の深さ以上の範囲で移動される。ミラー６５で反射された参照光（以下、「反射参照光」という）は、同一光路を通り、光ファイバーカプラ６２に入力する。
【００４１】
また、コリメータ６４ｂを透過した測定光は、分散補償ガラス６８を透過して、Ｘ−Ｙガルバノスキャナ７０により反射されて、レンズ７１を介して、農作物７２に照射される。すなわち、Ｘ−Ｙガルバノスキャナ７０の変動により、農作物７２への測定光の照射方向が農作物７２の横断面であるｘ−ｙ平面上で走査される。農作物７２で反射された測定光（以下、「反射測定光」という）は、同一光路を通り、光ファイバーカプラ６２に入力する。分散補償ガラス６８は、ミラー６５とコーナーキューブプリズム６７とに受ける参照光の分散と、測定光の分散とを同一に補償するための分散補償ガラスである。
【００４２】
レンズ７１の焦点距離は、１．５ｃｍ以上とするのがよい。深さ方向の位置ずれに対して農作物７２に照射される測定光のスポット径の変動が小さくなり、測定精度をより向上させることができる。ただし、焦点距離を長くしすぎると、十分に集光した測定光を農作物７２に照射することが難しくなるため、５０ｃｍ以下とするのがよい。
【００４３】
光ファイバーカプラ６２において、反射参照光と反射測定光とが合波されて干渉光が生成される。駆動装置６６によるコーナーキューブプリズム６７の移動により、光ファイバーカプラ６２では、反射測定光と反射参照光とのコヒーレント長における相互相関である干渉光が得られる。この干渉光の波形が図４に示すものとなり、包絡線の半値幅は６．２μｍとなる。この値は、ＳＣ光のコヒーレンス長に対応している。
【００４４】
光ファイバーカプラ６２からは干渉光が２分配されて、それぞれの分配された干渉光が光ファイバーカプラ６１ａ、６１ｂに入力する。分配された干渉光は、位相がπだけ異なる。光ファイバーカプラ６１ａ、６１ｂから出力されるそれぞれの干渉光は、コリメータ７３ａ、７３ｂを透過し、平行光線としてＰＤ７４ａ、７４ｂに入射する。ＰＤ７４ａ、７４ｂの出力は、多数のＳＣパルスについては積分値であり、コーナーキューブプリズムの移動に応じて変動する干渉光の包絡線である。その出力は差動増幅器７５に入力している。差動増幅器７５により同相成分が除去され、干渉成分のみを取り出すことができる。
【００４５】
差動増幅器７５の出力は所定時間間隔でサンプリングされ、Ａ／Ｄ変換されて、コンピュータ７６に入力されて画像処理され、農作物７２の３次元のグレー画像が得られる。すなわち、農作物７２のある照射位置（ｘ０，ｙ０）（ｘ軸、ｙ軸は農作物７２への測定光の照射方向に対して垂直な方向）に対して、図４のように干渉光の波形が得られる。この干渉波形の振幅は、測定光と参照光との光路長が等しくなる農作物７２の深さ方向（ｚ軸方向）の位置での反射光強度を示している。したがって、駆動装置６６によって参照光の光路長を変化させることで、深さ方向ｚ（時間軸）に関する干渉光の反射光強度Ｄ（ｘ０，ｙ０，ｚ）、つまり１次元のグレー画像が得られる。そして、Ｘ−Ｙガルバノスキャナ７０の変動により、農作物７２への測定光の照射位置（ｘ，ｙ）を走査することで、３次元の反射光強度Ｄ（ｘ，ｙ，ｚ）が得られ、これを画像処理することで農作物７２の３次元のグレー画像を得ることができる。たとえば、任意方向の縦断面画像（ｚ軸方向に平行な断面での画像）、任意の深さｚでの横断面画像（ｘｙ平面に平行な断面での画像）、任意の角度から見た立体画像などを得ることができる。これらの画像は、表示装置７７において、任意に表示される。なお、画像処理において積算処理を行い、より高精度の画像が得られるようにしてもよい。
【００４６】
なお、実施例１の農作物観測装置は、３次元のグレー画像を得るものであるが、照射位置（ｘ，ｙ）をある１方向の走査のみとして、縦断面画像のみを取得するものとしてもよいし、Ｘ−Ｙガルバノスキャナ７０を省略して照射位置（ｘ，ｙ）は固定し、一箇所または複数箇所の１次元のグレー画像を取得するものとしてもよい。
【００４７】
ミラー６５は、十分に薄いものか、反射面が全波長帯域で高反射率のものを用いることで、被測定物による波長分散を補償し、干渉ピークが広がってピーク位置の取得精度が劣化するのを防止している。駆動装置６６は、たとえば、ボイスコイルモータ型ディレイライン、ピエゾチューブ型ディレイライン、直動ステージ型ディレイライン、積層ピエゾ型ディレイライン、などである。
【００４８】
ＰＤ７４ａ、７４ｂの受光装置として、フォトダイオード、フォトトランジスタ、光電子増倍管などを用いることができる。なお、測定光系は、スーパーコンティニューム光源１００の出力端から、光ファイバーカプラ６１ａ、６２、偏波コントローラ６３ｂ、コリメータ６４ｂ、分散補償ガラス６８、Ｘ−Ｙガルバノスキャナ７０、レンズ７１と、それらを接続するファイバーと、農作物７２からの反射光を同一光路で光ファイバーカプラ６２に導く光路、光ファイバーカプラ６１ａ、６１ｂ、コリメータ７３ａ、７３ｂ、ＰＤ７４ａ、７４ｂ、差動増幅器７５と、コンピュータ７７、表示装置７７とから成る。また、参照光系は、スーパーコンティニューム光源１００の出力端から、光ファイバーカプラ６１ａ、６２、偏波コントローラ６３ａ、コリメータ６４ａ、コーナーキューブプリズム６７、ミラー６５、それらを接続するファイバーと、ミラー６５からの反射光を同一光路で光ファイバーカプラ６２に導く光路、光ファイバーカプラ６１ａ、６１ｂ、コリメータ７３ａ、７３ｂ、ＰＤ７４ａ、７４ｂ、差動増幅器７５と、コンピュータ７６、表示装置７７とから成る。
【００４９】
以上述べた実施例１の農作物観測装置は、時間領域での光断層測定（ＴＤ−ＯＣＴ）を行う装置であるが、波長領域での光断層測定（ＯＦＤＲ−ＯＣＴ）を行う装置であってもよい。ＯＦＤＲ−ＯＣＴは、参照光を反射させるミラー６５を機械的に動作させる必要がなく、測定精度向上や測定時間短縮の点で利点がある。
【００５０】
図５は、そのような変形例の農作物観測装置の構成を示した図である。図５のように、変形例の農作物観測装置は、主として、スーパーコンティニューム光源１００と、光ファイバーカプラ２６１と、ミラー６５と、Ｘ−Ｙガルバノスキャナ７０と、回折格子２７０と、ＣＣＤカメラ２７１と、干渉波形解析手段であるコンピュータ７６と、表示装置７７と、によって構成されている。なお、同一の符号のものは実施例１の農作物観測装置と同様の構成部分であり、説明を省略する。
【００５１】
光ファイバーカプラ２６１は、２×２ポートの構造であり、その入力ポートの一方にはスーパーコンティニューム光源１００からのＳＣ光が入力される。光ファイバーカプラ２６１に入力した光は、測定光と参照光に分波されて２つの出力ポートから出力される。２つの出力ポートには、偏波面の方位を一定方向に安定させる偏波コントローラ６３ａ、６３ｂが接続されており、それらの出力には、それぞれコリメータ６４ａ、６４ｂが接続されている。
【００５２】
コリメータ６４ａを透過した参照光は、フィルタ２６２、可変ＮＤフィルタ２６３を透過した後、ミラー６５に入射する。フィルタ２６２は、ミラー６５による分散を補償して、参照光の分散と測定光の分散を同一とするために設けたものである。また、可変ＮＤフィルタ２６３は、透過率を可変とすることで透過する参照光の光量を調整するものであり、検出器の飽和を防ぎながら干渉信号の大きさを最適にするために設けたものである。ミラー６５は、最も強い干渉が生じるような位置に固定されている。ミラー６５で反射された参照光（反射参照光）は、同一光路を通り、光ファイバーカプラ２６１に入力する。
【００５３】
また、コリメータ６４ｂを透過した測定光は、ＮＤフィルタ２６４、Ｘ−Ｙガルバノスキャナ７０により反射されて、レンズ７１を介して、農作物７２に照射される。農作物７２で反射された測定光（反射測定光）は、同一光路を通り、光ファイバーカプラ２６１に入力する。
【００５４】
反射参照光と反射測定光は、光ファイバーカプラ２６１において合波されて干渉光が生成され、スーパーコンティニューム光源１００が接続されているポートとは別のポートから出力される。そして、偏波コントローラ２７３によって偏波面を安定させた後、コリメータ２７４を透過する。
【００５５】
コリメータ２７４を透過した干渉光は、回折格子２７０に入射し、波長ごとに異なる方向に反射される。これをレンズ２７５によって波長ごとに平行光に変換した後、レンズ２７６によって波長ごとに集光する。その集光された干渉光はＣＣＤカメラ２７１に入力する。
【００５６】
ＣＣＤカメラ２７１は、リニアイメージセンサを有し、干渉光の強度に応じた受光信号を出力する。そして、受光信号をＰＣ７６にて画像処理する。干渉光は波長ごとにＣＣＤカメラ２７１の異なる位置に入射するため、波長ごとの干渉光の強度分布、つまり干渉光の波長スペクトルが得られ、これを逆フーリエ変換すれば、農作物７２のある照射位置（ｘ０，ｙ０）に対する深さ方向ｚ（時間軸）に関する干渉光の反射光強度Ｄ（ｘ０，ｙ０，ｚ）、つまり１次元のグレー画像が得られる。そして、実施例１と同様に、Ｘ−Ｙガルバノスキャナ７０によって農作物７２への測定光の照射位置（ｘ，ｙ）を走査することで、３次元の反射光強度Ｄ（ｘ，ｙ，ｚ）が得られ、これを画像処理することで農作物７２の３次元のグレー画像を得ることができる。このようにして得られた画像は、表示装置７７において表示される。この画像処理は、３次元の反射光強度Ｄ（ｘ，ｙ，ｚ）を複数回取得して積算処理を行うものであってもよい。より高精度な画像を得ることができる。
【００５７】
コリメータ２７４と回折格子２７０の間の光路上にＡＯＭ（音響光学変調器）を設け、干渉光の回折格子２７０への入射角度を制御できるようにしてもよい。ＣＣＤカメラ２７１による波長スペクトル測定において、１画素当たりの波長帯域幅が変わるため、測定波長範囲や波長分解能を制御でき、測定精度の向上を図ることができる。もちろん、回折格子２７０の姿勢を機械的に制御することで、干渉光の回折格子２７０への入射角度を制御してもよい。
【００５８】
画像の取得に加えて、農作物７２の分光情報（反射スペクトル）、つまり測定光の波長スペクトルを取得してもよい。これは、ＣＣＤカメラ２７１によって得られる受光信号をＰＣ７６にて解析することで抽出することができる。この分光情報と農作物７２の３次元のグレー画像を組み合わせることで、農作物７２内部の組成分布や、特定成分の濃度分布などを抽出することができる。
【００５９】
実施例１あるいはその変形例の農作物観測装置において、測定光系と参照光系とを一の筐体に納めて一体化し、手持ちが可能なハンディプローブとして構成してもよい。ハンディプローブとすることで、農作物７２への測定光の照射位置を使用者が自在に変更でき、農作物７２の断層画像の取得がより容易となる。たとえば、田畑において収穫前の農作物７２をその場でリアルタイム計測することが容易となる。また、測定光系と参照光系を一体としているため、測定光系に手振れなどが生じた場合に、参照光系も同様に手振れなどが生じ、測定光、参照光ともに手振れの影響を受ける。そのため、手振れによる測定誤差をキャンセルして軽減することができる。特に、ＯＦＤＲ−ＯＣＴである変形例の農作物観測装置において、このようなハンディプローブが有効である。変形例の農作物観測装置では、実施例１の農作物観測装置とは異なり参照光系のミラー６５を機械的に走査する必要がないため、ハンディプローブとしても測定誤差が少なくて済む。
【００６０】
図６は、変形例の農作物観測装置において、コリメータ６４ａを含む後段の参照光系と、コリメータ６４ｂを含む後段の測定光系を、手による持ち運びが可能な筐体に納めて一体化し、ハンディプローブ３００とした構成を示した図である。なお、コリメータ６４ａ、ｂよりも前段部分をハンディプローブ３００に納めてもよい。ただし、コリメータ６４ａ、ｂより前段では光は光ファイバー内を伝搬するため手振れの影響が少なく、コリメータ６４ａ、ｂより後段は光が空間を伝搬し、手振れの影響を大きく受けるため、ハンディプローブ３００にはコリメータ６４ａ、ｂよりも後段を納めれば十分である。
【００６１】
Ｘ−Ｙガルバノスキャナ７０に測定光を照射する際、そのミラーの中心から少しずらした位置に照射するとよい。Ｘ−Ｙガルバノスキャナ７０のミラー中心に測定光を照射する場合、干渉光強度の深さ方向（ｚ方向）依存性において、正側か負側の一方のみしか使用することができなかった。つまり、基準位置に対して参照光の光路長が増加する方向の走査と減少する方向の走査を区別できなかった。しかし、上記のように中心からずらした位置に照射することで、正・負両側を使用することができるため、深さ方向の測定レンジを２倍にすることができる。その結果、測定精度の向上を図ることができる。
【００６２】
また、２つのＳＣ光源を用意するなどして、ＳＣ光の中心波長を変えてそれぞれで農作物７２の３次元のグレー画像を取得するようにしてもよい。両者の画像を比較して差異点を抽出することで、より詳細な内部構造を解析することができる。
【００６３】
また、ＯＣＴ計測の際、農作物７２の表面にジェルを塗布してもよい。ジェルによって農作物７２の表面での反射を抑えることで、画像のノイズを低減することができ、より高精度な画像を取得することができる。
【００６４】
以上のようにして取得した農作物７２の３次元のグレー画像から、種々の情報を抽出することができる。
【００６５】
第１に、農作物７２の内部であって表面近傍の状態を、非破壊、非接触で高精度に観測することができる。光源として中心波長が０．７〜２．２μｍのＳＣ光を用いており、農作物７２内部へとある程度侵入させることができるためである。また、ＳＣ光であるため帯域幅が広く、高分解能であるため、内部構造も高精度に観測することができる。特に、細胞壁や葉脈などの内部構造の観測に適している。果実の成熟度が高いほど細胞が大きくなる傾向があるため、細胞の大きさを観測することで果実の成熟度の評価が可能となる。
【００６６】
第２に、農作物７２内部の糖類の濃度分布を抽出することができる。デンプン、セルロース、果糖、ショ糖などの糖類は、近赤外線領域に吸収を有している。そのため、取得した農作物７２の画像において、その糖類の吸収によって画像の濃淡が生じる。その濃淡は、糖類の濃度、種類、結晶化度などに依存している。したがって、農作物７２の３次元のグレー画像から、糖類の濃度分布を抽出することが可能である。抽出する方法としては、統計的な推定方法を用いることができ、たとえば、多変量解析などを用いることができる。特に、デンプンの濃度分布を高精度に求めることができる。農作物７２内部のデンプンの濃度分布が分かれば、農作物７２の成熟度などを非破壊、非接触で判断することができるようになる。成熟度が低いと、デンプンの結晶の角が取れて丸くなる傾向があり、デンプンの結晶間に隙間が生じる。その隙間によってクラックやデンプンの濃度分布が生じるため、成熟度の評価が可能となる。
【００６７】
農作物７２の内部構造や糖類の濃度分布は、一の断面画像から評価してもよいが、複数の断面画像から評価することでより精度を向上させることができる。
【００６８】
また、上記の農作物７２の内部構造や糖類の濃度分布の測定は、非破壊、非接触で行うことができるので、農作物７２の収穫後のみならず、農作物７２の育成中に行うこともでき、農作物７２の品質評価に利用することができる。たとえば、収穫後あるいは育成中の米、麦、豆などの穀物の抽出評価に利用したり、トマト、ブドウ、モモなどの野菜類、果物類のフィールドでの検査に利用したりすることができる。特に、米の登熟度（デンプンの充填度、成熟度）の評価に好適である。
【００６９】
特に、食用米、酒米、ワイン用ブドウなどの農作物７２に好適である。これらの農作物７２は高い品質を求められるが、その成熟度の指標として糖度や充填率を評価する為に、従来は破壊的に測定するか、また、光学的な機器を用い非破壊的に評価する場合は、プローブ等を接触させる必要があること、その果皮組織の発達程度によっては誤差が生じることが問題であった。しかし、本発明の農作物観測方法を用いれば、育成中の段階で非破壊、非接触で農作物の内部の状態を短時間で観測することができるので、品質評価へのフィードバックなども容易となり、たとえば育成中において品質管理が可能となる。
【００７０】
次に、実施例１の農作物観測装置を用いて農作物を観測した結果について説明する。
【００７１】
図７〜９は、米について観測した結果を示した図である。図７は、コシヒカリ（玄米、正常米）、図８は、コシヒカリ変異体（玄米、正常米）、図９は、コシヒカリ変異体（玄米、腹白米）である。コシヒカリ変異体は登熟度が低くなる変異体であり、コシヒカリに比べ、米のデンプン粒の発達程度が低い。また、一般的な外観品質基準に従い正常米と腹白米とを分類している。図７〜９において（ａ）は、米の外観を示した光学顕微鏡写真であり、（ｂ）は、ＯＣＴ画像（深さ方向の断面画像であり、（ａ）におけるｘｚ面での断面画像）である。ＯＣＴ画像取得に用いたＳＣ光は、中心波長１．７μｍ、半値幅２４２ｎｍ、分解能４．３μｍ、ビーム径は２０μｍである。
【００７２】
図７〜９のＯＣＴ画像を比較すると、コシヒカリとその変異体の識別、正常米と腹白米の識別が可能であることがわかる。図７〜９では、玄米のぬか層の部分が明瞭な黒色で示されており、コシヒカリ（玄米、正常米）では、図７のように、ぬか層直下の胚乳部分がほぼ均一な薄いグレーに表示されている。これに対し、コシヒカリ変異体（玄米、正常米）では、図８のように、ぬか層直下の胚乳部分は、濃い目のグレーで表示されている。また、コシヒカリ変異体（玄米、腹白米）では、図９のように、ぬか層直下の胚乳部分はやや薄いグレーであるが、ところどころ黒色の斑点が表示されている。このグレーの濃さと黒色の斑点は、米が未成熟のためデンプンの結晶化が不十分であり、それによりデンプン粒間にクラックや隙間が生じるためであると推察される。このＯＣＴ画像の濃淡は、デンプンの濃度や結晶化度を反映していると考えられる。このように、ＯＣＴ画像の濃淡によって、米の内部の状態を非破壊、非接触で観測することができ、コシヒカリとその変異体の識別、正常米と腹白米の識別が可能となっている。
【００７３】
図１０〜１２は、白米について観測した結果を示した図である。図１０は、コシヒカリ変異体（未熟米）、図１１は、コシヒカリ変異体（正常米）、図１２は、コシヒカリ（正常米）である。図１０〜１２において、（ａ）は白米を中央部で切断した断面のＳＥＭ画像であり、（ｂ）は、ＯＣＴ画像（深さ方向に垂直な断面画像）である。用いたＳＣ光は図７〜９の場合と同様である。
【００７４】
図１０〜１２のＳＥＭ画像を見ると、デンプンの結晶化度の違いや、細かな線状のクラックの有無などに違いがあることが認められ、特に図１２のコシヒカリ（正常米）では、非常に均一でなめらかな断面で、クラックが少ないことがわかる。
【００７５】
一方、図１０〜１２のＯＣＴ画像を見ると、胚乳部分の色の濃淡や、胚乳中の黒い斑点の数や大きさなどに違いがあることが認められる。特に、図１２のコシヒカリ（正常米）では、黒い斑点が少なく、胚乳部分は薄いグレーでほぼ均一である。
【００７６】
このように、図７〜１２のＯＣＴ画像によると、従来は米を切断してＳＥＭ画像を取得するなどしなければわからなかった米の内部の状態を、非破壊、非接触で観測できることがわかった。
【００７７】
図１３は、イチゴの葉のＯＣＴ画像（深さ方向に垂直な断面画像で深さ１５０μｍ、１ｍｍ角）である。用いたＳＣ光は図７〜９の場合と同様である。図１３のように、ツリー状の濃いグレーの部分は、葉脈を示している。このように、実施例１の農作物観測装置を用いたＯＣＴ計測により、葉の葉脈を非破壊、非接触で観測できることがわかる。このように葉脈を観測できるのは、葉脈を構成する細胞の配置が周辺と異なるため、葉脈の構造とその周辺構造との各境界から反射、散乱が生じており、それがＯＣＴ画像に濃淡として現れるためである。ＯＣＴ画像の取得に加えて分光情報を抽出することで、葉脈を構成する各成分の吸収の影響を観測することが可能となる。
【００７８】
図１４は、トマトの葉のＯＣＴ画像（深さ方向の断面画像）である。用いたＳＣ光は、図７〜９と同様である。ＯＣＴ画像は、葉の表面側の断面、葉の裏側の断面、葉の裏側に反射防止のジェルを塗布したときの葉の裏側の断面、葉の近くの茎の断面、葉の近くの茎の別の箇所の断面を示している。図１４のように、葉の場所によって異なる内部構造となっていることがわかる。また、葉の裏側に反射防止のジェルを塗布することで、ＯＣＴ画像のノイズが低減されていることがわかる。
【００７９】
図１５、１６は、トマトの葉の表側のＯＣＴ画像であり、（ａ）は深さ方向の断面画像、（ｂ）は深さ方向に垂直な断面画像である。また、図１５は、ＳＣ光として中心波長０．８μｍ、半値幅１２０ｎｍを用いており、図１６は、ＳＣ光として、図７〜９と同様に、中心波長１．７μｍ、半値幅２４２ｎｍを用いている。
【００８０】
図１５、１６のように、中心波長０．８μｍと中心波長１．７μｍとでは、ＯＣＴ画像に差異があることがわかる。ＳＣ光の中心波長を短くすれば、内部への進入度合は低下するものの、より高解像度なＯＣＴ画像を取得できる。また、ＳＣ光の中心波長を短くすることで、ＯＣＴ画像は、短波長側の吸収スペクトル、反射スペクトルを反映したものとなる。したがって、ＳＣ光の中心波長を変えて複数のＯＣＴ画像を取得し、それらのＯＣＴ画像を比較することで、より詳細な内部構造を解析することができる。
【００８１】
図１７、１８は、トマトの実のＯＣＴ画像（深さ方向の断面画像）を示した図である。用いたＳＣ光は、図１５と同様に、中心波長０．８μｍ、半値幅１２０ｎｍである。
【００８２】
図１７のＯＣＴ画像は、それぞれ以下の部分を計測したものである。１はトマトの実の赤い部分、２は同じく赤い部分であって表面に反射防止のジェルを塗布した場合、３は実のお尻部分、４は実の赤い部分であって皮を剥いた状態、５は実の赤い部分の皮を部分的に剥ぎ、皮を剥いだ部分と皮のある部分との境目部分、６は赤い皮の裏側、７は実の側面付近であって、ヘタから頂点に向かう方向での断面、を示している。
【００８３】
また、図１８のＯＣＴ画像は、それぞれ以下の部分を計測したものである。８は実の赤い部分であってヘタ付近、９は実の切り口の硬い部分、１０は実の切り口であって９よりは柔らかい部分、１１は実の切り口であって１０よりもさらに柔らかい部分、１２はトマトの切り口の種部分、１３はトマトの切り口であってヘタ付近の緑色部分、を示している。
【００８４】
図１７、１８のように、計測箇所の違いによって、トマトの実の内部の構造にも違いが生じていることがわかる。また、図１７、１８では、網目状に濃いグレーの部分が示されており、これはトマトの実の細胞壁を示している。この細胞の大きさは、トマトの実の成熟度を反映しており、成熟度が高いほど細胞が大きくなる傾向にある。そのため、細胞の大きさを観測して評価することで、トマトの成熟度を評価することができ、トマトの実の糖度などを推定することができる。
【実施例２】
【００８５】
実施例２の農作物の生産方法では、まず、実施例１の農作物観測装置を用いた農作物観測方法により、農作物を収穫することなく、育成途中の段階において、農作物の内部状態を観測する。たとえば、農作物の成熟度や糖度分布、細胞の大きさや形状、葉の葉脈の太さや形状などである。これは全量検査、全数検査でもよいし、抜き取り検査としてもよい。そして、その観測結果を育成方法にフィードバックする。たとえば、光量の調整、追肥の有無や量、水やりの量、温度、湿度、土の状態などを調整する。育成方法の調整は、株単位や農作物単位で行うこともできる。農作物の内部状態の観測と、育成方法へのフィードバックを繰り返すことで、育成中の農作物の品質を適正に評価、管理することができる。また、農作物の内部状態を観測することで、適正な収穫時期を判断することができる。その結果、農作物の品質のばらつきを低減したり、農作物の品質向上を図ったりすることができ、より高品質で高付加価値の農産物を生産することができる。
【産業上の利用可能性】
【００８６】
本発明は、農作物の表面近傍の断面の詳細な状態を非破壊、非接触で観測することができるので、農作物の成熟度などを非破壊、非接触で検査することができる。また、収穫後の農作物のみならず、育成中の農作物についても観測することができる。そのため、農作物を育成しながら、随時、農作物の成熟度などを判定して、それを育成方針にフィードバックすることにより、農作物の品質のばらつきを低減したり、農作物の品質向上を図ったりすることができる。
【符号の説明】
【００８７】
１００：スーパーコンティニューム光源
１０：ファイバーレーザ
１１：エルビウムドープ光ファイバー
１２：フィルム
２０：ファイバー増幅器
２１：エルビウムドープ光ファイバー
２２：パルス拡幅ファイバー
３０：波長シフトファイバー
１５、２３、２４：波長分割多重カプラ（ＷＤＭ）
１３、２８、３２：アイソレータ
５０：高非線形光ファイバー
１６、６１ａ、６１ｂ、６２、２６１：光ファイバーカプラ
２７０：回折格子
２７１：ＣＣＤカメラ
３００：ハンディプローブ
【図１】
2017181300.tif 000003
【図２】
2017181300.tif 000004
【図３】
2017181300.tif 000005
【図４】
2017181300.tif 000006
【図５】
2017181300.tif 000007
【図６】
2017181300.tif 000008
【図７】
2017181300.tif 000009
【図８】
2017181300.tif 000010
【図９】
2017181300.tif 000011
【図１０】
2017181300.tif 000012
【図１１】
2017181300.tif 000013
【図１２】
2017181300.tif 000014
【図１３】
2017181300.tif 000015
【図１４】
2017181300.tif 000016
【図１５】
2017181300.tif 000017
【図１６】
2017181300.tif 000018
【図１７】
2017181300.tif 000019
【図１８】
2017181300.tif 000020