CN115753712B - SIF anomaly detection method, device, system, terminal and medium based on PAR data - Google Patents
SIF anomaly detection method, device, system, terminal and medium based on PAR data Download PDFInfo
- Publication number
- CN115753712B CN115753712B CN202211441181.2A CN202211441181A CN115753712B CN 115753712 B CN115753712 B CN 115753712B CN 202211441181 A CN202211441181 A CN 202211441181A CN 115753712 B CN115753712 B CN 115753712B
- Authority
- CN
- China
- Prior art keywords
- par
- sif
- data
- spectrum data
- detection time
- Prior art date
- Legal status (The legal status 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 status listed.)
- Active
Links
Images
Landscapes
- Maintenance And Management Of Digital Transmission (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The invention discloses a SIF anomaly detection method, device, system, terminal and medium based on PAR data, and relates to the field of remote sensing data processing. Comprising the following steps: acquiring spectrum data; constructing a detection time group comprising at least one detection time unit based on the spectral data, the detection time unit comprising: an uplink spectrum data acquisition time t, and two downlink spectrum data acquisition times t-1 and t+1 adjacent to the acquisition time t; based on the pre-acquired PAR data, acquiring and calculating the absolute difference value PAR between the PAR at the t moment and the PAR at two adjacent moments in any detection time unit 1 And PAR 2 When PAR is 1 And PAR 2 Are all greater than or equal to a preset PAR Threshold value And judging the SIF calculated based on the corresponding spectrum data of the acquisition time of the detection time unit as an abnormal value and marking. According to the invention, through comparing the PAR data results monitored by the PAR monitor, whether the SIF is abnormal or not is checked, and SIF abnormal values are removed, so that the overall SIF result error is reduced.
Description
Technical Field
The invention relates to the field of remote sensing data processing, in particular to a SIF anomaly detection method, device, system, terminal and medium based on PAR data.
Background
Sunlight-induced chlorophyll fluorescence (Solar-Induced Chlorophyll Fluorescence/Sun-Induced Chlorophyll Fluorescence, SIF) is a fluorescent light generated by plants under natural light irradiation, is highly correlated with photosynthesis, and is a probe for detecting photosynthesis. Currently, the detection of SIF by researchers has been widely used for the tracking of photosynthesis and vegetation stress monitoring.
In the related study of SIF, hyperspectral measurement of the quality of up-and-down radiation data is the key for inverting vegetation SIF, and vegetation SIF brightness only accounts for about 5% of vegetation reflected radiation (i.e. up-radiation) radiance, so that some fine deviation of spectrum data will greatly influence SIF calculation accuracy. In the existing SIF observation and calculation method, the SIF observation and calculation method based on the double spectrometers is not applicable in most cases because of large overall uncertainty and high cost. The single spectrometer alternating SIF observation calculation method is a practical observation calculation method. However, under the condition of rapid change of light conditions such as overcast and rainy, cloudy and the like, huge deviation can occur in the spectrum data of a certain time acquired by adopting a single spectrometer, and the accuracy of the SIF calculated value acquired on the basis of the spectrum data containing the huge deviation is low.
In summary, there is no effective method for detecting SIF outliers in hyperspectral data alternately measured by a single spectrometer, and therefore, a method for eliminating the existing biased spectral data is needed to improve the accuracy of SIF calculation.
Disclosure of Invention
The invention aims to provide a SIF anomaly detection method, device, system, terminal and medium based on PAR data, so as to solve the problem of low SIF calculation accuracy caused by huge deviation spectrum data.
The first aspect of the present invention provides a SIF anomaly detection method based on PAR data, applied to a central control unit, comprising:
acquiring spectrum data;
constructing a detection time group comprising at least one detection time unit based on the spectral data, the detection time unit comprising: an uplink spectrum data acquisition time t, and two downlink spectrum data acquisition times t-1 and t+1 adjacent to the acquisition time t;
based on the pre-acquired PAR data and the acquisition time recorded by any one detection time unit, acquiring and calculating the PAR data of each acquisition time in any one detection time unit, and calculating the absolute difference value PAR 1 And PAR 2 Wherein PAR is t PAR data representing acquisition time t; PAR (PAR) t-1 PAR data representing acquisition time t-1; PAR (PAR) t+1 PAR data representing acquisition time t+1;
PAR 1 ＝|PAR t -PAR t-1 | (1)
PAR 2 ＝|PAR t -PAR t+1 | (2)
when PAR is 1 Greater than or equal to a predetermined PAR Threshold value And PAR 2 Greater than or equal to a predetermined PAR Threshold value Then determining based on the arbitrary detection time unitAnd (5) taking SIF obtained by calculating spectrum data corresponding to the acquisition time as an abnormal value and marking.
In the above embodiment of the present invention, it is preferable that correction calculation for calculating the SIF anomaly value is further included, specifically: when PAR is 1 And PAR 2 Any one of which is smaller than a preset PAR Threshold value If yes, adopting any detection time unit to be smaller than the preset PAR Threshold value And (3) calculating SIF according to the downlink spectrum data and the uplink spectrum data at the time corresponding to the PAR of the SIF, and finishing SIF abnormal value correction and marking.
In the above embodiment of the present invention, it is preferable that when PAR 1 And PAR 2 Are all smaller than the preset PAR Threshold value And judging that the SIF calculated by adopting the spectrum data corresponding to any one detection time unit is a normal value.
In the above embodiment of the present invention, preferably, the PAR data is collected by a PAR monitor, and the PAR monitor is disposed at a collecting end of a downlink radiation light path.
In the above embodiment of the present invention, it is preferable that the method further includes removing SIF outliers, specifically: the central control unit deletes the SIF outlier in the finally output SIF data based on the SIF outlier flag.
The second aspect of the present invention provides an apparatus for SIF anomaly detection based on PAR data, the apparatus comprising a central control unit, a spectrometer connected to the central control unit, and an optical path collection unit; the optical path collecting unit comprises a downlink radiation optical path unit and an uplink radiation optical path unit; the apparatus further comprises: and the photosynthetic active radiation monitor is connected with the central control unit and is arranged at the collecting end of the downlink radiation light path.
A third aspect of the present invention provides a system for SIF anomaly detection based on PAR data, the system comprising:
the data acquisition module acquires spectrum data and PAR data;
a detection time framework module that constructs a detection time group comprising at least one detection time cell based on the spectral data, the detection time cell comprising: an uplink spectrum data acquisition time t, and two downlink spectrum data acquisition times t-1 and t+1 adjacent to the acquisition time t;
the searching and calculating module is used for acquiring the PAR at the t moment in the detection time unit i to be detected based on the PAR data t PAR at time t-1 adjacent to said time t t-1 And PAR at time t+1 adjacent to said time t t+1 The method comprises the steps of carrying out a first treatment on the surface of the Calculating the absolute difference PAR according to the formula (3) and the formula (4) 1 And PAR 2 ；
PAR 1 ＝|PAR t -PAR t-1 | (3)
PAR 2 ＝|PAR t -PAR t+1 | (4)
Judgment modules for respectively judging PAR 1 And PAR 2 And preset PAR Threshold value When PAR is related to 1 And PAR 2 Are all greater than or equal to a preset PAR Threshold value If so, judging that SIF calculated based on the spectrum data corresponding to the detection time unit i is an abnormal value; when PAR is 1 And PAR 2 Any one of which is smaller than a preset PAR Threshold value When the SIF is detected, determining that SIF calculated based on the spectrum data corresponding to the detection time unit i needs to be corrected;
and a marking module, configured to mark the SIF outlier.
In the above embodiment of the present invention, preferably, the system further includes: the correction module is used for correcting the SIF required to be corrected based on the result calculated based on the spectrum data corresponding to the detection time unit i, specifically: acquiring a detection time unit corresponding to SIF to be corrected; based on less than a preset PAR Threshold value The PAR corresponding time of the time is obtained to obtain the downlink spectrum data A of the time; acquiring uplink spectrum data B based on the uplink spectrum data acquisition time in the detection time unit corresponding to the SIF to be corrected; and calculating SIF according to the downlink spectrum data A and the uplink spectrum data B, and finishing SIF correction and marking.
A third aspect of the present invention provides a terminal, comprising: one or more processors; and a memory for storing one or more programs that, when executed by the one or more processors, cause the one or more processors to implement the PAR data-based SIF anomaly detection method.
A fourth aspect of the present invention provides a computer-readable medium having stored thereon a computer program which, when executed by a processor, implements a SIF abnormality detection method based on PAR data.
The beneficial effects of the invention are as follows: according to the technical scheme, the SIF abnormality detection method based on the PAR data is used for checking whether SIF is abnormal or not by comparing the PAR data results monitored by the PAR monitor, eliminating SIF abnormal values and reducing overall SIF result errors. The method comprises the following steps: the technical scheme of the invention is that the SIF anomaly detection method based on the PAR data is applied to single spectrometer for alternately detecting SIF, whether the irradiance at a certain measurement time t of uplink spectrum data is greatly deviated from the irradiance at the measurement time of the downlink spectrum data adjacent to the front and rear is judged according to the PAR data measured in real time, and if so, the weather condition that the currently judged time t is in rapid light change is indicated. And eliminating SIF calculated values obtained by inversion of the spectrum measured value at the current judging moment t from the SIF calculated value final result. Compared with the prior SIF observation calculation method with single spectrometer alternation, which does not detect abnormal spectrum data, the technical scheme of the invention can accurately detect and reject SIF calculation values inverted by abnormal spectrum data, reduce SIF calculation value deviation which cannot be judged manually, improve the accuracy and reliability of final SIF calculation values, and provide support for further accurate detection and research of vegetation SIF in ecological field.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings that are needed in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and that other drawings can be obtained according to these drawings without inventive effort for a person skilled in the art.
Fig. 1 shows a flowchart of a SIF anomaly detection method based on PAR data according to one embodiment of the present invention;
fig. 2 shows a block diagram of an apparatus for SIF anomaly detection based on PAR data according to an embodiment of the present invention;
fig. 3 shows a control flow diagram of a central control unit in an apparatus for SIF anomaly detection based on PAR data according to an embodiment of the present invention;
fig. 4 shows a block diagram of a system for SIF anomaly detection based on PAR data according to one embodiment of the present invention.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings, in order to make the objects, technical solutions and advantages of the present invention more apparent. It should be understood that the detailed description is presented by way of example only and is not intended to limit the invention.
PAR means solar photosynthetically active radiation, and the solar photosynthetically active radiation PAR is rapidly measured by a photosynthetically active radiation monitor (hereinafter referred to as PAR monitor). In the application, the solar photosynthetically active radiation PAR provides abnormal data inspection for a spectrometer to measure uplink and downlink radiation paths. The method comprises the following steps: the PAR monitor is arranged beside the optical path collecting unit for measuring the downlink solar radiation upwards, is connected with the central control unit and monitors and collects solar radiation data in real time.
Examples
Referring to fig. 1, an embodiment of the present invention is a SIF anomaly detection method based on PAR data, applied to a central control unit, including:
s1, acquiring spectrum data;
s2, constructing a detection time group comprising at least one detection time unit based on the spectrum data, wherein the detection time unit comprises: an uplink spectrum data acquisition time t, and two downlink spectrum data acquisition times t-1 and t+1 adjacent to the acquisition time t;
s3, acquiring and calculating each acquisition time in any one detection time unit based on the PAR data acquired in advance and the acquisition time recorded in any one detection time unitAnd calculate the absolute value PAR of the difference 1 And PAR 2 Wherein PAR is t PAR data representing acquisition time t; PAR (PAR) t-1 PAR data representing acquisition time t-1; PAR (PAR) t+1 PAR data representing acquisition time t+1;
PAR 1 ＝|PAR t -PAR t-1 | (1)
PAR 2 ＝|PAR t -PAR t+1 | (2)
s4, when PAR 1 Greater than or equal to a predetermined PAR Threshold value And PAR 2 Greater than or equal to a predetermined PAR Threshold value And judging that SIF calculated based on the spectrum data corresponding to the acquisition time of any one detection time unit is an abnormal value and marking.
Explanation:
(one) with respect to spectral data and PAR data
Spectral data: the spectral data acquired from the spectrometer comprises at least one piece of upstream spectral data and two pieces of downstream spectral data at times adjacent to the time of acquisition of the at least one piece of upstream spectral data.
PAR data: the PAR data is collected by a PAR monitor, and the PAR monitor is arranged at a collecting end of a downlink radiation light path.
(II) S4 further comprises: when PAR is 1 And PAR 2 Are all smaller than the preset PAR Threshold value And when determining that the SIF calculated based on the spectrum data corresponding to the acquisition time in any one detection time unit is a normal value, specifically: when PAR is 1 And PAR 2 Are all smaller than the preset PAR Threshold value And acquiring two downlink spectrum data and one uplink spectrum data corresponding to the acquisition time based on the acquisition time in any one detection time unit, and calculating SIF through the average value of the two downlink spectrum data and the one uplink spectrum data.
S4 further includes: correction calculation of the calculated SIF outlier: when PAR is 1 And PAR 2 Any one of which is smaller than a preset PAR Threshold value When using any of the aboveLess than a predetermined PAR in a detection time unit Threshold value And (3) calculating SIF according to the downlink spectrum data and the uplink spectrum data at the corresponding acquisition time of PAR, and finishing SIF abnormal value correction and marking. The central control unit corrects the SIF that needs to be corrected in the finally output SIF data based on the correction flag of the SIF.
After S4, the SIF anomaly detection method based on PAR data further includes removing SIF anomaly values, specifically: the central control unit deletes the SIF outlier in the finally output SIF data based on the SIF outlier flag.
In S4, PAR 1 Greater than or equal to a predetermined PAR Threshold value And PAR 2 Greater than or equal to a predetermined PAR Threshold value The cases of (2) include the following: PAR (PAR) 1 Greater than a predetermined PAR Threshold value And PAR 2 Greater than a predetermined PAR Threshold value ；PAR 1 Greater than a predetermined PAR Threshold value And PAR 2 Equal to the preset PAR Threshold value ；PAR 1 Equal to the preset PAR Threshold value And PAR 2 Greater than a predetermined PAR Threshold value ；PAR 1 Equal to the preset PAR Threshold value And PAR 2 Equal to the preset PAR Threshold value 。
Referring to fig. 2, in one embodiment of the present invention, an apparatus for SIF anomaly detection based on PAR data, the apparatus includes: the system comprises a central control unit, a spectrometer connected with the central control unit, an optical path collecting unit and a PAR monitor; the optical path collecting unit comprises a downlink radiation optical path unit and an uplink radiation optical path unit; the PAR monitor is connected with the central control unit and is arranged at the collecting end of the downlink radiation light path.
The central control unit is an industrial personal computer and is responsible for data recording, control, acquisition, processing and storage.
And the spectrometer is responsible for measuring the light path collected by the light path collecting unit, and the measurement result is collected and processed by the central control unit.
And (III) the optical path collecting unit comprises a furcated optical fiber, a cosine corrector and two electronic shutters.
The tail end of the bifurcated optical fiber is connected with the spectrometer, a cosine corrector is arranged at the collecting end of the downlink radiation light path unit of the bifurcated optical fiber, and the cosine corrector is vertically and upwards installed and collects solar downlink radiation in a field of view (FOV) of 180 degrees. The collecting end of the uplink radiation light path unit of the bifurcated optical fiber is a bare optical fiber, points to the canopy and is responsible for collecting the ground vegetation uplink radiation in a field of view (FOV) of 25 degrees.
The two electronic shutters function to block the light path in order to achieve alternate measurements by the spectrometer. One electronic shutter is arranged at the optical fiber path of the downlink radiation optical path unit, the other electronic shutter is arranged at the optical fiber path of the uplink radiation optical path unit, and the two electronic shutters are respectively connected with the central control unit.
(IV) photosynthetic active radiation monitor
In this embodiment, the height of the PAR monitor is the same as the height of the collecting end of the downlink radiation optical path unit, and is connected to the central control unit, so as to monitor and collect PAR data in real time. The PAR monitor is used for realizing synchronous observation of solar downlink radiation and ground vegetation uplink radiation.
The PAR monitor rapidly measures the photosynthetically active radiation of the sun to obtain PAR data, and the central control unit detects whether SIF data calculated by any one spectrum data is abnormal or not based on the PAR data. The method comprises the following steps: based on the PAR data corresponding to the acquisition time in the detection time unit, the central control unit respectively judges the PAR of the acquisition time t t PAR with acquisition time t-1 t-1 Difference absolute value PAR 1 PAR at acquisition time t t PAR with acquisition time t+1 t+1 Difference absolute value PAR 2 And PAR (PAR) Threshold value When PAR is related to 1 、PAR 2 Are all greater than or equal to PAR Threshold value And when the SIF is in a weather condition of rapid light change in the period from t-1 to t+1, SIF calculated based on the spectral data of the period from t-1 to t+1 is an abnormal value and should be removed.
Referring to fig. 3, based on the device for detecting SIF abnormality based on PAR data in the present embodiment, the process of obtaining spectral data corresponding to the acquisition time in one detection time unit and calculating the SIF value reliability is as follows:
(1) The central control unit sends an instruction to control the uplink radiation electronic shutter to be closed, the downlink radiation light path collecting unit collects downlink radiation, the spectrometer measures and collects the 1 st solar downlink radiation spectrum data, the downlink radiation spectrum data is transmitted to the central control unit, and the measurement time t is recorded 1 。
(2) Then the downlink radiation electronic shutter is closed, the uplink radiation electronic shutter is opened, the uplink radiation light path collecting unit collects uplink radiation, the spectrometer measures and collects vegetation uplink radiation spectrum data, the uplink radiation spectrum data is transmitted to the central control unit, and the measuring time t is recorded 2 。
(3) Then the upward radiation electronic shutter is closed, the downward radiation electronic shutter is opened, the downward radiation light path collecting unit collects downward radiation, the spectrometer measures and collects the 2 nd solar downward radiation spectrum data, the downward radiation spectrum data is transmitted to the central control unit, and the measuring time t is recorded 3 。
(4) After the 3 times of spectrum data acquisition is completed, the central control unit performs data acquisition according to t 1 And t 3 Average value of two times of downlink spectrum data at moment and t 2 Calculating t from time uplink radiation spectrum data 1 To t 3 SIF of a time period.
(5) Real-time monitoring and recording PAR data by a PAR monitor, and calculating t by a central control unit 1 And t 2 Absolute value of PAR difference PAR at time instant 1 、t 2 And t 3 Absolute value of PAR difference PAR at time instant 2 。
(6) The central control unit judges SIF data quality conditions according to the following 3 conditions:
a. when PAR is 1 And PAR 2 Are all smaller than the preset PAR Threshold value At the time, then determine t 1 To t 3 The illumination change in the time period is small, and the SIF calculated value is not abnormal.
b. When PAR is 1 And PAR 2 When only 1 value is smaller than the PAR threshold value, judging t 2 Time of day and less than PAR Threshold value The illumination change between acquisition times corresponding to PAR is small, SIF needs correction, acquisitionBy t 2 Time-of-day uplink radiation spectral data and PAR less than Threshold value The SIF is recalculated in response to the PAR of the downlink radiation spectrum data at the acquisition time.
c. When PAR is 1 And PAR 2 Are all greater than or equal to PAR Threshold value Judging the light ray from t in the current detection time unit 1 To t 3 And (3) if the light condition is changed too much, eliminating the abnormal SIF value by taking the SIF calculated based on the spectrum data corresponding to the acquisition time of the detection time unit as the abnormal value.
The final output data of the central control unit comprises complete spectrum original data, complete real-time PAR data and SIF values of abnormal times.
Referring to fig. 4, the technical solution of this embodiment is a system for SIF anomaly detection based on PAR data, where the system includes:
the data acquisition module acquires spectrum data and PAR data;
a detection time framework module that constructs a detection time group comprising at least one detection time cell based on the spectral data, the detection time cell comprising: an uplink spectrum data acquisition time t, and two downlink spectrum data acquisition times t-1 and t+1 adjacent to the acquisition time t;
the searching and calculating module is used for acquiring the PAR at the t moment in the detection time unit i to be detected based on the PAR data t PAR at time t-1 adjacent to said time t t-1 And PAR at time t+1 adjacent to said time t t+1 The method comprises the steps of carrying out a first treatment on the surface of the Calculating the absolute difference PAR according to the formula (3) and the formula (4) 1 And PAR 2 ；
PAR 1 ＝|PAR t -PAR t-1 | (3)
PAR 2 ＝|PAR t -PAR t+1 | (4)
Judgment modules for respectively judging PAR 1 And PAR 2 And preset PAR Threshold value When PAR is related to 1 And PAR 2 Are all greater than or equal to a preset PAR Threshold value If so, determining the number of spectrums corresponding to the detection time unit iThe SIF obtained by calculation is an abnormal value; when PAR is 1 And PAR 2 Any one of which is smaller than a preset PAR Threshold value And if so, judging that SIF calculated based on the spectrum data corresponding to the detection time unit i needs to be corrected.
And the marking module is used for marking the SIF abnormal value based on the SIF abnormal value judging result of the judging module.
The system further comprises: the correction module is used for correcting the SIF required to be corrected based on the result calculated based on the spectrum data corresponding to the detection time unit i, specifically: acquiring a detection time unit corresponding to SIF to be corrected; based on less than a preset PAR Threshold value The PAR corresponding time of the time is obtained to obtain the downlink spectrum data A of the time; acquiring uplink spectrum data B based on the uplink spectrum data acquisition time in the detection time unit corresponding to the SIF to be corrected; and calculating SIF according to the downlink spectrum data A and the uplink spectrum data B, and finishing SIF correction and marking.
One technical solution of this embodiment is that a terminal includes: one or more processors; and a memory for storing one or more programs that, when executed by the one or more processors, cause the one or more processors to implement the method for SIF anomaly detection based on PAR data.
One aspect of this embodiment is a computer-readable medium having stored thereon a computer program that, when executed by a processor, implements the method for SIF anomaly detection based on PAR data.
By adopting the technical scheme disclosed by the invention, the following beneficial effects are obtained: according to the technical scheme, the SIF abnormality detection method based on the PAR data is used for checking whether SIF is abnormal or not by comparing the PAR data results monitored by the PAR monitor, eliminating SIF abnormal values and reducing overall SIF result errors. The method comprises the following steps: the technical scheme of the invention is that the SIF anomaly detection method based on the PAR data is applied to single spectrometer for alternately detecting SIF, whether the irradiance at a certain measurement time t of uplink spectrum data is greatly deviated from the irradiance at the measurement time of the downlink spectrum data adjacent to the front and rear is judged according to the PAR data measured in real time, and if so, the weather condition that the currently judged time t is in rapid light change is indicated. And eliminating SIF calculated values obtained by inversion of the spectrum measured value at the current judging moment t from the SIF calculated value final result. Compared with the prior SIF observation calculation method with single spectrometer alternation, which does not detect abnormal spectrum data, the technical scheme of the invention can accurately detect and reject SIF calculation values inverted by abnormal spectrum data, reduce SIF calculation value deviation which cannot be judged manually, improve the accuracy and reliability of final SIF calculation values, and provide support for further accurate detection and research of vegetation SIF in ecological field.
The foregoing is merely a preferred embodiment of the present invention and it should be noted that modifications and adaptations to those skilled in the art may be made without departing from the principles of the present invention, which is also intended to be covered by the present invention.
Claims (8)
1. A PAR data-based SIF anomaly detection method applied to a central control unit, comprising:
acquiring spectrum data;
constructing a detection time group comprising at least one detection time unit based on the spectral data, the detection time unit comprising: an uplink spectrum data acquisition time t, and two downlink spectrum data acquisition times t-1 and t+1 adjacent to the acquisition time t;
based on the pre-acquired PAR data and the acquisition time recorded by any one detection time unit, acquiring and calculating the PAR data of each acquisition time in any one detection time unit, and calculating the absolute difference value PAR 1 And PAR 2 Wherein PAR is t PAR data representing acquisition time t; PAR (PAR) t-1 PAR data representing acquisition time t-1; PAR (PAR) t+1 PAR data representing acquisition time t+1;
PAR 1 ＝|PAR t -PAR t-1 | (1)
PAR 2 ＝|PAR t -PAR t+1 | (2)
when PAR is 1 Greater than or equal to a predetermined PAR Threshold value And PAR 2 Greater than or equal to a predetermined PAR Threshold value If yes, determining SIF calculated based on spectrum data corresponding to the acquisition time of any one detection time unit as an abnormal value and marking;
when PAR is 1 And PAR 2 Any one of which is smaller than a preset PAR Threshold value If yes, adopting any detection time unit to be smaller than the preset PAR Threshold value Calculating SIF according to the downlink spectrum data and the uplink spectrum data at the moment corresponding to PAR, and finishing SIF abnormal value correction and marking;
when PAR is 1 And PAR 2 Are all smaller than the preset PAR Threshold value And judging that the SIF calculated by adopting the spectrum data corresponding to any one detection time unit is a normal value.
2. The SIF anomaly detection method based on PAR data of claim 1, wherein the PAR data is collected by a PAR monitor, the PAR monitor being disposed at a downstream radiation path collection end.
3. The method for detecting SIF anomalies based on PAR data according to claim 1, further comprising culling SIF anomalies, specifically: the central control unit deletes the SIF outlier in the finally output SIF data based on the SIF outlier flag.
4. An apparatus for implementing the PAR data-based SIF anomaly detection method of any one of claims 1 to 3, the apparatus comprising a central control unit, a spectrometer connected to the central control unit, and an optical path collection unit; the optical path collecting unit comprises a downlink radiation optical path unit and an uplink radiation optical path unit; characterized in that the device further comprises: and the photosynthetic active radiation monitor is connected with the central control unit and is arranged at the collecting end of the downlink radiation light path.
5. A system for implementing the PAR data-based SIF anomaly detection method of any one of claims 1 to 3, the system comprising:
the data acquisition module acquires spectrum data and PAR data;
a detection time framework module that constructs a detection time group comprising at least one detection time cell based on the spectral data, the detection time cell comprising: an uplink spectrum data acquisition time t, and two downlink spectrum data acquisition times t-1 and t+1 adjacent to the acquisition time t;
the searching and calculating module is used for acquiring the PAR at the t moment in the detection time unit i to be detected based on the PAR data t PAR at time t-1 adjacent to said time t t-1 And PAR at time t+1 adjacent to said time t t+1 The method comprises the steps of carrying out a first treatment on the surface of the Calculating the absolute difference PAR according to the formula (3) and the formula (4) 1 And PAR 2 ；
PAR 1 ＝|PAR t -PAR t-1 | (3)
PAR 2 ＝|PAR t -PAR t+1 | (4)
Judgment modules for respectively judging PAR 1 And PAR 2 And preset PAR Threshold value When PAR is related to 1 And PAR 2 Are all greater than or equal to a preset PAR Threshold value If so, judging that SIF calculated based on the spectrum data corresponding to the detection time unit i is an abnormal value; when PAR is 1 And PAR 2 Any one of which is smaller than a preset PAR Threshold value When the SIF is detected, determining that SIF calculated based on the spectrum data corresponding to the detection time unit i needs to be corrected; when PAR is 1 And PAR 2 Are all smaller than the preset PAR Threshold value When the SIF is detected, determining that SIF obtained by calculation by adopting spectrum data corresponding to any detection time unit is a normal value;
and the marking module marks the SIF abnormal value based on the SIF abnormal value judging result of the judging module.
6. The system of the SIF anomaly detection method based on PAR data of claim 5, the system further comprising: the correction module is used for correcting the SIF required to be corrected based on the result calculated based on the spectrum data corresponding to the detection time unit i, specifically:
acquiring a detection time unit corresponding to SIF to be corrected;
based on less than a preset PAR Threshold value The PAR corresponding time of the time is obtained to obtain the downlink spectrum data A of the time;
acquiring uplink spectrum data B based on the uplink spectrum data acquisition time in the detection time unit corresponding to the SIF to be corrected;
and calculating SIF according to the downlink spectrum data A and the uplink spectrum data B, and finishing SIF correction and marking.
7. A terminal, comprising: one or more processors; a memory for storing one or more programs that, when executed by the one or more processors, cause the one or more processors to implement the PAR data based SIF anomaly detection method of any one of claims 1-3.
8. A computer-readable medium, on which a computer program is stored, characterized in that the computer program, when executed by a processor, implements the SIF anomaly detection method based on PAR data as claimed in any one of claims 1 to 3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211441181.2A CN115753712B (en) | 2022-11-17 | 2022-11-17 | SIF anomaly detection method, device, system, terminal and medium based on PAR data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211441181.2A CN115753712B (en) | 2022-11-17 | 2022-11-17 | SIF anomaly detection method, device, system, terminal and medium based on PAR data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115753712A CN115753712A (en) | 2023-03-07 |
| CN115753712B true CN115753712B (en) | 2023-05-30 |
Family
ID=85372762
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211441181.2A Active CN115753712B (en) | 2022-11-17 | 2022-11-17 | SIF anomaly detection method, device, system, terminal and medium based on PAR data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115753712B (en) |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5410323B2 (en) * | 2009-11-27 | 2014-02-05 | 国立大学法人九州大学 | Optical vegetation index sensor |
| CN106770142A (en) * | 2017-02-28 | 2017-05-31 | 成都学院 | The Chlorophyll Fluorescence research method and device of a kind of Jiu Zhaigou water plant |
| US11287381B2 (en) * | 2017-06-30 | 2022-03-29 | Ut-Battelle, Llc | System and method for measuring sun-induced chlorophyll fluorescence |
| CN108593611B (en) * | 2018-04-25 | 2020-10-02 | 南京大学 | Double-channel automatic observation method for chlorophyll fluorescence of vegetation |
| CN111638195A (en) * | 2020-05-12 | 2020-09-08 | 中国科学院空天信息创新研究院 | Drought monitoring method and device |
| CN112557306A (en) * | 2020-12-07 | 2021-03-26 | 无锡钱荣分析仪器有限公司 | Abnormal spectrum eliminating method of full-spectrum type straightness spectrometer |
-
2022
- 2022-11-17 CN CN202211441181.2A patent/CN115753712B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN115753712A (en) | 2023-03-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Rahaman et al. | PV system defects identification using Remotely Piloted Aircraft (RPA) based infrared (IR) imaging: A review | |
| CN109034252B (en) | Automatic identification method for monitoring data abnormity of air quality station | |
| CN102589705B (en) | Single-blackbody temperature-controlled MRTD (Minimum Resolvable Temperature Difference) field online automatic detection device and method for thermal imager | |
| CN102680481B (en) | Detection method for cotton fiber impurities | |
| CN105739575A (en) | Greenhouse vegetable environmental parameter data fusion method, data fusion device, and data fusion system | |
| CN110567964A (en) | method and device for detecting defects of power transformation equipment and storage medium | |
| CN116429988B (en) | Dynamic remote sensing monitoring device for ocean plant carbon sink | |
| CN116046692B (en) | Soil heavy metal pollution monitoring method and device based on hyperspectrum | |
| Polonik et al. | Comparison of gas analyzers for eddy covariance: Effects of analyzer type and spectral corrections on fluxes | |
| US20230042106A1 (en) | System and method for the statistical analysis of images of photovoltaic panels | |
| CN114202675A (en) | Crop growth parameter determination method and system based on multispectral remote sensing image of unmanned aerial vehicle | |
| CN105890873A (en) | Apparatus and method for blind pixel detection of thermal infrared hyperspectral imager | |
| Deceglie et al. | Real-time series resistance monitoring in PV systems without the need for I–V curves | |
| CN114034826A (en) | Production environment monitoring system for freeze-drying cubilose process based on data analysis | |
| Payne et al. | Evaluating the accuracy of point spread function deconvolutions applied to luminescence images | |
| CN115753712B (en) | SIF anomaly detection method, device, system, terminal and medium based on PAR data | |
| CN112508995B (en) | Coal flow real-time dynamic metering method based on TOF camera | |
| Nollas et al. | Quality control procedure for 1-minute pyranometric measurements of global and shadowband-based diffuse solar irradiance | |
| CN116757332B (en) | Leaf vegetable yield prediction method, device, equipment and medium | |
| Puranik et al. | Progress in module level quantitative electroluminescence imaging of crystalline silicon PV module: A review | |
| Driesse et al. | Spectrally selective sensors for PV system performance monitoring | |
| EP3654013B1 (en) | Correction method of particle size distribution measuring apparatus, calibration method of particle size distribution measuring apparatus, particle size distribution measuring apparatus, and gas analysis system | |
| CN110608800A (en) | Method, device and system for measuring solar radiation illuminance of atmospheric absorption band | |
| Denz et al. | Quantitative assessment of the power loss of silicon PV modules by IR thermography and its practical application in the field | |
| Kilper et al. | A new generation of PV Monitoring system with high-grade remote diagnostics based on module level monitoring and integrated yield simulation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |