Search NASA⌕ Search

SEARCH · Search NASA

Results for “image simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,027 records · Page 57

SU-F-T-205: Effectiveness of Robust Treatment Planning to Account for Inter- Fractional Variation in Intensity Modulated Proton Therapy for Head Neck Cancer

Purpose: To evaluate the potential benefits of robust optimization in intensity modulated proton therapy(IMPT) treatment planning to account for inter-fractional variation for Head Neck Cancer(HNC). Methods: One patient with bilateral HNC previous treated at our institution was used in this study. Ten daily CBCTs were selected. The CT numbers of the CBCTs were corrected by mapping the CT numbers from simulation CT via Deformable Image Registration. The planning target volumes(PTVs) were defined by a 3mm expansion from clinical target volumes(CTVs). The prescription was 70Gy, 54Gy to CTV1, CTV2, and PTV1, PTV2 for robust optimized(RO) and conventionally optimized(CO) plans respectively. Both techniques were generated by RayStation with the same beam angles: two anterior oblique and two posterior oblique angles. The similar dose constraints were used to achieve 99% of CTV1 received 100% prescription dose while kept the hotspots less than 110% of the prescription. In order to evaluate the dosimetric result through the course of treatment, the contours were deformed from simulation CT to daily CBCTs, modified, and approved by a radiation oncologist. The initial plan on the simulation CT was re-replayed on the daily CBCTs followed the bony alignment. The target coverage was evaluated using the daily doses and the cumulative dose. Results: Eight of 10 daily deliveries with using RO plan achieved at least 95% prescription dose to CTV1 and CTV2, while still kept maximum hotspot less than 112% of prescription compared with only one of 10 for the CO plan to achieve the same standards. For the cumulative doses, the target coverage for both RO and CO plans was quite similar, which was due to the compensation of cold and hot spots. Conclusion: Robust optimization can be effectively applied to compensate for target dose deficit caused by inter-fractional target geometric variation in IMPT treatment planning.

60 APPLIED LIFE SCIENCES↗

IFSAR Simulation Using the Shooting and Bouncing Ray Technique

Interferometric Synthetic Aperture Radar (IFSAR) is a technique that allows an automated way to carry out terrain mapping. IFSAR is carried out by first generating a SAR image pair from two antennas that are spatially separated. The phase difference between the SAR image pair is proportional to the topography. After registering the SAR images, the difference in phase in each pixel is extracted to generate an interferogram. Since the phase can only be measured within 2pi radians, phase unwrapping is carried out to extract the absolute phase for each pixel that will be proportional to the local height. While IFSAR algorithm is typically applied to measurement data, it is useful to develop an IFSAR simulator to develop a better understanding of the IFSAR technique. The IFSAR simulator can be used in choosing system parameters, experimenting with processing procedures and mission planning. In this paper we will present an IFSAR simulation methodology to simulate the interferogram based on the shooting and bouncing ray (SBR) technique. SBR is a standard ray-tracing technique used to simulate scattering from large, complex targets. SBR is carried out by shooting rays at the target or scene. At the exit point of each ray, a ray-tube integration is done to find its contribution to the total field. A fast algorithm has been developed for the SBR for simulating SAR images of complex targets. In the IFSAR simulation, we build upon the fast SAR simulation technique. Given the antenna pair configuration, radar system parameters and the geometrical description of the scene, we first simulate two SAR images from each antenna. After post processing the two SAR images, we generate an interferogram. Phase unwrapping is then performed on the interferogram to arrive at the desired terrain map. We will present results from the SBR-based IFSAR simulator. The results will include terrain map reconstruction of urban environments. The reconstruction will be compared to the ground truth to examine the fidelity of the simulation. We will also investigate the effect of multi-bounce scattering in urban environments on phase unwrapping and reconstruction.

Houshmand, Bijan↗

Investigation of Unsteady Flow Behavior in Transonic Compressor Rotors with LES and PIV Measurements

In the present study, unsteady flow behavior in a modern transonic axial compressor rotor is studied in detail with large eddy simulation (LES) and particle image velocimetry (PIV). The main purpose of the study is to advance the current understanding of the flow field near the blade tip in an axial transonic compressor rotor near the stall and peak-efficiency conditions. Flow interaction between the tip leakage vortex and the passage shock is inherently unsteady in a transonic compressor. Casing-mounted unsteady pressure transducers have been widely applied to investigate steady and unsteady flow behavior near the casing. Although many aspects of flow have been revealed, flow structures below the casing cannot be studied with casing-mounted pressure transducers. In the present study, unsteady velocity fields are measured with a PIV system and the measured unsteady flow fields are compared with LES simulations. The currently applied PIV measurements indicate that the flow near the tip region is not steady even at the design condition. This self-induced unsteadiness increases significantly as the compressor rotor operates near the stall condition. Measured data from PIV show that the tip clearance vortex oscillates substantially near stall. The calculated unsteady characteristics of the flow from LES agree well with the PIV measurements. Calculated unsteady flow fields show that the formation of the tip clearance vortex is intermittent and the concept of vortex breakdown from steady flow analysis does not seem to apply in the current flow field. Fluid with low momentum near the pressure side of the blade close to the leading edge periodically spills over into the adjacent blade passage. The present study indicates that stall inception is heavily dependent on unsteady behavior of the flow field near the leading edge of the blade tip section for the present transonic compressor rotor.

Hah, Chunill↗

Spatial Signal Detection Using Continuous Shrinkage Priors

Motivated by the problem of detecting changes in two-dimensional X-ray diffraction data, we propose a Bayesian spatial model for sparse signal detection in image data. Our model places considerable mass near zero and has heavy tails to reflect the prior belief that the image signal is zero for most pixels and large for an important subset. We show that the spatial prior places mass on nearby locations simultaneously being zero, and also allows for nearby locations to simultaneously be large signals. The form of the prior also facilitates efficient computing for large images. We conduct a simulation study to evaluate the properties of the proposed prior and show that it outperforms other spatial models. As a result, we apply our method in the analysis of X-ray diffraction data from a two-dimensional area detector to detect changes in the pattern when the material is exposed to an electric field.

97 MATHEMATICS AND COMPUTING↗

X-Ray Computer Methods for Studying the Structural Integrity of Seeds and Their Importance in Modern Seed Science

X-ray computer methods of research (projection microfocus radiography and microtomography), which are used to study the problem of hidden defects of seeds and investigate its impact on sowing quality, have been considered. The description and main characteristics of technical means that were used to obtain digital two-dimensional and three-dimensional (tomographic) X-ray images of seeds have been given and the possible ways of their quantitative computer processing and analysis have been discussed. Conclusions about the abilities of the methods of projection microfocus radiography and microtomography to study the features of the internal structures of a seed that are related to the violation of its integrity have been formulated.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

SU-G-IeP3-09: Investigating the Interplay of Antiscatter Grids with Modern Detectors and Image Processing in Digital Radiography

Purpose: To use phantom and simulation experiments to relate technique factors, patient size and antiscatter grid use to image quality in portable digital radiography (DR), in light of advancements in detector design and image processing. Methods: Image contrast-to-noise ratio (CNR) on a portable DR system (MobileDaRt Evolution, Shimadzu) was measured by imaging four aluminum inserts of varying thickness, superimposed on a Lucite slab phantom using a pediatric abdominal protocol. Three thicknesses of Lucite were used: 6.1cm, 12cm, and 18.2cm, with both 55 and 65 kVp beams. The mAs was set so that detector entrance exposure (DEE) was matched between kVp values. Each technique and phantom was used with and without an antiscatter grid (focused linear grid embedded in aluminum with an 8:1 ratio). The CNR-improvement-factor was then used to determine the thickness- and technique-dependent appropriateness of grid use. Finally, the same experiment was performed via Monte Carlo simulation, integrating incident energy fluence at each detector pixel, so that effects of detector design and image processing could be isolated from physical factors upstream of the detector. Results: The physical phantom experiment demonstrated a clear improvement for the lower tube voltage (55kVp), along with substantial CNR benefits with grid use for 12–18cm phantoms. Neither trend was evident with Monte Carlo, suggesting that suboptimal quantum-detection-efficiency and automated grid-removal could explain trends in kVp and grid use, respectively. Conclusion: Physical experiments demonstrate marked improvement in CNR when using a grid for phantoms of 12 and 18cm Lucite thickness (above ∼10cm soft-tissue equivalent). This benefit is likely due to image processing, as this result was not seen with Monte Carlo. The impact of image processing on image resolution should also be investigated, and the CNR benefit of low kVp and grid use should be weighed against the increased exposure time necessary to achieve adequate DEE.

60 APPLIED LIFE SCIENCES↗

Large-eddy Simulations of Spray Variability Effects on Flow Variability in a Direct-injection Spark-ignition Engine Under Non-combusting Operating Conditions

Large-eddy Simulations (LES) have been carried out to investigate spray variability and its effect on cycle-to-cycle flow variability in a direct-injection spark-ignition (DISI) engine under non-reacting conditions. Initial simulations were performed of an injector in a constant volume spray chamber to validate the simulation spray set-up. Comparisons showed good agreement in global spray measures such as the penetration. Local mixing data and shot-to-shot variability were also compared using Rayleigh-scattering images and probability contours. The simulations were found to reasonably match the local mixing data and shot-to-shot variability using a random-seed perturbation methodology. After validation, the same spray set-up with only minor changes was used to simulate the same injector in an optically accessible DISI engine. Particle Image Velocimetry (PIV) measurements were used to quantify the flow velocity in a horizontal plane intersecting the spark plug gap. The engine was operated in a skip-fired operating mode and comparisons focused on cycles that included fuel injection, but no spark event and therefore no combustion. 105 total LES engine cycles were simulated using a parallel cycle simulation approach and 3 different perturbation methods in an attempt to isolate the effects of shot-to-shot spray variability and the initial turbulent flow field as well as their interaction effects on overall engine CCVs. The experimental mean and standard deviations were reasonably well matched by the simulations, though quantitative comparisons near the injection event during the intake stroke were difficult due to the high uncertainty in the PIV measurements at these crank angles. The 3 simulation perturbation methods resulted in very similar results, though further analysis found the current parallel cycle approach may be limiting the ability of the simulations to isolate the spray and flow effects.

Cycle-to-cycle Variability↗

Intrinsic anion diffusivity in lead halide perovskites is facilitated by a soft lattice

Significance Understanding and controlling the ion migration in a solid-state semiconductor is essential for device performance in modern photonics, electronics, and optoelectronics. Direct microstructural observation and atomic-scale mechanism of the intrinsic lattice diffusion are still highly challenging. Here, we report the direct visualization of halide anion interdiffusion in CsPbCl 3 –CsPbBr 3 single crystalline nanowire heterojunctions via nanoscale imaging, allowing for the direct measuring of intrinsic lattice diffusivities. The intrinsic halide diffusivities were found to be several orders of magnitudes lower than those reported in polycrystalline thin films. Together with multiscale modeling, we further demonstrate how anionic diffusion in lead halides is determined by low barriers to vacancy hopping and facilitated by the soft perovskite lattice.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental and theoretical insights into the effects of pH on catalysis of bond-cleavage by the lignin peroxidase isozyme H8 from Phanerochaete chrysosporium

Abstract Background Lignin peroxidases catalyze a variety of reactions, resulting in cleavage of both β-O-4′ ether bonds and C–C bonds in lignin, both of which are essential for depolymerizing lignin into fragments amendable to biological or chemical upgrading to valuable products. Studies of the specificity of lignin peroxidases to catalyze these various reactions and the role reaction conditions such as pH play have been limited by the lack of assays that allow quantification of specific bond-breaking events. The subsequent theoretical understanding of the underlying mechanisms by which pH modulates the activity of lignin peroxidases remains nascent. Here, we report on combined experimental and theoretical studies of the effect of pH on the enzyme-catalyzed cleavage of β-O-4′ ether bonds and of C–C bonds by a lignin peroxidase isozyme H8 from Phanerochaete chrysosporium and an acid stabilized variant of the same enzyme. Results Using a nanostructure initiator mass spectrometry assay that provides quantification of bond breaking in a phenolic model lignin dimer we found that catalysis of degradation of the dimer to products by an acid-stabilized variant of lignin peroxidase isozyme H8 increased from 38.4% at pH 5 to 92.5% at pH 2.6. At pH 2.6, the observed product distribution resulted from 65.5% β-O-4′ ether bond cleavage, 27.0% C α -C 1 carbon bond cleavage, and 3.6% C α -oxidation as by-product. Using ab initio molecular dynamic simulations and climbing-image Nudge Elastic Band based transition state searches, we suggest the effect of lower pH is via protonation of aliphatic hydroxyl groups under which extremely acidic conditions resulted in lower energetic barriers for bond-cleavages, particularly β-O-4′ bonds. Conclusion These coupled experimental results and theoretical explanations suggest pH is a key driving force for selective and efficient lignin peroxidase isozyme H8 catalyzed depolymerization of the phenolic lignin dimer and further suggest that engineering of lignin peroxidase isozyme H8 and other enzymes involved in lignin depolymerization should include targeting stability at low pH.

09 BIOMASS FUELS↗

Dynamics and energy coupling of gas puff Z-pinches on a fast linear transformer driver

Gas puff Z-pinch experiments with annular Ar and Ne gas shells have been conducted on the Compact Experimental System for Z-pinch and Ablation Research (CESZAR) linear transformer driver (LTD) with 500 kA current and 160 ns rise time. Here, we present results from the first systematic gas puff Z-pinch experiments using a fast (≤ 200 ns) LTD as a driver, in which we show that 7% of the stored energy in the capacitors is coupled to plasma kinetic energy as estimated via self-emission and laser schlieren images. 0D and 1D simulations—which do not allow instability growth and thus reach greater maximum average velocities—using initial conditions inferred from experimental implosion trajectories predict coupling in excess of 10% of the stored energy. Additionally, the Ar and Ne implosions were comparably massed and thus achieved similar maximum kinetic energies, though the Ne pinches were more stable and the x-ray pulses were longer and produced higher yield: 2–5 ns and 0.21–0.52 J (0.15–0.37 J/cm) of Ar K-shell and 12–25 ns and 2.2–3.9 J (1.6–2.6 J/cm) of Ne K-shell, respectively. The difference in stability is most likely attributed to variations in initial conditions such as density distribution and gas breakdown initiation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Simulating Remote Sensing Systems

The Application Research Toolbox (ART) is a collection of computer programs that implement algorithms and mathematical models for simulating remote sensing systems. The ART is intended to be especially useful for performing design-tradeoff studies and statistical analyses to support the rational development of design requirements for multispectral imaging systems. Among other things, the ART affords a capability to synthesize coarser-spatial-resolution image-data products. The ART also provides for simulations of image-degradation effects, including point-spread functions, misregistration of spectral images, and noise. The ART can utilize real or synthetic data sets, along with sensor specifications, to create simulated data sets. In one example of a particular application, simulated imagery of a coarse resolution system was created using high-resolution imagery from another system in order to perform a radiometric cross-comparison. In the case of a proposed sensor system, the simulated data can be used to conduct trade studies and statistical analyses to ensure that the sensor system will satisfy the requirements of potential scientific, academic, and commercial user communities.

Zanoni, Vicki↗

Design of teleoperation system with a force-reflecting real-time simulator

We developed a force-reflecting teleoperation system that uses a real-time graphic simulator. This system eliminates the effects of communication time delays in remote robot manipulation. The simulator provides the operator with predictive display and feedback of computed contact forces through a six-degree of freedom (6-DOF) master arm on a real-time basis. With this system, peg-in-hole tasks involving round-trip communication time delays of up to a few seconds were performed at three support levels: a real image alone, a predictive display with a real image, and a real-time graphic simulator with computed-contact-force reflection and a predictive display. The experimental results indicate the best teleoperation efficiency was achieved by using the force-reflecting simulator with two images. The shortest work time, lowest sensor maximum, and a 100 percent success rate were obtained. These results demonstrate the effectiveness of simulated-force-reflecting teleoperation efficiency.

Hirata, Mitsunori↗

Benefits of adaptive radiation therapy in lung cancer as a function of replanning frequency

Purpose: To quantify the potential benefit associated with daily replanning in lung cancer in terms of normal tissue dose sparing and to characterize the tradeoff between adaptive benefit and replanning frequency. Methods: A set of synthetic images and contours, derived from weekly active breathing control images of 12 patients who underwent radiation therapy treatment for nonsmall cell lung cancer, is generated for each fraction of treatment using principal component analysis in a way that preserves temporal anatomical trends (e.g., tumor regression). Daily synthetic images and contours are used to simulate four different treatment scenarios: (1) a “no-adapt” scenario that simulates delivery of an initial plan throughout treatment, (2) a “midadapt” scenario that implements a single replan for fraction 18, (3) a “weekly adapt” scenario that simulates weekly adaptations, and (4) a “full-adapt” scenario that simulates daily replanning. An initial intensity modulated radiation therapy plan is created for each patient and replanning is carried out in an automated fashion by reoptimizing beam apertures and weights. Dose is calculated on each image and accumulated to the first in the series using deformable mappings utilized in synthetic image creation for comparison between simulated treatments. Results: Target coverage was maintained and cord tolerance was not exceeded for any of the adaptive simulations. Average reductions in mean lung dose (MLD) and volume of lung receiving 20 Gy or more (V20{sub lung}) were 65 ± 49 cGy (p = 0.000 01) and 1.1% ± 1.2% (p = 0.0006), respectively, for all patients. The largest reduction in MLD for a single patient was 162 cGy, which allowed an isotoxic escalation of the target dose of 1668 cGy. Average reductions in cord max dose, mean esophageal dose (MED), dose received by 66% of the heart (D66{sub heart}), and dose received by 33% of the heart (D33{sub heart}), were 158 ± 280, 117 ± 121, 37 ± 77, and 99 ± 120 cGy, respectively. Average incremental reductions in MLD for the midadapt, weekly adapt, and full-adapt treatments were 38, 18, and 8 cGy, respectively. Incremental reductions in MED for the same treatments were 57, 37, and 23 cGy. Reductions in MLD and MED for the full-adapt treatment were correlated with the absolute decrease in the planning target volume (r = 0.34 and r = 0.26). Conclusions: Adaptive radiation therapy for lung cancer yields clinically relevant reductions in normal tissue doses for frequencies of adaptation ranging from a single replan up to daily replanning. Increased frequencies of adaptation result in additional benefit while magnitude of benefit decreases.

60 APPLIED LIFE SCIENCES↗

SU-G-TeP4-04: An Automated Monte Carlo Based QA Framework for Pencil Beam Scanning Treatments

Purpose: Prior to treating new PBS field, multiple (three) patient-field-specific QA measurements are performed: two 2D dose distributions at shallow depth (M1) and at the tumor depth (M2) with treatment hardware at zero gantry angle; one 2D dose distribution at iso-center (M3) without patient specific devices at the planned gantry angle. This patient-specific QA could be simplified by the use of MC model. The results of MC model commissioning for a spot-scanning system and the fully automated TOPAS/MC-based QA framework will be presented. Methods: We have developed in-house MC interface to access a TPS (Astroid) database from a computer cluster remotely. Once a plan is identified, the interface downloads information for the MC simulations, such as patient images, apertures points, and fluence maps and initiates calculations in both the patient and QA geometries. The resulting calculations are further analyzed to evaluate the TPS dose accuracy and the PBS delivery. Results: The Monte Carlo model of our system was validated within 2.0 % accuracy over the whole range of the dose distribution (proximal/shallow part, as well as target dose part) due to the location of the measurements. The averaged range difference after commissioning was 0.25 mm over entire treatment ranges, e.g., 6.5 cm to 31.6 cm. Conclusion: As M1 depths range typically from 1 cm to 4 cm from the phantom surface, The Monte Carlo model of our system was validated within +− 2.0 % in absolute dose level over a whole treatment range. The averaged range difference after commissioning was 0.25 mm over entire treatment ranges, e.g., 6.5 cm to 31.6 cm. This work was supported by NIH/NCI under CA U19 21239.

60 APPLIED LIFE SCIENCES↗

A finite mixtures algorithm for finding classes in images

The authors describe an algorithm for estimating the proportions of classes in a SAR (synthetic aperture radar) image without any user interaction. The method assumes that the image is a mixture of a known number of different pixel types. A maximum likelihood estimate of the parameters of the resulting mixture distribution is then used to find the proportions for the various classes. The technique was successfully applied to aircraft SAR images of sea ice. Computer simulations were used to determine the relative errors of the technique. The technique performs well even with extremely noisy images.

Samadani, Ramin↗

Evaluation of Techniques for Vegetation Removal from UAV-Based Photogrammetric DSM near WIPP Land Withdrawal Act Boundary, NM

The Waste Isolation Pilot Plant is located within New Mexico's karst landscape and is the United States' only deep geologic radioactive waste repository which isolates transuranic waste from defense activities underground in a bedded salt formation. Characteristic karst topography in this region is made up of material (such as limestone) that can be dissolved by water over a long period of time, and includes features such as springs, surface streams, sinking streams, caves, and sinkholes. Karst features can potentially impact the integrity of the waste repository in the future and compromise WIPP's performance, particularly when coupled with incompatible land-use activities within the Land Withdrawal Act (LWA) boundary as well as increased water withdrawals outside the LWA. In order to better capture the ground surface topography to support hydrology modeling efforts, the U.S. DOE has enlisted FIU to develop a high-resolution digital elevation model (DEM). In a previous field study, a high-resolution digital surface model (DSM) of an adjacent representative site, Basin 6, was developed by photogrammetry using aerial images captured with unmanned aerial vehicles. The DSM generated was based on point clouds representing the entire landscapes, including points of terrain, vegetation and infrastructure. These above-ground features need to be removed to create a DSM/DEM with good representation of the bare ground. In recent years, many techniques including both pre- and post-processing of DSMs, were developed for vegetation removal with varied success. This research focuses on applying and comparing various vegetation removal methods using point clouds generated of the Basin 6 pilot study area. The 3D RGB-based point cloud classification using Python scripts was a modification of the procedure by Themistocleous (2019). The method using Liblas classifies the point cloud three-dimensionally, whereas the method by Themistocleous (2019) is two-dimensional. The preliminary result of the method using Liblas, however, still had a large amount of vegetation remaining compared to the Pix4D Machine Learning method. The vegetation layer remains highly visible in the DEMs processed using the vegetation indices as opposed to the Pix4D Machine Learning method, which produces a DEM where mostly bare ground is seen. Further modifications to the input parameters are required to improve the vegetation vs. bare ground classification using this approach. The next step will be to improve the method using Liblas and to also test other vegetation removal methods found in the literature review to compare them against the two described here.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Four-Dimensional Dose Reconstruction for Scanned Proton Therapy Using Liver 4DCT-MRI

Purpose: Four-dimensional computed tomography-magnetic resonance imaging (4DCT-MRI) is an image-processing technique for simulating many 4DCT data sets from a static reference CT and motions extracted from 4DMRI studies performed using either volunteers or patients. In this work, different motion extraction approaches were tested using 6 liver cases, and a detailed comparison between 4DCT-MRI and 4DCT was performed. Methods and Materials: 4DCT-MRI has been generated using 2 approaches. The first approach used motion extracted from 4DMRI as being “most similar” to that of 4DCT from the same patient (subject-specific), and the second approach used the most similar motion obtained from a motion library derived from 4DMRI liver studies of 13 healthy volunteers (population-based). The resulting 4DCT-MRI and 4DCTs were compared using scanned proton 4D dose calculations (4DDC). Results: Dosimetric analysis showed that 93% ± 8% of points inside the clinical target volume (CTV) agreed between 4DCT and subject-specific 4DCT-MRI (gamma analysis: 3%/3 mm). The population-based approach however showed lower dosimetric agreement with only 79% ± 14% points in the CTV reaching the 3%/3 mm criteria. Conclusions: 4D CT-MRI extends the capabilities of motion modeling for dose calculations by accounting for realistic and variable motion patterns, which can be directly employed in clinical research studies. We have found that the subject-specific liver modeling appears more accurate than the population-based approach. The former is particularly interesting for clinical applications, such as improved target delineation and 4D dose reconstruction for patient-specific QA to allow for inter- and/or intra-fractional plan corrections.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Imaging surface topography with coherent x-ray reflectivity: Theory, kinematics, and simulations

A theoretical formalism is described for understanding coherent x-ray reflectivity (CXR) from the surface of a semi-infinite crystal having a variable surface topography, described by the height profile ℎ(𝑥,𝑦). The surface topography is imaged as a complex “effective density,” obtained from the phasing and inversion of the coherent x-ray reflectivity data, measured through a rocking scan centered at a vertical momentum transfer 𝑄$^{0}_{𝑧}$ and a vertical range Δ⁢𝑄 𝑧 . The formalism predicts that the effective density has an amplitude with a maximum located at the surface height for each position within the surface plane. The phase of the effective density has a lateral variation that is controlled by the surface height and a vertical variation that reflects a combination of the interfacial structure and specific choice of measurement conditions. This understanding enables direct observation of nanometer-scale interfacial topography, i.e., ℎ⁡(𝑥,𝑦)⁢𝑐 𝑠 (where 𝑐 𝑠 is the vertical substrate lattice parameter) with Å-scale sensitivity to surface height. Numerical simulations illustrate and confirm the theoretical results. These results show how the interpretation of the interfacial density phase obtained by CXR data inversion (i.e., surface topography with respect to a flat surface) is conceptually similar to that previously known for Bragg coherent diffraction imaging (BCDI) measurements of isolated nanoparticles (i.e., lattice displacements with respect to an ideal crystal lattice). This suggests that CXR can be thought of as a form of dark field imaging with respect to the bright field BCDI approach. An implication of these results is that interfacial imaging may bypass some of the significant challenges associated with BCDI imaging of multiple particles having different orientations.

X-ray imaging↗