# Analýza polních pokusů

Agronomové používají analytiku pokusů (Trial Analytics) k hodnocení výkonu různých odrůd plodin, technik pěstování a aplikací vstupů, včetně výsledků variabilních dávek (VRA) v precizním zemědělství. Sběrem, analýzou a interpretací dat z polních pokusů získávají výzkumníci poznatky o interakcích mezi genetikou, prostředím a postupy řízení. Tyto poznatky pak nasměrují vývoj strategie hospodaření na plodinách tak, aby se maximalizoval výnosový potenciál při minimalizaci spotřeby vstupů. Dále analytika pokusů nejen umožňuje zhodnotit účinnost postupů precizního zemědělství, ale také pomáhá identifikovat odrůdy odolné vůči různým a náročným podmínkám, čímž přispívá k potravinové bezpečnosti.

## Příprava dat

Pro efektivní analýzu pokusů jsou potřeba tyto základní datové sady:

1. **Sada dat výnosu**: \
   Tato sada zachycuje data o výnosu. \
   Můžeme ji importovat z [JohnDeere Operation Center](https://docs.geopard.tech/geopard-tutorials/cze/prohlidka-produktu-webova-aplikace/import/import-z-myjohndeere) nebo ji ručně nahrát jako [shapefile](https://docs.geopard.tech/geopard-tutorials/cze/prohlidka-produktu-webova-aplikace/import/import-dat-o-vynosu) nebo jako [vlastní formát strojů](https://docs.geopard.tech/geopard-tutorials/cze/prohlidka-produktu-webova-aplikace/import/proprietarni-formaty-stroju).
2. **Sada dat o aplikaci**: \
   Toto je klíčové pro pochopení skutečně provedené aplikace na poli. Minimálně obsahuje atributy jako TargetRate, AppliedRate a některé metriky související s technikou. \
   Stejně jako u sady dat výnosu máme možnost importovat ji z [JohnDeere Operation Center](https://docs.geopard.tech/geopard-tutorials/cze/prohlidka-produktu-webova-aplikace/import/import-z-myjohndeere) nebo ji ručně nahrát jako [shapefile](https://docs.geopard.tech/geopard-tutorials/cze/prohlidka-produktu-webova-aplikace/import/import-dat-as-applied-as-planted) nebo jako [vlastní formát strojů](https://docs.geopard.tech/geopard-tutorials/cze/prohlidka-produktu-webova-aplikace/import/proprietarni-formaty-stroju).
3. **zóny/parcely s pokusy/experimenty**: \
   Tyto ukazují plánované aplikační dávky pro naše pokusy a dávají přehled o experimentálním návrhu. \
   Pokud je taková datová vrstva k dispozici, nahrajeme ji jako [shapefile](https://docs.geopard.tech/geopard-tutorials/cze/prohlidka-produktu-webova-aplikace/import/import-dat-as-applied-as-planted) do kontroly AsApplied/AsPlanted nebo Yield. To zajišťuje kompatibilitu při tvorbě EquationMaps a zjednodušuje analýzu pokusů.\
   Pokud taková datová vrstva není k dispozici, atribut TargetRate ze sady dat o aplikaci může sloužit jako náhrada pro vyhodnocení pokusů.
4. **Historické zóny potenciálu pole:**\
   Tyto zóny generuje GeoPard (podrobnosti jsou [ZDE](https://docs.geopard.tech/geopard-tutorials/cze/prohlidka-produktu-webova-aplikace/mapy-zon-a-analyzy/viceleta-zonovani)). Jsou užitečné pro analýzu pokusů s konzistentní historickou produktivitou. To je zvláště přínosné, pokud jsou pokusy rozmístěny v regionech s rozdílnou historickou produktivitou.

Jakmile tyto datové sady sesbíráme, dalším krokem je zahájení procesu hodnocení pokusů.

## Přehled dat

Pro zemědělskou sezónu 2023 ozimé pšenice jsou k dispozici následující data:

* Sada dat výnosu zdůrazňující rozložení mokré hmoty *(Obr.1)*

<figure><img src="https://1171283066-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>Obr.1 Původní YieldDataset</p></figcaption></figure>

* Plán VRA dusíku (N34) (150 kg/ha) se 2 zkušebními parcelami (120 kg/ha a 180 kg/ha)*(Obr.2)*

<figure><img src="https://1171283066-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>Obr.2 Plán VRA dusíku (N34) se zkušebními parcelami </p></figcaption></figure>

* Sada dat o aplikaci ukazující aplikovanou statistiku *(Obr.3)*

<figure><img src="https://1171283066-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>Obr.3 Application Dataset</p></figcaption></figure>

* Historická produktivita pole (*Obr.4*)

<figure><img src="https://1171283066-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>Obr.4 Historická produktivita pole</p></figcaption></figure>

{% hint style="warning" %}
YieldDataset nebyla kalibrována: pracovaly tam více kombajnů, jsou patrné stopy po otáčení a chybějících datech a je viditelný šum. Doporučuje se aplikovat operace Kalibrace výnosu a Čištění pro optimální výsledky. \
Krok za krokem návod najdete na [ODKAZ](https://docs.geopard.tech/geopard-tutorials/cze/agronomie/kalibrace-a-cisteni-vynosu).
{% endhint %}

YieldDataset po kalibraci a čištění je zobrazena v *Obr.5*, spolu s aktualizovanými statistikami. Tato sada dat bude použita v následujících krocích.

<figure><img src="https://1171283066-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>Obr.5 Kalibrovaný a vyčištěný YieldDataset</p></figcaption></figure>

## Koncept

Cílem analytiky pokusů je zde určit nejúčinnější dávku dusíku (N34) pro pole. Jsou identifikovány oblasti s dávkami dusíku 120 kg/ha, 150 kg/ha a 180 kg/ha. Tato data pocházejí na jedné straně ze sady dat o aplikaci a na druhé straně z kalibrované sady dat výnosu.

Zaměřujeme svoji analýzu na tři odlišné zóny:

* 120 kg/ha (označeno jako zkušební zóna)
* 150 kg/ha (považováno za hlavní zónu)
* 180 kg/ha (další zkušební zóna)

Náš přístup bude zahrnovat následující hodnocení:

1. **Na základě plánu:** použití plánované variabilní aplikace (VRA) vázané na kalibrovaný výnos.
2. **Na základě aplikace:** srovnání skutečně aplikovaných dat s kalibrovaným výnosem.
3. **Na základě aplikace a historické produktivity:** srovnání skutečně aplikovaných dat s kalibrovaným výnosem překrytým historickými zónami potenciálu pole.

Tento metodický přístup umožní komplexní posouzení vlivu dusíku na výnos na základě jak plánovaných, tak skutečně aplikovaných dat.

## Na základě plánu

Vliv ~~aplikovaného~~ plánovaného dusíku (N34) na rozložení výnosu je vizuálně zachycen v následujících snímcích obrazovky *(Obr.6, Obr.7, Obr.8)*. Zde je stručné shrnutí zjištění:

* <mark style="color:modrá;">Hlavní zóna s dávkou dusíku 150 kg/ha zabírá 45,8 ha a průměrný výnos je 4,99 t/ha</mark> (*Obr.6*).
* <mark style="color:modrá;">První zkušební zóna, s aplikací 180 kg/ha dusíku, pokrývá 1,76 ha a průměrný výnos je 6,5 t/ha</mark> (*Obr.7*).
* <mark style="color:modrá;">Druhá zkušební zóna s 120 kg/ha dusíku má 1,86 ha a průměrný výnos 6,39 t/ha</mark> (*Obr.8*).

Výsledky vyvolávají zásadní otázku: <mark style="color:oranžová;">Proč se nižší aplikační dávka jeví jako efektivnější než vyšší?</mark> Pro hlubší porozumění zahrnuje další fáze[ hodnocení pokusů pomocí skutečně aplikovaných dat](#applied-based-evaluation).

<figure><img src="https://1171283066-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>Obr.6 Hlavní zóna s N34 150 kg/ha</p></figcaption></figure>

<figure><img src="https://1171283066-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>Obr.7 Zkušební zóna s N34 180 kg/ha</p></figcaption></figure>

<figure><img src="https://1171283066-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>Obr.8 Zkušební zóna s N34 120 kg/ha</p></figcaption></figure>

Níže najdete podrobnou diskusi o vzorcích a konfiguracích použitých při hodnocení.&#x20;

{% hint style="info" %}
Pro bližší seznámení s přístupem Equation a jeho provedením si prosím prostudujte naše návody jak pro [uživatelské rozhraní](https://docs.geopard.tech/geopard-tutorials/cze/prohlidka-produktu-webova-aplikace/analytika-zalozena-na-rovnicich) a [API](https://docs.geopard.tech/geopard-tutorials/cze/api-dokumentace/schemata-se-zakladnimi-toky/5.-spusteni-rovnic).
{% endhint %}

Zde jsou rovnice, které je třeba spustit pro reprodukci výpočtů.&#x20;

1. Hlavní s 150 kg/ha: \
   `Yield_Main = np.where(Zone==1, Yield_WetMass, np.nan)`
2. Zkušební s 120 kg/ha: \
   `Yield_Zone = np.where(Zone==3, Yield_WetMass, np.nan)`
3. Zkušební s 180 kg/ha: \
   `Yield_Zone = np.where(Zone==2, Yield_WetMass, np.nan)`

Je důležité aktivovat *Numpy* *(Obr.9)* a vypnout *Interpolaci* *(Obr.10)*.

<figure><img src="https://1171283066-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>Obr.9 Aktivovat "Numpy"</p></figcaption></figure>

<figure><img src="https://1171283066-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>Obr.10 Vypnout použití interpolovaných dat</p></figcaption></figure>

## Na základě aplikace

Významné zjištění je, že skutečná aplikovaná dávka během pokusu není vždy v souladu s plánovanou (cílovou) dávkou. Konkrétně se rozložení pohybuje od 120 kg/ha až do 189 kg/ha *(Obr.11)*. Vzhledem k této variabilitě bylo nutné stanovit mez chybovosti. Proto byla považována za přijatelnou tolerance ±5 %, aby byl pokus vhodný k vyhodnocení.

Prezentováno v následujících snímcích obrazovky (*Obr.12, Obr.13, Obr.14)* je statistické rozdělení výnosu se zaměřením na skutečně aplikovaná čísla dusíku (N34). Zde jsou shrnuté statistiky s ohledem na přijatelnou přesnost ±5 %:

* <mark style="color:modrá;">Hlavní zóna při 150 kg/ha měla aplikovanou plochu 43,5 ha a průměrný výnos 4,9 t/ha</mark> (*Obr.12*).&#x20;
* <mark style="color:modrá;">První zkušební zóna při 180 kg/ha pokryla plochu 1,47 ha a dosáhla průměrného výnosu 6,5 t/ha</mark> (*Obr.13*).&#x20;
* <mark style="color:modrá;">Druhá zkušební zóna nastavená na 120 kg/ha zabírala plochu 1,44 ha s průměrným výnosem 6,3 t/ha</mark> (*Obr.14*).

<figure><img src="https://1171283066-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>Obr.11 Skutečně aplikované dávky v pokusech</p></figcaption></figure>

<figure><img src="https://1171283066-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>Obr.12 Hlavní zóna s N34 150 kg/ha ±5%</p></figcaption></figure>

<figure><img src="https://1171283066-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>Obr.13 Zkušební zóna s N34 180 kg/ha ±5%</p></figcaption></figure>

<figure><img src="https://1171283066-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>Obr.14 Zkušební zóna s N34 120 kg/ha ±5%</p></figcaption></figure>

Pro hlubší pochopení metodiky a specifik těchto výsledků jsou níže použité rovnice:

1. Skutečně aplikovaný dusík v pokusu:\
   `Applied_Trial = np.where((Zone == 3) | (Zone == 2), Applied_Value, np.nan)`
2. Hlavní s 150 kg/ha s 5% akceptací: \
   `Yield_Main = np.where((Zone == 1) & (Applied_Value >= 142.5) & (Applied_Value <= 157.5), Yield_WetMass, np.nan)`
3. Zkušební s 120 kg/ha s 5% akceptací: \
   `Yield_Trial = np.where((Zone == 3) & (Applied_Value >= 114.0) & (Applied_Value <= 126.0), Yield_WetMass, np.nan)`
4. Zkušební s 180 kg/ha s 5% akceptací: \
   `Yield_Trial = np.where((Zone == 2) & (Applied_Value >= 171.0) & (Applied_Value <= 189.0), Yield_WetMass, np.nan)`

## **Na základě aplikace a historické produktivity**

Čísla výnosů z pokusů konzistentně převyšují průměrný výnos napříč celým polem. Klíčovým faktorem, který tento rozdíl způsobuje, se zdá být historicky vysoce produktivní zóna, kde byly pokusy prováděny, jak je vizualizováno v *Obr.15* a *Obr.16*. Pro nuance v hodnocení pokusů je důležité při analýze výsledků zohlednit produktivní zóny.

<figure><img src="https://1171283066-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>Obr.15 Historické zóny potenciálu pole</p></figcaption></figure>

<figure><img src="https://1171283066-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>Obr.16 Historické zóny potenciálu pole jako YieldDataset</p></figcaption></figure>

Prezentováno v následujících snímcích obrazovky (*Obr.17, Obr.18, Obr.19)* je statistické rozdělení výnosu se zaměřením na skutečně aplikovaná čísla dusíku (N34) překrytá historickými zónami produktivity (vytvořenými v GeoPard). Zde jsou shrnuté statistiky s ohledem na přijatelnou přesnost ±5 % pro aplikovaná čísla:

* <mark style="color:modrá;">Hlavní zóna při 150 kg/ha měla aplikovanou plochu 2,65 ha a průměrný výnos 6,34 t/ha</mark> (*Obr.17*).
* <mark style="color:modrá;">První zkušební zóna při 180 kg/ha pokryla plochu 1,08 ha a dosáhla průměrného výnosu 6,41 t/ha</mark> (*Obr.18*).
* <mark style="color:modrá;">Druhá zkušební zóna nastavená na 120 kg/ha zabírala plochu 1,78 ha s průměrným výnosem 6,33 t/ha</mark> (*Obr.19*).

<figure><img src="https://1171283066-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>Obr.17 Hlavní zóna s N34 150 kg/ha překrytá historickou produktivitou</p></figcaption></figure>

<figure><img src="https://1171283066-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>Obr.18 Zkušební zóna s N34 180 kg/ha ±5% překrytá historickou produktivitou</p></figcaption></figure>

<figure><img src="https://1171283066-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>Obr.19 Zkušební zóna s N34 120 kg/ha ±5% překrytá historickou produktivitou</p></figcaption></figure>

Pro hlubší pochopení metodiky a specifik těchto výsledků jsou níže použité rovnice:

1. Hlavní s 150 kg/ha s 5% akceptací překrytá historickou produktivitou: \
   `Yield_Main = np.where((Application_Zone == 1) & (Productivity_SubZone == 51) & (Applied_Value >= 142.5) & (Applied_Value <= 157.5), Yield_WetMass, np.nan)`
2. Zkušební s 120 kg/ha s 5% akceptací překrytá historickou produktivitou: \
   `Yield_Trial = np.where((Application_Zone == 3) & (Productivity_SubZone == 51) & (Applied_Value >= 114.0) & (Applied_Value <= 126.0), Yield_WetMass, np.nan)`
3. Zkušební s 180 kg/ha s 5% akceptací překrytá historickou produktivitou: \
   `Yield_Trial = np.where((Application_Zone == 2) & (Productivity_SubZone == 51) & (Applied_Value >= 171.0) & (Applied_Value <= 189.0), Yield_WetMass, np.nan)`

kde&#x20;

* část `Productivity_SubZone == 51` odráží vysoké produktivní zóny, ve kterých proběhly aplikované experimenty,
* části `(Applied_Value >= 142.5) & (Applied_Value <= 157.5)` , `(Applied_Value >= 114.0) & (Applied_Value <= 126.0)`,  `(Applied_Value >= 171.0) & (Applied_Value <= 189.0)` zahrnují přesnost ±5 % od hodnot dávek `150`, `120`, `180` kg/ha.

## Shrnutí

Výsledky výnosů z pokusů se úzce shodují s průměrným výnosem pozorovaným v rámci vysoce historicky produktivní zóny pole. Jinými slovy, experimentální aplikace produktu N34 při dávkách <mark style="color:modrá;">120 kg/ha - 150 kg/ha - 180 kg/ha</mark>, vedla k průměrným výnosům <mark style="color:modrá;">6,33 t/ha - 6,34 t/ha - 6,41 t/ha</mark> v uvedeném pořadí, a nemá v rámci vysoce produktivní zóny pole výrazný vliv na sklizený výnos.