KR20240007390A - A continuous leaf temperature prediction method using smart farm environmental data. - Google Patents

A continuous leaf temperature prediction method using smart farm environmental data. Download PDF

Info

Publication number
KR20240007390A
KR20240007390A KR1020220084210A KR20220084210A KR20240007390A KR 20240007390 A KR20240007390 A KR 20240007390A KR 1020220084210 A KR1020220084210 A KR 1020220084210A KR 20220084210 A KR20220084210 A KR 20220084210A KR 20240007390 A KR20240007390 A KR 20240007390A
Authority
KR
South Korea
Prior art keywords
temperature
smart farm
leaf
leaf temperature
measuring
Prior art date
Application number
KR1020220084210A
Other languages
Korean (ko)
Inventor
김준용
Original Assignee
주식회사 파모스
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 주식회사 파모스 filed Critical 주식회사 파모스
Priority to KR1020220084210A priority Critical patent/KR20240007390A/en
Publication of KR20240007390A publication Critical patent/KR20240007390A/en

Links

Classifications

    • AHUMAN NECESSITIES
    • A01AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
    • A01GHORTICULTURE; CULTIVATION OF VEGETABLES, FLOWERS, RICE, FRUIT, VINES, HOPS OR SEAWEED; FORESTRY; WATERING
    • A01G9/00Cultivation in receptacles, forcing-frames or greenhouses; Edging for beds, lawn or the like
    • A01G9/24Devices or systems for heating, ventilating, regulating temperature, illuminating, or watering, in greenhouses, forcing-frames, or the like
    • AHUMAN NECESSITIES
    • A01AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
    • A01GHORTICULTURE; CULTIVATION OF VEGETABLES, FLOWERS, RICE, FRUIT, VINES, HOPS OR SEAWEED; FORESTRY; WATERING
    • A01G7/00Botany in general
    • A01G7/06Treatment of growing trees or plants, e.g. for preventing decay of wood, for tingeing flowers or wood, for prolonging the life of plants
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Mining
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/10Services
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16YINFORMATION AND COMMUNICATION TECHNOLOGY SPECIALLY ADAPTED FOR THE INTERNET OF THINGS [IoT]
    • G16Y10/00Economic sectors
    • G16Y10/05Agriculture
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16YINFORMATION AND COMMUNICATION TECHNOLOGY SPECIALLY ADAPTED FOR THE INTERNET OF THINGS [IoT]
    • G16Y20/00Information sensed or collected by the things
    • G16Y20/10Information sensed or collected by the things relating to the environment, e.g. temperature; relating to location
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16YINFORMATION AND COMMUNICATION TECHNOLOGY SPECIALLY ADAPTED FOR THE INTERNET OF THINGS [IoT]
    • G16Y20/00Information sensed or collected by the things
    • G16Y20/20Information sensed or collected by the things relating to the thing itself
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16YINFORMATION AND COMMUNICATION TECHNOLOGY SPECIALLY ADAPTED FOR THE INTERNET OF THINGS [IoT]
    • G16Y40/00IoT characterised by the purpose of the information processing
    • G16Y40/10Detection; Monitoring
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16YINFORMATION AND COMMUNICATION TECHNOLOGY SPECIALLY ADAPTED FOR THE INTERNET OF THINGS [IoT]
    • G16Y40/00IoT characterised by the purpose of the information processing
    • G16Y40/20Analytics; Diagnosis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/02Indexing codes associated with the analysed material
    • G01N2291/028Material parameters
    • G01N2291/02845Humidity, wetness
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/02Indexing codes associated with the analysed material
    • G01N2291/028Material parameters
    • G01N2291/02881Temperature
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A40/00Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
    • Y02A40/10Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in agriculture
    • Y02A40/25Greenhouse technology, e.g. cooling systems therefor

Abstract

스마트팜 등에서 식물의 생육 상태를 측정하기 위하여 많은 센서들이 사용되고 있다. 특히, 식물의 현재 상태를 나타내는 지표로 식물 잎의 온도 즉, 엽온을 사용하고 있다. 스마트팜 내부의 엽온과 스마트팜 내부의 공기 온도의 비교를 통해 작물(식물)의 증산작용이 활발한지 여부를 측정할 수 있다. 또한 엽온을 이용하여 식물의 고온스트레스 발생 진단이 가능하다.
이렇게 작물(식물)의 생육상태를 반영하는 엽온을 측정하기 위해서는 열화상카메라나 적외선 온도센서가 추가 적으로 필요하며, 잎 끝에 맞춰 적외선센서를 설치하거나 열화상카메라로 촬영한 영상에서 잎의 이미지만을 구분하는 추가 작업이 필요하여, 엽온을 측정하는 것이 쉽지않는 문제가 있다. 상기와 같은 문제를 해결하기 위하여, 스마트팜 내부의 온도(intemp.)를 측정하는 내부온도측정단계; 및 스마트팜 외부의 온도(outtemp.)를 측정하는 외부온도측정단계; 및 스마트팜 내부로 입사되는 태양광의 일사량(rad.)을 측정하는 일사량측정단계; 및 엽온(canopy temp.)은 intemp.×A + outtemp.×B + rad.×C + D(식(1))로 계산되는 것을 특징으로 하는 스마트팜 내부환경 데이터를 이용한 엽온 예측방법을 제공한다. (상기 A, B, C, D는 각각 곱해지는 측정값에 대한 엽온의 증가에 기여하는 정도에 따른 상수이다)
본 출원 발명의 상기와 같은 구성에 의하여 스마트팜 내부에서 재배하는 식물의 엽온을 직접측정하지 않고도 식물의 엽온을 스마트팜 내부와 외부에 구비된 환경 변수 값으로부터 예측할 수 있는 수단을 제공하는 효과가 있다.
Many sensors are used to measure the growth status of plants in smart farms, etc. In particular, the temperature of plant leaves, that is, leaf temperature, is used as an indicator of the current state of the plant. By comparing the leaf temperature inside the smart farm and the air temperature inside the smart farm, it is possible to measure whether the transpiration of crops (plants) is active. Additionally, it is possible to diagnose high-temperature stress in plants using leaf temperature.
In order to measure leaf temperature, which reflects the growth state of crops (plants), an additional thermal imaging camera or infrared temperature sensor is required. An infrared sensor is installed at the tip of the leaf or only the image of the leaf is captured from the image taken with a thermal imaging camera. There is a problem in that it is not easy to measure leaf temperature because additional work is needed to differentiate. In order to solve the above problem, an internal temperature measurement step of measuring the temperature (intemp.) inside the smart farm; And an external temperature measurement step of measuring the temperature (outtemp.) outside the smart farm; And a solar radiation measurement step of measuring the solar radiation amount (rad.) of sunlight entering the smart farm; And canopy temperature (canopy temp.) is calculated as intemp. × A + outtemp. × B + rad. . (A, B, C, and D are constants depending on the degree to which they contribute to the increase in leaf temperature relative to the measured value being multiplied)
The above-mentioned configuration of the invention of this application has the effect of providing a means of predicting the leaf temperature of plants from the environmental variable values provided inside and outside the smart farm without directly measuring the leaf temperature of the plants grown inside the smart farm. .

Description

스마트팜 내부 환경데이터를 이용한 엽온 예측방법{.}Leaf temperature prediction method using smart farm internal environmental data{.}

본 출원 발명은 스마트팜에서 재배되는 식물의 엽온을 예측하는 기술에 관한 것이다. 더욱 자세하게는 직접 식물의 엽온을 측정하지 않고 정확히 예측하는 기술에 관한 것이다.The present application relates to technology for predicting leaf temperature of plants grown in smart farms. More specifically, it concerns technology that accurately predicts plant leaf temperature without directly measuring it.

본 발명의 출원 이전의 선행기술로 식물의 온도를 측정하여 식물생장 환경데이터를 수집하는 장치에 관한 기술이 개시되어 있다. 이 기술은 식물의 위치를 인식할 수 있는 영상처리 시스템과 이를 이용해 정확한 위치의 온도 값을 자동으로 측정할 수 있는 기술에 관한 것이다.Prior art prior to the application of the present invention discloses a device for collecting plant growth environment data by measuring the temperature of plants. This technology is about an image processing system that can recognize the location of plants and a technology that can automatically measure the temperature value of the exact location using this.

또 다른 선행기술로 식물을 촬영한 칼라 영상과 열 영상을 이용하여 비접촉 방식으로 잎의 온도를 측정하는 기술이다. 이 기술에서는 칼라 카메라와 열상 카메라를 통해 촬영한 동일 촬영에 대한 촬영 영상을 이용하여 칼라 영상과 열 영상 간의 매핑 파라미터를 산출하는 제1 단계와, 상기 칼라 카메라와 열상 카메라를 통해 촬영한 대상 식물에 대한 칼라 영상과 열 영상을 각각 획득하는 제2 단계, 상기 제2 단계에서 획득한 칼라 영상에서 식물의 잎 컬러를 근거로 해당 컬러를 갖는 잎 영역을 분리하는 제3 단계, 상기 제1 단계에서 산출된 매핑 파라미터를 근거로 열 영상에서 칼라 영상에서의 잎 영역에 대응되는 영역을 설정하는 제4 단계 및, 상기 제4 단계에서 설정된 열 영상에서의 잎 영역에 대해 화소별 온도 값에 대한 평균치를 산출하여 식물의 잎 온도를 획득하는 제5 단계를 포함하는 기술이 개시되어 있다.Another prior technology is a technology that measures the temperature of leaves in a non-contact manner using color images and thermal images taken of plants. In this technology, the first step is to calculate mapping parameters between color images and thermal images using captured images of the same shot taken through a color camera and a thermal camera, and to target plants photographed through the color camera and thermal camera. A second step of acquiring a color image and a thermal image respectively, a third step of separating a leaf area with the corresponding color based on the color of the plant's leaf in the color image obtained in the second step, calculated in the first step A fourth step of setting an area corresponding to the leaf area in the color image in the thermal image based on the mapping parameters, and calculating the average temperature value for each pixel for the leaf area in the thermal image set in the fourth step. Thus, a technology including a fifth step of obtaining the leaf temperature of the plant is disclosed.

공개특허공보 제10-2013-0035816호Public Patent Publication No. 10-2013-0035816 공개특허공보 제10-2016-0145319호Public Patent Publication No. 10-2016-0145319

스마트팜 등에서 식물의 생육 상태를 측정하기 위하여 많은 센서들이 사용되고 있다. 특히, 식물의 현재 상태를 나타내는 지표로 식물 잎의 온도 즉, 엽온을 사용하고 있다. 스마트팜 내부의 엽온과 스마트팜 내부의 공기 온도의 비교를 통해 작물(식물)의 증산작용이 활발한지 여부를 측정할 수 있다. 또한 엽온을 이용하여 식물의 고온스트레스 발생 진단이 가능하다.Many sensors are used to measure the growth status of plants in smart farms, etc. In particular, the temperature of plant leaves, that is, leaf temperature, is used as an indicator of the current state of the plant. By comparing the leaf temperature inside the smart farm and the air temperature inside the smart farm, it is possible to measure whether the transpiration of crops (plants) is active. Additionally, it is possible to diagnose high-temperature stress in plants using leaf temperature.

이렇게 작물(식물)의 생육상태를 반영하는 엽온을 측정하기 위해서는 열화상카메라나 적외선 온도센서가 추가 적으로 필요하며, 잎 끝에 맞춰 적외선센서를 설치하거나 열화상카메라로 촬영한 영상에서 잎의 이미지만을 구분하는 추가 작업이 필요하다. In order to measure leaf temperature, which reflects the growth state of crops (plants), an additional thermal imaging camera or infrared temperature sensor is required. An infrared sensor is installed at the tip of the leaf or only the image of the leaf is captured from the image taken with a thermal imaging camera. Additional work is needed to differentiate.

상기와 같은 문제를 해결하고자 본 출원 발명에서는 스마트팜 내부온도, 외부온도 및 일사량을 이용하여 회귀모델을 구하여 스마트팜 내부의 엽온을 예측하는 기술을 개발하고자 한다.In order to solve the above problems, the present application seeks to develop a technology to predict leaf temperature inside a smart farm by obtaining a regression model using the smart farm's internal temperature, external temperature, and solar radiation.

상기와 같은 문제를 해결하기 위하여 다음의 과제해결 수단을 제공한다.In order to solve the above problems, the following problem solving means are provided.

스마트팜 내부의 온도(intemp.)를 측정하는 내부온도측정단계; 및An internal temperature measurement step of measuring the temperature (intemp.) inside the smart farm; and

스마트팜 외부의 온도(outtemp.)를 측정하는 외부온도측정단계; 및An external temperature measurement step of measuring the temperature (outtemp.) outside the smart farm; and

스마트팜 내부로 입사되는 태양광의 일사량(rad.)을 측정하는 일사량측정단계; 및A solar radiation measurement step of measuring the solar radiation (rad.) of sunlight entering the smart farm; and

엽온(canopy temp.)은 intemp.×A + outtemp.×B + rad.×C + D(식(1))로 계산되는 것을 특징으로 하는 스마트팜 내부환경 데이터를 이용한 엽온 예측방법을 제공한다.We provide a canopy temperature prediction method using smart farm internal environment data, wherein canopy temperature is calculated as intemp.×A + outtemp.×B + rad.×C + D (Equation (1)).

(상기 A, B, C, D는 각각 곱해지는 측정값에 대한 엽온의 증가에 기여하는 정도에 따른 상수이다)(A, B, C, and D are constants depending on the degree to which they contribute to the increase in leaf temperature relative to the measured value being multiplied)

또한, 스마트팜 내부의 습도(inhum.)를 측정하는 내부습도측정단계; 및In addition, an internal humidity measurement step of measuring the humidity (inhum.) inside the smart farm; and

스마트팜 외부의 습도(outhum.)를 측정하는 외부습도측정단계; 및An external humidity measurement step of measuring the humidity (outhum.) outside the smart farm; and

스마트팜 내부로 입사되는 태양광의 일사량(rad.)을 측정하는 일사량측정단계; 및A solar radiation measurement step of measuring the solar radiation (rad.) of sunlight entering the smart farm; and

작물의 증산량은 일사량 700 w/m2 이상의 조건에서 (inhum. - outhum.) > 0 일때 (inhum.- outhum.)×V로 계산되는 것을 특징으로 하는 스마트팜 내부환경 데이터를 이용한 증산량 예측방법을 제공한다.Provides a transpiration prediction method using smart farm internal environment data, where the transpiration amount of crops is calculated as (inhum.- outhum.) do.

또한, 상기 엽온의 예측에 있어서, (inhum.-outhum.) > 0 인 경우에는 증산이 활발한 것에 해당되어 상기 식(1)을 사용하여 엽온을 예측하며, (inhum.-outhum.) <= 0 인 경우에는 증산이 활발하지 않아 엽온이 증산작용에 의하여 식혀지지 않는 것이므로, 엽온(canopy temp.)은 intemp.×A + outtemp.×B + rad.×C + D + t×E(식(2))로 계산되는 것을 특징으로 하는 스마트팜 내부환경 데이터를 이용한 엽온 예측방법을 제공한다. (상기 A, B, C, D, ㄸ는 각각 곱해지는 측정값에 대한 엽온의 증가에 기여하는 정도에 따른 상수이고, t는 (inhum.-outhum.) <= 0인 동안의 시간을 의미함(1시간 단위임))In addition, in predicting the leaf temperature, if (inhum.-outhum.) > 0, transpiration is active, so the leaf temperature is predicted using Equation (1), and (inhum.-outhum.) <= 0 In this case, transpiration is not active and the leaf temperature is not cooled by transpiration, so the leaf temperature (canopy temp.) is intemp.×A + outtemp.×B + rad.×C + D + t×E (equation (2) It provides a leaf temperature prediction method using smart farm internal environment data, which is calculated as )). (A, B, C, D, ㄸ are constants depending on the degree to which they contribute to the increase in leaf temperature for the measured value being multiplied, and t refers to the time during which (inhum.-outhum.) <= 0 (1 hour unit))

본 출원 발명의 상기와 같은 구성에 의하여 스마트팜 내부에서 재배하는 식물의 엽온을 직접 측정하지 않고도 식물의 엽온을 스마트팜 내부와 외부에 구비된 환경 변수 값으로부터 예측할 수 있는 수단을 제공하는 효과가 있다.The above-mentioned configuration of the invention of this application has the effect of providing a means of predicting the leaf temperature of plants from the environmental variable values provided inside and outside the smart farm without directly measuring the leaf temperature of the plants grown inside the smart farm. .

도 1은 본 발명의 엽온 예측에 사용한 데이터 셋
도 2는 본 발명의 예측 변수 설정을 위한 시뮬레이션 결과
도 3은 본 발명의 선형 예측 방법을 비교하기 위한 시뮬레이션 결과
도 4는 본 발명의 엽온 실측 및 예측 비교 그래프(여름)
도 5는 본 발명의 엽온 실측 및 예측 비교 그래프(겨울)
Figure 1 is a data set used to predict leaf temperature of the present invention.
Figure 2 is a simulation result for setting the predictor variable of the present invention
Figure 3 is a simulation result for comparing the linear prediction method of the present invention
Figure 4 is a graph comparing actual leaf temperature measurements and predictions of the present invention (summer)
Figure 5 is a graph comparing actual leaf temperature measurements and predictions of the present invention (winter)

작물의 수분스트레스는 작물이 토양으로부터 흡수하는 양보다 더 많은 수분을 증산작용으로 잃어버림으로써 체내 수분함량이 줄고, 이로 인해 작물 생장이 감소하는 현상을 의미한다. 작물이 수분 스트레스를 받게 되면 기공이 감소하여 증산작용이 저하되어 엽온(캐노피 온도, Canopy temperature, Tc)가 상승한다.Water stress in crops refers to a phenomenon in which the crop loses more moisture through transpiration than it absorbs from the soil, resulting in a decrease in moisture content in the body, which results in a decrease in crop growth. When crops are subjected to moisture stress, stomata are reduced, transpiration is reduced, and leaf temperature (canopy temperature, Tc) increases.

따라서, 작물의 엽온을 정확하게 측정하여 작물의 현재 생육상태를 확인하는 것이 중요하다. 그러나, 작물의 엽온을 측정하는 것이 물리적으로 간단하지 않은 문제가 있다. 하루에 1 내지 2회 정도 측정한다면, 인력으로 간단히 측정이 가능하다. 그러나, 하루 종일 엽온을 모니터링해야 한다면 쉽지 않다. 스마트팜 내부에는 환기 팬도 가동하고, 작업자도 이동하고 있기 때문에 엽온 센서가 잎에 정확하게 접촉하거나 원격으로 열화상카메라 등으로 정확한 위치에서 연속하여 측정하는 것이 쉽지 않다.Therefore, it is important to accurately measure the leaf temperature of the crop to check the current growth status of the crop. However, there is a problem in that measuring the leaf temperature of crops is not physically simple. If you measure once or twice a day, you can easily measure by manpower. However, it is not easy to monitor leaf temperature throughout the day. Because ventilation fans are running and workers are moving inside the smart farm, it is not easy for the leaf temperature sensor to accurately contact the leaves or to measure continuously at the exact location remotely using a thermal imaging camera.

본 출원 발명은 상기와 같은 직접 측정방법이 아닌, 실내외의 환경정보를 이용하여 간접적으로 엽온을 예측하는 방법을 제시하고자 한다.The present application seeks to propose a method of predicting leaf temperature indirectly using indoor and outdoor environmental information, rather than the direct measurement method as described above.

엽온과 관련된 생육 환경데이터는 도1에 도시된 바와 같다. Growth environment data related to leaf temperature are as shown in Figure 1.

본 발명에서는 스마트팜 내부 환경변수로 스마트팜 내부온도(intemp.), 스마트팜 내부습도(inhum.), 이산화탄소농도(co2), 태양복사량(rad)과, 스마트팜 외부 환경변수로 스마트팜 외부온도(outtemp.), 스마트팜 외부습도(outhum.), 스마트팜 외부풍속(windspeed)을 사용하였다.In the present invention, the smart farm internal environmental variables are smart farm internal temperature (intemp.), smart farm internal humidity (inhum.), carbon dioxide concentration (co2), and solar radiation (rad), and the smart farm external environmental variables are smart farm external temperature. (outtemp.), smart farm external humidity (outhum.), and smart farm external wind speed (windspeed) were used.

우선 이들 환경변수와 엽온과의 상관관계를 먼저 분석하였다.First, the correlation between these environmental variables and leaf temperature was analyzed.

3개의 변수 그룹을 만들고 각각에 대한 선형회귀 적합도를 분석하여 최적의 엽온 예측 변수를 선정하였다.Three variable groups were created and the linear regression fit for each was analyzed to select the optimal leaf temperature predictor.

제1 그룹에 사용한 input_data: 내부 온/습도, CO2, 일사량, 외부 온/습도, 풍속input_data used for the first group: internal temperature/humidity, CO2, solar radiation, external temperature/humidity, wind speed

제2 그룹에 사용한 input_data: 내/외부 온도, CO2, 일사량, 풍속input_data used for group 2: inside/outside temperature, CO2, solar radiation, wind speed

제3 그룹에 사용한 input_data: 내/외부 온도, 일사량input_data used for the third group: inside/outside temperature, solar radiation

이다. am.

분석결과 제3 그룹의 스마트팜의 내부 및 외부 온도와 일사량을 사용한 결과가 결정계수 0.702, RMSE 1.753으로 가장 값이 좋았다. 가장 좋게 나타났다. 실험결과는 도2에 도시하였다.As a result of the analysis, the result using the internal and external temperature and solar radiation of the third group's smart farm had the best value with a coefficient of determination of 0.702 and RMSE of 1.753. It turned out the best. The experimental results are shown in Figure 2.

상기 제3 그룹의 선정된 변수를 바탕으로 더 나은 엽온 예측 방법을 찾아보기 위하여 기존에 사용하였던 선형회귀 방법에 선형 모델(Linear model)의 예측도(accuracy) 및/ 또는 해석가능성(interpretability)을 더 높이기 위해 여러 정규화(regularization) 기능이 있는 Ridge와 Lasso 방법으로 엽온 예측을 추가로 수행하고 그 결과를 분석하였다.In order to find a better leaf temperature prediction method based on the selected variables of the third group, the accuracy and/or interpretability of the linear model was added to the previously used linear regression method. To improve this, additional leaf temperature prediction was performed using the Ridge and Lasso methods with various regularization functions, and the results were analyzed.

방법 1 : 선형회귀Method 1: Linear Regression

방법 2 : RidgeMethod 2: Ridge

방법 3 : LassoMethod 3: Lasso

그 결과 방법 3인 Lasso 모델을 사용한 경우에 결정계수 0.728, RMSE 1.676으로 가장 값이 좋았다. 도 3에 분석결과를 도시하고 있다.As a result, method 3, the Lasso model, had the best value with a coefficient of determination of 0.728 and RMSE of 1.676. Figure 3 shows the analysis results.

이때 예측 모델의 회귀식은 다음과 같다.At this time, the regression equation of the prediction model is as follows.

엽온= 9.651772156859273 + 0.46214678×내부온도(intemp.) + 0.08514499×외부온도(outtemp.)+0.00351328×일사량(rad.) -------- (식0)Leaf temperature = 9.651772156859273 + 0.46214678

으로 예측할 수 있다.can be predicted.

그러나, 이러한 예측은 실험에 사용한 데이터가 적용된 실증팜에 한정되는 것으로 해당 수식(식0)을 각각 사용하기 위해서는 스마트팜 내부온도, 외부온도, 일사량 등의 변화에 대한 가중치의 수정이 필요합니다. 또한 이러한 가중치의 수정은 해당 스마트팜의 구조와 위치 등과도 관련이 있습니다.However, these predictions are limited to the demonstration farm to which the data used in the experiment was applied, and in order to use each of the formulas (Equation 0), it is necessary to modify the weights for changes in the smart farm's internal temperature, external temperature, and solar radiation. Additionally, modification of these weights is related to the structure and location of the smart farm.

상기 환경 조건 변수와 관련하여 상기 가중치를 설정하는 방법은 하기와 같습니다.The method of setting the weight in relation to the environmental condition variable is as follows.

<실시예 1> <Example 1>

스마트팜의 내부에 작물의 증산에 필요한 정도의 바람이 항시 존재하고, 외부온도와 내부 온도가 5도 이상 차이가 있을 정도의 스마트팜 구조이며, 상대 습도가 70%를 넘지 않아 식물의 증산 작용이 무난한 경우에는 작물이 증산작용을 하여 태양 복사가 강한 경우에도 엽온을 정상적으로 발산하여 상기 식(0)으로 작물의 엽온의 예측이 가능하다.Inside the smart farm, there is always enough wind necessary for the transpiration of crops, and the smart farm structure is such that there is a difference of more than 5 degrees between the external and internal temperatures, and the relative humidity does not exceed 70%, so the transpiration function of the plants is maintained. In normal cases, the crop transpires and radiates leaf temperature normally even when solar radiation is strong, making it possible to predict the crop's leaf temperature using the equation (0) above.

한편 추가적으로 스마트팜의 시설의 차이에 의한 엽온 예측 차이를 보정하기 위하여 내부온도, 외부온도, 태양의 복사량과 엽온을 1회 이상 측정하여 식(0)의 계수들을 보정 또는 교정하여 사용할 수 있다.Meanwhile, in order to additionally correct for differences in leaf temperature prediction due to differences in smart farm facilities, the internal temperature, external temperature, solar radiation, and leaf temperature can be measured more than once and the coefficients in equation (0) can be corrected or calibrated.

<실시예 2><Example 2>

한 여름과 같이 스마트팜의 내부의 온도가 외부의 온도 보다 높고, 작물의 증산에 필요한 정도의 바람이 있으나, 상대 습도가 90%를 넘기 때문에 식물의 증산 작용이 쉽지 않은 경우 무난한 경우에는 작물의 증산작용이 원활하지 않기 때문에 증산 작용에 의한 작물의 엽온 발산이 불가능하여 작물의 엽온이 상기 식(0)으로 예측한 온도와 비교하여 더 높은 온도일 수 있다. In cases where the temperature inside the smart farm is higher than the outside temperature, such as in the middle of summer, and there is enough wind for transpiration of crops, but the relative humidity exceeds 90%, the transpiration of plants is not easy, in normal cases, transpiration of crops is possible. Because the operation is not smooth, it is impossible to dissipate the leaf temperature of the crop through transpiration, so the leaf temperature of the crop may be at a higher temperature compared to the temperature predicted by the above equation (0).

이러한 조건의 경우For these conditions

스마트팜 내부의 습도(inhum.)와 스마트팜 외부의 습도(outhum.) 와 한낮 조건을 확인하기 위하여 스마트팜 내부로 입사되는 태양광의 일사량(rad.)을 측정한 후 증산율을 계산하여 엽온의 예측 온도를 수정할 수 있다.In order to check the humidity (inhum.) inside the smart farm, the humidity (outhum.) outside the smart farm, and the midday conditions, the amount of solar radiation (rad.) entering the smart farm is measured, and then the transpiration rate is calculated to determine the leaf temperature. The predicted temperature can be modified.

즉, 작물의 증산량은 일사량 > 700 w/m2 이상의 조건에서 (inhum. - outhum.) > 0 일때 (inhum.- outhum.)*V로 계산할 수 있으며, 상기 V는 가중치 값이다. 스마트팜 내부의 습도가 높은 것은 증산에 의하여 발생한 수증기가 존재하는 것으로 작물이 증산작용을 하고 있다는 것으로 해석되어 지며, 따라서 식(0)의 엽온 예측 모델을 그대로 사용할 수 있다.In other words, the transpiration amount of the crop can be calculated as (inhum.-outhum.)*V when (inhum.-outhum.) > 0 under the condition of solar radiation > 700 w/m2, where V is a weight value. The high humidity inside the smart farm is interpreted as the presence of water vapor generated by transpiration, which means that the crops are transpiration. Therefore, the leaf temperature prediction model of equation (0) can be used as is.

그러나, 일사량 > 700 w/m2 이상의 조건에서 (inhum. - outhum.) <= 0 일때는 외부 습도가 높아 작물이 증산을 할 수 없는 조건으로 이때는 작물이 증산작용을 할 수 없어 엽온의 증가가 예측된다. 따라서, 상기 식(0)의 모델에 -(inhum.- outhum.)×V 항을 더 부가하여 증산이 100%가 되지 않는 경우의 엽온 예측모델을 제시하였다. 그러나, 이러한 엽온의 증가는 시간에 따라 누적되는 값이므로, 태양복사에 의한 일사량이 700 w/m2 이상인 시간을 누적하여 엽온을 예측하여야 한다. 그래서 일사량이 700 w/m2 이상인 시간을 곱하여 엽온을 예측하는 식은 하기와 같다.However, under conditions of solar irradiation > 700 w/m2 (inhum. - outhum.) <= 0, the external humidity is high and the crop cannot transpire. In this case, the crop cannot transpire, so an increase in leaf temperature is expected. do. Therefore, by adding -(inhum.-outhum.) However, since this increase in leaf temperature is a cumulative value over time, the leaf temperature must be predicted by accumulating the time when the amount of solar radiation is more than 700 w/m2. Therefore, the formula for predicting leaf temperature by multiplying the time when solar radiation is above 700 w/m2 is as follows.

여기서 시간(t)의 단위는 hour이며 분 단위는 소수점으로 표시한다.Here, the unit of time (t) is hour, and the unit of minutes is expressed as a decimal point.

엽온= 9.651772156859273 + 0.46214678×내부온도 + 0.08514499×외부온도+0.00351328×일사량 -(inhum.- outhum.)*V*t-------- (식0.1)Leaf temperature = 9.651772156859273 + 0.46214678

상기의 식으로부터 스마트팜에 따라 차이가 있을 수 있는 변수를 고려하면, 엽온(canopy temp.) = intemp.×A + outtemp.×B + rad.×C + D(식(1))로 표시되는 일반식을 구할 수 있으며, 스마트팜 내부의 습도가 높은 경우에는 습도를 고려한 (canopy temp.) = intemp.×A + outtemp.×B + rad.×C + D + t×E(식(2))를 일반화할 수 있었다.Considering the variables that may differ depending on the smart farm from the above equation, canopy temperature (canopy temp.) = intemp.×A + outtemp.×B + rad.×C + D (equation (1)) A general equation can be obtained, and if the humidity inside the smart farm is high, (canopy temp.) = intemp.×A + outtemp.×B + rad.×C + D + t×E (equation (2)) considering humidity ) could be generalized.

도4는 실제 엽온을 측정한 값과 본 발명에 의하여 예측한 엽온을 비교한 결과이다. 5월 중순부터 6월 초순의 실험결과이며, 실험에 사용한 작물은 토마토이다. 예측 결과보다 실제 측정 결과의 온도가 조금 높게 나온 것을 확인할 수 있다.Figure 4 shows the results of comparing the actual leaf temperature measured value and the leaf temperature predicted by the present invention. These are the results of an experiment from mid-May to early June, and the crop used in the experiment was tomatoes. It can be seen that the temperature of the actual measurement result is slightly higher than the predicted result.

도5는 12월 중순부터 2월말까지의 엽온의 실제측정 결과와 예측결과를 비교하여 도시한 그래프이다. 실제 엽온이 조금 아래쪽에 표시된 것을 볼 수 있다. Figure 5 is a graph comparing the actual measurement results and predicted results of leaf temperature from mid-December to the end of February. You can see that the actual leaf temperature is displayed slightly below.

그러나, 도4 및 도5의 그래프의 측정결과는 실제 엽온을 측정하는 노력과 시간과 비교하여 실제 사용할 수 있을 만큼의 예측력을 충족하는 실험결과이다.However, the measurement results in the graphs of Figures 4 and 5 are experimental results that meet the predictive power sufficient for actual use compared to the effort and time to measure actual leaf temperature.

그러나, 이러한 차이가 있는 원인을 우리나라 기후 사정을 고려하여 살펴보면, 겨울에는 습도가 낮고, 여름으로 갈수록 습도가 높은 기후 환경을 가지고 있다.However, if we look at the cause of this difference by considering our country's climate, we have a climate environment where humidity is low in winter and humidity is high in summer.

이러한 기후 사정을 고려하여 본다면, 습도가 낮은 겨울에는 상대습도가 낮기 때문에 증산이 원활하기 때문에 엽온이 쉽게 낮아질 수 있으며, 상대습도가 올라가는 여름에는 증산이 원활하지 않기 때문에 예측온도가 실측온도 보다 낮은 문제가 있었습니다. 이러한 문제를 해결하고자 스마트팜 내부습도와 외부습도를 추가로 수식에 추가하여 이를 보정하고자 하였습니다. 수식(2)가 이러한 내용입니다. 즉, 일사량이 많은 시간대에 스마트팜 내부의 습도가 높으면 증상이 어려워 일사량에 의하여 엽온이 올라가는 것을 일사량이 많은 시간대에 누적시간 개념으로 엽온을 보정하였습니다. 상기와 같은 보정에 의하여 좀 더 정확한 엽온 예측 결과를 구할 수 있었습니다.Taking these climate conditions into consideration, in the winter when relative humidity is low, transpiration is smooth, so the leaf temperature can easily drop, and in the summer when the relative humidity is high, because transpiration is not smooth, the predicted temperature is lower than the actual temperature. There was. To solve this problem, we attempted to correct this by adding the internal humidity and external humidity of the smart farm to the formula. This is what equation (2) contains. In other words, if the humidity inside the smart farm is high during times of high solar radiation, symptoms are difficult, and the leaf temperature rises due to solar radiation, so the leaf temperature was corrected using the concept of cumulative time during times of high solar radiation. By making the above corrections, we were able to obtain more accurate leaf temperature prediction results.

상기와 같은 발명의 작용효과를 나타내기 위한 발명의 구성은 다음과 같다.The structure of the invention to demonstrate the effects of the invention as described above is as follows.

스마트팜 내부의 온도(intemp.)를 측정하는 내부온도측정단계; 및An internal temperature measurement step of measuring the temperature (intemp.) inside the smart farm; and

스마트팜 외부의 온도(outtemp.)를 측정하는 외부온도측정단계; 및An external temperature measurement step of measuring the temperature (outtemp.) outside the smart farm; and

스마트팜 내부로 입사되는 태양광의 일사량(rad.)을 측정하는 일사량측정단계; 및A solar radiation measurement step of measuring the solar radiation (rad.) of sunlight entering the smart farm; and

엽온(canopy temp.)은 intemp.×A + outtemp.×B + rad.×C + D(식(1))로 계산되는 것을 특징으로 하는 스마트팜 내부환경 데이터를 이용한 엽온 예측방법을 제공한다.We provide a canopy temperature prediction method using smart farm internal environment data, wherein canopy temperature is calculated as intemp.×A + outtemp.×B + rad.×C + D (Equation (1)).

(상기 A, B, C, D는 각각 곱해지는 측정값에 대한 엽온의 증가에 기여하는 정도에 따른 상수이다)(A, B, C, and D are constants depending on the degree to which they contribute to the increase in leaf temperature relative to the measured value being multiplied)

또한, 스마트팜 내부의 습도(inhum.)를 측정하는 내부습도측정단계; 및In addition, an internal humidity measurement step of measuring the humidity (inhum.) inside the smart farm; and

스마트팜 외부의 습도(outhum.)를 측정하는 외부습도측정단계; 및An external humidity measurement step of measuring the humidity (outhum.) outside the smart farm; and

스마트팜 내부로 입사되는 태양광의 일사량(rad.)을 측정하는 일사량측정단계; 및A solar radiation measurement step of measuring the solar radiation (rad.) of sunlight entering the smart farm; and

작물의 증산량은 일사량 700 w/m2 이상의 조건에서 (inhum. - outhum.) > 0 일때 (inhum.- outhum.)×V로 계산되는 것을 특징으로 하는 스마트팜 내부환경 데이터를 이용한 증산량 예측방법을 제공한다.Provides a transpiration prediction method using smart farm internal environment data, where the transpiration amount of crops is calculated as (inhum.- outhum.) do.

또한, 상기 엽온의 예측에 있어서, (inhum.-outhum.) > 0 인 경우에는 증산이 활발한 것에 해당되어 상기 식(1)을 사용하여 엽온을 예측하며, (inhum.-outhum.) <= 0 인 경우에는 증산이 활발하지 않아 엽온이 증산작용에 의하여 식혀지지 않는 것이므로, 엽온(canopy temp.)은 intemp.×A + outtemp.×B + rad.×C + D + t×E(식(2))로 계산되는 것을 특징으로 하는 스마트팜 내부환경 데이터를 이용한 엽온 예측방법을 제공한다. (상기 A, B, C, D, ㄸ는 각각 곱해지는 측정값에 대한 엽온의 증가에 기여하는 정도에 따른 상수이고, t는 (inhum.-outhum.) <= 0인 동안의 시간을 의미함(1시간 단위임))In addition, in predicting the leaf temperature, if (inhum.-outhum.) > 0, transpiration is active, so the leaf temperature is predicted using Equation (1), and (inhum.-outhum.) <= 0 In this case, transpiration is not active and the leaf temperature is not cooled by transpiration, so the leaf temperature (canopy temp.) is intemp.×A + outtemp.×B + rad.×C + D + t×E (equation (2) It provides a leaf temperature prediction method using smart farm internal environment data, which is calculated as )). (A, B, C, D, ㄸ are constants depending on the degree to which they contribute to the increase in leaf temperature for the measured value being multiplied, and t refers to the time during which (inhum.-outhum.) <= 0 (1 hour unit))

Claims (3)

스마트팜 내부의 온도(intemp.)를 측정하는 내부온도측정단계; 및
스마트팜 외부의 온도(outtemp.)를 측정하는 외부온도측정단계; 및
스마트팜 내부로 입사되는 태양광의 일사량(rad.)을 측정하는 일사량측정단계; 및
엽온(canopy temp.)은 intemp.×A + outtemp.×B + rad.×C + D(식(1))로 계산되는 것을 특징으로 하는 스마트팜 내부환경 데이터를 이용한 엽온 예측방법.
(상기 A, B, C, D는 각각 곱해지는 측정값에 대한 엽온의 증가에 기여하는 정도에 따른 상수값)
An internal temperature measurement step of measuring the temperature (intemp.) inside the smart farm; and
An external temperature measurement step of measuring the temperature (outtemp.) outside the smart farm; and
A solar radiation measurement step of measuring the solar radiation (rad.) of sunlight entering the smart farm; and
Canopy temperature prediction method using smart farm internal environment data, characterized in that canopy temperature is calculated as intemp.×A + outtemp.×B + rad.×C + D (Equation (1)).
(A, B, C, and D are constant values depending on the degree to which they contribute to the increase in leaf temperature relative to the measured value being multiplied)
제1항에 있어서,
스마트팜 내부의 습도(inhum.)를 측정하는 내부습도측정단계; 및
스마트팜 외부의 습도(outhum.)를 측정하는 외부습도측정단계; 및
스마트팜 내부로 입사되는 태양광의 일사량(rad.)을 측정하는 일사량측정단계; 및작물의 증산량은 일사량 700 w/m2 이상의 조건에서 (inhum. - outhum.) > 0 일때 (inhum.- outhum.)×V로 계산되는 것을 특징으로 하는 특징으로 하는 스마트팜 내부환경 데이터를 이용한 엽온 예측방법.
According to paragraph 1,
An internal humidity measurement step of measuring the humidity (inhum.) inside the smart farm; and
An external humidity measurement step of measuring the humidity (outhum.) outside the smart farm; and
A solar radiation measurement step of measuring the solar radiation (rad.) of sunlight entering the smart farm; And the transpiration amount of the crop is calculated as (inhum.-outhum.) × V when (inhum.-outhum.) > 0 under conditions of solar radiation of 700 w/m2 or more. Prediction method.
제2항에 있어서,
상기 엽온의 예측에 있어서, (inhum.-outhum.) > 0 인 경우에는 증산이 활발한 것에 해당되어 상기 식(1)을 사용하여 엽온을 예측하며, (inhum.-outhum.) <= 0 인 경우에는 증산이 활발하지 않아 엽온이 증산작용에 의하여 식혀지지 않는 것이므로, 엽온(canopy temp.)은 intemp.×A + outtemp.×B + rad.×C + D + t×E(식(2))로 계산되는 것을 특징으로 하는 스마트팜 내부환경 데이터를 이용한 엽온 예측방법. (상기 A, B, C, D, ㄸ는 각각 곱해지는 측정값에 대한 엽온의 증가에 기여하는 정도에 따른 상수이고, t는 (inhum.-outhum.) <= 0인 동안의 시간을 의미함(1시간 단위임))
According to paragraph 2,
In predicting the leaf temperature, if (inhum.-outhum.) > 0, transpiration is active and the leaf temperature is predicted using Equation (1), and if (inhum.-outhum.) <= 0 Since transpiration is not active and the leaf temperature is not cooled by transpiration, the canopy temperature is intemp.×A + outtemp.×B + rad.×C + D + t×E (Equation (2)) A leaf temperature prediction method using smart farm internal environment data, which is calculated as . (A, B, C, D, ㄸ are constants depending on the degree to which they contribute to the increase in leaf temperature for the measured value being multiplied, and t refers to the time during which (inhum.-outhum.) <= 0 (1 hour unit))
KR1020220084210A 2022-07-08 2022-07-08 A continuous leaf temperature prediction method using smart farm environmental data. KR20240007390A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020220084210A KR20240007390A (en) 2022-07-08 2022-07-08 A continuous leaf temperature prediction method using smart farm environmental data.

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020220084210A KR20240007390A (en) 2022-07-08 2022-07-08 A continuous leaf temperature prediction method using smart farm environmental data.

Publications (1)

Publication Number Publication Date
KR20240007390A true KR20240007390A (en) 2024-01-16

Family

ID=89719809

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020220084210A KR20240007390A (en) 2022-07-08 2022-07-08 A continuous leaf temperature prediction method using smart farm environmental data.

Country Status (1)

Country Link
KR (1) KR20240007390A (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20130035816A (en) 2011-09-30 2013-04-09 순천대학교 산학협력단 Measuring temperature of the plant by using the image processing method
KR20160145319A (en) 2015-06-10 2016-12-20 주식회사 엘앤와이비젼 Method for temperature measuring of plant leaf using muti-image

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20130035816A (en) 2011-09-30 2013-04-09 순천대학교 산학협력단 Measuring temperature of the plant by using the image processing method
KR20160145319A (en) 2015-06-10 2016-12-20 주식회사 엘앤와이비젼 Method for temperature measuring of plant leaf using muti-image

Similar Documents

Publication Publication Date Title
Wang et al. Evaluation of ASTER and MODIS land surface temperature and emissivity products using long-term surface longwave radiation observations at SURFRAD sites
Cohen et al. Crop water status estimation using thermography: multi-year model development using ground-based thermal images
Papadavid et al. A modified SEBAL modeling approach for estimating crop evapotranspiration in semi-arid conditions
Marek et al. Post-processing techniques for reducing errors in weighing lysimeter evapotranspiration (ET) datasets
Mendiguren et al. Spatial pattern evaluation of a calibrated national hydrological model–a remote-sensing-based diagnostic approach
Cammalleri et al. Daily evapotranspiration assessment by means of residual surface energy balance modeling: A critical analysis under a wide range of water availability
Yang et al. Assessing light to moderate grazing effects on grassland production using satellite imagery
WO2020047579A1 (en) Method and system for plant stress determination and irrigation based thereon
Chen et al. A comparison of two adaptive multivariate analysis methods (PLSR and ANN) for winter wheat yield forecasting using Landsat-8 OLI images
KR20210045253A (en) Irrigation system and method using water shortage diagnosis of crops
Taki et al. Application of Neural Networks and multiple regression models in greenhouse climate estimation
Osroosh et al. Detecting fruit surface wetness using a custom-built low-resolution thermal-RGB imager
Soegaard Fluxes of carbon dioxide, water vapour and sensible heat in a boreal agricultural area of Sweden—scaled from canopy to landscape level
Dhillon et al. Development of a continuous leaf monitoring system to predict plant water status
KR20240007390A (en) A continuous leaf temperature prediction method using smart farm environmental data.
Liu et al. Fractional coverage rather than green chromatic coordinate is a robust indicator to track grassland phenology using smartphone photography
Sharabian et al. Evaluation of an active remote sensor for monitoring winter wheat growth status
Gutiérrez et al. On-the-go thermal imaging for water status assessment in commercial vineyards
Basit et al. Application of simplified surface energy balance index (S-SEBI) for crop evapotranspiration using Landsat 8
Nagy Thermographic evaluation of water stress in an apple orchard
CN116187859A (en) Method, device, equipment and medium for monitoring and diagnosing low-temperature cold injury of crops in irrigation area
Nguyen-Xuan et al. A dynamic model for temperature prediction in glass greenhouse
Ranjan et al. Internet of Things enabled crop physiology sensing system for abiotic crop stress management in apple and sweet cherry
Rougerie-Durocher et al. Measuring and modelling of apple flower stigma temperature as a step towards improved fire blight prediction
Kumar et al. Evapotranspiration mapping for agricultural water management: An overview

Legal Events

Date Code Title Description
J201 Request for trial against refusal decision