Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data segmentation”

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 955 records · Page 53

Radial widths, optical depths, and eccentricities of the Uranian rings

Observations of the stellar occultation by the Uranian rings of 15/16 August 1980 are used to estimate radial widths and normal optical depths for segments of rings 6, 5, 4, alpha, beta, eta, gamma, and delta. Synthetic occultation profiles are generated to match the observed light curves. A review of published data confirms the existence of width-radius relations for rings alpha and beta, and indicates that the optical depths of these two rings vary inversely with their radial widths. Masses are obtained for rings alpha and beta, on the assumption that differential precession is prevented by their self-gravity. A quantitative comparison of seven epsilon-ring occultation profiles obtained over a period of 3.4 yr reveals a consistent structure, which may reflect the presence of unresolved gaps and subrings.

Nicholson, P. D.↗

Fluoropolymer Films Deposited by Argon Ion-Beam Sputtering of Polytetrafluoroethylene

The FT-IR, XPS and UV spectra of fluoropolymer films (SPTFE-I) deposited by argon ion-beam sputtering of polytetrafluoroethylene (PTFE) were obtained and compared with prior corresponding spectra of fluoropolymer films (SPTFE-P) deposited by argon rf plasma sputtering of PTFE. Although the F/C ratios for SPTFE-I and -P (1.63 and 1.51) were similar, their structures were quite different in that there was a much higher concentration of CF2 groups in SPTFE-I than in SPTFE-P, ca. 61 and 33% of the total carbon contents, respectively. The FT-IR spectra reflect that difference, that for SPTFE-I showing a distinct doublet at 1210 and 1150 per centimeter while that for SPTFE-P presents a broad, featureless band at ca. 1250 per centimeter. The absorbance of the 1210-per centimeter band in SPTFE-I was proportional to the thickness of the film, in the range of 50-400 nanometers. The SPTFE-I was more transparent in the UV than SPTFE-P at comparable thickness. The mechanism for SPTFE-I formation likely involves "chopping off" of oligomeric segments of PTFE as an accompaniment to "plasma" polymerization of TFE monomer or other fluorocarbon fragments generated in situ from PTFE on impact with energetic Ar ions. Data are presented for SPTFE-I deposits and the associated Ar(+) bombarded PTFE targets where a fresh target was used for each run or a single target was used for a sequence of runs.

Golub, Morton A.↗

Aeroelastic Analysis of a Flexible Wing Wind Tunnel Model with Variable Camber Continuous Trailing Edge Flap Design

This paper presents data analysis of a flexible wing wind tunnel model with a variable camber continuous trailing edge flap (VCCTEF) design for drag minimization tested at the University of Washington Aeronautical Laboratory (UWAL). The wind tunnel test was designed to explore the relative merit of the VCCTEF concept for improved cruise efficiency through the use of low-cost aeroelastic model test techniques. The flexible wing model is a 10%-scale model of a typical transport wing and is constructed of woven fabric composites and foam core. The wing structural stiffness in bending is tailored to be half of the stiffness of a Boeing 757-era transport wing while the torsional stiffness is about the same. This stiffness reduction results in a wing tip deflection of about 10% of the wing semi-span. The VCCTEF is a multi-segment flap design having three chordwise camber segments and five spanwise flap sections for a total of 15 individual flap elements. The three chordwise camber segments can be positioned appropriately to create a desired trailing edge camber. Elastomeric material is used to cover the gaps in between the spanwise flap sections, thereby creating a continuous trailing edge. Wind tunnel data analysis conducted previously shows that the VCCTEF can achieve a drag reduction of up to 6.31% and an improvement in the lift-to-drag ratio (L=D) of up to 4.85%. A method for estimating the bending and torsional stiffnesses of the flexible wingUWAL wind tunnel model from static load test data is presented. The resulting estimation indicates that the stiffness of the flexible wing is significantly stiffer in torsion than in bending by as much as 9 to 1. The lift prediction for the flexible wing is computed by a coupled aerodynamic-structural model. The coupled model is developed by coupling a conceptual aerodynamic tool Vorlax with a finite-element model of the flexible wing via an automated geometry deformation tool. Based on the comparison of the lift curve slope, the lift prediction for the rigid wing is in good agreement with the estimated lift coefficients derived from the wind tunnel test data. Due to the movement of the VCCTEF during the wind tunnel test, uncertainty in the lift prediction due to the indicated variations of the VCCTEF deflection is studied. The results show a significant spread in the lift prediction which contradicts the consistency in the aerodynamic measurements, thus suggesting that the indicated variations as measured by the VICON system may not be reliable. The lift prediction of the flexible wing agrees very well with the measured lift curve for the baseline configuration. The computed bending deflection and wash-out twist of the flexible wing also match reasonably well with the aeroelastic deflection measurements. The results demonstrate the validity of the aerodynamic-structural tool for use to analyze aerodynamic performance of flexible wings.

Wind Tunnel↗

A Design Tool for Liquid Rocket Engine Injectors

A practical design tool for the analysis of flowfields near the injector face has been developed and used to analyze the Fastrac engine. The objective was to produce a computational design tool which was detailed enough to predict the interactive effects of injector element impingement angles and points and the momenta of the individual orifice flows. To obtain a model which could be used to simulate a significant number of individual orifices, a homogeneous computational fluid dynamics model was developed. To describe liquid and vapor sub- and super-critical flows, the model included thermal and caloric equations of state which were valid over a wide range of pressures and temperatures. A homogeneous model was constructed such that the local state of the flow was determined directly, i.e. the quality of the flow was calculated. Such a model does not identify drops or their distribution, but it does allow the flow along the injector face and into the acoustic cavity to be predicted. It also allows the film coolant flow to be accurately described. The initial evaluation of the injector code was made by simulating cold flow from an unlike injector element and from a like-on-like overlapping fan (LOL) injector element. The predicted mass flux distributions of these injector elements compared well to cold flow test results. These are the same cold flow tests which serve as the data base for the JANNAF performance prediction codes. The flux distributions 1 inch downstream of the injector face are very similar; the differences were somewhat larger at further distances from the faceplate. Since the cold flow testing did not achieve good mass balances when integrations across the entire fan were made, the CFD simulation appears to be reasonable alternative to future cold flow testing. To simulate the Fastrac, an RP-1/LOX combustion model must be chosen. This submodel must be relatively simple to accomplish three-dimensional, multiphase flow simulations. Single RP-1 pyrolysis and partial oxidation steps were chosen and the combustion was completed with the wet CO mechanism. Soot was also formed with a single global reaction. To validate the combustion submodel, global data from gas generator tests and from subscale motor test were used to predict qualitatively correct mean molecular weights, temperature, and soot levels. Because such tests do not provide general kinetics rates, the methodology is not necessarily appropriate for other than rocket type flows conditions. Soot predictions were made so that radiation heating to the motor walls can be made. These initial studies of the Fastrac were for a small region close to the injector face and chamber wall which included a segment of the acoustic cavity. The region analyzed includes 11 individual orifice holes to represent the LOL elements and the H2 film coolant holes. Typical results of this simulation are shown in Figure 1. At this point the only available test data to verify the predictions are temperatures measured in the acoustic cavity. These temperatures are in reasonable agreement at about 2000R (1111 K). Future work is expected to include improving the computational efficiency or the CFD model and/or using more computer capacity than the single Pentium PC with which these simulations were made.

Farmer, Richard C.↗

Augmented Method to Improve Thermal Data for the Figure Drift Thermal Distortion Predictions of the JWST OTIS Cryogenic Vacuum Test

The JWST Optical Telescope Element (OTE) assembly is the largest optically stable infrared-optimized telescope currently being manufactured and assembled, and is scheduled for launch in 2018. The JWST OTE, including the 18 segment primary mirror, secondary mirror, and the Aft Optics Subsystem (AOS) are designed to be passively cooled and operate near 45K. These optical elements are supported by a complex composite backplane structure. As a part of the structural distortion model validation efforts, a series of tests are planned during the cryogenic vacuum test of the fully integrated flight hardware at NASA JSC Chamber A. The successful ends to the thermal-distortion phases are heavily dependent on the accurate temperature knowledge of the OTE structural members. However, the current temperature sensor allocations during the cryo-vac test may not have sufficient fidelity to provide accurate knowledge of the temperature distributions within the composite structure. A method based on an inverse distance relationship among the sensors and thermal model nodes was developed to improve the thermal data provided for the nanometer scale WaveFront Error (WFE) predictions. The Linear Distance Weighted Interpolation (LDWI) method was developed to augment the thermal model predictions based on the sparse sensor information. This paper will encompass the development of the LDWI method using the test data from the earlier pathfinder cryo-vac tests, and the results of the notional and as tested WFE predictions from the structural finite element model cases to characterize the accuracies of this LDWI method.

Thermal Modeling↗

A High Performance Computing Approach to Tree Cover Delineation in 1-m NAIP Imagery Using a Probabilistic Learning Framework

Tree cover delineation is a useful instrument in deriving Above Ground Biomass (AGB) density estimates from Very High Resolution (VHR) airborne imagery data. Numerous algorithms have been designed to address this problem, but most of them do not scale to these datasets, which are of the order of terabytes. In this paper, we present a semi-automated probabilistic framework for the segmentation and classification of 1-m National Agriculture Imagery Program (NAIP) for tree-cover delineation for the whole of Continental United States, using a High Performance Computing Architecture. Classification is performed using a multi-layer Feedforward Backpropagation Neural Network and segmentation is performed using a Statistical Region Merging algorithm. The results from the classification and segmentation algorithms are then consolidated into a structured prediction framework using a discriminative undirected probabilistic graphical model based on Conditional Random Field, which helps in capturing the higher order contextual dependencies between neighboring pixels. Once the final probability maps are generated, the framework is updated and re-trained by relabeling misclassified image patches. This leads to a significant improvement in the true positive rates and reduction in false positive rates. The tree cover maps were generated for the whole state of California, spanning a total of 11,095 NAIP tiles covering a total geographical area of 163,696 sq. miles. The framework produced true positive rates of around 88% for fragmented forests and 74% for urban tree cover areas, with false positive rates lower than 2% for both landscapes. Comparative studies with the National Land Cover Data (NLCD) algorithm and the LiDAR canopy height model (CHM) showed the effectiveness of our framework for generating accurate high-resolution tree-cover maps.

Segments↗

Time Domain Tool Validation Using ARES I-X Flight Data

The ARES I-X vehicle was launched from NASA's Kennedy Space Center (KSC) on October 28, 2009 at approximately 11:30 EDT. ARES I-X was the first test flight for NASA s ARES I launch vehicle, and it was the first non-Shuttle launch vehicle designed and flown by NASA since Saturn. The ARES I-X had a 4-segment solid rocket booster (SRB) first stage and a dummy upper stage (US) to emulate the properties of the ARES I US. During ARES I-X pre-flight modeling and analysis, six (6) independent time domain simulation tools were developed and cross validated. Each tool represents an independent implementation of a common set of models and parameters in a different simulation framework and architecture. Post flight data and reconstructed models provide the means to validate a subset of the simulations against actual flight data and to assess the accuracy of pre-flight dispersion analysis. Post flight data consists of telemetered Operational Flight Instrumentation (OFI) data primarily focused on flight computer outputs and sensor measurements as well as Best Estimated Trajectory (BET) data that estimates vehicle state information from all available measurement sources. While pre-flight models were found to provide a reasonable prediction of the vehicle flight, reconstructed models were generated to better represent and simulate the ARES I-X flight. Post flight reconstructed models include: SRB propulsion model, thrust vector bias models, mass properties, base aerodynamics, and Meteorological Estimated Trajectory (wind and atmospheric data). The result of the effort is a set of independently developed, high fidelity, time-domain simulation tools that have been cross validated and validated against flight data. This paper presents the process and results of high fidelity aerospace modeling, simulation, analysis and tool validation in the time domain.

Hough, Steven↗

Cloud and Precipitation Analyses using Merged Datasets from Two Airborne Microwave Radiometers Covering 10–183 GHz

Microwave radiometers provide valuable insight into the structure and characteristics of clouds and precipitation. In NASA’s airborne remote-sensing arsenal, the Advanced Microwave Precipitation Radiometer (AMPR) and the Conical Scanning Millimeter-wave Imaging Radiometer (CoSMIR) have been used extensively in field campaigns throughout the world. AMPR operates with four channels between 10.7 and 85.5 GHz, while CoSMIR operates with nine channels ranging from 50.3 to 183.31 GHz. Although these datasets provide key information when used separately, the combination of these radiometers covers virtually the full range of frequencies used by the Global Precipitation Measurement (GPM) Microwave Imager (GMI), enabling suborbital observations to compare with GPM spaceborne measurements. This presentation will detail the merger of AMPR and CoSMIR data during two NASA airborne field campaigns: the Integrated Precipitation and Hydrology Experiment (IPHEx) and the Olympic Mountains Experiment / Radar Definition Experiment (OLYMPEX/RADEX). Other field campaigns, such as the Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS), may be included as well. Using the merged brightness temperature dataset, features of precipitating and non-precipitating clouds containing liquid and/or ice hydrometeors will be discussed from selected flight segments. Geophysical retrievals derived from these brightness temperatures using a one-dimensional variational (1DVAR) technique and/or multi-linear regression equations will also be employed in these analyses. Observed transitions between precipitating and non-precipitating systems will be explored in greater detail. Dropsonde data will be used to provide environmental contexts throughout each flight, and additional observations (e.g., from airborne and/or land-based radar) will be incorporated to supplement the radiometer-based results. The broader implications of these results and pathways for future work will also be discussed.

Corey G. Amiot↗

Wavefront Sensing of X-Ray Telescopes

Phase Retrieval analysis of off-axis or defocused focal-plane data from telescope optics has been proven effective in understanding misalignments and optical aberrations in normal incidence telescopes. The approach is used, e.g., in commissioning of the James Webb Space Telescope (JWST) segmented primary mirror. There is a similar need for evaluating low-order figure errors of grazing incidence mirrors and nested telescope assemblies. When implemented in these systems, phase retrieval does not depend on normal incidence access to each mirror (shell) surface and, therefore, provides an effective means for evaluating nested x-ray telescopes during integration and test. We have applied a well-known phase retrieval algorithm to grazing incidence telescopes. The algorithm uses the Levenberg-Marquardt optimization procedure to perform a non-linear least-squares fit of the telescope Point Spread Function (PSF). The algorithm can also retrieve low order figure errors at visible wavelengths where optical diffraction is the dominant defect in the PSF. In this paper we will present the analytical approach and its implementation for grazing incidence mirrors of the International X-Ray Observatory (IXO). We analyze the effects of low order axial surface errors individually, and in combination on the system PSF at 633 nanometers. We demonstrate via modeling that the wavefront sensing algorithm can recover axial errors (of the grazing incidence mirrors) to a small fraction of the known axial figure errors using simulated PSFs as input data to the algorithm.

Saha, Timo↗

Proposed methods for the reconstruction of the Pioneer/Venus probe entry trajectories

In connection with the Pioneer Venus Project, an orbiter and an array of atmospheric probes are sent to Venus in separate launches during 1978. The four probes, which are transported to the planet on an instrumented Bus are released approximately 20 days prior to encounters. There are three identical small probes and a large probe. Each is designed to make in-situ atmospheric measurements. One of the principal scientific objectives of the Pioneer Venus (P/V) multiprobe mission is the determination of the structure and composition of the Venusian atmosphere down to the planet surface. To enhance the science return it is of paramount importance to be able to accurately determine the spacecraft entry trajectories. In addition to the in-situ spacecraft measurements, earth-based Doppler data are available. The process felt to be most applicable for the P/V multiprobe entry trajectory reconstruction is considered, with emphasis placed on the highly dynamic upper entry segment of the trajectory.

Findlay, J. T.↗

Spatially Characterizing Effective Timber Supply

The structure of a computer-oriented cartographic model for assessing roundwood supply for generation of base load electricity is discussed. The model provides an analytical procedure for coupling spatial information of harvesting economics and owner willingness to sell stumpages. Supply is characterized in terms of standing timber; of accessibility considering various harvesting and hauling factors; and of availability as affected by ownership and residential patterns. Factors governing accessibility to timber include effective harvesting distance to haulic roads as modified by barriers and slopes. Haul distance is expressed in units that take into account the relative ease of travel along various road types to a central processing facility. Areas of accessible timber are grouped into spatial units, termed 'timbersheds', of common access to particular haul road segments that belong to unique 'transport zones'. Timber availability considerations include size of ownership parcels, housing density and excluded areas. The analysis techniques are demonstrated for a cartographic data base in western Massachusetts.

Berry, J. K.↗

Main rotor broadband noise study in the DNW

An acoustics test of a 2/5 scale model BO-105 helicopter main rotor was conducted in the Duits-Nederlandse Windtunnel (DNW). A range of operating conditions was tested from hover to moderately high flight speeds for various climb and descent rates at different thrust settings. Diagnostic tests including rotor speed and blade geometry changes were made to better isolate and study particular broadband self noise sources. Acoustic data in the form of acoustic pressure time histories and power spectra are used to demonstrate the regions of importance of the different broadband noise sources and their sensitivity to operating conditions. To help interpret the data, comparisons are made to predictions of rotor broadband noise. The predictions are based on self noise data previously obtained from isolated airfoil sections and the use of the NASA ROTONET program to define rotor performance and to sum contributions of noise from individual blade segments. An important result herein is the identification and articulation of a previously unheralded rotor broadband noise source. This source is blade-turbulent wake interaction (BWI) noise which dominates the spectra in the mid-frequencies for off-peak blade-vortex interaction (BVI) noise flight conditions.

Brooks, Thomas F.↗

An Imaging Spectrometric Observatory (ISO)

The Imaging Spectrometric Observatory (ISO) is designed for low light level spectroscopy of both the day and night side of the earth. The instrument is composed of five spectrometers, each of which covers part of the total wavelength range of 30 to 1300 nm spanned by the instrument. Wavelength resolution varies between 0.2 and 0.6 nm over the spectral range. The five spectrometers are each optimized for a portion of the spectrum by the choice of mirror reflective coatings and detector photocathode materials. The full spectral range for each spectrometer is covered in a total of 11 grating steps. The Imaging Spectrometric Observatory was flown for the first time on the Spacelab 1 mission during which it acquired almost 40 hr of observations. The ISO investigation to be flown on the Atmospheric Laboratory for Applications and Science (ATLAS 1) mission will draw on the experience gained from the data gathered on Spacelab 1. The detector system in each spectrometer was upgraded to provide both higher sensitivity at low light levels and simultaneous imaging over larger spectral segments than was achieved on Spacelab 1. In addition, the instrument and the observing sequences were modified to allow observation of the sun in the extreme ultraviolet. A summary of ISO parameters for ATLAS 1 (scheduled for late 1990) is given.

Torr, M. R.↗

An adaptive technique to maximize lossless image data compression of satellite images

Data compression will pay an increasingly important role in the storage and transmission of image data within NASA science programs as the Earth Observing System comes into operation. It is important that the science data be preserved at the fidelity the instrument and the satellite communication systems were designed to produce. Lossless compression must therefore be applied, at least, to archive the processed instrument data. In this paper, we present an analysis of the performance of lossless compression techniques and develop an adaptive approach which applied image remapping, feature-based image segmentation to determine regions of similar entropy and high-order arithmetic coding to obtain significant improvements over the use of conventional compression techniques alone. Image remapping is used to transform the original image into a lower entropy state. Several techniques were tested on satellite images including differential pulse code modulation, bi-linear interpolation, and block-based linear predictive coding. The results of these experiments are discussed and trade-offs between computation requirements and entropy reductions are used to identify the optimum approach for a variety of satellite images. Further entropy reduction can be achieved by segmenting the image based on local entropy properties then applying a coding technique which maximizes compression for the region. Experimental results are presented showing the effect of different coding techniques for regions of different entropy. A rule-base is developed through which the technique giving the best compression is selected. The paper concludes that maximum compression can be achieved cost effectively and at acceptable performance rates with a combination of techniques which are selected based on image contextual information.

Stewart, Robert J.↗

Polar Spin Axis Anomaly

The Global Geospace Science (GGS) Polar Plasma Laboratory (POLAR) spacecraft was launched on February 24, 1996, by a Delta 2. The spacecraft, a major axis spinner, appeared to function nominally throughout the early mission phase, which included several deployments, and orbit and attitude maneuvers. Of particular interest is the fact that the spacecraft was launched with a deliberate dynamic imbalance. During a segment of early orbit operations, a pair of Lanyard Deployed Booms (LDB) were extended. These booms were not identical; the intent was that the spacecraft would be nearly dynamically balanced after they were deployed. The spacecraft contained two dynamic balance mechanisms intended to fine tune the balance on orbit. However, subsequent images taken by the science instruments on the Despun Platform during the dynamic balancing segment indicated an offset of the principal spin axis from the geometric axis. This offset produced a sinusoidal blurring of the science images sufficiently large to degrade science data below mission requirement specifications. In the end, the imbalance encountered in flight was significantly outside the correction capability of the balances. The purpose of this paper is to examine the flight data during the various deployment and maneuver stages of the early orbit operations coupled with analytical simulations to discuss some of the potential causes of the resultant imbalance.

Crouse, Patrick L.↗

Improving Acoustic Models by Watching Television

Obtaining sufficient labelled training data is a persistent difficulty for speech recognition research. Although well transcribed data is expensive to produce, there is a constant stream of challenging speech data and poor transcription broadcast as closed-captioned television. We describe a reliable unsupervised method for identifying accurately transcribed sections of these broadcasts, and show how these segments can be used to train a recognition system. Starting from acoustic models trained on the Wall Street Journal database, a single iteration of our training method reduced the word error rate on an independent broadcast television news test set from 62.2% to 59.5%.

Witbrock, Michael J.↗

MOLA-derived Roughness Data Used to Predict Surface Scattering for Mars Subsurface Radar Sounding

The Mars Express orbiter, to be launched by the European Space Agency in 2003, will carry a low-frequency radar sounding instrument, MARSIS (Mars Advanced Radar for Subsurface and Ionospheric Sounding). The primary goal of MARSIS is to map the distribution of water, both solid and (if present) liquid, in the upper several km of the martian crust. Detecting discontinuities in the crust, such as an ice-water transition, presents many challenges for a Mars orbital radar sounder. One challenge that must be overcome is the presence of radar scattering (echoes) from the surface of Mars, expected to be detected by the sounder antennas at the same time as any echoes aris-ing from subsurface interfaces. As the transmitted spherical wavefront spreads within the crust of Mars,. it also interacts with surface topography at off-nadir positions, creating a "clutter" signal that can mask the subsurface echoes. The MARSIS instrument will utilize Doppler filtering to limit the off-nadir clutter in the along-track direction, and a nadir-null secondary antenna to identify strong off-nadir clutter from the cross-track direction. To evaluate the effects of off-nadir surface clutter and the capability of these schemes to reduce the clutter, it is necessary to predict the range of scattering behavior that may be expected from martian surface topography. In this paper, we utilize Mars Orbital Laser Altimeter (MOLA) data from the current Mars Global Surveyor mission to characterize the topographic roughness of a variety of martian terrain types, at scales relevant to the MARSIS clutter problem. Segments of MOLA altimetry profiles are reduced to the topographic parameters rms slope and fractal dimension, which then are used as inputs to a near-nadir radar scattering model to predict the strength of the clutter signal. Additional information is contained in the original extended abstract.

Plaut, J. J.↗

Spatial Analysis of Great Lakes Regional Icing Cloud Liquid Water Content

Abstract Clustering of cloud microphysical conditions, such as liquid water content (LWC) and drop size, can affect the rate and shape of ice accretion and the airworthiness of aircraft. Clustering may also degrade the accuracy of cloud LWC measurements from radars and microwave radiometers being developed by the government for remotely mapping icing conditions ahead of aircraft in flight. This paper evaluates spatial clustering of LWC in icing clouds using measurements collected during NASA research flights in the Great Lakes region. We used graphical and analytical approaches to describe clustering. The analytical approach involves determining the average size of clusters and computing a clustering intensity parameter. We analyzed flight data composed of 1-s-frequency LWC measurements for 12 periods ranging from 17.4 minutes (73 km) to 45.3 minutes (190 km) in duration. Graphically some flight segments showed evidence of consistency with regard to clustering patterns. Cluster intensity varied from 0.06, indicating little clustering, to a high of 2.42. Cluster lengths ranged from 0.1 minutes (0.6 km) to 4.1 minutes (17.3 km). Additional analyses will allow us to determine if clustering climatologies can be developed to characterize cluster conditions by region, time period, or weather condition. Introduction

Ryerson, Charles C.↗