# Analitika poskusov na polju

Agronomi uporabljajo analitiko poskusov za ocenjevanje uspešnosti različnih sort poljščin, tehnike pridelave in aplikacij vhodov, vključno z rezultati variabilnih aplikacij odmerka v preciznem kmetijstvu. Z zbiranjem, analizo in interpretacijo podatkov, ustvarjenih iz poljskih poskusov, raziskovalci pridobijo vpogled v interakcije med genetiko, okoljem in praksami upravljanja. To znanje usmerja razvoj strategij upravljanja posevkov, ki optimizirajo potencial pridelka ob hkratnem zmanjševanju porabe vhodov. Poleg tega analitika poskusov ne omogoča le ocenjevanja učinkovitosti praks preciznega kmetovanja, temveč tudi pomaga pri prepoznavanju odpornih sort poljščin, ki lahko uspevajo v raznolikih in zahtevnih razmerah, ter s tem prispeva k prehranski varnosti.

{% hint style="info" %}
GeoPard podpira tudi split-plot poskuse za dvofaktorske zasnove, kot so enake količine pri različnih hibridih. Ocenite tako glavne učinke kot interakcijo količina × hibrid na istem georeferenciranem sloju poskusa.
{% endhint %}

## Priprava podatkov

Za učinkovito analitiko poskusov so potrebni nekateri osnovni nabori podatkov:

1. **Nabor podatkov o pridelku**:\
   Ta nabor podatkov zajema podatke o pridelku.\
   Uvozimo ga lahko iz [JohnDeere Operation Center](https://docs.geopard.tech/geopard-tutorials/svn/predstavitev-izdelka-spletna-aplikacija/uvoz/uvoz-iz-myjohndeere) ali pa ga ročno naložimo kot [shapefile](https://docs.geopard.tech/geopard-tutorials/svn/predstavitev-izdelka-spletna-aplikacija/uvoz/uvoz-podatkov-o-pridelku) ali kot [proprietarni format stroja](https://docs.geopard.tech/geopard-tutorials/svn/predstavitev-izdelka-spletna-aplikacija/uvoz/lastniski-formati-strojev).
2. **Nabor podatkov o aplikaciji**:\
   To je ključnega pomena za razumevanje dejansko izvedene aplikacije na polju. Najmanj vsebuje atribute, kot so TargetRate, AppliedRate in nekatere metrike, povezane s strojem.\
   Tako kot pri naboru podatkov o pridelku ga lahko uvozimo iz [JohnDeere Operation Center](https://docs.geopard.tech/geopard-tutorials/svn/predstavitev-izdelka-spletna-aplikacija/uvoz/uvoz-iz-myjohndeere) ali pa ga ročno naložimo kot [shapefile](https://docs.geopard.tech/geopard-tutorials/svn/predstavitev-izdelka-spletna-aplikacija/uvoz/uvoz-podatkov-o-izvedeni-aplikaciji-setvi) ali kot [proprietarni format stroja](https://docs.geopard.tech/geopard-tutorials/svn/predstavitev-izdelka-spletna-aplikacija/uvoz/lastniski-formati-strojev).
3. **Cone/parcele s poskusi/eksperimenti**:\
   Te prikazujejo načrtovane količine aplikacije za naše poskuse in dajejo vpogled v poskusno zasnovo.\
   Če je tak podatkovni sloj na voljo, ga naložimo kot [shapefile](https://docs.geopard.tech/geopard-tutorials/svn/predstavitev-izdelka-spletna-aplikacija/uvoz/uvoz-podatkov-o-izvedeni-aplikaciji-setvi) v kontrolnik AsApplied/AsPlanted ali Yield. To zagotavlja združljivost pri gradnji EquationMaps in poenostavi izkušnjo analitike poskusov.\
   To je lahko zasnova z enim faktorjem ali split-plot zasnova z drugo dimenzijo obravnave, kot je hibrid ali sorta.\
   Če tak podatkovni sloj ni na voljo, lahko atribut TargetRate iz nabora podatkov o aplikaciji služi kot nadomestek za ocene poskusov.
4. **Cone zgodovinskega potenciala polja:**\
   Te cone generira GeoPard (podrobnosti so [TUKAJ](https://docs.geopard.tech/geopard-tutorials/svn/predstavitev-izdelka-spletna-aplikacija/karte-con-in-analitika/vecletne-cone)). Uporabne so za analizo poskusov z dosledno zgodovinsko produktivnostjo. To je še posebej koristno, kadar so poskusi razporejeni po regijah z različno zgodovinsko produktivnostjo.

Ko zberemo te nabore podatkov, je naslednji korak začetek procesa ocenjevanja poskusa.

## Pregled podatkov

Za kmetijsko sezono 2023 pri ozimni pšenici so na voljo naslednji podatki:

* Nabor podatkov o pridelku, ki prikazuje porazdelitev mokre mase *(Slika 1)*

<figure><img src="https://29491374-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>Slika 1 Izvorni YieldDataset</p></figcaption></figure>

* Načrt VRA za dušik (N34) (150 kg/ha) z 2 poskusnima parcelama (120 kg/ha in 180 kg/ha)*(Slika 2)*

<figure><img src="https://29491374-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>Slika 2 Načrt VRA za dušik (N34) s poskusnimi parcelami</p></figcaption></figure>

* Nabor podatkov o aplikaciji, ki prikazuje uporabljene statistike *(Slika 3)*

<figure><img src="https://29491374-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>Slika 3 Nabor podatkov o aplikaciji</p></figcaption></figure>

* Zgodovinska produktivnost polja (*Slika 4*)

<figure><img src="https://29491374-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>Slika 4 Zgodovinska produktivnost polja</p></figcaption></figure>

{% hint style="warning" %}
YieldDataset ni bil umerjen: tam je delovalo več kombajnov, prisotni so obrati na koncu vrste in sledi manjkajočih podatkov, vidna pa je tudi šumovitost. Za optimalne rezultate je priporočljivo uporabiti operaciji Yield Calibrate in Clean.\
Korak za korakom vodič je na voljo na [POVEZAVI](https://docs.geopard.tech/geopard-tutorials/svn/agronomija/kalibracija-in-ciscenje-pridelka).
{% endhint %}

YieldDataset je po umerjanju in čiščenju prikazan na *Sliki 5*, skupaj s posodobljenimi statistikami. Ta nabor podatkov bo uporabljen v naslednjih korakih.

<figure><img src="https://29491374-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>Slika 5 Umerjen in očiščen YieldDataset</p></figcaption></figure>

## Koncept

Cilj analitike poskusov je določiti najučinkovitejši odmerek dušika (N34) za polje. Označena so območja z odmerki dušika 120 kg/ha, 150 kg/ha in 180 kg/ha. Ti podatki izhajajo na eni strani iz ApplicationDataset, na drugi pa iz umerjenega YieldDataset.

Analizo osredotočamo na tri ločene cone:

* 120 kg/ha (označeno kot poskusna cona)
* 150 kg/ha (šteje kot glavna cona)
* 180 kg/ha (še ena poskusna cona)

Naš pristop bo vključeval naslednje ocene:

1. **Na podlagi načrta:** z uporabo načrtovane variabilne aplikacije odmerka (VRA), povezane z umerjenim pridelkom.
2. **Na podlagi dejanske aplikacije:** primerjava dejanskih podatkov o aplikaciji z umerjenim pridelkom.
3. **Na podlagi dejanske aplikacije in zgodovinske produktivnosti:** primerjava dejanskih podatkov o aplikaciji z umerjenim pridelkom, prekritim s conami zgodovinskega potenciala polja.

Ta metodičen pristop bo omogočil celovito oceno vpliva dušika na pridelek, na podlagi tako načrtovanih kot dejansko apliciranih podatkov o aplikaciji.

## Na podlagi načrta

Vpliv ~~dejansko apliciranega~~ načrtovanega dušika (N34) na porazdelitev pridelka je vizualno prikazan na naslednjih posnetkih zaslona *(Slika 6, Slika 7, Slika 8)*. Tukaj je kratek povzetek ugotovitev:

* <mark style="color:modra;">Glavna cona z odmerkom dušika 150 kg/ha zajema 45,8 ha in ima povprečen pridelek 4,99 t/ha</mark> (*Slika 6*).
* <mark style="color:modra;">Prva poskusna cona z uporabo aplikacije dušika 180 kg/ha pokriva 1,76 ha in daje povprečno 6,5 t/ha</mark> (*Slika 7*).
* <mark style="color:modra;">Druga poskusna cona z 120 kg/ha dušika obsega 1,86 ha in daje povprečen pridelek 6,39 t/ha</mark> (*Slika 8*).

Rezultati postavljajo pomembno vprašanje: <mark style="color:oranžna;">Zakaj se zdi, da je nižji odmerek učinkovitejši od višjega?</mark> Za globlji vpogled naslednja faza vključuje[ ocenjevanje poskusov z dejanskimi podatki o aplikaciji](#applied-based-evaluation).

<figure><img src="https://29491374-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>Slika 6 Glavna cona z N34 150kg/ha</p></figcaption></figure>

<figure><img src="https://29491374-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>Slika 7 Poskusna cona z N34 180kg/ha</p></figcaption></figure>

<figure><img src="https://29491374-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>Slika 8 Poskusna cona z N34 120kg/ha</p></figcaption></figure>

Spodaj je podrobna razprava o formulah in nastavitvah, uporabljenih med ocenjevanjem.

{% hint style="info" %}
Za podrobnejši vpogled v pristop Equation in njegovo izvedbo si oglejte naše vodiče za [uporabniški vmesnik](https://docs.geopard.tech/geopard-tutorials/svn/predstavitev-izdelka-spletna-aplikacija/analitika-na-podlagi-enacb) in [API](https://docs.geopard.tech/geopard-tutorials/svn/dokumentacija-api/diagrami-z-osnovnimi-tokovi/5.-izvedi-enacbe).
{% endhint %}

Tukaj so enačbe, ki jih je treba zagnati za ponovitev izračunov.

1. Glavna cona s 150 kg/ha:\
   `Yield_Main = np.where(Zone==1, Yield_WetMass, np.nan)`
2. Poskus s 120 kg/ha:\
   `Yield_Zone = np.where(Zone==3, Yield_WetMass, np.nan)`
3. Poskus s 180 kg/ha:\
   `Yield_Zone = np.where(Zone==2, Yield_WetMass, np.nan)`

Pomembno je aktivirati *Numpy* *(Slika 9)* in izklopiti *interpolacijo* *(Slika 10)*.

<figure><img src="https://29491374-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>Slika 9 Aktiviraj "Numpy"</p></figcaption></figure>

<figure><img src="https://29491374-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>Slika 10 Izklopi uporabo podatkov "Interpolated"</p></figcaption></figure>

## Na podlagi dejanske aplikacije

Pomembna ugotovitev je, da se dejanska stopnja aplikacije med poskusom ne ujema dosledno z načrtovano (ciljno) stopnjo. Natančneje, porazdelitev sega od 120 kg/ha do celo 189 kg/ha *(Slika 11)*. Glede na to spremenljivost je bilo nujno določiti referenčno mejo za toleranco napake. Tako se je ±5 % natančnost izkazala za sprejemljiv prag, da se poskus šteje za primernega za ocenjevanje.

Prikazano na naslednjih posnetkih zaslona (*Slika 12, Slika 13, Slika 14)* je statistična porazdelitev pridelka, s poudarkom na dejansko apliciranih količinah dušika (N34). Tukaj so povzete statistike, ob upoštevanju sprejemljive ±5 % natančnosti:

* <mark style="color:modra;">Glavna cona pri 150 kg/ha je imela aplicirano površino 43,5 ha in povprečen pridelek 4,9 t/ha</mark> (*Slika 12*).
* <mark style="color:modra;">Prva poskusna cona pri 180 kg/ha je pokrivala 1,47 ha in dala povprečen pridelek 6,5 t/ha</mark> (*Slika 13*).
* <mark style="color:modra;">Druga poskusna cona pri 120 kg/ha je obsegala 1,44 ha s povprečnim pridelkom 6,3 t/ha</mark> (*Slika 14*).

<figure><img src="https://29491374-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>Slika 11 Dejanske stopnje aplikacije v poskusih</p></figcaption></figure>

<figure><img src="https://29491374-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>Slika 12 Glavna cona z N34 150kg/ha ±5%</p></figcaption></figure>

<figure><img src="https://29491374-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>Slika 13 Poskusna cona z N34 180kg/ha ±5%</p></figcaption></figure>

<figure><img src="https://29491374-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>Slika 14 Poskusna cona z N34 120kg/ha ±5%</p></figcaption></figure>

Za globlje razumevanje metodologije in podrobnosti teh rezultatov so spodaj uporabljene enačbe:

1. Dejanski apliciran dušik v poskusu:\
   `Applied_Trial = np.where((Zone == 3) | (Zone == 2), Applied_Value, np.nan)`
2. Glavna cona s 150 kg/ha z upoštevanjem 5 % sprejemljivosti:\
   `Yield_Main = np.where((Zone == 1) & (Applied_Value >= 142.5) & (Applied_Value <= 157.5), Yield_WetMass, np.nan)`
3. Poskus s 120 kg/ha z upoštevanjem 5 % sprejemljivosti:\
   `Yield_Trial = np.where((Zone == 3) & (Applied_Value >= 114.0) & (Applied_Value <= 126.0), Yield_WetMass, np.nan)`
4. Poskus s 180 kg/ha z upoštevanjem 5 % sprejemljivosti:\
   `Yield_Trial = np.where((Zone == 2) & (Applied_Value >= 171.0) & (Applied_Value <= 189.0), Yield_WetMass, np.nan)`

## **Na podlagi dejanske aplikacije in zgodovinske produktivnosti**

Pridelek iz poskusov dosledno presega povprečni pridelek na celotnem polju. Ključni dejavnik za to razliko se zdi zgodovinsko visoko produktivna cona, kjer so potekali poskusi, kar je prikazano na *Sliki 15* in *Sliki 16*. Za bolj niansirano oceno poskusov je pri analizi rezultatov ključnega pomena upoštevati produktivne cone.

<figure><img src="https://29491374-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>Slika 15 Cone zgodovinskega potenciala polja</p></figcaption></figure>

<figure><img src="https://29491374-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>Slika 16 Cone zgodovinskega potenciala polja kot YieldDataset</p></figcaption></figure>

Prikazano na naslednjih posnetkih zaslona (*Slika 17, Slika 18, Slika 19)* je statistična porazdelitev pridelka, s poudarkom na dejansko apliciranih količinah dušika (N34), prekritih s conami zgodovinske produktivnosti (ustvarjenimi v GeoPardu). Tukaj so povzete statistike, ob upoštevanju sprejemljive ±5 % natančnosti za aplicirane količine:

* <mark style="color:modra;">Glavna cona pri 150 kg/ha je imela aplicirano površino 2,65 ha in povprečen pridelek 6,34 t/ha</mark> (*Slika 17*).
* <mark style="color:modra;">Prva poskusna cona pri 180 kg/ha je pokrivala 1,08 ha in dala povprečen pridelek 6,41 t/ha</mark> (*Slika 18*).
* <mark style="color:modra;">Druga poskusna cona pri 120 kg/ha je obsegala 1,78 ha s povprečnim pridelkom 6,33 t/ha</mark> (*Slika 19*).

<figure><img src="https://29491374-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>Slika 17 Glavna cona z N34 150kg/ha prekrita z zgodovinsko produktivnostjo</p></figcaption></figure>

<figure><img src="https://29491374-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>Slika 18 Poskusna cona z N34 180kg/ha ±5% prekrita z zgodovinsko produktivnostjo</p></figcaption></figure>

<figure><img src="https://29491374-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>Slika 19 Poskusna cona z N34 120kg/ha ±5% prekrita z zgodovinsko produktivnostjo</p></figcaption></figure>

Za globlje razumevanje metodologije in podrobnosti teh rezultatov so spodaj uporabljene enačbe:

1. Glavna cona s 150 kg/ha z upoštevanjem 5 % sprejemljivosti, prekrita z zgodovinsko produktivnostjo:\
   `Yield_Main = np.where((Application_Zone == 1) & (Productivity_SubZone == 51) & (Applied_Value >= 142.5) & (Applied_Value <= 157.5), Yield_WetMass, np.nan)`
2. Poskus s 120 kg/ha z upoštevanjem 5 % sprejemljivosti, prekrita z zgodovinsko produktivnostjo:\
   `Yield_Trial = np.where((Application_Zone == 3) & (Productivity_SubZone == 51) & (Applied_Value >= 114.0) & (Applied_Value <= 126.0), Yield_WetMass, np.nan)`
3. Poskus s 180 kg/ha z upoštevanjem 5 % sprejemljivosti, prekrita z zgodovinsko produktivnostjo:\
   `Yield_Trial = np.where((Application_Zone == 2) & (Productivity_SubZone == 51) & (Applied_Value >= 171.0) & (Applied_Value <= 189.0), Yield_WetMass, np.nan)`

kjer

* del `Productivity_SubZone == 51` odraža visoko produktivne cone z apliciranimi poskusi,
* deli `(Applied_Value >= 142.5) & (Applied_Value <= 157.5)` , `(Applied_Value >= 114.0) & (Applied_Value <= 126.0)`, `(Applied_Value >= 171.0) & (Applied_Value <= 189.0)` vključujejo ±5 % natančnost glede na odmerke `150`, `120`, `180` kg/ha.

## Povzetek

Rezultati pridelka iz poskusov se tesno ujemajo s povprečnim pridelkom, opaženim po celotni visoko zgodovinsko produktivni coni polja. Z drugimi besedami, poskusna aplikacija proizvoda N34 pri odmerkih <mark style="color:modra;">120 kg/ha - 150 kg/ha - 180 kg/ha</mark>, je prinesla povprečne pridelke <mark style="color:modra;">6,33 t/ha - 6,34 t/ha - 6,41 t/ha</mark> oziroma ni imela pomembnega vpliva na pobrani pridelek v visoko produktivni coni.