{"title":"An Operator Theory for Analyzing the Resolution of Multi-illumination Imaging Modalities","authors":"Ping Liu, Habib Ammari","doi":"10.1137/23m1551730","DOIUrl":"https://doi.org/10.1137/23m1551730","url":null,"abstract":"SIAM Journal on Imaging Sciences, Volume 16, Issue 4, Page 2105-2143, December 2023. <br/> Abstract. By introducing a new operator theory, we provide a unified mathematical theory for general source resolution in the multi-illumination imaging problem. Our main idea is to transform multi-illumination imaging into single-snapshot imaging with a new imaging kernel that depends on both the illumination patterns and the point spread function of the imaging system. We therefore prove that the resolution of multi-illumination imaging is approximately determined by the essential cutoff frequency of the new imaging kernel, which is roughly limited by the sum of the cutoff frequency of the point spread function and the maximum essential frequency in the illumination patterns. Our theory provides a unified way to estimate the resolution of various existing super-resolution modalities and results in the same estimates as those obtained in experiments. In addition, based on the reformulation of the multi-illumination imaging problem, we also estimate the resolution limits for resolving both complex and positive sources by sparsity-based approaches. We show that the resolution of multi-illumination imaging is approximately determined by the new imaging kernel from our operator theory and better resolution can be realized by sparsity-promoting techniques in practice but only for resolving very sparse sources. This explains experimentally observed phenomena in some sparsity-based super-resolution modalities.","PeriodicalId":49528,"journal":{"name":"SIAM Journal on Imaging Sciences","volume":"68 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2023-11-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138529000","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Marcelo Pereyra, Luis A. Vargas-Mieles, Konstantinos C. Zygalakis
{"title":"The Split Gibbs Sampler Revisited: Improvements to Its Algorithmic Structure and Augmented Target Distribution","authors":"Marcelo Pereyra, Luis A. Vargas-Mieles, Konstantinos C. Zygalakis","doi":"10.1137/22m1506122","DOIUrl":"https://doi.org/10.1137/22m1506122","url":null,"abstract":"Developing efficient Bayesian computation algorithms for imaging inverse problems is challenging due to the dimensionality involved and because Bayesian imaging models are often not smooth. Current state-of-the-art methods often address these difficulties by replacing the posterior density with a smooth approximation that is amenable to efficient exploration by using Langevin Markov chain Monte Carlo (MCMC) methods. An alternative approach is based on data augmentation and relaxation, where auxiliary variables are introduced in order to construct an approximate augmented posterior distribution that is amenable to efficient exploration by Gibbs sampling. This paper proposes a new accelerated proximal MCMC method called latent space SK-ROCK (ls SK-ROCK), which tightly combines the benefits of the two aforementioned strategies. Additionally, instead of viewing the augmented posterior distribution as an approximation of the original model, we propose to consider it as a generalisation of this model. Following on from this, we empirically show that there is a range of values for the relaxation parameter for which the accuracy of the model improves, and propose a stochastic optimisation algorithm to automatically identify the optimal amount of relaxation for a given problem. In this regime, ls SK-ROCK converges faster than competing approaches from the state of the art, and also achieves better accuracy since the underlying augmented Bayesian model has a higher Bayesian evidence. The proposed methodology is demonstrated with a range of numerical experiments related to image deblurring and inpainting, as well as with comparisons with alternative approaches from the state of the art. An open-source implementation of the proposed MCMC methods is available from https://github.com/luisvargasmieles/ls-MCMC.","PeriodicalId":49528,"journal":{"name":"SIAM Journal on Imaging Sciences","volume":" 409","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135185986","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Arttu Arjas, Mikko J. Sillanpää, Andreas S. Hauptmann
{"title":"Sequential Model Correction for Nonlinear Inverse Problems","authors":"Arttu Arjas, Mikko J. Sillanpää, Andreas S. Hauptmann","doi":"10.1137/23m1549286","DOIUrl":"https://doi.org/10.1137/23m1549286","url":null,"abstract":"","PeriodicalId":49528,"journal":{"name":"SIAM Journal on Imaging Sciences","volume":"193 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-10-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135731510","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A Common Lines Approach for Ab Initio Modeling of Molecules with Tetrahedral and Octahedral Symmetry","authors":"Adi Shasha Geva, Yoel Shkolnisky","doi":"10.1137/22m150383x","DOIUrl":"https://doi.org/10.1137/22m150383x","url":null,"abstract":"A main task in cryo-electron microscopy single particle reconstruction is to find a three-dimensional model of a molecule given a set of its randomly oriented and positioned noisy projection-images. In this work, we propose an algorithm for ab-initio reconstruction for molecules with tetrahedral or octahedral symmetry. The algorithm exploits the multiple common lines between each pair of projection-images as well as self common lines within each image. It is robust to noise in the input images as it integrates the information from all images at once. The applicability of the proposed algorithm is demonstrated using experimental cryo-electron microscopy data.","PeriodicalId":49528,"journal":{"name":"SIAM Journal on Imaging Sciences","volume":"44 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-10-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135884303","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"(boldsymbol{L_1-beta L_q}) Minimization for Signal and Image Recovery","authors":"Limei Huo, Wengu Chen, Huanmin Ge, Michael K. Ng","doi":"10.1137/22m1525363","DOIUrl":"https://doi.org/10.1137/22m1525363","url":null,"abstract":"The nonconvex optimization method has attracted increasing attention due to its excellent ability of promoting sparsity in signal processing, image restoration, and machine learning. In this paper, we consider a new minimization method and its applications in signal recovery and image reconstruction because minimization provides an effective way to solve the -ratio sparsity minimization model. Our main contributions are to establish a convex hull decomposition for and investigate RIP-based conditions for stable signal recovery and image reconstruction by minimization. For one-dimensional signal recovery, our derived RIP condition extends existing results. For two-dimensional image recovery under minimization of image gradients, we provide the error estimate of the resulting optimal solutions in terms of sparsity and noise level, which is missing in the literature. Numerical results of the limited angle problem in computed tomography imaging and image deblurring are presented to validate the efficiency and superiority of the proposed minimization method among the state-of-art image recovery methods.","PeriodicalId":49528,"journal":{"name":"SIAM Journal on Imaging Sciences","volume":"74 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-10-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"136210799","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}