# 圃場試験解析

農学者は試験解析（Trial Analytics）を用いて、さまざまな作物品種、栽培技術、および投入資材の適用（精密農業における可変施用率の結果を含む）の性能を評価します。圃場試験から得られるデータを収集、分析、解釈することで、研究者は遺伝、環境、管理実践の相互作用に関する洞察を得ます。この知見は、投入資材の使用を最小限に抑えつつ収量ポテンシャルを最適化する作物管理戦略の開発に役立ちます。さらに、Trial Analyticsは精密農業の有効性を評価するだけでなく、多様で厳しい条件下でも生育可能な耐性のある品種を特定するのにも役立ち、食料安全保障への貢献につながります。

## データ準備

効果的な試験解析のために、いくつかの必須データセットが必要です：

1. **収量データセット**: \
   このデータセットは収量データを記録します。 \
   これは以下からインポートできます、または手動で [JohnDeere Operation Center](https://docs.geopard.tech/geopard-tutorials/jpn/tsu-webapuri/inpto/myjohndeere-karainpto) からインポートするか、手動でアップロードして [シェイプファイル](https://docs.geopard.tech/geopard-tutorials/jpn/tsu-webapuri/inpto/dtanoinpto) として、または [機械メーカー独自形式](https://docs.geopard.tech/geopard-tutorials/jpn/tsu-webapuri/inpto/nofmatto).
2. **施用データセット**: \
   これは圃場で実際に実行された施用を理解するために重要です。最低限、TargetRate、AppliedRate、およびいくつかの機械関連指標の属性を含みます。 \
   収量データセットと同様に、これも [JohnDeere Operation Center](https://docs.geopard.tech/geopard-tutorials/jpn/tsu-webapuri/inpto/myjohndeere-karainpto) からインポートするか、手動でアップロードして [シェイプファイル](https://docs.geopard.tech/geopard-tutorials/jpn/tsu-webapuri/inpto/midtanoinpto) として、または [機械メーカー独自形式](https://docs.geopard.tech/geopard-tutorials/jpn/tsu-webapuri/inpto/nofmatto).
3. **試験/実験のゾーン/区画**: \
   これらは試験の計画された施用率を示し、実験設計の洞察を与えます。 \
   そのようなデータレイヤーが利用可能であれば、をとしてアップロードします [シェイプファイル](https://docs.geopard.tech/geopard-tutorials/jpn/tsu-webapuri/inpto/midtanoinpto) をAsApplied/AsPlantedまたはYieldコントロールに取り込みます。これによりEquationMaps作成時の互換性が確保され、試験解析が効率化されます。\
   そのようなデータレイヤーがない場合、施用データセットのTargetRate属性が試験評価の代替として利用できます。
4. **過去の圃場ポテンシャルゾーン：**\
   これらのゾーンはGeoPardによって生成されます（詳細は [ここ](https://docs.geopard.tech/geopard-tutorials/jpn/tsu-webapuri/znmapputo/zn)）。これらは一貫した過去の生産性を持つ試験を解析する際に有用です。試験が生産性の異なる地域に分散されている場合、特に有益です。

これらのデータセットを収集したら、次のステップは試験評価プロセスの開始です。

## データ概要

冬小麦の2023年作シーズンについて、次のデータがあります：

* 湿重量分布を示す収量データセット *（図1）*

<figure><img src="https://1993076779-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2Fg9Iv1AMpACkykujxvJ9s%2Fimage.png?alt=media&#x26;token=c2a61d7b-12e9-40b5-ad88-e6f820372638" alt=""><figcaption><p>図1 元の収量データセット</p></figcaption></figure>

* 窒素（N34）可変施用計画（150 kg/ha）と2つの試験区（120kg/haおよび180 kg/ha）*（図2）*

<figure><img src="https://1993076779-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FwkGO1Tp0Ji2KzUPhfQAE%2Fimage.png?alt=media&#x26;token=28a94d78-7360-43e2-865f-4a330028ab20" alt=""><figcaption><p>図2 窒素（N34）VRA計画と試験区 </p></figcaption></figure>

* 適用統計を示す施用データセット *（図3）*

<figure><img src="https://1993076779-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FF1YKbPqDw2HBDdnTqazX%2Fimage.png?alt=media&#x26;token=0f83b5b0-9220-4cea-984c-735338716048" alt=""><figcaption><p>図3 施用データセット</p></figcaption></figure>

* 過去の圃場生産性（*図4*)

<figure><img src="https://1993076779-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FCrFaBaZxtjd0xEXZ9TUd%2Fimage.png?alt=media&#x26;token=9cf68b49-a17c-4f2f-be15-1b0bfee4993e" alt=""><figcaption><p>図4 過去の圃場生産性</p></figcaption></figure>

{% hint style="warning" %}
収量データセットはキャリブレーションされていません：複数のコンバインが稼働しており、ターンアラウンドや欠損データの痕跡があり、ノイズが見られます。最適な結果を得るには、収量のキャリブレートおよびクリーン操作を適用することを推奨します。 \
ステップバイステップのチュートリアルは [リンク](https://docs.geopard.tech/geopard-tutorials/jpn/nong-xue/notokurningu).
{% endhint %}

収量データセットはキャリブレーションとクリーン後、 *図5*に更新された統計とともに表示されます。このデータセットは後続のステップで使用されます。

<figure><img src="https://1993076779-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2F62J6z5wt7ndZ4lWxZq58%2Fimage.png?alt=media&#x26;token=9e6ea416-81a0-4dfc-a94a-aa19f69c0ff3" alt=""><figcaption><p>図5 キャリブレーションおよびクリーン済み収量データセット</p></figcaption></figure>

## コンセプト

ここでのTrial Analyticsの目的は、圃場にとって最も有効な窒素（N34）施用率を見定めることです。窒素施用率120kg/ha、150kg/ha、180kg/haに設定された特定領域があります。このデータは一方で施用データセット、他方でキャリブレートされた収量データセットから得られます。

解析は以下の3つの異なるゾーンに焦点を当てます：

* 120kg/ha（試験区として指定）
* 150kg/ha（主要区と見なす）
* 180kg/ha（別の試験区）

我々のアプローチには次の評価が含まれます：

1. **計画ベース：** キャリブレートされた収量に紐づく計画された可変施用（VRA）を使用
2. **実施ベース：** 実際に施用されたデータセットをキャリブレートされた収量と比較
3. **実施ベースと過去の生産性：** 実際に施用されたデータセットを、過去の圃場ポテンシャルゾーンと重ね合わせたキャリブレート収量と比較

この体系的アプローチにより、計画および実際の施用データの両方に基づいて、窒素が収量に与える影響を包括的に評価できます。

## 計画ベース

の影響は、 ~~施用された~~ 計画された窒素（N34）が収量分布に与える影響は、次のスクリーンショットに視覚的に示されています *（図6、図7、図8）*。以下は発見の簡潔な概要です：

* <mark style="color:青;">主要区の窒素率150 kg/haは面積45.8 haで平均収量は4.99 t/ha</mark> (*図6*).
* <mark style="color:青;">最初の試験区（窒素180 kg/ha）は面積1.76 haで平均収量6.5 t/ha</mark> (*図7*).
* <mark style="color:青;">2番目の試験区（窒素120 kg/ha）は面積1.86 haで平均収量6.39 t/ha</mark> (*図8*).

この結果は重要な疑問を投げかけます： <mark style="color:オレンジ;">なぜ低い施用率の方が高い施用率より効率的に見えるのか？</mark> より深い洞察を得るために、次の段階では[ 実際に施用されたデータを用いて試験を評価することになります](#applied-based-evaluation).

<figure><img src="https://1993076779-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FraVRCgxuEVTGfTKK4c0s%2Fimage.png?alt=media&#x26;token=57602ad8-bd9f-46b4-adcf-5de53f99fe7b" alt=""><figcaption><p>図6 N34 150kg/haの主要区</p></figcaption></figure>

<figure><img src="https://1993076779-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FNM80FFM2TQULwV54LuqW%2Fimage.png?alt=media&#x26;token=1542336e-20a2-4a77-a343-9a92f67c90eb" alt=""><figcaption><p>図7 N34 180kg/haの試験区</p></figcaption></figure>

<figure><img src="https://1993076779-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FITg4CWAJb27X28jquqyO%2Fimage.png?alt=media&#x26;token=cff8d0b3-5896-4f34-bbdc-6a34a7352982" alt=""><figcaption><p>図8 N34 120kg/haの試験区</p></figcaption></figure>

以下に、評価で使用した数式と設定についての詳細な議論があります。&#x20;

{% hint style="info" %}
方程式アプローチとその実行方法の詳細については、次の両方のチュートリアルを参照してください： [ユーザーインターフェース](https://docs.geopard.tech/geopard-tutorials/jpn/tsu-webapuri/bsu) と [API](https://docs.geopard.tech/geopard-tutorials/jpn/apidokyumento/furno/5-wo).
{% endhint %}

計算を再現するために実行する方程式は以下の通りです。&#x20;

1. 主要（150 kg/ha）： \
   `Yield_Main = np.where(Zone==1, Yield_WetMass, np.nan)`
2. 試験（120 kg/ha）： \
   `Yield_Zone = np.where(Zone==3, Yield_WetMass, np.nan)`
3. 試験（180 kg/ha）： \
   `Yield_Zone = np.where(Zone==2, Yield_WetMass, np.nan)`

重要なのは、 *Numpy* *を有効にすることです* および *補間（Interpolation）をオフにすることです* *（図10）*.

<figure><img src="https://1993076779-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FixfesbJQxzVC93JW3H7b%2Fimage.png?alt=media&#x26;token=dcab1611-8a33-4450-9b06-5ba24760e8c3" alt=""><figcaption><p>図9 「Numpy」を有効にする</p></figcaption></figure>

<figure><img src="https://1993076779-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2F0V4AxWhE0zraNfH1axHQ%2Fimage.png?alt=media&#x26;token=d86b3921-c313-4946-9455-35fbef29cc13" alt=""><figcaption><p>図10 「補間（Interpolated）データの使用」をオフにする</p></figcaption></figure>

## 実施ベース

注目すべき観察点は、試験中の実際の施用率が計画（ターゲット）施用率と一貫して一致していないことです。具体的には、120 kg/haから最大189 kg/haまで分布しています *（図11）*。この変動性を考慮し、誤差許容範囲の基準を設定することが重要となりました。したがって、試験を評価に適すると見なすために±5%の精度を許容閾値としました。

次のスクリーンショットに示されるように（*図12、図13、図14）* は、実際に施用された窒素（N34）量に焦点を当てた収量の統計分布です。以下は±5%精度受け入れを考慮した要約統計です：

* <mark style="color:青;">主要区の150 kg/haでは施用面積が43.5 haで平均収量は4.9 t/ha</mark> (*図12*).&#x20;
* <mark style="color:青;">最初の試験区180 kg/haは面積1.47 haで平均収量6.5 t/ha</mark> (*図13*).&#x20;
* <mark style="color:青;">2番目の試験区120 kg/haは面積1.44 haで平均収量6.3 t/ha</mark> (*図14*).

<figure><img src="https://1993076779-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FQ5Aawg78wAVNsybfq2RB%2Fimage.png?alt=media&#x26;token=293f63d1-414b-4e23-be09-b54844f21988" alt=""><figcaption><p>図11 試験における実際の施用率</p></figcaption></figure>

<figure><img src="https://1993076779-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FyBHR1SLHA9ujXr5wk8xA%2Fimage.png?alt=media&#x26;token=78f67150-560b-4495-974e-30407e25aca9" alt=""><figcaption><p>図12 N34 150kg/haの主要区 ±5%</p></figcaption></figure>

<figure><img src="https://1993076779-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FeOMxE76sp3ugO1NffWPN%2Fimage.png?alt=media&#x26;token=9e3f47ae-250a-41dd-b9ce-f3d4bcaab85e" alt=""><figcaption><p>図13 N34 180kg/haの試験区 ±5%</p></figcaption></figure>

<figure><img src="https://1993076779-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FJUBlkogPAF2fTSYa0EOW%2Fimage.png?alt=media&#x26;token=48dd3f24-df15-40da-905c-9012b70971d2" alt=""><figcaption><p>図14 N34 120kg/haの試験区 ±5%</p></figcaption></figure>

方法論とこれらの結果の詳細を深く理解するために、使用した方程式は以下の通りです：

1. 試験の実際の施用窒素：\
   `Applied_Trial = np.where((Zone == 3) | (Zone == 2), Applied_Value, np.nan)`
2. 150 kg/haの主要区に5%受け入れを組み込む： \
   `Yield_Main = np.where((Zone == 1) & (Applied_Value >= 142.5) & (Applied_Value <= 157.5), Yield_WetMass, np.nan)`
3. 120 kg/haの試験に5%受け入れを組み込む： \
   `Yield_Trial = np.where((Zone == 3) & (Applied_Value >= 114.0) & (Applied_Value <= 126.0), Yield_WetMass, np.nan)`
4. 180 kg/haの試験に5%受け入れを組み込む： \
   `Yield_Trial = np.where((Zone == 2) & (Applied_Value >= 171.0) & (Applied_Value <= 189.0), Yield_WetMass, np.nan)`

## **実施ベースと過去の生産性**

試験からの収量は継続的に圃場全体の平均収量を上回っています。この差異を生む主要因の一つは、試験が実施された過去に生産性の高かったゾーンであるようで、これは *図15* と *図16*に可視化されています。試験結果をより詳細に評価するには、解析時に生産性ゾーンを考慮することが重要です。

<figure><img src="https://1993076779-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FH0aA9ZAe2IkSFtqcJe5B%2Fimage.png?alt=media&#x26;token=c58e81d6-b2df-4d6d-8682-3c750f38f2a1" alt=""><figcaption><p>図15 過去の圃場ポテンシャルゾーン</p></figcaption></figure>

<figure><img src="https://1993076779-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2Fi2agD6onOT4aVaW9LvHh%2Fimage.png?alt=media&#x26;token=feecbfd5-2433-4061-85b9-8233446572d8" alt=""><figcaption><p>図16 過去の圃場ポテンシャルゾーンを収量データセットとして表示</p></figcaption></figure>

次のスクリーンショットに示されるように（*図17、図18、図19）* は、GeoPardで作成した過去の生産性ゾーンと重ね合わせた、実際に施用された窒素（N34）量に焦点を当てた収量の統計分布です。以下は、実施値に対する±5%の精度受け入れを考慮した要約統計です：

* <mark style="color:青;">主要区の150 kg/haでは施用面積が2.65 haで平均収量は6.34 t/ha</mark> (*図17*).
* <mark style="color:青;">最初の試験区180 kg/haは面積1.08 haで平均収量6.41 t/ha</mark> (*図18*).
* <mark style="color:青;">2番目の試験区120 kg/haは面積1.78 haで平均収量6.33 t/ha</mark> (*図19*).

<figure><img src="https://1993076779-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FPqw7u62bJAezDn9UW5qs%2Fimage.png?alt=media&#x26;token=ff101b42-7882-4dc3-a2ad-b57053900768" alt=""><figcaption><p>図17 N34 150kg/haの主要区（過去の生産性と重ね合わせ）</p></figcaption></figure>

<figure><img src="https://1993076779-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FmzkQRoU3HMMOm1fLLjNY%2Fimage.png?alt=media&#x26;token=ba8105da-877b-4234-9e1d-f59826949bf0" alt=""><figcaption><p>図18 N34 180kg/haの試験区 ±5%（過去の生産性と重ね合わせ）</p></figcaption></figure>

<figure><img src="https://1993076779-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FzTOI3v14HdZLANjM3o36%2Fimage.png?alt=media&#x26;token=f83a50dd-585d-4b66-937d-eb0d2f1c78f0" alt=""><figcaption><p>図19 N34 120kg/haの試験区 ±5%（過去の生産性と重ね合わせ）</p></figcaption></figure>

方法論とこれらの結果の詳細を深く理解するために、使用した方程式は以下の通りです：

1. 150 kg/haの主要区に5%受け入れを組み込み、過去の生産性と重ね合わせた場合： \
   `Yield_Main = np.where((Application_Zone == 1) & (Productivity_SubZone == 51) & (Applied_Value >= 142.5) & (Applied_Value <= 157.5), Yield_WetMass, np.nan)`
2. 120 kg/haの試験に5%受け入れを組み込み、過去の生産性と重ね合わせた場合： \
   `Yield_Trial = np.where((Application_Zone == 3) & (Productivity_SubZone == 51) & (Applied_Value >= 114.0) & (Applied_Value <= 126.0), Yield_WetMass, np.nan)`
3. 180 kg/haの試験に5%受け入れを組み込み、過去の生産性と重ね合わせた場合： \
   `Yield_Trial = np.where((Application_Zone == 2) & (Productivity_SubZone == 51) & (Applied_Value >= 171.0) & (Applied_Value <= 189.0), Yield_WetMass, np.nan)`

ここで&#x20;

* の部分は `Productivity_SubZone == 51` 適用実験が行われた高生産性ゾーンを反映しています、
* の部分は `(Applied_Value >= 142.5) & (Applied_Value <= 157.5)` , `(Applied_Value >= 114.0) & (Applied_Value <= 126.0)`,  `(Applied_Value >= 171.0) & (Applied_Value <= 189.0)` これらは施用率からの±5%の精度を組み込んでいます `150`, `120`, `180` kg/ha。

## 要約

試験から得られた収量結果は、圃場の高い過去の生産性ゾーンで観察される平均収量と概ね一致しています。言い換えれば、N34製品を <mark style="color:青;">120 kg/ha - 150 kg/ha - 180 kg/ha</mark>の各施用率で実験的に適用した結果、平均収量が <mark style="color:青;">6.33 t/ha - 6.34 t/ha - 6.41 t/ha</mark> それぞれとなり、高生産性ゾーン内の収穫収量に対して有意な影響は見られません。