Détection et estimation par image en AEC — emplacement des capteurs et surface foliaire projetée

Essayer l'expérience de surface projetée dans votre navigateur

Le panneau ci-dessus exécute l'expérience d'image synthétique sans installer Python. Choisissez la scène propre, avec étiquette verte, sombre ou rognée ; modifiez la luminosité minimale V, la région d'intérêt et le mm/pixel. Neuf préréglages incluent les sept conditions publiées, une ROI sous-dimensionnée et une sélection vide. Les préréglages réinitialisent chaque entrée ; modifier le seuil manuellement change également le cas sombre récupéré.

Inspectez l'entrée, la référence dans le cadre, le masque extrait et la superposition d'erreur. Le rouge marque les faux positifs et le bleu les faux négatifs. Une bordure marron marque une ROI sélectionnée sur l'entrée. Enregistrez un résultat comme A et modifiez une entrée pour comparer images et métriques à la même échelle d'affichage. Le tableau de données donne des segments horizontaux de pixels (x/y indexés à zéro et longueur) ; ils reconstruisent exactement chaque région d'image, sans sous-échantillonnage. Partagez les conditions ou téléchargez le résultat JSON pour conserver l'expérience.

Les images restent les scènes géométriques originales de 320×240, sans téléversement de photo réelle ni entrée caméra. Les fonctions NumPy originales de géométrie et de mesure sont réutilisées sans modification. Le navigateur remplace uniquement la conversion/sélection de plage OpenCV par une recherche vérifiée de ses trois couleurs RGB fixes : fond (235,235,235), cible (50,160,65) et cible sombre (20,60,25). Leurs valeurs HSV uint8 OpenCV 4.6.0 sont (0,0,235), (64,175,160) et (64,170,60). Ce n'est pas un convertisseur HSV à usage général. H 35–85 et S 80–255 restent fixes ; le minimum V est réglable de 0 à 255, inclus. Aucune morphologie n'est appliquée.

La ROI fixe excluant l'étiquette est (50,20,230,190) ; la petite ROI est (140,60,20,30), en pixels x/y/largeur/hauteur. Aucune ne s'adapte à une plante. L'échelle varie de 0,1 à 1,0 mm/pixel comme conversion de plan de référence supposée, non un étalonnage. L'image reste identique lorsque l'échelle change. touches_frame apparaît comme 1 pour oui et 0 pour non ; une sélection vide a une précision indéfinie (—). Un IoU parfait pour une scène rognée ne dit toujours rien des portions hors de l'image. Ces sorties ne mesurent ni la biomasse, ni la croissance, ni le stress, ni la surface foliaire réelle.

Les téléchargements originaux et les références sauvegardées restent inchangés. Les instructions de reproduction locale ci-dessous sont facultatives. L'adaptateur navigateur utilise NumPy sans OpenCV ; enregistrez son point d'entrée, son noyau et sa palette ensemble pour l'exécuter localement avec NumPy installé.

leaf: experiment.py

leaf: leaf_kernel.py

leaf: palette.json

leaf: original.py

Détection et débit en AEC — où mesurer le CO₂ et la solution nutritive

Des affichages identiques peuvent renvoyer à des lieux et des instants différents

L'OD du réservoir diffère de l'OD à la sortie d'un canal de culture. Avant d'attribuer la différence à la consommation d'oxygène racinaire, vérifiez si les capteurs observent des instants comparables, si l'étalonnage est cohérent, et si l'échantillon change en transit.

Le guide de la zone racinaire séparait la CE, le pH, l'OD et le bilan hydrique. Cet article demande où et comment les comparer. Il reprend les concepts de retard et d'étalonnage du guide des capteurs de température/humidité, en ajoutant des enjeux spécifiques au CO₂, à la solution nutritive et au débit.

Le tableau des emplacements et la procédure d'investigation ci-dessous constituent notre approche de conception proposée. Ce ne sont pas des mesures d'installation, des nombres de capteurs requis ou des seuils d'alarme universels pour une culture.

Décider quelles différences vous devez observer

Utilisez des points permanents pour la surveillance et des mesures portables pour la cartographie. Conservez une référence fixe pour distinguer la variation temporelle pendant une enquête mobile. Attribuez des identifiants d'emplacement indépendants des identifiants de capteurs afin que les remplacements ne brisent pas l'historique de l'emplacement.

L'autocalibration du CO₂ correspond-elle aux conditions de culture ?

Une mesure immédiatement à côté d'une sortie de CO₂ ne devrait pas automatiquement représenter toute la canopée. Enregistrez les relations avec l'injection, l'air soufflé/d'extraction et la respiration humaine. Un point diagnostiquant l'équipement de dosage sert un objectif différent d'un point évaluant les conditions de la culture.

Vérifiez l'algorithme d'étalonnage. La note de calibration sur le terrain 2018 de Sensirion pour le SCD30 décrit l'ASC utilisant une concentration minimale correspondant à 400 ppm et nécessitant une exposition régulière à cette concentration.

Notre inférence est que les pièces enrichies en continu peuvent ne pas satisfaire cette hypothèse. Cela concerne l'algorithme documenté, et non chaque capteur de CO₂. Suivez les instructions du modèle et du micrologiciel sélectionnés et prévoyez une comparaison avec un gaz de référence ou un instrument étalonné. L'hypothèse historique de 400 ppm n'est pas une valeur d'étalonnage extérieure garantie actuellement.

Les sorties de température/HR intégrées à un appareil de CO₂ peuvent décrire son boîtier plutôt que l'air représentatif de la canopée. Distinguez les entrées de compensation environnementale des mesures utilisées pour évaluer l'environnement de la culture.

La CE, le pH et l'OD nécessitent des vérifications différentes

Le manuel des instruments multiparamètres de l'USGS relie stockage, maintenance, étalonnage, mesure et rapport. Sa portée sur les eaux naturelles n'établit pas de consignes de culture ni l'adéquation à une solution nutritive particulière.

La liste de contrôle proposée suivante applique cette approche générale ; les procédures et solutions réelles doivent correspondre à la documentation de l'instrument.

L'accord sur un étalon n'établit pas des lectures représentatives après installation. Étudiez les bulles, les dépôts, le liquide stagnant et les différences de concentration locales immédiatement après les ajouts.

La technologie de l'OD compte. Les recommandations de mesure d'YSI expliquent que les capteurs optiques ne consomment pas d'oxygène et n'ont pas d'exigence de brassage liée au principe de mesure, bien que le débit puisse améliorer le temps de réponse. Ne supposez pas des exigences de débit identiques pour les appareils électrochimiques et optiques. Même un capteur optique doit encore échantillonner un liquide représentatif de l'environnement racinaire visé.

Une ligne d'échantillonnage amène un liquide plus ancien au capteur

Une cellule de débit déportée sépare le temps d'affichage des conditions au point d'échantillonnage. Pour un volume de ligne synthétique de 0,50 L et un débit constant de 0,25 L/min, supposez un déplacement idéal sans mélange ni dispersion :

Diagram 1 · Use the button to switch views

Les deux minutes décrivent uniquement le transport idéal de la ligne. Le mélange dans la cellule, la réponse du capteur et le délai de communication sont des effets supplémentaires.

Doubler le débit à 0,50 L/min réduit ce temps idéal à une minute. Les lignes et cellules réelles se mélangent ; la température et les échanges gazeux peuvent aussi modifier l'échantillon. Décaler les horodatages de deux minutes ne reconstruit pas nécessairement les conditions d'origine.

Enregistrez séparément l'emplacement d'échantillonnage, la configuration de la ligne/cellule, le débit, le temps de mesure de l'instrument et le temps de réception. Lors du changement entre points d'échantillonnage, vérifiez aussi le report de l'échantillon précédent.

Comparer la commande de pompe, le débit, le niveau et la pression

Une pompe en marche n'est pas une mesure de débit. La livraison totale et la livraison à chaque branche sont des observations différentes. Suivez les exigences du débitmètre sélectionné concernant les sections droites, les tuyaux pleins, la plage et les bulles.

Convertissez le niveau en volume à l'aide de la géométrie du réservoir. Pour une section constante, la variation de volume est égale à la surface multipliée par la variation de niveau ; cette conversion n'est pas universelle pour les cuves à fond incliné. La mesure de niveau par pression nécessite également de prendre en compte la densité et la pression de référence.

La pression différentielle d'un filtre peut appuyer les comparaisons, mais elle dépend aussi du débit. Une hausse ne prouve pas à elle seule un colmatage lorsque les conditions de fonctionnement diffèrent. Enregistrez les altitudes de mesure, la pression relative ou absolue, et les unités.

Utiliser un résidu de bilan hydrique pour trouver les incohérences

Ici, la frontière est le réservoir seul. Q_{in} inclut chaque flux entrant et Q_{out} chaque flux sortant, y compris l'alimentation et le retour de circulation. Cela diffère de la frontière système entier de l'article précédent. L'exemple synthétique fixe d'autres termes, tels que l'évaporation, à zéro.

En utilisant les L/min et les minutes, on obtient un résidu r en L. Dans notre exemple original, l'entrée est de 1,0 L/min et la sortie de 0,2 L/min pendant cinq minutes. Le gain prévu est de 4 L, alors que le volume déduit du niveau ne passe que de 100 L à 101 L. Le résidu est −3 L.

Cela incite à enquêter sur des flux non enregistrés, un biais de débit/niveau, des différences de synchronisation ou des pertes omises. Ce n'est pas une preuve de fuite.

from math import isfinite


def transport_minutes(line_l, flow_l_min):
    if not all(isfinite(x) for x in (line_l, flow_l_min)):
        raise ValueError("finite values required")
    if line_l < 0 or flow_l_min <= 0:
        raise ValueError("nonnegative volume and positive flow required")
    return line_l / flow_l_min


def residual_l(v0, v1, intervals):
    if not all(isfinite(v) and v >= 0 for v in (v0, v1)):
        raise ValueError("invalid tank volume")
    if not intervals:
        raise ValueError("intervals required")
    predicted = 0.0
    for minutes, inflow, outflow in intervals:
        if not all(isfinite(x) for x in (minutes, inflow, outflow)):
            raise ValueError("finite interval values required")
        if minutes <= 0 or min(inflow, outflow) < 0:
            raise ValueError("invalid interval")
        predicted += (inflow - outflow) * minutes
    return (v1 - v0) - predicted


print(f"transport: {transport_minutes(0.50, 0.25):.1f} min")
print(f"residual: {residual_l(100, 101, [(5, 1.0, 0.2)]):.1f} L")

Enregistrez sous le nom measurement_check.py et exécutez python3 measurement_check.py pour obtenir 2.0 min et -3.0 L. Les débits sont constants dans chaque intervalle, et les intervalles sont supposés couvrir toute la période d'observation sans lacune ni chevauchement. Ce petit exemple n'accepte aucun horodatage, il ne peut donc pas vérifier automatiquement cette hypothèse.

Changez le volume final à 104 L et le résidu devient nul. Diviser cinq minutes en deux et trois minutes aux mêmes débits laisse le résultat inchangé. Pour des données mesurées, invalidez une intégrale incomplète plutôt que de traiter un débit manquant comme nul.

Supposons que chaque volume aux points terminaux ait une borne d'erreur de 0,5 L, que chaque débit ait une borne d'erreur de 0,1 L/min, et que la durée soit exacte. La contribution dans le pire des cas à l'erreur du résidu sur cinq minutes est 0.5+0.5+5(0.1+0.1)=2.0 L. Ces hypothèses seules ne peuvent pas expliquer un écart de 3 L. Il s'agit d'une borne illustrative, non de la précision d'un produit ni d'un intervalle de confiance statistique. Validez les bornes et les termes omis avant de diagnostiquer une défaillance.

Séparer la circulation de la ventilation

La circulation déplace et mélange l'air intérieur ; la ventilation échange l'air avec l'extérieur. Le guide des serres de NC State les traite séparément. Ne mettez pas le débit nominal d'un ventilateur de circulation dans un bilan de CO₂ comme un échange extérieur.

Pour les enquêtes de vitesse d'air de la canopée, enregistrez l'orientation de la sonde, l'emplacement, la durée de moyennage, l'éclairage et l'état des ventilateurs. Le débit à l'avant d'une étagère n'établit pas le débit à l'arrière. Conservez les différences spatiales et les fluctuations ainsi que les moyennes, et répétez l'enquête lorsque la croissance de la culture modifie les trajectoires d'écoulement. Une seule vitesse d'air ne prédit pas directement la transpiration ou une maladie.

Mettre en service la chaîne de mesure

1. Définir l'objectif et la frontière de l'observation ; aligner les identifiants d'emplacement, les unités et le temps.
2. Comparer à des étalons ou instruments de référence et conserver les résultats avant/après ajustement.
3. Étudier les différences d'instruments par des mesures colocalisées, séparément des enquêtes spatiales.
4. Comparer les journaux de fonctionnement aux observations, en distinguant transport, réponse et délais de communication.
5. Utiliser le bilan hydrique et des mesures multiples pour trouver des incohérences, puis étudier les causes concurrentes.

Ceci établit des observations pour le contrôle et l'estimation ultérieurs. L'estimation par image des plantes nécessite la même rigueur : conservez l'heure, l'éclairage et le point de vue de l'image avant d'interpréter des changements apparents comme de la croissance.

Images de plantes en AEC — mesurer la surface projetée, distinguer la surface foliaire

Plus de vert sur une photo signifie-t-il plus de feuillage ?

La même plante peut occuper une région visible différente lorsque l'orientation des feuilles, le chevauchement, l'éclairage ou la position de la caméra changent. Notre première tâche d'estimation se limite donc à la surface projetée d'une région sélectionnée dans une image. Elle n'estime pas directement la surface foliaire réelle, la masse, le rendement ou le stress.

Le guide de détection précédent séparait les différences spatiales, d'instrument et de synchronisation. Les images nécessitent aussi des enregistrements d'acquisition et de traitement. Ici, des ellipses qui se chevauchent forment des entrées synthétiques ; il n'y a ni photographies de plantes mesurées ni modèles entraînés. Des exemples contrôlés révèlent les causes d'échec avant de concevoir une évaluation sur de vraies plantes.

Définir d'abord la grandeur

Pour un masque binaire M(u,v) dont les pixels cibles valent un, le nombre de pixels est N=\sum_{u,v}M(u,v). Étant donné une échelle isotrope supposée s [mm/pixel], la surface équivalente au plan de référence est :

Notre référence synthétique contient 14 738 pixels. À 0,5 mm/pixel, cela devient 3 684,50 mm². L'échelle est fournie à titre d'illustration, non mesurée automatiquement.

L'analyse de taille de PlantCV calcule de même des propriétés à partir de régions étiquetées et prend en charge la conversion via les dimensions de pixel. Gardez l'extraction de région séparée de la conversion d'unités. Notre code compte directement les pixels et n'appelle pas PlantCV.

Garder l'acquisition et l'échelle cohérentes

Enregistrez la position et l'orientation de la caméra, la distance focale, les dimensions de l'image, l'exposition, la balance des blancs, l'éclairage, l'identifiant de la plante et l'horodatage. Si l'exposition ou la balance des blancs est automatique, enregistrez cet état. Un seuil fixe appliqué à des conditions d'image changeantes peut transformer des changements de couleur en changements apparents de surface.

Une longueur connue de 100 mm couvrant 200 pixels donne 0,5 mm/pixel à cet endroit. Placez la référence dans le plan de mesure et alignez la caméra aussi perpendiculairement que possible. Si la distorsion ou la perspective importe, utilisez l'étalonnage de caméra et une correction planaire appropriée. Des échelles horizontale et verticale différentes nécessitent Ns_xs_y.

Une référence au sol n'a pas besoin de partager l'échelle des feuilles poussant au-dessus. La correction planaire ne rend pas une canopée tridimensionnelle coplanaire. Lorsque ces hypothèses ne tiennent pas, rapportez le nombre de pixels sous acquisition fixe ou la surface équivalente au plan de référence, plutôt que de revendiquer une surface foliaire physique.

Utiliser 0,6 au lieu de 0,5 mm/pixel multiplie la surface du même masque par 1,44 : une augmentation de 44 %. Une extraction correcte peut coexister avec une conversion d'unités incorrecte.

Sélectionner les candidats de couleur

Le tutoriel OpenCV démontre la conversion en HSV et la sélection de plage de couleur avec inRange. Sélectionner le vert n'équivaut pas à reconnaître une feuille.

Les entrées font 320×240 pixels. Dans la représentation HSV 8 bits d'OpenCV, cet exemple sélectionne H 35–85, S 80–255 et V 100–255. H n'est pas saisi directement comme un angle de 0–360. Nos tableaux générés sont en RGB, donc la conversion utilise COLOR_RGB2HSV ; ne réutilisez pas aveuglément cette conversion sur les images BGR renvoyées par imread.

Le cœur tient en trois lignes. La génération d'image, la comparaison de conditions et la vérification se trouvent dans le téléchargement complet.

hsv = cv2.cvtColor(rgb, cv2.COLOR_RGB2HSV)
mask = cv2.inRange(hsv, (35, 80, min_v), (85, 255, 255)) != 0
pixels = int(np.count_nonzero(mask))

Aucun remplissage morphologique ni suppression de petits composants n'est appliqué. Si vous ajoutez un tel traitement, vérifiez s'il modifie les feuilles étroites ou les espaces entre feuilles par rapport aux masques de référence.

Comparer les échecs avec la même géométrie

La référence générée est l'union visible dans le cadre de trois ellipses. Le chevauchement est compté une seule fois, et les portions hors du cadre sont exclues.

Le rouge marque le fond sélectionné comme cible ; le bleu marque les pixels cibles manqués. Les régions vertes conservent leurs couleurs d'entrée. Ce sont des formes synthétiques, pas des photographies de plantes.

Le green_tag_roi fixe commence à (50,20), avec une largeur de 230 et une hauteur de 190 pixels. Il a été choisi pour éviter cette étiquette synthétique, et non comme une région d'intérêt indépendante de la culture. Une croissance au-delà de la ROI exclurait de véritables pixels cibles.

Abaisser le V minimum n'est pas une solution universelle. Cela récupère la cible sombre ici, mais peut admettre un fond vert sombre dans des scènes réelles. N'ajustez pas les seuils en fonction des réponses de chaque image et ne rapportez pas la précision uniquement sur ces mêmes images.

Évaluer surface et localisation séparément

Pour une prédiction P et une référence G, l'intersection sur union mesure le chevauchement spatial :

Un nombre égal de pixels en trop et manqués produit une erreur de surface nulle malgré des emplacements incorrects. Une référence [0,1,1,0] et une prédiction [1,1,0,0] contiennent chacune deux pixels sélectionnés, mais leur IoU est de 1/3.

Les sorties incluent aussi la précision, la fraction correcte des pixels sélectionnés, et le rappel, la fraction récupérée des pixels de référence. mask_area_error_percent est l'erreur de nombre de pixels par rapport à la référence, non une erreur physique due à l'étalonnage d'échelle. Elle reste nulle pour wrong_scale. Les ratios indéfinis, y compris l'IoU de deux masques vides, sont null en JSON plutôt qu'un score parfait supposé.

touches_frame signale les pixels sélectionnés touchant un bord d'image. C'est vrai pour clipped. Cela ne détecte pas chaque objet incomplet : la segmentation pourrait manquer les pixels de bord, ou une ROI interne pourrait retirer une partie de la cible. Inspectez l'image et la ROI même lorsque l'indicateur est faux.

Télécharger et reproduire

Téléchargez le ZIP complet, la source Python, les résultats attendus, le README et le manifeste. Le ZIP inclut les entrées, les références, les masques extraits et les superpositions d'erreur, tous régénérables à partir de la même source.

Environnement testé : Linux x86_64, Python 3.12.3, NumPy 1.26.4 et OpenCV 4.6.0. Les versions de dépendances sont figées pour la reproduction, non pour une certification d'autres systèmes d'exploitation ou versions de Python.

Après avoir extrait le ZIP, exécutez :

python3 -m venv .venv
.venv/bin/python -m pip install -r requirements.txt
.venv/bin/python experiment.py --self-test
.venv/bin/python experiment.py --verify expected.json --out reproduced

Sous Windows, remplacez .venv/bin/python par .venv\Scripts\python.exe. Le répertoire de sortie ne doit pas exister ; les sorties existantes ne sont pas écrasées. Le JSON de résultats et les fichiers PNG de chaque condition sont enregistrés. Les sorties de référence sont incluses sous samples/.

Modifiez le V minimum avec :

.venv/bin/python experiment.py --min-v 40 --out lower-v

Seul dim_lower_v utilise toujours 40 ; les autres conditions utilisent la valeur fournie. Les résultats modifiés échouent intentionnellement à la comparaison avec la référence expected.json. Ne vous attendez pas à ce que la vérification de référence réussisse après avoir modifié l'expérience.

Sept tests couvrent un rectangle connu, des masques de surface égale/emplacement différent, des masques vides, les relations entre sept cas, des entrées invalides, des résultats attendus modifiés et des omissions de ROI. Le téléchargement original n'a pas de fonctionnalité d'ingestion de photo réelle ; l'adaptateur navigateur est disponible ci-dessus.

Évaluer sur des plantes mesurées avant d'affirmer des conclusions biologiques

Ce qui suit est une proposition d'évaluation de l'étape suivante, non un travail accompli par cet exemple synthétique.

1. Conserver les identifiants de plante, les dates, l'illumination, la géométrie de la caméra et les images originales.
2. Définir l'appartenance à la cible, y compris si les tiges ou les feuilles sénescentes comptent, et créer des masques de référence humains.
3. Séparer les images de développement du seuil des images d'évaluation. Envisager des répartitions par plante ou par date d'acquisition plutôt que de simplement séparer des images vidéo adjacentes.
4. Rapporter les erreurs par éclairage, fond, chevauchement et rognage, en conservant les métriques de surface, l'IoU et les superpositions d'erreur.
5. Pour estimer la surface foliaire réelle ou la masse sèche, recueillir des mesures de référence indépendantes correspondantes et valider cette correspondance comme un modèle distinct.

Les changements de surface projetée dans une série temporelle peuvent résulter de la posture ou du chevauchement. Produire une métrique d'image n'établit pas un taux de croissance ou un stress. Rendez d'abord l'observation reproductible, puis utilisez des données de référence pour tester ce qu'elle représente.

Zone racinaire en AEC — séparer CE, pH, oxygène et bilan hydrique

Lorsque le niveau d'eau baisse, que faut-il remplacer ?

Une solution nutritive passe de 100 L à 90 L. Cette seule observation ne peut pas distinguer évaporation, absorption d'eau par la plante, fuite ou drainage. Elle n'établit pas non plus l'absorption de nutriments. Restaurer le volume et restaurer la composition en nutriments sont des décisions distinctes.

Faisant suite au guide de l'environnement lumineux, cet article examine l'environnement racinaire. Il introduit un système à solution recirculante ; ce n'est pas une recette de culture ni une implémentation de dosage automatique. Tous les exemples numériques sont synthétiques.

La zone racinaire est la région où les racines échangent eau, nutriments et gaz avec leur environnement. Une lecture du réservoir est utile, mais ne décrit pas nécessairement la surface des racines, l'extrémité d'un canal ou l'intérieur d'un substrat.

Garder cinq mesures distinctes

Ceci est un tableau de planification de mesure. Les conditions acceptables dépendent de la culture, du stade de croissance et du système de production.

La CE ne peut pas reconstituer la composition

1 mS/cm = 1 000 µS/cm = 1 dS/m. Enregistrez les unités. La conductivité change aussi avec la température du liquide : conservez la température mesurée, les réglages de compensation et la température de référence, et distinguez les lectures compensées des lectures non compensées.

La CE ne révèle pas les concentrations de sels individuels, comme l'explique le guide d'analyse d'eau de l'Oklahoma State University. Une seule lecture ne peut pas distinguer l'épuisement d'un nutriment de l'accumulation d'un autre ion.

Un affichage « ppm » dérivé de la CE n'est pas non plus une mesure de nutriments individuels. Enregistrez son facteur de conversion et distinguez-le d'un résultat de concentration de laboratoire. Utilisez l'analyse de l'eau source et de la solution lorsque la composition est nécessaire.

Le guide de gestion des nutriments de l'OSU identifie le test de l'eau source, la concentration en sels, le pH, l'alcalinité et les ratios de nutriments comme pertinents. Notre inférence opérationnelle est d'éviter de traiter un retour à la consigne de CE comme une preuve que la formulation originale a été restaurée.

Enregistrer pH et alcalinité séparément

Le pH décrit l'acidité ou la basicité actuelle. Parce qu'il est logarithmique, une moyenne pondérée par le volume de deux lectures de pH ne donne généralement pas le pH de la solution mélangée. USGS : le pH et l'eau

L'alcalinité mesure la capacité de neutralisation des acides. Deux eaux ayant le même pH initial peuvent réagir différemment au même ajout d'acide lorsque leur alcalinité diffère. USGS : l'alcalinité et l'eau

Une règle qui ajoute toujours la même dose pour une erreur de pH de 0,5 unité confond ces grandeurs. Pour évaluer une intervention, enregistrez l'eau source, le volume du réservoir, la substance et la quantité ajoutée, et l'heure de mesure après mélange. Le bilan hydrique ci-dessous ne peut pas calculer la dose requise.

L'oxygène nécessite une voie d'apport

La solubilité de l'oxygène à l'équilibre dépend de la température de l'eau, de la pression atmosphérique et de la salinité. USGS DOTABLES utilise ces conditions pour calculer la solubilité. À pression et salinité égales, une eau plus chaude retient généralement moins d'oxygène dissous à l'équilibre. Vérifiez les conditions de conversion lors de la comparaison entre mg/L et pourcentage de saturation.

Les racines ont besoin d'oxygène pour respirer. Le guide de Cornell sur les légumes-feuilles hydroponiques traite de l'apport d'oxygène aux racines et de la température de l'eau. Une lecture d'OD du réservoir apparemment adéquate n'établit pas les mêmes conditions à l'intérieur d'une masse racinaire dense. Enregistrez la circulation, l'aération et l'emplacement d'échantillonnage, y compris les différences entre l'entrée et la sortie du canal de culture.

Diagram 2 · Use the button to switch views

Pour une frontière englobant tout le système, la circulation interne n'est pas un apport d'eau externe. Le schéma n'implique pas un mélange complet ni un oxygène uniforme.

Avec un substrat, la rétention d'eau et les voies d'air comptent aussi. L'OD du réservoir seul n'établit pas l'aération du substrat. N'appliquez pas directement le calcul du compartiment liquide unique ci-dessous à l'intérieur d'un substrat.

Garder des comptes distincts pour l'eau et les solutés

Choisissez une frontière englobant tout le système recirculant, y compris le réservoir, les tuyaux et les canaux de culture. Sur un intervalle, soit V_{in} les ajouts externes, V_{out} le rejet liquide, et V_w l'évaporation plus l'eau transférée aux plantes, tous en L.

N'ajoutez pas la circulation interne de la pompe comme eau d'appoint. Une frontière limitée au réservoir seul compterait au contraire explicitement ses flux sortants et retournants.

Pour un soluté suivi, la masse est M [mg] et la concentration C=M/V [mg/L] :

U est le transfert vers les plantes ; R est le retrait net de la phase liquide, comme la précipitation. Le retrait du liquide ne signifie pas destruction de matière. Le rejet liquide emporte à la fois eau et soluté, contrairement à l'évaporation. Pour un rejet continu, utilisez M_{out}=\int C_{out}(t)Q_{out}(t)dt pour tenir compte de la concentration changeante.

L'exemple arithmétique original suivant vérifie cette comptabilité. Ce n'est pas un taux d'absorption mesuré par une plante, une recette d'engrais ou une prédiction de CE. Commencez avec 100 L contenant 10 000 mg d'un composant : 100 mg/L. Supposez aucun ajout ni précipitation.

La même baisse de volume peut s'accompagner d'une concentration croissante, décroissante ou inchangée. Ajouter 10 L d'eau ne contenant aucun soluté suivi ramène A à 100 mg/L, mais B n'atteint que 80 mg/L. Restaurer le volume ne restaure pas la masse de soluté.

Reproduire la comparaison en Python

Enregistrez le code suivant sous le nom root_zone.py et exécutez python3 root_zone.py. Seule la bibliothèque standard de Python est nécessaire.

from math import isfinite


def concentration(volume_l, mass_mg):
    if not (isfinite(volume_l) and isfinite(mass_mg)):
        raise ValueError("finite values required")
    if volume_l <= 0 or mass_mg < 0:
        raise ValueError("positive volume and nonnegative mass required")
    return mass_mg / volume_l


initial_v, initial_m = 100.0, 10000.0
initial_c = concentration(initial_v, initial_m)
cases = {
    "A": (initial_v - 10, initial_m),
    "B": (initial_v - 10, initial_m - 2000),
    "C": (initial_v - 10, initial_m - 10 * initial_c),
}
for name, (volume, mass) in cases.items():
    print(f"{name}: {concentration(volume, mass):.2f} mg/L")
print(f"B + water: {concentration(cases['B'][0] + 10, cases['B'][1]):.2f} mg/L")

Les sorties sont 111.11, 88.89, 100.00 et 80.00 mg/L. Changez le retrait de 2 000 mg dans B à 1 000 mg : après la perte de volume, la concentration devient 100 mg/L. L'absorption peut se produire même lorsque la concentration est inchangée. Un volume nul et une masse restante négative déclenchent des erreurs.

Ce code ne calcule pas la CE, le pH ni l'OD. Une comptabilité correcte pour un composant ne valide pas l'équilibre chimique multicomposant ni la réponse de la culture.

Décider quoi observer ensuite

Ce sont des combinaisons d'observation proposées, non des règles de dosage automatique.

Associez temps, emplacement, unités et étalonnage au volume d'eau, à l'appoint, au rejet et aux ajouts. L'emplacement détaillé des capteurs, la gestion individuelle des nutriments et le contrôle du dosage sont des étapes ultérieures. La base ici est de distinguer une concentration observée de la masse restant dans le système.