Andreas Alpers, Maximilian Fiedler, Peter Gritzmann, Fabian Klemm
{"title":"Turning Grain Maps into Diagrams","authors":"Andreas Alpers, Maximilian Fiedler, Peter Gritzmann, Fabian Klemm","doi":"10.1137/22m1491988","DOIUrl":"https://doi.org/10.1137/22m1491988","url":null,"abstract":"The present paper studies mathematical models for representing, imaging, and analyzing polycrystalline materials. We introduce various techniques for converting grain maps into diagram or tessellation representations that rely on constrained clustering. In particular, we show how to significantly accelerate the computation of generalized balanced power diagrams and how to extend it to allow for optimization over all relevant parameters. A comparison of the accuracy of the proposed approaches is given based on a three-dimensional real-world data set of voxels.","PeriodicalId":49528,"journal":{"name":"SIAM Journal on Imaging Sciences","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-02-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"136180941","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":"Steerable Near-Quadrature Filter Pairs in Three Dimensions.","authors":"Tommy M Tang, Hemant D Tagare","doi":"10.1137/21m143529x","DOIUrl":"https://doi.org/10.1137/21m143529x","url":null,"abstract":"<p><p>Steerable filter pairs that are near quadrature have many image processing applications. This paper proposes a new methodology for designing such filters. The key idea is to design steerable filters by minimizing a departure-from-quadrature function. These minimizing filter pairs are almost exactly in quadrature. The polar part of the filters is nonnegative, monotonic, and highly focused around an axis, and asymptotically the filters achieve exact quadrature. These results are established by exploiting a relation between the filters and generalized Hilbert matrices. These near-quadrature filters closely approximate three dimensional Gabor filters. We experimentally verify the asymptotic mathematical results and further demonstrate the use of these filter pairs by efficient calculation of local Fourier shell correlation of cryogenic electron microscopy.</p>","PeriodicalId":49528,"journal":{"name":"SIAM Journal on Imaging Sciences","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2022-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC9683347/pdf/nihms-1847751.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"40704943","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Linh Nguyen, Markus Haltmeier, Richard Kowar, Ngoc Do
{"title":"Analysis for Full-Field Photoacoustic Tomography with Variable Sound Speed.","authors":"Linh Nguyen, Markus Haltmeier, Richard Kowar, Ngoc Do","doi":"10.1137/21m1463409","DOIUrl":"10.1137/21m1463409","url":null,"abstract":"<p><p>Photoacoustic tomography (PAT) is a non-invasive imaging modality that requires recovering the initial data of the wave equation from certain measurements of the solution outside the object. In the standard PAT measurement setup, the used data consist of time-dependent signals measured on an observation surface. In contrast, the measured data from the recently invented full-field detection technique provide the solution of the wave equation on a spatial domain at a single instant in time. While reconstruction using classical PAT data has been extensively studied, not much is known for the full field PAT problem. In this paper, we build mathematical foundations of the latter problem for variable sound speed and settle its uniqueness and stability. Moreover, we introduce an exact inversion method using time-reversal and study its convergence. Our results demonstrate the suitability of both the full field approach and the proposed time-reversal technique for high resolution photoacoustic imaging.</p>","PeriodicalId":49528,"journal":{"name":"SIAM Journal on Imaging Sciences","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2022-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC10162777/pdf/nihms-1887591.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"9437974","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Centering Noisy Images with Application to Cryo-EM.","authors":"Ayelet Heimowitz, Nir Sharon, Amit Singer","doi":"10.1137/20m1365946","DOIUrl":"https://doi.org/10.1137/20m1365946","url":null,"abstract":"<p><p>We target the problem of estimating the center of mass of objects in noisy two-dimensional images. We assume that the noise dominates the image, and thus many standard approaches are vulnerable to estimation errors, e.g., the direct computation of the center of mass and the geometric median which is a robust alternative to the center of mass. In this paper, we define a novel surrogate function to the center of mass. We present a mathematical and numerical analysis of our method and show that it outperforms existing methods for estimating the center of mass of an object in various realistic scenarios. As a case study, we apply our centering method to data from single-particle cryo-electron microscopy (cryo-EM), where the goal is to reconstruct the three-dimensional structure of macromolecules. We show how to apply our approach for a better translational alignment of molecule images picked from experimental data. In this way, we facilitate the succeeding steps of reconstruction and streamline the entire cryo-EM pipeline, saving computational time and supporting resolution enhancement.</p>","PeriodicalId":49528,"journal":{"name":"SIAM Journal on Imaging Sciences","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2021-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC8813033/pdf/nihms-1739531.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"39591769","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Recovery of surfaces and functions in high dimensions: sampling theory and links to neural networks.","authors":"Qing Zou, Mathews Jacob","doi":"10.1137/20M1340654","DOIUrl":"10.1137/20M1340654","url":null,"abstract":"<p><p>Several imaging algorithms including patch-based image denoising, image time series recovery, and convolutional neural networks can be thought of as methods that exploit the manifold structure of signals. While the empirical performance of these algorithms is impressive, the understanding of recovery of the signals and functions that live on manifold is less understood. In this paper, we focus on the recovery of signals that live on a union of surfaces. In particular, we consider signals living on a union of smooth band-limited surfaces in high dimensions. We show that an exponential mapping transforms the data to a union of low-dimensional subspaces. Using this relation, we introduce a sampling theoretical framework for the recovery of smooth surfaces from few samples and the learning of functions living on smooth surfaces. The low-rank property of the features is used to determine the number of measurements needed to recover the surface. Moreover, the low-rank property of the features also provides an efficient approach, which resembles a neural network, for the local representation of multidimensional functions on the surface. The direct representation of such a function in high dimensions often suffers from the curse of dimensionality; the large number of parameters would translate to the need for extensive training data. The low-rank property of the features can significantly reduce the number of parameters, which makes the computational structure attractive for learning and inference from limited labeled training data.</p>","PeriodicalId":49528,"journal":{"name":"SIAM Journal on Imaging Sciences","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2021-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC8323788/pdf/nihms-1673810.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"39265387","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Spectral Embedding Norm: Looking Deep into the Spectrum of the Graph Laplacian.","authors":"Xiuyuan Cheng, Gal Mishne","doi":"10.1137/18m1283160","DOIUrl":"10.1137/18m1283160","url":null,"abstract":"<p><p>The extraction of clusters from a dataset which includes multiple clusters and a significant background component is a non-trivial task of practical importance. In image analysis this manifests for example in anomaly detection and target detection. The traditional spectral clustering algorithm, which relies on the leading <i>K</i> eigenvectors to detect <i>K</i> clusters, fails in such cases. In this paper we propose the <i>spectral embedding norm</i> which sums the squared values of the first <i>I</i> normalized eigenvectors, where <i>I</i> can be significantly larger than <i>K</i>. We prove that this quantity can be used to separate clusters from the background in unbalanced settings, including extreme cases such as outlier detection. The performance of the algorithm is not sensitive to the choice of <i>I</i>, and we demonstrate its application on synthetic and real-world remote sensing and neuroimaging datasets.</p>","PeriodicalId":49528,"journal":{"name":"SIAM Journal on Imaging Sciences","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2020-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC8204716/pdf/nihms-1594853.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"10320693","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Simplifying Transforms for General Elastic Metrics on the Space of Plane Curves.","authors":"Tom Needham, Sebastian Kurtek","doi":"10.1137/19m1265132","DOIUrl":"https://doi.org/10.1137/19m1265132","url":null,"abstract":"<p><p>In the shape analysis approach to computer vision problems, one treats shapes as points in an infinite-dimensional Riemannian manifold, thereby facilitating algorithms for statistical calculations such as geodesic distance between shapes and averaging of a collection of shapes. The performance of these algorithms depends heavily on the choice of the Riemannian metric. In the setting of plane curve shapes, attention has largely been focused on a two-parameter family of first order Sobolev metrics, referred to as elastic metrics. They are particularly useful due to the existence of simplifying coordinate transformations for particular parameter values, such as the well-known square-root velocity transform. In this paper, we extend the transformations appearing in the existing literature to a family of isometries, which take any elastic metric to the flat <i>L</i> <sup>2</sup> metric. We also extend the transforms to treat piecewise linear curves and demonstrate the existence of optimal matchings over the diffeomorphism group in this setting. We conclude the paper with multiple examples of shape geodesics for open and closed curves. We also show the benefits of our approach in a simple classification experiment.</p>","PeriodicalId":49528,"journal":{"name":"SIAM Journal on Imaging Sciences","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2020-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1137/19m1265132","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"39316946","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Anthony Yezzi, Ganesh Sundaramoorthi, Minas Benyamin
{"title":"Accelerated Optimization in the PDE Framework Formulations for the Active Contour Case.","authors":"Anthony Yezzi, Ganesh Sundaramoorthi, Minas Benyamin","doi":"10.1137/19m1304210","DOIUrl":"10.1137/19m1304210","url":null,"abstract":"<p><p>Following the seminal work of Nesterov, accelerated optimization methods have been used to powerfully boost the performance of first-order, gradient based parameter estimation in scenarios where second-order optimization strategies are either inapplicable or impractical. Not only does accelerated gradient descent converge considerably faster than traditional gradient descent, but it also performs a more robust local search of the parameter space by initially overshooting and then oscillating back as it settles into a final configuration, thereby selecting only local minimizers with a basis of attraction large enough to contain the initial overshoot. This behavior has made accelerated and stochastic gradient search methods particularly popular within the machine learning community. In their recent PNAS 2016 paper, <i>A Variational Perspective on Accelerated Methods in Optimization</i>, Wibisono, Wilson, and Jordan demonstrate how a broad class of accelerated schemes can be cast in a variational framework formulated around the Bregman divergence, leading to continuum limit ODEs. We show how their formulation may be further extended to infinite dimensional manifolds (starting here with the geometric space of curves and surfaces) by substituting the Bregman divergence with inner products on the tangent space and explicitly introducing a distributed mass model which evolves in conjunction with the object of interest during the optimization process. The coevolving mass model, which is introduced purely for the sake of endowing the optimization with helpful dynamics, also links the resulting class of accelerated PDE based optimization schemes to fluid dynamical formulations of optimal mass transport.</p>","PeriodicalId":49528,"journal":{"name":"SIAM Journal on Imaging Sciences","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2020-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC8320808/pdf/nihms-1676294.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"39265386","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Bin Li, Chenyang Shen, Yujie Chi, Ming Yang, Yifei Lou, Linghong Zhou, Xun Jia
{"title":"MULTI-ENERGY CONE-BEAM CT RECONSTRUCTION WITH A SPATIAL SPECTRAL NONLOCAL MEANS ALGORITHM.","authors":"Bin Li, Chenyang Shen, Yujie Chi, Ming Yang, Yifei Lou, Linghong Zhou, Xun Jia","doi":"10.1137/17M1123237","DOIUrl":"https://doi.org/10.1137/17M1123237","url":null,"abstract":"Multi-energy computed tomography (CT) is an emerging medical image modality with a number of potential applications in diagnosis and therapy. However, high system cost and technical barriers obstruct its step into routine clinical practice. In this study, we propose a framework to realize multi-energy cone beam CT (ME-CBCT) on the CBCT system that is widely available and has been routinely used for radiotherapy image guidance. In our method, a kVp switching technique is realized, which acquires x-ray projections with kVp levels cycling through a number of values. For this kVp-switching based ME-CBCT acquisition, x-ray projections of each energy channel are only a subset of all the acquired projections. This leads to an undersampling issue, posing challenges to the reconstruction problem. We propose a spatial spectral non-local means (ssNLM) method to reconstruct ME-CBCT, which employs image correlations along both spatial and spectral directions to suppress noisy and streak artifacts. To address the intensity scale difference at different energy channels, a histogram matching method is incorporated. Our method is different from conventionally used NLM methods in that spectral dimension is included, which helps to effectively remove streak artifacts appearing at different directions in images with different energy channels. Convergence analysis of our algorithm is provided. A comprehensive set of simulation and real experimental studies demonstrate feasibility of our ME-CBCT scheme and the capability of achieving superior image quality compared to conventional filtered backprojection-type (FBP) and NLM reconstruction methods.","PeriodicalId":49528,"journal":{"name":"SIAM Journal on Imaging Sciences","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2018-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1137/17M1123237","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"36568105","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Thomas P Matthews, Joemini Poudel, Lei Li, Lihong V Wang, Mark A Anastasio
{"title":"Parameterized joint reconstruction of the initial pressure and sound speed distributions for photoacoustic computed tomography.","authors":"Thomas P Matthews, Joemini Poudel, Lei Li, Lihong V Wang, Mark A Anastasio","doi":"10.1137/17M1153649","DOIUrl":"10.1137/17M1153649","url":null,"abstract":"<p><p>Accurate estimation of the initial pressure distribution in photoacoustic computed tomography (PACT) depends on knowledge of the sound speed distribution. However, the sound speed distribution is typically unknown. Further, the initial pressure and sound speed distributions cannot both, in general, be stably recovered from PACT measurements alone. In this work, a joint reconstruction (JR) method for the initial pressure distribution and a low-dimensional parameterized model of the sound speed distribution is proposed. By employing <i>a priori</i> information about the structure of the sound speed distribution, both the initial pressure and sound speed can be accurately recovered. The JR problem is solved by use of a proximal optimization method that allows constraints and non-smooth regularization functions for the initial pressure distribution. The gradients of the cost function with respect to the initial pressure and sound speed distributions are calculated by use of an adjoint state method that has the same per-iteration computational cost as calculating the gradient with respect to the initial pressure distribution alone. This approach is evaluated through 2D computer-simulation studies for a small animal imaging model and by application to experimental in vivo measurements of a mouse.</p>","PeriodicalId":49528,"journal":{"name":"SIAM Journal on Imaging Sciences","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2018-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6447310/pdf/nihms-1015338.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"37127613","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}