Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “signal processing algorithms”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Ultra-Stable Beacon Source for Laboratory Testing of Optical Tracking

The ultra-stable beacon source (USBS) provides a laser-beam output with a very low angular jitter and can be used as an absolute angular reference to simulate a beacon in the laboratory. The laser is mounted on the top of a very short (approximately equal to 1 m) inverted pendulum (IP) with its optical axis parallel to the carbon fiber pendulum leg. The 85-cm, carbon fiber rods making up the leg are very lightweight and rigid, and are supported by a flex-joint at the bottom (see figure). The gimbal-mounted laser is a weight-adjustable load of about 1.5 kg with its center of rotation co-located with the center of percussion of the inverted pendulum. This reduces the coupling of transverse motion at the base of the pendulum to angular motion of the laser at the top. The inverted pendulum is mounted on a gimbal with its center of rotation coinciding with the pivot position of the inverted pendulum flexure joint. This reduces coupling of ground tilt at the inverted pendulum base to motion of the laser mounted at the top. The mass of the top gimbal is adjusted to give the pendulum a very low resonant frequency (approximately equal to 10 mHz) that filters transverse seismic disturbances from the ground where the base is attached. The motion of the IP is monitored by an optical-lever sensor. The laser light is reflected by the mirror on the IP, and then is detected by a quadrant photo-detector (QPD). The position of the beam spot on the QPD corresponds to the tilt of the IP. Damping of this motion is provided by two coil and magnet pairs. The bottom gimbal mount consists of two plates. The IP is mounted on the second plate. The first plate is supported by two posts through needles and can be rotated about the axis connecting the tips of the needles. The second plate hangs from the first plate and can be rotated about the axis perpendicular to the first plate. As a result, the second plate acts as a two-axis rotation stage. Its center of rotation is located at the effective bending point of the flex-joint. The second plate is pressed against two screw actuators by the weight of the IP. The screw actuators are orthogonal to each other and are used to adjust the inclination of the second plate. The actuators are driven by stepper motors. The whole IP system is housed in a box made of Lexan plastic plates to provide isolation from air currents and temperature variations. The signals from the sensors are processed and recorded with a PC using the xPC Target realtime environment of Math- Works. The control algorithms are written using the Simulink package from The MathWorks.

Aso, Yoichi↗
NASA NTRS1979

A model for testing centerfinding algorithms for automated optical navigation

An efficient software simulation of the imaging process for optical navigation is presented, illustrating results using simple examples. The problems of image definition and optical system modeling, including ideal image containing features and realistic models of optical filtering performed by the entire camera system, are examined. A digital signal processing technique is applied to the problem of developing methods of automated optical navigation and the subsequent mathematical formulation is presented. Specific objectives such as an analysis of the effects of camera defocusing on centerfinding of planar targets, addition of noise filtering to the algorithm, and implementation of multiple frame capability were investigated.

Griffin, M. D.↗
NASA NTRS2004

Resolving the Issues with Flywheel Position Sensors

For the past few years, the Advanced Electrical Systems Branch here at NASA Glenn has been pursuing research in the area of flywheels. The purpose of these pursuits has been t o explore the potential for flywheels to replace current battery-powered systems in space. So far it has been learned that flywheels offer large momentum storage capacity, comparatively small volume, high durability, and near- complete discharge capabilities, all of which are advancements over the existing nickel hydrogen and nickel cadmium batteries. Another significant advantage of flywheels is the potential they offer for combining the function of attitude control with energy storage. During the summer of 2004, I worked with Dr. Barbara Kenny in the Advanced Electrical Systems Branch, supporting the work she is doing by analyzing and testing some new components for the new Generation-2 flywheel. To monitor the speed and angular position of the flywheel rotor, a once-around (OAR) signal along with a sensorless algorithm is used. The OAR signal is used for the magnetic bearings that keep the flywheel suspended for frictionless operation. The sensorless algorithm is used for the flywheel motor/generator control. The OAR is generated from position sensors that monitor a circular plate. The plate has a cut down the middle such that one half of the circle is on a slightly lower level than the other. Every half-turn, or 180, the sensors detect the "cut" on the plate, and trigger the OAR, telling the computer that the rotor has made half a revolution. This, however, doesn't provide needed detailed information about the angular position of the rotor, since it only provides a signal alert every half- revolution. This is enough information for the magnetic bearing control but is insufficient for the motor/generator control. A new resolver was designed such that it would give continuous angle information rather than the 180 degree information of the OAR. The new resolver has two separate observable pieces: a flat middle section to monitor vertical motion, and an angled section around the circumference, which, when observed from above, produces a sine-wave displacement through the entire 360" revolution. My first job when I arrived this summer was to calibrate the sensors that would be mounted on the inside of the flywheel casing to monitor the position (angular and vertical) of the shaft. After calibration, I used the sensors to evaluate voltage outputs created by position differences between two pairs of sensors on the angled portion of the resolver for eight different angular positions, moving the resolver vertically and laterally through its entire potential range of motion. The results of these tests will be used to determine the rotor angular (and axial) position from the sensor readings once the new flywheel unit is assembled. The sensorless algorighm mentioned above consists of two operations: the signal injection method and the back electro-motive force (EMF). The signal injection is meant to work at low speeds, while the back EMF algorithm is meant to work at higher speeds. Both work together to determine the correct estimate of rotor position and speed based on the measured motor/generator current. It was determined that we wanted to know exactly how accurate our estimation methods were, and so a resolver (a commercially available mechanical sensor mounted to the motor/generator shaft to measure rotor position and speed) and a "Resolver to Digital" (R2D) circuit board was purchased to make the comparison to the existing estimation. My work related to the R2D board has included the following: creating two connector cables (one to power the circuit and one to get readable output off the board), writing Simulink code to process the board's output, and building a dSpace panel to control and monitor the circuit. The next step in the process will be to perform tests to compare the estimated rotor position and speed from the sensorless algorithm to the actual rotor and speed from the resolver signal.

Fehrmann, Elizabeth A.↗
NASA NTRS2004

Software for Generating Strip Maps from SAR Data

Jurassicprok is a computer program that generates strip-map digital elevation models and other data products from raw data acquired by an airborne synthetic-aperture radar (SAR) system. This software can process data from a variety of airborne SAR systems but is designed especially for the GeoSAR system, which is a dual-frequency (P- and X-band), single-pass interferometric SAR system for measuring elevation both at the bare ground surface and top of the vegetation canopy. Jurassicprok is a modified version of software developed previously for airborne-interferometric- SAR applications. The modifications were made to accommodate P-band interferometric processing, remove approximations that are not generally valid, and reduce processor-induced mapping errors to the centimeter level. Major additions and other improvements over the prior software include the following: a) A new, highly efficient multi-stage-modified wave-domain processing algorithm for accurately motion compensating ultra-wideband data; b) Adaptive regridding algorithms based on estimated noise and actual measured topography to reduce noise while maintaining spatial resolution; c) Exact expressions for height determination from interferogram data; d) Fully calibrated volumetric correlation data based on rigorous removal of geometric and signal-to-noise decorrelation terms; e) Strip range-Doppler image output in user-specified Doppler coordinates; f) An improved phase-unwrapping and absolute-phase-determination algorithm; g) A more flexible user interface with many additional processing options; h) Increased interferogram filtering options; and i) Ability to use disk space instead of random- access memory for some processing steps.

Hensley, Scott↗
DOE OSTI2016

Applying the new gamma ray imager diagnostic to measurements of runaway electron Bremsstrahlung radiation in the DIII-D Tokamak (invited)

A new gamma ray imager (GRI) is developed to probe the electron distribution function with 2D spatial resolution during runaway electron (RE) experiments at the DIII-D tokamak. The diagnostic is sensitive to 0.5–100 MeV gamma rays, allowing characterization of the RE distribution function evolution during RE growth and dissipation. The GRI consists of a lead “pinhole camera” mounted on the DIII-D midplane with 123 honeycombed tangential chords 20 cm wide that span the vessel interior. Up to 30 bismuth germanate (BGO) scintillation detectors capture RE bremsstrahlung radiation for Pulse Height Analysis (PHA) capable of discriminating up to 20 000 pulses per second. Digital signal processing routines combining shaping filters are performed during PHA to reject noise and record gamma ray energy. The GRI setup and PHA algorithms will be described and initial data from experiments will be presented. A synthetic diagnostic is developed to generate the gamma ray spectrum of a GRI channel given the plasma information and a prescribed distribution function. Magnetic reconstructions of the plasma are used to calculate the angle between every GRI sightline and orient and discriminate gamma rays emitted by a field-aligned RE distribution function.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

The Octopus processor for the CMS L1 muon trigger for High Luminosity LHC

The upgraded L1 muon trigger system of the CMS experiment in the High Luminosity Large Hadron Collider is based on custom processors featuring large Field Programmable Gate Arrays (FPGAs) connected by large numbers of optical links. These provide the I/O bandwidth and power necessary to process the complex algorithms used during the collection of physics data. The design and performance requirements of these processors creates significant challenges in signal integrity, power delivery, and thermal management. In this paper we describe the Octopus processor, featuring a large Xilinx Virtex Ultrascale+ FPGA and up to 128 links interfaced to optics through high quality twin-ax copper cables. Results on signal integrity at 25 Gb/s and the first demonstration of 50+ Gb/s links with pluggable optics in CMS are also shown, demonstrating bit error rates below 10 –15 at a 95% confidence level. The thermal performance is measured inside an Advanced-TCA crate with acceptable thermal margins up to 200 W of chip power. Future improvements are mentioned, potentially allowing operation at up to 300 W.

Instruments & Instrumentation↗
DOE OSTI2022

Magnetic-Field Diffusion Effects in Beam Position Monitors II: Application to Calibration Single-Pulse Data

Beam position monitors (BPMs) provide time-resolved measurements of the current and centroid position of high-current electron beams in linear induction accelerators (LIAs). One of the types of detectors used in BPMs is the B-dot loop, which generates a signal from the EMF due to the time varying magnetic flux through the loop. If some of the boundaries of the loop are composed of thick metal walls with finite conductivity, the resulting signal must be corrected for the magnetic field diffusion into the metal. From first principles, we have derived the perturbation to BPM measurements due to this effect. The theoretical framework was used to design an algorithm for signal correction that does not require knowledge of the time history of the magnetic flux being measured. Corrected signals based on that process compared favorably with the know reference signals in a laboratory calibration test sequence.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Progress Toward the First Search for Bound Neutron Oscillation into Antineutron in a Liquid Argon TPC

This note presents current progress for a neutron-antineutron oscillation ($n–\overline{n}$) search in MicroBooNE paving the way for the first search analysis of such process in a Liquid Argon Time Projection Chamber (LArTPC). Convolutional Neural Network (CNN) and Boosted Decision Tree (BDT) algorithms were used to select signal $n–\overline{n}$ events over cosmogenic backgrounds. The CNN-only, BDT-only, and the combined (CNN+BDT) methods were demonstrated on the Monte-Carlo signal and background events. Validation of the CNNonly and the BDT-only methods was carried out on a small dataset of MicroBooNE Run1 off-beam data, setting the starting point toward further improvement of the analysis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Multivariable quantum signal processing (M-QSP): prophecies of the two-headed oracle

Recent work shows that quantum signal processing (QSP) and its multi-qubit lifted version, quantum singular value transformation (QSVT), unify and improve the presentation of most quantum algorithms. QSP/QSVT characterize the ability, by alternating ansätze, to obliviously transform the singular values of subsystems of unitary matrices by polynomial functions; these algorithms are numerically stable and analytically well-understood. That said, QSP/QSVT require consistent access to a single oracle, saying nothing about computing joint properties of two or more oracles; these can be far cheaper to determine given an ability to pit oracles against one another coherently. This work introduces a corresponding theory of QSP over multiple variables: M-QSP. Surprisingly, despite the non-existence of the fundamental theorem of algebra for multivariable polynomials, there exist necessary and sufficient conditions under which a desired stable multivariable polynomial transformation is possible. Moreover, the classical subroutines used by QSP protocols survive in the multivariable setting for non-obvious reasons, and remain numerically stable and efficient. Up to a well-defined conjecture, we give proof that the family of achievable multivariable transforms is as loosely constrained as could be expected. The unique ability of M-QSP to obliviously approximate joint functions of multiple variables coherently leads to novel speedups incommensurate with those of other quantum algorithms, and provides a bridge from quantum algorithms to algebraic geometry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1984

Eliminating "Hotspots" in Digital Image Processing

Signals from defective picture elements rejected. Image processing program for use with charge-coupled device (CCD) or other mosaic imager augmented with algorithm that compensates for common type of electronic defect. Algorithm prevents false interpretation of "hotspots". Used for robotics, image enhancement, image analysis and digital television.

Salomon, P. M.↗
NASA NTRS2003

Noninvasive identification of the total peripheral resistance baroreflex

We propose two identification algorithms for quantitating the total peripheral resistance (TPR) baroreflex, an important contributor to short-term arterial blood pressure (ABP) regulation. Each algorithm analyzes beat-to-beat fluctuations in ABP and cardiac output, which may both be obtained noninvasively in humans. For a theoretical evaluation, we applied both algorithms to a realistic cardiovascular model. The results contrasted with only one of the algorithms proving to be reliable. This algorithm was able to track changes in the static gains of both the arterial and cardiopulmonary TPR baroreflex. We then applied both algorithms to a preliminary set of human data and obtained contrasting results much like those obtained from the cardiovascular model, thereby making the theoretical evaluation results more meaningful. This study suggests that, with experimental testing, the reliable identification algorithm may provide a powerful, noninvasive means for quantitating the TPR baroreflex. This study also provides an example of the role that models can play in the development and initial evaluation of algorithms aimed at quantitating important physiological mechanisms.

Non-NASA Center↗
DOE OSTI2017

Adaptive OFDM Waveform Design for Spatio-Temporal-Sparsity Exploited STAP Radar

In this chapter, we describe a sparsity-based space-time adaptive processing (STAP) algorithm to detect a slowly moving target using an orthogonal frequency division multiplexing (OFDM) radar. The motivation of employing an OFDM signal is that it improves the target-detectability from the interfering signals by increasing the frequency diversity of the system. However, due to the addition of one extra dimension in terms of frequency, the adaptive degrees-of-freedom in an OFDM-STAP also increases. Therefore, to avoid the construction a fully adaptive OFDM-STAP, we develop a sparsity-based STAP algorithm. We observe that the interference spectrum is inherently sparse in the spatio-temporal domain, as the clutter responses occupy only a diagonal ridge on the spatio-temporal plane and the jammer signals interfere only from a few spatial directions. Hence, we exploit that sparsity to develop an efficient STAP technique that utilizes considerably lesser number of secondary data compared to the other existing STAP techniques, and produces nearly optimum STAP performance. In addition to designing the STAP filter, we optimally design the transmit OFDM signals by maximizing the output signal-to-interference-plus-noise ratio (SINR) in order to improve the STAP performance. The computation of output SINR depends on the estimated value of the interference covariance matrix, which we obtain by applying the sparse recovery algorithm. Therefore, we analytically assess the effects of the synthesized OFDM coefficients on the sparse recovery of the interference covariance matrix by computing the coherence measure of the sparse measurement matrix. Our numerical examples demonstrate the achieved STAP-performance due to sparsity-based technique and adaptive waveform design.

Sen, Satyabrata↗
DOE OSTI2015

Measurement of the Fast Neutron Response for {sup 4}He Scintillation Detectors Using a Coincidence Scattering Method

Previous work has measured the neutron response of pressurized {sup 4}He scintillation detectors, however these studies only examine the response as a function of incident neutron energy. Since the detection mechanism in {sup 4}He detectors is elastic scattering, and the interacting neutron will only deposit a fraction of its incident kinetic energy in the detector gas, an examination of the response of the detector output to deposited energy is necessary to transform these detectors into instruments for neutron spectrometry. Using a combined time-of-flight (TOF) and coincidence scattering method, this paper further characterizes the {sup 4}He light response to fast neutrons by examining the scintillation light yield as a function of deposited energy, measuring the light response up to 5 MeV. These {sup 4}He detectors are simple in design, and are manufactured by Arktis Radiation Detectors in several sizes. The specific model used in this experiment had an active volume 20 cm long with an inner diameter of 4.4 cm, giving a total active volume of 304 cm{sup 3}. The key components include the active volume, filled with 150 bar of helium-4 gas, and photomultiplier tubes (PMTs) mounted at either end of the active volume. The detector body is made of stainless steel. The detector response was experimentally measured using a two-detector coincidence arrangement with a {sup 252}Cf source. Two {sup 4}He detectors were vertically mounted, and the source was placed at a horizontal distance from the center of the bottom detector, forming a right angle. By requiring coincidence between the two detectors, it was confirmed that each neutron interacting in the second (top) detector must first have undergone a scattering interaction in the first (bottom) detector, and the time-of-flight (TOF) technique could then be used to determine the energy of the neutron as it traveled between the two detectors by the difference in time between the two detector events. More importantly, with the scattering angle known, the amount of energy deposited by the neutron in the bottom detector (ER) was also calculated using kinematic scattering equations. This deposited recoil energy was then compared to the corresponding light output for each event to form a deposited energy scintillation light response matrix. Similarly, the system's insensitivity to gammas and its ability to reject gammas by pulse shape discrimination (PSD) are often cited as an important advantage, although a detailed analysis of these capabilities has not yet been performed. This work therefore quantified these parameters in order to further characterize these detectors for future mixed radiation field measurements. Gamma sources were measured spanning a range of gamma-ray energies from 0.122 MeV to 1.332 MeV, including {sup 57}Co, {sup 137}Cs, {sup 54}Mn, and {sup 60}Co. Each source was counted by the {sup 4}He detector and the background subtracted. Taking the ratio of the number of events detected during the experimental source measurement to the number of gammas predicted by MCNPX to pass through the detector volume yields the detector's intrinsic gamma efficiency. The difference between this fraction and unity is therefore a measure of the detector's ability to ignore interfering gamma rays, defined as its inherent gamma rejection rate. The ability of post-processing PSD algorithms to further reduce the number of gammas is also investigated and quantified. Finally, it has been noted that the scintillation signal from a single neutron event can be separated in time into two components: the fast component is a sharp peak that exists on the order of nanoseconds; the slow component is a series of smaller pulses, stretched out over four microseconds. Whereas previous research has exclusively focused on the energy information contained in the slow component, this work demonstrates that the fast component is also sensitive to neutron energy, and the entire scintillation signal can therefore be used. In conclusion, the relationship of fast neutron {sup 4}He scintillation detectors to deposited neutron energy was explored, and will be combined with previous works that measured the scintillation response to incident neutron energy in order to develop a neutron spectrometer. Similarly, the ability of these {sup 4}He detectors to reject interfering gamma rays was also quantified, and so will enable this spectrometer to be deployed in mixed radiation field measurements. Finally, while previous works with these detectors have focused on an analysis of the slow scintillation component, it was demonstrated in this work that the fast component also contains significant energy information.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2001

Electronic Thermometer Readings

NASA Stennis' adaptive predictive algorithm for electronic thermometers uses sample readings during the initial rise in temperature and applies an algorithm that accurately and rapidly predicts the steady state temperature. The final steady state temperature of an object can be calculated based on the second-order logarithm of the temperature signals acquired by the sensor and predetermined variables from the sensor characteristics. These variables are calculated during tests of the sensor. Once the variables are determined, relatively little data acquisition and data processing time by the algorithm is required to provide a near-accurate approximation of the final temperature. This reduces the delay in the steady state response time of a temperature sensor. This advanced algorithm can be implemented in existing software or hardware with an erasable programmable read-only memory (EPROM). The capability for easy integration eliminates the expense of developing a whole new system that offers the benefits provided by NASA Stennis' technology.

Source record↗
DOE OSTI2015

Sparsity-based acoustic inversion in cross-sectional multiscale optoacoustic imaging

Purpose: With recent advancement in hardware of optoacoustic imaging systems, highly detailed cross-sectional images may be acquired at a single laser shot, thus eliminating motion artifacts. Nonetheless, other sources of artifacts remain due to signal distortion or out-of-plane signals. The purpose of image reconstruction algorithms is to obtain the most accurate images from noisy, distorted projection data. Methods: In this paper, the authors use the model-based approach for acoustic inversion, combined with a sparsity-based inversion procedure. Specifically, a cost function is used that includes the L1 norm of the image in sparse representation and a total variation (TV) term. The optimization problem is solved by a numerically efficient implementation of a nonlinear gradient descent algorithm. TV–L1 model-based inversion is tested in the cross section geometry for numerically generated data as well as for in vivo experimental data from an adult mouse. Results: In all cases, model-based TV–L1 inversion showed a better performance over the conventional Tikhonov regularization, TV inversion, and L1 inversion. In the numerical examples, the images reconstructed with TV–L1 inversion were quantitatively more similar to the originating images. In the experimental examples, TV–L1 inversion yielded sharper images and weaker streak artifact. Conclusions: The results herein show that TV–L1 inversion is capable of improving the quality of highly detailed, multiscale optoacoustic images obtained in vivo using cross-sectional imaging systems. As a result of its high fidelity, model-based TV–L1 inversion may be considered as the new standard for image reconstruction in cross-sectional imaging.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2002

Cloud Algorithm Design and Performance for the 2002 Geoscience Laser Altimeter System Mission

A satellite borne lidar instrument, the Geoscience Laser Altimeter System (GLAS), is to be launched in late 2002 and will provide continuous profiling of atmospheric clouds and aerosol on a global basis. Data processing algorithms have been developed to provide operational data products in near real time. Basic data products for cloud observations are the height of the top and bottom of single to multiple cloud layers and the lidar calibrated observed backscatter cross section up to the level of signal attenuation. In addition the optical depth and vertical profile of visible extinction cross section of many transmissive cloud layers and most haze layers are to be derived. The optical thickness is derivable in some cases from the attenuation of the molecular scattering below cloud base. In other cases an assumption of the scattering phase function is required. In both cases a estimated correction for multiple scattering is required. The data processing algorithms have been tested in part from aircraft measurements used to simulated satellite data. The GLAS lidar observations will be made from an orbit that will allow inter comparison with all other existing satellite cloud measurements.

Spinhirne, J. D.↗
DOE OSTI2016

SU-F-I-10: Spatially Local Statistics for Adaptive Image Filtering

Purpose: To facilitate adaptive image filtering operations, addressing spatial variations in both noise and signal. Such issues are prevalent in cone-beam projections, where physical effects such as X-ray scattering result in spatially variant noise, violating common assumptions of homogeneous noise and challenging conventional filtering approaches to signal extraction and noise suppression. Methods: We present a computational mechanism for probing into and quantifying the spatial variance of noise throughout an image. The mechanism builds a pyramid of local statistics at multiple spatial scales; local statistical information at each scale includes (weighted) mean, median, standard deviation, median absolute deviation, as well as histogram or dynamic range after local mean/median shifting. Based on inter-scale differences of local statistics, the spatial scope of distinguishable noise variation is detected in a semi- or un-supervised manner. Additionally, we propose and demonstrate the incorporation of such information in globally parametrized (i.e., non-adaptive) filters, effectively transforming the latter into spatially adaptive filters. The multi-scale mechanism is materialized by efficient algorithms and implemented in parallel CPU/GPU architectures. Results: We demonstrate the impact of local statistics for adaptive image processing and analysis using cone-beam projections of a Catphan phantom, fitted within an annulus to increase X-ray scattering. The effective spatial scope of local statistics calculations is shown to vary throughout the image domain, necessitating multi-scale noise and signal structure analysis. Filtering results with and without spatial filter adaptation are compared visually, illustrating improvements in imaging signal extraction and noise suppression, and in preserving information in low-contrast regions. Conclusion: Local image statistics can be incorporated in filtering operations to equip them with spatial adaptivity to spatial signal/noise variations. An efficient multi-scale computational mechanism is developed to curtail processing latency. Spatially adaptive filtering may impact subsequent processing tasks such as reconstruction and numerical gradient computations for deformable registration. NIH Grant No. R01-184173.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

WE-E-17A-06: Assessing the Scale of Tumor Heterogeneity by Complete Hierarchical Segmentation On MRI

Purpose: In many cancers, intratumoral heterogeneity exists in vascular and genetic structure. We developed an algorithm which uses clinical imaging to interrogate different scales of heterogeneity. We hypothesize that heterogeneity of perfusion at large distance scales may correlate with propensity for disease recurrence. We applied the algorithm to initial diagnosis MRI of rhabdomyosarcoma patients to predict recurrence. Methods: The Spatial Heterogeneity Analysis by Recursive Partitioning (SHARP) algorithm recursively segments the tumor image. The tumor is repeatedly subdivided, with each dividing line chosen to maximize signal intensity difference between the two subregions. This process continues to the voxel level, producing segments at multiple scales. Heterogeneity is measured by comparing signal intensity histograms between each segmented region and the adjacent region. We measured the scales of contrast enhancement heterogeneity of the primary tumor in 18 rhabdomyosarcoma patients. Using Cox proportional hazards regression, we explored the influence of heterogeneity parameters on relapse-free survival (RFS). To compare with existing methods, fractal and Haralick texture features were also calculated. Results: The complete segmentation produced by SHARP allows extraction of diverse features, including the amount of heterogeneity at various distance scales, the area of the tumor with the most heterogeneity at each scale, and for a given point in the tumor, the heterogeneity at different scales. 10/18 rhabdomyosarcoma patients suffered disease recurrence. On contrast-enhanced MRI, larger scale of maximum signal intensity heterogeneity, relative to tumor diameter, predicted for shorter RFS (p=0.05). Fractal dimension, fractal fit, and three Haralick features did not predict RFS (p=0.09-0.90). Conclusion: SHARP produces an automatic segmentation of tumor regions and reports the amount of heterogeneity at various distance scales. In rhabdomyosarcoma, RFS was shorter when the primary tumor exhibited larger scale of heterogeneity on contrast-enhanced MRI. If validated on a larger dataset, this imaging biomarker could be useful to help personalize treatment.

60 APPLIED LIFE SCIENCES↗