JP7464284B2 - プラズマ発光分光分析を用いた土壌診断方法 - Google Patents
プラズマ発光分光分析を用いた土壌診断方法 Download PDFInfo
- Publication number
- JP7464284B2 JP7464284B2 JP2021092618A JP2021092618A JP7464284B2 JP 7464284 B2 JP7464284 B2 JP 7464284B2 JP 2021092618 A JP2021092618 A JP 2021092618A JP 2021092618 A JP2021092618 A JP 2021092618A JP 7464284 B2 JP7464284 B2 JP 7464284B2
- Authority
- JP
- Japan
- Prior art keywords
- soil
- exchangeable
- samples
- data
- measured
- 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.)
- Active
Links
Images
Landscapes
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Description
本発明は、プラズマ発光分光分析を用いた土壌診断方法に関する。
土壌分析は作物の生産性向上を図る上で必須である。土壌診断に加えて、普及の進む精密農業・スマート農業においても各圃場での土壌特性の把握が必要であり、今後、土壌分析の需要はさらに高まるものと考えられる。
土壌試料の多項目を同時に診断する方法としては、下記の技術が報告されている。
（１）複数の抽出液を用いた土壌試料の多項目同時診断技術
馬場氏らによりディスクリート式自動化学分析装置による土壌診断技術が開発され（非特許文献１参照）、土壌試料の多項目同時診断が実施された。当該技術は、土壌試料の抽出および比色法による多項目同時測定を可能とするもので、ディスクリート式自動化学分析装置として市場化されている（非特許文献２参照）。この方法では、各試料に対して土壌分析の公定法に準じた抽出液ならびに発色液を用いるため、高精度の分析が可能であり、硝酸態窒素量、アンモニア態窒素量、有効態リン（トルオーグ法、Bray法）、交換性塩基（Ｃａ、Ｍｇ、Ｋ）、ＣＥＣ、鉄、ケイ酸、ホウ素量、腐植量を測定可能である。
（１）複数の抽出液を用いた土壌試料の多項目同時診断技術
馬場氏らによりディスクリート式自動化学分析装置による土壌診断技術が開発され（非特許文献１参照）、土壌試料の多項目同時診断が実施された。当該技術は、土壌試料の抽出および比色法による多項目同時測定を可能とするもので、ディスクリート式自動化学分析装置として市場化されている（非特許文献２参照）。この方法では、各試料に対して土壌分析の公定法に準じた抽出液ならびに発色液を用いるため、高精度の分析が可能であり、硝酸態窒素量、アンモニア態窒素量、有効態リン（トルオーグ法、Bray法）、交換性塩基（Ｃａ、Ｍｇ、Ｋ）、ＣＥＣ、鉄、ケイ酸、ホウ素量、腐植量を測定可能である。
（２）単一の土壌抽出液を用いた多項目同時診断技術
土壌試料からの単一抽出物を用いた多項目同時診断技術として、単一の土壌抽出液を用いた多項目同時診断技術の報告は少ないが、その中でエア・ウォーター・バイオデザイン株式会社では、独自の成分抽出液を用いて抽出し、ＬＥＤ／フォトディテクタによる吸光計測によって、硝酸態窒素量、アンモニア態窒素量、可給態リン酸、交換性塩基（Ｃａ、Ｍｇ、Ｋ）の６項目の同時測定が可能であるとしている（非特許文献３参照）。
土壌試料からの単一抽出物を用いた多項目同時診断技術として、単一の土壌抽出液を用いた多項目同時診断技術の報告は少ないが、その中でエア・ウォーター・バイオデザイン株式会社では、独自の成分抽出液を用いて抽出し、ＬＥＤ／フォトディテクタによる吸光計測によって、硝酸態窒素量、アンモニア態窒素量、可給態リン酸、交換性塩基（Ｃａ、Ｍｇ、Ｋ）の６項目の同時測定が可能であるとしている（非特許文献３参照）。
（３）土壌試料に対し近赤外線を照射して近赤外分光法により分析する方法
この方法は、多項目同時土壌診断技術として、土壌試料に対し近赤外線を照射し、近赤外域における成分別の分子振動による光吸収や拡散反射特性から、土壌中の炭素、窒素などの肥沃度関連指標を予測する技術である（特許文献１及び２、非特許文献４参照）。近赤外波長を用いた分析では、全炭素、全窒素、可給態窒素、ＣＥＣの分析が可能であることが示された。
この方法は、多項目同時土壌診断技術として、土壌試料に対し近赤外線を照射し、近赤外域における成分別の分子振動による光吸収や拡散反射特性から、土壌中の炭素、窒素などの肥沃度関連指標を予測する技術である（特許文献１及び２、非特許文献４参照）。近赤外波長を用いた分析では、全炭素、全窒素、可給態窒素、ＣＥＣの分析が可能であることが示された。
（４）プラズマ発光分光分析を用いた方法
誘導結合プラズマによる発光を分析するプラズマ発光分光分析(Inductively Coupled Plasma－Atomic Emission Spectrometry、ＩＣＰ－ＡＥＳと呼ぶ）は、液体試料中の多元素同時測定が可能な汎用元素分析装置である。土壌分析においては、各種溶媒で抽出された元素濃度の分析に用いられる。例えば、土壌中の交換性塩基についてはｐＨ７の１ＭＮＨ４ＯＡｃ（酢酸アンモニウム）で、また交換性Ａｌについては１ＭＫＣｌで抽出し、ＩＣＰ－ＡＥＳにより抽出液中の元素濃度を定量して、土壌あたりの含量として算出している。この際、ＩＣＰ－ＡＥＳでは対象元素の測定波長を指定して濃度を算出している。
誘導結合プラズマによる発光を分析するプラズマ発光分光分析(Inductively Coupled Plasma－Atomic Emission Spectrometry、ＩＣＰ－ＡＥＳと呼ぶ）は、液体試料中の多元素同時測定が可能な汎用元素分析装置である。土壌分析においては、各種溶媒で抽出された元素濃度の分析に用いられる。例えば、土壌中の交換性塩基についてはｐＨ７の１ＭＮＨ４ＯＡｃ（酢酸アンモニウム）で、また交換性Ａｌについては１ＭＫＣｌで抽出し、ＩＣＰ－ＡＥＳにより抽出液中の元素濃度を定量して、土壌あたりの含量として算出している。この際、ＩＣＰ－ＡＥＳでは対象元素の測定波長を指定して濃度を算出している。
馬場康尋、後藤逸男、「ディスクリート方式による自動化学分析装置の土壌診断への応用」、日本土壌科学雑誌、第８０巻、第６号、p.611-615（2009）
富士平工業株式会社、高速土壌養分自動分析装置、ＳＮＡ－３０ｉ、http://www.fujihira.co.jp/
エア・ウォーター・バイオデザイン株式会社、土壌分析装置、ＥＷ－ＴＨＡ１Ｊ、https://www.awbio.co.jp/soil-analyzer/
森次真一・鷲尾建紀・高原知佳子・大塚理哉・高野和夫、「近赤外分光法による土壌化学性診断の可能性」、岡山県農業研報６：41-48（２０１５）
上記のように土壌分析は種々の方法が検討されているが、公定法による土壌分析は依然として煩雑な作業を要し、一検体あたりの分析コストが高いという問題を抱えている。
上記（１）の複数の抽出液を用いた土壌試料の多項目同時診断技術では、高精度の分析が可能であるが、数種類の抽出液及び発色液を使用するので、抽出液及び発色液の準備等に時間が掛かるという課題がある。
上記（１）の複数の抽出液を用いた土壌試料の多項目同時診断技術では、高精度の分析が可能であるが、数種類の抽出液及び発色液を使用するので、抽出液及び発色液の準備等に時間が掛かるという課題がある。
上記（２）の単一の土壌抽出液を用いた多項目同時診断技術では、上記の硝酸態窒素量、アンモニア態窒素量、可給態リン酸、交換性塩基（Ｃａ、Ｍｇ、Ｋ）の６項目の同時測定が可能であるが、他のｐＨや粒径などの多くの分析項目の予測はできないという課題がある。
上記（３）の近赤外分光法により分析する方法では、全炭素、全窒素、可給態窒素、ＣＥＣの分析が可能であるが、その他の分析項目では予測精度が低く、予測出来ない項目も多いという課題がある。
上記（３）の近赤外分光法により分析する方法では、全炭素、全窒素、可給態窒素、ＣＥＣの分析が可能であるが、その他の分析項目では予測精度が低く、予測出来ない項目も多いという課題がある。
上記（４）のＩＣＰ－ＡＥＳにより分析する方法では、交換性塩基や交換性Ａｌのそれぞれを特定の抽出液を使用して分析をしているので、単一の土壌抽出液を用いた多項目同時診断ができないという課題がある。
以上説明したように、従来の土壌分析の何れの方法でも、単一の土壌抽出液を用い、例えば６項目を超える多項目の同時診断ができないという課題がある。
本発明者等は、ＩＣＰ－ＡＥＳによる土壌分析において、単一の土壌抽出液を用い、かつ、発光分析において、特定の波長だけではなく、全波長スペクトルデータを取得して、深層学習を用いた機械学習により、土壌試料の多項目の同時診断ができることを見出して本発明に到達した。
本発明は、上記課題に鑑み、ＩＣＰ－ＡＥＳにおいて、全波長スペクトルデータを用いて、交換性塩基、土壌ｐＨ(Ｈ２Ｏ抽出、ＫＣｌ抽出)、電気伝導率、有効態リン、全炭素、全窒素、交換性Ａｌ、ＣＥＣ、粒径組成（粘土画分量、砂画分量）等の土壌の種々の理化学特性に関する分析項目が取得できるプラズマ発光分光分析を用いた土壌診断方法を提供することを目的とする。なお、本発明において土壌の「診断」とは、必ずしも分析結果の評価を意味するものではなく、土壌の種々の理化学特性に関する各分析すべき項目を、分析し、予測し及び測定し又はこれを評価することの何れかを示す概念として使用しており、作物の生産性向上に向けた土壌診断や土壌分類のための基礎情報を取得することを総称している。
上記目的を達成するため、本発明のプラズマ発光分光分析を用いた土壌診断方法は、
学習用データ取得のための被測定土壌を所定の抽出液を用いて土壌抽出液とする第一のステップと、
プラズマ発光分光分析により土壌抽出液の全波長スペクトルデータを取得すると共に、プラズマ発光分光分析により被測定土壌のＣａ、Ｍｇ、Ｎａ、Ｋからなる交換性塩基の実測値を取得する第二のステップと、
交換性塩基以外の分析項目の実測値を公定法により取得する第三のステップと、
全波長スペクトルデータと交換性塩基及び前記交換性塩基以外の分析項目の実測値とを用いて、深層学習により学習を行い、前記学習の後で診断すべき被測定土壌の所定の分析項目の予測値を算出する第四のステップと、を、備える。
学習用データ取得のための被測定土壌を所定の抽出液を用いて土壌抽出液とする第一のステップと、
プラズマ発光分光分析により土壌抽出液の全波長スペクトルデータを取得すると共に、プラズマ発光分光分析により被測定土壌のＣａ、Ｍｇ、Ｎａ、Ｋからなる交換性塩基の実測値を取得する第二のステップと、
交換性塩基以外の分析項目の実測値を公定法により取得する第三のステップと、
全波長スペクトルデータと交換性塩基及び前記交換性塩基以外の分析項目の実測値とを用いて、深層学習により学習を行い、前記学習の後で診断すべき被測定土壌の所定の分析項目の予測値を算出する第四のステップと、を、備える。
上記構成において、抽出液は、好ましくは酢酸アンモニウム（ＮＨ４ＯＡｃ）溶液である。
被測定土壌の分析項目は、好ましくは、交換性塩基（Ｃａ、Ｍｇ、Ｎａ、Ｋ）、ｐＨ、電気伝導率、有効態リン、全炭素、全窒素、交換性Ａｌ、ＣＥＣ、粒径組成、可溶性Ｆｅ、交換性Ｍｎ、可溶性Ｚｎ、リン酸吸収係数、硝酸態窒素量、アンモニア態窒素量を含む。粒径組成は、好ましくは粘土画分量及び砂画分量である。
第二のステップにおいて、好ましくは、全波長スペクトルデータの強度を自然対数に変換し、かつ、全波長スペクトルデータの強度の内、飽和した強度とマイナスの強度とをそれぞれ所定の値に設定して自然対数に変換する。
さらに本発明は、深層学習を終了してよいか否かの判定を行う第五のステップを有し、
この第五のステップにおいて、深層学習を終了してよいと判定された場合には、深層学習に用いた学習用データをデータセットとして格納する第六のステップと、
新規の被測定土壌を所定の抽出液を用いて新規の土壌抽出液とする第七のステップと、
プラズマ発光分光分析により新規の土壌抽出液の全波長スペクトルデータを取得する第八のステップと、
第八のステップで取得した全波長スペクトルデータと、第五のステップで取得した学習用データとを用いて、新規の被測定土壌の所定の分析項目の予測値を算出する第九のステップとを、備えてもよい。
好ましくは、第八のステップで取得した全波長スペクトルデータと第九のステップで取得した新規の被測定土壌の所定の分析項目の予測値とを新規な学習用データとする。
被測定土壌の分析項目は、好ましくは、交換性塩基（Ｃａ、Ｍｇ、Ｎａ、Ｋ）、ｐＨ、電気伝導率、有効態リン、全炭素、全窒素、交換性Ａｌ、ＣＥＣ、粒径組成、可溶性Ｆｅ、交換性Ｍｎ、可溶性Ｚｎ、リン酸吸収係数、硝酸態窒素量、アンモニア態窒素量を含む。粒径組成は、好ましくは粘土画分量及び砂画分量である。
第二のステップにおいて、好ましくは、全波長スペクトルデータの強度を自然対数に変換し、かつ、全波長スペクトルデータの強度の内、飽和した強度とマイナスの強度とをそれぞれ所定の値に設定して自然対数に変換する。
さらに本発明は、深層学習を終了してよいか否かの判定を行う第五のステップを有し、
この第五のステップにおいて、深層学習を終了してよいと判定された場合には、深層学習に用いた学習用データをデータセットとして格納する第六のステップと、
新規の被測定土壌を所定の抽出液を用いて新規の土壌抽出液とする第七のステップと、
プラズマ発光分光分析により新規の土壌抽出液の全波長スペクトルデータを取得する第八のステップと、
第八のステップで取得した全波長スペクトルデータと、第五のステップで取得した学習用データとを用いて、新規の被測定土壌の所定の分析項目の予測値を算出する第九のステップとを、備えてもよい。
好ましくは、第八のステップで取得した全波長スペクトルデータと第九のステップで取得した新規の被測定土壌の所定の分析項目の予測値とを新規な学習用データとする。
本発明によれば、プラズマ発光分光分析によって全波長スペクトルデータを用いた分析を行うことにより、精度良く迅速に各種土壌分析項目を診断でき、優れたプラズマ発光分光分析を用いた土壌診断方法を提供することができる。
以下、図面を参照しながら本発明を実施するための形態について説明するが、本発明の範囲はここに記載する実施形態に限定されることなく適宜変更することができる。特に、図面に記載した各部の名称などについては概念的な事項を示すに過ぎず、その適用場面に応じて適宜変更することができる。
図１は、本発明の一実施形態に係る土壌診断方法に用いる土壌析装置１の構成を示すブロック図である。土壌分析装置１は、ＩＣＰ分析部１０とＩＣＰ分析部１０に接続される土壌分析部２０とから構成される。
ＩＣＰ分析部１０は、誘導結合プラズマによる発光を分析するプラズマ発光分光分析が取得できる装置であり、所謂ＩＣＰ発光分光分析装置と、例えば７１元素２３２波長等の後述する全波長スペクトルデータを取得できるソフトウェアが格納されたＰＣやワークステーションからなるコンピュータと、で構成されている。
ＩＣＰ分析部１０は、誘導結合プラズマによる発光を分析するプラズマ発光分光分析が取得できる装置であり、所謂ＩＣＰ発光分光分析装置と、例えば７１元素２３２波長等の後述する全波長スペクトルデータを取得できるソフトウェアが格納されたＰＣやワークステーションからなるコンピュータと、で構成されている。
土壌分析部２０は、ＩＣＰ分析部１０からのＩＣＰ発光スペクトルデータを処理するＩＣＰ発光スペクトル処理部２１と、ＩＣＰ発光スペクトル処理部２１で処理された処理済ＩＣＰ発光スペクトルデータが入力される制御部２２と、制御部２２に接続されるデータ処理部３０とから構成されている。
データ処理部３０は、制御部２２とインターフェースをする入出力インターフェース部３１と、主記憶装置及び補助記憶装置を備える記憶装置３２と、四則演算等の演算処理を行う演算装置と、記憶装置及び演算装置を制御する制御装置とを備えるコンピュータで構成され、データ処理プログラムが記憶装置３２に格納されており、データ処理プログラムが演算装置に展開されて実行される。これにより、データ処理部３０は図１に示すような記憶部３３及び予測部３４を機能的に備える。
土壌分析部２０は具体的には、ＣＰＵ（Central Processing Unit）、揮発性メモリのＤＲＡＭ（ダイナミックランダムアクセスモリ）、不揮発性のメモリであるＨＤＤ（ハードディスクドライブ）等からなるＰＣやワークステーション等のコンピュータで構成される。ＣＰＵはマイクロプロセッサにより構成することができる。さらに、データ処理部３０は、後述する予測部３１で実行する深層学習のために多数のデータ処理を高速に行うよう、ＧＰＵ（Graphics Processing Unit）と、記憶装置３２として大容量のＤＲＡＭを備えて構成してもよい。
土壌分析部２０において、ＩＣＰ発光スペクトル処理部２１で処理された処理済ＩＣＰ発光スペクトルデータのデジタルデータから、予測部３４において診断すべき土壌の土壌地力形質の分析項目を予測（推定とも呼ぶ）する。
図２は、本発明の土壌析装置１を用いて土壌診断方法を行うときの手順を示すフロー図である。
第一のステップ（ＳＴ１）において、学習用データ取得のための被測定土壌を所定の抽出液を用いて土壌抽出液とする。つまり、土壌抽出液を調製する。
第二のステップ（ＳＴ２）において、ＩＣＰ分析部１０でプラズマ発光分光分析により土壌抽出液の全波長スペクトルデータを取得すると共に、プラズマ発光分光分析により被測定土壌のＣａ、Ｍｇ、Ｎａ、Ｋからなる交換性塩基の実測値を取得する。
第三のステップ（ＳＴ３）において、被測定土壌の交換性塩基（Ｃａ、Ｍｇ、Ｎａ、Ｋ）以外の分析項目の実測値を公定法により取得する。
第四のステップ（ＳＴ４）において、ＩＣＰ分析部１０で取得した全波長スペクトルデータと、被測定土壌の交換性塩基及び交換性塩基以外の分析項目を含む所定の分析項目の実測値を用いて、深層学習により学習を行い、学習の後で診断すべき被測定土壌の所定の分析項目の予測値を算出する。
第一のステップ（ＳＴ１）において、学習用データ取得のための被測定土壌を所定の抽出液を用いて土壌抽出液とする。つまり、土壌抽出液を調製する。
第二のステップ（ＳＴ２）において、ＩＣＰ分析部１０でプラズマ発光分光分析により土壌抽出液の全波長スペクトルデータを取得すると共に、プラズマ発光分光分析により被測定土壌のＣａ、Ｍｇ、Ｎａ、Ｋからなる交換性塩基の実測値を取得する。
第三のステップ（ＳＴ３）において、被測定土壌の交換性塩基（Ｃａ、Ｍｇ、Ｎａ、Ｋ）以外の分析項目の実測値を公定法により取得する。
第四のステップ（ＳＴ４）において、ＩＣＰ分析部１０で取得した全波長スペクトルデータと、被測定土壌の交換性塩基及び交換性塩基以外の分析項目を含む所定の分析項目の実測値を用いて、深層学習により学習を行い、学習の後で診断すべき被測定土壌の所定の分析項目の予測値を算出する。
（第一のステップ）
抽出液は、土壌から交換性塩基（Ｃａ,Ｍｇ,Ｋ,Ｎａ）等の陽イオンを抽出できる液が好ましく、例えば所定濃度で所定のｐＨを有している酢酸アンモニウム溶液を用いることができる。酢酸アンモニウム溶液の濃度は、例えば０．０５Ｍから２Ｍの濃度範囲で土壌抽出液を得ることができる。さらに好ましくは、０．９Ｍ（モル）から１Ｍ程度とすることができる。酢酸アンモニウム溶液の濃度は、１Ｍ程度とすることが好ましい。酢酸アンモニウム溶液のｐＨは、例えば６．５～７とすることができ、好ましくはｐＨは７である。
抽出液は、土壌から交換性塩基（Ｃａ,Ｍｇ,Ｋ,Ｎａ）等の陽イオンを抽出できる液が好ましく、例えば所定濃度で所定のｐＨを有している酢酸アンモニウム溶液を用いることができる。酢酸アンモニウム溶液の濃度は、例えば０．０５Ｍから２Ｍの濃度範囲で土壌抽出液を得ることができる。さらに好ましくは、０．９Ｍ（モル）から１Ｍ程度とすることができる。酢酸アンモニウム溶液の濃度は、１Ｍ程度とすることが好ましい。酢酸アンモニウム溶液のｐＨは、例えば６．５～７とすることができ、好ましくはｐＨは７である。
（第二のステップ）
後述するＩＣＰ分析部１０で用いるソフト（ＩＣＰＥ－９０００ソフト（ＩＣＰＥ ｓｏｌｕｔｉｏｎｓ）のデフォルトで設定されている７１元素２３２波長に加えて、炭素の感度情報のある２波長を加えた２３４波長を対象とした。ここで、ＩＣＰ分析部１０で取得される２３４波長のスペクトルデータを全波長スペクトルデータと呼ぶ。発光を検出するＣＣＤ検知器においては、例えば各測定波長を中央とした、３１ピクセル分のピークプロファイルデータが得られる。各ピークプロファイルデータに対して、ピーク強度計算の積分範囲を１ピクセルとなるように指定し、中央ピクセルからの相対位置が０,３,６,９,１２,１５,－３,－６,－９,－１２,－１５となる１１ピクセルにおける個別の強度値を求め、得られた強度値の自然対数を算出し説明変数（ｘ）とした（２５７４ピクセル/試料）。従属変数（ｙ）は土壌地力形質の各分析項目の実測値を用いた。
後述するＩＣＰ分析部１０で用いるソフト（ＩＣＰＥ－９０００ソフト（ＩＣＰＥ ｓｏｌｕｔｉｏｎｓ）のデフォルトで設定されている７１元素２３２波長に加えて、炭素の感度情報のある２波長を加えた２３４波長を対象とした。ここで、ＩＣＰ分析部１０で取得される２３４波長のスペクトルデータを全波長スペクトルデータと呼ぶ。発光を検出するＣＣＤ検知器においては、例えば各測定波長を中央とした、３１ピクセル分のピークプロファイルデータが得られる。各ピークプロファイルデータに対して、ピーク強度計算の積分範囲を１ピクセルとなるように指定し、中央ピクセルからの相対位置が０,３,６,９,１２,１５,－３,－６,－９,－１２,－１５となる１１ピクセルにおける個別の強度値を求め、得られた強度値の自然対数を算出し説明変数（ｘ）とした（２５７４ピクセル/試料）。従属変数（ｙ）は土壌地力形質の各分析項目の実測値を用いた。
ＩＣＰ分析部１０では、プラズマ発光分光分析により土壌地力形質の内、被測定土壌の交換性塩基（Ｃａ、Ｍｇ、Ｎａ、Ｋ）の実測値を取得して、ＩＣＰ発光スペクトル処理部２１に出力する。
次に、ＩＣＰ分析部１０からＩＣＰ発光スペクトル処理部２１に入力された全波長スペクトルデータは、ＩＣＰ発光スペクトル処理部２１で処理が施されて処理済ＩＣＰ発光スペクトルデータとされる。全波長スペクトルデータは、その強度のダイナミックレンジが大きいので、処理を高速に行うために全波長スペクトルデータの強度を全て自然対数に変換した。自然対数には、飽和した強度とマイナスの強度をそれぞれ所定の値に設定してから変換することができる。具体例を以下に示す。
（ａ）ＩＣＰ発光スペクトルの各波長の強度値は全て自然対数に変換される。
（ｂ）ＩＣＰ発光スペクトルの発光強度飽和の場合、例えば強度値１００００００と仮定して、自然対数に変換する。
（ｃ）得られた強度値がマイナスの試料は強度値を１として、自然対数に変換する。
次に、ＩＣＰ分析部１０からＩＣＰ発光スペクトル処理部２１に入力された全波長スペクトルデータは、ＩＣＰ発光スペクトル処理部２１で処理が施されて処理済ＩＣＰ発光スペクトルデータとされる。全波長スペクトルデータは、その強度のダイナミックレンジが大きいので、処理を高速に行うために全波長スペクトルデータの強度を全て自然対数に変換した。自然対数には、飽和した強度とマイナスの強度をそれぞれ所定の値に設定してから変換することができる。具体例を以下に示す。
（ａ）ＩＣＰ発光スペクトルの各波長の強度値は全て自然対数に変換される。
（ｂ）ＩＣＰ発光スペクトルの発光強度飽和の場合、例えば強度値１００００００と仮定して、自然対数に変換する。
（ｃ）得られた強度値がマイナスの試料は強度値を１として、自然対数に変換する。
（第三のステップ）
被測定土壌の分析項目は、土壌地力形質とも呼ばれ、従属変数とも呼ばれている。交換性塩基（Ｃａ、Ｍｇ、Ｎａ、Ｋ）以外の分析項目は、土壌ｐＨ（Ｈ２Ｏ抽出、ＫＣｌ抽出)、土壌ＥＣ、有効態リン（Bｒａｙ－Ｉ法）、全炭素、全窒素、交換性Ａｌ、ＣＥＣ、粒径組成（粘土画分量、砂画分量）等である。
被測定土壌の分析項目は、土壌地力形質とも呼ばれ、従属変数とも呼ばれている。交換性塩基（Ｃａ、Ｍｇ、Ｎａ、Ｋ）以外の分析項目は、土壌ｐＨ（Ｈ２Ｏ抽出、ＫＣｌ抽出)、土壌ＥＣ、有効態リン（Bｒａｙ－Ｉ法）、全炭素、全窒素、交換性Ａｌ、ＣＥＣ、粒径組成（粘土画分量、砂画分量）等である。
ここで、土壌ＥＣは、土壌の電気伝導率である。ＣＥＣ（Cation Exchange Capacity、交換性塩基容量と呼ぶ）の値は、電気的にマイナスの土壌が最大限どの程度の塩基（Ｃａ、Ｍｇ、Ｎａ、Ｋ、アンモニウム、Ｈ等）の陽イオンを吸着することができるかの指標である。ＣＥＣは、土壌に肥効養分を蓄えておける量及び緩衝力を示す指標となる数値である。
交換性塩基は、第二のステップにおいて、ＩＣＰ分析部１０で発光分析されるが、上記の土壌ｐＨから粒径組成の各分析項目は、後述する公定法により実測値を取得する。
交換性塩基は、第二のステップにおいて、ＩＣＰ分析部１０で発光分析されるが、上記の土壌ｐＨから粒径組成の各分析項目は、後述する公定法により実測値を取得する。
公定法による実測値としては、所謂従来の化学分析を用いた定量分析（定量法とも呼ぶ）又は、ＩＣＰ分析部１０に用いるＩＣＰ分析装置やエネルギー分散型蛍光Ｘ線分析等のような定量分析が可能な装置により取得できる。
（第四のステップ）
第四のステップの深層学習について説明する。図３は、予測部３４において土壌地力形質の推定を行う深層学習を模式的に説明する図である。
予測部３４で実行される深層学習は、図３に示すように、処理済ＩＣＰ発光スペクトルデータを入力データ４１とする入力層４２、中間層４３及び出力層４４を備えたニューラルネットワークから構成されており、例えばコンピュータ上で動作する機械学習や深層学習のためのニューラルネットワーク用のプログラムにより動作する。このようなニューラルネットワーク用のプログラムを作成するソフトウェアとして、Neural Network Console version 1.7.7352.44102（Sony Network Communications Inc.製）等を使用することができる。
第四のステップの深層学習について説明する。図３は、予測部３４において土壌地力形質の推定を行う深層学習を模式的に説明する図である。
予測部３４で実行される深層学習は、図３に示すように、処理済ＩＣＰ発光スペクトルデータを入力データ４１とする入力層４２、中間層４３及び出力層４４を備えたニューラルネットワークから構成されており、例えばコンピュータ上で動作する機械学習や深層学習のためのニューラルネットワーク用のプログラムにより動作する。このようなニューラルネットワーク用のプログラムを作成するソフトウェアとして、Neural Network Console version 1.7.7352.44102（Sony Network Communications Inc.製）等を使用することができる。
中間層４３は入力層４２と出力層４４との間に挿入され、入力層４２から予測用の学習用データが入力されると、出力層４４から出力データ４５が出力される。学習用データは教師データとも呼ばれている。
予測部３４においては、深層学習による学習と予測とにより、各分析項目の予測が行われ、以下のステップ２０の演算により学習が行われ、ステップ３０の演算により評価用データから各分析項目の予測が行われる。
（ステップ２０）
ＩＣＰ発光スペクトル処理部２１で処理された処理済ＩＣＰ発光スペクトルデータ４１が入力層４２に入力される。入力層４２は、２３４波長に対応する２３４個の入力ユニット４２ａを有しており、各入力ユニット４２ａのそれぞれに、各波長で得られた発光強度値の自然対数が説明変数（ｘ）として入力される。
ＩＣＰ発光スペクトル処理部２１で処理された処理済ＩＣＰ発光スペクトルデータ４１が入力層４２に入力される。入力層４２は、２３４波長に対応する２３４個の入力ユニット４２ａを有しており、各入力ユニット４２ａのそれぞれに、各波長で得られた発光強度値の自然対数が説明変数（ｘ）として入力される。
予測部３４は、入力層４２の各入力ユニット４２ａに入力された入力データ４１から畳み込み演算を実行して、重み係数が計算される。この重み係数は、後述するニューラルネットワークの機械学習の結果、つまり学習結果によって決定される。
中間層４３は、各中間値に対応するｎ個の中間ユニット４３ａを有しており、各中間ユニット４３ａには、それぞれ個々の中間値ｙｉｊ（１≦ｉ≦ｎ、及び１≦ｊ≦ｍ）が入力される。
予測部３４は、中間層４３の各中間ユニット４３ａに入力された中間値ｙｉｊから演算し、重み係数を計算して、ニューラルネットワークの機械学習によって出力層４４から出力データ４５が出力される。
予測部３４は、中間層４３の各中間ユニット４３ａに入力された中間値ｙｉｊから演算し、重み係数を計算して、ニューラルネットワークの機械学習によって出力層４４から出力データ４５が出力される。
ニューラルネットワークの機械学習において、学習用データとして、ｐＨ７の１Ｍの酢酸アンモニウム溶液による土壌抽出液をＩＣＰ発光分光分析装置で測定して取得した全波長スペクトルデータと各分析値の実測値とを用いることができる。
上記ステップ２０により、ｐＨが７で、１Ｍ酢酸アンモニウム溶液による土壌抽出液を、マルチタイプＩＣＰ発光分光分析装置を用いて取得した全波長スペクトルと、各実測値を、人工知能（ＡＩ（Artificial Intelligence））による深層学習（Deep Learning）により学習する。即ち、予測部３４は、学習のための既知のデータとして、処理済ＩＣＰ発光スペクトルデータ４１が入力層４２に入力されると共に、教師信号として、出力データ４５が得られるように重み係数を演算し、決定する。そして、出力データ４５として得られた各土壌分析値の予測値との比較により、予測精度を検証する。ここで、入力される既知のデータ数、つまり学習用データが多いほど、決定される重み係数の精度が高まる。
（ステップ３０）
このようなニューラルネットワークの機械学習の結果、診断すべき被測定土壌を測定したデータ、つまり評価用データが処理済ＩＣＰ発光スペクトルデータ４１として入力されると、予測部３３は、中間データそして出力データ４５を二段階で重み係数を演算することにより、出力データ４５を出力する。即ち、出力データ４５としては、評価用データの各分析項目を出力する。このようなニューラルネットワークの学習結果に基づく各分析項目の出力データ４５は、予測値又は推定値とも呼ばれている。
従って、予測部３４の出力層４４からの出力データ４５から、ｘ軸に実測値、ｙ軸に予測値をプロットして、相関関係（ｙ＝ａｘ＋ｂ）と決定係数Ｒ２を求めることにより、予測精度が容易に判断される。
このようなニューラルネットワークの機械学習の結果、診断すべき被測定土壌を測定したデータ、つまり評価用データが処理済ＩＣＰ発光スペクトルデータ４１として入力されると、予測部３３は、中間データそして出力データ４５を二段階で重み係数を演算することにより、出力データ４５を出力する。即ち、出力データ４５としては、評価用データの各分析項目を出力する。このようなニューラルネットワークの学習結果に基づく各分析項目の出力データ４５は、予測値又は推定値とも呼ばれている。
従って、予測部３４の出力層４４からの出力データ４５から、ｘ軸に実測値、ｙ軸に予測値をプロットして、相関関係（ｙ＝ａｘ＋ｂ）と決定係数Ｒ２を求めることにより、予測精度が容易に判断される。
このように、学習のための既知のデータとして、処理済ＩＣＰ発光スペクトルデータ４１を入力層４２に入力すると共に、学習用データとして、出力データ４５が得られるように、予測部４３は重み係数を演算し決定する。
土壌分析装置１においては、例えば、土壌の１Ｍ酢酸アンモニウム溶液（ＮＨ４ＯＡｃ）による抽出液の処理済ＩＣＰ発光スペクトルデータについて、一例として、完全無作為抽出法により、８０％を学習用データとして選択し、残り２０％を評価用データとして使用してもよい。
本発明の土壌分析装置１によれば、一般的な試薬である１Ｍ酢酸アンモニウム溶液で抽出した土壌抽出液をＩＣＰ発光分光分析装置により分析し、得られた全波長スペクトルデータから、交換性塩基（Ｃａ、Ｋ、Ｍｇ、Ｎa）、交換性Ａｌ、有効態Ｐ、交換性塩基容量（ＣＥＣ）、ｐＨ（Ｈ２Ｏ抽出）、ｐＨ（ＫＣｌ抽出）、電気伝導度（ＥＣ）、全炭素、全窒素、粒径組成として粘土画分量及び砂画分量の全ての項目を同時に予測可能となる。
（深層学習後の新規の被測定土壌の診断方法）
図４は、本発明の土壌析装置１を用いて土壌診断方法を行う際に深層学習が済んだ後に新規の被測定土壌を診断する手順を示すフロー図である。図３において、第一のステップから第四のステップは、図２と同様に深層学習が行われ、第五のステップから第九のステップにおいて新規の被測定土壌の診断が行われる。
第五のステップ（ＳＴ５）において、深層学習を終了してよいか否かの判定を行う。
第六のステップ（ＳＴ６）において、第五のステップにて深層学習を終了してよいと判定された場合には、深層学習に用いた学習用データをデータセットとして格納する。
第七のステップ（ＳＴ７）において、新規の被測定土壌を所定の抽出液を用いて新規の土壌抽出液とする。
第八のステップ（ＳＴ８）において、プラズマ発光分光分析により新規の土壌抽出液の全波長スペクトルデータを取得する。
第九のステップ（ＳＴ９）において、第八のステップで取得した新規の土壌抽出液による新たな全波長スペクトルデータと第五のステップで取得した深層学習に用いた学習用データとを用いて、新規の被測定土壌の所定の分析項目の予測値を算出する。各ステップの詳細を次に説明する。
図４は、本発明の土壌析装置１を用いて土壌診断方法を行う際に深層学習が済んだ後に新規の被測定土壌を診断する手順を示すフロー図である。図３において、第一のステップから第四のステップは、図２と同様に深層学習が行われ、第五のステップから第九のステップにおいて新規の被測定土壌の診断が行われる。
第五のステップ（ＳＴ５）において、深層学習を終了してよいか否かの判定を行う。
第六のステップ（ＳＴ６）において、第五のステップにて深層学習を終了してよいと判定された場合には、深層学習に用いた学習用データをデータセットとして格納する。
第七のステップ（ＳＴ７）において、新規の被測定土壌を所定の抽出液を用いて新規の土壌抽出液とする。
第八のステップ（ＳＴ８）において、プラズマ発光分光分析により新規の土壌抽出液の全波長スペクトルデータを取得する。
第九のステップ（ＳＴ９）において、第八のステップで取得した新規の土壌抽出液による新たな全波長スペクトルデータと第五のステップで取得した深層学習に用いた学習用データとを用いて、新規の被測定土壌の所定の分析項目の予測値を算出する。各ステップの詳細を次に説明する。
（第五のステップ）
第五のステップ（ＳＴ５）において、深層学習を終了してよいか否かの判定を行う。この第五のステップ（ＳＴ５）では、上述した実測値と予測値の相関関係（ｙ＝ａｘ＋ｂ）と決定係数Ｒ２が所定の値以内であることを判定基準としてもよい。第五のステップにおいて、深層学習を終了してよい（ＹＥＳ）と判定された場合には第六のステップに進む。一方、深層学習が十分ではない（ＮＯ）場合には第四のステップに戻る。
第五のステップ（ＳＴ５）において、深層学習を終了してよいか否かの判定を行う。この第五のステップ（ＳＴ５）では、上述した実測値と予測値の相関関係（ｙ＝ａｘ＋ｂ）と決定係数Ｒ２が所定の値以内であることを判定基準としてもよい。第五のステップにおいて、深層学習を終了してよい（ＹＥＳ）と判定された場合には第六のステップに進む。一方、深層学習が十分ではない（ＮＯ）場合には第四のステップに戻る。
（第六のステップ）
第五のステップにおいて、深層学習を終了してよい（ＹＥＳ）と判定された場合には、第六のステップ（ＳＴ６）にて、上記の深層学習に用いた学習用データを記憶部３３にデータセットとして格納する。データセットは、学習用データに用いた各被測定土壌の採取情報と、各被測定土壌の土壌抽出液により取得したプラズマ発光分光分析により全波長スペクトルデータ、Ｃａ、Ｍｇ、Ｎａ、Ｋからなる交換性塩基の実測値、交換性塩基以外の分析項目の実測値、交換性塩基以外の分析項目の予測値等を含むように構成される。
第五のステップにおいて、深層学習を終了してよい（ＹＥＳ）と判定された場合には、第六のステップ（ＳＴ６）にて、上記の深層学習に用いた学習用データを記憶部３３にデータセットとして格納する。データセットは、学習用データに用いた各被測定土壌の採取情報と、各被測定土壌の土壌抽出液により取得したプラズマ発光分光分析により全波長スペクトルデータ、Ｃａ、Ｍｇ、Ｎａ、Ｋからなる交換性塩基の実測値、交換性塩基以外の分析項目の実測値、交換性塩基以外の分析項目の予測値等を含むように構成される。
（第七のステップ）
第七のステップ（ＳＴ７）において、新規の被測定土壌を所定の抽出液を用いて新規の土壌抽出液とする。第七のステップは、新規の被測定土壌を使用する以外は、第一のステップと同様に実施できる。
第七のステップ（ＳＴ７）において、新規の被測定土壌を所定の抽出液を用いて新規の土壌抽出液とする。第七のステップは、新規の被測定土壌を使用する以外は、第一のステップと同様に実施できる。
（第八のステップ）
第八のステップ（ＳＴ８）において、新規の土壌抽出液をＩＣＰ分析部１０でプラズマ発光分光分析により新たに全波長スペクトルデータを取得する。
第八のステップ（ＳＴ８）において、新規の土壌抽出液をＩＣＰ分析部１０でプラズマ発光分光分析により新たに全波長スペクトルデータを取得する。
（第九のステップ）
第九のステップ（ＳＴ９）において、第八のステップで取得した新たな全波長スペクトルデータと、第五のステップで取得した深層学習に用いた学習用データとから、診断すべき新規の被測定土壌の所定の分析項目の予測値を算出する。予測値を算出した後、被測定土壌の採取情報と、全波長スペクトルデータと被測定土壌の所定の分析項目の予測値とを、記憶部３３に新規の被測定土壌による新規な学習用データ、つまり、データセットとして格納する。新規の被測定土壌のデータセットを、データ処理部３０の図示しないディスプレイに表示してもよい。
第九のステップ（ＳＴ９）において、第八のステップで取得した新たな全波長スペクトルデータと、第五のステップで取得した深層学習に用いた学習用データとから、診断すべき新規の被測定土壌の所定の分析項目の予測値を算出する。予測値を算出した後、被測定土壌の採取情報と、全波長スペクトルデータと被測定土壌の所定の分析項目の予測値とを、記憶部３３に新規の被測定土壌による新規な学習用データ、つまり、データセットとして格納する。新規の被測定土壌のデータセットを、データ処理部３０の図示しないディスプレイに表示してもよい。
なお、図４に示す被測定土壌の診断方法を機能させるプログラムは、例えばＣ＋＋（登録商標）言語やＰｙｔｈｏｎ（登録商標）言語等のプログラム言語で記述され、当該プログラムを実行するデータ処理部３０の記憶装置３２に格納されている。当該プログラムは、記憶媒体に格納してから記憶装置３２に格納されてもよい。記憶媒体としては、ＣＤ－ＲＯＭ、ＤＶＤ－ＲＯＭやＵＳＢメモリ等を用いることができる。また、プログラムは、外部のサーバーからネットワークを介してＩＣＰ発光スペクトル処理部２０にダウンロードされてから記憶装置３２に格納されてもよい。
深層学習後の新規の被測定土壌の診断方法によれば、十分な深層学習を行った後で、新規の被測定土壌について、所定の分析項目の予測を精度良く迅速に行うことができる。
次に、土壌抽出液から土壌地力形質の各分析について、実施例により詳細に説明する。
次に、土壌抽出液から土壌地力形質の各分析について、実施例により詳細に説明する。
（土壌の詳細）
分析した土壌の詳細について説明する。
国際農林水産業研究センターが２０１６年から２０２０年の間に収集した熱帯・亜熱帯環境における土壌試料のうち、ＩＣＰによる交換性塩基の定量を実施した合計１９４２試料を対象とした。
国別の内訳はブルキナファソ（１０７４試料)、ラオス（３３４試料）、日本（２９８試料）、モザンビーク（１１５試料）、パラオ（５９試料）、マダガスカル（４９試料）、フィリピン（１３試料）である。交換性塩基と他の分析項目の国別の試料数を表１に示す。
分析した土壌の詳細について説明する。
国際農林水産業研究センターが２０１６年から２０２０年の間に収集した熱帯・亜熱帯環境における土壌試料のうち、ＩＣＰによる交換性塩基の定量を実施した合計１９４２試料を対象とした。
国別の内訳はブルキナファソ（１０７４試料)、ラオス（３３４試料）、日本（２９８試料）、モザンビーク（１１５試料）、パラオ（５９試料）、マダガスカル（４９試料）、フィリピン（１３試料）である。交換性塩基と他の分析項目の国別の試料数を表１に示す。
各国の農耕地および林野土壌は多様な目的で採取しているため、表層試料に加えて深度別試料を含み、施肥履歴も多用な試料セットとなっている。当該試料について、完全無作為抽出法により８０％（１５５３件）をステップ１０の学習用データとして選択し、残り２０％（３８８件）をステップ２０の評価用データとして使用した。
（土壌の処理方法）
採取した土壌を室温環境にて約１週間の風乾後、目開き２ｍｍで篩別し、植物残渣を取り除いたものを土壌試料とした。
全波長スペクトルは、ｐＨが７の１Ｍ酢酸アンモニウム溶液を用いて抽出した土壌抽出液により取得した。土壌抽出液の調製及び測定は以下の通り実施した。
（１）土壌試料２．５ｇ（±０．０１ｇ）を５０ｍＬの蓋つき試験管に秤取し、１Ｍの酢酸アンモニウム溶液１２．５ｍＬを加えて振とう器（YAMATO社製、型番、SA300、以下振とう器と呼ぶ）を用いて３０分間振とうした。
（２）その後、３０００ｒｐｍで３から５分間遠心分離し、上清をフィルター（アドバンテック社製、以下、ADVANTEC No.５Cと呼ぶ）により、ろ過した。
次に、試験管に残った土壌に、再び１Ｍの酢酸アンモニウム１２．５ｍＬを加えて、振とう器で３０分間振とうした。その後３０００ｒｐｍで３から５分間遠心分離し、上清をフィルター（ADVANTEC No.5C）によりろ過した。これをもう一度繰り返すことにより合計３回の交換性塩基の抽出を行った。
最後に、１Ｍの酢酸アンモニウム１２．５ｍＬを加えて撹拌機（サイエンティフィックインダストリーズ社製、型番、VOLTEX、以下撹拌機と呼ぶ）で１０秒間振とうし、３０００ｒｐｍで３から５分間遠心分離し、上清をフィルター（ADVANTEC 5C 110mm）を使ってろ過した。計５０ｍＬ採れた抽出液を分析試料とした。
採取した土壌を室温環境にて約１週間の風乾後、目開き２ｍｍで篩別し、植物残渣を取り除いたものを土壌試料とした。
全波長スペクトルは、ｐＨが７の１Ｍ酢酸アンモニウム溶液を用いて抽出した土壌抽出液により取得した。土壌抽出液の調製及び測定は以下の通り実施した。
（１）土壌試料２．５ｇ（±０．０１ｇ）を５０ｍＬの蓋つき試験管に秤取し、１Ｍの酢酸アンモニウム溶液１２．５ｍＬを加えて振とう器（YAMATO社製、型番、SA300、以下振とう器と呼ぶ）を用いて３０分間振とうした。
（２）その後、３０００ｒｐｍで３から５分間遠心分離し、上清をフィルター（アドバンテック社製、以下、ADVANTEC No.５Cと呼ぶ）により、ろ過した。
次に、試験管に残った土壌に、再び１Ｍの酢酸アンモニウム１２．５ｍＬを加えて、振とう器で３０分間振とうした。その後３０００ｒｐｍで３から５分間遠心分離し、上清をフィルター（ADVANTEC No.5C）によりろ過した。これをもう一度繰り返すことにより合計３回の交換性塩基の抽出を行った。
最後に、１Ｍの酢酸アンモニウム１２．５ｍＬを加えて撹拌機（サイエンティフィックインダストリーズ社製、型番、VOLTEX、以下撹拌機と呼ぶ）で１０秒間振とうし、３０００ｒｐｍで３から５分間遠心分離し、上清をフィルター（ADVANTEC 5C 110mm）を使ってろ過した。計５０ｍＬ採れた抽出液を分析試料とした。
（土壌分析装置）
土壌分析装置１として、以下の構成のＩＣＰ分析部１０と土壌分析部２０とを用い、ソフトウェアとして、Neural Network Console version 1.7.7352.44102（Sony Network Communications Inc.）を用いたプログラムにより実行し、土壌分析を実施した。
ＩＣＰ分析部：島津製作所製（ＩＣＰ発光分光分析装置、ＩＣＰＥ９０００）
土壌分析部：
ワークステーション：アプライド社製（型番CERVO-Deep for windows）
ＣＰＵ: インテル社（登録商標）製（Xeon（登録商標）8-Core W-2245 3.90GHz）
ＧＰＵ:NVIDIA社(登録商標)製（GeForce（登録商標）（RTX 2080Ti（DRAM：11（GDDR6×２））
ＤＲＡＭメモリの容量：１５２ＧＢ
ＨＤＤ容量：２ＴＢ
土壌分析装置１として、以下の構成のＩＣＰ分析部１０と土壌分析部２０とを用い、ソフトウェアとして、Neural Network Console version 1.7.7352.44102（Sony Network Communications Inc.）を用いたプログラムにより実行し、土壌分析を実施した。
ＩＣＰ分析部：島津製作所製（ＩＣＰ発光分光分析装置、ＩＣＰＥ９０００）
土壌分析部：
ワークステーション：アプライド社製（型番CERVO-Deep for windows）
ＣＰＵ: インテル社（登録商標）製（Xeon（登録商標）8-Core W-2245 3.90GHz）
ＧＰＵ:NVIDIA社(登録商標)製（GeForce（登録商標）（RTX 2080Ti（DRAM：11（GDDR6×２））
ＤＲＡＭメモリの容量：１５２ＧＢ
ＨＤＤ容量：２ＴＢ
学習用データにおいて、実測した土壌分析値は、以下の方法により定量分析したものを使用した。
交換性塩基：
土壌をｐＨ７の１Ｍ酢酸アンモニウム溶液による抽出液を、マルチタイプＩＣＰ発光分光分析装置を用いて、全波長スペクトルを分析することにより、交換性塩基（Ｃａ,Ｍｇ,Ｋ,Ｎａ）を定量した。
ＣＥＣ：
ＣＥＣの実測値を以下の方法により定量した。交換性塩基抽出残渣(２.５g土壌)を蒸留水及び８０％メタノールで洗浄した後、１０％ＫＣｌを１２．５ｍＬ添加し、３０分間振とう機で振とうして３０００ｒｐｍで遠心し上澄みをろ過した。１０％ＫＣｌによる抽出処理を合計３回繰り返した後、さらに１０％ＫＣｌを１２．５ｍＬ添加して、撹拌機でよく攪拌した後、遠心分離したものを同様にろ過することで５０ｍＬの抽出液を得た。得られた抽出液のアンモニア態窒素量についてオートアナライザー３型(ビーエルテック社製（BL-tec社、AAIII)）を用いてサリチル酸法で定量してＣＥＣの実測値を定量した。
ｐＨ（Ｈ２Ｏ抽出）：ガラス電極法（固液比は１：５）により定量した。
ｐＨ（ＫＣｌ抽出）：ガラス電極法（固液比は１：５）により定量した。
電気伝導度（ＥＣ）：ガラス電極法（固液比は１：５）により定量した。
有効態リン：Ｂｒａｙ１法を用いた。モリブデンブルー発色後、紫外可視分光光度計により波長７１０ｎｍで定量した。
交換性アルミニウム：１Ｍ塩化カリウム溶液で抽出後、ＩＣＰで定量分析した。
全炭素：ＮＣアナライザーを用い、乾式燃焼法により定量した。
全窒素：ＮＣアナライザーを用い、乾式燃焼法により定量した。
粒径組成：粘土画分及び砂画分は、沈降法により測定した。過酸化水素水による有機物分解後、ヘキサメタリン酸ナトリウムを分散材として添加し、１６時間振とうした後、ピペット法により定量した。
交換性塩基：
土壌をｐＨ７の１Ｍ酢酸アンモニウム溶液による抽出液を、マルチタイプＩＣＰ発光分光分析装置を用いて、全波長スペクトルを分析することにより、交換性塩基（Ｃａ,Ｍｇ,Ｋ,Ｎａ）を定量した。
ＣＥＣ：
ＣＥＣの実測値を以下の方法により定量した。交換性塩基抽出残渣(２.５g土壌)を蒸留水及び８０％メタノールで洗浄した後、１０％ＫＣｌを１２．５ｍＬ添加し、３０分間振とう機で振とうして３０００ｒｐｍで遠心し上澄みをろ過した。１０％ＫＣｌによる抽出処理を合計３回繰り返した後、さらに１０％ＫＣｌを１２．５ｍＬ添加して、撹拌機でよく攪拌した後、遠心分離したものを同様にろ過することで５０ｍＬの抽出液を得た。得られた抽出液のアンモニア態窒素量についてオートアナライザー３型(ビーエルテック社製（BL-tec社、AAIII)）を用いてサリチル酸法で定量してＣＥＣの実測値を定量した。
ｐＨ（Ｈ２Ｏ抽出）：ガラス電極法（固液比は１：５）により定量した。
ｐＨ（ＫＣｌ抽出）：ガラス電極法（固液比は１：５）により定量した。
電気伝導度（ＥＣ）：ガラス電極法（固液比は１：５）により定量した。
有効態リン：Ｂｒａｙ１法を用いた。モリブデンブルー発色後、紫外可視分光光度計により波長７１０ｎｍで定量した。
交換性アルミニウム：１Ｍ塩化カリウム溶液で抽出後、ＩＣＰで定量分析した。
全炭素：ＮＣアナライザーを用い、乾式燃焼法により定量した。
全窒素：ＮＣアナライザーを用い、乾式燃焼法により定量した。
粒径組成：粘土画分及び砂画分は、沈降法により測定した。過酸化水素水による有機物分解後、ヘキサメタリン酸ナトリウムを分散材として添加し、１６時間振とうした後、ピペット法により定量した。
（学習用データ及び評価用データ）
上記のように完全無作為抽出法により全試料から８０％の試料を選択し学習用データとした。なお、残りの２０％を評価用データとした。
１Ｍ酢酸アンモニウム溶液で抽出した抽出液を希釈し、ＩＣＰ分析部１０により全波長スペクトルデータを取得し、Ｃａ、Ｋ、Ｍｇ、Ｎａを定量した。また、その際に得られる全波長スペクトルデータを取得し、上述したＩＣＰ発光スペクトル処理部２１で処理された処理済ＩＣＰ発光スペクトルデータを取得し、学習用データを取得した。学習の後で診断すべき被測定土壌の所定の分析項目の予測値を算出した。具体的には、全試料の８０％の試料により学習用データを取得した後で、全試料の２０％の評価用データにより土壌分析を実施した。プラズマ発光分光分析による土壌抽出液の全波長スペクトルデータを評価用データとし、この全波長スペクトルデータと学習用データとから予測値を取得した。
上記のように完全無作為抽出法により全試料から８０％の試料を選択し学習用データとした。なお、残りの２０％を評価用データとした。
１Ｍ酢酸アンモニウム溶液で抽出した抽出液を希釈し、ＩＣＰ分析部１０により全波長スペクトルデータを取得し、Ｃａ、Ｋ、Ｍｇ、Ｎａを定量した。また、その際に得られる全波長スペクトルデータを取得し、上述したＩＣＰ発光スペクトル処理部２１で処理された処理済ＩＣＰ発光スペクトルデータを取得し、学習用データを取得した。学習の後で診断すべき被測定土壌の所定の分析項目の予測値を算出した。具体的には、全試料の８０％の試料により学習用データを取得した後で、全試料の２０％の評価用データにより土壌分析を実施した。プラズマ発光分光分析による土壌抽出液の全波長スペクトルデータを評価用データとし、この全波長スペクトルデータと学習用データとから予測値を取得した。
（交換性塩基（Ｃａ、Ｍｇ、Ｋ、Ｎａ）の予測結果）
図５は、土壌の交換性塩基の予測値と実測値との関係図であり、（ａ）はＣａ、（ｂ）はＭｇ、（ｃ）はＫ、（ｄ）はＮａを示す。
図５（ａ）のＣａでは、相関関係はｙ＝０．９６７ｘ＋０．２３１であり、決定係数Ｒ２が０．９９５であった。図５（ｂ）のＭｇでは、相関関係はｙ＝１．０２８ｘ-０．０６０であり、決定係数Ｒ２が０．９８３であった。図５（ｃ）のＫでは、相関関係はｙ＝０．９８９ｘ-０．０８であり、決定係数Ｒ２が０．９９３であった。図５（ｄ）のＮａでは、相関関係はｙ＝１．０５１ｘ+０．０４８であり、決定係数Ｒ２が０．９４６であった。図５に示す交換性塩基（Ｃａ、Ｍｇ、Ｎａ、Ｋ）の予測については、ＩＣＰ分析部１０を用いた公定法による分析であり、ほぼ完全に実測値と予測値が一致することを確認した。
図５は、土壌の交換性塩基の予測値と実測値との関係図であり、（ａ）はＣａ、（ｂ）はＭｇ、（ｃ）はＫ、（ｄ）はＮａを示す。
図５（ａ）のＣａでは、相関関係はｙ＝０．９６７ｘ＋０．２３１であり、決定係数Ｒ２が０．９９５であった。図５（ｂ）のＭｇでは、相関関係はｙ＝１．０２８ｘ-０．０６０であり、決定係数Ｒ２が０．９８３であった。図５（ｃ）のＫでは、相関関係はｙ＝０．９８９ｘ-０．０８であり、決定係数Ｒ２が０．９９３であった。図５（ｄ）のＮａでは、相関関係はｙ＝１．０５１ｘ+０．０４８であり、決定係数Ｒ２が０．９４６であった。図５に示す交換性塩基（Ｃａ、Ｍｇ、Ｎａ、Ｋ）の予測については、ＩＣＰ分析部１０を用いた公定法による分析であり、ほぼ完全に実測値と予測値が一致することを確認した。
（交換性Ａｌの予測結果）
表１に示すように、交換性Ａｌの実測は合計１０４２試料を対象とし、国別の内訳はブルキナファソ（５１６試料)、ラオス（３３４試料）、日本（０試料）、モザンビーク（７１試料）、パラオ（５９試料）、マダガスカル（４９試料）、フィリピン（１３試料）である。
図６は土壌の交換性Ａｌの予測値と実測値との関係を示す図である。交換性Ａｌでは、相関関係はｙ＝０．９６４ｘ+０．０２４であり、決定係数Ｒ２が０．９６４であった。図６に示す交換性Ａｌについては、公定法では１ＭのＫＣｌによる抽出後、ＩＣＰで測定するが、抽出液として公定法で用いた１ＭＫＣｌとは異なる１Ｍ酢酸アンモニウム溶液とした本実施例においても、高精度の予測が可能であることが判明した。
表１に示すように、交換性Ａｌの実測は合計１０４２試料を対象とし、国別の内訳はブルキナファソ（５１６試料)、ラオス（３３４試料）、日本（０試料）、モザンビーク（７１試料）、パラオ（５９試料）、マダガスカル（４９試料）、フィリピン（１３試料）である。
図６は土壌の交換性Ａｌの予測値と実測値との関係を示す図である。交換性Ａｌでは、相関関係はｙ＝０．９６４ｘ+０．０２４であり、決定係数Ｒ２が０．９６４であった。図６に示す交換性Ａｌについては、公定法では１ＭのＫＣｌによる抽出後、ＩＣＰで測定するが、抽出液として公定法で用いた１ＭＫＣｌとは異なる１Ｍ酢酸アンモニウム溶液とした本実施例においても、高精度の予測が可能であることが判明した。
（有効態Ｐの予測結果）
表１に示すように、有効態Ｐの実測は合計１９４２試料を対象とし、国別の内訳はブルキナファソ（１０７４試料)、ラオス（３３４試料）、日本（２９８試料）、モザンビーク（１１５試料）、パラオ（５９試料）、マダガスカル（４９試料）、フィリピン（１３試料）である。
図７は有効態Ｐの予測値と実測値との関係を示す図である。有効態Ｐｌでは、相関関係はｙ＝０．９８５ｘ+０．８３５であり、決定係数Ｒ２が０．９６４であった。
実施例の有効態Ｐの予測方法は、従来のＢｒａｙ法を用いた有効態Ｐの公定法とは、抽出液及び測定方法も異なるが、本実施例の予測結果により高精度の予測が可能であることが判明した。
このように、１Ｍ酢酸アンモニウム溶液により抽出される交換性Ａｌ及び有効態Ｐの予測値は、従来の公定法と一定の相関があったものと考えられる。
表１に示すように、有効態Ｐの実測は合計１９４２試料を対象とし、国別の内訳はブルキナファソ（１０７４試料)、ラオス（３３４試料）、日本（２９８試料）、モザンビーク（１１５試料）、パラオ（５９試料）、マダガスカル（４９試料）、フィリピン（１３試料）である。
図７は有効態Ｐの予測値と実測値との関係を示す図である。有効態Ｐｌでは、相関関係はｙ＝０．９８５ｘ+０．８３５であり、決定係数Ｒ２が０．９６４であった。
実施例の有効態Ｐの予測方法は、従来のＢｒａｙ法を用いた有効態Ｐの公定法とは、抽出液及び測定方法も異なるが、本実施例の予測結果により高精度の予測が可能であることが判明した。
このように、１Ｍ酢酸アンモニウム溶液により抽出される交換性Ａｌ及び有効態Ｐの予測値は、従来の公定法と一定の相関があったものと考えられる。
（ｐＨの予測結果)
表１に示すように、ｐＨの実測はＨ２Ｏ抽出及びＫＣｌ抽出のそれぞれにおいて合計１９４１試料を対象とし、国別の内訳はブルキナファソ（１０７３試料)、ラオス（３３４試料）、日本（２９８試料）、モザンビーク（１１５試料）、パラオ（５９試料）、マダガスカル（４９試料）、フィリピン（１３試料）である。
図８は、土壌ｐＨの予測値と実測値との関係図であり、（ａ）はＨ２Ｏ抽出、（ｂ）はＫＣｌ抽出の場合を示す。
図８（ａ）のＨ２Ｏ抽出における相関関係はｙ＝０．９３６ｘ+０．３７２であり、決定係数Ｒ２が０．９５１であった。図８（ｂ）のＫＣｌ抽出における相関関係はｙ＝０．９９６６ｘ+０．１４１であり、決定係数Ｒ２が０．９５９であった。
表１に示すように、ｐＨの実測はＨ２Ｏ抽出及びＫＣｌ抽出のそれぞれにおいて合計１９４１試料を対象とし、国別の内訳はブルキナファソ（１０７３試料)、ラオス（３３４試料）、日本（２９８試料）、モザンビーク（１１５試料）、パラオ（５９試料）、マダガスカル（４９試料）、フィリピン（１３試料）である。
図８は、土壌ｐＨの予測値と実測値との関係図であり、（ａ）はＨ２Ｏ抽出、（ｂ）はＫＣｌ抽出の場合を示す。
図８（ａ）のＨ２Ｏ抽出における相関関係はｙ＝０．９３６ｘ+０．３７２であり、決定係数Ｒ２が０．９５１であった。図８（ｂ）のＫＣｌ抽出における相関関係はｙ＝０．９９６６ｘ+０．１４１であり、決定係数Ｒ２が０．９５９であった。
土壌ｐＨは、従来の公定法ではガラス電極法により測定するため、処理済ＩＣＰ発光スペクトルデータ４１とは直接的な関係性は無いものと考えられる。しかしながら、上記実施例においては、土壌診断に向けた予測値としては十分な精度で予測が可能であることが判明した。これは、土壌ｐＨは、処理済ＩＣＰ発光スペクトルデータ４１によって検出された１Ｍ酢酸アンモニウム溶液の抽出液のイオンのバランスによって規定されるので、ｐＨとの相関があることにより十分な精度で予測できたと推定される。
（電気伝導度の予測結果)
表１に示すように、電気伝導度（ＥＣ）の実測は合計１９４１試料を対象とし、国別の内訳はブルキナファソ（１０７３試料)、ラオス（３３４試料）、日本（２９８試料）、モザンビーク（１１５試料）、パラオ（５９試料）、マダガスカル（４９試料）、フィリピン（１３試料）である。
図９は、電気伝導度（ｍＳ／ｍ）の予測値と実測値との関係図である。電気伝導度の相関関係はｙ＝０．９８９ｘ+０．７０７であり、決定係数Ｒ２が０．９６０であった。
電気伝導度は、従来の公定法では土壌ｐＨと同様にガラス電極法により測定するため、処理済ＩＣＰ発光スペクトルデータ４１とは直接的な関係性は無いものと考えられる。しかしながら、上記実施例においては、土壌診断に向けた電気伝導度の予測値としては十分な精度で予測が可能であることが判明した。電気伝導度は、処理済ＩＣＰ発光スペクトルデータ４１によって検出された１Ｍ酢酸アンモニウム溶液の抽出液中の溶出イオン総量と相関があることにより十分な精度で予測できたと推定される。
表１に示すように、電気伝導度（ＥＣ）の実測は合計１９４１試料を対象とし、国別の内訳はブルキナファソ（１０７３試料)、ラオス（３３４試料）、日本（２９８試料）、モザンビーク（１１５試料）、パラオ（５９試料）、マダガスカル（４９試料）、フィリピン（１３試料）である。
図９は、電気伝導度（ｍＳ／ｍ）の予測値と実測値との関係図である。電気伝導度の相関関係はｙ＝０．９８９ｘ+０．７０７であり、決定係数Ｒ２が０．９６０であった。
電気伝導度は、従来の公定法では土壌ｐＨと同様にガラス電極法により測定するため、処理済ＩＣＰ発光スペクトルデータ４１とは直接的な関係性は無いものと考えられる。しかしながら、上記実施例においては、土壌診断に向けた電気伝導度の予測値としては十分な精度で予測が可能であることが判明した。電気伝導度は、処理済ＩＣＰ発光スペクトルデータ４１によって検出された１Ｍ酢酸アンモニウム溶液の抽出液中の溶出イオン総量と相関があることにより十分な精度で予測できたと推定される。
（ＣＥＣ（交換性塩基容量）の予測結果）
表１に示すように、ＣＥＣの実測は合計１９４１試料を対象とし、国別の内訳はブルキナファソ（１０７４試料)、ラオス（３３４試料）、日本（２９８試料）、モザンビーク（１１５試料）、パラオ（５９試料）、マダガスカル（４９試料）、フィリピン（１３試料）である。
図１０は、ＣＥＣの予測値と実測値との関係図である。ＣＥＣの相関関係はｙ＝０．９２４ｘ+０．９８４であり、決定係数Ｒ２が０．９１１であった。これにより、上記実施例においては、土壌診断に向けたＣＥＣの予測値としては十分な精度で予測が可能であることが判明した。これは、ＣＥＣが図５に示した交換性塩基（Ｃａ、Ｍｇ、Ｋ、Ｎa）と図６に示した交換性Ａｌの和によって、一定の予測が出来ることが知られており、交換性塩基の予測値と交換性Ａｌの予測値との相互関係から精度良く予測されていると推定される。
表１に示すように、ＣＥＣの実測は合計１９４１試料を対象とし、国別の内訳はブルキナファソ（１０７４試料)、ラオス（３３４試料）、日本（２９８試料）、モザンビーク（１１５試料）、パラオ（５９試料）、マダガスカル（４９試料）、フィリピン（１３試料）である。
図１０は、ＣＥＣの予測値と実測値との関係図である。ＣＥＣの相関関係はｙ＝０．９２４ｘ+０．９８４であり、決定係数Ｒ２が０．９１１であった。これにより、上記実施例においては、土壌診断に向けたＣＥＣの予測値としては十分な精度で予測が可能であることが判明した。これは、ＣＥＣが図５に示した交換性塩基（Ｃａ、Ｍｇ、Ｋ、Ｎa）と図６に示した交換性Ａｌの和によって、一定の予測が出来ることが知られており、交換性塩基の予測値と交換性Ａｌの予測値との相互関係から精度良く予測されていると推定される。
（全炭素（ＴＣ）の予測結果)
表１に示すように、全炭素の実測は合計１８１５試料を対象とし、国別の内訳はブルキナファソ（９４７試料)、ラオス（３３４試料）、日本（２９８試料）、モザンビーク（１１５試料）、パラオ（５９試料）、マダガスカル（４９試料）、フィリピン（１３試料）である。
図１１は全炭素（ＴＣ）の予測値と実測値との関係図である。全炭素の相関関係はｙ＝０．８６６ｘ+０．１６４であり、決定係数Ｒ２が０．８１２であった。
実施例の抽出液であるＮＨ４ＯＡｃは、炭素を多く含む１Ｍ酢酸アンモニウム溶液であり、関連するＩＣＰ発光スペクトルをマスクして予測の妨げとなることが考えられたが、図１１に示すように全炭素についても一定程度の予測が可能であることが判明した。
表１に示すように、全炭素の実測は合計１８１５試料を対象とし、国別の内訳はブルキナファソ（９４７試料)、ラオス（３３４試料）、日本（２９８試料）、モザンビーク（１１５試料）、パラオ（５９試料）、マダガスカル（４９試料）、フィリピン（１３試料）である。
図１１は全炭素（ＴＣ）の予測値と実測値との関係図である。全炭素の相関関係はｙ＝０．８６６ｘ+０．１６４であり、決定係数Ｒ２が０．８１２であった。
実施例の抽出液であるＮＨ４ＯＡｃは、炭素を多く含む１Ｍ酢酸アンモニウム溶液であり、関連するＩＣＰ発光スペクトルをマスクして予測の妨げとなることが考えられたが、図１１に示すように全炭素についても一定程度の予測が可能であることが判明した。
（全窒素（ＴＮ）の予測結果)
表１に示すように、全窒素の実測は合計１８１５試料を対象とし、国別の内訳はブルキナファソ（９４７試料)、ラオス（３３４試料）、日本（２９８試料）、モザンビーク（１１５試料）、パラオ（５９試料）、マダガスカル（４９試料）、フィリピン（１３試料）である。
図１２は全窒素（ＴＮ）の予測値と実測値との関係図である。全窒素の相関関係はｙ＝０．９０５ｘ+０．０１０であり、決定係数Ｒ２が０．９５６であった。実施例の抽出液であるＮＨ４ＯＡｃは、窒素を多く含む１Ｍ酢酸アンモニウム溶液であり、関連するＩＣＰ発光スペクトルをマスクして予測の妨げとなることが考えられたが、図１２に示すように全窒素についても一定程度の予測が可能であることが判明した。
表１に示すように、全窒素の実測は合計１８１５試料を対象とし、国別の内訳はブルキナファソ（９４７試料)、ラオス（３３４試料）、日本（２９８試料）、モザンビーク（１１５試料）、パラオ（５９試料）、マダガスカル（４９試料）、フィリピン（１３試料）である。
図１２は全窒素（ＴＮ）の予測値と実測値との関係図である。全窒素の相関関係はｙ＝０．９０５ｘ+０．０１０であり、決定係数Ｒ２が０．９５６であった。実施例の抽出液であるＮＨ４ＯＡｃは、窒素を多く含む１Ｍ酢酸アンモニウム溶液であり、関連するＩＣＰ発光スペクトルをマスクして予測の妨げとなることが考えられたが、図１２に示すように全窒素についても一定程度の予測が可能であることが判明した。
全炭素及び全窒素は、土壌が含有する有機物量を反映しており、従来の公定法では乾式燃焼法によって定量され、ＩＣＰ発光スペクトルとの関係性は低いものであった。しかしながら、上記実施例で説明したように全炭素及び全窒素が所定の精度で予測できた理由として、ＩＣＰ発光スペクトルにおいて炭素に特徴的な波長が発光するとの報告があり、当該波長と有機物の安定化に寄与する種々のイオン濃度から予測された可能性があると推定される。
（粒径組成の予測結果)
表１に示すように、粒径組成の内砂画分量の実測は合計８５８試料を対象とし、国別の内訳はブルキナファソ（３４８試料)、ラオス（３３４試料）、日本（４８試料）、モザンビーク（１１５試料）、パラオ（０試料）、マダガスカル（０試料）、フィリピン（１３試料）である。粒径組成の内粘土画分量の実測は合計８５２試料を対象とし、国別の内訳はブルキナファソ（３４２試料)とした以外の他の国の試料数は、砂画分量と同じである。
図１３は、土壌の粒径組成の予測値と実測値との関係図であり、（ａ）は砂画分量（％）、（ｂ）は粘土画分量（％）を示す。砂画分量の相関関係はｙ＝０．９６２ｘ+１．０４３であり、決定係数Ｒ２が０．８７１であった。図１３（ｂ）に示すように、粘土画分量の相関関係はｙ＝０．８６１ｘ+４．１６６であり、決定係数Ｒ２が０．８４６であった。
図１３に示すように土壌の粒径組成として、砂画分量及び粘土画分量についても一定程度の予測が可能であることが判明した。
表１に示すように、粒径組成の内砂画分量の実測は合計８５８試料を対象とし、国別の内訳はブルキナファソ（３４８試料)、ラオス（３３４試料）、日本（４８試料）、モザンビーク（１１５試料）、パラオ（０試料）、マダガスカル（０試料）、フィリピン（１３試料）である。粒径組成の内粘土画分量の実測は合計８５２試料を対象とし、国別の内訳はブルキナファソ（３４２試料)とした以外の他の国の試料数は、砂画分量と同じである。
図１３は、土壌の粒径組成の予測値と実測値との関係図であり、（ａ）は砂画分量（％）、（ｂ）は粘土画分量（％）を示す。砂画分量の相関関係はｙ＝０．９６２ｘ+１．０４３であり、決定係数Ｒ２が０．８７１であった。図１３（ｂ）に示すように、粘土画分量の相関関係はｙ＝０．８６１ｘ+４．１６６であり、決定係数Ｒ２が０．８４６であった。
図１３に示すように土壌の粒径組成として、砂画分量及び粘土画分量についても一定程度の予測が可能であることが判明した。
土壌が含有する砂画分量及び粘土画分量は、土壌の物理的性質を示し、これらの従来の定量分析は、粒子径によって水中沈降速度が変化することを利用した物理的手法によって定量されており、ＩＣＰ発光スペクトルとの関係性は低いものであった。しかしながら、上記実施例で説明したように砂画分量及び粘土画分量が所定の精度で予測できた理由として、砂画分量および粘土画分量の差異により、土壌中のＦｅやＡｌ等の１Ｍ酢酸アンモニウム溶液への溶出特性が変化し、これらのＩＣＰ発光スペクトルから予測がされたものと推測される。
本発明はその趣旨を逸脱しない範囲において様々な形態で実施することができる。
例えば、上述した実施形態においては、対象としなかった、可溶性Ｆｅ（鉄）、交換性Ｍｎ（マンガン）、可溶性Ｚｎ（亜鉛）、リン酸吸収係数、硝酸態窒素量、アンモニア態窒素量などの分析項目についても、Ｆｅ、Ｍｎ、ＺｎはＩＣＰ発光スペクトルで分析できる元素であり、リン酸吸収係数、硝酸態窒素量、アンモニア態窒素量について土壌ｐＨや交換性Ａｌ、全窒素等の関連の分析項目が高精度で予測されているので、深層学習で実測値と共に学習をすれば、同様に予測が可能であることが強く推定される。
例えば、上述した実施形態においては、対象としなかった、可溶性Ｆｅ（鉄）、交換性Ｍｎ（マンガン）、可溶性Ｚｎ（亜鉛）、リン酸吸収係数、硝酸態窒素量、アンモニア態窒素量などの分析項目についても、Ｆｅ、Ｍｎ、ＺｎはＩＣＰ発光スペクトルで分析できる元素であり、リン酸吸収係数、硝酸態窒素量、アンモニア態窒素量について土壌ｐＨや交換性Ａｌ、全窒素等の関連の分析項目が高精度で予測されているので、深層学習で実測値と共に学習をすれば、同様に予測が可能であることが強く推定される。
例えば、Ｆｅの予測は、以下に示す公定法によりＦｅの実測値を取得し、取得した全波長スペクトルデータとＦｅの実測値とにより学習データの作成を行い、学習をした後で評価用データの分析を行うことにより、交換性塩基（Ｃａ、Ｍｇ、Ｎａ、Ｋ）と同様にＦｅの予測が可能となる。
可溶性Ｆｅ、交換性Ｍｎ、可溶性Ｚｎ、リン酸吸収係数、硝酸態窒素量、アンモニア態窒素量などの実測値の取得方法を以下に示す。
可溶性Ｆｅ：
土壌試料１．００ｇに対して、ｐＨ３．０の０．２Ｍ酸性シュウ酸溶液１００ｍＬを添加し、室温（約２０℃）且つ暗条件で４時間振とうする。この試料液約２５ｍＬを容積５０ｍＬの遠沈管に取り、高分子凝集剤（例えばアコフロック（登録商標），ＭＴアクアポリマー株式会社製）を１滴加えて振り混ぜ、２０００から５０００ｒｐｍで１０分間遠心分離する。得られた上澄み液のＦｅ濃度をＩＣＰ発光分析装置で定量する。
交換性Ｍｎ：
土壌試料５ｇに対して５０ｍＬのｐＨ７で１ＭのＮＨ４ＯＡｃ溶液を添加し、１時間振とうした後、放置または遠心分離して得られた上澄み液をろ過する。得られたろ液中のＭｎ濃度をＩＣＰ発光分光分析装置により定量する。
可溶性Ｚｎ：
土壌試料１０ｇに対して５０ｍＬの０．１ＭＨＣｌ溶液を添加し、恒温水平振り混ぜ機で３０℃に保ち、１時間振り混ぜた後、乾燥ろ紙でろ過し、得られたろ液中のＺｎ濃度をＩＣＰ発光分光分析装置で定量する。
リン酸吸収係数：
一定濃度のリン酸溶液を土壌に添加し、２４時間後の抽出溶液中のリン酸量を定量することで、土壌に吸着されたリン酸量を測定する。例えば、１Ｌ（リットル）当たり１３．４４ｇのＰ２Ｏ５（１３．４４ｇ／Ｌ）のリン酸アンモニウム溶液を土壌に添加した後、適宜、振り混ぜながら２４時間放置し、乾燥ろ紙を用いてろ過し、得られた抽出液のリン酸量をバナドモリブデン酸法などのリン酸定量法で定量することで、土壌に吸着したリン酸量を計算する。
硝酸態窒素量：
土壌試料に２ＭＫＣｌを添加し、１時間振とうした後、上澄み液をろ過し、抽出液中の硝酸態窒素量を、銅・カドミウム還元－ナフチルエチレンジアミン吸光光度法によって定量する。
アンモニア態窒素量：
土壌試料に２ＭＫＣｌを添加し、１時間振とうした後、上澄み液をろ過し、抽出液中の硝酸態窒素量を、サリチル酸吸光光度法で定量する。
可溶性Ｆｅ：
土壌試料１．００ｇに対して、ｐＨ３．０の０．２Ｍ酸性シュウ酸溶液１００ｍＬを添加し、室温（約２０℃）且つ暗条件で４時間振とうする。この試料液約２５ｍＬを容積５０ｍＬの遠沈管に取り、高分子凝集剤（例えばアコフロック（登録商標），ＭＴアクアポリマー株式会社製）を１滴加えて振り混ぜ、２０００から５０００ｒｐｍで１０分間遠心分離する。得られた上澄み液のＦｅ濃度をＩＣＰ発光分析装置で定量する。
交換性Ｍｎ：
土壌試料５ｇに対して５０ｍＬのｐＨ７で１ＭのＮＨ４ＯＡｃ溶液を添加し、１時間振とうした後、放置または遠心分離して得られた上澄み液をろ過する。得られたろ液中のＭｎ濃度をＩＣＰ発光分光分析装置により定量する。
可溶性Ｚｎ：
土壌試料１０ｇに対して５０ｍＬの０．１ＭＨＣｌ溶液を添加し、恒温水平振り混ぜ機で３０℃に保ち、１時間振り混ぜた後、乾燥ろ紙でろ過し、得られたろ液中のＺｎ濃度をＩＣＰ発光分光分析装置で定量する。
リン酸吸収係数：
一定濃度のリン酸溶液を土壌に添加し、２４時間後の抽出溶液中のリン酸量を定量することで、土壌に吸着されたリン酸量を測定する。例えば、１Ｌ（リットル）当たり１３．４４ｇのＰ２Ｏ５（１３．４４ｇ／Ｌ）のリン酸アンモニウム溶液を土壌に添加した後、適宜、振り混ぜながら２４時間放置し、乾燥ろ紙を用いてろ過し、得られた抽出液のリン酸量をバナドモリブデン酸法などのリン酸定量法で定量することで、土壌に吸着したリン酸量を計算する。
硝酸態窒素量：
土壌試料に２ＭＫＣｌを添加し、１時間振とうした後、上澄み液をろ過し、抽出液中の硝酸態窒素量を、銅・カドミウム還元－ナフチルエチレンジアミン吸光光度法によって定量する。
アンモニア態窒素量：
土壌試料に２ＭＫＣｌを添加し、１時間振とうした後、上澄み液をろ過し、抽出液中の硝酸態窒素量を、サリチル酸吸光光度法で定量する。
さらに、深層学習後の新規の被測定土壌の測定においては、第八のステップで取得した全波長スペクトルデータと、第九のステップで取得した所定の分析項目の予測値とを学習用データとしてもよい。
１ 土壌分析装置
１０ ＩＣＰ分析部
２０ 土壌分析部
２１ ＩＣＰ発光スペクトル処理部
２２ 制御部
３０ データ処理部
３１ 入出力インターフェース部
３２ 記憶装置
３３ 記憶部
３４ 予測部
４１ 入力データ
４２ 入力層
４２ａ 入力ユニット
４３ 中間層
４３ａ 中間ユニット
４４ 出力層
４５ 出力データ
１０ ＩＣＰ分析部
２０ 土壌分析部
２１ ＩＣＰ発光スペクトル処理部
２２ 制御部
３０ データ処理部
３１ 入出力インターフェース部
３２ 記憶装置
３３ 記憶部
３４ 予測部
４１ 入力データ
４２ 入力層
４２ａ 入力ユニット
４３ 中間層
４３ａ 中間ユニット
４４ 出力層
４５ 出力データ
Claims (6)
- 学習用データ取得のための被測定土壌を所定の抽出液を用いて土壌抽出液とする第一のステップと、
プラズマ発光分光分析により前記土壌抽出液の全波長スペクトルデータを取得すると共に、前記プラズマ発光分光分析により被測定土壌のＣａ、Ｍｇ、Ｎａ、Ｋからなる交換性塩基の実測値を取得する第二のステップと、
前記交換性塩基以外の分析項目の実測値を公定法により取得する第三のステップと、
前記全波長スペクトルデータと前記交換性塩基及び前記交換性塩基以外の分析項目の実測値とを用いて、深層学習により学習を行う第四のステップと、
新規の被測定土壌を所定の抽出液を用いて新規の土壌抽出液とする第七のステップと、
前記プラズマ発光分光分析により前記新規の土壌抽出液の全波長スペクトルデータを取得する第八のステップと、
前記第八のステップで取得した全波長スペクトルデータと、前記深層学習で取得した学習用データとを用いて、前記新規の被測定土壌の所定の分析項目の予測値を算出する第九のステップと、
を、備え、
前記交換性塩基以外の分析項目は、ｐＨ、電気伝導率、有効態リン、全炭素、全窒素、交換性Ａｌ、ＣＥＣ、粒径組成、可溶性Ｆｅ、交換性Ｍｎ、可溶性Ｚｎ、リン酸吸収係数、硝酸態窒素量及びアンモニア態窒素量を含み、
前記所定の分析項目は、交換性塩基（Ｃａ、Ｍｇ、Ｎａ、Ｋ）、ｐＨ、電気伝導率、有効態リン、全炭素、全窒素、交換性Ａｌ、ＣＥＣ、粒径組成、可溶性Ｆｅ、交換性Ｍｎ、可溶性Ｚｎ、リン酸吸収係数、硝酸態窒素量及びアンモニア態窒素量を含む、プラズマ発光分光分析を用いた土壌診断方法。 - 前記抽出液が、酢酸アンモニウム（ＮＨ４ＯＡｃ）溶液である、請求項１に記載のプラズマ発光分光分析を用いた土壌診断方法。
- 前記粒径組成が、粘土画分量及び砂画分量である、請求項１に記載のプラズマ発光分光分析を用いた土壌診断方法。
- 前記第二のステップにおいて、前記全波長スペクトルデータの強度を自然対数に変換し、かつ、前記全波長スペクトルデータの強度の内、飽和した強度とマイナスの強度とをそれぞれ所定の値として自然対数に変換する、請求項１に記載のプラズマ発光分光分析を用いた土壌診断方法。
- 前記第四のステップの後において、前記深層学習を終了してよいか否かの判定を行う第五のステップと、
前記第五のステップにおいて、深層学習を終了してよいと判定された場合には、前記深層学習に用いた学習用データをデータセットとして格納する第六のステップと、
を備え、
前記第九のステップでは、前記第八のステップで取得した全波長スペクトルデータと、前記第五のステップで取得した学習用データとを用いて、前記新規の被測定土壌の所定の分析項目の予測値を算出する、請求項１から４の何れかに記載のプラズマ発光分光分析を用いた土壌診断方法。 - 前記第八のステップで取得した新規の被測定土壌の全波長スペクトルデータと、前記第九のステップで取得した前記新規の被測定土壌の所定の分析項目の予測値とを新規な学習用データとする、請求項５に記載のプラズマ発光分光分析を用いた土壌診断方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2021092618A JP7464284B2 (ja) | 2021-06-01 | 2021-06-01 | プラズマ発光分光分析を用いた土壌診断方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2021092618A JP7464284B2 (ja) | 2021-06-01 | 2021-06-01 | プラズマ発光分光分析を用いた土壌診断方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2022184647A JP2022184647A (ja) | 2022-12-13 |
| JP7464284B2 true JP7464284B2 (ja) | 2024-04-09 |
Family
ID=84437567
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2021092618A Active JP7464284B2 (ja) | 2021-06-01 | 2021-06-01 | プラズマ発光分光分析を用いた土壌診断方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP7464284B2 (ja) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR102675928B1 (ko) * | 2023-03-14 | 2024-06-21 | 메타오가닉테크놀러지(주) | 초미세 기포를 이용한 비료 공급 방법 및 장치 |
| CN117147526B (zh) * | 2023-10-31 | 2024-02-06 | 西北师范大学 | 大面积土壤pH值快速检测方法 |
| WO2025105976A1 (en) * | 2023-11-15 | 2025-05-22 | Aramco Innovations LLC | Method for detecting trace metal anomalies for geochemical exploration |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2013064603A (ja) | 2011-09-15 | 2013-04-11 | Gunma Univ | 重金属分析装置及び重金属の分析方法 |
| JP2013524237A (ja) | 2010-04-05 | 2013-06-17 | アプライド・リサーチ・アソシエイツ，インコーポレーテッド | レーザ誘起ブレークダウン分光のための認識アルゴリズムを形成するための方法 |
| CN107044976A (zh) | 2017-05-10 | 2017-08-15 | 中国科学院合肥物质科学研究院 | 基于libs与堆叠式rbm深度学习技术的土壤重金属含量分析预测方法 |
| CN108596246A (zh) | 2018-04-23 | 2018-09-28 | 浙江科技学院 | 基于深度神经网络的土壤重金属含量检测模型的建立方法 |
| CN109374860A (zh) | 2018-11-13 | 2019-02-22 | 西北大学 | 一种基于机器学习算法的土壤养分预测及综合评价方法 |
| JP2022134096A (ja) | 2021-03-02 | 2022-09-14 | 国立研究開発法人農業・食品産業技術総合研究機構 | 土壌の可給態養分濃度の推定方法及び推定装置、並びに、推定モデルの生成方法及び生成装置 |
-
2021
- 2021-06-01 JP JP2021092618A patent/JP7464284B2/ja active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2013524237A (ja) | 2010-04-05 | 2013-06-17 | アプライド・リサーチ・アソシエイツ，インコーポレーテッド | レーザ誘起ブレークダウン分光のための認識アルゴリズムを形成するための方法 |
| JP2013064603A (ja) | 2011-09-15 | 2013-04-11 | Gunma Univ | 重金属分析装置及び重金属の分析方法 |
| CN107044976A (zh) | 2017-05-10 | 2017-08-15 | 中国科学院合肥物质科学研究院 | 基于libs与堆叠式rbm深度学习技术的土壤重金属含量分析预测方法 |
| CN108596246A (zh) | 2018-04-23 | 2018-09-28 | 浙江科技学院 | 基于深度神经网络的土壤重金属含量检测模型的建立方法 |
| CN109374860A (zh) | 2018-11-13 | 2019-02-22 | 西北大学 | 一种基于机器学习算法的土壤养分预测及综合评价方法 |
| JP2022134096A (ja) | 2021-03-02 | 2022-09-14 | 国立研究開発法人農業・食品産業技術総合研究機構 | 土壌の可給態養分濃度の推定方法及び推定装置、並びに、推定モデルの生成方法及び生成装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2022184647A (ja) | 2022-12-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP7464284B2 (ja) | プラズマ発光分光分析を用いた土壌診断方法 | |
| Guezenoc et al. | Critical review and advices on spectral-based normalization methods for LIBS quantitative analysis | |
| Zhang et al. | Predicting cadmium concentration in soils using laboratory and field reflectance spectroscopy | |
| Riedel et al. | Prediction of soil parameters using the spectral range between 350 and 15,000 nm: A case study based on the Permanent Soil Monitoring Program in Saxony, Germany | |
| Andrade et al. | Micronutrients prediction via pXRF spectrometry in Brazil: Influence of weathering degree | |
| Lebret et al. | High abundances of the nuisance raphidophyte Gonyostomum semen in brown water lakes are associated with high concentrations of iron | |
| Goranov et al. | TEnvR: MATLAB-based toolbox for environmental research | |
| CN115931828B (zh) | 适应复杂土壤基体的成分分析与预测方法、单元及系统 | |
| Araújo et al. | Soil contaminated with chromium by tannery sludge and identified by vis-NIR-mid spectroscopy techniques | |
| dos Santos et al. | Comparison between energy dispersive X-ray fluorescence spectral data and elemental data for soil attributes modelling | |
| Slonecker | Analysis of the effects of heavy metals on vegetation hyperspectral reflectance properties | |
| Wang et al. | Prediction of N, P, and K contents in sugarcane leaves by VIS-NIR spectroscopy and modeling of NPK interaction effects | |
| Luna et al. | Classification of soil samples based on Raman spectroscopy and X-ray fluorescence spectrometry combined with chemometric methods and variable selection | |
| Breure et al. | Transfer functions for phosphorus and potassium soil tests and implications for the QUEFTS model | |
| Farrar et al. | Rapid assessment of soil carbon and nutrients following application of organic amendments | |
| Srivastava et al. | Soil chemical properties estimation using hyperspectral remote sensing: A review | |
| CN121027016A (zh) | 一种用于河流环境保护的河道水质检测分析方法及装置 | |
| CN118169106B (zh) | 一种基于免定标libs技术的数据筛选与多元素定量分析方法 | |
| CN119438105A (zh) | 一种基于航天高光谱遥感的土壤盐碱化监测方法及系统 | |
| Zarei et al. | Development and application of novel risk indices for assessing heavy metal pollution in aquatic sediments | |
| Ichsan et al. | The best performing color space and machine learning regression algorithm for the accurate estimation of chromium (VI) and iron (III) in aqueous samples using low-cost and portable flatbed scanner colorimetry | |
| Hartono et al. | Development of quick test method for soil pH, nitrate, phosphorus, and potassium combining chemicals and phone cellular application | |
| Dibi et al. | Use of fluorescence and reflectance spectra for predicting okra (Abelmoschus esculentus) yield and macronutrient contents of leaves | |
| Mikkelsen et al. | Advances in assessing nutrient availability in soils | |
| Dasgupta et al. | Contemporary use of sensors for soil qualitative and quantitative assessment in the context of climate change |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20231030 |
|
| A871 | Explanation of circumstances concerning accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A871 Effective date: 20231030 |
|
| RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20231030 |
|
| RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20231101 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20231129 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20240116 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20240226 |
|
| 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: 20240312 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20240321 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7464284 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |