CN108734596A - A kind of time series interpolating method of aquaculture water quality monitoring missing data - Google Patents
A kind of time series interpolating method of aquaculture water quality monitoring missing data Download PDFInfo
- Publication number
- CN108734596A CN108734596A CN201810311429.0A CN201810311429A CN108734596A CN 108734596 A CN108734596 A CN 108734596A CN 201810311429 A CN201810311429 A CN 201810311429A CN 108734596 A CN108734596 A CN 108734596A
- Authority
- CN
- China
- Prior art keywords
- time series
- principal component
- water quality
- sequence
- monitoring data
- 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.)
- Pending
Links
- 238000012544 monitoring process Methods 0.000 title claims abstract description 46
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 27
- 238000009360 aquaculture Methods 0.000 title claims abstract description 17
- 244000144974 aquaculture Species 0.000 title claims abstract description 17
- 238000000034 method Methods 0.000 title claims abstract description 14
- 239000011159 matrix material Substances 0.000 claims abstract description 16
- 238000010276 construction Methods 0.000 claims abstract description 4
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 4
- 230000000694 effects Effects 0.000 abstract description 2
- 230000007704 transition Effects 0.000 abstract description 2
- 230000007613 environmental effect Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000036760 body temperature Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A20/00—Water conservation; Efficient water supply; Efficient water use
- Y02A20/152—Water filtration
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Business, Economics & Management (AREA)
- Data Mining & Analysis (AREA)
- Agronomy & Crop Science (AREA)
- Algebra (AREA)
- General Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computing Systems (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Farming Of Fish And Shellfish (AREA)
Abstract
The present invention relates to the time series interpolating methods that a kind of aquaculture water quality monitors missing data, including step to have:The first step:Choose the identical water quality monitoring time series containing missing data of two groups of length;Second step:Establish delay matrix and calculating matrix inner product;Third walks:Time series Toeplitz of the construction with missing data postpones autocorrelation matrix;4th step:Singular value decomposition is made to autocorrelation matrix, obtains feature vector;5th step:Calculate principal component vector;6th step:Time series is reconstructed using each principal component, obtains interpolation sequence.The present invention decomposes the time series containing missing values, reconstructs to obtain the time series after interpolation, the interpolating method can hold the transition information of time series comprehensively, achieve good comprehensive interpolation effect, the continuity that ensure that data time series, to realize that the prediction of water monitoring data time series provides effective means.
Description
Technical field
The invention belongs to aquaculture water quality monitoring technical fields, especially set and a kind of aquaculture water quality monitors missing number
According to time series interpolating method.
Background technology
Aquatic environment monitoring technology is used widely in intensive aquaculture industry, and it is strong to be greatly promoted aquatic products
The development of health aquaculture.Aquatic environment information monitoring is the key link for realizing high intensity aquaculture, by monitoring water body temperature
Degree, pH, dissolved oxygen etc. have aquatic products growing environment the water quality factor of significant impact, and best growth ring is provided for aquatic products
Border.Using automatic on-line monitoring system, water quality monitoring sensor automatic collection aquatic environment data are carried out, realizes and covers extensively
Lid, automatic measurement, real-time Transmission represent the developing direction of aquaculture water quality environmental monitoring;But there is measurements at this stage
The disadvantage that precision is low, data easily lack becomes the bottleneck of such method development, hinders the extensive use in aquaculture.
It is insufficient due to sampling in aquaculture water quality environmental monitoring, system mistake is monitored, or other
Reason frequently can lead to error or the missing of monitoring data.In this case, it is handled in progress aquatic environment big data
With analysis when, it is common practice to by with missing or mistake observation data collective delete, although do so ensure that it is used
The reliability of data, but cause the loss of a large amount of useful information.It realizes to the pre- of water monitoring data time series
Survey, it is necessary to the continuity for ensureing time series, for this problem, need combine water monitoring data the characteristics of and difference
Missing degree, be specially designed for the interpolating method of aquaculture water quality environmental monitoring data.
Invention content
It is an object of the invention to overcome the deficiencies in the prior art, and provide a kind of aquaculture water quality monitoring missing data
Time series interpolating method,
The present invention solves its technical problem and following technical scheme is taken to realize:
A kind of time series interpolating method of aquaculture water quality monitoring missing data, including steps are as follows:
The first step:Choose the identical water quality monitoring time series containing missing data of two groups of length, time series vector
Respectively A=(a1,a2,…aN), B=(b1,b2,…bN), sequence length N is established corresponding with nested dimension M, M < N/2
Delay matrix is respectively:
Second step:The missing data component set of above-mentioned two groups of water monitoring data time serieses is:And haveCalculating matrix inner product
Third walks:Time series Toeplitz delay autocorrelation matrix S of the construction with missing data
4th step:Singular value decomposition, S=E λ E are made to autocorrelation matrix ST, obtain eigenvalue λkAnd corresponding feature vector Ek,
Wherein 1≤k≤M, characteristic value are arranged as λ in descending order1≥λ2≥…λM≥0;
5th step:Calculate principal component vector, wherein i-th of value of k-th of principal component of monitoring data time series A be:
I-th of value of k-th of principal component of monitoring data time series B be:
6th step:Time series is reconstructed using each principal component,
I-th of element representation of sequence of k-th of principal component reconstruct of monitoring data time series A be:Monitoring data
I-th of element representation of sequence of k-th of principal component reconstruct of time series B be:
As M≤i≤N-M+1,
As 1≤i≤M-1,
As N-M+2≤i≤N,
I-th of the element x of sequence reconstructed using the preceding K principal component of monitoring data time series AiIt is expressed as:
I-th of the element yi of sequence reconstructed using the preceding K principal component of monitoring data time series B is expressed as:
The advantages and positive effects of the present invention are:
1, the present invention decomposes the time series containing missing values, reconstructs to obtain the time series after interpolation, this is inserted
Compensating method can hold the transition information of time series comprehensively, achieve good comprehensive interpolation effect.
2, the present invention overcomes in aquaculture water quality environmental monitoring, since sampling is insufficient, system mistake is monitored, or
Person is the loss of a large amount of useful informations caused by some other reasons, ensure that the continuity of data time series, is real
The prediction of existing water monitoring data time series provides effective means.
Description of the drawings
Fig. 1 is the logical procedure diagram of the method for the present invention.
Specific implementation mode
The embodiment of the present invention is further described below in conjunction with attached drawing, it is emphasized that, following implementation is to say
Bright property, without being restrictive, it cannot be used as limitation of the invention in this embodiment.
A kind of time series interpolating method of aquaculture water quality monitoring missing data, as shown in Figure 1, such as including step
Under:
The first step:Choose the identical water quality monitoring time series containing missing data of two groups of length, time series vector
Respectively A=(a1,a2,…aN), B=(b1,b2,…bN), sequence length N is established corresponding with nested dimension M, M < N/2
Delay matrix is respectively:
Second step:The missing data component set of above-mentioned two groups of water monitoring data time serieses is:And haveCalculating matrix inner product
Third walks:Time series Toeplitz delay autocorrelation matrix S of the construction with missing data
4th step:Singular value decomposition, S=E λ E are made to autocorrelation matrix ST, obtain eigenvalue λkAnd corresponding feature vector Ek,
Wherein 1≤k≤M, characteristic value are arranged as λ in descending order1≥λ2≥…λM≥0;
5th step:Calculate principal component vector, wherein i-th of value of k-th of principal component of monitoring data time series A be:
I-th of value of k-th of principal component of monitoring data time series B be:
6th step:Time series is reconstructed using each principal component,
I-th of element representation of sequence of k-th of principal component reconstruct of monitoring data time series A be:Monitoring data
I-th of element representation of sequence of k-th of principal component reconstruct of time series B be:
As M≤i≤N-M+1,
As 1≤i≤M-1,
As N-M+2≤i≤N,
I-th of the element x of sequence reconstructed using the preceding K principal component of monitoring data time series AiIt is expressed as:
I-th of the element y of sequence reconstructed using the preceding K principal component of monitoring data time series BiIt is expressed as:
Claims (1)
1. a kind of time series interpolating method of aquaculture water quality monitoring missing data, it is characterised in that including steps are as follows:
The first step:Choose the identical water quality monitoring time series containing missing data of two groups of length, time series vector difference
For A=(a1,a2,…aN), B=(b1,b2,…bN), sequence length N establishes corresponding time lag with nested dimension M, M < N/2
Matrix is respectively:
Second step:The missing data component set of above-mentioned two groups of water monitoring data time serieses is:And haveCalculating matrix inner product
Third walks:Time series Toeplitz delay autocorrelation matrix S of the construction with missing data
4th step:Singular value decomposition, S=E λ E are made to autocorrelation matrix ST, obtain eigenvalue λkAnd corresponding feature vector Ek, wherein
1≤k≤M, characteristic value are arranged as λ in descending order1≥λ2≥…λM≥0;
5th step:Calculate principal component vector, wherein i-th of value of k-th of principal component of monitoring data time series A be:
I-th of value of k-th of principal component of monitoring data time series B be:
6th step:Time series is reconstructed using each principal component,
I-th of element representation of sequence of k-th of principal component reconstruct of monitoring data time series A be:The monitoring data time
I-th of element representation of sequence of k-th of principal component reconstruct of sequence B be:
As M≤i≤N-M+1,
As 1≤i≤M-1,
As N-M+2≤i≤N,It utilizes
I-th of element x of sequence of the preceding K principal component reconstruct of monitoring data time series AiIt is expressed as:
I-th of the element y of sequence reconstructed using the preceding K principal component of monitoring data time series BiIt is expressed as:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810311429.0A CN108734596A (en) | 2018-04-09 | 2018-04-09 | A kind of time series interpolating method of aquaculture water quality monitoring missing data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810311429.0A CN108734596A (en) | 2018-04-09 | 2018-04-09 | A kind of time series interpolating method of aquaculture water quality monitoring missing data |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108734596A true CN108734596A (en) | 2018-11-02 |
Family
ID=63940687
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810311429.0A Pending CN108734596A (en) | 2018-04-09 | 2018-04-09 | A kind of time series interpolating method of aquaculture water quality monitoring missing data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108734596A (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002039254A1 (en) * | 2000-11-09 | 2002-05-16 | Spss Inc. | System and method for building a time series model |
CN102446302A (en) * | 2011-12-31 | 2012-05-09 | 浙江大学 | Data preprocessing method of water quality prediction system |
CN103577694A (en) * | 2013-11-07 | 2014-02-12 | 广东海洋大学 | Aquaculture water quality short-time combination forecast method on basis of multi-scale analysis |
CN107577649A (en) * | 2017-09-26 | 2018-01-12 | 广州供电局有限公司 | The interpolation processing method and device of missing data |
-
2018
- 2018-04-09 CN CN201810311429.0A patent/CN108734596A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002039254A1 (en) * | 2000-11-09 | 2002-05-16 | Spss Inc. | System and method for building a time series model |
CN102446302A (en) * | 2011-12-31 | 2012-05-09 | 浙江大学 | Data preprocessing method of water quality prediction system |
CN103577694A (en) * | 2013-11-07 | 2014-02-12 | 广东海洋大学 | Aquaculture water quality short-time combination forecast method on basis of multi-scale analysis |
CN107577649A (en) * | 2017-09-26 | 2018-01-12 | 广州供电局有限公司 | The interpolation processing method and device of missing data |
Non-Patent Citations (1)
Title |
---|
孙东: "赤潮多源监测数据处理与综合预测预报方法研究", 《万方数据》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bakun et al. | Environmental inputs to fishery population models for eastern boundary current regions | |
Mohan et al. | Two new indices for stomach content analysis of fishes | |
Munoz-Benavent et al. | Automatic Bluefin Tuna sizing using a stereoscopic vision system | |
CN106442233A (en) | Novel mixed method for estimating particle size of seawater suspended particulate matter by means of satellite data | |
CN116091944A (en) | Aquatic ecology investigation method for aquatic vegetation coverage | |
Aron et al. | A DESCRIPTION OF A DISCRETE DEPTH PLANKTON SAMPLER WITH SOME NOTES ON THE TOWING BEHAVIOR OF AB‐FOOT ISAACS‐KIDD MID‐WATER TRAWL AND A ONE‐METER RING NET | |
Ionescu et al. | DeePS at: A deep learning model for prediction of satellite images for nowcasting purposes | |
CN108734596A (en) | A kind of time series interpolating method of aquaculture water quality monitoring missing data | |
Nouchi et al. | Effects of non-uniform vertical constituent profiles on remote sensing reflectance of oligo-to mesotrophic lakes | |
CN111723479B (en) | Process-oriented eutrophic lake algae total quantity remote sensing estimation method | |
CN105677770B (en) | Self-adaptive sampling method for monitoring offshore marine environment data | |
CN110308254B (en) | Accurate monitoring method for low dissolved oxygen | |
Salminen et al. | Varve distribution reveals spatiotemporal hypolimnetic hypoxia oscillations during the past 200 years in Lake Lehmilampi, eastern Finland | |
Shelton et al. | Among-year variation in growth of Antarctic krill Euphausia superba based on length-frequency data | |
Fulanda et al. | Determination of effective tow duration for estimation of fish and shellfish abundance in deeper-water bottom trawl surveys | |
Dawson et al. | Electrochemical noise--some new developments in corrosion monitoring | |
CN206282358U (en) | A kind of cross section of river two-dimensional display system | |
CN214122000U (en) | Floating algae detection equipment | |
Evangelista et al. | Comparison of diversity indices applied to macrophyte incidence-based data | |
Gerdeaux et al. | Inferring occurrence of growth checks in pike (Esox lucius) scales by using sequential isotopic analysis of otoliths | |
Sidik et al. | The identification of outliers in wrapped normal data by using Ga statistics | |
Koranteng et al. | Population parameters, biannual cohorts, and assessment in the Pagellus bellottii (Sparidae) fishery off Ghana | |
Komeyama et al. | Body size estimation method for seasonally growing farmed yellowtail Seriola quinqueradiata in an aquaculture net cage using a stereo camera | |
Haotian et al. | An augmented reality visualization algorithm of the lightweight ROV | |
Herman et al. | Meso-scale spatial distribution of plankton: co-evolution of concepts and instrumentation |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181102 |
|
RJ01 | Rejection of invention patent application after publication |