Search NASA⌕ Search

SEARCH · Search NASA

Results for “statistical image 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 973 records · Page 54

Combining human and computer interpretation capabilities to analyze ERTS imagery

The human photointerpreter and the computer have complementary capabilities that are exploited in a computer-based data analysis system developed at the Forestry Remote Sensing Laboratory, University of California. This system is designed to optimize the process of extracting resource information from ERTS images. The human has the ability to quickly delineate gross differences in land classes, such as wildland, urban, and agriculture on appropriate ERTS images, and to further break these gross classes into meaningful subclasses. The computer, however, can more efficiently analyze point-by-point spectral information and localized textural information which can result in a much more detailed agricultural or wildland classification based on species composition and/or plant association. These human and computer capabilities have been integrated through the use of an inexpensive small scale computer dedicated to the interactive preprocessing of the human inputs and the display of raw ERTS images and computer classified images. The small computer is linked to a large scale computer system wherein the bulk of the statistical work and the automatic point-by-point classification is done.

Nichols, J. D.↗

Space vehicle approach velocity judgments under simulated visual space conditions

There were 35 volunteers who responded when they first perceived an increase in apparent size of a collimated, two-dimensional perspective image of an Orbiter vehicle. The variables of interest included the presence (or absence) of a fixed reticle within the field of view (FOV), background starfield velocity, initial range to the vehicle and vehicle closure velocity. It was found that: 1) increasing vehicle approach velocity yielded a very small (but significant) effect of faster detection of vehicle movement, nevertheless, response variability was relatively large; 2) including the fixed reticle in the FOV produced significantly slower detection of vehicle radial movement, however this occurred only at the largest range and the magnitude of the effect was only about 15% of the one sigma value; and 3) increasing background star velocity during this judgment led to slower detection of vehicle movement. While statistically significant, this effect was small and occurred primarily at the largest range. A possible explanation for the last two findings is that other static and dynamic objects within the visual field may compete for available attention which otherwise would be available for judging image expansion; thus, the target's retinal image has to expand more than otherwise for its movement to be detected. This study also showed that the Proximity Operations Research Mockup at NASA/Ames can be used effectively to investigate a variety of visual judgment questions related to future space operations. These findings are discussed in relation to previous research and possible underlying mechanisms.

Space Flight↗

Robotic NDA of Holdup Deposits in Gaseous Diffusion Piping by Gamma Assay of Uranium-235 - 18331

This research exploits in-pipe access to develop robotically deployed, high-cadence, accurate, certain, paperless assay of U-235 in holdup deposits within process piping. Prior radiometric assay from outside the pipes suffered from manual deployment challenges, attenuation of detection through pipe walls, long counting times, approximate modeling, and shortfalls associated with transcription and human interpretation. These downsides limited the speed, quality, and economy of assay resulting in vast cost and budget consequences to D and D. D and D of outdated facilities cuts into these pipes as part of the demolition process, creating the unique, previously unexploited possibility of robotic in-pipe assay. The assay collimates radiation emanating from all but a short segment of pipe wall from reaching a gamma detector. Hence, the detector only views and measures source from a short segment of pipe at a given time. This is achieved by an innovative pair of collimating discs that are coaxial with the pipe and positioned fore and aft of the detector. The detector assembly is translated through pipes by an autonomous mobile robot. Beyond radiometric assay of U-235, the robot images and geometrically models the pipe interior and deposit appearance. The robot is recovered from the same pipe opening from which it is launched, hence it drives the same distance out and back, measuring the same deposits twice. This achieves redundant radiometric and odometric measurements which adds further to statistical significance of the method. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Experimental and Computational Evaluation of Lipidomic In-Source Fragmentation as a Result of Postionization with Matrix-Assisted Laser Desorption/Ionization

Matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI) can provide spatially resolved molecular information about a sample. Recently, a postionization approach (MALDI-2) has been commercially integrated with MALDI-MSI, allowing for bettered sensitivity and consequent improved spatial resolution. While advantages of MALDI-2 have previously been established, we demonstrate here statistically increased in-source fragmentation (ISF) results from postionization with a commercial instrument. Via lipid standard analyses, known MALDI ISF pathways (e.g., loss of trimethylamine) were statistically increased in MALDI-2 compared to MALDI-1 (65–172% increase in fragmentation). Gas phase molecular modeling with density functional theory estimated that the most-weighted virtual orbitals to excite within lipids involve ester and phosphate bonds. Protonated lipid excitation energies are furthermore red-shifted compared to those of other adduct types [e.g., 254 nm for protonated PC(16:0/18:1)] and approach the MALDI-2 laser energy (266 nm). Analysis of rat brain homogenate detected statistically more positive-ion mode peaks with MALDI-2 (1090) than that with MALDI-1 (719), where Kernel density estimations showed that the majority of this enhancement occurs with low m/z ions (i.e., m/z 75–500). Taken together with the lipid standard data, these observations may indicate ISF due to postionization. Finally, while artifact contributions from matrix blanks were also noted, both experimental and computational data sets suggest that the overall extent of ISF is statistically increased in MALDI-2 compared to MALDI-1.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Solar X-Ray Limb

We describe a new technique to measure the height of the X-ray limb with observations from occulted X-ray flare sources as observed by the RHESSI (the Reuven Ramaty High-Energy Spectroscopic Imager) satellite. This method has model dependencies different from those present in traditional observations at optical wavelengths, which depend upon detailed modeling involving radiative transfer in a medium with complicated geometry and flows. It thus provides an independent and more rigorous measurement of the “true” solar radius, which means that of the mass distribution. RHESSI’s measurement makes use of the flare X-ray source’s spatial Fourier components (the visibilities), which are sensitive to the presence of the sharp edge at the lower boundary of the occulted source. We have found a suitable flare event for analysis, SOL2011-10-20T03:25 (M1.7), and report a first result from this novel technique here. Using a four-minute integration over the 3–25 keV photon energy range, we find R{sub X-ray}=960.11±0.15±0.29 arcsec, at 1 au, where the uncertainties include statistical uncertainties from the method and a systematic error. The standard VAL-C model predicts a value of 959.94 arcsec, which is about 1σ below our value.

79 ASTRONOMY AND ASTROPHYSICS↗

Application of Support Vector Regression to Derive Crater Depth/Diameter From Satellite Images

Through the study of impact crater shapes, one can draw important conclusions about the nature and evolution of planetary surfaces [e.g., 1-4].In particular, studying the depth (d) to diameter (D)ratio (d/D) of a population of impact craters, in combination with crater count statistics, can yield valuable insights regarding rates of erosion and burial[5]. Motivated by the great abundance of available planetary surface image data, the goal of this project is to develop an efficient way to estimate d/D from satellite images of impact craters for which stereo information is not available [6]. We set out to develop and train a machine learning algorithm to extract d/D from a dataset of synthetic impact crater images for which model d/D is known. The applications of machine learning to planetary science are numerous and diverse [7], including automatic planetary surface mapping [8] and the detection of impact craters [9]. Our algorithm makes use of Support Vector Regression (SVR), which is a type of Support Vector Machine (SVM) [10, 11].SVMs are a branch of supervised machine learning valued for their straightforward implementation and versatility in solving both classification and regression problems. In regression analysis, an SVR algorithm produces a hyperplane function to fit the training data points, as well as an ε-tube that surrounds the hyperplane. Tunable hyperparameters include the width of the ε-tube (ε) and the amount an algorithm is penalized for points which fall outside the ε-tube.

L R Chin↗

SU-G-BRC-13: Model Based Classification for Optimal Position Selection for Left-Sided Breast Radiotherapy: Free Breathing, DIBH, Or Prone

Purpose: There are clinical decision challenges to select optimal treatment positions for left-sided breast cancer patients—supine free breathing (FB), supine Deep Inspiration Breath Hold (DIBH) and prone free breathing (prone). Physicians often make the decision based on experiences and trials, which might not always result optimal OAR doses. We herein propose a mathematical model to predict the lowest OAR doses among these three positions, providing a quantitative tool for corresponding clinical decision. Methods: Patients were scanned in FB, DIBH, and prone positions under an IRB approved protocol. Tangential beam plans were generated for each position, and OAR doses were calculated. The position with least OAR doses is defined as the optimal position. The following features were extracted from each scan to build the model: heart, ipsilateral lung, breast volume, in-field heart, ipsilateral lung volume, distance between heart and target, laterality of heart, and dose to heart and ipsilateral lung. Principal Components Analysis (PCA) was applied to remove the co-linearity of the input data and also to lower the data dimensionality. Feature selection, another method to reduce dimensionality, was applied as a comparison. Support Vector Machine (SVM) was then used for classification. Thirtyseven patient data were acquired; up to now, five patient plans were available. K-fold cross validation was used to validate the accuracy of the classifier model with small training size. Results: The classification results and K-fold cross validation demonstrated the model is capable of predicting the optimal position for patients. The accuracy of K-fold cross validations has reached 80%. Compared to PCA, feature selection allows causal features of dose to be determined. This provides more clinical insights. Conclusion: The proposed classification system appeared to be feasible. We are generating plans for the rest of the 37 patient images, and more statistically significant results are to be presented.

60 APPLIED LIFE SCIENCES↗

Exploration of thermal counterflow in He II using particle tracking velocimetry

Flow visualization using particle image velocimetry (PIV) and particularly particle tracking velocimetry (PTV) has been applied to thermal counterflow in He II for nearly two decades now, but the results remain difficult to interpret because tracer particle motion can be influenced by both the normal fluid and superfluid components of He II as well as the quantized vortex tangle. For instance, in one early experiment it was observed (using PTV) that tracer particles move at the normal fluid velocity v n , while in another it was observed (using PIV) that particles move at v n /2. Besides the different visualization methods, the range of applied heat flux investigated by these experiments differed by an order of magnitude. To resolve this apparent discrepancy and explore the statistics of particle motion in thermal counterflow, we apply the PTV method to a wide range of heat flux at a number of different fluid temperatures. In our analysis, we introduce a scheme for analyzing the velocity of particles presumably moving with the normal fluid separately from those presumably influenced by the quantized vortex tangle. Our results show that for lower heat flux there are two distinct peaks in the streamwise particle velocity probability density function (PDF), with one centered at the normal fluid velocity v n (named G2 for convenience) while the other is centered near v n /2 (G1). For higher heat flux there is a single peak centered near v n /2 (G3). Using our separation scheme, we show quantitatively that there is no size difference between the particles contributing to G1 and G2. We also show that nonclassical features of the transverse particle velocity PDF arise entirely from G1, while the corresponding PDF for G2 exhibits the classical Gaussian form. The G2 transverse velocity fluctuation, backed up by second sound attenuation in decaying counterflow, suggests that large-scale turbulence in the normal fluid is absent from the two-peak region. We offer a brief discussion of the physical mechanisms that may be responsible for our observations, revealing that G1 velocity fluctuations may be linked to fluctuations of quantized vortex line velocity, and suggest a number of numerical simulations that may reveal the underlying physics in detail.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Depth-diameter ratios for Martian impact craters: Implications for target properties and episodes of degradation

This study determines crater depth through use of photoclinometric profiles. Random checks of the photoclinometric results are performed using shadow estimation techniques. The images are Viking Orbiter digital format frames; in cases where the digital image is unusable for photoclinometric analysis, shadow estimation is used to determine crater depths. The two techniques provide depth results within 2 percent of each other. Crater diameters are obtained from the photoclinometric profiles and checked against the diameters measured from the hard-copy images using a digitizer. All images used in this analysis are of approximately 40 m/pixel resolution. The sites that have been analyzed to date include areas within Arabia, Maja Valles, Memnonia, Acidalia, and Elysium. Only results for simple craters (craters less than 5 km in diameter) are discussed here because of the low numbers of complex craters presently measured in the analysis. General results indicate that impact craters are deeper than average. A single d/D relationship for fresh impact craters on Mars does not exist due to changes in target properties across the planet's surface. Within regions where target properties are approximately constant, however, d/D ratios for fresh craters can be determined. In these regions, the d/D ratios of nonpristine craters can be compared with the fresh crater d/D relationship to obtain information on relative degrees of crater degradation. This technique reveals that regional episodes of enhanced degradation have occurred. However, the lack of statistically reliable size-frequency distribution data prevents comparison of the relative ages of these events between different regions, and thus determination of a large-scale episode (or perhaps several episodes) cannot be made at this time.

Barlow, N. G.↗

Preliminary Optical And Electric Field Pulse Statistics From Storm Overflights During The Altus Cumulus Electrification Study

The Altus Cumulus Electrification Study (ACES) was conducted during the month of August, 2002 in an area near Key West, Florida. One of the goals of this uninhabited aerial vehicle (UAV) study was to collect high resolution optical pulse and electric field data from thunderstorms. During the month long campaign, we acquired 5294 lightning generated optical pulses with associated electric field changes. Most of these observations were made while close to the top of the storms. We found filtered mean and median 10-10% optical pulse widths of 875 and 830 microns respectively while the 50-50% mean and median optical pulse widths are 422 and 365 microns respectively. These values are similar to previous results as are the 10-90% mean and median rise times of 327 and 265 microns. The peak electrical to optical pulse delay mean and median were 209 and 145 microns which is longer than one would expect from theoretical results. The results of the pulse analysis will contribute to further validation of the Optical Transient Detector (OTD) and the Lightning Imaging Sensor (LIS) satellites. Pre-launch estimates of the flash detection efficiency were based on a small sample of optical pulse measurements associated with less than 350 lightning discharges collected by NASA U-2 aircraft in the early 1980s. Preliminary analyses of the ACES measurements show that we have greatly increased the number of optical pulses available for validation of the LIS and other orbital lightning optical sensors. Since the Altus was often close to the cloud tops, many of the optical pulses are from low-energy pulses. From these low-energy pulses, we can determine the fraction of optical lightning pulses below the thresholds of LIS, OTD, and any future satellite-based optical sensors such as the geostationary Lightning Mapping Sensor.

Mach, D. A.↗

Bhutan Agriculture: Developing a Crop Mask for Rice and Creating a Data Collection Protocol Utilizing Remotely Sensed Data in Bhutan

Rice cultivation in Bhutan has been increasingly threatened by deteriorating soil health and outbreaks of diseases and pests associated with the global change in climate patterns. Field surveys, which the national government of Bhutan has relied on to monitor remote agricultural lands, are becoming increasingly overwhelmed by growing threats to agricultural health. To address these concerns, NASA DEVELOP partnered with the Department of Agriculture of Bhutan, the Bhutan Foundation, and the Ugyen Wangchuck Institute of Conservation and Environmental Research (UWICER) and worked to increase the government of Bhutan’s agricultural monitoring capacity. Utilizing Earth observations including Landsat 8 Operational Land Imager (OLI), Sentinel-1 C-band Synthetic Aperture Radar (C-SAR), Shuttle Radar Topography Mission (SRTM), and Planet imagery, the DEVELOP team worked with NASA SERVIR and created a sampling protocol to identify rice plantations and supplement field surveys for more efficient agriculture monitoring. The analysis focused on districts Paro, Punakha, Samtse, Sarpang, Trongsa, Zhemgang, Wangdue Phodrang, and Samdrup Jongkhar in the year 2020 during the period of transplantation (June) to harvesting of rice (November). The team provided the partners with a sampling protocol for integrating NASA Earth observations into their crop monitoring methods, as well as a crop mask for rice identification and to aid crop management. The crop mask for rice was developed using the Random Forest (RF) classifier for the eight districts of Bhutan. Visually, the random forest model has proved to be more accurate and precise than the classification and Regression Tree model. Statistically, the Random Forest model was 91.8% accurate in identifying rice in Bhutan.

Yeshey Seldon↗

The dilemma of parotid gland and pharyngeal constrictor muscles preservation—Is daily online image guidance required? A dosimetric analysis

With margin reduction common in head and neck radiotherapy, it is critical that the dosimetric effects of setup deviations are quantified. With past studies focusing on the quantification of positional and volumetric changes of organs at risk (OARs), this study aimed to measure the dose delivered to these the parotid gland (PG) and pharyngeal constrictor muscles (PCMs) using cone beam computed tomography (CBCT). Furthermore, this investigation sought to establish a potential time trend of change in dose delivered to target volumes secondary to ascertaining the need for daily image guidance (IG) to reduce the dose burden to these important OARs. Intensity modulated radiotherapy (IMRT) plans for 5 locally advanced head and neck patients' plans were created and mapped to weekly CBCTs. Each plan was recalculated without heterogeneity correction allowing for dosimetric comparison. Dosimetric endpoints recorded to assess the effect of positional variation were as per ICRU 83 and included D{sub 95} and D{sub 98} for the target volumes, mean dose (MD) and V{sub 30} {sub Gy} for the PGs, and V{sub 50} {sub Gy} and MD for the PCMs. Results were deemed statistically significant if p < 0.05. No significant time trends were established for these OARs. A significant decrease in V{sub 50} {sub Gy} was observed for all PCMs (p < 0.001) on all CBCTs relative to the original plan. Regarding target volumes, a highly significant decrease in MD (MD = 20 Gy, CI: −20.310 to −19.820) in D{sub 98} of the high-dose planning target volume (PTV [70 Gy]; PTVD{sub 98%} = 70 Gy) for case 3 was found (p ≤ 0.001). A nonpredictable, yet significant dosimetric effect was found. A clinically acceptable balance must be achieved between OAR dosimetry and target coverage as can be achieved by frequent IG.

61 RADIATION PROTECTION AND DOSIMETRY↗

Temperature and magnetic-field driven dynamics in artificial magnetic square ice

Artificial spin ices are often spoken of as being realisations of some of the celebrated vertex models of statistical mechanics, where the exact microstate of the system can be imaged using advanced magnetic microscopy methods. The fact that a stable image can be formed means that the system is in fact athermal and not undergoing the usual finite-temperature fluctuations of a statistical mechanical system. In this paper we report on the preparation of artificial spin ices with islands that are thermally fluctuating due to their very small size. The relaxation rate of these islands was determined using variable frequency focused magneto-optic Kerr measurements. We performed magnetic imaging of artificial spin ice under varied temperature and magnetic field using X-ray transmission microscopy which uses X-ray magnetic circular dichroism to generate magnetic contrast. Furthermore, we have developed an on-membrane heater in order to apply temperatures in excess of 700 K and have shown increased dynamics due to higher temperature. Due to the ‘photon-in, photon-out' method employed here, it is the first report where it is possible to image the microstates of an ASI system under the simultaneous application of temperature and magnetic field, enabling the determination of relaxation rates, coercivties, and the analysis of vertex population during reversal.

36 MATERIALS SCIENCE↗

High-count-rate effects in event processing for the XRISM/Resolve X-ray microcalorimeter. II. Energy scale and resolution in orbit

The Resolve instrument on the X-ray Imaging and Spectroscopy Mission (XRISM) uses a 36 pixel microcalorimeter designed to deliver high-resolution, non-dispersive X-ray spectroscopy. Although it is optimized for extended sources with low count rates, Resolve observations of bright point sources are still able to provide unique insights into the physics of these objects, as long as high-count-rate effects are addressed in the analysis. These effects include the loss of exposure time for each pixel, changes in the energy scale, and changes in the energy resolution. To investigate these effects under realistic observational conditions, we observed the bright X-ray source, the Crab Nebula, with XRISM at several offset positions with respect to the Resolve field of view and with continuous illumination from 55 Fe sources on the filter wheel. For the spectral analysis, we excluded data where exposure-time loss was too significant to ensure reliable spectral statistics. The energy scale at 6 keV shows a slight negative shift in the high-count-rate regime. The energy resolution at 6 keV worsens as the count rate in electrically neighboring pixels increases, but can be restored by applying a nearest-neighbor coincidence cut (“cross-talk cut”). We examined how these effects influence the observation of bright point sources, using GX 13+1 as a test case, and identified an eV-scale energy offset at 6 keV between the inner (brighter) and outer (fainter) pixels. Users who seek to analyze velocity structures on the order of tens of km s–1 should account for such high-count-rate effects. These findings will aid in the interpretation of Resolve data from bright sources and provide valuable considerations for designing and planning for future microcalorimeter missions.

X-rays: general↗

Workshop to Identify Opportunities in Biological and Environmental Research Uniquely Enabled by the APS Upgrade

The development of synchrotron-based x-ray science has enabled tremendous scientific advances, driven both by discoveries of new fundamental ways to derive insight from the interaction of x-rays with matter and by the evolution of x-ray light sources themselves. Over the last 20-plus years, the U.S. Department of Energy’s (DOE’s) Advanced Photon Source (APS) at Argonne has yielded a multitude of transformational insights into science questions across numerous areas relevant to the biological, geological, geochemical, biogeochemical, and environmental sciences. Recent developments in accelerator technology, particularly the multi-bend achromat lattice, have led to plans to dramatically upgrade the APS, allowing an upgraded APS to offer orders-of-magnitude improvements in key parameters such as coherence and brightness, while retaining high-energy capabilities. These enhanced parameters are directly relevant to almost all hard x-ray imaging, scattering, and spectroscopy experiments (Eriksson et al., 2014). The APS Upgrade Project exceeds the capabilities of today’s synchrotrons by 2 to 3 orders of magnitude in brightness and coherent flux in the hard x-ray range, enabling a transformational range of new probes for the structure, properties, and functionality of matter. With vastly improved coherent flux, nearly all techniques can become microscopies. For example, frontiers inaccessible today will be opened by x-ray instruments that provide three-dimensional (3-D) real-space resolution down to nanometers across, up to millimeters of field of view, providing the ability to follow in real time relevant processes and trace element sensitivity down to a few atoms. In addition to the unprecedented spatial resolution across a large field of view, and with appropriate time resolution, the ability for high-throughput, four-dimension/five-dimension data collection enables unprecedented statistical analysis of complex biological and environmental systems in changing environments.

54 ENVIRONMENTAL SCIENCES↗

Segmentation and Classification of Fission as Pores in Reactor Irradiated Annular U–10Zr Metallic Fuel Using Machine Learning Models

Metallic fuels, particularly U—10Zr, are promising candidates for next-generation sodium-cooled fast reactors. Irradiation of nuclear fuels in reactors can lead to the formation of solid and gas fission product which subsequently forms microstructural pores, deteriorating fuel performance. Due to the massive amount of pores and complex phases formed, a quantitative description of fission gas pores is not yet available, preventing the development of microstructure-informed fuel performance modeling for fuel qualification. This paper applied a pre-trained deep learning model to ~10,260 high magnification scanning electron microscopy images. This method increased the accuracy of fission gas pore segmentation and allows statistical features to be extracted which cannot be achieved manually. A pre-trained decision tree model worked on the segemenation results and further classified the pores into different categories to produce a correlation between the pores, movement of lanthanides, and temperature gradient during irradiation. Finally, this paper emphasizes the potentials of machine learning models to accelerate fuel research, development, and qualification for advanced reactors.

36 MATERIALS SCIENCE↗

Blob-hole correlation model for edge turbulence and comparisons with NSTX GPI data

Gas puff imaging (GPI) observations made in NSTX [Zweben S J, et al., 2017 Phys. Plasmas 24 102509] have revealed two-point spatial correlations of edge and scrape-off layer turbulence in the plane perpendicular to the magnetic field. A common feature is the occurrence of dipole-like patterns with significant regions of negative correlation. In this paper, we explore the possibility that these dipole patterns may be due to blob-hole pairs. Statistical methods are applied to determine the two-point spatial correlation that results from a model of blob-hole pair formation. It is shown that the model produces dipole correlation patterns that are qualitatively similar to the GPI data in several respects. Effects of the reference location (confined surfaces or scrape-off layer), a superimposed random background, hole velocity and lifetime, and background sheared flows are explored and discussed with respect to experimental observations. Additional analysis of the experimental GPI dataset is performed to further test this blob-hole correlation model. A time delay two-point spatial correlation study did not reveal inward propagation of the negative correlation structures that were postulated to correspond to holes in the data nor did it suggest that the negative correlation structures are due to neutral shadowing. However, tracing of the highest and lowest values (extrema) of the normalized GPI fluctuations shows strong evidence for mean inward propagation of minima and outward propagation of maxima, in qualitative agreement with theoretical expectations. Other properties of the experimentally observed extrema are discussed.

Myra, J R↗

Experimental philosophy leading to a small scale digital data base of the conterminous United States for designing experiments with remotely sensed data

Research using satellite remotely sensed data, even within any single scientific discipline, often lacked a unifying principle or strategy with which to plan or integrate studies conducted over an area so large that exhaustive examination is infeasible, e.g., the U.S.A. However, such a series of studies would seem to be at the heart of what makes satellite remote sensing unique, that is the ability to select for study from among remotely sensed data sets distributed widely over the U.S., over time, where the resources do not exist to examine all of them. Using this philosophical underpinning and the concept of a unifying principle, an operational procedure for developing a sampling strategy and formal testable hypotheses was constructed. The procedure is applicable across disciplines, when the investigator restates the research question in symbolic form, i.e., quantifies it. The procedure is set within the statistical framework of general linear models. The dependent variable is any arbitrary function of remotely sensed data and the independent variables are values or levels of factors which represent regional climatic conditions and/or properties of the Earth's surface. These factors are operationally defined as maps from the U.S. National Atlas (U.S.G.S., 1970). Eighty-five maps from the National Atlas, representing climatic and surface attributes, were automated by point counting at an effective resolution of one observation every 17.6 km (11 miles) yielding 22,505 observations per map. The maps were registered to one another in a two step procedure producing a coarse, then fine scale registration. After registration, the maps were iteratively checked for errors using manual and automated procedures. The error free maps were annotated with identification and legend information and then stored as card images, one map to a file. A sampling design will be accomplished through a regionalization analysis of the National Atlas data base (presently being conducted). From this analysis a map of homogeneous regions of the U.S.A. will be created and samples (LANDSAT scenes) assigned by region.

Labovitz, M. L.↗