Showing 1–40 of 40 results for author: Buzzard, G T

ResSR: A Residual Approach to Super-Resolving Multispectral Images

Authors: Haley Duba-Sullivan, Emma J. Reid, Sophie Voisin, Charles A. Bouman, Gregery T. Buzzard

Abstract: Multispectral imaging sensors typically have wavelength-dependent resolution, which reduces the ability to distinguish small features in some spectral bands. Existing super-resolution methods upsample a multispectral image (MSI) to achieve a common resolution across all bands but are typically sensor-specific, computationally expensive, and may assume invariant image statistics across multiple len… ▽ More Multispectral imaging sensors typically have wavelength-dependent resolution, which reduces the ability to distinguish small features in some spectral bands. Existing super-resolution methods upsample a multispectral image (MSI) to achieve a common resolution across all bands but are typically sensor-specific, computationally expensive, and may assume invariant image statistics across multiple length scales. In this paper, we introduce ResSR, an efficient and modular residual-based method for super-resolving the lower-resolution bands of a multispectral image. ResSR uses singular value decomposition (SVD) to identify correlations across spectral bands and then applies a residual correction process that corrects only the high-spatial frequency components of the upsampled bands. The SVD formulation improves the conditioning and simplifies the super-resolution problem, and the residual method retains accurate low-spatial frequencies from the measured data while incorporating high-spatial frequency detail from the SVD solution. While ResSR is formulated as the solution to an optimization problem, we derive an approximate closed-form solution that is fast and accurate. We formulate ResSR for any number of distinct resolutions, enabling easy application to any MSI. In a series of experiments on simulated and measured Sentinel-2 MSIs, ResSR is shown to produce image quality comparable to or better than alternative algorithms. However, it is computationally faster and can run on larger images, making it useful for processing large data sets. △ Less

Submitted 23 August, 2024; originally announced August 2024.

Comments: Submitted to IEEE Transactions on Geoscience and Remote Sensing

Pixel-weighted Multi-pose Fusion for Metal Artifact Reduction in X-ray Computed Tomography

Authors: Diyu Yang, Craig A. J. Kemp, Soumendu Majee, Gregery T. Buzzard, Charles A. Bouman

Abstract: X-ray computed tomography (CT) reconstructs the internal morphology of a three dimensional object from a collection of projection images, most commonly using a single rotation axis. However, for objects containing dense materials like metal, the use of a single rotation axis may leave some regions of the object obscured by the metal, even though projections from other rotation axes (or poses) migh… ▽ More X-ray computed tomography (CT) reconstructs the internal morphology of a three dimensional object from a collection of projection images, most commonly using a single rotation axis. However, for objects containing dense materials like metal, the use of a single rotation axis may leave some regions of the object obscured by the metal, even though projections from other rotation axes (or poses) might contain complementary information that would better resolve these obscured regions. In this paper, we propose pixel-weighted Multi-pose Fusion to reduce metal artifacts by fusing the information from complementary measurement poses into a single reconstruction. Our method uses Multi-Agent Consensus Equilibrium (MACE), an extension of Plug-and-Play, as a framework for integrating projection data from different poses. A primary novelty of the proposed method is that the output of different MACE agents are fused in a pixel-weighted manner to minimize the effects of metal throughout the reconstruction. Using real CT data on an object with and without metal inserts, we demonstrate that the proposed pixel-weighted Multi-pose Fusion method significantly reduces metal artifacts relative to single-pose reconstructions. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: Submitted to IEEE MMSP 2024. arXiv admin note: substantial text overlap with arXiv:2209.07561

CLAMP: Majorized Plug-and-Play for Coherent 3D LIDAR Imaging

Authors: Tony G. Allen, David J. Rabb, Gregery T. Buzzard, Charles A. Bouman

Abstract: Coherent LIDAR uses a chirped laser pulse for 3D imaging of distant targets. However, existing coherent LIDAR image reconstruction methods do not account for the system's aperture, resulting in sub-optimal resolution. Moreover, these methods use majorization-minimization for computational efficiency, but do so without a theoretical treatment of convergence. In this paper, we present Coherent LID… ▽ More Coherent LIDAR uses a chirped laser pulse for 3D imaging of distant targets. However, existing coherent LIDAR image reconstruction methods do not account for the system's aperture, resulting in sub-optimal resolution. Moreover, these methods use majorization-minimization for computational efficiency, but do so without a theoretical treatment of convergence. In this paper, we present Coherent LIDAR Aperture Modeled Plug-and-Play (CLAMP) for multi-look coherent LIDAR image reconstruction. CLAMP uses multi-agent consensus equilibrium (a form of PnP) to combine a neural network denoiser with an accurate physics-based forward model. CLAMP introduces an FFT-based method to account for the effects of the aperture and uses majorization of the forward model for computational efficiency. We also formalize the use of majorization-minimization in consensus optimization problems and prove convergence to the exact consensus equilibrium solution. Finally, we apply CLAMP to synthetic and measured data to demonstrate its effectiveness in producing high-resolution, speckle-free, 3D imagery. △ Less

Submitted 19 June, 2024; originally announced June 2024.

MACE CT Reconstruction for Modular Material Decomposition from Energy Resolving Photon-Counting Data

Authors: Natalie M. Jadue, Madhuri Nagare, Jonathan S. Maltz, Gregery T. Buzzard, Charles A. Bouman

Abstract: X-ray computed tomography (CT) based on photon counting detectors (PCD) extends standard CT by counting detected photons in multiple energy bins. PCD data can be used to increase the contrast-to-noise ratio (CNR), increase spatial resolution, reduce radiation dose, reduce injected contrast dose, and compute a material decomposition using a specified set of basis materials. Current commercial and p… ▽ More X-ray computed tomography (CT) based on photon counting detectors (PCD) extends standard CT by counting detected photons in multiple energy bins. PCD data can be used to increase the contrast-to-noise ratio (CNR), increase spatial resolution, reduce radiation dose, reduce injected contrast dose, and compute a material decomposition using a specified set of basis materials. Current commercial and prototype clinical photon counting CT systems utilize PCD-CT reconstruction methods that either reconstruct from each spectral bin separately, or first create an estimate of a material sinogram using a specified set of basis materials and then reconstruct from these material sinograms. However, existing methods are not able to utilize simultaneously and in a modular fashion both the measured spectral information and advanced prior models in order to produce a material decomposition. We describe an efficient, modular framework for PCD-based CT reconstruction and material decomposition using on Multi-Agent Consensus Equilibrium (MACE). Our method employs a detector proximal map or agent that uses PCD measurements to update an estimate of the pathlength sinogram. We also create a prior agent in the form of a sinogram denoiser that enforces both physical and empirical knowledge about the material-decomposed sinogram. The sinogram reconstruction is computed using the MACE algorithm, which finds an equilibrium solution between the two agents, and the final image is reconstructed from the estimated sinogram. Importantly, the modularity of our method allows the two agents to be designed, implemented, and optimized independently. Our results on simulated data show a substantial (450%) CNR boost vs conventional maximum likelihood reconstruction when applied to a phantom used to evaluate low contrast detectability. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 11 pages, 10 figures, submitted to SPIE Medical Imaging 2024

Texture Matching GAN for CT Image Enhancement

Authors: Madhuri Nagare, Gregery T. Buzzard, Charles A. Bouman

Abstract: Deep neural networks (DNN) are commonly used to denoise and sharpen X-ray computed tomography (CT) images with the goal of reducing patient X-ray dosage while maintaining reconstruction quality. However, naive application of DNN-based methods can result in image texture that is undesirable in clinical applications. Alternatively, generative adversarial network (GAN) based methods can produce appro… ▽ More Deep neural networks (DNN) are commonly used to denoise and sharpen X-ray computed tomography (CT) images with the goal of reducing patient X-ray dosage while maintaining reconstruction quality. However, naive application of DNN-based methods can result in image texture that is undesirable in clinical applications. Alternatively, generative adversarial network (GAN) based methods can produce appropriate texture, but naive application of GANs can introduce inaccurate or even unreal image detail. In this paper, we propose a texture matching generative adversarial network (TMGAN) that enhances CT images while generating an image texture that can be matched to a target texture. We use parallel generators to separate anatomical features from the generated texture, which allows the GAN to be trained to match the desired texture without directly affecting the underlying CT image. We demonstrate that TMGAN generates enhanced image quality while also producing image texture that is desirable for clinical application. △ Less

Submitted 20 December, 2023; originally announced December 2023.

Comments: Submitted to IEEE Transactions on Medical Imaging

Generative Plug and Play: Posterior Sampling for Inverse Problems

Authors: Charles A. Bouman, Gregery T. Buzzard

Abstract: Over the past decade, Plug-and-Play (PnP) has become a popular method for reconstructing images using a modular framework consisting of a forward and prior model. The great strength of PnP is that an image denoiser can be used as a prior model while the forward model can be implemented using more traditional physics-based approaches. However, a limitation of PnP is that it reconstructs only a sing… ▽ More Over the past decade, Plug-and-Play (PnP) has become a popular method for reconstructing images using a modular framework consisting of a forward and prior model. The great strength of PnP is that an image denoiser can be used as a prior model while the forward model can be implemented using more traditional physics-based approaches. However, a limitation of PnP is that it reconstructs only a single deterministic image. In this paper, we introduce Generative Plug-and-Play (GPnP), a generalization of PnP to sample from the posterior distribution. As with PnP, GPnP has a modular framework using a physics-based forward model and an image denoising prior model. However, in GPnP these models are extended to become proximal generators, which sample from associated distributions. GPnP applies these proximal generators in alternation to produce samples from the posterior. We present experimental simulations using the well-known BM3D denoiser. Our results demonstrate that the GPnP method is robust, easy to implement, and produces intuitively reasonable samples from the posterior for sparse interpolation and tomographic reconstruction. Code to accompany this paper is available at https://github.com/gbuzzard/generative-pnp-allerton . △ Less

Submitted 12 June, 2023; originally announced June 2023.

Comments: 8 pages, submitted to 2023 IEEE Allerton Conference

MSC Class: 94A08; 68U10; 60J22

Dynamic DH-MBIR for Phase-Error Estimation from Streaming Digital-Holography Data

Authors: Ali G. Sheikh, Casey J. Pellizzari, Sherman J. Kisner, Gregery T. Buzzard, Charles A. Bouman

Abstract: Directed energy applications require the estimation of digital-holographic (DH) phase errors due to atmospheric turbulence in order to accurately focus the outgoing beam. These phase error estimates must be computed with very low latency to keep pace with changing atmospheric parameters, which requires that phase errors be estimated in a single shot of DH data. The digital holography model-based i… ▽ More Directed energy applications require the estimation of digital-holographic (DH) phase errors due to atmospheric turbulence in order to accurately focus the outgoing beam. These phase error estimates must be computed with very low latency to keep pace with changing atmospheric parameters, which requires that phase errors be estimated in a single shot of DH data. The digital holography model-based iterative reconstruction (DH-MBIR) algorithm is capable of accurately estimating phase errors in a single shot using the expectation maximization (EM) algorithm. However, existing implementations of DH-MBIR require hundreds of iterations, which is not practical for real-time applications. In this paper, we present the Dynamic DH-MBIR (DDH-MBIR) algorithm for estimating isoplanatic phase errors from streaming single-shot data with extremely low latency. The Dynamic DH-MBIR algorithm reduces the computation and latency by orders of magnitude relative to conventional DH-MBIR, making real-time throughput and latency feasible in applications. Using simulated data that models frozen flow of atmospheric turbulence, we show that our algorithm can achieve a consistently high Strehl ratio with realistic simulation parameters using only 1 iteration per timestep. △ Less

Submitted 5 May, 2023; originally announced May 2023.

Comments: Submitted to 2023 IEEE Asilomar Conference on Signals, Systems, and Computers

Projected Multi-Agent Consensus Equilibrium (PMACE) with Application to Ptychography

Authors: Qiuchen Zhai, Gregery T. Buzzard, Kevin Mertes, Brendt Wohlberg, Charles A. Bouman

Abstract: Multi-Agent Consensus Equilibrium (MACE) formulates an inverse imaging problem as a balance among multiple update agents such as data-fitting terms and denoisers. However, each such agent operates on a separate copy of the full image, leading to redundant memory use and slow convergence when each agent affects only a small subset of the full image. In this paper, we extend MACE to Projected Multi-… ▽ More Multi-Agent Consensus Equilibrium (MACE) formulates an inverse imaging problem as a balance among multiple update agents such as data-fitting terms and denoisers. However, each such agent operates on a separate copy of the full image, leading to redundant memory use and slow convergence when each agent affects only a small subset of the full image. In this paper, we extend MACE to Projected Multi-Agent Consensus Equilibrium (PMACE), in which each agent updates only a projected component of the full image, thus greatly reducing memory use for some applications.We describe PMACE in terms of an equilibrium problem and an equivalent fixed point problem and show that in most cases the PMACE equilibrium is not the solution of an optimization problem. To demonstrate the value of PMACE, we apply it to the problem of ptychography, in which a sample is reconstructed from the diffraction patterns resulting from coherent X-ray illumination at multiple overlapping spots. In our PMACE formulation, each spot corresponds to a separate data-fitting agent, with the final solution found as an equilibrium among all the agents. Our results demonstrate that the PMACE reconstruction algorithm generates more accurate reconstructions at a lower computational cost than existing ptychography algorithms when the spots are sparsely sampled. △ Less

Submitted 5 October, 2023; v1 submitted 27 March, 2023; originally announced March 2023.

Autonomous Polycrystalline Material Decomposition for Hyperspectral Neutron Tomography

Authors: Mohammad Samin Nur Chowdhury, Diyu Yang, Shimin Tang, Singanallur V. Venkatakrishnan, Hassina Z. Bilheux, Gregery T. Buzzard, Charles A. Bouman

Abstract: Hyperspectral neutron tomography is an effective method for analyzing crystalline material samples with complex compositions in a non-destructive manner. Since the counts in the hyperspectral neutron radiographs directly depend on the neutron cross-sections, materials may exhibit contrasting neutron responses across wavelengths. Therefore, it is possible to extract the unique signatures associated… ▽ More Hyperspectral neutron tomography is an effective method for analyzing crystalline material samples with complex compositions in a non-destructive manner. Since the counts in the hyperspectral neutron radiographs directly depend on the neutron cross-sections, materials may exhibit contrasting neutron responses across wavelengths. Therefore, it is possible to extract the unique signatures associated with each material and use them to separate the crystalline phases simultaneously. We introduce an autonomous material decomposition (AMD) algorithm to automatically characterize and localize polycrystalline structures using Bragg edges with contrasting neutron responses from hyperspectral data. The algorithm estimates the linear attenuation coefficient spectra from the measured radiographs and then uses these spectra to perform polycrystalline material decomposition and reconstructs 3D material volumes to localize materials in the spatial domain. Our results demonstrate that the method can accurately estimate both the linear attenuation coefficient spectra and associated reconstructions on both simulated and experimental neutron data. △ Less

Submitted 21 August, 2023; v1 submitted 27 February, 2023; originally announced February 2023.

X-ray Spectral Estimation using Dictionary Learning

Authors: Wenrui Li, Venkatesh Sridhar, K. Aditya Mohan, Saransh Singh, Jean-Baptiste Forien, Xin Liu, Gregery T. Buzzard, Charles A. Bouman

Abstract: As computational tools for X-ray computed tomography (CT) become more quantitatively accurate, knowledge of the source-detector spectral response is critical for quantitative system-independent reconstruction and material characterization capabilities. Directly measuring the spectral response of a CT system is hard, which motivates spectral estimation using transmission data obtained from a collec… ▽ More As computational tools for X-ray computed tomography (CT) become more quantitatively accurate, knowledge of the source-detector spectral response is critical for quantitative system-independent reconstruction and material characterization capabilities. Directly measuring the spectral response of a CT system is hard, which motivates spectral estimation using transmission data obtained from a collection of known homogeneous objects. However, the associated inverse problem is ill-conditioned, making accurate estimation of the spectrum challenging, particularly in the absence of a close initial guess. In this paper, we describe a dictionary-based spectral estimation method that yields accurate results without the need for any initial estimate of the spectral response. Our method utilizes a MAP estimation framework that combines a physics-based forward model along with an $L_0$ sparsity constraint and a simplex constraint on the dictionary coefficients. Our method uses a greedy support selection method and a new pair-wise iterated coordinate descent method to compute the above estimate. We demonstrate that our dictionary-based method outperforms a state-of-the-art method as shown in a cross-validation experiment on four real datasets collected at beamline 8.3.2 of the Advanced Light Source (ALS). △ Less

Submitted 26 February, 2023; originally announced February 2023.

Comments: Document Release Number: LLNL-CONF-845171 Submitted to 2023 ICIP conference

Ringing Artifact Reduction Method for Ultrasound Reconstruction Using Multi-Agent Consensus Equilibrium

Authors: Abdulrahman M. Alanazi, Singanallur Venkatakrishnan, Gregery T. Buzzard, Charles A. Bouman

Abstract: Non-destructive characterization of multi-layered structures that can be accessed from only a single side is important for applications such as well-bore integrity inspection. Existing methods related to Synthetic Aperture Focusing Technique (SAFT) rapidly produce acceptable results but with significant artifacts. Recently, ultrasound model-based iterative reconstruction (UMBIR) approaches have sh… ▽ More Non-destructive characterization of multi-layered structures that can be accessed from only a single side is important for applications such as well-bore integrity inspection. Existing methods related to Synthetic Aperture Focusing Technique (SAFT) rapidly produce acceptable results but with significant artifacts. Recently, ultrasound model-based iterative reconstruction (UMBIR) approaches have shown significant improvements over SAFT. However, even these methods produce ringing artifacts due to the high fractional-bandwidth of the excitation signal. In this paper, we propose a ringing artifact reduction method for ultrasound image reconstruction that uses a multi-agent consensus equilibrium (RARE-MACE) framework. Our approach integrates a physics-based forward model that accounts for the propagation of a collimated ultrasonic beam in multi-layered media, a spatially varying image prior, and a denoiser designed to suppress the ringing artifacts that are characteristic of reconstructions from high-fractional bandwidth ultrasound sensor data. We test our method on simulated and experimental measurements and show substantial improvements in image quality compared to SAFT and UMBIR. △ Less

Submitted 9 February, 2023; originally announced February 2023.

Comments: arXiv admin note: text overlap with arXiv:2211.15859, arXiv:2202.09703

An Edge Alignment-based Orientation Selection Method for Neutron Tomography

Authors: Diyu Yang, Shimin Tang, Singanallur V. Venkatakrishnan, Mohammad S. N. Chowdhury, Yuxuan Zhang, Hassina Z. Bilheux, Gregery T. Buzzard, Charles A. Bouman

Abstract: Neutron computed tomography (nCT) is a 3D characterization technique used to image the internal morphology or chemical composition of samples in biology and materials sciences. A typical workflow involves placing the sample in the path of a neutron beam, acquiring projection data at a predefined set of orientations, and processing the resulting data using an analytic reconstruction algorithm. Typi… ▽ More Neutron computed tomography (nCT) is a 3D characterization technique used to image the internal morphology or chemical composition of samples in biology and materials sciences. A typical workflow involves placing the sample in the path of a neutron beam, acquiring projection data at a predefined set of orientations, and processing the resulting data using an analytic reconstruction algorithm. Typical nCT scans require hours to days to complete and are then processed using conventional filtered back-projection (FBP), which performs poorly with sparse views or noisy data. Hence, the main methods in order to reduce overall acquisition time are the use of an improved sampling strategy combined with the use of advanced reconstruction methods such as model-based iterative reconstruction (MBIR). In this paper, we propose an adaptive orientation selection method in which an MBIR reconstruction on previously-acquired measurements is used to define an objective function on orientations that balances a data-fitting term promoting edge alignment and a regularization term promoting orientation diversity. Using simulated and experimental data, we demonstrate that our method produces high-quality reconstructions using significantly fewer total measurements than the conventional approach. △ Less

Submitted 8 March, 2023; v1 submitted 1 December, 2022; originally announced December 2022.

Model-based Reconstruction for Multi-Frequency Collimated Beam Ultrasound Systems

Authors: Abdulrahman M. Alanazi, Singanallur Venkatakrishnan, Hector Santos-Villalobos, Gregery T. Buzzard, Charles Bouman

Abstract: Collimated beam ultrasound systems are a technology for imaging inside multi-layered structures such as geothermal wells. These systems work by using a collimated narrow-band ultrasound transmitter that can penetrate through multiple layers of heterogeneous material. A series of measurements can then be made at multiple transmit frequencies. However, commonly used reconstruction algorithms such as… ▽ More Collimated beam ultrasound systems are a technology for imaging inside multi-layered structures such as geothermal wells. These systems work by using a collimated narrow-band ultrasound transmitter that can penetrate through multiple layers of heterogeneous material. A series of measurements can then be made at multiple transmit frequencies. However, commonly used reconstruction algorithms such as Synthetic Aperture Focusing Technique (SAFT) tend to produce poor quality reconstructions for these systems both because they do not model collimated beam systems and they do not jointly reconstruct the multiple frequencies. In this paper, we propose a multi-frequency ultrasound model-based iterative reconstruction (UMBIR) algorithm designed for multi-frequency collimated beam ultrasound systems. The combined system targets reflective imaging of heterogeneous, multi-layered structures. For each transmitted frequency band, we introduce a physics-based forward model to accurately account for the propagation of the collimated narrow-band ultrasonic beam through the multi-layered media. We then show how the joint multi-frequency UMBIR reconstruction can be computed by modeling the direct arrival signals, detector noise, and incorporating a spatially varying image prior. Results using both simulated and experimental data indicate that multi-frequency UMBIR reconstruction yields much higher reconstruction quality than either single frequency UMBIR or SAFT. △ Less

Submitted 28 November, 2022; originally announced November 2022.

Multi-Pose Fusion for Sparse-View CT Reconstruction Using Consensus Equilibrium

Authors: Diyu Yang, Craig A. J. Kemp, Gregery T. Buzzard, Charles A. Bouman

Abstract: CT imaging works by reconstructing an object of interest from a collection of projections. Traditional methods such as filtered-back projection (FBP) work on projection images acquired around a fixed rotation axis. However, for some CT problems, it is desirable to perform a joint reconstruction from projection data acquired from multiple rotation axes. In this paper, we present Multi-Pose Fusion… ▽ More CT imaging works by reconstructing an object of interest from a collection of projections. Traditional methods such as filtered-back projection (FBP) work on projection images acquired around a fixed rotation axis. However, for some CT problems, it is desirable to perform a joint reconstruction from projection data acquired from multiple rotation axes. In this paper, we present Multi-Pose Fusion, a novel algorithm that performs a joint tomographic reconstruction from CT scans acquired from multiple poses of a single object, where each pose has a distinct rotation axis. Our approach uses multi-agent consensus equilibrium (MACE), an extension of plug-and-play, as a framework for integrating projection data from different poses. We apply our method on simulated data and demonstrate that Multi-Pose Fusion can achieve a better reconstruction result than single pose reconstruction. △ Less

Submitted 15 September, 2022; originally announced September 2022.

Comments: To appear in 58th Annual Allerton Conference on Communication, Control, and Computing

Plug-and-Play Methods for Integrating Physical and Learned Models in Computational Imaging

Authors: Ulugbek S. Kamilov, Charles A. Bouman, Gregery T. Buzzard, Brendt Wohlberg

Abstract: Plug-and-Play Priors (PnP) is one of the most widely-used frameworks for solving computational imaging problems through the integration of physical models and learned models. PnP leverages high-fidelity physical sensor models and powerful machine learning methods for prior modeling of data to provide state-of-the-art reconstruction algorithms. PnP algorithms alternate between minimizing a data-fid… ▽ More Plug-and-Play Priors (PnP) is one of the most widely-used frameworks for solving computational imaging problems through the integration of physical models and learned models. PnP leverages high-fidelity physical sensor models and powerful machine learning methods for prior modeling of data to provide state-of-the-art reconstruction algorithms. PnP algorithms alternate between minimizing a data-fidelity term to promote data consistency and imposing a learned regularizer in the form of an image denoiser. Recent highly-successful applications of PnP algorithms include bio-microscopy, computerized tomography, magnetic resonance imaging, and joint ptycho-tomography. This article presents a unified and principled review of PnP by tracing its roots, describing its major variations, summarizing main results, and discussing applications in computational imaging. We also point the way towards further developments by discussing recent results on equilibrium equations that formulate the problem associated with PnP algorithms. △ Less

Submitted 12 August, 2022; v1 submitted 31 March, 2022; originally announced March 2022.

Sparse-View CT Reconstruction using Recurrent Stacked Back Projection

Authors: Wenrui Li, Gregery T. Buzzard, Charles A. Bouman

Abstract: Sparse-view CT reconstruction is important in a wide range of applications due to limitations on cost, acquisition time, or dosage. However, traditional direct reconstruction methods such as filtered back-projection (FBP) lead to low-quality reconstructions in the sub-Nyquist regime. In contrast, deep neural networks (DNNs) can produce high-quality reconstructions from sparse and noisy data, e.g.… ▽ More Sparse-view CT reconstruction is important in a wide range of applications due to limitations on cost, acquisition time, or dosage. However, traditional direct reconstruction methods such as filtered back-projection (FBP) lead to low-quality reconstructions in the sub-Nyquist regime. In contrast, deep neural networks (DNNs) can produce high-quality reconstructions from sparse and noisy data, e.g. through post-processing of FBP reconstructions, as can model-based iterative reconstruction (MBIR), albeit at a higher computational cost. In this paper, we introduce a direct-reconstruction DNN method called Recurrent Stacked Back Projection (RSBP) that uses sequentially-acquired backprojections of individual views as input to a recurrent convolutional LSTM network. The SBP structure maintains all information in the sinogram, while the recurrent processing exploits the correlations between adjacent views and produces an updated reconstruction after each new view. We train our network on simulated data and test on both simulated and real data and demonstrate that RSBP outperforms both DNN post-processing of FBP images and basic MBIR, with a lower computational cost than MBIR. △ Less

Submitted 9 December, 2021; originally announced December 2021.

Comments: 5 pages, 5 pages, 2021 Asilomar Conference on Signals, Systems, and Computers

Projected Multi-Agent Consensus Equilibrium for Ptychographic Image Reconstruction

Authors: Qiuchen Zhai, Brendt Wohlberg, Gregery T. Buzzard, Charles A. Bouman

Abstract: Ptychography is a computational imaging technique using multiple, overlapping, coherently illuminated snapshots to achieve nanometer resolution by solving a nonlinear phase-field recovery problem. Ptychography is vital for imaging of manufactured nanomaterials, but existing algorithms have computational shortcomings that limit large-scale application. In this paper, we present the Projected Multi-… ▽ More Ptychography is a computational imaging technique using multiple, overlapping, coherently illuminated snapshots to achieve nanometer resolution by solving a nonlinear phase-field recovery problem. Ptychography is vital for imaging of manufactured nanomaterials, but existing algorithms have computational shortcomings that limit large-scale application. In this paper, we present the Projected Multi-Agent Consensus Equilibrium (PMACE) approach for solving the ptychography inversion problem. This approach extends earlier work on MACE, which formulates an inversion problem as an equilibrium among multiple agents, each acting independently to update a full reconstruction. In PMACE, each agent acts on a portion (projection) corresponding to one of the snapshots, and these updates to projections are then combined to give an update to the full reconstruction. The resulting algorithm is easily parallelized, with convergence properties inherited from convergence results associated with MACE. We apply our method on simulated data and demonstrate that it outperforms competing algorithms in both reconstruction quality and convergence speed. △ Less

Submitted 8 December, 2021; v1 submitted 28 November, 2021; originally announced November 2021.

Comments: To be published in Asilomar Conference on Signals, Systems, and Computers 2021

Multi-Resolution Data Fusion for Super Resolution Imaging

Authors: Emma J Reid, Lawrence F Drummy, Charles A Bouman, Gregery T Buzzard

Abstract: Applications in materials and biological imaging are limited by the ability to collect high-resolution data over large areas in practical amounts of time. One solution to this problem is to collect low-resolution data and interpolate to produce a high-resolution image. However, most existing super-resolution algorithms are designed for natural images, often require aligned pairing of high and low-… ▽ More Applications in materials and biological imaging are limited by the ability to collect high-resolution data over large areas in practical amounts of time. One solution to this problem is to collect low-resolution data and interpolate to produce a high-resolution image. However, most existing super-resolution algorithms are designed for natural images, often require aligned pairing of high and low-resolution training data, and may not directly incorporate a model of the imaging sensor. In this paper, we present a Multi-resolution Data Fusion (MDF) algorithm for accurate interpolation of low-resolution data at multiple resolutions up to 8x. Our approach uses small quantities of unpaired high-resolution data to train a neural network prior model denoiser and then uses the Multi-Agent Consensus Equilibrium (MACE) problem formulation to balance this denoiser with a forward model agent that promotes fidelity to measured data. A key theoretical novelty is the analysis of mismatched back-projectors, which modify typical forward model updates for computational efficiency or improved image quality. We use MACE to prove that using a mismatched back-projector is equivalent to using a standard back-projector and an appropriately modified prior model. We present electron microscopy results at 4x and 8x interpolation factors that exhibit reduced artifacts relative to existing methods while maintaining fidelity to acquired data and accurately resolving sub-pixel-scale features. △ Less

Submitted 1 January, 2022; v1 submitted 13 May, 2021; originally announced May 2021.

Convergence of Griddy Gibbs Sampling and other perturbed Markov chains

Authors: Vu Dinh, Ann E. Rundell, Gregery T. Buzzard

Abstract: The Griddy Gibbs sampling was proposed by Ritter and Tanner (1992) as a computationally efficient approximation of the well-known Gibbs sampling method. The algorithm is simple and effective and has been used successfully to address problems in various fields of applied science. However, the approximate nature of the algorithm has prevented it from being widely used: the Markov chains generated by… ▽ More The Griddy Gibbs sampling was proposed by Ritter and Tanner (1992) as a computationally efficient approximation of the well-known Gibbs sampling method. The algorithm is simple and effective and has been used successfully to address problems in various fields of applied science. However, the approximate nature of the algorithm has prevented it from being widely used: the Markov chains generated by the Griddy Gibbs sampling method are not reversible in general, so the existence and uniqueness of its invariant measure is not guaranteed. Even when such an invariant measure uniquely exists, there was no estimate of the distance between it and the probability distribution of interest, hence no means to ensure the validity of the algorithm as a means to sample from the true distribution. In this paper, we show, subject to some fairly natural conditions, that the Griddy Gibbs method has a unique, invariant measure. Moreover, we provide $L^p$ estimates on the distance between this invariant measure and the corresponding measure obtained from Gibbs sampling. These results provide a theoretical foundation for the use of the Griddy Gibbs sampling method. We also address a more general result about the sensitivity of invariant measures under small perturbations on the transition probability. That is, if we replace the transition probability $P$ of any Monte Carlo Markov chain by another transition probability $Q$ where $Q$ is close to $P$, we can still estimate the distance between the two invariant measures. The distinguishing feature between our approach and previous work on convergence of perturbed Markov chain is that by considering the invariant measures as fixed points of linear operators on function spaces, we don't need to impose any further conditions on the rate of convergence of the Markov chain. △ Less

Submitted 29 March, 2021; originally announced March 2021.

Comments: 28 pages, 5 figures

MSC Class: 60J20; 65C40; 92B05; 32A70

Journal ref: Journal of Statistical Computation and Simulation 87.7 (2017): 1379-1400

Multi-Slice Fusion for Sparse-View and Limited-Angle 4D CT Reconstruction

Authors: Soumendu Majee, Thilo Balke, Craig A. J. Kemp, Gregery T. Buzzard, Charles A. Bouman

Abstract: Inverse problems spanning four or more dimensions such as space, time and other independent parameters have become increasingly important. State-of-the-art 4D reconstruction methods use model based iterative reconstruction (MBIR), but depend critically on the quality of the prior modeling. Recently, plug-and-play (PnP) methods have been shown to be an effective way to incorporate advanced prior mo… ▽ More Inverse problems spanning four or more dimensions such as space, time and other independent parameters have become increasingly important. State-of-the-art 4D reconstruction methods use model based iterative reconstruction (MBIR), but depend critically on the quality of the prior modeling. Recently, plug-and-play (PnP) methods have been shown to be an effective way to incorporate advanced prior models using state-of-the-art denoising algorithms. However, state-of-the-art denoisers such as BM4D and deep convolutional neural networks (CNNs) are primarily available for 2D or 3D images and extending them to higher dimensions is difficult due to algorithmic complexity and the increased difficulty of effective training. In this paper, we present multi-slice fusion, a novel algorithm for 4D reconstruction, based on the fusion of multiple low-dimensional denoisers. Our approach uses multi-agent consensus equilibrium (MACE), an extension of plug-and-play, as a framework for integrating the multiple lower-dimensional models. We apply our method to 4D cone-beam X-ray CT reconstruction for non destructive evaluation (NDE) of samples that are dynamically moving during acquisition. We implement multi-slice fusion on distributed, heterogeneous clusters in order to reconstruct large 4D volumes in reasonable time and demonstrate the inherent parallelizable nature of the algorithm. We present simulated and real experimental results on sparse-view and limited-angle CT data to demonstrate that multi-slice fusion can substantially improve the quality of reconstructions relative to traditional methods, while also being practical to implement and train. △ Less

Submitted 19 February, 2021; v1 submitted 31 July, 2020; originally announced August 2020.

Comments: arXiv admin note: substantial text overlap with arXiv:1906.06601

Distributed Iterative CT Reconstruction using Multi-Agent Consensus Equilibrium

Authors: Venkatesh Sridhar, Xiao Wang, Gregery T. Buzzard, Charles A. Bouman

Abstract: Model-Based Image Reconstruction (MBIR) methods significantly enhance the quality of computed tomographic (CT) reconstructions relative to analytical techniques, but are limited by high computational cost. In this paper, we propose a multi-agent consensus equilibrium (MACE) algorithm for distributing both the computation and memory of MBIR reconstruction across a large number of parallel nodes. In… ▽ More Model-Based Image Reconstruction (MBIR) methods significantly enhance the quality of computed tomographic (CT) reconstructions relative to analytical techniques, but are limited by high computational cost. In this paper, we propose a multi-agent consensus equilibrium (MACE) algorithm for distributing both the computation and memory of MBIR reconstruction across a large number of parallel nodes. In MACE, each node stores only a sparse subset of views and a small portion of the system matrix, and each parallel node performs a local sparse-view reconstruction, which based on repeated feedback from other nodes, converges to the global optimum. Our distributed approach can also incorporate advanced denoisers as priors to enhance reconstruction quality. In this case, we obtain a parallel solution to the serial framework of Plug-n-play (PnP) priors, which we call MACE-PnP. In order to make MACE practical, we introduce a partial update method that eliminates nested iterations and prove that it converges to the same global solution. Finally, we validate our approach on a distributed memory system with real CT data. We also demonstrate an implementation of our approach on a massive supercomputer that can perform large-scale reconstruction in real-time. △ Less

Submitted 20 November, 2019; originally announced November 2019.

Safe Approximate Dynamic Programming Via Kernelized Lipschitz Estimation

Authors: Ankush Chakrabarty, Devesh K. Jha, Gregery T. Buzzard, Yebin Wang, Kyriakos Vamvoudakis

Abstract: We develop a method for obtaining safe initial policies for reinforcement learning via approximate dynamic programming (ADP) techniques for uncertain systems evolving with discrete-time dynamics. We employ kernelized Lipschitz estimation and semidefinite programming for computing admissible initial control policies with provably high probability. Such admissible controllers enable safe initializat… ▽ More We develop a method for obtaining safe initial policies for reinforcement learning via approximate dynamic programming (ADP) techniques for uncertain systems evolving with discrete-time dynamics. We employ kernelized Lipschitz estimation and semidefinite programming for computing admissible initial control policies with provably high probability. Such admissible controllers enable safe initialization and constraint enforcement while providing exponential stability of the equilibrium of the closed-loop system. △ Less

Submitted 3 July, 2019; originally announced July 2019.

4D X-Ray CT Reconstruction using Multi-Slice Fusion

Authors: Soumendu Majee, Thilo Balke, Craig A. J. Kemp, Gregery T. Buzzard, Charles A. Bouman

Abstract: There is an increasing need to reconstruct objects in four or more dimensions corresponding to space, time and other independent parameters. The best 4D reconstruction algorithms use regularized iterative reconstruction approaches such as model based iterative reconstruction (MBIR), which depends critically on the quality of the prior modeling. Recently, Plug-and-Play methods have been shown to be… ▽ More There is an increasing need to reconstruct objects in four or more dimensions corresponding to space, time and other independent parameters. The best 4D reconstruction algorithms use regularized iterative reconstruction approaches such as model based iterative reconstruction (MBIR), which depends critically on the quality of the prior modeling. Recently, Plug-and-Play methods have been shown to be an effective way to incorporate advanced prior models using state-of-the-art denoising algorithms designed to remove additive white Gaussian noise (AWGN). However, state-of-the-art denoising algorithms such as BM4D and deep convolutional neural networks (CNNs) are primarily available for 2D and sometimes 3D images. In particular, CNNs are difficult and computationally expensive to implement in four or more dimensions, and training may be impossible if there is no associated high-dimensional training data. In this paper, we present Multi-Slice Fusion, a novel algorithm for 4D and higher-dimensional reconstruction, based on the fusion of multiple low-dimensional denoisers. Our approach uses multi-agent consensus equilibrium (MACE), an extension of Plug-and-Play, as a framework for integrating the multiple lower-dimensional prior models. We apply our method to the problem of 4D cone-beam X-ray CT reconstruction for Non Destructive Evaluation (NDE) of moving parts. This is done by solving the MACE equations using lower-dimensional CNN denoisers implemented in parallel on a heterogeneous cluster. Results on experimental CT data demonstrate that Multi-Slice Fusion can substantially improve the quality of reconstructions relative to traditional 4D priors, while also being practical to implement and train. △ Less

Submitted 15 June, 2019; originally announced June 2019.

Comments: 8 pages, 8 figures, IEEE International Conference on Computational Photography 2019, Tokyo

Plug and play methods for magnetic resonance imaging (long version)

Authors: Rizwan Ahmad, Charles A. Bouman, Gregery T. Buzzard, Stanley Chan, Sizhou Liu, Edward T. Reehorst, Philip Schniter

Abstract: Magnetic Resonance Imaging (MRI) is a non-invasive diagnostic tool that provides excellent soft-tissue contrast without the use of ionizing radiation. Compared to other clinical imaging modalities (e.g., CT or ultrasound), however, the data acquisition process for MRI is inherently slow, which motivates undersampling and thus drives the need for accurate, efficient reconstruction methods from unde… ▽ More Magnetic Resonance Imaging (MRI) is a non-invasive diagnostic tool that provides excellent soft-tissue contrast without the use of ionizing radiation. Compared to other clinical imaging modalities (e.g., CT or ultrasound), however, the data acquisition process for MRI is inherently slow, which motivates undersampling and thus drives the need for accurate, efficient reconstruction methods from undersampled datasets. In this article, we describe the use of "plug-and-play" (PnP) algorithms for MRI image recovery. We first describe the linearly approximated inverse problem encountered in MRI. Then we review several PnP methods, where the unifying commonality is to iteratively call a denoising subroutine as one step of a larger optimization-inspired algorithm. Next, we describe how the result of the PnP method can be interpreted as a solution to an equilibrium equation, allowing convergence analysis from the equilibrium perspective. Finally, we present illustrative examples of PnP methods applied to MRI image recovery. △ Less

Submitted 10 December, 2019; v1 submitted 20 March, 2019; originally announced March 2019.

Deep Back Projection for Sparse-View CT Reconstruction

Authors: Dong Hye Ye, Gregery T. Buzzard, Max Ruby, Charles A. Bouman

Abstract: Filtered back projection (FBP) is a classical method for image reconstruction from sinogram CT data. FBP is computationally efficient but produces lower quality reconstructions than more sophisticated iterative methods, particularly when the number of views is lower than the number required by the Nyquist rate. In this paper, we use a deep convolutional neural network (CNN) to produce high-quality… ▽ More Filtered back projection (FBP) is a classical method for image reconstruction from sinogram CT data. FBP is computationally efficient but produces lower quality reconstructions than more sophisticated iterative methods, particularly when the number of views is lower than the number required by the Nyquist rate. In this paper, we use a deep convolutional neural network (CNN) to produce high-quality reconstructions directly from sinogram data. A primary novelty of our approach is that we first back project each view separately to form a stack of back projections and then feed this stack as input into the convolutional neural network. These single-view back projections convert the encoding of sinogram data into the appropriate spatial location, which can then be leveraged by the spatial invariance of the CNN to learn the reconstruction effectively. We demonstrate the benefit of our CNN based back projection on simulated sparse-view CT data over classical FBP. △ Less

Submitted 6 July, 2018; originally announced July 2018.

Comments: GlobalSIP 2018

Deep neural networks for non-linear model-based ultrasound reconstruction

Authors: Hani Almansouri, S. V. Venkatakrishnan, Gregery T. Buzzard, Charles A. Bouman, Hector Santos-Villalobos

Abstract: Ultrasound reflection tomography is widely used to image large complex specimens that are only accessible from a single side, such as well systems and nuclear power plant containment walls. Typical methods for inverting the measurement rely on delay-and-sum algorithms that rapidly produce reconstructions but with significant artifacts. Recently, model-based reconstruction approaches using a linear… ▽ More Ultrasound reflection tomography is widely used to image large complex specimens that are only accessible from a single side, such as well systems and nuclear power plant containment walls. Typical methods for inverting the measurement rely on delay-and-sum algorithms that rapidly produce reconstructions but with significant artifacts. Recently, model-based reconstruction approaches using a linear forward model have been shown to significantly improve image quality compared to the conventional approach. However, even these techniques result in artifacts for complex objects because of the inherent non-linearity of the ultrasound forward model. In this paper, we propose a non-iterative model-based reconstruction method for inverting measurements that are based on non-linear forward models for ultrasound imaging. Our approach involves obtaining an approximate estimate of the reconstruction using a simple linear back-projection and training a deep neural network to refine this to the actual reconstruction. We apply our method to simulated ultrasound data and demonstrate dramatic improvements in image quality compared to the delay-and-sum approach and the linear model-based reconstruction approach. △ Less

Submitted 28 September, 2018; v1 submitted 3 July, 2018; originally announced July 2018.

Plug-and-Play Unplugged: Optimization Free Reconstruction using Consensus Equilibrium

Authors: Gregery T. Buzzard, Stanley H. Chan, Suhas Sreehari, Charles A. Bouman

Abstract: Regularized inversion methods for image reconstruction are used widely due to their tractability and ability to combine complex physical sensor models with useful regularity criteria. Such methods motivated the recently developed Plug-and-Play prior method, which provides a framework to use advanced denoising algorithms as regularizers in inversion. However, the need to formulate regularized inver… ▽ More Regularized inversion methods for image reconstruction are used widely due to their tractability and ability to combine complex physical sensor models with useful regularity criteria. Such methods motivated the recently developed Plug-and-Play prior method, which provides a framework to use advanced denoising algorithms as regularizers in inversion. However, the need to formulate regularized inversion as the solution to an optimization problem limits the possible regularity conditions and physical sensor models. In this paper, we introduce Consensus Equilibrium (CE), which generalizes regularized inversion to include a much wider variety of both forward components and prior components without the need for either to be expressed with a cost function. CE is based on the solution of a set of equilibrium equations that balance data fit and regularity. In this framework, the problem of MAP estimation in regularized inversion is replaced by the problem of solving these equilibrium equations, which can be approached in multiple ways. The key contribution of CE is to provide a novel framework for fusing multiple heterogeneous models of physical sensors or models learned from data. We describe the derivation of the CE equations and prove that the solution of the CE equations generalizes the standard MAP estimate under appropriate circumstances. We also discuss algorithms for solving the CE equations, including ADMM with a novel form of preconditioning and Newton's method. We give examples to illustrate consensus equilibrium and the convergence properties of these algorithms and demonstrate this method on some toy problems and on a denoising example in which we use an array of convolutional neural network denoisers, none of which is tuned to match the noise level in a noisy image but which in consensus can achieve a better result than any of them individually. △ Less

Submitted 13 June, 2018; v1 submitted 24 May, 2017; originally announced May 2017.

Comments: 17 pages

MSC Class: 94A08; 68U10

A Framework for Dynamic Image Sampling Based on Supervised Learning (SLADS)

Authors: G. M. Dilshan P. Godaliyadda, Dong Hye Ye, Michael D. Uchic, Michael A. Groeber, Gregery T. Buzzard, Charles A. Bouman

Abstract: Sparse sampling schemes have the potential to dramatically reduce image acquisition time while simultaneously reducing radiation damage to samples. However, for a sparse sampling scheme to be useful it is important that we are able to reconstruct the underlying object with sufficient clarity using the sparse measurements. In dynamic sampling, each new measurement location is selected based on info… ▽ More Sparse sampling schemes have the potential to dramatically reduce image acquisition time while simultaneously reducing radiation damage to samples. However, for a sparse sampling scheme to be useful it is important that we are able to reconstruct the underlying object with sufficient clarity using the sparse measurements. In dynamic sampling, each new measurement location is selected based on information obtained from previous measurements. Therefore, dynamic sampling schemes have the potential to dramatically reduce the number of measurements needed for high fidelity reconstructions. However, most existing dynamic sampling methods for point-wise measurement acquisition tend to be computationally expensive and are therefore too slow for practical applications. In this paper, we present a framework for dynamic sampling based on machine learning techniques, which we call a supervised learning approach for dynamic sampling (SLADS). In each step of SLADS, the objective is to find the pixel that maximizes the expected reduction in distortion (ERD) given previous measurements. SLADS is fast because we use a simple regression function to compute the ERD, and it is accurate because the regression function is trained using data sets that are representative of the specific application. In addition, we introduce a method to terminate dynamic sampling at a desired level of distortion, and we extended the SLADS methodology to sample groups of pixels at each step. Finally, we present results on computationally-generated synthetic data and experimentally-collected data to demonstrate a dramatic improvement over state-of-the-art static sampling methods. △ Less

Submitted 14 March, 2017; originally announced March 2017.

Simultaneous state and exogenous input estimation for nonlinear systems using boundary-layer sliding mode observers

Authors: Ankush Chakrabarty, Gregery T. Buzzard, Stanislaw H. Zak, Fanglai Zhu, Ann E. Rundell

Abstract: While sliding mode observers (SMOs) using discontinuous relays are widely analyzed, most SMOs are implemented computationally using a continuous approximation of the discontinuous relays. This approximation results in the formation of a boundary layer in a neighborhood of the sliding manifold in the observer error space. Therefore, it becomes necessary to develop methods for attenuating the effect… ▽ More While sliding mode observers (SMOs) using discontinuous relays are widely analyzed, most SMOs are implemented computationally using a continuous approximation of the discontinuous relays. This approximation results in the formation of a boundary layer in a neighborhood of the sliding manifold in the observer error space. Therefore, it becomes necessary to develop methods for attenuating the effect of the boundary layer and guaranteeing performance bounds on the resulting state estimation error. In this paper, a method is proposed for constructing boundary-layer SMOs (BL-SMOs) with prescribed state estimation error bounds. The BL-SMO formulation is then extended to simultaneously estimate exogenous inputs (disturbance signals in the state and output vector fields), along with the system state. Two numerical examples are presented to illustrate the effectiveness of the proposed approach. △ Less

Submitted 17 December, 2015; v1 submitted 14 July, 2015; originally announced July 2015.

Journal ref: SIAM J. Control Optim 2018

Experimental Design for Dynamics Identification of Cellular Processes

Authors: Vu Dinh, Ann E. Rundell, Gregery T. Buzzard

Abstract: We address the problem of using nonlinear models to design experiments to characterize the dynamics of cellular processes by using the approach of the Maximally Informative Next Experiment (MINE), which was introduced in [W. Dong, et al. Systems biology of the clock in neurospora crassa. {\em PLoS ONE}, page e3105, 2008] and independently in [M. M. Donahue, et al. Experiment design through dynamic… ▽ More We address the problem of using nonlinear models to design experiments to characterize the dynamics of cellular processes by using the approach of the Maximally Informative Next Experiment (MINE), which was introduced in [W. Dong, et al. Systems biology of the clock in neurospora crassa. {\em PLoS ONE}, page e3105, 2008] and independently in [M. M. Donahue, et al. Experiment design through dynamical characterization of non-linear systems biology models utilising sparse grids. {\em IET System Biology}, 4:249--262, 2010]. In this approach, existing data is used to define a probability distribution on the parameters; the next measurement point is the one that yields the largest model output variance with this distribution. Building upon this approach, we introduce the Expected Dynamics Estimator (EDE), which is the expected value using this distribution of the output as a function of time. We prove the consistency of this estimator (uniform convergence to true dynamics) even when the chosen experiments cluster in a finite set of points. We extend this proof of consistency to various practical assumptions on noisy data and moderate levels of model mismatch. Through the derivation and proof, we develop a relaxed version of MINE that is more computationally tractable and robust than the original formulation. The results are illustrated with numerical examples on two nonlinear ordinary differential equation models of biomolecular and cellular processes. △ Less

Submitted 13 November, 2013; originally announced November 2013.

MSC Class: 62K20; 65D15; 62L05; 92C45

Maps conjugating holomorphic maps in C^n

Authors: Gregery T. Buzzard, Sergei Merenkov

Abstract: If f is a bijection from C^n onto a complex manifold M, which conjugates every holomorphic map in C^n to an endomorphism in M, then we prove that f is necessarily biholomorphic or antibiholomorphic. This extends a result of A. Hinkkanen to higher dimensions. As a corollary, we prove that if there is an epimorphism from the semigroup of all holomorphic endomorphisms of C^n to the semigroup of hol… ▽ More If f is a bijection from C^n onto a complex manifold M, which conjugates every holomorphic map in C^n to an endomorphism in M, then we prove that f is necessarily biholomorphic or antibiholomorphic. This extends a result of A. Hinkkanen to higher dimensions. As a corollary, we prove that if there is an epimorphism from the semigroup of all holomorphic endomorphisms of C^n to the semigroup of holomorphic endomorphisms in M, or an epimorphism in the opposite direction for a doubly-transitive M, then it is given by conjugation by some biholomorphic or antibiholomorphic map. We show also that there are two unbounded domains in C^n with isomorphic endomorphism semigroups but which are neither biholomorphically nor antibiholomorphically equivalent. △ Less

Submitted 10 November, 2003; v1 submitted 12 January, 2003; originally announced January 2003.

Comments: 10 pages

MSC Class: 32H02; 32H50

Double Sections and Dominating Maps

Authors: Steven Shin-Yi Lu, Gregery T. Buzzard

Abstract: As is well-known, there exist nonconstant holomorphic maps from the plane into the Riemann sphere $\PP^1$ minus two points, the simplest example of which is an explicit realization of the uniformization map given by applying the exponential map and then composing with a Mobius transformation taking 0 and $\infty$ to the two given punctures. Likewise, we can map the plane into the sphere minus on… ▽ More As is well-known, there exist nonconstant holomorphic maps from the plane into the Riemann sphere $\PP^1$ minus two points, the simplest example of which is an explicit realization of the uniformization map given by applying the exponential map and then composing with a Mobius transformation taking 0 and $\infty$ to the two given punctures. Likewise, we can map the plane into the sphere minus one point by simply applying directly a Mobius transformation taking $\infty$ to this puncture. In this paper we prove two parametrized versions of this uniformization result using different methods. We first present a slightly easier version using complex analysis, which allows us to construct the uniformizing maps more or less explicitly in terms of the initial coordinates given. Then we give the more general and coordinate invariant version, which we prove by extending Kodaira's theory of the Jacobian fibration to a family of singular curves constructed via algebraic geometry. Finally, using the results obtained with the Jacobian fibration, we obtain two equivalent conditions for an complex surface nonhyperbolically fibered over a complex curve to be holomorphically dominable by $\CC^2$. In theorem \ref{thm:surf} we show that this dominability is equivalent to the apparently weaker condition of the existence of a Zariski dense image of $\CC$ and equivalent to the quasiprojectivity of the base curve together with the nonnegativity of the orbifold Euler characteristic. △ Less

Submitted 5 June, 2000; originally announced June 2000.

Comments: 21 pages

Report number: FI-NP1999-003

Algebraic Surfaces Holomorphically Dominable by C2

Authors: Stephen S. Y. Lu, Gregery T. Buzzard

Abstract: An n-dimensional complex manifold M is said to be (holomorphically) dominable by $\CC^n$ if there is a map $F:\CC^n \ra M$ which is holomorphic such that the Jacobian determinant $\det(DF)$ is not identically zero. Such a map F is called a dominating map. In this paper, we attempt to classify algebraic surfaces X which are dominable by $\CC^2$ using a combination of techniques from algebraic top… ▽ More An n-dimensional complex manifold M is said to be (holomorphically) dominable by $\CC^n$ if there is a map $F:\CC^n \ra M$ which is holomorphic such that the Jacobian determinant $\det(DF)$ is not identically zero. Such a map F is called a dominating map. In this paper, we attempt to classify algebraic surfaces X which are dominable by $\CC^2$ using a combination of techniques from algebraic topology, complex geometry and analysis. One of the key tools in the study of algebraic surfaces is the notion of Kodaira dimension (defined in section 2). By Kodaira's pioneering work and its extensions, an algebraic surface which is dominable by $\CC^2$ must have Kodaira dimension less than two. Using the Kodaira dimension and the fundamental group of X, we succeed in classifying algebraic surfaces which are dominable by $\CC^2$ except for certain cases in which X is an algebraic surface of Kodaira dimension zero and the case when X is rational without any logarithmic 1-form. More specifically, in the case when X is compact (namely projective), we need to exclude only the case when X is birationally equivalent to a K3 surface (a simply connected compact complex surface which admits a globally non-vanishing holomorphic 2-form) that is neither elliptic nor Kummer (see sections 3 and 4 for the definition of these types of surfaces). △ Less

Submitted 23 May, 2000; originally announced May 2000.

Comments: 39 pages

Report number: FI-NP1999-002

Double sections and dominating maps

Authors: Gregery T. Buzzard, Steven Lu

Abstract: As is well-known, given the complex sphere P^1 minus two points, there exist nonconstant holomorphic maps from the plane into this set, the simplest example of which is given by applying the exponential map and then composing with a Möbius transformation taking 0 and 1 to the two given punctures. Likewise, given the sphere minus one point, we can map the plane into this set by simply applying di… ▽ More As is well-known, given the complex sphere P^1 minus two points, there exist nonconstant holomorphic maps from the plane into this set, the simplest example of which is given by applying the exponential map and then composing with a Möbius transformation taking 0 and 1 to the two given punctures. Likewise, given the sphere minus one point, we can map the plane into this set by simply applying directly a Möbius transformation taking 1 to this puncture. In this paper we prove a parametrized version of this result. △ Less

Submitted 18 April, 1999; originally announced April 1999.

Report number: FI-NP1999-003

Algebraic surfaces holomorphically dominable by C^2

Authors: Gregery T. Buzzard, Stephen Lu

Abstract: Using the Kodaira dimension and the fundamental group of X, we succeed in classifying algebraic surfaces which are dominable by C^2 except for certain cases in which X is an algebraic surface of Kodaira dimension zero and the case when X is rational without any logarithmic 1-form. More specifically, in the case when X is compact (namely projective), we need to exclude only the case when X is bir… ▽ More Using the Kodaira dimension and the fundamental group of X, we succeed in classifying algebraic surfaces which are dominable by C^2 except for certain cases in which X is an algebraic surface of Kodaira dimension zero and the case when X is rational without any logarithmic 1-form. More specifically, in the case when X is compact (namely projective), we need to exclude only the case when X is birationally equivalent to a K3 surface (a simply connected compact complex surface which admits a globally non-vanishing holomorphic 2-form) that is neither elliptic nor Kummer. With the exceptions noted above, we show that for any algebraic surface of Kodaira dimension less than 2, dominability by C^2 is equivalent to the apparently weaker requirement of the existence of a holomorphic image of C which is Zariski dense in the surface. With the same exceptions, we will also show the very interesting and revealing fact that dominability by C^2 is preserved even if a sufficiently small neighborhood of any finite set of points is removed from the surface. In fact, we will provide a complete classification in the more general category of (not necessarily algebraic) compact complex surfaces before tackling the problem in the case of non-compact algebraic surfaces. △ Less

Submitted 13 March, 1999; originally announced March 1999.

Report number: FI-NP1999-002

A Carleman type theorem for proper holomorphic embeddings

Authors: Gregery T. Buzzard, Franc Forstneric

Abstract: In 1927, Carleman showed that a continuous, complex-valued function on the real line can be approximated in the Whitney topology by an entire function restricted to the real line. In this paper, we prove a similar result for proper holomorphic embeddings. Namely, we show that a proper $\cC^r$ embedding of the real line into $\C^n$ can be approximated in the strong $\cC^r$ topology by a proper ho… ▽ More In 1927, Carleman showed that a continuous, complex-valued function on the real line can be approximated in the Whitney topology by an entire function restricted to the real line. In this paper, we prove a similar result for proper holomorphic embeddings. Namely, we show that a proper $\cC^r$ embedding of the real line into $\C^n$ can be approximated in the strong $\cC^r$ topology by a proper holomorphic embedding of $\C$ into $\C^n$. △ Less

Submitted 3 April, 1996; originally announced April 1996.

Report number: MSRI 1996-021

An interpolation theorem for holomorphic automorphisms of {\bf C}$^n$

Authors: Gregery T. Buzzard, Franc Forstneric

Abstract: We construct automorphisms of $\C^n$ which map certain discrete sequences one onto another with prescribed finite jet at each point, thus solving a general Mittag-Leffler interpolation problem for automorphisms. Under certain circumstances, this can be done while also approximating a given automorphism on a compact set. We construct automorphisms of $\C^n$ which map certain discrete sequences one onto another with prescribed finite jet at each point, thus solving a general Mittag-Leffler interpolation problem for automorphisms. Under certain circumstances, this can be done while also approximating a given automorphism on a compact set. △ Less

Submitted 3 March, 1996; originally announced March 1996.

Report number: MSRI 1996-015

Complex dynamics in several variables

Authors: John Smillie, Gregery T. Buzzard

Abstract: These notes are the outgrowth of a series of lectures given at MSRI in January 1995 at the beginning of the special semester in complex dynamics and hyperbolic geometry. In these notes, the primary aim is to motivate the study of complex dynamics in two variables and to introduce the major ideas in the field. Hence the treatment of the subject is mostly expository. These notes are the outgrowth of a series of lectures given at MSRI in January 1995 at the beginning of the special semester in complex dynamics and hyperbolic geometry. In these notes, the primary aim is to motivate the study of complex dynamics in two variables and to introduce the major ideas in the field. Hence the treatment of the subject is mostly expository. △ Less

Submitted 7 February, 1996; originally announced February 1996.

Report number: MSRI 1996-011

Compositional roots of Hénon maps

Authors: Gregery T. Buzzard, John Erik Fornaess

Abstract: Let $H$ denote a composition of complex Hénon maps in {\bf C}$^2$. In this paper we show that the only possible compositional roots of $H$ are also compositions of Hénon maps, and that $H$ can have compositional roots of only finitely many distinct orders. Let $H$ denote a composition of complex Hénon maps in {\bf C}$^2$. In this paper we show that the only possible compositional roots of $H$ are also compositions of Hénon maps, and that $H$ can have compositional roots of only finitely many distinct orders. △ Less

Submitted 19 June, 1995; originally announced June 1995.

Report number: MSRI 1995-052

An embedding of {\bf C} in {\bf C}$^2$ with hyperbolic complement

Authors: Gregery T. Buzzard, John Erik Fornaess

Abstract: Let $X$ be a closed, $1$-dimensional, complex subvariety of $\CC^2$ and let $\ol{\BB}$ be a closed ball in $\CC^2 - X$. Then there exists a Fatou-Bieberbach domain $Ω$ with $X \subseteq Ω\subseteq \CC^2 - \ol{\BB}$ and a biholomorphic map $Φ: Ω\ra \CC^2$ such that $\CC^2 - Φ(X)$ is Kobayashi hyperbolic. As corollaries, there is an embedding of the plane in $\CC^2$ whose complement is hyperbolic,… ▽ More Let $X$ be a closed, $1$-dimensional, complex subvariety of $\CC^2$ and let $\ol{\BB}$ be a closed ball in $\CC^2 - X$. Then there exists a Fatou-Bieberbach domain $Ω$ with $X \subseteq Ω\subseteq \CC^2 - \ol{\BB}$ and a biholomorphic map $Φ: Ω\ra \CC^2$ such that $\CC^2 - Φ(X)$ is Kobayashi hyperbolic. As corollaries, there is an embedding of the plane in $\CC^2$ whose complement is hyperbolic, and there is a nontrivial Fatou-Bieberbach domain containing any finite collection of complex lines. △ Less

Submitted 19 June, 1995; originally announced June 1995.

Report number: MSRI 1995-051