Search NASA⌕ Search

SEARCH · Search NASA

Results for “Geoscience”

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 289 records · Page 16

Recent Results in ICESat Geolocation Parameter Calibration from Range Residual Analysis

The Ice, Cloud and land Elevation Satellite (ICESat), launched on January 12, 2003, carries the Geoscience Laser Altimeter System (GLAS), which has a primary measurement of short-pulse laser-ranging to the Earths surface at 1064 nm wavelength at a rate of 40 pulses per second. The instrument has collected (at the time of this writing) over 540 million observations of the Earth's surface including precise elevation measurements of the ice sheets, sea ice roughness and thickness, ocean and land surface elevations and surface reflectivity. The accurate geolocation of GLASs surface returns, the spots from which the laser energy reflects on the Earths surface, is a critical issue in the scientific application of these data. Pointing, ranging, timing and orbit errors must be compensated to accurately geolocate the laser altimeter surface returns. Towards this end, the laser range observations can be fully exploited in an integrated residual analysis to accurately calibrate these geolocation/instrument parameters. The ICESat laser altimeter data have been simultaneously processed as direct altimetry from ocean sweeps and around-the-world scans along with dynamic crossovers in order to calibrate pointing ranging and timing. The analysis has resolved complex time varying signal resulting from thermal variations of instrument pointing and boresite shadowing. The calibration methodology and current calibration results are discussed along with future efforts and current geolocation accuracies.

Luthcke, Scott B.↗

Aerosol and Cloud Observations and Data Products by the GLAS Polar Orbiting Lidar Instrument

The Geoscience Laser Altimeter System (GLAS) launched in 2003 is the first polar orbiting satellite lidar. The instrument was designed for high performance observations of the distribution and optical scattering cross sections of clouds and aerosol. The backscatter lidar operates at two wavelengths, 532 and 1064 nm. Both receiver channels meet and exceed their design goals, and beginning with a two month period through October and November 2003, an excellent global lidar data set now exists. The data products for atmospheric observations include the calibrated, attenuated backscatter cross section for cloud and aerosol; height detection for multiple cloud layers; planetary boundary layer height; cirrus and aerosol optical depth and the height distribution of aerosol and cloud scattering cross section profiles. The data sets are now in open release through the NASA data distribution system. The initial results on global statistics for cloud and aerosol distribution has been produced and in some cases compared to other satellite observations. The sensitivity of the cloud measurements is such that the 70% global cloud coverage result should be the most accurate to date. Results on the global distribution of aerosol are the first that produce the true height distribution for model inter-comparison.

Spinhirne, J. D.↗

Antarctica Cloud Cover for October 2003 from GLAS Satellite Lidar Profiling

Seeing clouds in polar regions has been a problem for the imagers used on satellites. Both clouds and snow and ice are white, which makes clouds over snow hard to see. And for thermal infrared imaging both the surface and the clouds cold. The Geoscience Laser Altimeter System (GLAS) launched in 2003 gives an entirely new way to see clouds from space. Pulses of laser light scatter from clouds giving a signal that is separated in time from the signal from the surface. The scattering from clouds is thus a sensitive and direct measure of the presence and height of clouds. The GLAS instrument orbits over Antarctica 16 times a day. All of the cloud observations for October 2003 were summarized and compared to the results from the MODIS imager for the same month. There are two basic cloud types that are observed, low stratus with tops below 3 km and high cirrus form clouds with cloud top altitude and thickness tending at 12 km and 1.3 km respectively. The average cloud cover varies from over 93 % for ocean and coastal regions to an average of 40% over the East Antarctic plateau and 60-90% over West Antarctica. When the GLAS monthly average cloud fractions are compared to the MODIS cloud fraction data product, differences in the amount of cloud cover are as much as 40% over the continent. The results will be used to improve the way clouds are detected from the imager observations. These measurements give a much improved understanding of distribution of clouds over Antarctica and may show how they are changing as a result of global warming.

Spinhirne, J. D.↗

ICESat Observations of Topographic Change in the Northern Segment of the 2004 Sumatra-Andaman Islands Earthquake Rupture Zone

The Andaman Islands are located 120 km east of the Sunda trench in the northern quarter of the 1300 km long rupture zone of the 2004 Sumatra-Andaman Islands earthquake inferred from the distribution of aftershocks. Initial field reports indicate that several meters of uplift and up to a meter of submergence occurred on the western and eastern shorelines of the Andaman Islands, respectively, associated with the earthquake (Bilham, 2005). Satellite images also document uplift of western shoreline coral reef platforms above sea level. Body-wave (Ji, 2005; Yamamaka, 2005) and tide-gauge (Ortiz, 2005) slip inversions only resolve coseismic slip in the southern one-third to one-half of the rupture zone. The amount of coseismic slip in the Andaman Islands region is poorly constrained by these inversions. The Ice, Cloud, and land Elevation Satellite (ICESat), a part of the NASA Earth Observing System, is being used to document the spatial pattern of Andaman Islands vertical displacements in order to constrain models of slip distribution in the northern part of the rupture zone. ICESat carries the Geoscience Laser Altimeter System (GLAS) that obtains elevation measurements from 80 m diameter footprints spaced 175 m apart along profiles. For surfaces of low slope, single-footprint absolute elevation and horizontal accuracies of 10 cm and 6 m (1 sigma), respectively, referenced to the ITRF 2002 TOPEX/Poseidon ellipsoid are being obtained. Laser pulse backscatter waveforms enable separation of ground topography and overlying vegetation cover. During each 33-day observing period ICESat acquires three profiles crossing the Andaman Islands. A NNE-SSW oriented track consists of 1600 laser footprints along the western side of North, Middle, and South Andaman Islands and 240 laser footprints across the center of Great Andaman Island. Two NNW-SSE tracks consist of 440 footprints across Middle Andaman Island and 25 footprints across the west side of Sentinel Island. Cloud-free profiles were acquired in the fall of 2003 and 2004. During February-March, 2005 ICESat's precise pointing capability will be used to exactly repeat these three profiles, with a cross-track accuracy of better than 100 m, providing trench- parallel and -perpendicular observations of topographic change of the Andaman Islands that will compliment geodetic field surveys. The observed elevation changes will be compared to models of coseismic deformation associated with the mainshock and large aftershocks in the Andaman Islands region.

Harding, David↗

Height Distribution Between Cloud and Aerosol Layers from the GLAS Spaceborne Lidar in the Indian Ocean Region

The Geoscience Laser Altimeter System (GLAS), a nadir pointing lidar on the Ice Cloud and land Elevation Satellite (ICESat) launched in 2003, now provides important new global measurements of the relationship between the height distribution of cloud and aerosol layers. GLAS data have the capability to detect, locate, and distinguish between cloud and aerosol layers in the atmosphere up to 40 km altitude. The data product algorithm tests the product of the maximum attenuated backscatter coefficient b'(r) and the vertical gradient of b'(r) within a layer against a predetermined threshold. An initial case result for the critical Indian Ocean region is presented. From the results the relative height distribution between collocated aerosol and cloud shows extensive regions where cloud formation is well within dense aerosol scattering layers at the surface. Citation: Hart, W. D., J. D. Spinhime, S. P. Palm, and D. L. Hlavka (2005), Height distribution between cloud and aerosol layers from the GLAS spaceborne lidar in the Indian Ocean region,

Hart, William D.↗

The Beauty and Complexity of the Brunt Ice Shelf from MOA and ICESat

Beginning in February 2003, NASA's Ice, Cloud, and land Elevation Satellite (ICESat) has determined surface elevations from approx. 86degN to 86degS latitude. To date, altimetry data have been acquired in a series of observation periods in repeated track patterns using all three Geoscience Laser Altimeter System (GLAS) lasers. This paper will focus on ice shelf elevation data that were obtained in 2003 across the Brunt Ice Shelf and the Stancomb-Wills Ice Tongue. Integrating the altimetry with the recently available MODIS Mosaic of Antarctica (MOA), quantifies the relative accuracy and precision of the resulting ice shelf elevations. Furthermore, the elevation data was processed onto an elevation grid, by regional interpolation across the area s complex glacial features only. Ice thickness estimation from the altimetry of the floating ice is discussed. ICESat operates at 40Hz and its elevation data is obtained every 172m along track. These elevations have a relative accuracy of about 14cm based on the standard deviation of low-slope crossover differences and a precision of close to 2cm for the Laser 2a, Release 21, GLA12 data used here.

Humbert, Angelika↗

Improved Thin Cirrus and Terminator Cloud Detection in CERES Cloud Mask

Thin cirrus clouds account for about 20-30% of the total cloud coverage and affect the global radiation budget by increasing the Earth's albedo and reducing infrared emissions. Thin cirrus, however, are often underestimated by traditional satellite cloud detection algorithms. This difficulty is caused by the lack of spectral contrast between optically thin cirrus and the surface in techniques that use visible (0.65 micron ) and infrared (11 micron ) channels. In the Clouds and the Earth s Radiant Energy System (CERES) Aqua Edition 1 (AEd1) and Terra Edition 3 (TEd3) Cloud Masks, thin cirrus detection is significantly improved over both land and ocean using a technique that combines MODIS high-resolution measurements from the 1.38 and 11 micron channels and brightness temperature differences (BTDs) of 11-12, 8.5-11, and 3.7-11 micron channels. To account for humidity and view angle dependencies, empirical relationships were derived with observations from the 1.38 micron reflectance and the 11-12 and 8.5-11 micron BTDs using 70 granules of MODIS data in 2002 and 2003. Another challenge in global cloud detection algorithms occurs near the day/night terminator where information from the visible 0.65 micron channel and the estimated solar component of 3.7 micron channel becomes less reliable. As a result, clouds are often underestimated or misidentified near the terminator over land and ocean. Comparisons between the CLAVR-x (Clouds from Advanced Very High Resolution Radiometer [AVHRR]) cloud coverage and Geoscience Laser Altimeter System (GLAS) measurements north of 60 N indicate significant amounts of missing clouds from CLAVR-x because this part of the world was near the day/night terminator viewed by AVHRR. Comparisons between MODIS cloud products (MOD06) and GLAS in the same region also show similar difficulties with MODIS cloud retrievals. The consistent detection of clouds through out the day is needed to provide reliable cloud and radiation products for CERES and other research efforts involving the modeling of clouds and their interaction with the radiation budget.

Trepte, Qing↗

Improvements in Near-Terminator and Nocturnal Cloud Masks using Satellite Image Data over the Atmospheric Radiation Measurement Sites

Cloud detection using satellite measurements presents a big challenge near the terminator where the visible (VIS; 0.65 (micro)m) channel becomes less reliable and the reflected solar component of the solar infrared 3.9-(micro)m channel reaches very low signal-to-noise ratio levels. As a result, clouds are underestimated near the terminator and at night over land and ocean in previous Atmospheric Radiation Measurement (ARM) Program cloud retrievals using Geostationary Operational Environmental Satellite (GOES) imager data. Cloud detection near the terminator has always been a challenge. For example, comparisons between the CLAVR-x (Clouds from Advanced Very High Resolution Radiometer (AVHRR)) cloud coverage and Geoscience Laser Altimeter System (GLAS) measurements north of 60 degrees N indicate significant amounts of missing clouds from AVHRR because this part of the world was near the day/night terminator viewed by AVHRR. Comparisons between MODIS cloud products and GLAS at the same regions also shows the same difficulty in the MODIS cloud retrieval (Pavolonis and Heidinger 2005). Consistent detection of clouds at all times of day is needed to provide reliable cloud and radiation products for ARM and other research efforts involving the modeling of clouds and their interaction with the radiation budget. To minimize inconsistencies between daytime and nighttime retrievals, this paper develops an improved twilight and nighttime cloud mask using GOES-9, 10, and 12 imager data over the ARM sites and the continental United States (CONUS).

Detection↗

Meeting Report: Long Term Monitoring of Global Vegetation using Moderate Resolution Satellites

The international community has long recognized the need to coordinate observations of Earth from space. In 1984, this situation provided the impetus for creating the Committee on Earth Observation Satellites (CEOS), an international coordinating mechanism charged with coordinating international civil spaceborne missions designed to observe and study planet Earth. Within CEOS, its Working Group on Calibration and Validation (WGCV) is tasked with coordinating satellite-based global observations of vegetation. Currently, several international organizations are focusing on the requirements for Earth observation from space to address key science questions and societal benefits related to our terrestrial environment. The Global Vegetation Workshop, sponsored by the WGCV and held in Missoula, Montana, 7-10 August, 2006, was organized to establish a framework to understand the inter-relationships among multiple, global vegetation products and identify opportunities for: 1) Increasing knowledge through combined products, 2) Realizing efficiency by avoiding redundancy, and 3) Developing near- and long-term plans to avoid gaps in our understanding of critical global vegetation information. The Global Vegetation Workshop brought together 135 researchers from 25 states and 14 countries to advance these themes and formulate recommendations for CEOS members and the Global Earth Observation System of Systems (GEOSS). The eighteen oral presentations and most of the 74 posters presented at the meeting can be downloaded from the meeting website (www.ntsg.umt.edu/VEGMTG/). Meeting attendees were given a copy of the July 2006 IEEE Transactions on Geoscience and Remote Sensing Special Issue on Global Land Product Validation, coordinated by the CEOS Working Group on Calibration and Validation (WGCV). This issue contains 29 articles focusing on validation products from several of the sensors discussed during the workshop.

Morisette, Jeffrey↗

Initial Verification of GEOS-4 Aerosols Using CALIPSO and MODIS: Scene Classification

A-train sensors such as MODIS and MISR provide column aerosol properties, and in the process a means of estimating aerosol type (e.g. smoke vs. dust). Correct classification of aerosol type is important because retrievals are often dependent upon selection of the right aerosol model. In addition, aerosol scene classification helps place the retrieved products in context for comparisons and analysis with aerosol transport models. The recent addition of CALIPSO to the A-train now provides a means of classifying aerosol distribution with altitude. CALIPSO level 1 products include profiles of attenuated backscatter at 532 and 1064 nm, and depolarization at 532 nm. Backscatter intensity, wavelength ratio, and depolarization provide information on the vertical profile of aerosol concentration, size, and shape. Thus similar estimates of aerosol type using MODIS or MISR are possible with CALIPSO, and the combination of data from all sensors provides a means of 3D aerosol scene classification. The NASA Goddard Earth Observing System general circulation model and data assimilation system (GEOS-4) provides global 3D aerosol mass for sulfate, sea salt, dust, and black and organic carbon. A GEOS-4 aerosol scene classification algorithm has been developed to provide estimates of aerosol mixtures along the flight track for NASA's Geoscience Laser Altimeter System (GLAS) satellite lidar. GLAS launched in 2003 and did not have the benefit of depolarization measurements or other sensors from the A-train. Aerosol typing from GLAS data alone was not possible, and the GEOS-4 aerosol classifier has been used to identify aerosol type and improve the retrieval of GLAS products. Here we compare 3D aerosol scene classification using CALIPSO and MODIS with the GEOS-4 aerosol classifier. Dust, smoke, and pollution examples will be discussed in the context of providing an initial verification of the 3D GEOS-4 aerosol products. Prior model verification has only been attempted with surface mass comparisons and column optical depth from AERONET and MODIS.

Welton, Ellsworth J.↗

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↗

Timber Volume and Biomass Estimates in Central Siberia from Satellite Data

Mapping of boreal forest's type, structure parameters and biomass are critical for understanding the boreal forest's significance in the carbon cycle, its response to and impact on global climate change. The biggest deficiency of the existing ground based forest inventories is the uncertainty in the inventory data, particularly in remote areas of Siberia where sampling is sparse, lacking, and often decades old. Remote sensing methods can help overcome these problems. In this joint US and Russian study, we used the moderate resolution imaging spectroradiometer (MODIS) and unique waveform data of the geoscience laser altimeter system (GLAS) and produced a map of timber volume for a 10degx12deg area in Central Siberia. Using these methods, the mean timber volume for the forested area in the total study area was 203 m3/ ha. The new remote sensing methods used in this study provide a truly independent estimate of forest structure, which is not dependent on traditional ground forest inventory methods.

Ranson, K. Jon↗

Efficient Kriging via Fast Matrix-Vector Products

Interpolating scattered data points is a problem of wide ranging interest. Ordinary kriging is an optimal scattered data estimator, widely used in geosciences and remote sensing. A generalized version of this technique, called cokriging, can be used for image fusion of remotely sensed data. However, it is computationally very expensive for large data sets. We demonstrate the time efficiency and accuracy of approximating ordinary kriging through the use of fast matrixvector products combined with iterative methods. We used methods based on the fast Multipole methods and nearest neighbor searching techniques for implementations of the fast matrix-vector products.

Memarsadeghi, Nargess↗

Performance of the GLAS Laser Transmitter in Space

The Geoscience Laser Altimeter System (GLAS), launched in January 2003, is a laser altimeter and lidar for the Earth Observing System's (EOS) ICESat mission. The laser transmitter requirements, design and qualification test results and in-flight performance for this space-based remote sensing instrument is summarized and presented.

Yu, Anthony W.↗

Performance of the Lidar Design and Data Algorithms for the GLAS Global Cloud and Aerosol Measurements

The Geoscience Laser Altimeter System (GLAS) launched in early 2003 is the first polar orbiting satellite lidar. The instrument design includes high performance observations of the distribution and optical scattering cross sections of atmospheric clouds and aerosol. The backscatter lidar operates at two wavelengths, 532 and 1064 nm. For the atmospheric cloud and aerosol measurements, the 532 nm channel was designed for ultra high efficiency with solid state photon counting detectors and etalon filtering. Data processing algorithms were developed to calibrate and normalize the signals and produce global scale data products of the height distribution of cloud and aerosol layers and their optical depths and particulate scattering cross sections up to the limit of optical attenuation. The paper will concentrate on the effectiveness and limitations of the lidar channel design and data product algorithms. Both atmospheric receiver channels meet and exceed their design goals. Geiger Mode Avalanche Photodiode modules are used for the 532 nm signal. The operational experience is that some signal artifacts and non-linearity require correction in data processing. As with all photon counting detectors, a pulse-pile-up calibration is an important aspect of the measurement. Additional signal corrections were found to be necessary relating to correction of a saturation signal-run-on effect and also for daytime data, a small range dependent variation in the responsivity. It was possible to correct for these signal errors in data processing and achieve the requirement to accurately profile aerosol and cloud cross section down to 10-7 llm-sr. The analysis procedure employs a precise calibration against molecular scattering in the mid-stratosphere. The 1064 nm channel detection employs a high-speed analog APD for surface and atmospheric measurements where the detection sensitivity is limited by detector noise and is over an order of magnitude less than at 532 nm. A unique feature of the GLAS is a full acquisition of the surface return pulse, which has important application to the atmospheric transmission retrieval.

Spinhirne, James D.↗

Abstracts of the Annual Meeting of Planetary Geologic Mappers, Flagstaff, AZ, 2008

Topics discussed include: Merging of the USGS Atlas of Mercury 1:5,000,000 Geologic Series; Geologic Mapping of the V-36 Thetis Regio Quadrangle: 2008 Progress Report; Structural Maps of the V-17 Beta Regio Quadrangle, Venus; Geologic Mapping of Isabella Quadrangle (V-50) and Helen Planitia, Venus; Renewed Mapping of the Nepthys Mons Quadrangle (V-54), Venus; Mapping the Sedna-Lavinia Region of Venus; Geologic Mapping of the Guinevere Planitia Quadrangle of Venus; Geological Mapping of Fortuna Tessera (V-2): Venus and Earth's Archean Process Comparisons; Geological Mapping of the North Polar Region of Venus (V-1 Snegurochka Planitia): Significant Problems and Comparisons to the Earth's Archean; Venus Quadrangle Geological Mapping: Use of Geoscience Data Visualization Systems in Mapping and Training; Geologic Map of the V-1 Snegurochka Planitia Quadrangle: Progress Report; The Fredegonde (V-57) Quadrangle, Venus: Characterization of the Venus Midlands; Formation and Evolution of Lakshmi Planum (V-7), Venus: Assessment of Models using Observations from Geological Mapping; Geologic Map of the Meskhent Tessera Quadrangle (V-3), Venus: Evidence for Early Formation and Preservation of Regional Topography; Geological Mapping of the Lada Terra (V-56) Quadrangle, Venus: A Progress Report; Geology of the Lachesis Tessera Quadrangle (V-18), Venus; Geologic Mapping of the Juno Chasma Quadrangle, Venus: Establishing the Relation Between Rifting and Volcanism; Geologic Mapping of V-19, V-28, and V-53; Lunar Geologic Mapping Program: 2008 Update; Geologic Mapping of the Marius Quadrangle, the Moon; Geologic Mapping along the Arabia Terra Dichotomy Boundary: Mawrth Vallis and Nili Fossae, Mars: Introductory Report; New Geologic Map of the Argyre Region of Mars; Geologic Evolution of the Martian Highlands: MTMs -20002, -20007, -25002, and -25007; Mapping Hesperia Planum, Mars; Geologic Mapping of the Meridiani Region, Mars; Geology of Holden Crater and the Holden and Ladon Multi-Ring Impact Basins, Margaritifer Terra, Mars; Geologic Mapping of Athabasca Valles; Geologic Mapping of MTM -30247, -35247 and -40247 Quadrangles, Reull Vallis Region of Mars; Geologic Mapping of the Martian Impact Crater Tooting; Geology of the Southern Utopia Planitia Highland-Lowland Boundary Plain: First Year Results and Second Year Plan; Mars Global Geologic Mapping: Amazonian Results; Recent Geologic Mapping Results for the Polar Regions of Mars; Geologic Mapping of the Medusae Fossae Formation on Mars (MC-8 SE and MC-23 NW) and the Northern Lowlands of Venus (V-16 and V-15); Geologic Mapping of the Zal, Hi'iaka, and Shamshu Regions of Io; Global Geologic Map of Europa; Material Units, Structures/Landforms, and Stratigraphy for the Global Geologic Map of Ganymede (1:15M); and Global Geologic Mapping of Io: Preliminary Results.

Bleamaster, Leslie F., III↗

Using MODIS and GLAS Data to Develop Timber Volume Estimates in Central Siberia

The boreal forest is the Earth's largest terrestrial biome, covering some 12 million km2 and accounting for about one third of this planet's total forest area. Mapping of boreal forest's type, structure parameters and biomass are critical for understanding the boreal forest's significance in the carbon cycle, its response to and impact on global climate change. Ground based forest inventories, have much uncertainty in the inventory data, particularly in remote areas of Siberia where sampling is sparse and/or lacking. In addition, many of the forest inventories that do exist for Siberia are now a decade or more old. Thus, available forest inventories fail to capture the current conditions. Changes in forest structure in a particular forest-type and region can change significantly due to changing environment conditions, and natural and anthropogenic disturbance. Remote sensing methods can potentially overcome these problems. Multispectral sensors can be used to provide vegetation cover maps that show a timely and accurate geographic distribution of vegetation types rather than decade old ground based maps. Lidar sensors can be used to directly obtain measurements that can be used to derive critical forest structure information (e.g., height, density, and volume). These in turn can used to estimate biomass components using allometric equations without having to use out dated forest inventory. Finally, remote sensing data is ideally suited to provide a sampling basis for a rigorous statistical estimate of the variance and error bound on forest structure measures. In this study, new remote sensing methods were applied to develop estimates timber volume using NASA's MODerate resolution Imaging Spectroradiometer (MODIS) and unique waveform data of the geoscience laser altimeter system (GLAS) for a 10 deg x 10 deg area in central Siberia. Using MODIS and GLAS data, maps were produced for cover type and timber volume for 2003, and a realistic variance (error bound) for timber volume was calculated for the study area. In this 'study we used only GLAS footprints that had a slope value of less than 10 deg. This was done to avoid large errors due to the effect of slope on the GLAS models. The method requires the integration of new remote sensing methods with available ground studies of forest timber volume conducted in Russian forests. The results were compared to traditional ground forest inventory methods reported in the literature and to ground truth collected in the study area.

Ranson, K. Jon↗

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↗