Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Image processing”

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

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

At least 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Graphics Processing Unit Assisted Thermographic Compositing

Objective: To develop a software application utilizing general purpose graphics processing units (GPUs) for the analysis of large sets of thermographic data. Background: Over the past few years, an increasing effort among scientists and engineers to utilize the GPU in a more general purpose fashion is allowing for supercomputer level results at individual workstations. As data sets grow, the methods to work them grow at an equal, and often greater, pace. Certain common computations can take advantage of the massively parallel and optimized hardware constructs of the GPU to allow for throughput that was previously reserved for compute clusters. These common computations have high degrees of data parallelism, that is, they are the same computation applied to a large set of data where the result does not depend on other data elements. Signal (image) processing is one area were GPUs are being used to greatly increase the performance of certain algorithms and analysis techniques.

Ragasa, Scott↗
NASA NTRS2002

On the nature and evolution of the neural bases of human language

The traditional theory equating the brain bases of language with Broca's and Wernicke's neocortical areas is wrong. Neural circuits linking activity in anatomically segregated populations of neurons in subcortical structures and the neocortex throughout the human brain regulate complex behaviors such as walking, talking, and comprehending the meaning of sentences. When we hear or read a word, neural structures involved in the perception or real-world associations of the word are activated as well as posterior cortical regions adjacent to Wernicke's area. Many areas of the neocortex and subcortical structures support the cortical-striatal-cortical circuits that confer complex syntactic ability, speech production, and a large vocabulary. However, many of these structures also form part of the neural circuits regulating other aspects of behavior. For example, the basal ganglia, which regulate motor control, are also crucial elements in the circuits that confer human linguistic ability and abstract reasoning. The cerebellum, traditionally associated with motor control, is active in motor learning. The basal ganglia are also key elements in reward-based learning. Data from studies of Broca's aphasia, Parkinson's disease, hypoxia, focal brain damage, and a genetically transmitted brain anomaly (the putative "language gene," family KE), and from comparative studies of the brains and behavior of other species, demonstrate that the basal ganglia sequence the discrete elements that constitute a complete motor act, syntactic process, or thought process. Imaging studies of intact human subjects and electrophysiologic and tracer studies of the brains and behavior of other species confirm these findings. As Dobzansky put it, "Nothing in biology makes sense except in the light of evolution" (cited in Mayr, 1982). That applies with as much force to the human brain and the neural bases of language as it does to the human foot or jaw. The converse follows: the mark of evolution on the brains of human beings and other species provides insight into the evolution of the brain bases of human language. The neural substrate that regulated motor control in the common ancestor of apes and humans most likely was modified to enhance cognitive and linguistic ability. Speech communication played a central role in this process. However, the process that ultimately resulted in the human brain may have started when our earliest hominid ancestors began to walk.

Review, Tutorial↗
NASA NTRS1994

Interactions between nested sunspots. 1: The formation and breakup of a delta-type sunspot

We investigate a nest of sunspots in which three ordinary bipolar pairs of sunspots are aligned collinearly. The usual spreading action of the growing regions brings two spots of leading polarity together (p-p collision) and forces the leading and trailing spots of the two interior regions to overlap inot a single penumbra (p-f collision), thus forming a delta-spot. We examine digitally processed images from the Ottawa River Solar Observatory of two related events inside the delta-spot 5 days after the p-f collision begins: the violent disruption of the f-umbra, and the formation in less than a day of an hydrogen-alpha filament. The evolutionary changes in shape, area, relative motions, and brightness that we measure for each spot in the elongated nest are more compatible with Parker's (1979a) hypothesis of a sunspot as a cluster of flux tubes held together by downdrafts than with the notion of a sunspot as a monolithic plug of magnetic flux. From chromospheric developments over the delta-spot, we show that a shearing motion along a polarity inversion is more effective than convergence for creating a chromospheric filament. We invoke the release of an instability, triggered by a sequence of processes lasting 1 day or more, to explain the disruption of the f-umbra in this delta-spot. We show that the sequence is initiated when the colliding p-f umbrae reach a critical separation around 3200 +/- 200 km. We present a descriptive model in which the reconnected magnetic fields block vertical transport of convective heat flux just beneath the photosphere. We observe the formation of an unusual type of penumbra adjacent to the f-polarity portion of this delta-spot just before its disruption. A tangential penumbral band grows out of disordered matter connected to the f-umbra. We present this as evidence for the extrusion of umbral magnetic flux by thermal plumes rising through a loosely bound umbra.

Gaizauskas, V.↗
NASA NTRS2005

A Fast Implementation of the ISODATA Clustering Algorithm

Clustering is central to many image processing and remote sensing applications. ISODATA is one of the most popular and widely used clustering methods in geoscience applications, but it can run slowly, particularly with large data sets. We present a more efficient approach to ISODATA clustering, which achieves better running times by storing the points in a kd-tree and through a modification of the way in which the algorithm estimates the dispersion of each cluster. We also present an approximate version of the algorithm which allows the user to further improve the running time, at the expense of lower fidelity in computing the nearest cluster center to each point. We provide both theoretical and empirical justification that our modified approach produces clusterings that are very similar to those produced by the standard ISODATA approach. We also provide empirical studies on both synthetic data and remotely sensed Landsat and MODIS images that show that our approach has significantly lower running times.

Memarsadeghi, Nargess↗
NASA NTRS2007

A Fast Implementation of the Isodata Clustering Algorithm

Clustering is central to many image processing and remote sensing applications. ISODATA is one of the most popular and widely used clustering methods in geoscience applications, but it can run slowly, particularly with large data sets. We present a more efficient approach to IsoDATA clustering, which achieves better running times by storing the points in a kd-tree and through a modification of the way in which the algorithm estimates the dispersion of each cluster. We also present an approximate version of the algorithm which allows the user to further improve the running time, at the expense of lower fidelity in computing the nearest cluster center to each point. We provide both theoretical and empirical justification that our modified approach produces clusterings that are very similar to those produced by the standard ISODATA approach. We also provide empirical studies on both synthetic data and remotely sensed Landsat and MODIS images that show that our approach has significantly lower running times.

Memarsadeghi, Nargess↗
NASA NTRS2019

Overview of the New Version 3 NASA Micro-Pulse Lidar Network (MPLNET) Automatic Precipitation Detection Algorithm

Precipitation modifies atmospheric column thermodynamics through the process of evaporation and serves as a proxy for latent heat modulation. For this reason, a correct precipitation parameterization (especially for low-intensity precipitation) within global scale models is crucial. In addition to improving our modeling of the hydrological cycle, this will reduce the associated uncertainty of global climate models in correctly forecasting future scenarios, and will enable the application of mitigation strategies. In this manuscript we present a proof of concept algorithm to automatically detect precipitation from lidar measurements obtained from the National Aeronautics and Space Administration Micropulse lidar network (MPLNET). The algorithm, once tested and validated against other remote sensing instruments, will be operationally implemented into the network to deliver a near real time (latency <1.5 h) rain masking variable that will be publicly available on MPLNET website as part of the new Version 3 data products. The methodology, based on an image processing technique, detects only light precipitation events (defined by intensity and duration) such as light rain, drizzle, and virga. During heavy rain events, the lidar signal is completely extinguished after a few meters in the precipitation or it is unusable because of water accumulated on the receiver optics. Results from the algorithm, in addition to filling a gap in light rain, drizzle, and virga detection by radars, are of particular interest for the scientific community as they help to fully characterize the aerosol cycle, from emission to deposition, as precipitation is a crucial meteorological phenomenon accelerating atmospheric aerosol removal through the scavenging effect. Algorithm results will also help the understanding of long term aerosol–cloud interactions, exploiting the multi-year database from several MPLNET permanent observational sites across the globe. The algorithm is also applicable to other lidar and/or ceilometer network infrastructures in the framework of the Global Aerosol Watch (GAW) aerosol lidar observation network (GALION).

Simone Lolli↗
NASA NTRS2003

Structural Anomalies Detected in Ceramic Matrix Composites Using Combined Nondestructive Evaluation and Finite Element Analysis (NDE and FEA)

Most reverse engineering approaches involve imaging or digitizing an object and then creating a computerized reconstruction that can be integrated, in three dimensions, into a particular design environment. The rapid prototyping technique builds high-quality physical prototypes directly from computer-aided design files. This fundamental technique for interpreting and interacting with large data sets is being used here via Velocity2 (an integrated image-processing software, ref. 1) using computed tomography (CT) data to produce a prototype three-dimensional test specimen model for analyses. A study at the NASA Glenn Research Center proposes to use these capabilities to conduct a combined nondestructive evaluation (NDE) and finite element analysis (FEA) to screen pretest and posttest structural anomalies in structural components. A tensile specimen made of silicon nitrite (Si3N4) ceramic matrix composite was considered to evaluate structural durability and deformity. Ceramic matrix composites are being sought as candidate materials to replace nickel-base superalloys for turbine engine applications. They have the unique characteristics of being able to withstand higher operating temperatures and harsh combustion environments. In addition, their low densities relative to metals help reduce component mass (ref. 2). Detailed three-dimensional volume rendering of the tensile test specimen was successfully carried out with Velocity2 (ref. 1) using two-dimensional images that were generated via computed tomography. Subsequent, three-dimensional finite element analyses were performed, and the results obtained were compared with those predicted by NDE-based calculations and experimental tests. It was shown that Velocity2 software can be used to render a three-dimensional object from a series of CT scan images with a minimum level of complexity. The analytical results (ref. 3) show that the high-stress regions correlated well with the damage sites identified by the CT scans and the experimental data. Furthermore, modeling of the voids collected via NDE offered an analytical advantage that resulted in more accurate assessments of the material s structural strength. The top figure shows a CT scan image of the specimen test section illustrating various hidden structural entities in the material and an optical image of the test specimen considered in this study. The bottom figure represents the stress response predicted from the finite element analyses (ref .3 ) for a selected CT slice where it clearly illustrates the correspondence of the high stress risers due to voids in the material with those predicted by the NDE. This study is continuing, and efforts are concentrated on improving the modeling capabilities to imitate the structural anomalies as detected.

Abdul-Aziz, Ali↗
NASA NTRS1998

Mesogranulation as A Distinct Scale of Convection in the Sun

We present evidence for the existence of mesogranulation as a scale of convection distinct from granulation and supergranulation through analysis of full-disk Doppler velocity images of the Sun collected by the Michelson Doppler Imager (MDI) aboard the NASA/ESA Solar and Heliospheric Observatory (SOHO). Our analysis procedures isolate nearly steady flows in the solar photosphere and yield power spectra of convection for spherical harmonic degrees up to I = 1000. Each spectrum exhibits an obvious supergranulation peak at I approximately 130 and a broad secondary peak at I approximately 600 with a distinct break in the spectrum between these peaks at I approximately 300. We believe that this secondary peak is a signature of mesogranulation with typical cell diameters of about 7 Mm. Our standard analysis procedure is to first remove the p-mode oscillation signal by averaging individual Dopplergrams over 17-minute intervals. Next, by fitting to standard functional forms we remove Doppler signals due to the motion of the spacecraft, the convective blueshift, solar rotation including differential rotation, and the meridional circulation in order to produce Dopplergrams dominated by convective motions. By mapping these processed images onto heliographic coordinates and projecting onto spherical harmonics, we produce a power spectrum of solar convection for each 17-minute period. We construct synthetic images and pass them through the same analysis procedure in order to determine the actual solar convection spectrum that reproduces the analyzed results. We find that a small but increasing percentage of high-degree convective power is lost in the analysis as we approach the limit of resolution of the detector but'that the broad, mesogranulation peak at I approximately 600 must be included in the convection spectrum of the synthetic images.

Bachmann, Kurt T.↗
NASA NTRSDate not supplied

Pushing the Limits of Aquatic Remote Sensing: Synthetic Data and Deep Learning for Fast Inverse Emulation of A Coupled Ocean-Atmosphere Radiative Transfer Model

The inversion of electromagnetic information to physical and biological properties of the water column is a notoriously difficult problem, yet fundamental to our ability of understanding aquatic processes on large time and space scales. There is now a growing necessity to develop pragmatic approaches that allow timely and effective extrapolation of local processes, to spatially resolved global products, and to promote operational and sustainable resource policy management. This presentation will discuss research integrating advanced biological and radiative modeling, high-end computation, and machine learning to develop a portable global processor for simultaneous retrieval of atmosphere and water optics for diverse aquatic systems from the open and coastal ocean to optically extreme inland waters and harmful algal blooms. We will discuss some of the basic concepts behind the forward modeling approach including DEAP, the novel Distributed Equivalent Algal Populations model, for developing large spectral libraries of aquatic particle optics to aid in our ability to distinguish phytoplankton functional types (PFTs) and inorganic material, as well as other factors which enable comprehensive modeling from the benthos to top-of-atmosphere (TOA). This information is being used to understand how we can leverage next-generation deep learning methods for maximum information retrieval and rapid image processing, while also providing capabilities to identify minimum sensor spectral requirements necessary for certain aquatic applications. Further, I will touch on how we envision this research to enable the aquatic community for science discovery and how we are moving closer towards the capability for high-fidelity global analysis of aquatic ecosystems.

Jeremy Alan Kravitz↗
NASA NTRS1991

A cost-effective intelligent robotic system with dual-arm dexterous coordination and real-time vision

Dexterous coordination of manipulators based on the use of redundant degrees of freedom, multiple sensors, and built-in robot intelligence represents a critical breakthrough in development of advanced manufacturing technology. A cost-effective approach for achieving this new generation of robotics has been made possible by the unprecedented growth of the latest microcomputer and network systems. The resulting flexible automation offers the opportunity to improve the product quality, increase the reliability of the manufacturing process, and augment the production procedures for optimizing the utilization of the robotic system. Moreover, the Advanced Robotic System (ARS) is modular in design and can be upgraded by closely following technological advancements as they occur in various fields. This approach to manufacturing automation enhances the financial justification and ensures the long-term profitability and most efficient implementation of robotic technology. The new system also addresses a broad spectrum of manufacturing demand and has the potential to address both complex jobs as well as highly labor-intensive tasks. The ARS prototype employs the decomposed optimization technique in spatial planning. This technique is implemented to the framework of the sensor-actuator network to establish the general-purpose geometric reasoning system. The development computer system is a multiple microcomputer network system, which provides the architecture for executing the modular network computing algorithms. The knowledge-based approach used in both the robot vision subsystem and the manipulation control subsystems results in the real-time image processing vision-based capability. The vision-based task environment analysis capability and the responsive motion capability are under the command of the local intelligence centers. An array of ultrasonic, proximity, and optoelectronic sensors is used for path planning. The ARS currently has 18 degrees of freedom made up by two articulated arms, one movable robot head, and two charged coupled device (CCD) cameras for producing the stereoscopic views, and articulated cylindrical-type lower body, and an optional mobile base. A functional prototype is demonstrated.

Marzwell, Neville I.↗
NASA NTRS2002

Demonstrating Acquisition of Real-Time Thermal Data over Fires Utilizing UAVs

NASA-Ames Research Center, in collaboration with General Atomics Aeronautical Systems, Inc. has been developing real-time data acquisition and information delivery systems employing uninhabited aerial vehicle (UAV) technology for disaster mitigation and assessment demonstrations. Working in conjunction with the US Forest Service, a disaster community agency responsible for wildfire management and mitigation, we developed a large-scale wildfire demonstration called the First Response Experiment (FIRE). During that experiment in late summer 2001, the participants demonstrated the melding of innovative technologies such as UAV platforms, real-time data processing, and data telemetry for quick analysis of a disaster event. The General Atomics ALTUS UAV, the Airborne Infrared Disaster Assessment System (AIRDAS) and Over-The-Horizon (OTH) satellite data telemetry equipment were employed over a controlled burn to test the feasibility of a disaster monitoring and mitigation platform for hazardous duty. The ALTUS UAV was employed to demonstrate the long duration, altitude, and payload capability of unmanned platforms for acquiring disaster related data. The ALTUS has an operational altitude to 45,000 feet (13,700 in), with a flight duration of twenty-four hours and a payload capacity of over 300 lbs. (148.5 kg). This allows the platform to operate under the conditions that would be necessary for monitoring and mitigating disaster events throughout the Unites States. The four channel AIRDAS data (calibrated thermal infrared digital imagery of the fire event) was sent from the ALTUS UAV via a satellite communications system (NERA transponder and INMARSAT satellite) to a data archive server and an image processing work station at NASA-Ames Research Center, 400 miles away.

Wegener, Steven S.↗
NASA NTRS2004

Process for thermal imaging scanning of a swaged heater for an anode subassembly of a hollow cathode assembly

A process for thermal imaging scanning of a swaged heater of an anode subassembly of a hollow cathode assembly, comprising scanning a swaged heater with a thermal imaging radiometer to measure a temperature distribution of the heater; raising the current in a power supply to increase the temperature of the swaged heater; and measuring the swaged heater temperature using the radiometer, whereupon the temperature distribution along the length of the heater shall be less than plus or minus 5 degrees C.

Patterson, Michael J.↗
NASA NTRSDate not supplied

SWIPE: Spectral Water Inversion Processor and Emulator

Degradation of Earth’s inland water resources due to anthropogenic perturbations and climate anomalies at both local and global scales continues to place human health at substantial risk. There is now a growing necessity to develop pragmatic approaches that allow timely and effective extrapolation of local processes, to spatially resolved global products, and to promote operational and sustainable resource policy management. This presentation will be discussing the progress made developing SWIPE: Spectral Water Inversion Processor and Emulator. SWIPE is a platform for advanced modeling of coastal and inland aquatic habitats. The goal is create a comprehensive and cohesive system to leverage recent advancements in computation and machine learning to develop a synthetic training ground for sensitivity studies and algorithm development. The four principal facets of SWIPE include: 1. Advanced two-layer coated sphere bio-optical modeling and GPU radiative transfer modeling, 2. Big Data involving massive synthetic spectral libraries of optical properties of various global aquatic particles, surface reflectance, and top-of-atmosphere reflectance, all at hyperspectral resolution leveraging high-end computing systems at NASA Ames Research Center, 3. Deep Learning for algorithm development for water quality inversion of concentrations of common biogeophysical variables as well as optics, full uncertainty characterization by water type, and forward emulation, and lastly, 4. Image Processing for application of developed retrieval algorithms for both hyperspectral and multispectral sensors with experimental corrections for global adjacency, noise, sunglint, and benthic reflectance. This presentation will demonstrate the Equivalent Algal Populations (EAP) two-layer coated sphere scattering model which has been used develop spectral libraries of hyperspectral inherent optical properties of roughly 80 species of phytoplankton, covering 15 different classes and nine taxonomic functional types. The EAP model was also used to derive spectral properties of 10 different non-algal particle functional types. Examples of how the SMART-G (Speed-up Monte-carlo Advanced Radiative Transfer using GPU) radiative transfer code is used to model optically complex aquatic signals will be presented and discussed in the context of creating a massive synthetic database which can leverage the full power of next generation machine learning techniques and high end computing for water quality inversion. We will discuss our active investigation in things like appropriate model architectures, dimensionality reduction techniques such as PCA and autoencoders, uncertainty quantification and abstaining, and which variables actually benefit most from hyperspectral information versus multispectral resolution. We are also curious about questions relating to cost/benefit analysis in terms of computation resources, neural network complexity, and data volumes. Answers to these questions will hopefully elaborate on cost efficiency for potential future sensor design considerations.

SWIPE↗
NASA NTRS1988

Spatial vision processes: From the optical image to the symbolic structures of contour information

The significance of machine and natural vision is discussed together with the need for a general approach to image acquisition and processing aimed at recognition. An exploratory scheme is proposed which encompasses the definition of spatial primitives, intrinsic image properties and sampling, 2-D edge detection at the smallest scale, the construction of spatial primitives from edges, and the isolation of contour information from textural information. Concepts drawn from or suggested by natural vision at both perceptual and physiological levels are relied upon heavily to guide the development of the overall scheme. The scheme is intended to provide a larger context in which to place the emerging technology of detector array focal-plane processors. The approach differs from many recent efforts in edge detection and image coding by emphasizing smallest scale edge detection as a foundation for multi-scale symbolic processing while diminishing somewhat the importance of image convolutions with multi-scale edge operators. Cursory treatments of information theory illustrate that the direct application of this theory to structural information in images could not be realized.

Jobson, Daniel J.↗
NASA NTRS1992

Fortran for the nineties

Fortran has largely enjoyed prominence for the past few decades as the computer programming language of choice for numerically intensive scientific, engineering, and process control applications. Fortran's well understood static language syntax has allowed resulting parsers and compiler optimizing technologies to often generate among the most efficient and fastest run-time executables, particularly on high-end scalar and vector supercomputers. Computing architectures and paradigms have changed considerably since the last ANSI/ISO Fortran release in 1978, and while FORTRAN 77 has more than survived, it's aged features provide only partial functionality for today's demanding computing environments. The simple block procedural languages have been necessarily evolving, or giving way, to specialized supercomputing, network resource, and object-oriented paradigms. To address these new computing demands, ANSI has worked for the last 12-years with three international public reviews to deliver Fortran 90. Fortran 90 has superseded and replaced ISO FORTRAN 77 internationally as the sole Fortran standard; while in the US, Fortran 90 is expected to be adopted as the ANSI standard this summer, coexisting with ANSI FORTRAN 77 until at least 1996. The development path and current state of Fortran will be briefly described highlighting the many new Fortran 90 syntactic and semantic additions which support (among others): free form source; array syntax; new control structures; modules and interfaces; pointers; derived data types; dynamic memory; enhanced I/O; operator overloading; data abstraction; user optional arguments; new intrinsics for array, bit manipulation, and system inquiry; and enhanced portability through better generic control of underlying system arithmetic models. Examples from dynamical astronomy, signal and image processing will attempt to illustrate Fortran 90's applicability to today's general scalar, vector, and parallel scientific and engineering requirements and object oriented programming paradigms. Time permitting, current work proceeding on the future development of Fortran 2000 and collateral standards will be introduced.

Himer, J. T.↗
NASA NTRS2013

A System for Measuring the Sway of the Vehicle Assembly Building

A system was developed to measure the sway of the Vehicle Assembly Building (VAB) at Kennedy Space Center. This system was installed in the VAB and gathered more than one total year of data. The building movement was correlated with measurements provided by three wind towers in order to determine the maximum deflection of the building during high-wind events. The VAB owners were in the process of obtaining new platforms for use in assembling very tall rockets when analysis of the VAB showed that a high wind could move the building sufficiently that an upper platform might impact a rocket. The problem arises because safety requires a relatively small gap between the platform and the rocket, while a large enough gap is needed to ensure that stacking tolerances prevent contact between the rocket and the platform. This only leaves an inch or two (approximately 2 to 5 cm) of total clearance, so when the analysis showed that more than a couple of inches of motion could occur in a high wind, there was a potential for damaging the rocket. The KSC Applied Physics Laboratory was asked to install a system in the VAB that could measure the motion of the building in high winds to determine the actual building sway. The motion of the VAB roof under wind load was measured optically, and under analysis, it was determined that a relatively large-aperture optical system would be required to reduce diffraction effects to less than a small fraction of an inch (approximately mm) at a distance of 500 ft (≈150 m). A 10-in. (approximately 250 mm) telescope was placed on the floor of the building, looking at the ceiling. On the ceiling, a flat plate with three white LEDs was mounted in an "L" shape, such that the telescope was essentially looking at three stars. Software was written to track the motion of these three points using an image processing system. This provided a better than 1/10-in. (approximately 2.5-mm) 2D measurement faster than once a second. Data was downloaded once a month for comparison with the wind tower data. The system was fully operational and provided enough data to show that the VAB will only move 1 in. (approximately 2.5 cm) at the ceiling under 70-knot winds. Adjustable platforms are not required.

Youngquist, Robert↗
NASA NTRS2003

A Fast Implementation of the ISOCLUS Algorithm

Unsupervised clustering is a fundamental tool in numerous image processing and remote sensing applications. For example, unsupervised clustering is often used to obtain vegetation maps of an area of interest. This approach is useful when reliable training data are either scarce or expensive, and when relatively little a priori information about the data is available. Unsupervised clustering methods play a significant role in the pursuit of unsupervised classification. One of the most popular and widely used clustering schemes for remote sensing applications is the ISOCLUS algorithm, which is based on the ISODATA method. The algorithm is given a set of n data points (or samples) in d-dimensional space, an integer k indicating the initial number of clusters, and a number of additional parameters. The general goal is to compute a set of cluster centers in d-space. Although there is no specific optimization criterion, the algorithm is similar in spirit to the well known k-means clustering method in which the objective is to minimize the average squared distance of each point to its nearest center, called the average distortion. One significant feature of ISOCLUS over k-means is that clusters may be merged or split, and so the final number of clusters may be different from the number k supplied as part of the input. This algorithm will be described in later in this paper. The ISOCLUS algorithm can run very slowly, particularly on large data sets. Given its wide use in remote sensing, its efficient computation is an important goal. We have developed a fast implementation of the ISOCLUS algorithm. Our improvement is based on a recent acceleration to the k-means algorithm, the filtering algorithm, by Kanungo et al.. They showed that, by storing the data in a kd-tree, it was possible to significantly reduce the running time of k-means. We have adapted this method for the ISOCLUS algorithm. For technical reasons, which are explained later, it is necessary to make a minor modification to the ISOCLUS specification. We provide empirical evidence, on both synthetic and Landsat image data sets, that our algorithm's performance is essentially the same as that of ISOCLUS, but with significantly lower running times. We show that our algorithm runs from 3 to 30 times faster than a straightforward implementation of ISOCLUS. Our adaptation of the filtering algorithm involves the efficient computation of a number of cluster statistics that are needed for ISOCLUS, but not for k-means.

Memarsadeghi, Nargess↗
NASA NTRS2013

A Combined Very Large Telescope and Gemini Study of the Atmosphere of the Directly Imaged Planet, Beta Pictoris b

We analyze new/archival VLT/NaCo and Gemini/NICI high-contrast imaging of the young, self-luminous planet Beta Pictoris b in seven near-to-mid IR photometric filters, using advanced image processing methods to achieve high signal-to-noise, high precision measurements. While Beta Pic b's near-IR colors mimic those of a standard, cloudy early-to-mid L dwarf, it is overluminous in the mid-infrared compared to the field L/T dwarf sequence. Few substellar/planet-mass objects-i.e., κ And b and 1RXJ 1609B-match Beta Pic b's JHKsL photometry and its 3.1 micron and 5 micron photometry are particularly difficult to reproduce. Atmosphere models adopting cloud prescriptions and large (approx. 60 micron)dust grains fail to reproduce the Beta Pic b spectrum. However, models incorporating thick clouds similar to those found forHR8799 bcde, but also with small (a fewmicrons) modal particle sizes, yield fits consistent with the data within the uncertainties. Assuming solar abundance models, thick clouds, and small dust particles (a = 4 micron), we derive atmosphere parameters of log(g) = 3.8 +/- 0.2 and Teff = 1575-1650 K, an inferred mass of 7+4 −3 MJ, and a luminosity of log(L/L) approx. −3.80 +/- 0.02. The best-estimated planet radius, is approx. equal to 1.65 +/- 0.06 RJ, is near the upper end of allowable planet radii for hot-start models given the host star's age and likely reflects challenges constructing accurate atmospheric models. Alternatively, these radii are comfortably consistent with hot-start model predictions if Beta Pic b is younger than is approx. equal to 7 Myr, consistent with a late formation well after its host star's birth approx. 12+8 −4 Myr ago.

overluminous↗