Xueqi Guo - Academia.edu (original) (raw)

Papers by Xueqi Guo

Research paper thumbnail of SEAM-STRESS: A Weakly Supervised Framework for Interstitial Lung Disease Segmentation in Chest CT

2023 IEEE 20th International Symposium on Biomedical Imaging (ISBI)

Research paper thumbnail of FedFTN: Personalized Federated Learning with Deep Feature Transformation Network for Multi-institutional Low-count PET Denoising

arXiv (Cornell University), Apr 2, 2023

Research paper thumbnail of Characterization of Early Stage Parkinson's Disease From Resting-State fMRI Data Using a Long Short-Term Memory Network

Frontiers in Neuroimaging

Parkinson's disease (PD) is a common and complex neurodegenerative disorder with five stages ... more Parkinson's disease (PD) is a common and complex neurodegenerative disorder with five stages on the Hoehn and Yahr scaling. Characterizing brain function alterations with progression of early stage disease would support accurate disease staging, development of new therapies, and objective monitoring of disease progression or treatment response. Functional magnetic resonance imaging (fMRI) is a promising tool in revealing functional connectivity (FC) differences and developing biomarkers in PD. While fMRI and FC data have been utilized for diagnosis of PD through application of machine learning approaches such as support vector machine and logistic regression, the characterization of FC changes in early-stage PD has not been investigated. Given the complexity and non-linearity of fMRI data, we propose the use of a long short-term memory (LSTM) network to distinguish the early stages of PD and understand related functional brain changes. The study included 84 subjects (56 in stage...

Research paper thumbnail of Fast-MC-PET: A Novel Deep Learning-aided Motion Correction and Reconstruction Framework for Accelerated PET

arXiv (Cornell University), Feb 14, 2023

Patient motion during PET is inevitable. Its long acquisition time not only increases the motion ... more Patient motion during PET is inevitable. Its long acquisition time not only increases the motion and the associated artifacts but also the patient's discomfort, thus PET acceleration is desirable. However, accelerating PET acquisition will result in reconstructed images with low SNR, and the image quality will still be degraded by motion-induced artifacts. Most of the previous PET motion correction methods are motion type specific that require motion modeling, thus may fail when multiple types of motion present together. Also, those methods are customized for standard long acquisition and could not be directly applied to accelerated PET. To this end, modeling-free universal motion correction reconstruction for accelerated PET is still highly under-explored. In this work, we propose a novel deep learning-aided motion correction and reconstruction framework for accelerated PET, called Fast-MC-PET. Our framework consists of a universal motion correction (UMC) and a shortto-long acquisition reconstruction (SL-Reon) module. The UMC enables modeling-free motion correction by estimating quasi-continuous motion from ultra-short frame reconstructions and using this information for motion-compensated reconstruction. Then, the SL-Recon converts the accelerated UMC image with low counts to a high-quality image with high counts for our final reconstruction output. Our experimental results on human studies show that our Fast-MC-PET can enable 7-fold acceleration and use only 2 minutes acquisition to generate high-quality reconstruction images that outperform/match previous motion correction reconstruction methods using standard 15 minutes long acquisition data.

Research paper thumbnail of Generation of Whole-Body FDG Parametric K i Images From Static PET Images Using Deep Learning

IEEE Transactions on Radiation and Plasma Medical Sciences

F-fluorodeoxyglucose parametric K i images show a great advantage over static standard uptake val... more F-fluorodeoxyglucose parametric K i images show a great advantage over static standard uptake value (SUV) images, due to the higher contrast and better accuracy in tracer uptake rate estimation. In this study, we explored the feasibility of generating synthetic K i images from static SUV ratio (SUVR) images using three configurations of U-Nets with different sets of input and output image patches, which were the U-Nets with single input and single output (SISO), multiple inputs and single output (MISO), and single input and multiple outputs (SIMO). SUVR images were generated by averaging three 5-min dynamic SUV frames starting at 60-min post-injection, and then normalized by the mean SUV values in the blood pool. The corresponding ground-truth K i images were derived using Patlak graphical analysis with input functions from the measurement of arterial blood samples. Even though the synthetic K i values were not quantitatively accurate compared with ground truth, the linear regression analysis of joint histograms in the voxels of body regions showed that the mean R 2 values were higher between U-Net prediction and ground truth (0.596, 0.580, and 0.576 in SISO, MISO, and SIMO), than that between SUVR and ground truth K i (0.571). In terms of similarity metrics, the synthetic K i images were closer to the ground-truth K i images (mean SSIM = 0.729, 0.704, and 0.704 in SISO, MISO, and MISO) than the input SUVR images (mean SSIM = 0.691). Therefore, it is feasible to use deep learn

Research paper thumbnail of Federated Transfer Learning for Low-Dose PET Denoising: A Pilot Study With Simulated Heterogeneous Data

IEEE Transactions on Radiation and Plasma Medical Sciences

Positron emission tomography (PET) with a reduced injection dose, i.e., low-dose PET, is an effic... more Positron emission tomography (PET) with a reduced injection dose, i.e., low-dose PET, is an efficient way to reduce radiation dose. However, low-dose PET reconstruction suffers from a low signal-to-noise ratio (SNR), affecting diagnosis and other PET-related applications. Recently, deep learning-based PET denoising methods have demonstrated superior performance in generating high-quality reconstruction. However, these methods require a large amount of representative data for training, which can be difficult to collect and share due to medical data privacy regulations. Moreover, low-dose PET data at different institutions may use different low-dose protocols, leading to nonidentical data distribution. While previous federated learning (FL) algorithms enable multiinstitution collaborative training without the need of aggregating local data, it is challenging for previous methods to address the large domain shift caused by different low-dose PET settings, and the application of FL to PET is still underexplored. In this work, we propose a federated transfer learning (FTL) framework for low-dose PET denoising using heterogeneous low-dose data. Our experimental results on simulated multiinstitutional data demonstrate that our method can efficiently utilize heterogeneous low-dose data without compromising data privacy for achieving superior low-dose PET

Research paper thumbnail of Characterization of tunnel oxide passivated contact with n-type poly-Si on p-type c-Si wafer substrate

Current Applied Physics, 2019

The junction properties of tunnel silicon oxide (SiO x) passivated contact (TOPCon) with n-type p... more The junction properties of tunnel silicon oxide (SiO x) passivated contact (TOPCon) with n-type poly-Si on p-type c-Si wafer are characterized using current-voltage (J-V) and capacitance-voltage (C-V) measurements. The dark J-V curves show a standard diode characteristic with a turn-on voltage of ∼0.63V, indicating a p-n junction is formed. While the C-V curve displays an irregular shape with features of 1) a slow C increase with the decrease of the magnitude of reverse bias voltage, being used to estimate the built-in potential (V bi), 2) a significant increase at a given positive bias voltage, corresponding to the geometric capacitance crossing the ultrathin SiO x , and 3) a sharp decrease to

Research paper thumbnail of Inter-Pass Motion Correction for Whole-Body Dynamic PET and Parametric Imaging

IEEE Transactions on Radiation and Plasma Medical Sciences

Whole-body dynamic fluoro-D-glucose (FDG)positron emission tomography (PET) imaging through conti... more Whole-body dynamic fluoro-D-glucose (FDG)positron emission tomography (PET) imaging through continuous-bed-motion (CBM) mode multi-pass acquisition protocol is a promising metabolism measurement. However, inter-pass misalignment originating from body movement could degrade parametric quantification. We aim to apply a non-rigid registration method for inter-pass motion correction in whole-body dynamic PET. 27 subjects underwent a 90-min whole-body FDG CBM PET scan on a Biograph mCT (Siemens Healthineers), acquiring 9 over-the-heart single-bed passes and subsequently 19 CBM passes (frames). The inter-pass motion correction was executed using non-rigid image registration with multi-resolution, B-spline free-form deformations. The parametric images were then generated by Patlak analysis. The overlaid Patlak slope K i and y-intercept V b images were visualized to qualitatively evaluate motion impact and correction effect. The normalized weighted mean-squared Patlak fitting errors (NFEs) were compared in the whole body, head, and hypermetabolic regions of interest (ROIs). In K i images, ROI statistics were collected and malignancy discrimination capacity was estimated by the area under the receiver operating characteristic curve

Research paper thumbnail of Assessment of lower extremities flow using dynamic Rb-82 PET: Acquisition protocols and quantification methods

The Journal of Nuclear Medicine, May 1, 2021

Research paper thumbnail of Early Disease Stage Characterization in Parkinson’s Disease from Resting-state fMRI Data Using a Long Short-term Memory Network

Parkinson’s disease (PD) is a common and complex neurodegenerative disorder with 5 stages in the ... more Parkinson’s disease (PD) is a common and complex neurodegenerative disorder with 5 stages in the Hoehn and Yahr scaling. Given the heterogeneity of PD, it is challenging to classify early stages 1 and 2 and detect brain function alterations. Functional magnetic resonance imaging (fMRI) is a promising tool in revealing functional connectivity (FC) differences and developing biomarkers in PD. Some machine learning approaches like support vector machine and logistic regression have been successfully applied in the early diagnosis of PD using fMRI data, which outperform classifiers based on manually selected morphological features. However, the early-stage characterization in FC changes has not been fully investigated. Given the complexity and non-linearity of fMRI data, we propose the use of a long short-term memory (LSTM) network to characterize the early stages of PD. The study included 84 subjects (56 in stage 2 and 28 in stage 1) from the Parkinson’s Progression Markers Initiative ...

Research paper thumbnail of Generalized prediction framework for reconstructed image properties using neural networks

Medical Imaging 2019: Physics of Medical Imaging, Mar 1, 2019

Model-based reconstruction (MBR) algorithms in CT have demonstrated superior dose-image quality t... more Model-based reconstruction (MBR) algorithms in CT have demonstrated superior dose-image quality tradeoffs compared to traditional analytical methods. However, the nonlinear and data-dependent nature of these algorithms pose significant challenges for performance evaluation and parameter optimization. To address these challenges, this work presents an analysis framework for quantitative and predictive modeling of image properties in general nonlinear MBR algorithms. We propose to characterize the reconstructed appearance of arbitrary stimuli by the generalized system response function that accounts for dependence on the imaging conditions, reconstruction parameters, object, and the stimulus itself (size, contrast, location). We estimate this nonlinear function using a multilayer perceptron neural network by providing input and output pairs that samples the range of imaging parameters of interest. The feasibility of this approach was demonstrated for predicting the appearance of a spiculated lesion reconstructed by a penalized-likelihood objective with a Huber penalty in a physical phantom as a function of its location and reconstruction parameters β and δ. The generalized system response functions predicted from the trained neural network show good agreement with those computed from mean reconstructions, proving the ability of the framework in mapping out the nonlinear function for combinations of imaging parameters not present in the training data. We demonstrated utility of the framework to achieve desirable (e.g., non-blocky) lesion appearance in arbitrary locations in the phantom without the need for performing actual reconstructions. The proposed prediction framework permits efficient and quantifiable performance evaluations to provide robust control and understanding of image properties for general classes of nonlinear MBR algorithms.

Research paper thumbnail of Performance analysis for nonlinear tomographic data processing

15th International Meeting on Fully Three-Dimensional Image Reconstruction in Radiology and Nuclear Medicine, 2019

Image quality analysis of nonlinear algorithms is challenging due to numerous dependencies on the... more Image quality analysis of nonlinear algorithms is challenging due to numerous dependencies on the imaging system, algorithmic parameters, object, and stimulus. In particular, traditional notions of linearity and local linearity are of limited utility when the system response is dependent on the stimulus itself. In this work, we analyze the performance of nonlinear systems using perturbation response - the difference between the mean output with and without a stimulus, and introduce a new metric to examine variation of the responses in individual images. We applied the analysis to four algorithms with different degrees of nonlinearity for a spherical stimulus of varying contrast. For model-based reconstruction methods [penalized-likelihood (PL) reconstruction with a quadratic penalty and a Huber penalty], perturbation response analysis reaffirmed known trends in terms of object- and location-dependence. For a CNN denoising network, the response exhibits highly nonlinear behavior as the contrast increases - from the stimulus completely disappearing, to appearing at the right contrast but smaller in size, to being fully admitted by the algorithm. Furthermore, the variation metric for PL reconstruction with a Huber penalty and the CNN network reveals high variation at the edge of the stimulus, i.e., perturbation response computed from the mean images is a smoothed version of individual responses due to "jitter" in edges. This behavior suggests that the mean response alone may not be representative of performance in individual images and image quality metrics traditionally defined based on the mean response may be inappropriate for certain nonlinear algorithms. This work demonstrates the potential utility of perturbation response and response variation in the analysis and optimization of nonlinear imaging algorithms.

Research paper thumbnail of SEAM-STRESS: A Weakly Supervised Framework for Interstitial Lung Disease Segmentation in Chest CT

2023 IEEE 20th International Symposium on Biomedical Imaging (ISBI)

Research paper thumbnail of FedFTN: Personalized Federated Learning with Deep Feature Transformation Network for Multi-institutional Low-count PET Denoising

arXiv (Cornell University), Apr 2, 2023

Research paper thumbnail of Characterization of Early Stage Parkinson's Disease From Resting-State fMRI Data Using a Long Short-Term Memory Network

Frontiers in Neuroimaging

Parkinson's disease (PD) is a common and complex neurodegenerative disorder with five stages ... more Parkinson's disease (PD) is a common and complex neurodegenerative disorder with five stages on the Hoehn and Yahr scaling. Characterizing brain function alterations with progression of early stage disease would support accurate disease staging, development of new therapies, and objective monitoring of disease progression or treatment response. Functional magnetic resonance imaging (fMRI) is a promising tool in revealing functional connectivity (FC) differences and developing biomarkers in PD. While fMRI and FC data have been utilized for diagnosis of PD through application of machine learning approaches such as support vector machine and logistic regression, the characterization of FC changes in early-stage PD has not been investigated. Given the complexity and non-linearity of fMRI data, we propose the use of a long short-term memory (LSTM) network to distinguish the early stages of PD and understand related functional brain changes. The study included 84 subjects (56 in stage...

Research paper thumbnail of Fast-MC-PET: A Novel Deep Learning-aided Motion Correction and Reconstruction Framework for Accelerated PET

arXiv (Cornell University), Feb 14, 2023

Patient motion during PET is inevitable. Its long acquisition time not only increases the motion ... more Patient motion during PET is inevitable. Its long acquisition time not only increases the motion and the associated artifacts but also the patient's discomfort, thus PET acceleration is desirable. However, accelerating PET acquisition will result in reconstructed images with low SNR, and the image quality will still be degraded by motion-induced artifacts. Most of the previous PET motion correction methods are motion type specific that require motion modeling, thus may fail when multiple types of motion present together. Also, those methods are customized for standard long acquisition and could not be directly applied to accelerated PET. To this end, modeling-free universal motion correction reconstruction for accelerated PET is still highly under-explored. In this work, we propose a novel deep learning-aided motion correction and reconstruction framework for accelerated PET, called Fast-MC-PET. Our framework consists of a universal motion correction (UMC) and a shortto-long acquisition reconstruction (SL-Reon) module. The UMC enables modeling-free motion correction by estimating quasi-continuous motion from ultra-short frame reconstructions and using this information for motion-compensated reconstruction. Then, the SL-Recon converts the accelerated UMC image with low counts to a high-quality image with high counts for our final reconstruction output. Our experimental results on human studies show that our Fast-MC-PET can enable 7-fold acceleration and use only 2 minutes acquisition to generate high-quality reconstruction images that outperform/match previous motion correction reconstruction methods using standard 15 minutes long acquisition data.

Research paper thumbnail of Generation of Whole-Body FDG Parametric K i Images From Static PET Images Using Deep Learning

IEEE Transactions on Radiation and Plasma Medical Sciences

F-fluorodeoxyglucose parametric K i images show a great advantage over static standard uptake val... more F-fluorodeoxyglucose parametric K i images show a great advantage over static standard uptake value (SUV) images, due to the higher contrast and better accuracy in tracer uptake rate estimation. In this study, we explored the feasibility of generating synthetic K i images from static SUV ratio (SUVR) images using three configurations of U-Nets with different sets of input and output image patches, which were the U-Nets with single input and single output (SISO), multiple inputs and single output (MISO), and single input and multiple outputs (SIMO). SUVR images were generated by averaging three 5-min dynamic SUV frames starting at 60-min post-injection, and then normalized by the mean SUV values in the blood pool. The corresponding ground-truth K i images were derived using Patlak graphical analysis with input functions from the measurement of arterial blood samples. Even though the synthetic K i values were not quantitatively accurate compared with ground truth, the linear regression analysis of joint histograms in the voxels of body regions showed that the mean R 2 values were higher between U-Net prediction and ground truth (0.596, 0.580, and 0.576 in SISO, MISO, and SIMO), than that between SUVR and ground truth K i (0.571). In terms of similarity metrics, the synthetic K i images were closer to the ground-truth K i images (mean SSIM = 0.729, 0.704, and 0.704 in SISO, MISO, and MISO) than the input SUVR images (mean SSIM = 0.691). Therefore, it is feasible to use deep learn

Research paper thumbnail of Federated Transfer Learning for Low-Dose PET Denoising: A Pilot Study With Simulated Heterogeneous Data

IEEE Transactions on Radiation and Plasma Medical Sciences

Positron emission tomography (PET) with a reduced injection dose, i.e., low-dose PET, is an effic... more Positron emission tomography (PET) with a reduced injection dose, i.e., low-dose PET, is an efficient way to reduce radiation dose. However, low-dose PET reconstruction suffers from a low signal-to-noise ratio (SNR), affecting diagnosis and other PET-related applications. Recently, deep learning-based PET denoising methods have demonstrated superior performance in generating high-quality reconstruction. However, these methods require a large amount of representative data for training, which can be difficult to collect and share due to medical data privacy regulations. Moreover, low-dose PET data at different institutions may use different low-dose protocols, leading to nonidentical data distribution. While previous federated learning (FL) algorithms enable multiinstitution collaborative training without the need of aggregating local data, it is challenging for previous methods to address the large domain shift caused by different low-dose PET settings, and the application of FL to PET is still underexplored. In this work, we propose a federated transfer learning (FTL) framework for low-dose PET denoising using heterogeneous low-dose data. Our experimental results on simulated multiinstitutional data demonstrate that our method can efficiently utilize heterogeneous low-dose data without compromising data privacy for achieving superior low-dose PET

Research paper thumbnail of Characterization of tunnel oxide passivated contact with n-type poly-Si on p-type c-Si wafer substrate

Current Applied Physics, 2019

The junction properties of tunnel silicon oxide (SiO x) passivated contact (TOPCon) with n-type p... more The junction properties of tunnel silicon oxide (SiO x) passivated contact (TOPCon) with n-type poly-Si on p-type c-Si wafer are characterized using current-voltage (J-V) and capacitance-voltage (C-V) measurements. The dark J-V curves show a standard diode characteristic with a turn-on voltage of ∼0.63V, indicating a p-n junction is formed. While the C-V curve displays an irregular shape with features of 1) a slow C increase with the decrease of the magnitude of reverse bias voltage, being used to estimate the built-in potential (V bi), 2) a significant increase at a given positive bias voltage, corresponding to the geometric capacitance crossing the ultrathin SiO x , and 3) a sharp decrease to

Research paper thumbnail of Inter-Pass Motion Correction for Whole-Body Dynamic PET and Parametric Imaging

IEEE Transactions on Radiation and Plasma Medical Sciences

Whole-body dynamic fluoro-D-glucose (FDG)positron emission tomography (PET) imaging through conti... more Whole-body dynamic fluoro-D-glucose (FDG)positron emission tomography (PET) imaging through continuous-bed-motion (CBM) mode multi-pass acquisition protocol is a promising metabolism measurement. However, inter-pass misalignment originating from body movement could degrade parametric quantification. We aim to apply a non-rigid registration method for inter-pass motion correction in whole-body dynamic PET. 27 subjects underwent a 90-min whole-body FDG CBM PET scan on a Biograph mCT (Siemens Healthineers), acquiring 9 over-the-heart single-bed passes and subsequently 19 CBM passes (frames). The inter-pass motion correction was executed using non-rigid image registration with multi-resolution, B-spline free-form deformations. The parametric images were then generated by Patlak analysis. The overlaid Patlak slope K i and y-intercept V b images were visualized to qualitatively evaluate motion impact and correction effect. The normalized weighted mean-squared Patlak fitting errors (NFEs) were compared in the whole body, head, and hypermetabolic regions of interest (ROIs). In K i images, ROI statistics were collected and malignancy discrimination capacity was estimated by the area under the receiver operating characteristic curve

Research paper thumbnail of Assessment of lower extremities flow using dynamic Rb-82 PET: Acquisition protocols and quantification methods

The Journal of Nuclear Medicine, May 1, 2021

Research paper thumbnail of Early Disease Stage Characterization in Parkinson’s Disease from Resting-state fMRI Data Using a Long Short-term Memory Network

Parkinson’s disease (PD) is a common and complex neurodegenerative disorder with 5 stages in the ... more Parkinson’s disease (PD) is a common and complex neurodegenerative disorder with 5 stages in the Hoehn and Yahr scaling. Given the heterogeneity of PD, it is challenging to classify early stages 1 and 2 and detect brain function alterations. Functional magnetic resonance imaging (fMRI) is a promising tool in revealing functional connectivity (FC) differences and developing biomarkers in PD. Some machine learning approaches like support vector machine and logistic regression have been successfully applied in the early diagnosis of PD using fMRI data, which outperform classifiers based on manually selected morphological features. However, the early-stage characterization in FC changes has not been fully investigated. Given the complexity and non-linearity of fMRI data, we propose the use of a long short-term memory (LSTM) network to characterize the early stages of PD. The study included 84 subjects (56 in stage 2 and 28 in stage 1) from the Parkinson’s Progression Markers Initiative ...

Research paper thumbnail of Generalized prediction framework for reconstructed image properties using neural networks

Medical Imaging 2019: Physics of Medical Imaging, Mar 1, 2019

Model-based reconstruction (MBR) algorithms in CT have demonstrated superior dose-image quality t... more Model-based reconstruction (MBR) algorithms in CT have demonstrated superior dose-image quality tradeoffs compared to traditional analytical methods. However, the nonlinear and data-dependent nature of these algorithms pose significant challenges for performance evaluation and parameter optimization. To address these challenges, this work presents an analysis framework for quantitative and predictive modeling of image properties in general nonlinear MBR algorithms. We propose to characterize the reconstructed appearance of arbitrary stimuli by the generalized system response function that accounts for dependence on the imaging conditions, reconstruction parameters, object, and the stimulus itself (size, contrast, location). We estimate this nonlinear function using a multilayer perceptron neural network by providing input and output pairs that samples the range of imaging parameters of interest. The feasibility of this approach was demonstrated for predicting the appearance of a spiculated lesion reconstructed by a penalized-likelihood objective with a Huber penalty in a physical phantom as a function of its location and reconstruction parameters β and δ. The generalized system response functions predicted from the trained neural network show good agreement with those computed from mean reconstructions, proving the ability of the framework in mapping out the nonlinear function for combinations of imaging parameters not present in the training data. We demonstrated utility of the framework to achieve desirable (e.g., non-blocky) lesion appearance in arbitrary locations in the phantom without the need for performing actual reconstructions. The proposed prediction framework permits efficient and quantifiable performance evaluations to provide robust control and understanding of image properties for general classes of nonlinear MBR algorithms.

Research paper thumbnail of Performance analysis for nonlinear tomographic data processing

15th International Meeting on Fully Three-Dimensional Image Reconstruction in Radiology and Nuclear Medicine, 2019

Image quality analysis of nonlinear algorithms is challenging due to numerous dependencies on the... more Image quality analysis of nonlinear algorithms is challenging due to numerous dependencies on the imaging system, algorithmic parameters, object, and stimulus. In particular, traditional notions of linearity and local linearity are of limited utility when the system response is dependent on the stimulus itself. In this work, we analyze the performance of nonlinear systems using perturbation response - the difference between the mean output with and without a stimulus, and introduce a new metric to examine variation of the responses in individual images. We applied the analysis to four algorithms with different degrees of nonlinearity for a spherical stimulus of varying contrast. For model-based reconstruction methods [penalized-likelihood (PL) reconstruction with a quadratic penalty and a Huber penalty], perturbation response analysis reaffirmed known trends in terms of object- and location-dependence. For a CNN denoising network, the response exhibits highly nonlinear behavior as the contrast increases - from the stimulus completely disappearing, to appearing at the right contrast but smaller in size, to being fully admitted by the algorithm. Furthermore, the variation metric for PL reconstruction with a Huber penalty and the CNN network reveals high variation at the edge of the stimulus, i.e., perturbation response computed from the mean images is a smoothed version of individual responses due to "jitter" in edges. This behavior suggests that the mean response alone may not be representative of performance in individual images and image quality metrics traditionally defined based on the mean response may be inappropriate for certain nonlinear algorithms. This work demonstrates the potential utility of perturbation response and response variation in the analysis and optimization of nonlinear imaging algorithms.