CN112579687B - Marine environment monitoring data compressed sensing online reconstruction method - Google Patents
Marine environment monitoring data compressed sensing online reconstruction method Download PDFInfo
- Publication number
- CN112579687B CN112579687B CN202011409107.3A CN202011409107A CN112579687B CN 112579687 B CN112579687 B CN 112579687B CN 202011409107 A CN202011409107 A CN 202011409107A CN 112579687 B CN112579687 B CN 112579687B
- Authority
- CN
- China
- Prior art keywords
- data
- matrix
- reconstruction
- solving
- compressed sensing
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/26—Visual data mining; Browsing structured data
Landscapes
- Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
The invention discloses an online compressed sensing reconstruction method for marine environment monitoring data. Aiming at the problem that the reconstruction performance of time-varying ocean monitoring data is not high in the traditional compressed sensing method, the ocean measurement data online reconstruction method based on the low-rank regularized sequential compressed sensing is provided. The method includes the steps of firstly adding low-rank performance of data on a spatial structure into an optimization algorithm, constructing a low-rank regular term by combining existing historical data through a sliding window mechanism, establishing a data fidelity term according to the condition that data in overlapping regions of previous and next moments are equal, and finally solving the reconstruction optimization algorithm based on an alternating direction multiplier method. The method can improve the compressed sensing reconstruction performance of the time-varying ocean monitoring data.
Description
Technical Field
The invention relates to a compressed sensing technology, in particular to an online compressed sensing reconstruction method for marine environment monitoring data, and belongs to the technical field of compressed sensing reconstruction.
Background
The ocean is the foundation on which the human beings rely for survival and development, along with the development of productivity and the progress of social civilization, the world economy, science and technology and military strength are improved, the resource demand of the human society is increased day by day, and because 71 percent of the surface of the ocean floor ball is rich in resources and is not fully developed, the vision of various countries is changed from land to ocean, the enthusiasm of developing and utilizing the ocean is high and the conflict of the ocean rights of the ocean is increased. In this context, the problem of monitoring marine environments is beginning to be of widespread concern. Because the scale range of data acquired by means such as satellite remote sensing and radar is too large and the resolution is too low, the sensor network is called as a main means for acquiring marine environment monitoring data.
In a real scene, a sensor network is formed by deploying a plurality of sensors in a certain ocean monitoring area, the same object is continuously sensed, various environment monitoring data are acquired and transmitted to a monitoring center through a certain communication means, and different types of sensing applications are supported. Compared with the land environment, the marine environment has advantages in depth and breadth, but the advantages place higher demands on the deployment of the marine sensor network: on one hand, the ocean range is too large, so that the number of deployed nodes is too large, and the construction cost is too high; on the other hand, how to transmit the collected mass information timely and effectively is also limited by node energy consumption and communication cost. And the sensor network has the characteristics of more sensors, dynamic change of a topological structure, easiness in damage of sensor nodes, limited storage capacity and the like. This determines that a reasonable solution must be found to solve the contradiction between limited energy and data transmission during the operation of the sensor network.
In recent years, Compressed Sensing (CS) provides a new idea for solving the above problems. According to the CS theory, if a signal is sparse in a certain domain, it can be reconstructed from linear measurements of very few raw data as long as certain conditions are met. Bajwa et al first proposed the introduction of the CS theory into the efficient data collection problem of sensor networks, after which, researchers began conducting a great deal of research into the use of CS theory in sensor networks. The existing algorithm can greatly reduce the communication energy consumption overhead and realize the acquisition, transmission and application of data. None of these algorithms take into account correlations created by dynamic changes in the sensor network data collection process. The sequential compression acquisition and progressive reconstruction of the sensor data stream, however, in the process of collecting the sensor time series data, only the dynamic evolution relation of the time domain stream signal before and after is considered, the spatial structure correlation among multiple sensors is not considered, and the reconstruction performance needs to be improved.
Disclosure of Invention
The invention aims to design a compressed sensing-based marine environment monitoring data online reconstruction method, and aims to solve the problem that the reconstruction performance of the traditional compressed sensing method for processing time-varying marine monitoring data is not high. Firstly, designing a sliding window mechanism, combining the existing historical data, and fully utilizing the low-rank property of a matrix to construct a low-rank matrix; then, according to the rule that the data of the overlapping area are equal, an error correction term of the overlapping area is added; and finally, designing a solving method based on an alternating direction multiplier (ADMM) in the process of solving the reconstruction optimization algorithm.
In order to achieve the purpose, the technical scheme adopted by the invention is a compressed sensing online reconstruction method of marine environment monitoring data, which comprises the following steps:
step 2, establishing a compressed sensing online reconstruction model as follows:
s.t.S(t)=Ψ(t)-1vec(X(t)),Z(t)=X(t)
wherein Y (t) ═ Y(t-W):(t-1),yt],Representing an observation matrix, m representing the number of observations, N representing the length of the signal, and m < N,representing all sensor node data at time t, λ1And λ2Representing regularization parameters, | Z (t) | non-woven calculation*Represents the nuclear norm of the matrix z (t),new matrix of first (W-1) rows representing the spatio-temporal basis Ψ (t), ΨT(W-1) is the time domain basis ΨTOf the first (W-1) row ofSA spatial basis is represented by a number of spatial groups,showing the estimated values of the overlapping areas of the two moments before and after the moment (t-1), Z (t) showing an intermediate variable, S (t) being a sparse representation matrix of X (t),estimates of X (t), S (t), and Z (t), respectively, vec (-) indicating vectorization of the matrix; step 3, utilizing the input ytRapid reconstruction of data x at time ttAnd with input X(t-w):(t-1)Build the data matrix X in the window1(t) initializing parameters as initial data;
step 3, utilizing the input ytRapid reconstruction of data x at time ttAnd with input X(t-w):(t-1)Build the data matrix X in the window1(t) initializing parameters as initial data;
step 4, solving the reconstruction model by adopting an alternating direction multiplier method, rewriting the reconstruction model into an augmented Lagrange form,
a and B represent Lagrange multipliers, alpha and beta are penalty parameters, the problem is decomposed into a plurality of subproblems, and iterative solution is respectively carried out on the subproblems to obtain an original data estimation value at the time t
And 5, when t is t +1, the sliding window moves forward, data at the time t moves in, data at the time t-w moves out, and an observation data matrix Y (t +1) at the time t +1 is constructed. If T is equal to T, outputting, otherwise, returning to the step 2;
compared with the prior art, the invention has the following advantages:
1. the low-rank constraint is added into an optimization algorithm, a low-rank regular term is constructed by combining the existing historical data, and the data reconstruction can be optimally constrained by fully utilizing the internal structural features of the data;
2. the condition that data in the overlapped area of a matrix formed by the time-varying data in the sliding window are equal is fully utilized, a data fidelity item is established, the characteristics of the time-varying data are added into an optimization constraint, and the reconstruction precision is further improved.
Drawings
FIG. 1 is a schematic flow chart of an online compressed sensing reconstruction method for marine environmental data according to the present invention;
FIG. 2 is a schematic view of a sliding window data acquisition;
fig. 3 is a schematic diagram of an observation process of sequential compressed sensing.
Detailed Description
The present invention will be described in more detail with reference to the following examples and the accompanying drawings.
As shown in fig. 1, an embodiment of the invention discloses a method for compressed sensing online reconstruction of marine environmental monitoring data, which comprises the following specific steps:
step 2, establishing a compressed sensing online reconstruction model as follows:
s.t.S(t)=Ψ-1vec(X(t)),Z(t)=X(t)
wherein, yN(t)∈Rm×1,Φ(t)∈Rm×NM represents the number of observations, N represents the magnitude of the signal, and m < N, xN(t)∈RN×1Represents all sensor node data at time t, λ1And λ2Representing regularization parameters, | X | | non-woven phosphor*Representing the kernel norm of matrix X, w represents the window size,ΨT(w-1) represents taking the time domain base psiTOf the first (w-1) row ofSA spatial basis is represented by a number of spatial groups,an estimated value representing the overlapping area calculated at time (t-1);
step 3, utilizing the transfusionIn ytRapid reconstruction of data x at time ttAnd with input X(t-w):(t-1)Build data matrix X in window1(t) initializing parameters as initial data;
step 4, solving the reconstruction model by adopting an alternating direction multiplier method, rewriting the reconstruction model into an augmented Lagrange form,
a and B represent Lagrange multipliers, alpha and beta are penalty parameters, the problem is decomposed into a plurality of subproblems, and iterative solution is respectively carried out on the subproblems to obtain an original data estimated value at the time tThe optimization problem is solved by using an alternating direction multiplier method, and an original reconstruction model is divided into a plurality of subproblems to carry out iterative solution. Knowing the data for the kth iteration, the kth (k +1) iteration is updated
Wherein the iteration parameter is Ak+1=Ak+α(Sk+1(t)-Ψvec(Xk+1(t))),Bk+1=Bk+β(Zk+1(t)-Xk+1(t)); all matrices with subscript k (or k +1) represent the values of the original matrix without subscript at the kth (or k +1) iteration. If K > K (K represents the preset maximum iteration number), stopping iteration, and stopping iteration at the momentIs equal to matrix XkThe last column of (t);
the solving method of the neutron problem X (t) in the iterative process comprises the following specific steps:
(1) first, the gradient of the first sub-problem is calculated
Derivative and let the derivative equal zero to obtain
WhereinIn order to further reduce the calculated amount of matrix inverse operation, a conjugate gradient algorithm is adopted to solve an optimal solution;
(2) the solving method of the sub-problem S (t) in the iteration process comprises the following specific steps:
this equation is a standard linear least squares problem, the optimal solution being
Wherein G ═ 2 λ1Ψ(w-1)TΨ(w-1)+αI);
(3) The solving method of the neutron problem Z (t) in the iterative process comprises the following specific steps:
solving Z by singular value shrinkage operator methodk+1(t)
Zk+1(t)=shrink(Xk+1(t)-Bk/β,λ2/β)
Where shrink (H, τ) refers to a function that soft-thresholds the singular values of matrix H at threshold τ.
And step 5, when t is t +1, the sliding window moves forward, the data at the time t moves in, and the data at the time t-w moves out, so as to construct an observation data matrix Y (t +1) at the time t +1, as shown in fig. 2 and fig. 3. If T is equal to T, outputting, otherwise, returning to the step 2;
the above is a detailed description of the present invention in the following detailed description, which is not intended to limit the present invention in form, so that any person skilled in the art can make changes, modifications and substitutions based on the principle and spirit of the present invention after reading the present invention, and the present invention shall fall into the protection scope of the present invention.
Claims (5)
1. A marine environment monitoring data online reconstruction method based on sequential compressed sensing is characterized by comprising the following steps:
step 1, inputting observation data y at t momentt(t-W) raw data X in time period (t-1)(t-W):(t-1)And its observation data Y(t-W):(t-1)Constructing an observation data matrix Y (t), wherein W represents the size of a sliding window;
step 2, establishing a compressed sensing online reconstruction model as follows:
s.t.S(t)=Ψ(t)-1vec(X(t)),Z(t)=X(t)
wherein Y (t) ═ Y(t-W):(t-1),yt],Representing an observation matrix, m representing the number of observations, N representing the length of the signal, and m < N,represents all sensor node data at time t, λ1And λ2Representing a regularization parameter, | Z (t) | ventilated computing*Represents the kernel norm of the matrix z (t),new matrix of first (W-1) rows representing the space-time basis Ψ (t), ΨT(W-1) is the time domain basis ΨTOf the first (W-1) row ofSA spatial basis is represented by a number of spatial groups,showing the estimated values of the overlapping areas of the two moments before and after the moment (t-1), Z (t) showing an intermediate variable, S (t) being a sparse representation matrix of X (t),estimated values of x (t), s (t) and z (t), respectively, vec (-) indicating vectorization of the matrix;
step 3, utilizing the input ytRapid reconstruction of data x at time ttAnd with input X(t-w):(t-1)Build the data matrix X in the window1(t) initializing parameters as initial data;
step 4, solving the reconstruction model by adopting an alternating direction multiplier method, rewriting the reconstruction model into an augmented Lagrange form,
a and B represent Lagrange multipliers, alpha and beta are penalty parameters, the problem is decomposed into a plurality of subproblems, and iteration solving is carried out respectively to obtain x at t momenttEstimated value
And 5, if T is equal to T, outputting, otherwise, returning to the step 2.
2. The method according to claim 1, wherein the method for solving the reconstructed model in step 4 is specifically: dividing the original reconstruction model into a plurality of subproblems to carry out iterative solution, knowing the data of the kth iteration, updating the kth (k +1) iteration process
3. The method according to claim 2, wherein the solving method of the neutron problem x (t) in the iterative process is specifically:
first, the gradient of the first sub-problem is calculated
Derivative and make the derivative equal to zero to obtain
5. The method according to claim 2, wherein the solving method of the neutron problem z (t) in the iterative process is specifically:
solving Z by singular value shrinkage operator methodk+1(t)
Zk+1(t)=shrink(Xk+1(t)-Bk/β,λ2/β)
Where shrnk (H, τ) refers to a function that soft-thresholds singular values of matrix H at threshold τ.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011409107.3A CN112579687B (en) | 2020-12-04 | 2020-12-04 | Marine environment monitoring data compressed sensing online reconstruction method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011409107.3A CN112579687B (en) | 2020-12-04 | 2020-12-04 | Marine environment monitoring data compressed sensing online reconstruction method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112579687A CN112579687A (en) | 2021-03-30 |
CN112579687B true CN112579687B (en) | 2022-07-15 |
Family
ID=75127270
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011409107.3A Active CN112579687B (en) | 2020-12-04 | 2020-12-04 | Marine environment monitoring data compressed sensing online reconstruction method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112579687B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113408644B (en) * | 2021-07-02 | 2023-07-14 | 南京信息工程大学 | Satellite data reconstruction method and method for detecting response of upper ocean to typhoon |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107945129A (en) * | 2017-11-28 | 2018-04-20 | 电子科技大学 | A kind of MRI image reconstructing method |
CN110827198A (en) * | 2019-10-14 | 2020-02-21 | 唐山学院 | Multi-camera panoramic image construction method based on compressed sensing and super-resolution reconstruction |
CN111047661A (en) * | 2019-12-12 | 2020-04-21 | 重庆大学 | CS-MRI image reconstruction method based on sparse manifold joint constraint |
WO2020157754A1 (en) * | 2019-01-30 | 2020-08-06 | Technology Innovation Momentum Fund (Israel) Limited Partnership | System and method for reconstruction of compressed signal data using artificial neural networking |
-
2020
- 2020-12-04 CN CN202011409107.3A patent/CN112579687B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107945129A (en) * | 2017-11-28 | 2018-04-20 | 电子科技大学 | A kind of MRI image reconstructing method |
WO2020157754A1 (en) * | 2019-01-30 | 2020-08-06 | Technology Innovation Momentum Fund (Israel) Limited Partnership | System and method for reconstruction of compressed signal data using artificial neural networking |
CN110827198A (en) * | 2019-10-14 | 2020-02-21 | 唐山学院 | Multi-camera panoramic image construction method based on compressed sensing and super-resolution reconstruction |
CN111047661A (en) * | 2019-12-12 | 2020-04-21 | 重庆大学 | CS-MRI image reconstruction method based on sparse manifold joint constraint |
Non-Patent Citations (2)
Title |
---|
基于Jacobi ADMM的传感网分布式压缩感知数据重构算法;李国瑞等;《计算机研究与发展》;20200607(第06期);164-171 * |
时变水声信道的动态压缩感知估计;江伟华等;《声学学报》;20190515(第03期);90-98 * |
Also Published As
Publication number | Publication date |
---|---|
CN112579687A (en) | 2021-03-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103810755B (en) | Compressed sensing spectrum picture method for reconstructing based on documents structured Cluster rarefaction representation | |
Zhu et al. | Wind speed behaviors feather analysis and its utilization on wind speed prediction using 3D-CNN | |
CN103020939B (en) | Method for removing large-area thick clouds for optical remote sensing images through multi-temporal data | |
CN110827198B (en) | Multi-camera panoramic image construction method based on compressed sensing and super-resolution reconstruction | |
CN105513026A (en) | Compressed sensing reconstruction method based on image nonlocal similarity | |
Xie et al. | Adaptive measurement network for cs image reconstruction | |
CN111046737B (en) | Efficient intelligent sensing acquisition method for microseism signal detection | |
CN102075749B (en) | Image compression reconstruction method under compressed sensing frame based on non-convex model | |
CN109684314A (en) | A kind of wireless sensor network missing value estimation method based on space structure | |
CN112579687B (en) | Marine environment monitoring data compressed sensing online reconstruction method | |
CN107705342A (en) | A kind of infrared image reconstructing method based on self-adapting generalized orthogonal matching pursuit | |
CN105894547A (en) | Image processing method based on group-wave transformation compressed sensing | |
CN103258134B (en) | A kind of dimension-reduction treatment method of vibration signal of higher-dimension | |
CN104318046A (en) | System and method for incrementally converting high dimensional data into low dimensional data | |
Liu et al. | Detail injection-based spatio-temporal fusion for remote sensing images with land cover changes | |
CN114140357A (en) | Multi-temporal remote sensing image cloud region reconstruction method based on cooperative attention mechanism | |
CN102355752A (en) | Wireless sensor network compressed sensing measurement matrix based on expander graph and reconfiguring method | |
CN114387531A (en) | Ground surface temperature downscaling method based on improved geographic weighted regression model | |
CN112463877A (en) | Compressed sensing observation and reconstruction method of three-dimensional marine environment monitoring data | |
CN103295249B (en) | A kind of method for reconstructing of various visual angles compressed sensing image | |
CN112785684B (en) | Three-dimensional model reconstruction method based on local information weighting mechanism | |
Sadiq et al. | Recovering defective Landsat 7 Enhanced Thematic Mapper Plus images via multiple linear regression model | |
CN103442236B (en) | A kind of multistage remote sensing signal compaction coding method with point Dimension Vector Quantization of Linear Prediction | |
CN112150566A (en) | Dense residual error network image compressed sensing reconstruction method based on feature fusion | |
CN110139233A (en) | Wireless sensor network data restorative procedure based on space-time characteristic fusion |
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 |