Voorbeeldresultaten worden gepresenteerd voor wortelsegmentatie op basis van RGB- en HS-beeldvorming. Voor spectrale beeldvorming wordt een voorbeeld van hoge-resolutie watermapping gegeven. Ten slotte worden resultaten getoond die de wetenschappelijke context demonstreren waarin beeldgebaseerde wortelgegevens worden toegepast.
RGB-gebaseerde wortelmeting
Figuur 8 toont een RGB-tijdreeks van wortelbeelden van suikerbieten cultivar Ferrara. De beelden vertonen enkele artefacten door inhomogene verlichting van de rhizoboxen, met lichtere gebieden langs de linkerkant en een verschillende helderheid in het overlappende gebied tussen de bovenste en onderste beelden.

Figuur 8: Tijdserie van wortelgroei op basis van RGB-beeldvorming.
De afbeeldingen tonen de suikerbietcultivar Ferrara op verschillende dagen na zaaien (DNS). De beelden vertonen enkele artefacten als gevolg van niet-homogene verlichting aan de linkerzijde van de afbeelding en tussen de bovenste en onderste beelden. Schaalbalken, 2 cm. Klik hier om een grotere versie van deze figuur te bekijken.
Figuur 9 geeft details over wortelsegmentatie op basis van kleurthresholding voor het ras Ferrara op 35 dagen na zaaien (DAS). Als referentie (Figuur 9A) is een binair beeld gebruikt waarbij alle wortels handmatig met een grafisch tablet zijn gevolgd. De tijd die nodig was voor het handmatig volgen van het volledige, volledig ontwikkelde, dichte wortelstelsel van de suikerbiet was ongeveer vier uur. Figuur 9B biedt een gedetailleerd beeld van een geselecteerd gebied aan de bovenkant van de afbeelding waar oude zijwortels aanwezig zijn. Hier worden verschillende wortelassen niet geclassificeerd door de kleurthreshold. Onderaan (Figuur 9C) daarentegen, waar witte jonge wortels overheersen, classificeert de op kleur gebaseerde segmentatie alle wortelassen correct. Het gebinariseerde wortelstelsel (Figuur 9D) vertoont aan de linkerkant een zwart gebied door een belichtingsartefact, dat als uitsluitingsgebied is gedefinieerd vóór het uitvoeren van de kwantitatieve analyse. Figuur 9E toont de bijbehorende pixelhistogrammen van geselecteerde kenmerken (wortels versus bodem) voor het rode kanaal van de RGB-afbeelding van Ferrara op DAS 35. De wortelpixels (blauwe kleur) vertonen duidelijk drie pieken die overeenkomen met heldere jonge zijwortels, donkere oude zijwortels en de penwortel. De overlap tussen de oude zijwortels en de bodemachtergrond is zeer sterk, wat leidt tot ongeclassificeerde wortelassen (cf. Figuur 9B).

Figuur 9: Wortelsegmentatie met behulp van een kleurthreshold.
(A) Handmatig gesegmenteerd wortelstelsel met behulp van een grafisch tablet, (B) gebied met slecht gesegmenteerde oude wortelasjes in de boven- en (C) correct gesegmenteerde jonge assen onderaan de afbeelding, en (D) binair beeld verkregen uit kleurgebaseerde drempelwaardeanalyse. (E) Pixelhistogrammen voor geselecteerde kenmerken van de RGB-afbeelding. Wortels worden weergegeven door de blauwe balken, waarbij verschillende worteltypen zijn aangegeven; bodem wordt weergegeven door de rode balken. Schaalbalken in A en D, 2 cm; schaalbalken in B en C, 1 cm. Klik hier om een grotere versie van deze figuur te bekijken.
De resulterende totale zichtbare wortellengte gekwantificeerd voor de handmatig gesegmenteerde referentieafbeelding is 1534.1 cm, terwijl de geautomatiseerde, op kleur gebaseerde segmentatie een totale wortellengte van 1427.6 cm geeft.
Grijswaardenbeelden verkregen via UV-verlichting bieden in het hier getoonde geval geen voordeel en presteerden slechter in vergelijking met kleur-thresholding (wortellengte: 1679,7 cm). Oude wortels konden niet worden gesegmenteerd en er was meer ruis in het beeld, waarschijnlijk door de lagere lichtintensiteit van de UV-lampen. Echter, in het geval van jonge wortels met een hoge autofluorescentie en een helder achtergrondsubstraat, kan UV-verlichting nog steeds een optie zijn, zoals blijkt uit een beeld verkregen uit een ander experiment waarbij zand als achtergrondsubstraat werd gebruikt (Figuur 10).

Figuur 10: UV-verlichting om wortels te visualiseren op een heldere achtergrond.
Voorbeeld van een wortelstelsel van harde tarwe dat groeit in een rhizobox gevuld met kwartszand. De rhizobox is gefotografeerd met verlichting voor (A) UV-licht en (B) fluorescerend (dag)licht. Schaalbalken, 2 cm. Klik hier om een grotere versie van deze figuur te bekijken.
HSI-gebaseerde wortelmetingen
Figuur 1 toont de gemiddelde spectra voor drie root ROI's (oude en jonge zijwortels, penwortel) en twee bodem ROI's (boven- en onderkant van de rhizobox).

Figuur 11: Gemiddelde spectra van wortel en bodem.
Spectra van interessegebieden (ROIs) op de wortel (drie worteltypen) en in de bodem (boven- en onderkant van de rhizobox). De ROI's zijn geselecteerd om een optimaal segmentatiecriterium tussen wortel en bodem te bepalen. Klik hier om een grotere versie van deze figuur te bekijken.
Het is duidelijk dat de penwortel en jonge zijwortels aanzienlijk verschillen van de achtergrond in de intensiteit van de meeste spectrale banden. Voor de oude zijwortels zijn de intensiteitsverschillen veel kleiner. Een kenmerk dat visueel kan worden afgeleid, is de verschillende helling van het spectrum rond het waterabsorptiegebied (1450 nm). Hier is de helling van de wortelspectra hoger in vergelijking met de bodemspectra. Bovendien kan een verandering in de spectra van de penwortel en jonge zijwortels worden geïdentificeerd in het gebied rond 10 nm, die niet voorkomt bij de oude zijwortels.
Figuur 12A toont het resultaat van het zoekalgoritme dat een spectrale ratio met het sterkste contrast tussen voorgrond en achtergrond identificeert. De ratio van spectra bij 1476 nm ten opzichte van 1076 nm biedt de beste scheiding tussen wortels en bodem. Het resulterende histogram van wortelvoorgrond- en bodemachtergrondpixels is weergegeven in Figuur 12B>. Hoewel er enige overlap is, zijn de meeste pixels duidelijk gescheiden van de bodemachtergrond. Door een bimodale Gaussische curve door het histogram te fitten en de afstand van Bhattacharyya voor kwantificering te gebruiken, wordt een waarde van 7,80 verkregen. Een waarde hoger dan 3,0 duidt op een sterk beeldcontrast, waardoor een betrouwbare scheiding mogelijk is28.

Figuur 12: Verschil in reflectie tussen de wortel op de voorgrond en de bodem op de achtergrond voor verschillende spectrale bandverhoudingen en pixelhistogram bij de spectrale verhouding die is gebruikt voor segmentatie.
(A) Helle kleuren (geel) tonen een hoog contrast tussen voorgrond en achtergrond, donkere kleuren (blauw) tonen een laag contrast. De eerste 15 banden zijn verwijderd vanwege ruis. De rode lijnen geven de bandratio aan met het hoogste contrast. (B) Pixelhistogram van wortels (blauw) en bodem (rood) bij de spectrale ratio voor segmentatie. Blauwe balken representeren de wortel en rode balken de bodem. De intensiteitswaarde komt overeen met de ratio van spectrale band 160 ten opzichte van spectrale band 49. Klik hier om een grotere versie van deze figuur te bekijken.
Het binaire beeld (Figuur 13) is gemaakt door een globale intensiteitsdrempel toe te passen op de geïdentificeerde spectrale ratio met een waarde van 1,08, berekend op basis van de histogramafstand27. Analyse van de wortellengte van dit beeld geeft een totale lengte van 1557,3 cm, wat een fout van slechts 1,5% vertegenwoordigt in vergelijking met het handmatig gevolgde referentiebeeld.

Figuur 13: Binair beeld van het wortelstelsel van suikerbietcultivar Ferrara.
De afbeelding is verkregen door een globale spectrale drempelwaarde toe te passen. Schaalbalk (linksonder), 2 cm. Klik hier om een grotere versie van deze figuur te bekijken.
Hoewel wortelsegmentatie is verbeterd door het gebruik van spectrale informatie vergeleken met informatie op basis van kleur, is het hoofddoel van HS-beeldvorming de analyse van chemometrische beeldeigenschappen. Dit wordt geïllustreerd door het in kaart brengen van het watergehalte van een rhizobox-beeld.
Figuur 14A toont de gemiddelde spectra van de compartimenten in de kalibratierhizobox (vgl. Figuur 7) gevuld met bodem met verschillende watergehaltes. De vorm van de spectra is vergelijkbaar tussen de compartimenten; een spectrale ratio biedt hier dus niet noodzakelijkerwijs een stabieler classificatiecriterium. Daarom is de intensiteit bij een enkele spectrale band (1680 nm), waar het gemiddelde verschil tussen aangrenzende watergehaltes maximaal is, geïdentificeerd als het beste scheidingscriterium. De resulterende pixelhistogrammen voor deze spectrale golflengte worden getoond in Figuur 14B.

Figuur 14: Spectrale kenmerken voor kalibratie van het watergehalte.
(A) Gemiddelde spectra van negen watercompartimenten uit de kalibratierhizobox met verschillende watergehalten; (B) Pixelhistogrammen voor de watercompartimenten bij band 216, waar de gemiddelde afstand tussen naburige compartimenten maximaal is. Klik hier om een grotere versie van deze figuur te bekijken.
De relatie tussen de gemiddelde pixelintensiteit bij 1680 nm en het gemeten watergehalte wordt getoond in Figuur 15.

Figuur 15: Relatie tussen spectrale reflectie en volumetrisch watergehalte.
De figuur toont dataparen van het gemeten watergehalte en de spectrale reflectie, met empirische curven (lineair en exponentieel) die zijn aangepast aan de gegevens met uitzondering van de hoogste watergehaltes (rode driehoeken). Klik hier om een grotere versie van deze figuur te bekijken.
Het onderscheid tussen hogere watergehaltes op basis van spectrale intensiteit wordt moeilijk. Een significante regressie (lineair of exponentieel) met een hoge R2 kan worden gefit voor watergehaltes tot ongeveer 0,30 cm3 cm-3. Nattere bodemcondities kunnen niet betrouwbaar worden voorspeld aan de hand van de intensiteitswaarde. Vergelijkbaar gedrag van een exponentiële relatie tussen reflectiviteit en watergehalte, met een afnemende respons bij watergehaltes hoger dan 0,30 cm3 cm-3, werd ook in andere studies gevonden30.
Een afbeelding van een rhizobox met een gedetailleerde mapping van het watergehalte is weergegeven in Figuur 16. Vier aspecten dienen te worden opgemerkt. Ten eerste is een gebied met een lager watergehalte zichtbaar in de bewortelde delen van de rhizobox. Ten tweede is de sterkste uitputting geconcentreerd in de nabijheid van individuele wortelasen. Ten derde komen uitputtingszones ook voor waar geen wortelasen aan het oppervlak zichtbaar zijn, wat wijst op gebieden waar wortels in de bodem verborgen zijn. Ten vierde resulteert watermapping zonder verdere beeldbewerking in een vlekkerig uiterlijk als gevolg van de geaggregeerde bodem. Dit kan duiden op een inhomogene watergehalteverdeling op aggregaatschaal, maar ook op het effect van de oppervlaktemorfologie op de beeldkwaliteit. Chemometrische beeldbewerkingstechnieken zijn een optie om dergelijke morfologische effecten in spectrale beelden te overwinnen31, maar zijn tot nu toe niet geïmplementeerd in de hier gebruikte Matlab-scripts.

Figuur 16: Watergehaltekaarten van een rhizobox.
De donkerblauwe kleuren vertegenwoordigen gebieden met een hoog watergehalte, groen tot rode gebieden tonen gebieden met een laag watergehalte. De plantenwortel is in het zwart over de afbeelding geplaatst. Schaalbalk (linksonder), 2 cm. Klik hier om een grotere versie van deze figuur te bekijken.
Toepassingsvoorbeelden
Figuur 17 relateert kwantitatieve worteleigenschappen uit beeldanalyse aan bovengrondse metingen.

Figuur 17: Typische toepassingsvoorbeelden voor wortelgegevens.
(A) en (B) toon wortelinformatie die wordt gebruikt voor de karakterisering van bovengrondse en ondergrondse plantonderdelen in een fenotyperingscontext. (A) stelt de wortelgroei voor van de suikerbietcultivar Ferrara, (B) vergelijkt zes in rhizoboxen gekweekte suikerbietcultivars met behulp van de blad-worteloppervlakteverhouding (gegevens van één replica). (C) en (D) zijn functionele relaties tussen kenmerken zoals aangetroffen in fysiologisch plantonderzoek. (C) toont de invloed van de blad-worteloppervlakteverhouding op de productie van droge stof en (D) de relatie tussen het worteloppervlak en de stomatale conductantie. Klik hier om een grotere versie van deze figuur te bekijken.
Figuur 17A en 17B zijn relevant voor fenotypering gericht op een uitgebreide karakterisering van de plant boven- en ondergronds. Figuur 17A toont de wortelgroei van suikerbietcultivar Ferrara (cf. Figuur 8 voor afbeeldingen). De expansie van het wortelstelsel geeft de capaciteit van een cultivar aan om het bodemvolume te exploiteren binnen een bepaalde tijdsspanne van de vegetatieperiode. Figuur 17B toont de verhouding tussen blad- en worteloppervlak van zes suikerbietcultivars, wat een descriptor vormt voor de balans tussen het aanbod van de plant (wortel) en de vraag (blad).
Figuren 17C en 17D geven voorbeelden van functionele relaties die relevant zijn in fysiologisch onderzoek. In Figuur 17C is de ratio van het bladoppervlak ten opzichte van het worteloppervlak gerelateerd aan de droge stof die tijdens het experiment is gevormd, wat wijst op de overheersende rol van het assimilerende oppervlak als beperkende factor voor de accumulatie van droge stof. Het gebrek aan significantie, ondanks een relatief hoge R2, is te wijten aan het lage aantal gepaarde gegevens (n=6) dat hier is gebruikt. Figuur 17D laat zien dat cultivars met een groter worteloppervlak (verbeterde opname) gedurende het experiment een gemiddeld hogere stomatale conductantie vertonen. Het grotere worteloppervlak ondersteunt blijkbaar de waterextractie, waardoor de opening van de huidmondjes wordt verlengd.