In this paper, we deal with (nonlinear) ill-posed problems that are regularized by minimizing Tikhonov-type functionals. If the minimization is tedious for some penalty term begin{document}$ P_0 $end{document}, we approximate it by a family of penalty terms begin{document}$ ({P_beta}) $end{document} having nicer properties and analyze what happens as begin{document}$ betato 0 $end{document}.
We investigate the discrepancy principle for the choice of the regularization parameter and apply all results to linear problems with sparsity constraints. Numerical results show that the proposed method yields good results.
In this paper, we deal with (nonlinear) ill-posed problems that are regularized by minimizing Tikhonov-type functionals. If the minimization is tedious for some penalty term begin{document}$ P_0 $end{document}, we approximate it by a family of penalty terms begin{document}$ ({P_beta}) $end{document} having nicer properties and analyze what happens as begin{document}$ betato 0 $end{document}.We investigate the discrepancy principle for the choice of the regularization parameter and apply all results to linear problems with sparsity constraints. Numerical results show that the proposed method yields good results.
{"title":"On Tikhonov-type regularization with approximated penalty terms","authors":"A. Neubauer","doi":"10.3934/IPI.2021027","DOIUrl":"https://doi.org/10.3934/IPI.2021027","url":null,"abstract":"<p style='text-indent:20px;'>In this paper, we deal with (nonlinear) ill-posed problems that are regularized by minimizing Tikhonov-type functionals. If the minimization is tedious for some penalty term <inline-formula><tex-math id=\"M1\">begin{document}$ P_0 $end{document}</tex-math></inline-formula>, we approximate it by a family of penalty terms <inline-formula><tex-math id=\"M2\">begin{document}$ ({P_beta}) $end{document}</tex-math></inline-formula> having nicer properties and analyze what happens as <inline-formula><tex-math id=\"M3\">begin{document}$ betato 0 $end{document}</tex-math></inline-formula>.</p><p style='text-indent:20px;'>We investigate the discrepancy principle for the choice of the regularization parameter and apply all results to linear problems with sparsity constraints. Numerical results show that the proposed method yields good results.</p>","PeriodicalId":50274,"journal":{"name":"Inverse Problems and Imaging","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2021-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"76107262","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Lingfeng Li, Shousheng Luo, X. Tai, Jiang Yang, kaifeng China Algorithms
Seeking the convex hull of an object (or point set) is a very fundamental problem arising from various tasks. In this work, we propose a variational approach based on the level-set representation for convex hulls of 2-dimensional objects. This method can adapt to exact and inexact convex hull problems. In addition, this method can compute multiple convex hulls simultaneously. In this model, the convex hull is characterized by the zero sublevel-set of a level-set function. For the exact case, we require the zero sublevel-set to be convex and contain the whole given object, where the convexity is characterized by the non-negativity of Laplacian of the level-set function. Then, the convex hull can be obtained by minimizing the area of the zero sublevel-set. For the inexact case, instead of requiring all the given points are included, we penalize the distance from all given points to the zero sublevel-set. Especially, the inexact model can handle the convex hull problem of the given set with outliers very well, while most of the existing methods fail. An efficient numerical scheme using the alternating direction method of multipliers is developed. Numerical examples are given to demonstrate the advantages of the proposed methods.
{"title":"A new variational approach based on level-set function for convex hull problem with outliers","authors":"Lingfeng Li, Shousheng Luo, X. Tai, Jiang Yang, kaifeng China Algorithms","doi":"10.3934/ipi.2020070","DOIUrl":"https://doi.org/10.3934/ipi.2020070","url":null,"abstract":"Seeking the convex hull of an object (or point set) is a very fundamental problem arising from various tasks. In this work, we propose a variational approach based on the level-set representation for convex hulls of 2-dimensional objects. This method can adapt to exact and inexact convex hull problems. In addition, this method can compute multiple convex hulls simultaneously. In this model, the convex hull is characterized by the zero sublevel-set of a level-set function. For the exact case, we require the zero sublevel-set to be convex and contain the whole given object, where the convexity is characterized by the non-negativity of Laplacian of the level-set function. Then, the convex hull can be obtained by minimizing the area of the zero sublevel-set. For the inexact case, instead of requiring all the given points are included, we penalize the distance from all given points to the zero sublevel-set. Especially, the inexact model can handle the convex hull problem of the given set with outliers very well, while most of the existing methods fail. An efficient numerical scheme using the alternating direction method of multipliers is developed. Numerical examples are given to demonstrate the advantages of the proposed methods.","PeriodicalId":50274,"journal":{"name":"Inverse Problems and Imaging","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2021-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"89906835","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Bao Wang, A. Lin, Penghang Yin, Wei Zhu, A. Bertozzi, S. Osher
We improve the robustness of Deep Neural Net (DNN) to adversarial attacks by using an interpolating function as the output activation. This data-dependent activation remarkably improves both the generalization and robustness of DNN. In the CIFAR10 benchmark, we raise the robust accuracy of the adversarially trained ResNet20 from begin{document}$ sim 46% $end{document} to begin{document}$ sim 69% $end{document} under the state-of-the-art Iterative Fast Gradient Sign Method (IFGSM) based adversarial attack. When we combine this data-dependent activation with total variation minimization on adversarial images and training data augmentation, we achieve an improvement in robust accuracy by 38.9 begin{document}$ % $end{document} for ResNet56 under the strongest IFGSM attack. Furthermore, We provide an intuitive explanation of our defense by analyzing the geometry of the feature space.
We improve the robustness of Deep Neural Net (DNN) to adversarial attacks by using an interpolating function as the output activation. This data-dependent activation remarkably improves both the generalization and robustness of DNN. In the CIFAR10 benchmark, we raise the robust accuracy of the adversarially trained ResNet20 from begin{document}$ sim 46% $end{document} to begin{document}$ sim 69% $end{document} under the state-of-the-art Iterative Fast Gradient Sign Method (IFGSM) based adversarial attack. When we combine this data-dependent activation with total variation minimization on adversarial images and training data augmentation, we achieve an improvement in robust accuracy by 38.9 begin{document}$ % $end{document} for ResNet56 under the strongest IFGSM attack. Furthermore, We provide an intuitive explanation of our defense by analyzing the geometry of the feature space.
{"title":"Adversarial defense via the data-dependent activation, total variation minimization, and adversarial training","authors":"Bao Wang, A. Lin, Penghang Yin, Wei Zhu, A. Bertozzi, S. Osher","doi":"10.3934/ipi.2020046","DOIUrl":"https://doi.org/10.3934/ipi.2020046","url":null,"abstract":"We improve the robustness of Deep Neural Net (DNN) to adversarial attacks by using an interpolating function as the output activation. This data-dependent activation remarkably improves both the generalization and robustness of DNN. In the CIFAR10 benchmark, we raise the robust accuracy of the adversarially trained ResNet20 from begin{document}$ sim 46% $end{document} to begin{document}$ sim 69% $end{document} under the state-of-the-art Iterative Fast Gradient Sign Method (IFGSM) based adversarial attack. When we combine this data-dependent activation with total variation minimization on adversarial images and training data augmentation, we achieve an improvement in robust accuracy by 38.9 begin{document}$ % $end{document} for ResNet56 under the strongest IFGSM attack. Furthermore, We provide an intuitive explanation of our defense by analyzing the geometry of the feature space.","PeriodicalId":50274,"journal":{"name":"Inverse Problems and Imaging","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2021-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"86454372","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
This paper considers the inverse elastic wave scattering by a bounded penetrable or impenetrable scatterer. We propose a novel technique to show that the elastic obstacle can be uniquely determined by its far-field pattern associated with all incident plane waves at a fixed frequency. In the first part of this paper, we establish the mixed reciprocity relation between the far-field pattern corresponding to special point sources and the scattered field corresponding to plane waves, and the mixed reciprocity relation is the key point to show the uniqueness results. In the second part, besides the mixed reciprocity relation, a priori estimates of solution to the transmission problem with boundary data in begin{document}$ [L^{mathrm{q}}(partialOmega)]^{3} $end{document} ( begin{document}$ 1 ) is deeply investigated by the integral equation method and also we have constructed a well-posed modified static interior transmission problem on a small domain to obtain the uniqueness result.
{"title":"The uniqueness of the inverse elastic wave scattering problem based on the mixed reciprocity relation","authors":"Jianlin Xiang, G. Yan","doi":"10.3934/ipi.2021004","DOIUrl":"https://doi.org/10.3934/ipi.2021004","url":null,"abstract":"This paper considers the inverse elastic wave scattering by a bounded penetrable or impenetrable scatterer. We propose a novel technique to show that the elastic obstacle can be uniquely determined by its far-field pattern associated with all incident plane waves at a fixed frequency. In the first part of this paper, we establish the mixed reciprocity relation between the far-field pattern corresponding to special point sources and the scattered field corresponding to plane waves, and the mixed reciprocity relation is the key point to show the uniqueness results. In the second part, besides the mixed reciprocity relation, a priori estimates of solution to the transmission problem with boundary data in begin{document}$ [L^{mathrm{q}}(partialOmega)]^{3} $end{document} ( begin{document}$ 1 ) is deeply investigated by the integral equation method and also we have constructed a well-posed modified static interior transmission problem on a small domain to obtain the uniqueness result.","PeriodicalId":50274,"journal":{"name":"Inverse Problems and Imaging","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2021-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"91073438","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
In this paper, a Cauchy problem of non-homogenous stochastic heat equation is considered together with its inverse source problem, where the source term is assumed to be driven by an additive white noise. The Cauchy problem (direct problem) is to determine the displacement of random temperature field, while the considered inverse problem is to reconstruct the statistical properties of the random source, i.e. the mean and variance of the random source. It is proved constructively that the Cauchy problem has a unique mild solution, which is expressed an integral form. Then the inverse random source problem is formulated into two Fredholm integral equations of the first kind, which are typically ill-posed. To obtain stable inverse solutions, the regularized block Kaczmarz method is introduced to solve the two Fredholm integral equations. Finally, numerical experiments are given to show that the proposed method is efficient and robust for reconstructing the statistical properties of the random source.
{"title":"Cauchy problem of non-homogenous stochastic heat equation and application to inverse random source problem","authors":"Shuli Chen, Zewen Wang, Guo-Xin Chen","doi":"10.3934/IPI.2021008","DOIUrl":"https://doi.org/10.3934/IPI.2021008","url":null,"abstract":"In this paper, a Cauchy problem of non-homogenous stochastic heat equation is considered together with its inverse source problem, where the source term is assumed to be driven by an additive white noise. The Cauchy problem (direct problem) is to determine the displacement of random temperature field, while the considered inverse problem is to reconstruct the statistical properties of the random source, i.e. the mean and variance of the random source. It is proved constructively that the Cauchy problem has a unique mild solution, which is expressed an integral form. Then the inverse random source problem is formulated into two Fredholm integral equations of the first kind, which are typically ill-posed. To obtain stable inverse solutions, the regularized block Kaczmarz method is introduced to solve the two Fredholm integral equations. Finally, numerical experiments are given to show that the proposed method is efficient and robust for reconstructing the statistical properties of the random source.","PeriodicalId":50274,"journal":{"name":"Inverse Problems and Imaging","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2021-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"76341180","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Variational source conditions for inverse Robin and flux problems by partial measurements","authors":"De-Han Chen, Daijun Jiang, I. Yousept, J. Zou","doi":"10.3934/ipi.2021050","DOIUrl":"https://doi.org/10.3934/ipi.2021050","url":null,"abstract":"<jats:p xml:lang=\"fr\" />","PeriodicalId":50274,"journal":{"name":"Inverse Problems and Imaging","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2021-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"81321815","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
In this paper we study the inverse problem of determining an electrical inclusion in a multi-layer composite from boundary measurements in 2D. We assume the conductivities in different layers are different and derive a stability estimate for the linearized map with explicit formulae on the conductivity and the thickness of each layer. Intuitively, if an inclusion is surrounded by a highly conductive layer, then, in view of "the principle of the least work", the current will take a path in the highly conductive layer and disregard the existence of the inclusion. Consequently, a worse stability of identifying the hidden inclusion is expected in this case. Our estimates indeed show that the ill-posedness of the problem increases as long as the conductivity of some layer becomes large. This work is an extension of the previous result by Nagayasu-Uhlmann-Wang[15], where a depth-dependent estimate is derived when an inclusion is deeply hidden in a conductor. Estimates in this work also show the influence of the depth of the inclusion.
{"title":"Refined stability estimates in electrical impedance tomography with multi-layer structure","authors":"Haigang Li, Jenn-Nan Wang, Ling Wang","doi":"10.3934/ipi.2021048","DOIUrl":"https://doi.org/10.3934/ipi.2021048","url":null,"abstract":"In this paper we study the inverse problem of determining an electrical inclusion in a multi-layer composite from boundary measurements in 2D. We assume the conductivities in different layers are different and derive a stability estimate for the linearized map with explicit formulae on the conductivity and the thickness of each layer. Intuitively, if an inclusion is surrounded by a highly conductive layer, then, in view of \"the principle of the least work\", the current will take a path in the highly conductive layer and disregard the existence of the inclusion. Consequently, a worse stability of identifying the hidden inclusion is expected in this case. Our estimates indeed show that the ill-posedness of the problem increases as long as the conductivity of some layer becomes large. This work is an extension of the previous result by Nagayasu-Uhlmann-Wang[15], where a depth-dependent estimate is derived when an inclusion is deeply hidden in a conductor. Estimates in this work also show the influence of the depth of the inclusion.","PeriodicalId":50274,"journal":{"name":"Inverse Problems and Imaging","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2021-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"91086780","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Ying Zhang, Xuhua Ren, B. Clifford, Qian Wang, Xiaoqun Zhang
In recent years multi-modal data processing methods have gained considerable research interest as technological advancements in imaging, computing, and data storage have made the collection of redundant, multi-modal data more commonplace. In this work we present an image restoration method tailored for scenarios where pre-existing, high-quality images from different modalities or contrasts are available in addition to the target image. Our method is based on a novel network architecture which combines the benefits of traditional multi-scale signal representation, such as wavelets, with more recent concepts from data fusion methods. Results from numerical simulations in which T1-weighted MRI images are used to restore noisy and undersampled T2-weighted images demonstrate that the proposed network successfully utilizes information from high-quality reference images to improve the restoration quality of the target image beyond that of existing popular methods.
{"title":"Image fusion network for dual-modal restoration","authors":"Ying Zhang, Xuhua Ren, B. Clifford, Qian Wang, Xiaoqun Zhang","doi":"10.3934/ipi.2021067","DOIUrl":"https://doi.org/10.3934/ipi.2021067","url":null,"abstract":"In recent years multi-modal data processing methods have gained considerable research interest as technological advancements in imaging, computing, and data storage have made the collection of redundant, multi-modal data more commonplace. In this work we present an image restoration method tailored for scenarios where pre-existing, high-quality images from different modalities or contrasts are available in addition to the target image. Our method is based on a novel network architecture which combines the benefits of traditional multi-scale signal representation, such as wavelets, with more recent concepts from data fusion methods. Results from numerical simulations in which T1-weighted MRI images are used to restore noisy and undersampled T2-weighted images demonstrate that the proposed network successfully utilizes information from high-quality reference images to improve the restoration quality of the target image beyond that of existing popular methods.","PeriodicalId":50274,"journal":{"name":"Inverse Problems and Imaging","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2021-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"91367630","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
A spectral problem occurring in description of small transverse vibrations of a star graph of Stieltjes strings is considered. At all but one pendant vertices Dirichlet conditions are imposed which mean that these vertices are clamped. One vertex (the root) can move with damping in the direction orthogonal to the equilibrium position of the strings. We describe the spectrum of such spectral problem. The corresponding inverse problem lies in recovering the values of point masses and the lengths of the intervals between the masses using the spectrum and some other parameters. We propose conditions on a sequence of complex numbers and a collection of real numbers to be the spectrum of a problem we consider and the lengths of the edges, correspondingly.
{"title":"Direct and inverse spectral problems for a star graph of Stieltjes strings damped at a pendant vertex","authors":"Lu Yang, Guangsheng Wei, V. Pivovarchik","doi":"10.3934/ipi.2020063","DOIUrl":"https://doi.org/10.3934/ipi.2020063","url":null,"abstract":"A spectral problem occurring in description of small transverse vibrations of a star graph of Stieltjes strings is considered. At all but one pendant vertices Dirichlet conditions are imposed which mean that these vertices are clamped. One vertex (the root) can move with damping in the direction orthogonal to the equilibrium position of the strings. We describe the spectrum of such spectral problem. The corresponding inverse problem lies in recovering the values of point masses and the lengths of the intervals between the masses using the spectrum and some other parameters. We propose conditions on a sequence of complex numbers and a collection of real numbers to be the spectrum of a problem we consider and the lengths of the edges, correspondingly.","PeriodicalId":50274,"journal":{"name":"Inverse Problems and Imaging","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2021-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"76969476","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
This short note was motivated by our efforts to investigate whether there exists a half plane free of transmission eigenvalues for Maxwell's equations. This question is related to solvability of the time domain interior transmission problem which plays a fundamental role in the justification of linear sampling and factorization methods with time dependent data. Our original goal was to adapt semiclassical analysis techniques developed in [21,23] to prove that for some combination of electromagnetic parameters, the transmission eigenvalues lie in a strip around the real axis. Unfortunately we failed. To try to understand why, we looked at the particular example of spherically symmetric media, which provided us with some insight on why we couldn't prove the above result. Hence this paper reports our findings on the location of all transmission eigenvalues and the existence of complex transmission eigenvalues for Maxwell's equations for spherically stratified media. We hope that these results can provide reasonable conjectures for general electromagnetic media.
{"title":"A note on transmission eigenvalues in electromagnetic scattering theory","authors":"F. Cakoni, S. Meng, Jingni Xiao","doi":"10.3934/IPI.2021025","DOIUrl":"https://doi.org/10.3934/IPI.2021025","url":null,"abstract":"This short note was motivated by our efforts to investigate whether there exists a half plane free of transmission eigenvalues for Maxwell's equations. This question is related to solvability of the time domain interior transmission problem which plays a fundamental role in the justification of linear sampling and factorization methods with time dependent data. Our original goal was to adapt semiclassical analysis techniques developed in [21,23] to prove that for some combination of electromagnetic parameters, the transmission eigenvalues lie in a strip around the real axis. Unfortunately we failed. To try to understand why, we looked at the particular example of spherically symmetric media, which provided us with some insight on why we couldn't prove the above result. Hence this paper reports our findings on the location of all transmission eigenvalues and the existence of complex transmission eigenvalues for Maxwell's equations for spherically stratified media. We hope that these results can provide reasonable conjectures for general electromagnetic media.","PeriodicalId":50274,"journal":{"name":"Inverse Problems and Imaging","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2021-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"74860289","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}