In this paper,the application of an algorithm for precipitation retrieval based on Himawari-8 (H8) satellite infrared data is studied.Based on GPM precipitation data and H8 Infrared spectrum channel brightness tempera...In this paper,the application of an algorithm for precipitation retrieval based on Himawari-8 (H8) satellite infrared data is studied.Based on GPM precipitation data and H8 Infrared spectrum channel brightness temperature data,corresponding "precipitation field dictionary" and "channel brightness temperature dictionary" are formed.The retrieval of precipitation field based on brightness temperature data is studied through the classification rule of k-nearest neighbor domain (KNN) and regularization constraint.Firstly,the corresponding "dictionary" is constructed according to the training sample database of the matched GPM precipitation data and H8 brightness temperature data.Secondly,according to the fact that precipitation characteristics in small organizations in different storm environments are often repeated,KNN is used to identify the spectral brightness temperature signal of "precipitation" and "non-precipitation" based on "the dictionary".Finally,the precipitation field retrieval is carried out in the precipitation signal "subspace" based on the regular term constraint method.In the process of retrieval,the contribution rate of brightness temperature retrieval of different channels was determined by Bayesian model averaging (BMA) model.The preliminary experimental results based on the "quantitative" evaluation indexes show that the precipitation of H8 retrieval has a good correlation with the GPM truth value,with a small error and similar structure.展开更多
The classical elastic impedance(EI) inversion method,however,is based on the L2-norm misfit function and considerably sensitive to outliers,assuming the noise of the seismic data to be the Guassian-distribution.So we ...The classical elastic impedance(EI) inversion method,however,is based on the L2-norm misfit function and considerably sensitive to outliers,assuming the noise of the seismic data to be the Guassian-distribution.So we have developed a more robust elastic impedance inversion based on the L1-norm misfit function,and the noise is assumed to be non-Gaussian.Meanwhile,some regularization methods including the sparse constraint regularization and elastic impedance point constraint regularization are incorporated to improve the ill-posed characteristics of the seismic inversion problem.Firstly,we create the Ll-norm misfit objective function of pre-stack inversion problem based on the Bayesian scheme within the sparse constraint regularization and elastic impedance point constraint regularization.And then,we obtain more robust elastic impedances of different angles which are less sensitive to outliers in seismic data by using the IRLS strategy.Finally,we extract the P-wave and S-wave velocity and density by using the more stable parameter extraction method.Tests on synthetic data show that the P-wave and S-wave velocity and density parameters are still estimated reasonable with moderate noise.A test on the real data set shows that compared to the results of the classical elastic impedance inversion method,the estimated results using the proposed method can get better lateral continuity and more distinct show of the gas,verifying the feasibility and stability of the method.展开更多
Some remarks are made on the use of the Abadie constraint qualification, the Guignard constraint qualifications and the Guignard regularity condition in obtaining weak and strong Kuhn-Tucker type optimality conditions...Some remarks are made on the use of the Abadie constraint qualification, the Guignard constraint qualifications and the Guignard regularity condition in obtaining weak and strong Kuhn-Tucker type optimality conditions in differentiable vector optimization problems.展开更多
Surface mass anomalies estimated by mass concentration(mascon)approach using Gravity Recovery and Climate Experiment(GRACE)observations with regularization constraints generally present higher spatial resolution than ...Surface mass anomalies estimated by mass concentration(mascon)approach using Gravity Recovery and Climate Experiment(GRACE)observations with regularization constraints generally present higher spatial resolution than the spheric harmonic(SH)solutions.To analyze the influence of different types of constraints on the estimation of mascon solutions,we carried out a closed-loop simulation experiment to estimate surface mass anomalies over South America based on simulated GRACE intersatellite geopotential differences.Tikhonov regularization with spatial constraint(SC),uniform weighting constraint(UWC),and a prior information constraint(APC)were employed to stabilize the mascon solutions,and the corresponding optimal regularization parameters were determined based on the minimum residual root-mean-square(RMS)criterion.The results show that mascon solutions estimated under different types of constraints are consistent and equivalent when the optimal regularization parameters are selected.The spatial distributions and main characteristics of regional surface mass anomalies estimated by the three types of constraints agree well,and the values of residual RMS with different constraints are very close.But due to the smoothing effect of regularization,the signal strength of mascon solutions is a bit weaker than that of original true signal,especially in the regions with strong signals.In addition,due to the ill-conditioned problem is more serious for higher grid resolution,the relative contribution of the three types of constraints to the final mascon solutions would be stronger.The results show that the averages of relative contribution percentages of these constraints for 2°×2° mascon grids are 80%-90%,while the corresponding values for 4°×4° mascon grids are 30%-60%.However,based on the minimum residual RMS criterion,the accuracy of estimation results is not affected by the type of constraints and their relative contribution to the final mascon solutions.展开更多
A new normalized least mean square(NLMS) adaptive filter is first derived from a cost function, which incorporates the conventional one of the NLMS with a minimum-disturbance(MD)constraint. A variable regularization f...A new normalized least mean square(NLMS) adaptive filter is first derived from a cost function, which incorporates the conventional one of the NLMS with a minimum-disturbance(MD)constraint. A variable regularization factor(RF) is then employed to control the contribution made by the MD constraint in the cost function. Analysis results show that the RF can be taken as a combination of the step size and regularization parameter in the conventional NLMS. This implies that these parameters can be jointly controlled by simply tuning the RF as the proposed algorithm does. It also demonstrates that the RF can accelerate the convergence rate of the proposed algorithm and its optimal value can be obtained by minimizing the squared noise-free posteriori error. A method for automatically determining the value of the RF is also presented, which is free of any prior knowledge of the noise. While simulation results verify the analytical ones, it is also illustrated that the performance of the proposed algorithm is superior to the state-of-art ones in both the steady-state misalignment and the convergence rate. A novel algorithm is proposed to solve some problems. Simulation results show the effectiveness of the proposed algorithm.展开更多
Image restoration is often solved by minimizing an energy function consisting of a data-fidelity term and a regularization term. A regularized convex term can usually preserve the image edges well in the restored imag...Image restoration is often solved by minimizing an energy function consisting of a data-fidelity term and a regularization term. A regularized convex term can usually preserve the image edges well in the restored image. In this paper, we consider a class of convex and edge-preserving regularization functions, I.e., multiplicative half-quadratic regularizations, and we use the Newton method to solve the correspondingly reduced systems of nonlinear equations. At each Newton iterate, the preconditioned conjugate gradient method, incorporated with a constraint preconditioner, is employed to solve the structured Newton equation that has a symmetric positive definite coefficient matrix.The igenvalue bounds of the preconditioned matrix are deliberately derived, which can be used to estimate the convergence speed of the preconditioned conjugate gradient method. We use experimental results to demonstrate that this new approach is efficient,and the effect of image restoration is r0easonably well.展开更多
高光谱图像(HSI)具有良好的光谱识别能力,但在采集过程中易受到混合噪声的污染,严重影响了后续任务的精度,因此HSI去噪是重要的预处理过程。针对现有去噪方法对空间-光谱先验信息利用不足、条纹噪声建模不合理的问题,提出一种新的基于...高光谱图像(HSI)具有良好的光谱识别能力,但在采集过程中易受到混合噪声的污染,严重影响了后续任务的精度,因此HSI去噪是重要的预处理过程。针对现有去噪方法对空间-光谱先验信息利用不足、条纹噪声建模不合理的问题,提出一种新的基于群稀疏正则化的高光谱图像去噪算法。该算法将干净HSI的空间-光谱低秩特性和各波段上条纹噪声的低秩结构融入一个新框架,实现了干净HSI与高强度结构化条纹噪声的分离;同时为了有效保持图像的边缘信息,在去噪模型中引入新的群稀疏正则化,即基于L_(2,1)范数的增强型三维全变分正则化(enhanced 3D total variation, E3DTV),充分挖掘HSI差分图像的稀疏先验信息,进一步提升了图像的分段平滑性。采用交替方向乘子法对变量优化求解,在仿真和真实数据集上进行数值实验表明,所提模型具有更好的去噪和去条纹性能,在视觉效果和定量评价结果上都明显优于其他对比算法。展开更多
基金Supported by National Natural Science Foundation of China(41805080)Natural Science Foundation of Anhui Province,China(1708085QD89)+1 种基金Key Research and Development Program Projects of Anhui Province,China(201904a07020099)Open Foundation Project Shenyang Institute of Atmospheric Environment,China Meteorological Administration(2016SYIAE14)
文摘In this paper,the application of an algorithm for precipitation retrieval based on Himawari-8 (H8) satellite infrared data is studied.Based on GPM precipitation data and H8 Infrared spectrum channel brightness temperature data,corresponding "precipitation field dictionary" and "channel brightness temperature dictionary" are formed.The retrieval of precipitation field based on brightness temperature data is studied through the classification rule of k-nearest neighbor domain (KNN) and regularization constraint.Firstly,the corresponding "dictionary" is constructed according to the training sample database of the matched GPM precipitation data and H8 brightness temperature data.Secondly,according to the fact that precipitation characteristics in small organizations in different storm environments are often repeated,KNN is used to identify the spectral brightness temperature signal of "precipitation" and "non-precipitation" based on "the dictionary".Finally,the precipitation field retrieval is carried out in the precipitation signal "subspace" based on the regular term constraint method.In the process of retrieval,the contribution rate of brightness temperature retrieval of different channels was determined by Bayesian model averaging (BMA) model.The preliminary experimental results based on the "quantitative" evaluation indexes show that the precipitation of H8 retrieval has a good correlation with the GPM truth value,with a small error and similar structure.
基金Projects(U1562215,41674130,41404088)supported by the National Natural Science Foundation of ChinaProjects(2013CB228604,2014CB239201)supported by the National Basic Research Program of China+1 种基金Projects(2016ZX05027004-001,2016ZX05002006-009)supported by the National Oil and Gas Major Projects of ChinaProject(15CX08002A)supported by the Fundamental Research Funds for the Central Universities,China
文摘The classical elastic impedance(EI) inversion method,however,is based on the L2-norm misfit function and considerably sensitive to outliers,assuming the noise of the seismic data to be the Guassian-distribution.So we have developed a more robust elastic impedance inversion based on the L1-norm misfit function,and the noise is assumed to be non-Gaussian.Meanwhile,some regularization methods including the sparse constraint regularization and elastic impedance point constraint regularization are incorporated to improve the ill-posed characteristics of the seismic inversion problem.Firstly,we create the Ll-norm misfit objective function of pre-stack inversion problem based on the Bayesian scheme within the sparse constraint regularization and elastic impedance point constraint regularization.And then,we obtain more robust elastic impedances of different angles which are less sensitive to outliers in seismic data by using the IRLS strategy.Finally,we extract the P-wave and S-wave velocity and density by using the more stable parameter extraction method.Tests on synthetic data show that the P-wave and S-wave velocity and density parameters are still estimated reasonable with moderate noise.A test on the real data set shows that compared to the results of the classical elastic impedance inversion method,the estimated results using the proposed method can get better lateral continuity and more distinct show of the gas,verifying the feasibility and stability of the method.
文摘Some remarks are made on the use of the Abadie constraint qualification, the Guignard constraint qualifications and the Guignard regularity condition in obtaining weak and strong Kuhn-Tucker type optimality conditions in differentiable vector optimization problems.
基金funded by the National Key Research and Development Program of China(Grant No.2018YFC1503503)the National Natural Science Foundation of China(Grant Nos.41974015,42061134007,41474019)。
文摘Surface mass anomalies estimated by mass concentration(mascon)approach using Gravity Recovery and Climate Experiment(GRACE)observations with regularization constraints generally present higher spatial resolution than the spheric harmonic(SH)solutions.To analyze the influence of different types of constraints on the estimation of mascon solutions,we carried out a closed-loop simulation experiment to estimate surface mass anomalies over South America based on simulated GRACE intersatellite geopotential differences.Tikhonov regularization with spatial constraint(SC),uniform weighting constraint(UWC),and a prior information constraint(APC)were employed to stabilize the mascon solutions,and the corresponding optimal regularization parameters were determined based on the minimum residual root-mean-square(RMS)criterion.The results show that mascon solutions estimated under different types of constraints are consistent and equivalent when the optimal regularization parameters are selected.The spatial distributions and main characteristics of regional surface mass anomalies estimated by the three types of constraints agree well,and the values of residual RMS with different constraints are very close.But due to the smoothing effect of regularization,the signal strength of mascon solutions is a bit weaker than that of original true signal,especially in the regions with strong signals.In addition,due to the ill-conditioned problem is more serious for higher grid resolution,the relative contribution of the three types of constraints to the final mascon solutions would be stronger.The results show that the averages of relative contribution percentages of these constraints for 2°×2° mascon grids are 80%-90%,while the corresponding values for 4°×4° mascon grids are 30%-60%.However,based on the minimum residual RMS criterion,the accuracy of estimation results is not affected by the type of constraints and their relative contribution to the final mascon solutions.
基金supported by the National Natural Science Foundation of China(61571131 11604055)
文摘A new normalized least mean square(NLMS) adaptive filter is first derived from a cost function, which incorporates the conventional one of the NLMS with a minimum-disturbance(MD)constraint. A variable regularization factor(RF) is then employed to control the contribution made by the MD constraint in the cost function. Analysis results show that the RF can be taken as a combination of the step size and regularization parameter in the conventional NLMS. This implies that these parameters can be jointly controlled by simply tuning the RF as the proposed algorithm does. It also demonstrates that the RF can accelerate the convergence rate of the proposed algorithm and its optimal value can be obtained by minimizing the squared noise-free posteriori error. A method for automatically determining the value of the RF is also presented, which is free of any prior knowledge of the noise. While simulation results verify the analytical ones, it is also illustrated that the performance of the proposed algorithm is superior to the state-of-art ones in both the steady-state misalignment and the convergence rate. A novel algorithm is proposed to solve some problems. Simulation results show the effectiveness of the proposed algorithm.
基金supported by the National Basic Research Program (No.2005CB321702)the National Outstanding Young Scientist Foundation(No. 10525102)the Specialized Research Grant for High Educational Doctoral Program(Nos. 20090211120011 and LZULL200909),Hong Kong RGC grants and HKBU FRGs
文摘Image restoration is often solved by minimizing an energy function consisting of a data-fidelity term and a regularization term. A regularized convex term can usually preserve the image edges well in the restored image. In this paper, we consider a class of convex and edge-preserving regularization functions, I.e., multiplicative half-quadratic regularizations, and we use the Newton method to solve the correspondingly reduced systems of nonlinear equations. At each Newton iterate, the preconditioned conjugate gradient method, incorporated with a constraint preconditioner, is employed to solve the structured Newton equation that has a symmetric positive definite coefficient matrix.The igenvalue bounds of the preconditioned matrix are deliberately derived, which can be used to estimate the convergence speed of the preconditioned conjugate gradient method. We use experimental results to demonstrate that this new approach is efficient,and the effect of image restoration is r0easonably well.
文摘高光谱图像(HSI)具有良好的光谱识别能力,但在采集过程中易受到混合噪声的污染,严重影响了后续任务的精度,因此HSI去噪是重要的预处理过程。针对现有去噪方法对空间-光谱先验信息利用不足、条纹噪声建模不合理的问题,提出一种新的基于群稀疏正则化的高光谱图像去噪算法。该算法将干净HSI的空间-光谱低秩特性和各波段上条纹噪声的低秩结构融入一个新框架,实现了干净HSI与高强度结构化条纹噪声的分离;同时为了有效保持图像的边缘信息,在去噪模型中引入新的群稀疏正则化,即基于L_(2,1)范数的增强型三维全变分正则化(enhanced 3D total variation, E3DTV),充分挖掘HSI差分图像的稀疏先验信息,进一步提升了图像的分段平滑性。采用交替方向乘子法对变量优化求解,在仿真和真实数据集上进行数值实验表明,所提模型具有更好的去噪和去条纹性能,在视觉效果和定量评价结果上都明显优于其他对比算法。