Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time-series analysis”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Tropical Tropospheric Ozone (TTO) Maps from Nimbus 7 and Earth-Probe TOMS by the Modified-Residual Method: Validation, Evaluation and Trends based on Atlantic Regional Time Series - 1

The well-known wave-one pattern seen in tropical total ozone [Shiotani, 1992; Ziemke et al., 1996, 1998] has been used to develop a modified-residual (MR) method for retrieving time-averaged stratospheric ozone and tropospheric ozone column amount from TOMS (Total Ozone Mapping Spectrometer) over the 14 complete calendar years of Nimbus 7 observations (1979-1992) and from TOMS on the Earth-Probe (1996-present) and ADEOS platforms (1996- 1997). Nine- to sixteen-day averaged tropical tropospheric ozone (TTO) maps, validated with ozonesondes, show a seasonality expected from dynamical and chemical influences. The maps may be viewed on a homepage: http://metosrv2.umd.edu/~tropo. Stratospheric column ozone, which is also derived by the modified-residual method, compares well with sondes (to within 6-7 DU) and with stratospheric ozone column derived from other satellites (within 8-10 DU). Validation of the TTO time-series is presently limited to ozonesonde comparisons with Atlantic stations and sites on the adjacent continents (Ascension Island, Natal, Brazil; Brazzaville); for the sounding periods, TTO at all locations agrees with the sonde record to +/-7 DU. TTO time-series and the magnitude of the wave-one pattern show ENSO signals in the strongest El Nifio periods from 1979-1998. From 12degN and 12degS, zonally averaged tropospheric ozone shows no significant trend from 1980-1990. Trends are also not significant during this period in localized regions, e.g. from just west of South America across to southern Africa. This is consistent with the ozonesonde record at Natal, Brazil (the only tropical ozone data publicly available for the 1980's), which shows a not statistically significant increase. The lack of trend in tropospheric ozone agrees with a statistical analysis based on another method for deriving TTO from TOMS, the so-called Convective-Cloud-Differential approach of Ziemke et al. [1998].

Thompson, Anne M.↗
NASA NTRS2012

Demystifying Kepler Data: A Primer for Systematic Artifact Mitigation

The Kepler spacecraft has collected data of high photometric precision and cadence almost continuously since operations began on 2009 May 2. Primarily designed to detect planetary transits and asteroseismological signals from solar-like stars, Kepler has provided high quality data for many areas of investigation. Unconditioned simple aperture time-series photometry are however affected by systematic structure. Examples of these systematics are differential velocity aberration, thermal gradients across the spacecraft, and pointing variations. While exhibiting some impact on Kepler's primary science, these systematics can critically handicap potentially ground-breaking scientific gains in other astrophysical areas, especially over long timescales greater than 10 days. As the data archive grows to provide light curves for 10(exp 5) stars of many years in length, Kepler will only fulfill its broad potential for stellar astrophysics if these systematics are understood and mitigated. Post-launch developments in the Kepler archive, data reduction pipeline and open source data analysis software have occurred to remove or reduce systematic artifacts. This paper provides a conceptual primer for users of the Kepler data archive to understand and recognize systematic artifacts within light curves and some methods for their removal. Specific examples of artifact mitigation are provided using data available within the archive. Through the methods defined here, the Kepler community will find a road map to maximizing the quality and employment of the Kepler legacy archive.

Kinemuchi, K.↗
NASA NTRS1999

Psychological changes in hundred-day remote Antarctic field groups

Psychological adaptation to extreme environments has been examined from several perspectives. In this study, two Australian teams, each consisting of six male crew members, completed computer-administered questionnaires twice weekly during 100-day traverses around the Lambert Glacier Basin, Antarctica. Only small trends were noted when data were aggregated at the group level, which is consistent with the findings of others. Data were then analyzed using pooled time-series regression. These analyses incorporated personality characteristics, environmental factors, and interpersonal factors as predictors of Group Tensions, Personal Morale, Emotional State, Cognitive Readiness, and the Team's Work Life. Most of the psychological discomfort and problems that occurred appeared to be within the individual or between individuals. They did not affect all members of the group equally.

NASA Discipline Space Human Factors↗
NASA NTRS2001

Dynamical Interpolation of Mesoscale Flows in the TOPEX/Poseidon Diamond Surrounding the U.S. Joint Global Ocean Flux Study Bermuda Atlantic Time-Series Study Site

An open boundary ocean model is configured in a domain bounded by the four TOPEX/Poseidon (T/P) ground tracks surrounding the US Joint Global Ocean Flux Study Bermuda Atlantic Time-Series Study (BATS) site. This implementation facilitates prescription of model boundary conditions directly from altimetric measurements (both TIP and ERS-2). The expected error characteristics for a domain of this size with periodically updated boundary conditions are established with idealized numerical experiments using simulated data. A hindcast simulation is then constructed using actual altimetric observations during the period October 1992 through September 1998. Quantitative evaluation of the simulation suggests significant skill. The correlation coefficient between predicted sea level anomaly and ERS observations in the model interior is 0.89; that for predicted versus observed dynamic height anomaly based on hydrography at the BATS site is 0.73. Comparison with the idealized experiments suggests that the main source of error in the hindcast is temporal undersampling of the boundary conditions. The hindcast simulation described herein provides a basis for retrospective analysis of BATS observations in the context of the mesoscale eddy field.

McGillicuddy, Dennis J., Jr.↗
NASA NTRS2001

Dynamical Interpolation of Mesoscale Flows in the TOPEX/ Poseidon Diamond Surrounding the U.S. Joint Global Ocean Flux Study Bermuda Atlantic Time-Series Study Site

An open boundary ocean model is configured in a domain bounded by the four TOPEX/Poseidon (TIP) ground tracks surrounding the U.S. Joint Global Ocean Flux Study Bermuda Atlantic Time-series Study (BATS) site. This implementation facilitates prescription of model boundary conditions directly from altimetric measurements (both TIP and ERS-2). The expected error characteristics for a domain of this size with periodically updated boundary conditions are established with idealized numerical experiments using simulated data. A hindcast simulation is then constructed using actual altimetric observations during the period October 1992 through September 1998. Quantitative evaluation of the simulation suggests significant skill. The correlation coefficient between predicted sea level anomaly and ERS observations in the model interior is 0.89; that for predicted versus observed dynamic height anomaly based on hydrography at the BATS site is 0.73. Comparison with the idealized experiments suggests that the main source of error in the hindcast is temporal undersampling of the boundary conditions. The hindcast simulation described herein provides a basis for retrospective analysis of BATS observations in the context of the mesoscale eddy field.

McGillicuddy, D. J.↗
NASA NTRS2011

Modeling the Height of Young Forests Regenerating from Recent Disturbances in Mississippi using Landsat and ICESat data

Many forestry and earth science applications require spatially detailed forest height data sets. Among the various remote sensing technologies, lidar offers the most potential for obtaining reliable height measurement. However, existing and planned spaceborne lidar systems do not have the capability to produce spatially contiguous, fine resolution forest height maps over large areas. This paper describes a Landsat-lidar fusion approach for modeling the height of young forests by integrating historical Landsat observations with lidar data acquired by the Geoscience Laser Altimeter System (GLAS) instrument onboard the Ice, Cloud, and land Elevation (ICESat) satellite. In this approach, "young" forests refer to forests reestablished following recent disturbances mapped using Landsat time-series stacks (LTSS) and a vegetation change tracker (VCT) algorithm. The GLAS lidar data is used to retrieve forest height at sample locations represented by the footprints of the lidar data. These samples are used to establish relationships between lidar-based forest height measurements and LTSS-VCT disturbance products. The height of "young" forest is then mapped based on the derived relationships and the LTSS-VCT disturbance products. This approach was developed and tested over the state of Mississippi. Of the various models evaluated, a regression tree model predicting forest height from age since disturbance and three cumulative indices produced by the LTSS-VCT method yielded the lowest cross validation error. The R(exp 2) and root mean square difference (RMSD) between predicted and GLAS-based height measurements were 0.91 and 1.97 m, respectively. Predictions of this model had much higher errors than indicated by cross validation analysis when evaluated using field plot data collected through the Forest Inventory and Analysis Program of USDA Forest Service. Much of these errors were due to a lack of separation between stand clearing and non-stand clearing disturbances in current LTSS-VCT products and difficulty in deriving reliable forest height measurements using GLAS samples when terrain relief was present within their footprints. In addition, a systematic underestimation of about 5 m by the developed model was also observed, half of which could be explained by forest growth that occurred between field measurement year and model target year. The remaining difference suggests that tree height measurements derived using waveform lidar data could be significantly underestimated, especially for young pine forests. Options for improving the height modeling approach developed in this study were discussed.

Li, Ainong↗
NASA NTRS1998

Ice Surface Elevation Changes in East Antarctica from Satellite Altimetry

Estimates of the overall mass balance and seasonal and inter-annual variations in the surface mass balance are obtainable from time-series of ice surface elevations measured by satellite altimetry. Beginning in 2001, NASA's ICESat laser altimeter and lidar mission will significantly improve the range accuracy, the orbit accuracy, and the spatial coverage for measurement of ice sheet elevations (to 86 S) , as compared to previous radar altimeters designed for ocean measurements The radar altimeters on Seasat and Geosat provided ice sheet measurements to 72 S, and on ERS-1 and ERS-2 to 81 S. Although radar altimetry has significant limitations in coverage (due to loss of tracking) and accuracy over sloping surfaces, information on ice-sheet surface-elevation changes has been derived for parts of Antarctica. Recently, the accuracy of the ice measurements by Seasat (3 months of 1978) and Geosat (1985 to 1989) have been improved by new calculations of the satellite orbit heights and other altimeter corrections. Residual orbit errors and inter-satellite biases are evaluated by crossover analysis and by global adjustments to an ocean surface derived from altimeter data. The standard deviation of the orbit error is less than 9 cm, and the long-term trend in the error appears to be less than 1 cm/yr. Orbit errors can be further reduced by adjustment to the ocean surface, but false signals of several cm/yr may be also introduced by the adjustments. These false signals are caused mainly by residual errors in the altimeter corrections over the ocean, and secondary by real changes in the ocean surface elevation. Maps of ice sheet elevation changes north of 72 S are derived from Seasat-Geosat crossovers and from 4.5 years of Geosat crossovers. A notable ice thinning rate of about 50 cm/yr is found at elevations below 2200 meters between 70 and 72 S to the East of the Amery ice shelf, in both the Seasat-Geosat and Geosat-Geosat time intervals Above 2200 meters, to the ridge around 95 E, the thinning is about 10 to 20 cm/yr. For the Seasat/Geosat interval, a thickening is noted on the area East of the ridge between about 95 E and 115 E. Preliminary results from analysis of ERS-1 and ERS-2 radar altimeter data are also described.

Zwally, H. Jay↗
NASA NTRS2012

Remote Sensing of Tropical Ecosystems: Atmospheric Correction and Cloud Masking Matter

Tropical rainforests are significant contributors to the global cycles of energy, water and carbon. As a result, monitoring of the vegetation status over regions such as Amazonia has been a long standing interest of Earth scientists trying to determine the effect of climate change and anthropogenic disturbance on the tropical ecosystems and its feedback on the Earth's climate. Satellite-based remote sensing is the only practical approach for observing the vegetation dynamics of regions like the Amazon over useful spatial and temporal scales, but recent years have seen much controversy over satellite-derived vegetation states in Amazônia, with studies predicting opposite feedbacks depending on data processing technique and interpretation. Recent results suggest that some of this uncertainty could stem from a lack of quality in atmospheric correction and cloud screening. In this paper, we assess these uncertainties by comparing the current standard surface reflectance products (MYD09, MYD09GA) and derived composites (MYD09A1, MCD43A4 and MYD13A2 - Vegetation Index) from the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard the Aqua satellite to results obtained from the Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm. MAIAC uses a new cloud screening technique, and novel aerosol retrieval and atmospheric correction procedures which are based on time-series and spatial analyses. Our results show considerable improvements of MAIAC processed surface reflectance compared to MYD09/MYD13 with noise levels reduced by a factor of up to 10. Uncertainties in the current MODIS surface reflectance product were mainly due to residual cloud and aerosol contamination which affected the Normalized Difference Vegetation Index (NDVI): During the wet season, with cloud cover ranging between 90 percent and 99 percent, conventionally processed NDVI was significantly depressed due to undetected clouds. A smaller reduction in NDVI due to increased aerosol levels was observed during the dry season, with an inverse dependence of NDVI on aerosol optical thickness (AOT). NDVI observations processed with MAIAC showed highly reproducible and stable inter-annual patterns with little or no dependence on cloud cover, and no significant dependence on AOT (p less than 0.05). In addition to a better detection of cloudy pixels, MAIAC obtained about 20-80 percent more cloud free pixels, depending on season, a considerable amount for land analysis given the very high cloud cover (75-99 percent) observed at any given time in the area. We conclude that a new generation of atmospheric correction algorithms, such as MAIAC, can help to dramatically improve vegetation estimates over tropical rain forest, ultimately leading to reduced uncertainties in satellite-derived vegetation products globally.

tropical ecosystems↗
NASA NTRS2003

NASA Tech Briefs, November 2003

Topics covered include: Computer Program Recognizes Patterns in Time-Series Data; Program for User-Friendly Management of Input and Output Data Sets; Noncoherent Tracking of a Source of a Data-Modulated Signal; Software for Acquiring Image Data for PIV; Detecting Edges in Images by Use of Fuzzy Reasoning; A Timer for Synchronous Digital Systems; Prototype Parts of a Digital Beam-Forming Wide-Band Receiver; High-Voltage Droplet Dispenser; Network Extender for MIL-STD-1553 Bus; MMIC HEMT Power Amplifier for 140 to 170 GHz; Piezoelectric Diffraction-Based Optical Switches; Numerical Modeling of Nanoelectronic Devices; Organizing Diverse, Distributed Project Information; Eigensolver for a Sparse, Large Hermitian Matrix; Modified Polar-Format Software for Processing SAR Data; e-Stars Template Builder; Software for Acoustic Rendering; Functionally Graded Nanophase Beryllium/Carbon Composites; Thin Thermal-Insulation Blankets for Very High Temperatures; Prolonging Microgravity on Parabolic Airplane Flights; Device for Locking a Control Knob; Cable-Dispensing Cart; Foam Sensor Structures Would be Self-Deployable and Survive Hard Landings; Real-Gas Effects on Binary Mixing Layers; Earth-Space Link Attenuation Estimation via Ground Radar Kdp; Wedge Heat-Flux Indicators for Flash Thermography; Measuring Diffusion of Liquids by Common-Path Interferometry; Zero-Shear, Low-Disturbance Optical Delay Line; Whispering-Gallery Mode-Locked Lasers; Spatial Light Modulators as Optical Crossbar Switches; Update on EMD and Hilbert-Spectra Analysis of Time Series; Quad-Tree Visual-Calculus Analysis of Satellite Coverage; Dyakonov-Perel Effect on Spin Dephasing in n-Type GaAs; Update on Area Production in Mixing of Supercritical Fluids; and Quasi-Sun-Pointing of Spacecraft Using Radiation Pressure.

Source record↗
NASA NTRS2004

Advances in Measuring Antarctic Sea-Ice Thickness and Ice-Sheet Elevations with ICESat Laser Altimetry

NASA's Ice, Cloud and Land Elevation Satellite (ICESat) has been measuring elevations of the Antarctic ice sheet and sea-ice freeboard elevations with unprecedented accuracy. Since February 20,2003, data has been acquired during three periods of laser operation varying from 36 to 54 days, which is less than the continuous operation of 3 to 5 years planned for the mission. The primary purpose of ICESat is to measure time-series of ice-sheet elevation changes for determination of the present-day mass balance of the ice sheets, study of associations between observed ice changes and polar climate, and estimation of the present and future contributions of the ice sheets to global sea level rise. ICESat data will continue to be acquired for approximately 33 days periods at 3 to 6 month intervals with the second of ICESat's three lasers, and eventually with the third laser. The laser footprints are about 70 m on the surface and are spaced at 172 m along-track. The on-board GPS receiver enables radial orbit determinations to an accuracy better than 5 cm. The orbital altitude is around 600 km at an inclination of 94 degrees with a 8-day repeat pattern for the calibration and validation period, followed by a 91 -day repeat period for the rest of the mission. The expected range precision of single footprint measurements was 10 cm, but the actual range precision of the data has been shown to be much better at 2 to 3 cm. The star-tracking attitude-determination system should enable footprints to be located to 6 m horizontally when attitude calibrations are completed. With the present attitude calibration, the elevation accuracy over the ice sheets ranges from about 30 cm over the low-slope areas to about 80 cm over areas with slopes of 1 to 2 degrees, which is much better than radar altimetry. After the first period of data collection, the spacecraft attitude was controlled to point the laser beam to within 50 m of reference surface tracks over the ice sheets. Detection of ice elevation changes over select areas of the ice sheet is demonstrated with using both crossover analysis and precise-repeat track analysis. Sea ice freeboard-height distributions over the Antarctic sea pack are derived over distances of 50 km and converted into maps of average freeboard thickness and sea-ice thickness.

Zwally, H. Jay↗
NASA NTRS1971

A real-time statistical time-series analyzer

Device extracts average frequency of human speech and produces second, third, and fourth moments of instantaneous frequency about this average. It operates on electrical time representation of input signal, performs statistical analysis on zero-crossing of almost any signal, and does not require specialized personnel to operate it.

Stewart, C. H.↗
NASA NTRS1989

An introduction to chaotic and random time series analysis

The origin of chaotic behavior and the relation of chaos to randomness are explained. Two mathematical results are described: (1) a representation theorem guarantees the existence of a specific time-domain model for chaos and addresses the relation between chaotic, random, and strictly deterministic processes; (2) a theorem assures that information on the behavior of a physical system in its complete state space can be extracted from time-series data on a single observable. Focus is placed on an important connection between the dynamical state space and an observable time series. These two results lead to a practical deconvolution technique combining standard random process modeling methods with new embedded techniques.

Scargle, Jeffrey D.↗
NASA NTRS1989

Wavelet analysis introduction and application to radar scattering from water waves

Double-parameter expansions, leading to an instantaneous multiscale analysis, are used to obtain time-frequency graphs from time-series data for interpretation of scatterometer response to sea surface geometry. This method is used to define, as well as possible, one-dimensional signal input (elevation and slope). However, some important fundamental parameters, such as the asymmetric development of short scales in crests and troughs, the long-wave amplitude modulation, the different zero-crossings, and so forth, are easily extracted from the wavelet images.

Chapron, B.↗
NASA NTRS1994

Time-series photometric spot modeling. 2: Fifteen years of photometry of the bright RS CVn binary HR 7275

We present a time-dependent spot modeling analysis of 15 consecutive years of V-band photometry of the long-period (P(sub orb) = 28.6 days) RS CVn binary HR 7275. This baseline in time is one of the longest, uninterrupted intervals a spotted star has been observed. The spot modeling analysis yields a total of 20 different spots throughout the time span of our observations. The distribution of the observed spot migration rates is consistent with solar-type differential rotation and suggests a lower limit of the differential-rotation coefficient of 0.022 +/-0.004. The observed, maximum lifetime of a single spot (or spot group) is 4.5 years, the minimum lifetime is approximately one year, but an average spot lives for 2.2 years. If we assume that the mechanical shear by differential rotation sets the upper limit to the spot lifetime, the observed maximum lifetime in turn sets an upper limit to the differential-rotation coefficient, namely 0.04 +/- 0.01. This would be differential rotation just 5 to 8 times less than the solar value and one of the strongest among active binaries. We found no conclusive evidence for the existence of a periodic phenomenon that could be attributed to a stellar magnetic cycle.

Strassmeier, K. G.↗
NASA NTRS1994

(abstract) Mount Rainier: New Remote Sensing Observations of a Decade Volcano

Mount Rainier was selected as a Decade Volcano by the International Association of Volcanology and Chemistry of the Earth's Interior. The purpose of this selection is to focus scientific and public attention on Mount Rainier during the current decade, the United Nations-designated International Decade of Natural Hazard Reduction. The Mount Rainier science plan calls for remote sensing surveys to monitor the volcano. To date, we have conducted airborne surveys with visible and near-infrared, thermal infrared, and interferometric radar instruments. Our preliminary analysis of some night-time time-series thermal infrared survey data sets of the summit suggests that, aside from seasonal variations in snow cover, there have been no qualitative changes in the size or pattern of the summit hot spots. Day-time airborne surveys were done to record the current surface appearance of the volcano and map hydrothermal alteration in the summit region. An interferometric radar survey yielded a high-resolution digital elevation model (DEM) which serves as a base for the registration of the other remote sensing data sets. More importantly, the DEM documents the current topography of glaciers and valleys. Planned biannual radar survey of mount rainier will produce a data set from which seasonal changes in glacier and valley topography can be characterized. Such characterization is essential if we are to recognize geothermally induced changes in snow and ice cover.

infrared radar interferometry glaciers seasonal va↗
NASA NTRS2007

Development of a Rayleigh Scattering Diagnostic for Time-Resolved Gas Flow Velocity, Temperature, and Density Measurements in Aerodynamic Test Facilities

A molecular Rayleigh scattering technique is developed to measure time-resolved gas velocity, temperature, and density in unseeded turbulent flows at sampling rates up to 32 kHz. A high power continuous-wave laser beam is focused at a point in an air flow field and Rayleigh scattered light is collected and fiber-optically transmitted to the spectral analysis and detection equipment. The spectrum of the light, which contains information about the temperature and velocity of the flow, is analyzed using a Fabry-Perot interferometer. Photomultiplier tubes operated in the photon counting mode allow high frequency sampling of the circular interference pattern to provide time-resolved flow property measurements. An acoustically driven nozzle flow is studied to validate velocity fluctuation measurements, and an asymmetric oscillating counterflow with unequal enthalpies is studied to validate the measurement of temperature fluctuations. Velocity fluctuations are compared with constant temperature anemometry measurements and temperature fluctuations are compared with constant current anemometry measurements at the same locations. Time-series and power spectra of the temperature and velocity measurements are presented. A numerical simulation of the light scattering and detection process was developed and compared with experimental data for future use as an experiment design tool.

Mielke, Amy F.↗
NASA NTRS2008

Remote-Sensing Time Series Analysis, a Vegetation Monitoring Tool

The Time Series Product Tool (TSPT) is software, developed in MATLAB , which creates and displays high signal-to- noise Vegetation Indices imagery and other higher-level products derived from remotely sensed data. This tool enables automated, rapid, large-scale regional surveillance of crops, forests, and other vegetation. TSPT temporally processes high-revisit-rate satellite imagery produced by the Moderate Resolution Imaging Spectroradiometer (MODIS) and by other remote-sensing systems. Although MODIS imagery is acquired daily, cloudiness and other sources of noise can greatly reduce the effective temporal resolution. To improve cloud statistics, the TSPT combines MODIS data from multiple satellites (Aqua and Terra). The TSPT produces MODIS products as single time-frame and multitemporal change images, as time-series plots at a selected location, or as temporally processed image videos. Using the TSPT program, MODIS metadata is used to remove and/or correct bad and suspect data. Bad pixel removal, multiple satellite data fusion, and temporal processing techniques create high-quality plots and animated image video sequences that depict changes in vegetation greenness. This tool provides several temporal processing options not found in other comparable imaging software tools. Because the framework to generate and use other algorithms is established, small modifications to this tool will enable the use of a large range of remotely sensed data types. An effective remote-sensing crop monitoring system must be able to detect subtle changes in plant health in the earliest stages, before the effects of a disease outbreak or other adverse environmental conditions can become widespread and devastating. The integration of the time series analysis tool with ground-based information, soil types, crop types, meteorological data, and crop growth models in a Geographic Information System, could provide the foundation for a large-area crop-surveillance system that could identify a variety of plant phenomena and improve monitoring capabilities.

McKellip, Rodney↗
NASA NTRS2003

Update on EMD and Hilbert-Spectra Analysis of Time Series

This method is especially well suited for analyzing time-series data that represent nonstationary and nonlinear physical phenomena. The method is based principally on the concept of empirical mode decomposition (EMD), according to which any complicated signal (as represented by digital samples) can be decomposed into a finite number of functions, called "intrinsic mode functions" (IMFs), that admit well-behaved Hilbert transforms. The local energies and the instantaneous frequencies derived from the IMFs through Hilbert transforms can be used to construct an energy-frequency-time distribution, denoted a Hilbert spectrum.

Huang, Norden E.↗