Search NASA⌕ Search

SEARCH · Search NASA

Results for “matched field processing”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Sea surface topography fields of the tropical Pacific from data assimilation

Time series of maps of monthly tropical Pacific dynamic topography anomalies from 1979 through 1985 were constructed by means of assimilation of the tide gauge and expendable bathythermograph (XBT) data into a linear model driven by observed winds. Estimates of error statistics were calculated and compared to actual differences between hindcasts and observations. Four experiements were performed as follows: one with no assimilation, one with assimiation of sea level anomaly data from eight selected island tide gauge stations, one with assimilation of dynamic height anomalies derived from XBT data, and one with both XBT and tide gauge data assimilated. Data from seven additional tide gauge stations were withheld from the assimilation process and used for verification in all four experiements. Statistical objective maps based on data alone were also constructed for comparison purposes. The dynamic response of the model without assimilation was, in general, weaker than the observed response. Assimilation resulted in enhanced signal amplitude in all three assimilation experiments. With few exceptions the error estimated generated by the Kalman filter appeared quite reasonable. Since the error processes cannot be assumed to be white or stationary, we could find no straightforward way to test the formal statistical hypothesis that the time series of differences between the filter ouput and the actual observations were drawn from a population with statistics given by the Kalman filter estimaes. We found that the best filter we could devise was still short of the goal of producing a white innovation sequence. Extensive changes in the assumed error statistics make only marginal differences. The same is true for long time and space scale behavior of different models with richer physics and finer resolution. Better data assimilation results will probably require relaxation of the assumptions of stationarity and serial independence of the errors. Formation of such detailed noise models will require longer time series, with the attendant problems of matching very different data sets.

Miller, Robert N.↗
NASA NTRS2013

Support Routines for In Situ Image Processing

This software consists of a set of application programs that support ground-based image processing for in situ missions. These programs represent a collection of utility routines that perform miscellaneous functions in the context of the ground data system. Each one fulfills some specific need as determined via operational experience. The most unique aspect to these programs is that they are integrated into the large, in situ image processing system via the PIG (Planetary Image Geometry) library. They work directly with space in situ data, understanding the appropriate image meta-data fields and updating them properly. The programs themselves are completely multimission; all mission dependencies are handled by PIG. This suite of programs consists of: (1)marscahv: Generates a linearized, epi-polar aligned image given a stereo pair of images. These images are optimized for 1-D stereo correlations, (2) marscheckcm: Compares the camera model in an image label with one derived via kinematics modeling on the ground, (3) marschkovl: Checks the overlaps between a list of images in order to determine which might be stereo pairs. This is useful for non-traditional stereo images like long-baseline or those from an articulating arm camera, (4) marscoordtrans: Translates mosaic coordinates from one form into another, (5) marsdispcompare: Checks a Left Right stereo disparity image against a Right Left disparity image to ensure they are consistent with each other, (6) marsdispwarp: Takes one image of a stereo pair and warps it through a disparity map to create a synthetic opposite- eye image. For example, a right eye image could be transformed to look like it was taken from the left eye via this program, (7) marsfidfinder: Finds fiducial markers in an image by projecting their approximate location and then using correlation to locate the markers to subpixel accuracy. These fiducial markets are small targets attached to the spacecraft surface. This helps verify, or improve, the pointing of in situ cameras, (8) marsinvrange: Inverse of marsrange . given a range file, re-computes an XYZ file that closely matches the original. . marsproj: Projects an XYZ coordinate through the camera model, and reports the line/sample coordinates of the point in the image, (9) marsprojfid: Given the output of marsfidfinder, projects the XYZ locations and compares them to the found locations, creating a report showing the fiducial errors in each image. marsrad: Radiometrically corrects an image, (10) marsrelabel: Updates coordinate system or camera model labels in an image, (11) marstiexyz: Given a stereo pair, allows the user to interactively pick a point in each image and reports the XYZ value corresponding to that pair of locations. marsunmosaic: Extracts a single frame from a mosaic, which will be created such that it could have been an input to the original mosaic. Useful for creating simulated input frames using different camera models than the original mosaic used, and (12) merinverter: Uses an inverse lookup table to convert 8-bit telemetered data to its 12-bit original form. Can be used in other missions despite the name.

Deen, Robert G.↗
NASA NTRS2024

Oculometric Biomarkers of Visuomotor Deficits in Clinically Asymptomatic Patients With Systemic Lupus Erythematosus Undergoing Long-Term Hydroxychloroquine Treatment

Introduction: This study examines a set of oculomotor measurements, or “oculometric” biomarkers, as potential early indicators of visual and visuomotor deficits due to retinal toxicity in asymptomatic Systemic Lupus Erythematosus (SLE) patients on long-term hydroxychloroquine (HCQ) treatment. The aim is to identify subclinical functional impairments that are otherwise undetectable by standard clinical tests and to link them to structural retinal changes. Methods: We measured oculomotor responses in a cohort of SLE patients on chronic HCQ therapy using a previously established behavioral task and analysis technique. We also examined the relationship between oculometrics, OCT measures of retinal thickness, and standard clinical perimetry measures of visual function in our patient group using Bivariate Pearson Correlation and a Linear Mixed-Effects Model (LMM). Results: Significant visual and visuomotor deficits were found in 12 asymptomatic SLE patients on long-term HCQ therapy compared to a cohort of 17 age-matched healthy controls. Notably, six oculometrics were significantly different. The median initial pursuit acceleration was 22%, steady-state pursuit gain 16%, proportion smooth 7%, and target speed responsiveness 31% lower, while catch-up saccade amplitude was 46% and fixation error 46% larger. Excluding the two patients with diagnosed mild toxicity, four oculometrics, all but fixation error and proportion smooth, remained significantly impaired compared to controls. Across our population of 12 patients (24 retinae), we found that pursuit latency, initial acceleration, steady-state gain, and fixation error were linearly related to retinal thickness even when age was accounted for, while standard measures of clinical function (Mean Deviation and Pattern Standard Deviation) were not. Discussion: Our data show that specific oculometrics are sensitive early biomarkers of functional deficits in SLE patients on HCQ that could be harnessed to assist in the early detection of HCQ-induced retinal toxicity and other visual pathologies, potentially providing early diagnostic value beyond standard visual field and OCT evaluations.

eye movements↗
DOE OSTI2016

SU-F-I-31: Reproducibility of An Automatic Exposure Control Technique in the Low-Dose CT Scan of Cardiac PET/CT Exams

Purpose: Two separate low-dose CT scans are usually performed for attenuation correction of rest and stress N-13 ammonia PET/CT myocardial perfusion imaging (PET/CT). We utilize an automatic exposure control (AEC) technique to reduce CT radiation dose while maintaining perfusion image quality. Our goal is to assess the reproducibility of displayed CT dose index (CTDI) on same-day repeat CT scans (CT1 and CT2). Methods: Retrospectively, we reviewed CT images of PET/CT studies performed on the same day. Low-dose CT utilized AEC technique based on tube current modulation called Smart-mA. The scan parameters were 64 × 0.625mm collimation, 5mm slice thickness, 0.984 pitch, 1-sec rotation time, 120 kVp, and noise index 50 with a range of 10–200 mA. The scan length matched with PET field of view (FOV) with the heart near the middle of axial FOV. We identified the reference slice number (RS) for an anatomical landmark (carina) and used it to estimate axial shift between two CTs. For patient size, we measured an effective diameter on the reference slice. The effect of patient positioning to CTDI was evaluated using the table height. We calculated the absolute percent difference of the CTDI (%diff) for estimation of the reproducibility. Results: The study included 168 adults with an average body-mass index of 31.72 ± 9.10 (kg/m{sup 2}) and effective diameter was 32.72 ± 4.60 cm. The average CTDI was 1.95 ± 1.40 mGy for CT1 and 1.97 ± 1.42mGy for CT2. The mean %diff was 7.8 ± 6.8%. Linear regression analysis showed a significant correlation between the table height and %diff CTDI. (r=0.82, p<0.001) Conclusion: We have shown for the first time in human subjects, using two same-day CT images, that the AEC technique in low-dose CT is reproducible within 10% and significantly depends on the patient centering.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2011

Preliminary Study on Two Possible Bunch Compression Schemes at NLCTA

In this paper, two possible bunch compression configurations are proposed and evaluated by numerical simulation in the Next Linear Collider Test Accelerator (NLCTA) at SLAC. A bunch compression ratio up to 20 could be achieved under a perfect condition, without consideration for the timing jitter and other error sources. The NLCTA is a test accelerator built at SLAC, which is approximately 42 meters long and composed of X-band acceleration structures. The main aim of building NLCTA is to develop and demonstrate the X-band rf acceleration technologies for the next generation linear collider, with a relatively high acceleration gradient between 50 MV/m and 100 MV/m. The current operation configuration of NLCTA features a thermionic-cathode electron gun at its starting point which generates an electron beam with an energy of 5 MeV. This is followed by a roughly 1.5 meter long X-band acceleration structure which boosts the electron beam energy to 60 MeV. Then there is a four-dipole magnetic chicane which is 6 meters long and provides a first order longitudinal dispersion of R{sub 56} = -73mm. Next the electron beam passes by several matching quadrupoles and can be accelerated further to 120 MeV through another one-meter-long X-band acceleration structure. After that, there are three small chicanes downstream, with a total first order longitudinal dispersion of R{sub 56} = -10mm. A sketch of the main components of NLCTA is shown in Figure 1, where the total length of this accelerator is 45 meters. Free Electron Lasers (FELs), proposed by J. Madey and demonstrated for the first time at Stanford University in 1970s [2] [3], use the lasing of relativistic electron beam traveling through a magnetic undulator, which can reach high power and can be widely tunable in wavelength. Linac based FEL source can provide sufficient brightness, and a short X-ray wavelength down to angstrom scale, which promises in supporting wide range of research experiments. In order to have an electron beam lasing coherently in an undulator, one needs a very bright beam in all three dimensions. In other words, one needs an electron beam with very short bunch length (high intensity), very small transverse emittance and very small energy spread. Most FELs currently being operated, commissioned, constructed or proposed are based on RF acceleration in a frequency range from L-band ( 1 GHz) to C-band ( 6 GHz). As RF frequency goes higher, wake fields effects tend to be much stronger and jitter tolerances are tighter. To demonstrate that X-band acceleration structures can be applied in constructing an FEL, one could perform bunch compression experiments at NLCTA as a first step, and investigate tolerances on timing jitter, misalignments etc.. Another important point is to evaluate the transverse emittance growth in this bunch compression process. In the following sections, two possible bunch compression schemes are proposed to be tested at NLCTA. Elegant [4] 3-D simulation is performed to evaluate these two schemes, with wake fields, space charge and coherent synchrotron radiation (CSR) effects included. One million macro particles are adopted in the numerical simulations. The simulation starts with an electron beam of 20 pC at a beam energy of 5 MeV. The initial RMS bunch length is taken as 0.5 ps at such a low bunch charge, and the RMS energy spread is 5 x 10{sup -3}. The normalized transverse emittance is 1 mm.mrad.

43 PARTICLE ACCELERATORS↗
NASA NTRS2004

Winds, Water Budgets and Stable Isotopes in Tropical Cyclones using TRMM and QUICKSCAT

Water vapor is the most abundant greenhouse gas in the atmosphere. Changes in its concentration and distribution are controlled by the hydrologic cycle. Because of its capacity to absorb and emit long wave radiation, release latent heat during condensation in storms and reflect short wave radiation when clouds form it has a major impact on Global climate change. The stable isotope ratios of water are H20 H2l6O and H0 H2l6O. These ratios change whenever water undergoes a phase change. They also change in both rain and water vapor whenever an air parcel is exposed to rain. In addition the relative changes in the two ratios differ as a &nction of the relative humidity. In short, the stable isotope ratios in water vapor in the atmosphere contain an integrated history of the processes affecting the concentration and distribution of water vapor in the atmosphere. Therefore the measurement and interpretation of changes in these stable isotope ratios are a powerful tool matched by no other method in tracing the transport history of water in the atmosphere. Our initial studies under this grant focused on the changes of the stable isotope ratios of precipitation and water vapor in tropical cyclones. The changes in time and space were found to be very large and to trace the transport of water in the storms reflecting changes in basic structural features. Because the stable isotope ratios of rains from tropical cyclones are so low flooding associated with land falling tropical cyclones introduces a negative isotopic spike into the coastal surface waters. In addition the stable isotope ratios of water vapor in the vicinity of tropical cyclones is anomalously low. This suggests that carbonate shelled organisms such as ostracoda living in coastal waters have the potential to record the isotopic spike and thereby provide a long term record of tropical storm activity in sediment cores containing fossil shells. Likewise, tree rings in coastal environments offer a similar potential. We have analyzed the oxygen isotopic composition of ostrcoda shells formed in the floodwaters of Tropical Storm Allison (2001) and discovered the negative isotopic 1 16 spike. Because we had learned that storm activity has a major impact on the stable isotope ratios of water vapor in the tropics and sub-tropics we decided to analyze the isotopic compositions of water vapor in different locations in the tropics. We did this in Puerto Escondido, Mexico in July 1998, near Kwajalein Island in the Pacific in 1999 as part of a TRMM summer field program and in 2001 in Key West, Florida as part of the CAMEX 4 summer field program. Our isotopic studies along with our earlier tropical cyclone studies showed that the low isotopic ratios in water vapor induced by exposure to rains the storms persisted for 48 hours often far away from the original storm site. We also noted that positive isotopic spikes were introduced into atmospheric water vapor if winds were high and extensive sea spray was present. These findings have a significant impact on the interpretation of the stable isotope studies of tropical ice cores found in the high mountain regions of the tropics. The assumption made in interpreting the ice core record is that the source water vapor evaporated from the sea surface is in near isotopic equilibrium with the seawater and undergoes a decrease during its transport that reflects the change in temperature from the sea surface to the site of the ice core. Because an additional isotopic depletion occurs at the sea surface source area that depends on the intensity, duration and size of the tropical rain system the isotopic variations found in the ice cores must take into account changes in past storm activity in the tropics. These systems must be an important source of water vapor to the ice cores because they charge the troposphere with water vapor to a far greater vertical height than evaporation in quiescent regions. Finally, an interest in increased heat transfer in thnterior of tropical cyclones resulting from greater amounts of sea spray is a topic of considerable interest to the research community. Increases in sea spray may be related to rapid changes in the intensity of hurricanes, a feature of hurricanes that currently is very poorly forecast. Project CBLAST of the Hurricane Research Division of NOAA is an active program that uses P3 research aircraft to evaluate this problem. An instrument has been designed and built at the University of Houston that will be placed on the P3 research aircraft during the 2004 hurricane season. It continuously measures the salt content of rain in hurricanes. Changes in the salt content of the rains should reflect the abundance of sea spray at the sea surface. In this way maps of sea spray intensity in hurricanes may be forthcoming. This should help computer modelers who simulate hurricanes to better understand the potential of changes in sea spray to change the intensity of hurricanes. The cost of designing, building and installing this instrument was borne largely by funding from this NASA grant. A list of presentations at national meetings and publications that were as the result of funding from this NASA grant are found in the report.

Lawrence, James R.↗
DOE OSTI2021

Factorization at subleading power and endpoint divergences in h → γγ decay. Part II. Renormalization and scale evolution

Building on the recent derivation of a bare factorization theorem for the b-quark induced contribution to the h → γγ decay amplitude based on soft-collinear effective theory, we derive the first renormalized factorization theorem for a process described at subleading power in scale ratios, where λ = m b /M h « 1 in our case. We prove two refactorization conditions for a matching coefficient and an operator matrix element in the endpoint region, where they exhibit singularities giving rise to divergent convolution integrals. The refactorization conditions ensure that the dependence of the decay amplitude on the rapidity regulator, which regularizes the endpoint singularities, cancels out to all orders of perturbation theory. We establish the renormalized form of the factorization formula, proving that extra contributions arising from the fact that “endpoint regularization” does not commute with renormalization can be absorbed, to all orders, by a redefinition of one of the matching coefficients. We derive the renormalization-group evolution equation satisfied by all quantities in the factorization formula and use them to predict the large logarithms of order $αα^2_sL^k$ in the three-loop decay amplitude, where $L=ln(-M^2_h/m^2_b)$ and k = 6, 5, 4, 3. We find perfect agreement with existing numerical results for the amplitude and analytical results for the three-loop contributions involving a massless quark loop. On the other hand, we disagree with the results of previous attempts to predict the series of subleading logarithms $~αα^n_sL^{2n+1}$

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2011

Multi-View Budgeted Learning under Label and Feature Constraints Using Label-Guided Graph-Based Regularization

Budgeted learning under constraints on both the amount of labeled information and the availability of features at test time pertains to a large number of real world problems. Ideas from multi-view learning, semi-supervised learning, and even active learning have applicability, but a common framework whose assumptions fit these problem spaces is non-trivial to construct. We leverage ideas from these fields based on graph regularizers to construct a robust framework for learning from labeled and unlabeled samples in multiple views that are non-independent and include features that are inaccessible at the time the model would need to be applied. We describe examples of applications that fit this scenario, and we provide experimental results to demonstrate the effectiveness of knowledge carryover from training-only views. As learning algorithms are applied to more complex applications, relevant information can be found in a wider variety of forms, and the relationships between these information sources are often quite complex. The assumptions that underlie most learning algorithms do not readily or realistically permit the incorporation of many of the data sources that are available, despite an implicit understanding that useful information exists in these sources. When multiple information sources are available, they are often partially redundant, highly interdependent, and contain noise as well as other information that is irrelevant to the problem under study. In this paper, we are focused on a framework whose assumptions match this reality, as well as the reality that labeled information is usually sparse. Most significantly, we are interested in a framework that can also leverage information in scenarios where many features that would be useful for learning a model are not available when the resulting model will be applied. As with constraints on labels, there are many practical limitations on the acquisition of potentially useful features. A key difference in the case of feature acquisition is that the same constraints often don't pertain to the training samples. This difference provides an opportunity to allow features that are impractical in an applied setting to nevertheless add value during the model-building process. Unfortunately, there are few machine learning frameworks built on assumptions that allow effective utilization of features that are only available at training time. In this paper we formulate a knowledge carryover framework for the budgeted learning scenario with constraints on features and labels. The approach is based on multi-view and semi-supervised learning methods that use graph-encoded regularization. Our main contributions are the following: (1) we propose and provide justification for a methodology for ensuring that changes in the graph regularizer using alternate views are performed in a manner that is target-concept specific, allowing value to be obtained from noisy views; and (2) we demonstrate how this general set-up can be used to effectively improve models by leveraging features unavailable at test time. The rest of the paper is structured as follows. In Section 2, we outline real-world problems to motivate the approach and describe relevant prior work. Section 3 describes the graph construction process and the learning methodologies that are employed. Section 4 provides preliminary discussion regarding theoretical motivation for the method. In Section 5, effectiveness of the approach is demonstrated in a series of experiments employing modified versions of two well-known semi-supervised learning algorithms. Section 6 concludes the paper.

99 GENERAL AND MISCELLANEOUS↗
NASA NTRSDate not supplied

Stability Analysis of the Flow Over a Swept Forward-Facing Step Using PIV Base Flows

Step excrescences are a type of surface imperfection encountered on swept wings of commercial aircraft, commonly because accessibility requirements prevent creating the wing’s surface from a single panel. If the step height is too large, super-critical, the flow can undergo an early transition to turbulence, which can be highly detrimental for the performance of the wing. In the case of a forward-facing step on a swept wing in a low-disturbance environment, stationary crossflow vortices can develop a significant amplitude upstream of the step and hence dominate the structure of the boundary-layer flow over the step. The main goal of the present investigation is to illuminate the path to transition supported by the fascinatingly complex flow field in the direct downstream vicinity of a step with a super critical height. The high-resolution, stereographic Particle Image Velocimetry (PIV) measurement dataset presently available for this flow field provides a complete description of the laminar flow for the execution of BiGlobal stability analysis in a plane parallel to the step. Although the notorious sensitivity of stability results to the description of the base flow demands a very careful uncertainty analysis of those results, it is argued that this very fact can be leveraged to produce new insight into the supported perturbation dynamics. In performing the analysis, several unsteady mode families are discovered that display the explosive perturbation expected for early transition to be induced. In considering domain widths equal to an integer-multiple of the incident crossflow-vortex wavelength and analyzing an extent of 5 crossflow-vortex wavelengths parallel to the step, it is found that the stability results converge while increasing the domain width. It is demonstrated, moreover, that the results for the wider domains can be approximated by appropriately averaging the results on neighboring single-crossflow-vortex-wavelength domains covering the same region. Besides being useful for computational purposes, this observed property suggests interpreting the instability mechanism as a distorted primary mechanism rather than a “proper” secondary mechanism. This follows in the context of the secondary in-stability analysis of three-dimensional boundary layers, because the secondary mechanism is usually characterized by being localized in a pocket of strong shear, while the distorted primary mechanism typically has an infinite support in the direction parallel to the step. Even though the growth rates are found to be sensitive to the interrogation-window size inherent to the PIV post-processing procedure, the spatial structure of the eigenfunctions is found to be relatively insensitive. Lastly, the spatial structure of the eigenfunctions corresponding to all velocity components are matched with the shape functions determined by computing the Spectral Proper Orthogonal Decomposition (SPOD) of a time-resolved measurement of the perturbation content.

forward-facing step↗
DOE OSTI2015

The morphologies of massive galaxies from z ∼ 3—witnessing the two channels of bulge growth

We quantify the morphological evolution of z∼0 massive galaxies (M{sub ∗}/M{sub ⊙}∼10{sup 11.2±0.3}) from z∼3 in the 5 CANDELS fields. The progenitors are selected using abundance matching techniques to account for the mass growth. At z<1, the population matches the massive end of the Hubble sequence, with 30% of pure spheroids, 50% of galaxies with equally dominant disk and bulge components, and 20% of disks. At z∼2−3 however, there is a majority of irregular systems (∼60%−70%) with still 30% of pure spheroids. We then analyze the stellar populations, star formation rates (SFRs), gas fractions, and structural properties for the different morphologies independently. Our results suggest two distinct channels for the growth of bulges in massive galaxies. Around ∼30%−40% were already bulges at z∼2.5, with low average SFRs and gas fractions (10%−15%), high Sérsic indices (n>3−4), and small effective radii (R{sub e}∼1 kpc), pointing toward an even earlier formation through gas-rich mergers or violent disk instabilities. Between z∼2.5 and z∼0, they rapidly increase their size by a factor of ∼4−5, are quenched, and slightly increase their Sérsic indices (n∼5) but their global morphology remains unaltered. The structural evolution is independent of the gas fractions, suggesting that it is driven by ex situ events. The remaining 60% experience a gradual morphological transformation, from clumpy disks to more regular bulge+disk systems, essentially happening at z>1. This results in the growth of a significant bulge component (n∼3) for 2/3 of the systems, possibly through the migration of clumps, while the remaining 1/3 retain a rather small bulge (n∼1.5−2). The transition phase between disturbed and relaxed systems and the emergence of the bulge is correlated with a decrease in the star formation activity and the gas fractions, suggesting a morphological quenching process as a plausible mechanism for the formation of these bulges.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

WE-AB-204-09: Respiratory Motion Correction in 4D-PET by Simultaneous Motion Estimation and Image Reconstruction (SMEIR)

Purpose: In conventional 4D-PET, images from different frames are reconstructed individually and aligned by registration methods. Two issues with these approaches are: 1) Reconstruction algorithms do not make full use of all projections statistics; and 2) Image registration between noisy images can Result in poor alignment. In this study we investigated the use of simultaneous motion estimation and image reconstruction (SMEIR) method for cone beam CT for motion estimation/correction in 4D-PET. Methods: Modified ordered-subset expectation maximization algorithm coupled with total variation minimization (OSEM- TV) is used to obtain a primary motion-compensated PET (pmc-PET) from all projection data using Demons derived deformation vector fields (DVFs) as initial. Motion model update is done to obtain an optimal set of DVFs between the pmc-PET and other phases by matching the forward projection of the deformed pmc-PET and measured projections of other phases. Using updated DVFs, OSEM- TV image reconstruction is repeated and new DVFs are estimated based on updated images. 4D XCAT phantom with typical FDG biodistribution and a 10mm diameter tumor was used to evaluate the performance of the SMEIR algorithm. Results: Image quality of 4D-PET is greatly improved by the SMEIR algorithm. When all projections are used to reconstruct a 3D-PET, motion blurring artifacts are present, leading to a more than 5 times overestimation of the tumor size and 54% tumor to lung contrast ratio underestimation. This error reduced to 37% and 20% for post reconstruction registration methods and SMEIR respectively. Conclusion: SMEIR method can be used for motion estimation/correction in 4D-PET. The statistics is greatly improved since all projection data are combined together to update the image. The performance of the SMEIR algorithm for 4D-PET is sensitive to smoothness control parameters in the DVF estimation step.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

RESOLVING THE FAR-IR LINE DEFICIT: PHOTOELECTRIC HEATING AND FAR-IR LINE COOLING IN NGC 1097 AND NGC 4559

The physical state of interstellar gas and dust is dependent on the processes which heat and cool this medium. To probe heating and cooling of the interstellar medium over a large range of infrared surface brightness, on sub-kiloparsec scales, we employ line maps of [C II] 158 {mu}m, [O I] 63 {mu}m, and [N II] 122 {mu}m in NGC 1097 and NGC 4559, obtained with the Photodetector Array Camera and Spectrometer on board Herschel. We matched new observations to existing Spitzer Infrared Spectrograph data that trace the total emission of polycyclic aromatic hydrocarbons (PAHs). We confirm at small scales in these galaxies that the canonical measure of photoelectric heating efficiency, ([C II] + [O I])/TIR, decreases as the far-infrared (far-IR) color, {nu}f{sub {nu}}(70 {mu}m) {nu}f{sub {nu}}(100 {mu}m), increases. In contrast, the ratio of far-IR cooling to total PAH emission, ([C II] + [O I])/PAH, is a near constant {approx}6% over a wide range of far-IR color, 0.5 < {nu}f{sub {nu}}(70 {mu}m) {nu}f{sub {nu}}(100 {mu}m) {approx}< 0.95. In the warmest regions, where {nu}f{sub {nu}}(70 {mu}m) {nu}f{sub {nu}}(100 {mu}m) {approx}> 0.95, the ratio ([C II] + [O I])/PAH drops rapidly to 4%. We derived representative values of the local ultraviolet radiation density, G{sub 0}, and the gas density, n{sub H}, by comparing our observations to models of photodissociation regions. The ratio G{sub 0}/n{sub H}, derived from fine-structure lines, is found to correlate with the mean dust-weighted starlight intensity, (U), derived from models of the IR spectral energy distribution. Emission from regions that exhibit a line deficit is characterized by an intense radiation field, indicating that small grains are susceptible to ionization effects. We note that there is a shift in the 7.7/11.3 {mu}m PAH ratio in regions that exhibit a deficit in ([C II] + [O I])/PAH, suggesting that small grains are ionized in these environments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

V‐Groove Si Nanopatterning for the Direct Epitaxy of Orientation‐Patterned III–V Nonlinear Optical Crystals

Orientation‐patterned (OP) III–V semiconductors—used as quasi‐phase‐matched crystals for nonlinear optics applications—are typically epitaxially grown on expensive III–V substrates using a complex process requiring three separate epitaxy steps. In this work, a method is demonstrated for growing orientation‐patterned III–V crystals on Si substrates through V‐groove nanopatterning and a single epitaxial growth. V‐groove Si allows for suppression of random antiphase domain formation that is typical of III–V growth on (001)‐oriented Si through the use of (111)‐faceted trenches patterned on (001)‐oriented Si substrates. By alternating the directions of the V‐groove trenches between [110] and [10], antiphase boundaries can be selectively induced at the boundaries between the two directions of trenches due to the difference in symmetry between the III–V material and Si. This approach allows for a greatly simplified process for growing OP‐III‐Vs and more broadly is a new, robust approach for precisely patterning antiphase boundaries of arbitrary shape and length scales down to 50 nm.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Active Measurement of Plasma Parameters in Over-Dense MFE Reactor Plasma Using High Frequency Beat-Wave Resonantly Generated Thermal Bernstein Waves

UC Davis was funded to explore the combination of beat-wave excitation and laser scattering to make completely non-invasive determination of local plasma conditions. This involves the use of two high frequency mm-wave sources of differing frequencies (the pump waves) to excite plasma waves at the difference of the frequencies, which is known as the beat frequency. The difference frequency or the frequency of the beat wave is chosen so that we produce plasma waves whose fluctuating electron density is a sensitive function of plasma parameters at the location where the source waves overlap. This fluctuation can be detected from outside the plasma, by measuring the deflection of a probe beam directed at the overlap region between the two pump waves (i.e., collective scattering). This work is an extension of mm-wave frequencies of a 10 GHz beat-wave experiment conducted on the DDT toroidal device (see Fig. 1) in a magnetized plasma where the excited Langmuir wave is replaced by the magnetized upper hybrid wave. Here, two high power (> 1 kW) ~94 GHz pump sources will be employed to match the upper hybrid resonance. The primary envisioned experiment is to drive the electron Bernstein wave (EBW) in the interior of a magnetized plasma. The EBW is a hot-plasma electrostatic wave, which cannot propagate in vacuum. Roots of the EBW dispersion relation occur at frequencies slightly higher than harmonics of the electron cyclotron frequency. Use of the EBW has particular advantages for measuring plasma magnetic field in plasmas with relatively high density and low magnetic field, as is the case in spherical tokamaks. This is the “over dense” condition of future efficient high beta magnetic confinement fusion reactors such as the spherical tokamak (STs). The overarching goal of the proposed project is to demonstrate beat wave EBW excitation in low density and temperature parameters that are directly scalable to over-dense spherical tokamaks such as NSTX or MAST, where measurement of magnetic field by electron cyclotron emission (ECE) is experimentally complex, and requires models of plasma profiles in order to interpret results. In such experiments, availability of point measurements of magnetic field strength provided by a beat-wave diagnostic would serve as a powerful independent check on the ECE results, serving to validate the assumed profiles and wave-propagation models used. The first step was to assemble/modify the target chamber in which the beat-wave plasma experiments will take place. Here, the target chamber of the CTIX compact-toroid (spheromak) injector was chosen. At the time that the proposal was submitted, the aim was to conduct the experiments in the CTIX experimental facilities located at the Livermore campus of UC Davis using microwave/mm-wave equipment brought from the main campus., However, by the time DoE funds for this project were awarded, the CTIX experimental facilities located at the Livermore campus of UC Davis were in the process of being shut down due to the closure of the Livermore campus by Lawrence Livermore National Laboratory (LLNL). Consequently, the equipment needed for this project was crated up and shipped from Livermore to Davis as part of the laboratory shutdown after space on campus was secured and College of Engineering approval secured. Once it arrived in Davis on December 20, 2018, work on the project began in earnest with the initial efforts directed toward developing a new plasma target source. Because of the shutdown of the CTIX facility, the scope of the project thus had to be reduced due to the abbreviated amount of time allocated to the project once the shutdown/move had been completed. Chamber modifications have been completed, including the fabrication and installation of a 38 GHz interferometer, as described in Sec. 3. The thermionic plasma source has been tested and characterized. The magnetic field required for the EBW experiments will be generated by permanent magnets mounted on plates positioned above and below the chamber. Magnetic simulations have been performed, and show that a 3-ring magnet arrangement on each plate, with a plate separation of 50 cm, is required to achieve the ~100 G uniform magnetic field in the interaction region that we feel is necessary to ensure good EBW coupling. We proved unable to successfully revive either of Prof. Luhmann’s two high power (≥20 kW) 94 GHz millimeter-wave sources after their extended time in storage and decided to postpone further efforts and concentrate on slightly lower power sources. Consequently, we were successful in reviving a ~2 kW extended klystron amplifier (EIK), which together with a ~2 kW extended interaction oscillator (EIO) will provide the two pump sources required for the project. Safe operation of the EIO requires extensive modifications to one of Prof. Luhmann’s pulse modulators - these modulator modifications are ongoing, which has unfortunately delayed the planned beat-wave experiments. Details of this work are provided in Sec. 3. Even though project funding has now ended, we still plan to proceed to complete the EIO pulse modulator modifications, install the two ~94 GHz pump wave sources on the target chamber, and conduct unmagnetized tests of the beat-wave concept which involve coupling to unmagnetized electron plasma waves rather than magnetized EBWs. Should these experiments prove successful, we plan to submit a follow-up proposal for additional funding which would allow us to undertake the magnetized EBW experiments.

42 ENGINEERING↗
DOE OSTI2012

Relationship between x-ray illumination field size and flat field intensity and its impacts on x-ray imaging

Purpose: X-ray cone-beam CT (CBCT) is being increasingly used for various clinical applications, while its performance is still hindered by image artifacts. This work investigates a new source of reconstruction error, which is often overlooked in the current CBCT imaging. The authors find that the x-ray flat field intensity (I{sub 0}) varies significantly as the illumination volume size changes at different collimator settings. A wrong I{sub 0} value leads to inaccurate CT numbers of reconstructed images as well as wrong scatter measurements in the CBCT research. Methods: The authors argue that the finite size of x-ray focal spot together with the detector glare effect cause the I{sub 0} variation at different illumination sizes. Although the focal spot of commercial x-ray tubes typically has a nominal size of less than 1 mm, the off-focal-spot radiation covers an area of several millimeters on the tungsten target. Due to the large magnification factor from the field collimator to the detector, the penumbra effects of the collimator blades result in different I{sub 0} values for different illumination field sizes. Detector glare further increases the variation, since one pencil beam of incident x-ray is scattered into an area of several centimeters on the detector. In this paper, the authors study these two effects by measuring the focal spot distribution with a pinhole assembly and the detector point spread function (PSF) with an edge-spread function method. The authors then derive a formula to estimate the I{sub 0} value for different illumination field sizes, using the measured focal spot distribution and the detector PSF. Phantom studies are carried out to investigate the accuracy of scatter measurements and CT images with and without considering the I{sub 0} variation effects. Results: On our tabletop system with a Varian Paxscan 4030CB flat-panel detector and a Varian RAD-94 x-ray tube as used on a clinical CBCT system, the focal spot distribution has a measured full-width-at-half-maximum (FWHM) of around 0.4 mm, while non-negligible off-focal-spot radiation is observed at a distance of over 2 mm from the center. The measured detector PSF has an FWHM of 0.510 mm, with a shape close to Gaussian. From these two distributions, the author calculate the estimated I{sub 0} values at different collimator settings. The I{sub 0} variation mainly comes from the focal spot effect. The estimation matches well with the measurements at different collimator widths in both horizontal and vertical directions, with an average error of less than 3%. Our method improves the accuracy of conventional scatter measurements, where the scatter is measured as the difference between fan-beam and cone-beam projections. On a uniform water cylinder phantom, more accurate I{sub 0} suppresses the unfaithful high-frequency signals at the object boundaries of the measured scatter, and the SPR estimation error is reduced from 0.158 to 0.014. The proposed I{sub 0} estimation also reduces the reconstruction error from about 20 HU on the Catphan Copyright-Sign 600 phantom in the selected regions of interest to less than 4 HU. Conclusions: The I{sub 0} variation is identified as one additional error source in x-ray imaging. By measuring the focal-spot distribution and detector PSF, the authors propose an accurate method of estimating the I{sub 0} value for different illumination field sizes. The method obtains more accurate scatter measurements and therefore facilitates scatter correction algorithm designs. As correction methods for other CBCT artifacts become more successful, our research is significant in further improving the CBCT imaging accuracy.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

SU-E-T-389: Effect of Interfractional Shoulder Motion On Low Neck Nodal Targets for Patients Treated Using Volume Modulated Arc Therapy (VMAT)

Purpose: To quantify the dosimetric impact of interfractional shoulder motion on targets in the low neck for head and neck patients treated with volume modulated arc therapy (VMAT). Methods: Three patients with head and neck cancer were selected. All three required treatment to nodal regions in the low neck in addition to the primary tumor. The patients were immobilized during simulation and treatment with a custom thermoplastic mask covering the head and shoulders. One VMAT plan was created for each patient utilizing two full 360° arcs. A second plan was created consisting of two superior VMAT arcs matched to an inferior static AP supraclavicular field. A CT-on-rails alignment verification was performed weekly during each patient's treatment course. The weekly CT images were registered to the simulation CT and the target contours were deformed and applied to the weekly CT. The two VMAT plans were copied to the weekly CT datasets and recalculated to obtain the dose to the low neck contours. Results: The average observed shoulder position shift in any single dimension relative to simulation was 2.5 mm. The maximum shoulder shift observed in a single dimension was 25.7 mm. Low neck target mean doses, normalized to simulation and averaged across all weekly recalculations were 0.996, 0.991, and 1.033 (Full VMAT plan) and 0.986, 0.995, and 0.990 (Half-Beam VMAT plan) for the three patients, respectively. The maximum observed deviation in target mean dose for any individual weekly recalculation was 6.5%, occurring with the Full VMAT plan for Patient 3. Conclusion: Interfractional variation in dose to low neck nodal regions was quantified for three head and neck patients treated with VMAT. Mean dose was 3.3% higher than planned for one patient using a Full VMAT plan. A Half-Beam technique is likely a safer choice when treating the supraclavicular region with VMAT.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1999

Solar-B X-ray Telescope (XRT) Concept Study Report

The X-ray observations from the Yohkoh SXT provided the greatest step forward in our understanding of the solar corona in nearly two decades. Expanding on the accomplishments of Yohkoh, we believe that the scientific objectives of the Solar-B mission are achieved with a significantly improved X-ray telescope (XRT) similar to the SXT. The Solar-B XRT will have twice the spatial resolution and a broader temperature response, while building on the knowledge gained from the successful Yohkoh mission. We present the scientific justification for this view, discuss the instrumental requirements that flow from the scientific objectives, and describe the instrumentation to meet these requirements. We then provide a detailed discussion of the design activities carried out during Phase A, noting the conclusions that were reached in terms of their implications for the detailed design activities which are now commencing. Details of the instrument that have changed as a result of the Phase A studied are specifically noted, and areas of concern going into Phase B are highlighted. XRT is a grazing-incidence (GI) modified Wolter I X-ray telescope, of 35cm inner diameter and 2.7m focal length. The 2048x2048 back-illuminated CCD (now an ISAS responsibility) has 13.5 micron pixels, corresponding to 1.0 arcsec and giving full Sun field of view. This will be the highest resolution GI X-ray telescope ever flown for Solar coronal studies, and it has been designed specifically to observe both the high and low temperature coronal plasma. A small optical telescope provides visible light images for co-alignment with the Solar-B optical and EUV instruments. The XRT science team is working in close cooperation with our Japanese colleagues in the design and construction of this instrument. All of the expertise and resources of the High Energy and Solar/Stellar Divisions of the Center for Astrophysics are being made available to this program, and our team will carry its full share of responsibility for mission operations, data reduction and education and public outreach. All aspects of the XRT design were reviewed during Phase A. The study focussed particularly on those aspects that have the greatest affect on instrument performance and extended lifetime, on the image quality error budget, and on the camera (mechanical and electrical) interface and the instrument mounting interfaces. The present instrument design differs in some details from that originally proposed. Selection of the XRT for Phase A study was contingent upon the removal of the camera and its associated electronics, and the acceptance of a stringent cost cap. The removal of the electronics left the XRT without control electronics for the instrument mechanisms. A mechanism controller was therefore added. The removal of the camera resulted in major complications to the integration and test plan. After many discussions, it was decided that the system would be less expensive, and the risk of unacceptable performance lower, if we include a focus mechanism. The remainder of the XRT design baseline matches the proposed configuration. Data requirements for the XRT are driven by the science plans, which are based on the physical processes in the solar outer atmosphere. Discussions to date of the XRT observing plan, both alone and in conjunction with the other Solar-B instruments, shows that the XRT needs 2 Gbits of on-board storage, at least one circulating buffer of 640 Mbits, and twelve 10- minute downlinks per day in order to carry out its required programs.

Golub, Leon↗
DOE OSTI2011

Deformable Image Registration of Liver With Consideration of Lung Sliding Motion

Purpose: A feature based deformable registration model with sliding transformation was developed in the upper abdominal region for liver cancer. Methods: A two-step thin-plate spline (bi-TPS) algorithm was implemented to deformably register the liver organ. The first TPS registration was performed to exclusively quantify the sliding displacement component. A manual segmentation of the thoracic and abdominal cavity was performed as a priori knowledge. Tissue feature points were automatically identified inside the segmented contour on the images. The scale invariant feature transform method was utilized to match feature points that served as landmarks for the subsequent TPS registration to derive the sliding displacement vector field. To a good approximation, only motion along superior/inferior (SI) direction of voxels on each slice was averaged to obtain the sliding displacement for each slice. A second TPS transformation, as the last step, was carried out to obtain the local deformation field. Manual identification of bifurcation on liver, together with the manual segmentation of liver organ, was employed as a ''ground truth'' for assessing the algorithm's performance. Results: The proposed two-step TPS was assessed with six liver patients. The average error of liver bifurcation between manual identification and calculation for these patients was less than 1.8 mm. The residual errors between manual contour and propagated contour of liver organ using the algorithm fell in the range between 2.1 and 2.8 mm. An index of Dice similarity coefficient (DSC) between manual contour and calculated contour for liver tumor was 93.6% compared with 71.2% from the conventional TPS calculation. Conclusions: A high accuracy ({approx}2 mm) of the two-step feature based TPS registration algorithm was achievable for registering the liver organ. The discontinuous motion in the upper abdominal region was properly taken into consideration. Clinical implementation of the algorithm will find broad application in radiation therapy of liver cancer.

62 RADIOLOGY AND NUCLEAR MEDICINE↗