(57)【要約】  本発明は、以下の工程を含む、対象である土地に樹木を植えた際の該樹木の成長度合いを予測する方法を与える：  （ａ）上記土地の近傍に栽植された複数の樹木についての成長に関するデータを取得する工程  （ｂ）前記樹木のそれぞれの下の土壌について、土壌の色、土壌の特性及び／又は土壌の養分を分析する工程  （ｃ）該分析の結果から説明変数を得る工程  （ｄ）該説明変数を用いて、前記樹木についての成長に関するデータとの検量式を得る工程  （ｅ）上記土地の土壌について土壌の色、土壌の特性及び／又は土壌の養分を分析し、当該分析結果を前記検量式に当てはめて、上記土地に樹木を植えた際の該樹木の成長度合いを計算し、予測する工程。
2016043007.tif 000002
【特許請求の範囲】
【請求項１】
以下の工程を含む、対象である土地に樹木を植えた際の該樹木の成長度合いを予測する方法：
（ａ）上記土地の近傍に栽植された複数の樹木についての成長に関するデータを取得する工程
（ｂ）前記樹木のそれぞれの下の土壌について、土壌の色、土壌の特性及び／又は土壌の養分を分析する工程
（ｃ）該分析の結果から説明変数を得る工程
（ｄ）該説明変数を用いて、前記樹木についての成長に関するデータとの検量式を得る工程
（ｅ）上記土地の土壌について土壌の色、土壌の特性及び／又は土壌の養分を分析し、当該分析結果を前記検量式に当てはめて、上記土地に樹木を植えた際の該樹木の成長度合いを計算し、予測する工程。
【請求項２】
土壌の色を分析することを含み、分析の結果得られる説明変数として少なくともＬ*、ａ*、ｂ*を含む請求項１に記載の方法。
【請求項３】
土壌の特性を赤外分光法により分析することを含み、分析の結果得られる説明変数として少なくとも赤外分光スペクトルを含む請求項１に記載の方法。
【請求項４】
土壌の養分を分析することを含み、分析の結果得られる説明変数として少なくとも全炭素、全窒素、交換性Ｃａ、交換性Ｋ及び交換性Ｍｇの量を含む請求項１に記載の方法。
【請求項５】
土壌の特性及び土壌の養分を分析することを含む請求項１に記載の方法。
【請求項６】
土壌の色を分析することをさらに含む請求項５に記載の方法。
【請求項７】
請求項１〜６のいずれかに記載の方法により、植林予定地における樹木の成長度合いを植林のための樹木の植栽前及び／又は植栽後に予測することを含む、植林方法。
【請求項８】
以下を含む、請求項１〜６のいずれかに記載の方法を行うためのシステム：
（ア）土壌の色及び／又は土壌の特性を測定する測定装置；及び
（イ）前記装置を搭載した移動機具。
【請求項９】
測定装置を用いて得られた測定値を樹高に換算して樹木の成長度合いを予測する機能を備えた機材をさらに具備する、請求項８に記載のシステム。
【請求項１０】
移動機具が地上移動型機具又は空中移動型機具である、請求項８又は９に記載のシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、樹木成長の簡易予測法に関する。
【背景技術】
【０００２】
ある土地への植林前に、植栽が予定される樹木について植栽後の成長予測を行うことは、植林の効率を確保するうえで必要である。しかしながら、これまでかかる成長の予測が行われることは実質的になかった。予測のための方法が存在しないことが、その一因である。林業は農業とは異なり、施肥が難しいため、土壌環境が重要であるにもかかわらず、適切な土壌の管理の方法さえ明らかではないといった背景もある。【０００３】
樹木の成長量自体を評価する方法（特許文献１）、樹木の健全度を評価する方法（特許文献２）や、樹木の木材成分を定量する方法（特許文献３）については、報告がなされている。
また土壌に着目し、土壌分析方法として光スペクトル（近赤外分光スペクトル）から土壌の特性及び成分等を予測する方法についても報告がなされている（特許文献４）。
さらに非特許文献１には林分成長量の予測のために、過去の成長量を推定するための方法について記載されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００８−７９５４９号公報
【特許文献２】特開２００７−１６６９６７号公報
【特許文献３】特開平１１−２５８１６１号公報
【特許文献４】特開２００６−３８５１１号公報
【非特許文献】
【０００５】
【非特許文献１】「林分成長量の推定および予測方法に関する研究」西沢正久、林業試験研究報告、第１２９号、p.１〜p.１６１
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、特許文献１〜４に開示されている技術的事項は、いずれも植林前に樹木の成長予測を可能とするものではない。特許文献４に記載の方法は土壌の特性及び成分等を予測する方法に留まり、同文献にも樹木の成長あるいはその予測については記載されていない。また非特許文献１は過去の成長量を推定するための方法について記載するものであって、樹木の成長予測の方法を与えるものではない。
すなわち、植林前に樹木の成長予測を可能とする技術はこれまで知られていない。
上記背景に鑑み、本発明は、植林前に樹木の成長を予測する方法を提供することを目的とした。
【課題を解決するための手段】
【０００７】
そこで本発明者らは樹木の成長予測に資する指標を探索したところ、これまで前記成長予測に用いられていなかった特定の指標により上記課題が解決される可能性があることを見出し、さらに鋭意研究を進めた結果本発明を完成するに至った。
すなわち、本発明は、少なくとも以下の発明に関する：
［１］
以下の工程を含む、対象である土地に樹木を植えた際の該樹木の成長度合いを予測する方法：
（ａ）上記土地の近傍に栽植された複数の樹木についての成長に関するデータを取得する工程
（ｂ）前記樹木のそれぞれの下の土壌について、土壌の色、土壌の特性及び／又は土壌の養分を分析する工程
（ｃ）該分析の結果から説明変数を得る工程
（ｄ）該説明変数を用いて、前記樹木についての成長に関するデータとの検量式を得る工程
（ｅ）上記土地の土壌について土壌の色、土壌の特性及び／又は土壌の養分を分析し、当該分析結果を前記検量式に当てはめて、上記土地に樹木を植えた際の該樹木の成長度合いを計算し、予測する工程。
［２］
土壌の色を分析することを含み、分析の結果得られる説明変数として少なくともＬ*、ａ*、ｂ*を含む上記［１］に記載の方法。
［３］
土壌の特性を赤外分光法により分析することを含み、分析の結果得られる説明変数として少なくとも赤外分光スペクトルを含む上記［１］に記載の方法。
［４］
土壌の養分を分析することを含み、分析の結果得られる説明変数として少なくとも全炭素、全窒素、交換性Ｃａ、交換性Ｋ及び交換性Ｍｇの量を含む上記［１］に記載の方法。
［５］
土壌の特性及び土壌の養分を分析することを含む上記［１］に記載の方法。
［６］
土壌の色を分析することをさらに含む上記［５］に記載の方法。
［７］
上記［１］〜上記［６］のいずれかに記載の方法により、植林予定地における樹木の成長度合いを植林のための樹木の植栽前及び／又は植栽後に予測することを含む、植林方法。
［８］
以下を含む、上記［１］〜［６］のいずれかに記載の方法を行うためのシステム：
（ア）土壌の色及び／又は土壌の特性を測定する測定装置
（イ）前記装置を搭載した移動機具。
［９］
測定装置を用いて得られた測定値を樹高に換算して樹木の成長度合いを予測する機能を備えた機材をさらに具備する、［８］に記載のシステム。
［１０］
移動機具が地上移動型機具又は空中移動型機具である、［８］又は［９］に記載のシステム。
【発明の効果】
【０００８】
本発明によれば、少なくとも以下の効果が奏される：
（ア）土壌色や赤外分光スペクトル、土壌養分から成長予測を行うことが可能となる。
（イ）土壌において不足する養分が明らかになるため、施肥などによる土壌改良の際の指針を提供することができる。
【図面の簡単な説明】
【０００９】
【図１】本発明の概略を示すスキームである。
【図２】土壌の色を用いる本発明の方法において得られた樹木成長の予測値と実測値との相関を示す図である。  図２ＡはｐＨ及びＥＣを考慮せずＬ*、ａ*及びｂ*のみを用いて一般化線形モデルにより１次式で作成した場合であり、図２Ｂは一般化線形モデルにより２次までの交互作用を含めて作成した予測式から得られた予測高さと実測高さとの相関を示す。図２ＣはさらにｐＨ及びＥＣの測定値を考慮した結果を示す。
【図３】土壌の特性を用いる本発明の方法において得られた樹木成長の予測値と実測値との相関を示す図である。
【図４】赤外スペクトルデータ（図４Ａ）及びそのスムージング処理後のデータ（図４Ｂ）を示す図である。
【図５】養分１法（図５Ａ）及び養分２法（図５Ｂ）による結果を示す図である。
【図６】ユーカリについての、土壌の特性（赤外分光スペクトル）を用いる本発明の方法による樹木成長の予測の結果を示す図である。  縦軸は地位指数（平均樹高を表す指標）の実測値を、横軸は地位指数の予測値を、それぞれ表す。
【発明を実施するための形態】
【００１０】
以下に本発明をより詳細に説明する。
なお、本明細書において「土壌を分析（する）」とは、対象の土壌の色、土壌の特性及び／又は土壌に含まれる養分（土壌の養分）についての情報を得ることを意味する。
また。本明細書において「（樹木の）成長」とは、例えば樹木の樹高により定量される成長を意味するが、他の定量項目（例えば枝下高、太さ（胸高直径）であってもよい。
【００１１】
上記のとおり、少なくとも本発明は以下の工程を含む、対象である土地に樹木を植えた際の該樹木の成長度合いを予測する方法に関する：
（ａ）上記土地の近傍に栽植された複数の樹木についての成長に関するデータを取得する工程
（ｂ）前記樹木のそれぞれの下の土壌について、土壌の色、土壌の特性及び／又は土壌の養分を分析する工程
（ｃ）該分析の結果から説明変数を得る工程
（ｄ）該説明変数を用いて、前記樹木についての成長に関するデータとの検量式を得る工程
（ｅ）上記土地の土壌について土壌の色、土壌の特性及び／又は土壌の養分を分析し、当該分析結果を前記検量式に当てはめて、上記土地に樹木を植えた際の該樹木の成長度合いを計算し、予測する工程。
本発明の例についての手順を概念的に図１に示した。すなわち、本発明は土壌の色、土壌の特性及び土壌の養分等についての分析を行い、当該分析結果から検量線モデル及び／又は検量式を作製し、該検量線モデル及び／又は検量式を用いて、樹木成長の推定（予測）を行うことに基礎を置くものである。
【００１２】
本発明の方法について、以下により詳細に説明する。
なお本発明の方法は、いずれの方法においても
（１）当該樹木の樹高など樹木の成長に関するデータを取得する工程、及び
（２）当該樹木下の土壌を採取し分析に供するか、又はその場で分析データを得る工程は共通して行われる。
【００１３】
土壌の色、土壌の特性及び／又は土壌の養分を用いる本発明の方法のそれぞれについて、以下に説明する。
Ａ．土壌の色を用いる方法
土壌の色の分析は比較的容易に行うことができるため、当該分析結果を用いる方法は簡便さの点において好ましい。
土壌の色を用いる場合、Ｌａｂの情報を入手し、高さおよび養分との相関を解析する。
ＣＩＥＬＡＢの３つの座標は、色の明度（Ｌ*＝０は黒、Ｌ*＝１００は白の拡散色で、白の反射色はさらに高い）、赤／マゼンタと緑の間の位置（ａ*、負の値は緑寄りで、正の値はマゼンタ寄り）、黄色と青の間の位置（ｂ*、負の値は青寄り、正の値は黄色寄り）に対応している。
土壌の色を用いる場合、本発明の方法は少なくとも以下の各ステップを含む：  
（ａ１）スキャナー、デジタルカメラや土色計などの色情報取得装置により土壌の画像を取得し、画像解析ソフトによりサンプルの平均的な色情報（ＲＧＢなど）を得る。
（ｂ１）上記色情報をＣＩＥ（Commission Internationale de l’Eclairage：国際照明委員会）により策定されたＣＩＥ  Ｌａｂ色空間情報に変換する（ＲＧＢからＣＩＥ  Ｌａｂなど各色情報間の変換式は公開されている。http://www.cs.rit.edu/~ncs/color/t_convert.html）。
（ｃ１）得られた色情報についての変数であるＬ*、ａ*、ｂ*を説明変数とし、樹木成長データを目的変数として、一般化線形モデルにより、検量線モデル及び／又は検量式を作成する。
（ｄ１）前記検量線モデル及び／又は検量式を用いて、樹木成長の推定（予測）を行う。
変数の個数は限定されないが、当該個数を増やすことによってより正確な予測が可能になる場合がある。例えば変数を３個又は４個以上用いる当該方法は好ましく、６個以上用いる方法はより好ましい。
用いられる変数は、典型的には上記Ｌ*、ａ*及びｂ*であり、この場合検量式は一次式である。またこれら３つの変数に加えてＬａ、Ｌｂ及びａｂをさらに用いる方法は好ましい。各変数の交互作用を考慮して説明パラメーターの個数が増えるため、予測の精度が向上するからである。
また、土壌の色の分析結果に加重して、ｐＨ及び／又はＥＣの測定データを考慮すると一層精度の高い予測結果が得られるため好ましい。この傾向はＬ*、ａ*及びｂ*のみに加えてｐＨ及びＥＣの測定データを考慮した場合に顕著である。Ｌ*、ａ*及びｂ*は全炭素などの有機物の量や交換性養分の量をより強く反映する一方、予測に寄与するｐＨ及びＥＣについては反映しづらいためであると考えられる。
上記スキャナー、デジタルカメラや土色計といった色情報取得装置は、土壌の色情報を与えるものであれば限定されない。色情報取得装置の形態としては、携帯可能な小型のものが挙げられる。当該携帯可能な小型色情報取得装置を移動機具に搭載して、地こしらえの際に当該植林地全体の土壌についての色情報を取得するようにしたものは好ましい。移動機具を用いることにより、より広範な領域についての色情報を正確かつ効率的に取得することができるばかりでなく、地上移動型機具には表面より下方に存する土壌についての色情報の取得も可能になるからである。
色情報取得装置に加えて、得られた色情報を変換・解析し、解析結果を樹高に換算して樹木の成長度合いの推定・予測を行う機能をさらに具備した移動機具を含むシステムは、本発明の方法の実施のためにより好ましく用いられる。解析結果を樹高に換算する際に用いる換算式を環境条件に応じて改変できるシステムは、一層より好ましい。環境条件には、土壌型、気候、季節、地形及び微地形といった立地条件が、天候などの条件とともに包含される。
本発明において用いられる色情報取得装置を搭載した移動機具として、トラクターや林業機械（ハーベスタ及びスキッダ等）といった地上移動型機具、ならびに有人航空機や無人航空機（ドローン及びラジコンヘリ等）といった空中移動型機具が例示される。移動機具として、移動パターンが規則的である、トラクター及びハーベスタは好ましく、無人航空機も操作が簡便であるため好ましい。
【００１４】
Ｂ．土壌の特性を用いる方法
土壌の養分が樹木の成長に影響すると考えられるが、樹木の成長に必要な養分は多種類あり、成長との関連を明らかにするための方法は未だ確立されていないばかりでなく、全ての養分を測定することは、時間的費用的に制約がある場合が多い。また、樹木の成長には様々な養分が関連しあい、必要な養分量とともに、過剰である場合には却って成長制限要因となる場合もあるため、養分による樹木の成長予測は、必ずしも容易ではないといった事情がある。
一方、土壌の特性（全炭素量、ＥＣ（電気伝導度）、交換性Ｃａ量、交換性Ｍｇ量及び交換性Ｋ量等）を反映する土壌の赤外スペクトル情報を得て樹木成長を推定するための検量線モデルを作成すると、精度が高い樹木成長の推定（予測）を行うことができる。したがって当該方法は、簡便さと精度のバランスに優れている。
この場合、本発明の方法は少なくとも以下の各ステップを含む：  
（ａ２）赤外分光測定装置により、土壌サンプルの赤外分光スペクトルデータを得る。
（ｂ２）前記赤外分光スペクトルデータの標準化処理を行う。
（ｃ２）樹木高さを説明変数とし、多変量解析ＰＬＳ解析を行い、検量線モデル及び／又は検量式を作成する。
（ｄ２）前記検量線モデル及び／又は検量式により、樹木成長の推定（予測）を行う。
より具体的には、例えば赤外分光装置を用いて複数の標準土壌サンプルの３５０〜７８００ｃｍ−１の赤外スペクトルデータを取得し、最適な波数範囲として８００〜１２００ｃｍ−１を選択し、４ｃｍ−１ごとにスムージング処理を行い、スペクトルデータの標準化処理を行ってよい。樹高など樹木成長のデータ結果を用いて、ＰＬＳ回帰により検量線モデルを作成する。その検量線モデルを用いて、樹木成長が未知である場所において樹木を植栽した場合の樹高成長を推定（予測）する。
なお４ｃｍ−１ごとにスムージング処理を行うことにより、２ｃｍ−１〜１６ｃｍ−１の範囲のスムージング処理のうち最も正確な予測がなされる。
当該方法において反映される土壌の特性として、全炭素量を含むものは好ましく、交換性Ｃａ、Ｍｇ及びＫ等の微量成分の量を含むものはより好ましい。また、全Ｆｅ及び／又は全Ｓｉを反映してもよい。
上記赤外分光測定装置の形態としては室内据え置き型のもののほか、デジタルカメラや赤外スペクトルカメラといった携帯可能な小型のものを含むものが挙げられる。
携帯可能な小型のものを含む赤外分光測定装置の場合、当該赤外分光測定装置を移動機具に搭載して、地こしらえの際に当該植林地全体の土壌についての色情報を取得するようにしたものは好ましい。移動機具を用いることにより、より広範な領域についての赤外スペクトルデータ情報を正確かつ効率的に取得することができるばかりでなく、地上移動型機具の場合には表面より下方に存する土壌についての赤外スペクトルデータ情報の取得も可能になるからである。
赤外スペクトルデータ情報取得装置に加えて、得られた赤外スペクトルデータ情報を変換・解析し、解析結果を樹高に換算して樹木の成長度合いの推定・予測を行う機能を備えた機材をさらに具備した移動機具を含むシステムは、本発明の方法の実施のためにより好ましく用いられる。解析結果を樹高に換算する際に用いる換算式を環境条件に応じて改変できるシステムは、一層より好ましい。当該機材は上記測定装置の一部を構成していてもよく、あるいは上記測定装置と独立した別個の機材であってもよい。
本発明において用いられる赤外分光測定装置を搭載した移動機具として、トラクターや林業機械（ハーベスタ及びスキッダ等）といった地上移動型機具、ならびに有人航空機や無人航空機（ドローン及びラジコンヘリ等）といった空中移動型機具が例示される。移動機具として、移動パターンが規則的である、トラクター及びハーベスタは好ましく、無人航空機も操作が簡便であるため好ましい。  赤外分光測定装置に加えて色情報取得装置を搭載した本発明のシステムによれば、樹木の成長予測をより精密に行うことができるため好ましい。
【００１５】
Ｃ．土壌の養分を用いる方法
上記のとおり養分による樹木の成長予測は必ずしも容易ではないといった事情はあるにせよ、有効な栄養成分（養分）を特定し、樹木の成長予測を行うことは、一層精度の高い予測方法を与えると推測された。本発明者らはかかる推測の下、土壌の養分を用いる樹木の成長予測の方法について探究し、その結果本発明の他の態様を創出するに至った。
土壌の養分を分析することにより樹木の成長を予測する本発明の方法によれば、一般に上記Ａ．土壌の色及びＢ．土壌の特性のみを用いる方法より精度の高い予測結果を与える点において好ましい。
土壌の養分を分析することを含む本発明の方法は、少なくとも以下の各ステップを含む：  
（ａ３）全炭素、全窒素、酢酸アンモニア抽出による交換性微量金属元素（Ｃａ、Ｋ、Ｍｇ）の量（割合）を測定し、さらに土壌ｐＨ及び電気伝導度（ＥＣ）も併せて測定する。
（ｂ３）これら土壌の養分についての測定値を説明変数とし、樹高などの樹木成長の測定結果を用いて、一般化線形モデルにより、樹木成長を推定するための検量線モデル及び／又は検量式を作成する。
交換性微量金属元素として上記Ｃａ、Ｋ、Ｍｇのほか、ＡｌやＮａといった成分についての分析結果を加味することにより、一層精度の高い予測が可能となるため好ましい。ＡｌやＮａといった成分は、樹木の成長に対する影響は比較的小さいことを併せ考えれば、かかる態様が好ましいことは予想外のことである。
【００１６】
また土壌の特性及び土壌の養分を分析することの両方を含む本発明の方法も好ましい。さらにまた、土壌の色を分析することをさらに含むことにより、一層精度の高い予測が可能になる場合があり好ましい。
土壌の色、特性及び養分についての各説明変数と各種特性項目（ｐＨ、ＥＣ、元素量）についての相関を調査し、前記項目を適宜用いて検量式／検量線の作成に用いてよい。
【００１７】
本発明の方法が適用される面積はとくに限定されず、１植林地に相当する面積である数１００〜数１０００ｈａから数ｈａ以下まで適用可能である。面積に応じてサンプリング方法・数を改変してよい。
土壌のサンプリング数、サンプリングを行う場所の数（密度）及びサンプリングする土壌の量は、調査を行う土地の面請及び調査において必要な精度に応じて適宜決定してよい。
サンプルとして事前に成長調査を行う樹木の本数も、調査を行う土地の面請及び調査において必要な精度に応じて決定してよい。
【００１８】
地形の影響を低減するために、サンプリングは各形状の部位から均等に行うことは好ましい。さらに地域・地形、樹種、季節及び期間（年数）等の影響を排除するために、比較的狭いプロットでサンプリング等を行えばよい。比較的広いエリアにおいて予測を行う場合には、これらの潜在的な変動要因についてのデータを検量線に用いて行われる解析の中に入れることにより、これらの変動要因を入れ込んだ予測も可能である。
【００１９】
本発明は前記いずれかの方法により、植林予定地における樹木の成長度合いを植林のための樹木の植栽前及び／又は植栽後に予測することを含む植林方法にも関する。
当該予測に用いられる本発明の予測方法は限定されず、時間、場所、コスト等を勘案して適宜決定してよい。簡便さの観点からは土壌の色を用いる方法が有利であり、予測の正確さの観点からは土壌の養分の分析結果を用いる方法が有利である。また土壌の特性を用いる方法は、土壌の色を用いる方法及び土壌の養分の分析結果を用いる方法の両方の利点を相当程度併せて有するという優位性がある。
【００２０】
本発明の方法を行うに際して用いられる機材は限定されないところ、以下を含む、本発明の方法を行うためのシステムは好ましい：
（ア）土壌の色及び／又は土壌の特性を測定する測定装置；及び
（イ）前記装置を搭載した移動機具。
かかるシステムのうち、測定装置を用いて得られた測定値を樹高に換算して樹木の成長度合いを予測する機能を備えた機具をさらに具備するシステムはより好ましく、移動機具が地上移動型機具又は空中移動型機具であるシステムはより一層好ましい。
上記システムの例として、色情報取得装置及び／又は赤外分光測定装置を搭載した移動機具を移動手段として含むものが挙げられる。当該移動機具として、トラクターや林業機械（ハーベスタ及びスキッダ等）といった地上移動型機具、ならびに有人航空機や無人航空機（ドローン及びラジコンヘリ等）といった空中移動型機具が例示され、トラクター及びハーベスタは好ましく、無人航空機も好ましい。
【実施例】
【００２１】
次に実施例に基づいて本発明を更により詳細に説明するが、本発明はこれらの実施例によって何等制限されるものではない。【００２２】
（実施例１）
樹木の成長（樹高）を予測する方法についての試験（１）
[試験方法] 
・試験地  ベトナム国ロンアン省内のアカシア・ハイブリッド植林地１年生。
・同一プロット内で樹木の成長にバラつき（３０ｃｍ〜４３０ｃｍ）が生じていた。土壌成分の相違が原因と思われたが、明確には判じられなかった。 
・樹木下の土壌０〜５ｃｍをサンプリングした（合計２２地点）
・土壌の色（以下「Ａ法」ということがある）、土壌の特性（以下「Ｂ法」ということがある）又は土壌の養分（以下「Ｃ法」ということがある）を用いる方法のそれぞれにつき、説明変数と樹高との相関を調査した。Ａ法〜Ｃ法のそれぞれにおいて用いられた手法は本明細書本文において上記したとおりである。【００２３】
[結果]
Ａ法、Ｂ法及びＣ法を用いた結果を以下に示す。
（１）Ａ法
ｐＨ及びＥＣを考慮せずＬ*、ａ*及びｂ*のみを用いた場合、予測高さと実測高さとの決定係数（Ｒ２）は０．４４０７であった（図２Ａ）。
また、一般化線形モデルにより2次までの交互作用を含めて作成した予測式から得られた予測高さと実測高さとの決定係数は０．５４７８に向上し（図２Ｂ）、さらにｐＨ及びＥＣの測定値を考慮すると決定係数は０．７９８７にさらに向上した（図２Ｃ）。
（２）Ｂ法
予測高さと実測高さとの決定係数は０．８４１９であった（図３）。なお５つのサンプルについての赤外スペクトルデータ及びそのスムージング処理後のデータ（チャート）を図４Ａ及び図４Ｂに示す。
（３）Ｃ法
土壌養分として通常用いられることの多い、ｐＨ、ＥＣ、全炭素、全窒素、交換性Ｃａ、Ｍｇ、Ｋを用いて一般化線形モデルにより、予測式を作成した。交互作用を含まない１次式により作成した場合、決定係数は０．６５９９であった（「養分１法」、図５Ａ）。
次に土壌養分として通常用いられることの多い、ｐＨ、ＥＣ、全炭素、全窒素、交換性Ｃａ、Ｍｇ、Ｋに加えて、通常では植物の成長に影響が少なかったり、有害である成分も含めて、一般化線形モデルにより、予測式を作成した。交互作用を含まない１次式により作成した。養分１法に用いられた養分種に加えて、Ａｌ、Ｂ、Ｆｅ、Ｎａ、Ｓｉ、Ｚｎを説明変数に加えたところ、決定係数は０．８８１６となりより予測精度の高い測定が可能となった（「養分２法」、図５Ｂ）。
なお、養分２法における回帰直線の作成に用いられた各養分種に割り当てられた係数は下表のとおりである（表１）。切片についても同様に下表に挿入して数値を示した。【００２４】
【表１】
2016043007.tif 000003
【００２５】
（実施例２）樹木の成長（樹高）を予測する方法についての試験（２）
実施例１と同様な試験をパプア・ニューギニアのユーカリ（７年生から２０年生）について行った。
[試験方法] 
・試験地：パプア・ニューギニア国のユーカリ植林地７年生から２０年生。
・１５年生時の主林木の平均樹高について地位指数４を平均とし、２ｍごとに１〜７までの７段階の樹高に区分し、樹高の実測値と赤外分光法による予測値とを比較することとした。なお、「地位指数」は、あるプロットにおける主林木の平均樹高をより簡便に表すために本試験において採用された指標である。
・土壌供試サンプルとして、樹木下の土壌０〜１０ｃｍの深さの地点からプロット当たり５点採取し、混合したものを各プロットにおける土壌供試サンプルとして用いた。
・土壌供試サンプルについて赤外分光スペクトルデータを得た後、前記赤外分光スペクトルデータの標準化処理を行った。さらに、樹木高さを説明変数とし、多変量解析ＰＬＳ解析を行い、検量式を作成した。前記検量式により、説明変数と樹高との相関を表す指標としての決定係数を計算し、樹木成長の推定（予測）の可否を判断した。【００２６】
[結果]
樹高についての説明変数（予測値）と実測値との決定係数は０．４８５２であり（図６）、ユーカリについても成長が予測できることが明らかになった。
土壌供試サンプルのサンプリングをより深い地点の土壌も併せて行うことにより、より精度の高い成長の予測が可能になることも示唆された。
【産業上の利用可能性】
【００２７】
本発明によれば、従来の方法が達成し得なかった植林前に樹木の成長を予測する方法が提供される。したがって、本発明は、植林産業および関連産業の発展に寄与するところ大である。
【図１】
2016043007.tif 000004
【図２Ａ】
2016043007.tif 000005
【図２Ｂ】
2016043007.tif 000006
【図２Ｃ】
2016043007.tif 000007
【図３】
2016043007.tif 000008
【図４Ａ】
2016043007.tif 000009
【図４Ｂ】
2016043007.tif 000010
【図５Ａ】
2016043007.tif 000011
【図５Ｂ】
2016043007.tif 000012
【図６】
2016043007.tif 000013
【国際調査報告】
2016043007.tif 000014
2016043007.tif 000015
2016043007.tif 000016
2016043007.tif 000017