Jesienią 2013 roku przeprowadzono eksperyment polowy w celu porównania rozwoju rozłogów i kłączy czterech odmian trawy bermudzkiej, w tym dwóch typów wysiewnych („LaPaloma” i „Yukon”) oraz dwóch sterylnych hybryd wegetatywnych („Patriot” i „Tifway”). Zastosowano układ bloków randomizowanych z trzema powtórzeniami, co łącznie dało 12 poletek (2 x 2 m).
Z każdej odmiany typu trawnikowego oraz z dzikiej bermudagrass randomly zebrano czternaście rozłogów i czternaście kłączy z poletek, a także z roślin dzikiej bermudagrass rosnących w pobliżu poletek, co dało łącznie 70 rozłogów i 70 kłączy. Wszystkie rozłogi i kłącza zostały oczyszczone zgodnie z opisem w protokole (krok 2) przed dalszymi pomiarami. Średnicę i długość międzywęźli zmierzono odpowiednio za pomocą suwmiarki i linijki, a dla każdego rozłogu lub kłącza policzono liczbę węzłów. Zapisano również czas niezbędny do oczyszczenia i pomiaru próbek rozłogów i kłączy za pomocą linijki i suwmiarki. Średnice rozłogów i kłączy obliczono jako średnią wszystkich zmierzonych średnic międzywęźli. Całkowite długości rozłogów i kłączy obliczono jako sumę wszystkich długości międzywęźli. Ponadto, za pomocą systemu cyfrowej analizy obrazu, zgodnie z opisem w krokach 3 i 5, zmierzono całkowite długości skanowane oraz średnice skanowane każdego rozłogu i kłącza. Zapisano czas niezbędny do pomiaru cech rozłogów i kłączy za pomocą systemu cyfrowej analizy. Następnie każdy rozłóg i kłącze przecięto nożyczkami, aby oddzielić międzywęźla od węzłów, a międzywęźla wykorzystano do oszacowania skanowanej średnicy międzywęźla, zgodnie z opisem w krokach 3 i 5. Współczynniki korelacji Pearsona obliczono dla rozłogów i kłączy (n = 70 rozłogów, n = 70 kłączy) pomiędzy długością mierzoną a skanowaną, średnicą mierzoną a skanowaną, liczbą węzłów a wartością bezwzględną różnicy między średnicą mierzoną a skanowaną oraz średnicami mierzonymi a skanowanymi średnicami międzywęźli. Długości zmierzone linijką posłużyły do kalibracji długości oszacowanych za pomocą systemu cyfrowej analizy obrazu.
Analiza regresji wykazała wysoką korelację między skanowaną a zmierzoną długością rozłogów (Rycina 1a), ze współczynnikiem nachylenia 1,03 i punktem przecięcia -4,22, a także między skanowaną a zmierzoną długością kłączy (Rycina 1b), ze współczynnikiem nachylenia 1,03 i punktem przecięcia 4,22. Ręczne oczyszczanie 14 rozłogów i 14 kłączy zajęło średnio odpowiednio 21 min i 24 s oraz 11 min i 12 s. Średni czas pomiaru długości i średnicy za pomocą linijki i suwmiarki wyniósł 14 min i 6 s dla rozłogów oraz 13 min i 35 s dla kłączy. Skanowanie i analiza programowa próbek przy użyciu WinRHIZO zajęły średnio 11 min dla rozłogów i 12 min i 4 s dla kłączy.
Zmierzone i zeskanowane średnice były również istotnie skorelowane zarówno w przypadku rozłogów, jak i kłączy. Zależności między średnicą zmierzoną a zeskanowaną były bliskie stosunku 1:1, co wskazuje na dobre dopasowanie danych (Rycina 2a i 2b). Jednak wartość punktu przecięcia sugerowała, że system analizy obrazu cyfrowego zawyżał zmierzoną średnicę, szczególnie przy niższych wartościach, natomiast wyższe wartości średnicy kłączy były niedoszacowane. To zawyżenie może wynikać z obecności węzłów rozłogów, które są skanowane przez oprogramowanie, wpływając na całkowitą powierzchnię rzutu wykorzystywaną do obliczania średnicy (stosunek całkowitej powierzchni rzutu do całkowitej długości), podczas gdy są one pomijane podczas pomiarów suwmiarką. Korelacja między liczbą węzłów a różnicą między wartościami średnicy uzyskanymi obiema metodami (pomiarową i skanowaną) była istotna tylko w przypadku rozłogów (Rycina 3a); ponadto zmienność liczby węzłów wyjaśniała jedynie niewielką część zmienności tej różnicy (R2 = 14%). Istotna korelacja stwierdzona między zeskanowaną średnicą międzywęźla a średnicą zmierzoną (nachylenia wynoszące odpowiednio 1,01 dla rozłogów i 0,98 dla kłączy; punkty przecięcia bliskie zeru) (Rycina 4a i 4b) wykazuje, że średnicę międzywęźla można dokładnie oszacować za pomocą systemu analizy obrazu cyfrowego, pod warunkiem usunięcia węzłów. W związku z tym całkowitą długość rozłogów oraz średnią średnicę próbek składających się z licznych rozłogów lub kłączy można w łatwy i dokładny sposób określić ilościowo za pomocą systemu analizy obrazu cyfrowego.
W ramach trwającego eksperymentu w każdej polecenie sezonowo pobierano jedną próbkę darni (20 x 20 x 15 cm głębokości) od jesieni 2013 r. do lata 2015 r., a następnie postępowano z nią zgodnie z opisem w protokole. Długość rozłogów i kłączy na jednostkę powierzchni (gęstość długości) oraz ich masa na jednostkę powierzchni (gęstość masy) dla próbek pobranych w lipcu 2014 r. są przedstawione na Rysunku 5. Różnice w gęstości długości rozłogów zaobserwowano pomiędzy odmianami rozmnażanymi wegetatywnie („Patriot” i „Tifway”) a tymi rozmnażanymi z nasion („La Paloma” i „Yukon”). Odmiana „Patriot” wykazała najwyższą gęstość długości kłączy, a następnie „Tifway” oraz odmiany rozsiewne. Gęstość masy rozłogów różniła się dla wszystkich odmian, przy czym najwyższą wartość odnotowano dla „Patriot”, a następnie dla „Tifway”, „La Paloma” i „Yukon”. Odmiany rozmnażane wegetatywnie wykazały również wyższe gęstości masy kłączy niż odmiany rozsiewne. Zmiany gęstości długości i gęstości masy rozłogów oraz kłączy odmiany Patriot w całym okresie badania przedstawiono na Rysunku 6. Gęstość długości rozłogów wzrosła od marca 2014 r. do lipca 2014 r., a następnie nie zmieniała się od lipca 2014 r. do lipca 2015 r. W próbkach pobranych w październiku 2013 r. i marcu 2014 r. znaleziono jedynie kilka kłączy. Gęstość długości kłączy wzrosła w lipcu 2104 r., osiągając najwyższe wartości, lecz ponownie spadła w październiku 2014 r. Gęstość masy rozłogów nieznacznie wzrosła od marca do lipca 2014 r.; jednak szybszy wzrost zaobserwowano od lipca do października 2014 r., z późniejszym spadkiem w marcu 2015 r. Gęstość masy kłączy wykazywała trend podobny do gęstości długości kłączy, osiągając najwyższą wartość w lipcu 2014 r.
Oprogramowanie uwzględnia w analizie wszystkie obiekty znajdujące się na zeskanowanym obrazie. Przedstawiono przykład układu analizy obrazu cyfrowego z oprogramowania WinRHIZO (Rycina 7), w którym linie o różnych kolorach nakładają się na obiekty (stolony) o różnych średnicach w celu obliczenia całkowitej długości dla poszczególnych klas średnic. Można zauważyć, że analiza bierze pod uwagę fragmenty korzeni lub liści. Jak opisano w kroku 3.9, istnieje możliwość ograniczenia szerokości i liczby analizowanych klas średnic. Histogram przedstawia rozkład długości w wybranych klasach średnic (Rycina 7). Histogram ten może być wykorzystany do oceny minimalnych klas średnic, które należy wykluczyć. Wizualna obserwacja tego wykresu w górnej części obrazu ekranu wskazuje, że długość ma rozkład normalny wokół średniej klasy średnicy, z wyjątkiem dwóch pierwszych klas, które wykazują wartości wyższe niż te pasujące do rozkładu normalnego. Nawet jeśli próbki zostały starannie oczyszczone, uwzględnienie tych mniejszych klas w analizie danych może wpłynąć na wyniki, prowadząc do zawyżenia gęstości długości i zaniżenia średniej średnicy. Nasze wyniki pokazują, że długość mniejszych klas (średnica < 0.2 mm) stanowiła 13-32% całkowitych wartości długości kłącza wynikających z analizy programowej (Tabela 1). Ponadto, średnia średnica została zaniżona o 2-17% (Tabela 1).

Rysunek 1: Analiza regresji wartości długości zmierzonych linijką w stosunku do wartości oszacowanych za pomocą cyfrowego systemu analizy obrazu dla stolonów24 (a) oraz kłączy (b) trawy bermuda. Linia przerywana reprezentuje stosunek 1:1. Panel A został zmodyfikowany na podstawie pracy Pornaro et al.24. Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.

Rysunek 2: Analiza regresji wartości średnic zmierzonych suwmiarką w stosunku do wartości oszacowanych za pomocą systemu cyfrowej analizy obrazu dla stolonów trawy bermudzkiej24 (a) oraz kłączy (b). Linia przerywana reprezentuje stosunek 1:1. Panel A został zmodyfikowany na podstawie pracy Pornaro et al.24. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 3: Analiza regresji liczby węzłów rozłogów trawy bermudzkiej24 (a) i kłączy (b) w stosunku do wartości bezwzględnych różnicy między średnicą oszacowaną za pomocą systemu cyfrowej analizy obrazu a zmierzoną suwmiarką. Linia przerywana reprezentuje stosunek 1:1. Panel A został zmodyfikowany na podstawie pracy Pornaro et al.24. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 4: Analiza regresji wartości średnicy zmierzonych za pomocą suwmiarki w stosunku do wartości oszacowanych za pomocą cyfrowego systemu analizy obrazu dla stolonów trawy bermudzkiej24 (a) i kłączy (b) wyłącznie dla stawek. Linia przerywana reprezentuje stosunek 1:1. Panel A został zmodyfikowany na podstawie pracy Pornaro et al.24. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 5: Przykładowe wyniki gęstości długości i masy rozłogów oraz kłączy z badania polowego porównującego cztery odmiany trawy bermudzkiej (Patriot, Tifway, La Paloma, Yukon). Gęstość długości rozłogów (a), gęstość długości kłączy (b), gęstość masy rozłogów (c) oraz gęstość masy kłączy (d). Pionowe słupki reprezentują błędy standardowe dla sześciu powtórzeń. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 6: Przykładowe wyniki gęstości długości i masy rozłogów oraz kłączy z badania polowego, przedstawiające rozwój rozłogów i kłączy odmiany trawy bermudzkiej Patriot. Gęstość długości rozłogów (a), gęstość długości kłączy (b), gęstość masy rozłogów (c) oraz gęstość masy kłączy (d). Słupki pionowe reprezentują błędy standardowe dla sześciu powtórzeń. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 7: Przykładowy układ analizy obrazu cyfrowego z oprogramowania WinRHIZO. Zeskanowany obraz na pierwszym planie oraz wykresy słupkowe w górnej części obrazu ekranu przedstawiają rozkład długości w wybranych klasach średnic. Kolorowe linie wskazują analizę obrazu, a każdy kolor odpowiada kolorom klas średnic przedstawionych na wykresach słupkowych. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
| Odmiana uprawna | Blokada | Długość kłącza (cm/dm)2) | Średnia średnica (mm) |
| < 0,2 mm | całkowity | stosuneka | < 0,2 mm | łącznie | stosunekb |
| Patriota | 1 | 231 | 278 | 16.9 | 1.637846 | 1.5994 | 97.7 |
| Patriota | 2 | 304 | 349 | 12.8 | 1.620667 | 1.588371 | 98.0 |
| Patriota | 3 | 304 | 366 | 16.8 | 1.649918 | 1.621367 | 98.3 |
| Tifway | 1 | 184 | 231 | 20.6 | 2.149745 | 1.9951 | 92.8 |
| Tifway | 2 | 155 | 193 | 19.9 | 1.866253 | 1.76605 | 94.6 |
| Tifway | 3 | 119 | 150 | 20.9 | 1.877386 | 1.75865 | 93.7 |
| La Paloma | 1 | 17 | 23 | 24.4 | 2.139019 | 1.8904 | 88.4 |
| La Paloma | 2 | 26 | 38 | 31.6 | 2.101385 | 1.7455 | 83.1 |
| La Paloma | 3 | 34 | 47 | 27.5 | 2.033729 | 1.7354 | 85.3 |
| Yukon | 1 | 32 | 44 | 28.0 | 1.700155 | 1.4945 | 87.9 |
| Yukon | 2 | 17 | 25 | 33.2 | 1.68339 | 1.4284 | 84.9 |
| Yukon | 3 | 67 | 87 | 23.6 | 1.844721 | 1.6774 | 90.9 |
| a Długość klas ≤0,2 mm/długość całkowita | | | | |
| b Całkowita średnica/średnica klas ≤0,2 mm | | | |
Tabela 1: Gęstość długości kłączy oraz średnia średnica kłączy z uwzględnieniem i bez uwzględnienia klas o mniejszej średnicy. Gęstość długości z uwzględnieniem i bez uwzględnienia klas średnicy mniejszych niż 0,2 mm oraz ich stosunki (długość klas ≤ 0,2 mm/długość całkowita); oraz średnia średnica z uwzględnieniem i bez uwzględnienia klas średnicy mniejszych niż 0,2 mm oraz ich stosunki (średnica z uwzględnieniem klas < 0,2/średnica bez uwzględnienia klas < 0,2 mm).