SEARCH · Engineering Papers

Results for “Drift vector”

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.

SourceResearch connectionFrom yearThrough year

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

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Reduced Common-Mode Voltage PWM Scheme for Full-SiC Three-Level Uninterruptible Power Supply With Small DC-Link Capacitors

In this article, a pulsewidth modulation (PWM) scheme for three-level full-SiC uninterruptible power supplies is developed to achieve a high power density. Two key passive components are selected for size reduction of the ac-ac stage: common-mode (CM) EMI filter, and dc-link capacitors. To reduce the CM noise, a new vector combination is proposed based on synchronous switching among three-phases. The proposed combinations align CM voltage (CMV) to be a single pulse per switching period. Owing to the simple shape, CMV cancellation between a three-level rectifier and inverter can be maximally utilized. A transition between the three combinations can control the drift of neutral point voltage. An equivalent carrier-based implementation is developed. Second, a simple algorithm to compensate neutral point voltage fluctuation is proposed both for differential mode (DM) and CM output voltage. Low-order harmonics on three-phase currents and an additional high-frequency CM noise by misaligned switching instant can be eliminated. The proposed compensation can be implemented by a simple correction on carrier slopes and injected zero-sequence voltage. The proposed PWM scheme is verified with 20-kW full-SiC UPS switching at 60 kHz with 140 μF dc-link capacitors.

42 ENGINEERING↗

DOE OSTI2015

SU-E-J-237: Real-Time 3D Anatomy Estimation From Undersampled MR Acquisitions

Recent developments made MRI guided radiotherapy feasible. Performing simultaneous imaging during fractions can provide information about changing anatomy by means of deformable image registration for either immediate plan adaptations or accurate dose accumulation on the changing anatomy. In 3D MRI, however, acquisition time is considerable and scales with resolution. Furthermore, intra-scan motion degrades image quality.In this work, we investigate the sensitivity of registration quality on imageresolution: potentially, by employing spatial undersampling, the acquisition timeof MR images for the purpose of deformable image registration can be reducedsignificantly.On a volunteer, 3D-MR imaging data was sampled in a navigator-gated manner, acquiring one axial volume (360×260×100mm{sup 3}) per 3s during exhale phase. A T1-weighted FFE sequence was used with an acquired voxel size of (2.5mm{sup 3}) for a duration of 17min. Deformation vector fields were evaluated for 100 imaging cycles with respect to the initial anatomy using deformable image registration based on optical flow. Subsequently, the imaging data was downsampled by a factor of 2, simulating a fourfold acquisition speed. Displacements of the downsampled volumes were then calculated by the same process.In kidneyliver boundaries and the region around stomach/duodenum, prominent organ drifts could be observed in both the original and the downsampled imaging data. An increasing displacement of approximately 2mm was observed for the kidney, while an area around the stomach showed sudden displacements of 4mm. Comparison of the motile points over time showed high reproducibility between the displacements of high-resolution and downsampled volumes: over a 17min acquisition, the componentwise RMS error was not more than 0.38mm.Based on the synthetic experiments, 3D nonrigid image registration shows little sensitivity to image resolution and the displacement information is preserved even when halving the resolution. This can be employed to greatly reduce image acquisition times for interventional applications in real-time. This work was funded by the SoRTS consortium, which includes the industry partners Elekta, Philips and Technolution.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2024

Self-Supervised Cloud Classification

Abstract Low-level marine clouds play a pivotal role in Earth’s weather and climate through their interactions with radiation, heat and moisture transport, and the hydrological cycle. These interactions depend on a range of dynamical and microphysical processes that result in a broad diversity of cloud types and spatial structures, and a comprehensive understanding of cloud morphology is critical for continued improvement of our atmospheric modeling and prediction capabilities moving forward. Deep learning has recently accelerated our ability to study clouds using satellite remote sensing, and machine learning classifiers have enabled detailed studies of cloud morphology. A major limitation of deep learning approaches to this problem, however, is the large number of hand-labeled samples that are required for training. This work applies a recently developed self-supervised learning scheme to train a deep convolutional neural network (CNN) to map marine cloud imagery to vector embeddings that capture information about mesoscale cloud morphology and can be used for satellite image classification. The model is evaluated against existing cloud classification datasets and several use cases are demonstrated, including training cloud classifiers with very few labeled samples, interrogation of the CNN’s learned internal feature representations, cross-instrument application, and resilience against sensor calibration drift and changing scene brightness. The self-supervised approach learns meaningful internal representations of cloud structures and achieves comparable classification accuracy to supervised deep learning methods without the expense of creating large hand-annotated training datasets. Significance Statement Marine clouds heavily influence Earth’s weather and climate, and improved understanding of marine clouds is required to improve our atmospheric modeling capabilities and physical understanding of the atmosphere. Recently, deep learning has emerged as a powerful research tool that can be used to identify and study specific marine cloud types in the vast number of images collected by Earth-observing satellites. While powerful, these approaches require hand-labeling of training data, which is prohibitively time intensive. This study evaluates a recently developed self-supervised deep learning method that does not require human-labeled training data for processing images of clouds. We show that the trained algorithm performs competitively with algorithms trained on hand-labeled data for image classification tasks. We also discuss potential downstream uses and demonstrate some exciting features of the approach including application to multiple satellite instruments, resilience against changing image brightness, and its learned internal representations of cloud types. The self-supervised technique removes one of the major hurdles for applying deep learning to very large atmospheric datasets.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1979

A digital scene matching technique for geometric image correction and autonomous navigation

A technique is described for precise registration of two images of the same area, taken under different conditions. The technique, called AUTO-MATCH, involves digital preprocessing of the images to extract edge contours, followed by a correlation of corresponding edge end points. The availability of an array of endpoints makes possible subpixed registration accuracy. The technique is applied to a system for assessment of geometric image quality, to be installed at NASA-Goddard. This system measures the registration vectors over an array of window pairs from LANDSAT images, so as to determine the distortion between them. Similar measurements could be used to infer relative positioning of the spacecraft. The technique is considered for updating the inertial navigation systems of missiles or aircraft. Scene matching between an image obtained aboard the vehicle and a stored reference can eliminate the drift of an inertial platform in order to demonstrate this capability in real time. A laboratory demonstration of AUTO-MATCH is described.

Tisdale, G. E.↗

NASA NTRS1990

Attitude determination and parameter estimation using vector observations - Application

This paper presents tests of a new method for the simultaneous estimation of spacecraft attitude and sensor biases, based on a quaternion estimation algorithm minimizing Wahba's loss function. The new method is compared with a conventional batch least-squares differential correction algorithm. The estimates are based on data from strapdown gyros and star trackers, simulated with varying levels of Gaussian noise for both inertially-fixed and earth-pointing attitudes. Both algorithms solve for the spacecraft attitude and the gyro drift rate biases. In the majority of tests performed, the two methods converge to the same estimates in the same number of iterations, but the new algorithm requires about 60 percent more computational effort. Some cases were found in which the new method converges in fewer iterations than the differential correction, and some for which the differential correction requires fewer iterations.

Markley, F. Landis↗

DOE OSTI2026

Stability and Characteristics of Lower-hybrid Drift Waves: Dependence on Electron Beta and Cross-field Relative Drift

Lower-hybrid drift waves (LHDWs) are frequently observed microinstabilities in both space and laboratory plasmas. Despite decades of study, the relationship between electrostatic (ES-LHDW) and electromagnetic (EM-LHDW) variants and the plasma parameters controlling their stability remains unclear. Here, we systematically examine LHDW behavior by solving the local linear dispersion relation over a wide range of plasma and field conditions. Our results demonstrate that ES-LHDWs and EM-LHDWs are not distinct modes but rather two different regimes of the same drift wave whose character evolves smoothly with electron beta (β e ) and the cross-field electron drift velocity relative to ions, normalized to the ion sound speed (u 0x /C s ). The nature of the waves changes from electrostatic to electromagnetic when β e increases. Growth rates increase with u 0x /C s but decrease with β e , while the most unstable wavelength remains nearly universal, with kρ e ∼ 0.8 (k is the magnitude of the wave vector and ρ e is the electron gyroradius). We further present quasi-linear estimates of nonlinear saturation properties, including energy partition among electric fields, magnetic fields, and particle kinetic responses. We show that ES-LHDWs reach higher electric-field saturation amplitudes, whereas EM-LHDWs generate strong magnetic perturbations and parallel electric fields that may enable efficient particle heating. Comparisons with the classical model reveal that retaining electromagnetic effects is essential for accurate predictions of frequency, growth rate, and the propagation angle. These findings provide a unified framework for understanding LHDWs across diverse collisionless plasma environments, including current sheets of magnetic reconnection, shear layers, collisionless shocks, and boundary regions.

Solar coronal waves↗

NASA NTRS2013

Improving Photometry and Stellar Signal Preservation with Pixel-Level Systematic Error Correction

The Kepler Mission has demonstrated that excellent stellar photometric performance can be achieved using apertures constructed from optimally selected CCD pixels. The clever methods used to correct for systematic errors, while very successful, still have some limitations in their ability to extract long-term trends in stellar flux. They also leave poorly correlated bias sources, such as drifting moiré pattern, uncorrected. We will illustrate several approaches where applying systematic error correction algorithms to the pixel time series, rather than the co-added raw flux time series, provide significant advantages. Examples include, spatially localized determination of time varying moiré pattern biases, greater sensitivity to radiation-induced pixel sensitivity drops (SPSDs), improved precision of co-trending basis vectors (CBV), and a means of distinguishing the stellar variability from co-trending terms even when they are correlated. For the last item, the approach enables physical interpretation of appropriately scaled coefficients derived in the fit of pixel time series to the CBV as linear combinations of various spatial derivatives of the pixel response function (PRF). We demonstrate that the residuals of a fit of soderived pixel coefficients to various PRF-related components can be deterministically interpreted in terms of physically meaningful quantities, such as the component of the stellar flux time series which is correlated with the CBV, as well as, relative pixel gain, proper motion and parallax. The approach also enables us to parameterize and assess the limiting factors in the uncertainties in these quantities.

Kolodzijczak, Jeffrey J.↗

DOE OSTI2024

Investigating instabilities in magnetized low-pressure capacitively coupled RF plasma using particle-in-cell (PIC) simulations

The effect of a uniform magnetic field on particle transport in low-pressure radio frequency (RF) capacitively coupled plasma (CCP) has been studied using a particle-in-cell model. Three distinct regimes of plasma behavior can be identified as a function of the magnetic field. In the first regime at low magnetic fields, asymmetric plasma profiles are observed within the CCP chamber due to the effect of $\overrightarrow{E}$ × $\overrightarrow{B}$ drift. As the magnetic field increases, instabilities develop and form self-organized spoke-shaped structures that are distinctly seen within the bulk plasma closer to the sheath. In this second regime, the spoke-shaped coherent structures rotate inside the plasma chamber in the -$\overrightarrow{E}$ × $\overrightarrow{B}$ direction, where $\overrightarrow{E}$ × $\overrightarrow{B}$ are the DC electric and magnetic field vectors, respectively, and the DC electric field exists in the sheath and pre-sheath regions. The spoke rotation frequency is in the megahertz range. As the magnetic field strength increases further, the rotating coherent spokes continue to exist near the sheath. The coherent structures are, however, accompanied by new small-scale incoherent structures originating and moving within the bulk plasma region away from the sheath. This is the third regime of plasma behavior. The threshold values of the magnetic field between these regimes were found not to vary with changing plasma reactor geometry (e.g., area ratio between ground and powered electrodes) or the use of an external capacitor between the RF-powered electrode and the RF source. The threshold values of the magnetic field between these regimes shift toward higher values with increasing gas pressure. Analysis of the results indicates that the rotating structures are due to the lower hybrid instability driven by density gradients and electron-neutral collisions. This paper provides guidance on the upper limit of the magnetic field for instability-free operation in low-pressure CCP-based semiconductor deposition and etch systems that use the external magnetic field for plasma uniformity control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1993

MeV Gamma Ray Detection Algorithms for Stacked Silicon Detectors

By making use of the signature of a gamma ray event as it appears in N = 5 to 20 lithium-drifted silicon detectors and applying smart selection algorithms, gamma rays in the energy range of 1 to 8 MeV can be detected with good efficiency and selectivity. Examples of the types of algorithms used for different energy regions include the simple sum mode, the sum-coincidence mode used in segmented detectors, unique variations on sum-coincidence for an N-dimensional vector event, and a new and extremely useful mode for double escape peak spectroscopy at pair-production energies. The latter algorithm yields a spectrum similar to that of the pair spectrometer, but without the need of the dual external segments for double escape coincidence, and without the large loss in efficiency of double escape events. Background events due to Compton scattering are largely suppressed. Monte Carlo calculations were used to model the gamma ray interactions in the silicon, in order to enable testing of a wide array of different algorithms on the event N-vectors for a large-N stack.

McMurray, Robert E. Jr.↗

DOE OSTI2015

ION KINETIC ENERGY CONSERVATION AND MAGNETIC FIELD STRENGTH CONSTANCY IN MULTI-FLUID SOLAR WIND ALFVÉNIC TURBULENCE

We investigate the properties of plasma fluid motion in the large-amplitude, low-frequency fluctuations of highly Alfvénic fast solar wind. We show that protons locally conserve total kinetic energy when observed from an effective frame of reference comoving with the fluctuations. For typical properties of the fast wind, this frame can be reasonably identified by alpha particles which, due to their drift with respect to protons at about the Alfvén speed along the magnetic field, do not partake in the fluid low-frequency fluctuations. Using their velocity to transform the proton velocity into the frame of Alfvénic turbulence, we demonstrate that the resulting plasma motion is characterized by a constant absolute value of the velocity, zero electric fields, and aligned velocity and magnetic field vectors as expected for unidirectional Alfvénic fluctuations in equilibrium. We propose that this constraint, via the correlation between velocity and magnetic field in Alfvénic turbulence, is the origin of the observed constancy of the magnetic field; while the constant velocity corresponding to constant energy can only be observed in the frame of the fluctuations, the corresponding constant total magnetic field, invariant for Galilean transformations, remains the observational signature in the spacecraft frame of the constant total energy in the Alfvén turbulence frame.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Reconstruction for Liquid Argon TPC Neutrino Detectors Using Parallel Architectures.

Neutrinos are particles that interact rarely, so identifying them re-quires large detectors which produce lots of data. Processing this data with thecomputing power available is becoming more dicult as the detectors increasein size to reach their physics goals. In liquid argon time projection chambers(TPCs) the charged particles from neutrino interactions produce ionization elec-trons which drift in an electric field towards a series of collection wires, and thesignal on the wires is used to reconstruct the interaction. The MicroBooNEdetector currently collecting data at Fermilab has 8000 wires, and planned fu-ture experiments like DUNE will have 100 times more, which means that thetime required to reconstruct an event will scale accordingly. Modernization ofliquid argon TPC reconstruction code, including vectorization, parallelizationand code portability to GPUs, will help to mitigate these challenges. The liquidargon TPC hit finding algorithm within theLArSoftframework used acrossmultiple experiments has been vectorized and parallelized. This increases thespeed of the algorithm on the order of ten times within a standalone version onIntel architectures. This new version has been incorporated back intoLArSoftso that it can be generally used. These methods will also be applied to otherlow-level reconstruction algorithms of the wire signals such as the deconvolu-tion. The applications and performance of this modernized liquid argon TPCwire reconstruction will be presented.

Berkman, Sophie↗

DOE OSTI2019

Reconstruction for Liquid Argon TPC Neutrino Detectors Using Parallel Architectures [Slides]

Neutrinos are particles that interact rarely, so identifying them requires large detectors which produce lots of data. Processing this data with the computing power available is becoming more difficult as the detectors increase in size to reach their physics goals. In liquid argon time projection chambers (TPCs) the charged particles from neutrino interactions produce ionization electrons which drift in an electric field towards a series of collection wires, and the signal on the wires is used to reconstruct the interaction. The MicroBooNE detector currently collecting data at Fermilab has 8000 wires, and planned future experiments like DUNE will have 100 times more, which means that the time required to reconstruct an event will scale accordingly. Modernization of liquid argon TPC reconstruction code, including vectorization, parallelization and code portability to GPUs, will help to mitigate these challenges. The liquid argon TPC hit finding algorithm within the L arSoft framework used across multiple experiments has been vectorized and parallelized. This increases the speed of the algorithm on the order of ten times within a stand alone version on Intel architectures. This new version has been incorporated back into LarSoft so that it can be generally used. These methods will also be applied to other low level reconstruction of the wire signals such as the deconvolution. The applications and performance of this modernized liquid argon TPC wire reconstruction will be presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2021

Powering the Blue Economy: A Survey of Station-Keeping Methods for Mooringless Platforms

The term “ocean platform” is used to reference everything from stationary, typically moored, buoys to mobile water vehicles, whether they operate on the ocean’s surface or underwater. For certain applications for which relatively stationary station-keeping conditions are desired, the use of mooring systems is not always a viable alternative either for economic, environmental, regulatory, or otherwise practical reasons, or a combination thereof, (e.g., short deployments, sensitive ecosystems, very deep project sites). Maintaining a platform at a single waypoint or reference location without being moored would require additional control systems and a power source to counteract the drift forces that would naturally displace it. Mobile platforms, which are usually untethered except for remotely operated vehicles, typically require energy input to power their station-keeping capabilities so that they hold or control their location in the ocean. Currently, most of these platforms use combustion engines or batteries for this purpose, which, depending on the specific systems, may be costly, pollute the environment, or create limitations on the length of the deployment. However, powering this kind of platforms with surrounding renewable resources (waves, currents, winds, or sun) has been identified as a promising solution to expand their application. The intent of this report is to investigate station-keeping methods for various ocean platforms that are not moored or otherwise anchored to the ocean floor, or another platform or vessel, paying particular interest to technologies that use marine renewable resources to power their operation, because that is of particular interest to the U.S. Department of Energy’s Powering the Blue Economy (PBE) initiative. As a first step, 72 articles and technical reports related to mooringless station-keeping methods were collected for review. The preliminary literature review provided a broad overview of common themes across the literature from which a descriptive methodology for analyzing various platforms was developed. That is, station-keeping methods were categorized based on their predominant energy source and consumption (renewable, nonrenewable, or hybrid if the platform uses renewable and nonrenewable resources equally), and their localization strategy (drift reduction, “path-planning or “waypoint-holding”). In addition, platform types were segregated into the following groups: buoys, surface drifters, and unoccupied surface vehicles (USVs); offshore renewable energy systems; and unoccupied underwater vehicles (UUVs). The main types of station-keeping methods encountered in this report achieve their intended localization strategy by means of drift mitigation, steering, and/or propulsion. Drift mitigation is commonly accomplished via drogues and sea anchors. Stand-along steering subsystems use control surfaces (e.g., ship rudder, wing sail, etc.) that react to ocean currents, waves, or winds to provide varying-degrees of course adjustments. Combined steering and propulsion subsystems include differential thrusters, directional thrusters separate from a primary thruster that cause the platform to pitch up/down or yaw clockwise/counterclockwise, or vectored thrusters that direct the propulsion in a range of directions relative to the platform’s local coordinate system. Propulsion is often achieved by running a motor and applying active control strategies but can also involve buoyancy shifts and using sails to generate lifting forces that propel a platform in a desired direction. Future research is primarily expected to take place in the form of a technoeconomic analysis that would aim to determine the technological viability, cost, and added value of mooringless station-keeping use cases identified through this research, including docking for UUV recharging or for georeferencing drifter buoys, deep-sea floating wind farms, U.S. Navy sonar arrays, and a Pacific Ocean wave buoy network.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NASA NTRS1987

A survey of the cosmic ray diurnal variation during 1973-1979. I - Persistence of solar diurnal variation. II - Application of diffusion-convection model to diurnal anisotropy data

An analysis of data obtained with the vertical underground muon telescope at Embudo, NM shows that the solar diurnal variation in cosmic ray intensity is a persistent phenomenon over the 1973-1979 period. Assuming that the daily fluctuations in the amplitude and phase of the diurnal variation are random perturbations about the mean vector, the relative magnitude of the random component is determined. In the second part, the Diffusion-Convection model of cosmic ray transport is applied to high rigidity particles detected at the earth in order to deduce the behavior of the heliospheric transport parameters between 1973 and 1979. It is suggested that the diurnal variation observed at Embudo during 1979 may require a contribution from the charged particle drifts.

Riker, J. F.↗

DOE OSTI2013

Ion temperature gradient instability at sub-Larmor radius scales with non-zero ballooning angle

Linear gyro-kinetic stability calculations predict unstable toroidal ion temperature gradient modes (ITGs) with normalised poloidal wave vectors well above one (k{sub {theta}}{rho}{sub i}>1) for standard tokamak parameters with adiabatic electron response. These modes have a maximum amplitude at a poloidal angle {theta} that is shifted away from the low field side ({theta}{ne}0). The physical mechanism is clarified through the use of a fluid model. It is shown that the shift of the mode away from the low field side reduces the effective drift frequency which allows for the instability to develop. Numerical tests using the gyro-kinetic model confirm this physical mechanism. Furthermore, it is shown that modes localized away from the low field side can be important also for k{sub {theta}}{rho}{sub i}<1 close to the threshold of the ITG. In fact, modes with maximum amplitude at {theta}{ne}0 can exist for normalised temperature gradient lengths below the threshold of the ITG obtained for the case with the maximum at {theta}=0.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Implementation of a Doppler-Free Saturation Spectroscopy (DFSS) Diagnostic for Helicon Wave Electric Field Vector Measurement in Edge Plasma in DIII-D

A laser-based technique known as Doppler-free saturation spectroscopy (DFSS) has been designed, fabricated, and installed on the DIII-D National Fusion Facility to measure the helicon wave electric field vector in the edge plasma. These experimental measurements quantify phenomena resulting in decreased current drive efficiency due to wave/edge plasma interactions. This implementation of DFSS on DIII-D is the first of its kind on a tokamak and thus presents unique engineering challenges, including integration of the system onto an existing multidiagnostic port flange without impacts to system serviceability, as well as maintaining precise laser alignment over a 2-m distance during disruptions and thermal drift of the vessel. Further, these challenges were resolved using innovative design approaches such as a novel decoupled shutter system to facilitate serviceability of the in-vessel mirror assemblies without the need for personnel vessel entry, as well as an ex-vessel piezo-mirror-based optical system for laser beam shaping and real-time steering of the measurement location. The solutions to these engineering challenges were demonstrated during the successful installation and operation of these diagnostic components during the 2022 DIII-D vent and subsequent experimental campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1993

Voyager planetary radio astronomy studies

Analysis of nonthermal radio emission data obtained by the Planetary Radio Astronomy (PRA) spectrometers on the Voyager 1 and 2 spacecraft was performed. This PRA data provided unique insights into the radio emission characteristics of the outer planets because of PRA's unique spectral response below the terrestrial ionospheric plasma frequency and its unprecedented proximity to the source. Of those results which were documented or published, this final report surveys only the highlights and cites references for more complete discussions. Unpublished results for Uranus, Neptune, and theoretical Ionian current distributions are presented at greater length. The most important conclusion to be drawn from these observations is that banded spectral emission is common to the radio emission below 1-2 MHz observed from all four Jovian planets. In every case multiple spectral features evolve on time scales of seconds to minutes. To the extent these features drift in frequency, they appear never to cross one another. The Neptunian spectral features appear to drift little or not at all, their evolution consisting principally of waxing and waning. Since other evidence strongly suggests that most or all of this radio emission is occurring near the local magnetospheric electron cyclotron frequency, this implies that this emission preferentially occurs at certain continually changing planetary radii. It remains unknown why certain radii might be favored, unless radial electric field components or other means serve to differentiate radially the magnetospheric plasma density, particle energy vectors, or particle coherence. Calculation of the spatial distribution and intensity of the Io-generated magnetospheric currents are also presented; these currents may be limited principally by wave impedance and local field strengths.

Staelin, David H.↗

DOE OSTI2014

SU-C-18A-04: 3D Markerless Registration of Lung Based On Coherent Point Drift: Application in Image Guided Radiotherapy

Purpose: This study evaluated a new probabilistic non-rigid registration method called coherent point drift for real time 3D markerless registration of the lung motion during radiotherapy. Method: 4DCT image datasets Dir-lab (www.dir-lab.com) have been used for creating 3D boundary element model of the lungs. For the first step, the 3D surface of the lungs in respiration phases T0 and T50 were segmented and divided into a finite number of linear triangular elements. Each triangle is a two dimensional object which has three vertices (each vertex has three degree of freedom). One of the main features of the lungs motion is velocity coherence so the vertices that creating the mesh of the lungs should also have features and degree of freedom of lung structure. This means that the vertices close to each other tend to move coherently. In the next step, we implemented a probabilistic non-rigid registration method called coherent point drift to calculate nonlinear displacement of vertices between different expiratory phases. Results: The method has been applied to images of 10-patients in Dir-lab dataset. The normal distribution of vertices to the origin for each expiratory stage were calculated. The results shows that the maximum error of registration between different expiratory phases is less than 0.4 mm (0.38 SI, 0.33 mm AP, 0.29 mm RL direction). This method is a reliable method for calculating the vector of displacement, and the degrees of freedom (DOFs) of lung structure in radiotherapy. Conclusions: We evaluated a new 3D registration method for distribution set of vertices inside lungs mesh. In this technique, lungs motion considering velocity coherence are inserted as a penalty in regularization function. The results indicate that high registration accuracy is achievable with CPD. This method is helpful for calculating of displacement vector and analyzing possible physiological and anatomical changes during treatment.

60 APPLIED LIFE SCIENCES↗