Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ORTHOGONAL FUNCTION”

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 973 records · Page 54

Planck 2013 results. XXIV. Constraints on primordial non-Gaussianity

The Planck nominal mission cosmic microwave background (CMB) maps yield unprecedented constraints on primordial non-Gaussianity (NG). Using three optimal bispectrum estimators, separable template-fitting (KSW), binned, and modal, we obtain consistent values for the primordial local, equilateral, and orthogonal bispectrum amplitudes, quoting as our final result f NL local = 2.7 ± 5.8, f NL equil = -42 ± 75, and f NL orth = -25 ± 39 (68% CL statistical). Non-Gaussianity is detected in the data; using skew-C &ell; statistics we find a nonzero bispectrum from residual point sources, and the integrated-Sachs-Wolfe-lensing bispectrum at a level expected in the ΛCDM scenario. The results are based on comprehensive cross-validation of these estimators on Gaussian and non-Gaussian simulations, are stable across component separation techniques, pass an extensive suite of tests, and are confirmed by skew-C &ell; , wavelet bispectrum and Minkowski functional estimators. Beyond estimates of individual shape amplitudes, we present model-independent, three-dimensional reconstructions of the Planck CMB bispectrum and thus derive constraints on early-Universe scenarios that generate primordial NG, including general single-field models of inflation, excited initial states (non-Bunch-Davies vacua), and directionally-dependent vector models. We provide an initial survey of scale-dependent feature and resonance models. These results bound both general single-field and multi-field model parameter ranges, such as the speed of sound, c s ≥ 0.02 (95% CL), in an effective field theory parametrization, and the curvaton decay fraction r D ≥ 0.15 (95% CL). The Planck data significantly limit the viable parameter space of the ekpyrotic/cyclic scenarios. The amplitude of the four-point function in the local model τ NL < 2800 (95% CL). Taken together, these constraints represent the highest precision tests to date of physical mechanisms for the origin of cosmic structure.

79 ASTRONOMY AND ASTROPHYSICS↗

On implementation of the extended interior penalty function

The extended interior penalty function formulation is implemented. A rational method for determining the transition between the interior and extended parts is set forth. The formulation includes a straightforward method for avoiding design points with some negative components, which are physically meaningless in structural analysis. The technique, when extended to problems involving parametric constraints, can facilitate closed form integration of the penalty terms over the most important parts of the parameter interval. The method lends itself well to the use of approximation concepts, such as design variable linking, constraint deletion and Taylor series expansions of response quantities in terms of design variables. Examples demonstrating the algorithm, in the context of planar orthogonal frames subjected to ground motion, are included.

Cassis, J. H.↗

An Exploratory Exercise in Taguchi Analysis of Design Parameters: Application to a Shuttle-to-space Station Automated Approach Control System

The chief goals of the summer project have been twofold - first, for my host group and myself to learn as much of the working details of Taguchi analysis as possible in the time allotted, and, secondly, to apply the methodology to a design problem with the intention of establishing a preliminary set of near-optimal (in the sense of producing a desired response) design parameter values from among a large number of candidate factor combinations. The selected problem is concerned with determining design factor settings for an automated approach program which is to have the capability of guiding the Shuttle into the docking port of the Space Station under controlled conditions so as to meet and/or optimize certain target criteria. The candidate design parameters under study were glide path (i.e., approach) angle, path intercept and approach gains, and minimum impulse bit mode (a parameter which defines how Shuttle jets shall be fired). Several performance criteria were of concern: terminal relative velocity at the instant the two spacecraft are mated; docking offset; number of Shuttle jet firings in certain specified directions (of interest due to possible plume impingement on the Station's solar arrays), and total RCS (a measure of the energy expended in performing the approach/docking maneuver). In the material discussed here, we have focused on single performance criteria - total RCS. An analysis of the possibility of employing a multiobjective function composed of a weighted sum of the various individual criteria has been undertaken, but is, at this writing, incomplete. Results from the Taguchi statistical analysis indicate that only three of the original four posited factors are significant in affecting RCS response. A comparison of model simulation output (via Monte Carlo) with predictions based on estimated factor effects inferred through the Taguchi experiment array data suggested acceptable or close agreement between the two except at the predicted optimum point, where a difference outside a rule-of-thumb bound was observed. We have concluded that there is most likely an interaction effect not provided for in the original orthogonal array selected as the basis for our experimental design. However, we feel that the data indicates that this interaction is a mild one and that inclusion of its effect will not alter the location of the optimum.

Deal, Don E.↗

Convergence of Weak-SINDy Surrogate Models

In this paper, we give an in-depth error analysis for surrogate models generated by a variant of the Sparse Identification of Nonlinear Dynamics (SINDy) method. We start with an overview of a variety of nonlinear system identification techniques, namely SINDy, weak-SINDy, and the occupation kernel method. Under the assumption that the dynamics are a finite linear combination of a set of basis functions, these methods establish a linear system to recover coefficients. We illuminate the structural similarities between these techniques and establish a projection property for the weak-SINDy technique. Following the overview, we analyze the error of surrogate models generated by a simplified version of weak-SINDy. In particular, under the assumption of boundedness of a composition operator given by the solution, we show that (i) the surrogate dynamics converges towards the true dynamics and (ii) the solution of the surrogate model is reasonably close to the true solution. Finally, as an application, we discuss the use of a combination of weak-SINDy surrogate modeling and proper orthogonal decomposition (POD) to build a surrogate model for partial differential equations (PDEs).

97 MATHEMATICS AND COMPUTING↗

Analytical gradients for molecular-orbital-based machine learning

We report molecular-orbital-based machine learning (MOB-ML) enables the prediction of accurate correlation energies at the cost of obtaining molecular orbitals. Here, we present the derivation, implementation, and numerical demonstration of MOB-ML analytical nuclear gradients, which are formulated in a general Lagrangian framework to enforce orthogonality, localization, and Brillouin constraints on the molecular orbitals. The MOB-ML gradient framework is general with respect to the regression technique (e.g., Gaussian process regression or neural networks) and the MOB feature design. We show that MOB-ML gradients are highly accurate compared to other ML methods on the ISO17 dataset while only being trained on energies for hundreds of molecules compared to energies and gradients for hundreds of thousands of molecules for the other ML methods. The MOB-ML gradients are also shown to yield accurate optimized structures at a computational cost for the gradient evaluation that is comparable to a density-corrected density functional theory calculation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Design, characterization, and fabrication of solar-retroreflective cool-wall materials

Raising urban albedo increases the fraction of incident sunlight returned to outer space, cooling cities and their buildings. We evaluated the angular distribution of solar radiation incident on exterior walls in 17 U S. climates to develop performance parameters for solar-retroreflective walls, then applied first-principle physics and ray-tracing simulations to explore designs. Our analysis indicates that retroreflective walls must function at large incidence angles to reflect a substantial portion of summer sunlight, and that this will be difficult to attain with materials that rely on total internal reflection. Gonio-spectrophotometer measurements of the solar spectral bi-directional reflectivity of a bicycle reflector showed little to no retroreflection at large incidence angles. Visual comparisons of retroreflection to specular first-surface reflection for four different retroreflective safety films using violet and green lasers suggest their retroreflection to be no greater than 0.09 at incidence angles up to 45°, and no greater than 0.30 at incidence angles of up to 70°. Attempts to produce a two-surface retroreflector with orthogonal mirror grooves by cutting and polishing an aluminum block indicate that residual surface roughness impedes retroreflection. Ongoing efforts focus on forming orthogonal surfaces with aluminized Mylar film, a material with very high specular reflectance across the solar spectrum. We investigated (1) folding or stamping a free film; (2) adhering the film to a pre-shaped substrate; or (3) attaching the film to a flat ductile substrate, then shaping. The latter two methods were more successful but yielded imperfect right angles.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Optically and thermally driven huge lattice orbital and spin angular momenta from spinning fullerenes

Lattice vibration in solids may carry angular momentum. But unlike the intrinsic spin of electrons, the lattice vibration is rarely rotational. To induce angular momentum, one needs to find a material that can accommodate a twisted normal mode, two orthogonal modes, or excitation of magnons. If excitation is too strong, one may exceed the Lindemann limit, so the material melts. Therefore these methods are not ideal. Here, we theoretically propose a route to phonon angular momentum in a molecular crystal C 60 . We find that a single laser pulse is able to inject a significant amount of angular momentum to C 60 , and the momentum transfer is helicity dependent. Changing from right-circularly polarized light to left-circularly polarized light switches the direction of phonon angular momentum. On the ultrafast timescale, the orbital angular momentum change closely resembles the displacive excitation of coherent phonons, with a cosine-function dependence on time, different from the spin counterpart. Atomic displacements, even under strong laser excitation, remain far below the Lindemann criterion. Under thermal excitation, spinning C 60 even at room temperature generates a huge angular momentum close to several hundred h¯. Our finding opens the door to a large group of fullerenes, from C 60 , C 70 to their endohedral derivatives, where angular momentum can be generated through light or temperature. Furthermore, this paves the way to the phononic control electronic spin and harvesting thermal energy through phonon angular momentum.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Combining glomerular basement membrane and tubular basement membrane assessment improves the prediction of diabetic end‐stage renal disease

Abstract Background To address the prognostic value of combining tubular basement membrane (TBM) and glomerular basement membrane (GBM) thickness in diabetic nephropathy (DN). Methods This retrospective study enrolled 110 patients with type 2 diabetes and biopsy‐proven DN from 2011 to 2018. The pathological findings were confirmed according to the Renal Pathology Society classifications. GBM and TBM thicknesses were determined using the Haas' direct measurement/arithmetic mean method and orthogonal intercept method, respectively. Cox proportional hazard models were used to investigate the hazard ratios (HRs) for the influence of combined GBM and TBM thickness for predicting end‐stage renal disease (ESRD). Results Patients were assigned to three groups according to the median GBM and TBM thickness: GBM lo TBM lo (GBM < 681 nm and TBM < 1200 nm), GBM hi TBM lo /GBM lo TBM hi (GBM ≥ 681 nm and TBM < 1200 nm, or GBM < 681 nm and TBM ≥ 1200 nm), and GBM hi TBM hi (GBM ≥ 681 nm and TBM ≥ 1200 nm). The GBM hi TBM lo /GBM lo TBM hi and GBM hi TBM hi groups displayed poorer renal function, a more severe glomerular classification and interstitial inflammation, and poorer renal survival rates than the GBM lo TBM lo group The GBM hi TBM lo /GBM lo TBM hi and GBM hi TBM hi groups had adjusted HRs of 1.49 (95% confidence interval [CI], 1.21‐9.75) and 3.07 (95% CI, 2.87‐12.78), respectively, compared with the GBM lo TBM lo group. Conclusions TBM thickness enhanced GBM thickness for renal prognosis in patients with type 2 diabetes.

Zhao, Lijun↗

Spectral fluctuations in the Sachdev-Ye-Kitaev model

We present a detailed quantitative analysis of spectral correlations in the Sachdev-Ye-Kitaev (SYK) model. We find that the deviations from universal Random Matrix Theory (RMT) behavior are due to a small number of long-wavelength fluctuations (of the order of the number of Majorana fermions N ) from one realization of the ensemble to the next one. These modes can be parameterized effectively in terms of Q-Hermite orthogonal polynomials, the main contribution being due to scale fluctuations for which we give a simple analytical estimate. Our numerical results for N = 32 show that only the lowest eight polynomials are needed to eliminate the nonuniversal part of the spectral fluctuations. The covariance matrix of the coefficients of this expansion can be obtained analytically from low-order double-trace moments. We evaluate the covariance matrix of the first six moments and find that it agrees with the numerics. We also analyze the spectral correlations in terms of a nonlinear σ-model, which is derived through a Fierz transformation, and evaluate the one and two-point spectral correlation functions to two-loop order. The wide correlator is given by the sum of the universal RMT result and corrections whose lowest-order term corresponds to scale fluctuations. However, the loop expansion of the σ-model results in an ill-behaved expansion of the resolvent, and it gives universal RMT fluctuations not only for q= 4 or higher even q-body interactions, but also for the q= 2 SYK model albeit with a much smaller Thouless energy while the correct result in this case should have been Poisson statistics. In our numerical studies we analyze the number variance and spectral form factor for N = 32 and q = 4. We show that the quadratic deviation of the number variance for large energies appears as a peak for small times in the spectral form factor. After eliminating the long-wavelength fluctuations, we find quantitative agreement with RMT up to an exponentially large number of level spacings for the number variance or exponentially short times in the case of the spectral form factor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Spectra-Orthogonal Optical Anisotropy in Wafer-Scale Molecular Crystal Monolayers

Controlling the spectral and polarization response of two-dimensional (2D) crystals is vital for developing ultrathin platforms for compact optoelectronic devices. However, independently tuning optical anisotropy and spectral response remains challenging in conventional semiconductors due to the intertwined nature of their lattice and electronic structures. Here, we report spectra-orthogonal optical anisotropy─where polarization anisotropy is tuned independently of spectral response─in wafer-scale, one-atom-thick 2D molecular crystal (2DMC) monolayers synthesized on monolayer transition-metal dichalcogenide (TMD) crystals. Utilizing the concomitant spectral consistency and structural tunability of perylene derivatives, we demonstrate tunable optical polarization anisotropy in 2DMCs with similar spectral profiles, as confirmed by room-temperature scanning tunneling microscopy and cross-polarized reflectance microscopy. Additional angle-dependent analysis of the single-crystal and polycrystalline molecular domains reveals an epitaxial relationship between the 2DMC and TMD. In conclusion, our results establish a scalable, molecule-based 2D crystalline platform for unique and tunable functionalities unattainable in covalent 2D solids.

2D materials↗

Skew-Symmetric adjacency matrices for clustering directed graphs

Graph clustering methods often critically rely on the symmetry of graph matrices. Developing analogous methods for digraphs often proves more challenging, because digraph matrices are typically asymmetric and not orthogonally diagonalizable. However, researchers have recently proposed several complex-valued Hermitian digraph matrices. In particular, one such representation has been utilized as an input to an algorithm for finding imbalanced cuts. In this work, we establish an algebraic relationship between this matrix and an associated real-valued matrix. We show using this real-valued matrix for imbalanced cut-finding algorithms is not only sufficient but advantageous. Our algorithm uses less memory and asymptotically less computation while provably preserving solution quality. We also show our method can be easily implemented using standing computational building blocks, possesses better numerical properties, and loans itself to a natural interpretation via an objective function relaxation argument. We empirically demonstrate these advantages on real world data sets and show how our algorithm can uncover meaningful cluster structure.

Hayashi, Koby↗

Data-Driven Reduced-Order Modeling of Convective Heat Transfer in Porous Media

This work presents a data-driven Reduced-Order Model (ROM) for parametric convective heat transfer problems in porous media. The intrusive Proper Orthogonal Decomposition aided Reduced-Basis (POD-RB) technique is employed to reduce the porous medium formulation of the incompressible Reynolds-Averaged Navier–Stokes (RANS) equations coupled with heat transfer. Instead of resolving the exact flow configuration with high fidelity, the porous medium formulation solves a homogenized flow in which the fluid-structure interactions are captured via volumetric flow resistances with nonlinear, semi-empirical friction correlations. A supremizer approach is implemented for the stabilization of the reduced fluid dynamics equations. The reduced nonlinear flow resistances are treated using the Discrete Empirical Interpolation Method (DEIM), while the turbulent eddy viscosity and diffusivity are approximated by adopting a Radial Basis Function (RBF) interpolation-based approach. The proposed method is tested using a 2D numerical model of the Molten Salt Fast Reactor (MSFR), which involves the simulation of both clean and porous medium regions in the same domain. For the steady-state example, five model parameters are considered to be uncertain: the magnitude of the pumping force, the external coolant temperature, the heat transfer coefficient, the thermal expansion coefficient, and the Prandtl number. For transient scenarios, on the other hand, the coastdown-time of the pump is the only uncertain parameter. The results indicate that the POD-RB-ROMs are suitable for the reduction of similar problems. The relative L2 errors are below 3.34% for every field of interest for all cases analyzed, while the speedup factors vary between 54 (transient) and 40,000 (steady-state).

German, Péter (ORCID:0000000307285283)↗

cMonkey2: Automated, systematic, integrated detection of co-regulated gene modules for any organism

The cMonkey integrated biclustering algorithm identifies conditionally co-regulated modules of genes (biclusters). cMonkey integrates various orthogonal pieces of information which support evidence of gene co-regulation, and optimizes biclusters to be supported simultaneously by one or more of these prior constraints. The algorithm served as the cornerstone for constructing the first global, predictive Environmental Gene Regulatory Influence Network (EGRIN) model for a free-living cell, and has now been applied to many more organisms. However, due to its computational inefficiencies, long run-time and complexity of various input data types, cMonkey was not readily usable by the wider community. To address these primary concerns, we have significantly updated the cMonkey algorithm and refactored its implementation, improving its usability and extendibility. These improvements provide a fully functioning and user-friendly platform for building co-regulated gene modules and the tools necessary for their exploration and interpretation. We show, via three separate analyses of data for E. coli, M. tuberculosis and H. sapiens, that the updated algorithm and inclusion of novel scoring functions for new data types (e.g. ChIP-seq and transcription factor over-expression [TFOE]) improve discovery of biologically informative co-regulated modules. The complete cMonkey2 software package, including source code, is available at https://github.com/baliga-lab/cmonkey2.

59 BASIC BIOLOGICAL SCIENCES↗

Vestibular convergence patterns in vestibular nuclei neurons of alert primates

Sensory signal convergence is a fundamental and important aspect of brain function. Such convergence may often involve complex multidimensional interactions as those proposed for the processing of otolith and semicircular canal (SCC) information for the detection of translational head movements and the effective discrimination from physically congruent gravity signals. In the present study, we have examined the responses of primate rostral vestibular nuclei (VN) neurons that do not exhibit any eye movement-related activity using 0.5-Hz translational and three-dimensional (3D) rotational motion. Three distinct neural populations were identified. Approximately one-fourth of the cells exclusively encoded rotational movements (canal-only neurons) and were unresponsive to translation. The canal-only central neurons encoded head rotation in SCC coordinates, exhibited little orthogonal canal convergence, and were characterized with significantly higher sensitivities to rotation as compared to primary SCC afferents. Another fourth of the neurons modulated their firing rates during translation (otolith-only cells). During rotations, these neurons only responded when the axis of rotation was earth-horizontal and the head was changing orientation relative to gravity. The remaining one-half of VN neurons were sensitive to both rotations and translations (otolith + canal neurons). Unlike primary otolith afferents, however, central neurons often exhibited significant spatiotemporal (noncosine) tuning properties and a wide variety of response dynamics to translation. To characterize the pattern of SCC inputs to otolith + canal neurons, their rotational maximum sensitivity vectors were computed using exclusively responses during earth-vertical axis rotations (EVA). Maximum sensitivity vectors were distributed throughout the 3D space, suggesting strong convergence from multiple SCCs. These neurons were also tested with earth-horizontal axis rotations (EHA), which would activate both vertical canals and otolith organs. However, the recorded responses could not be predicted from a linear combination of EVA rotational and translational responses. In contrast, one-third of the neurons responded similarly during EVA and EHA rotations, although a significant response modulation was present during translation. Thus this subpopulation of otolith + canal cells, which included neurons with either high- or low-pass dynamics to translation, appear to selectively ignore the component of otolith-selective activation that is due to changes in the orientation of the head relative to gravity. Thus contrary to primary otolith afferents and otolith-only central neurons that respond equivalently to tilts relative to gravity and translational movements, approximately one-third of the otolith + canal cells seem to encode a true estimate of the translational component of the imposed passive head and body movement.

Non-NASA Center↗

Charge transfer excitations from excited state Hartree-Fock subsequent minimization scheme

Photoinduced charge-transfer processes play a key role for novel photovoltaic phenomena and devices. Thus, the development of ab initio methods that allow for an accurate and computationally inexpensive treatment of charge-transfer excitations is a topic that nowadays attracts a lot of scientific attention. In this paper we extend an approach recently introduced for the description of single and double excitations [M. Tassi, I. Theophilou, and S. Thanos, Int. J. Quantum Chem. 113, 690 (2013); M. Tassi, I. Theophilou, and S. Thanos, J. Chem. Phys. 138, 124107 (2013)] to allow for the description of intermolecular charge-transfer excitations. We describe an excitation where an electron is transferred from a donor system to an acceptor one, keeping the excited state orthogonal to the ground state and avoiding variational collapse. These conditions are achieved by decomposing the space spanned by the Hartree-Fock (HF) ground state orbitals into four subspaces: The subspace spanned by the occupied orbitals that are localized in the region of the donor molecule, the corresponding for the acceptor ones and two more subspaces containing the virtual orbitals that are localized in the neighborhood of the donor and the acceptor, respectively. Next, we create a Slater determinant with a hole in the subspace of occupied orbitals of the donor and a particle in the virtual subspace of the acceptor. Subsequently we optimize both the hole and the particle by minimizing the HF energy functional in the corresponding subspaces. Finally, we test our approach by calculating the lowest charge-transfer excitation energies for a set of tetracyanoethylene-hydrocarbon complexes that have been used earlier as a test set for such kind of excitations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fluorescence-Enabled Self-Reporting for Redox Flow Batteries

Monitoring battery health is challenging. Self-reporting enables the rapid health assessment of redox flow batteries (RFBs) and provides insight into degradation mechanisms of electrochemically active molecules (redoxmers). In this study, we introduce fluorescence as an orthogonal property to probe this chemistry in situ. An established anolyte redoxmer, 2,1,3-benzothiadiazole, is rendered an efficient blue-green fluorophore through minimalistic derivatization. One of the derivatives is electrochemically reversible with a long lifetime and cycling stability in the charged state. Spectroscopic measurements on this new redoxmer reveal strong effects of the electrolyte cation on the fluorescence that are useful for probing the solvent microenvironment. Using this probe, we demonstrate proof-of-concept in situ crossover sensing and characterize the effects of electrolyte composition on this crossover. In this way, real-time tracing of redoxmers in a flow cell is demonstrated, paving the way to include still more self-reporting functions into the redoxmers.

25 ENERGY STORAGE↗

Bottom-up synthesis of mesoscale nanomeshes of graphene nanoribbons on germanium

The synthesis of functional graphene nanostructures on Ge(001) provides an attractive route toward integrating graphene-based electronic devices onto complementary metal oxide semiconductor-compatible platforms. In this study, we leverage the phenomenon of the anisotropic growth of graphene nanoribbons from rationally placed graphene nanoseeds and their rotational self-alignment during chemical vapor deposition to synthesize mesoscale graphene nanomeshes over areas spanning several hundred square micrometers. Lithographically patterned nanoseeds are defined on a Ge(001) surface at pitches ranging from 50 to 100 nm, which serve as starting sites for subsequent nanoribbon growth. Rotational self-alignment of the nanoseeds followed by anisotropic growth kinetics causes the resulting nanoribbons to be oriented along each of the equivalent, orthogonal Ge < 110 > directions with equal probability. As the nanoribbons grow, they fuse, creating a continuous nanomesh. In contrast to nanomesh synthesis via top-down approaches, this technique yields nanomeshes with atomically faceted edges and covalently bonded junctions, which are important for maximizing charge transport properties. Additionally, we simulate the electrical characteristics of nanomeshes synthesized from different initial nanoseed-sizes, size-polydispersities, pitches, and device channel lengths to identify a parameter-space for acceptable on/off ratios and on-conductance in semiconductor electronics. The simulations show that decreasing seed diameter and pitch are critical to increasing nanomesh on/off ratio and on-conductance, respectively. With further refinements in lithography, nanomeshes obtained via seeded synthesis and anisotropic growth are likely to have superior electronic properties with tremendous potential in a multitude of applications, such as radio frequency communications, sensing, thin-film electronics, and plasmonics.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

TU-FG-209-12: Treatment Site and View Recognition in X-Ray Images with Hierarchical Multiclass Recognition Models

Purpose: To investigate an approach of automatically recognizing anatomical sites and imaging views (the orientation of the image acquisition) in 2D X-ray images. Methods: A hierarchical (binary tree) multiclass recognition model was developed to recognize the treatment sites and views in x-ray images. From top to bottom of the tree, the treatment sites are grouped hierarchically from more general to more specific. Each node in the hierarchical model was designed to assign images to one of two categories of anatomical sites. The binary image classification function of each node in the hierarchical model is implemented by using a PCA transformation and a support vector machine (SVM) model. The optimal PCA transformation matrices and SVM models are obtained by learning from a set of sample images. Alternatives of the hierarchical model were developed to support three scenarios of site recognition that may happen in radiotherapy clinics, including two or one X-ray images with or without view information. The performance of the approach was tested with images of 120 patients from six treatment sites – brain, head-neck, breast, lung, abdomen and pelvis – with 20 patients per site and two views (AP and RT) per patient. Results: Given two images in known orthogonal views (AP and RT), the hierarchical model achieved a 99% average F1 score to recognize the six sites. Site specific view recognition models have 100 percent accuracy. The computation time to process a new patient case (preprocessing, site and view recognition) is 0.02 seconds. Conclusion: The proposed hierarchical model of site and view recognition is effective and computationally efficient. It could be useful to automatically and independently confirm the treatment sites and views in daily setup x-ray 2D images. It could also be applied to guide subsequent image processing tasks, e.g. site and view dependent contrast enhancement and image registration. The senior author received research grants from ViewRay Inc. and Varian Medical System.

60 APPLIED LIFE SCIENCES↗