Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “volume estimation”

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

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

At least 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

Bistatic scatter from rain

An experimental investigation of bistatic scatter from rain was conducted, using a 143-km scatter path at frequencies of 4.5 and 7.7 GHz. The ratio of transmitted to received power (transmission loss) was measured for scattering angles ranging from 6 to 130 deg. Simultaneous weather radar observations were made at a frequency of 1.3 GHz. Transmission loss estimates for the bistatic scatter path were computed, using the weather radar data, the bistatic radar equation, and a model for the scattering cross section per unit volume of rain based upon Rayleigh scattering by an ensemble of water spheres. The measured and estimated transmission loss values were compared to test the use of the scattering model for the estimation of interference. The averaged ratio of measured-to-calculated transmission loss for the 4.5 GHz data is 1.2 plus or minus 0.4 dB. The averaged ratio for the 7.7 GHz data is -1.6 plus or minus 0.5 dB.

Crane, R. K.↗
NASA NTRS1991

Measures of large-scale structure in the CfA redshift survey slices

Variations of the counts-in-cells with cell size are used here to define two statistical measures of large-scale clustering in three 6 deg slices of the CfA redshift survey. A percolation criterion is used to estimate the filling factor which measures the fraction of the total volume in the survey occupied by the large-scale structures. For the full 18 deg slice of the CfA redshift survey, f is about 0.25 + or - 0.05. After removing groups with more than five members from two of the slices, variations of the counts in occupied cells with cell size have a power-law behavior with a slope beta about 2.2 on scales from 1-10/h Mpc. Application of both this statistic and the percolation analysis to simulations suggests that a network of two-dimensional structures is a better description of the geometry of the clustering in the CfA slices than a network of one-dimensional structures. Counts-in-cells are also used to estimate at 0.3 galaxy h-squared/Mpc the average galaxy surface density in sheets like the Great Wall.

De Lapparent, Valerie↗
NASA NTRS1999

Coherent Doppler Lidar Data Products from Space-Based Platforms

Coherent Doppler lidar is a promising technique for the global measurements of winds using a space-based platform. Doppler lidar produces estimates of the radial component of the velocity vector averaged over the resolution volume of the measurement. Profiles of the horizontal vector winds are produced by scanning the lidar beam or stepping the lidar beam through a sequence of different angles (step-stare). The first design for space-based measurements proposed a conical scan which requires a high power laser to produce acceptable signal levels for every laser pulse. Performance is improved by fixing the laser beam and accumulating the signal from many lidar pulses for each range-gate. This also improves the spatial averaging of the wind estimates and reduces the threshold signal energy required for a good estimate. Coherent Doppler lidar performance for space-based operation is determined using computer simulations and including the wind variability over the measurement volume as well as the variations of the atmospheric aerosol backscatter.

Frehlich, Rod↗
NASA NTRS2011

Method to Estimate the Dissolved Air Content in Hydraulic Fluid

In order to verify the air content in hydraulic fluid, an instrument was needed to measure the dissolved air content before the fluid was loaded into the system. The instrument also needed to measure the dissolved air content in situ and in real time during the de-aeration process. The current methods used to measure the dissolved air content require the fluid to be drawn from the hydraulic system, and additional offline laboratory processing time is involved. During laboratory processing, there is a potential for contamination to occur, especially when subsaturated fluid is to be analyzed. A new method measures the amount of dissolved air in hydraulic fluid through the use of a dissolved oxygen meter. The device measures the dissolved air content through an in situ, real-time process that requires no additional offline laboratory processing time. The method utilizes an instrument that measures the partial pressure of oxygen in the hydraulic fluid. By using a standardized calculation procedure that relates the oxygen partial pressure to the volume of dissolved air in solution, the dissolved air content is estimated. The technique employs luminescent quenching technology to determine the partial pressure of oxygen in the hydraulic fluid. An estimated Henry s law coefficient for oxygen and nitrogen in hydraulic fluid is calculated using a standard method to estimate the solubility of gases in lubricants. The amount of dissolved oxygen in the hydraulic fluid is estimated using the Henry s solubility coefficient and the measured partial pressure of oxygen in solution. The amount of dissolved nitrogen that is in solution is estimated by assuming that the ratio of dissolved nitrogen to dissolved oxygen is equal to the ratio of the gas solubility of nitrogen to oxygen at atmospheric pressure and temperature. The technique was performed at atmospheric pressure and room temperature. The technique could be theoretically carried out at higher pressures and elevated temperatures.

Hauser, Daniel M.↗
DOE OSTI2014

A multifractal approach to space-filling recovery for PET quantification

Purpose: A new image-based methodology is developed for estimating the apparent space-filling properties of an object of interest in PET imaging without need for a robust segmentation step and used to recover accurate estimates of total lesion activity (TLA). Methods: A multifractal approach and the fractal dimension are proposed to recover the apparent space-filling index of a lesion (tumor volume, TV) embedded in nonzero background. A practical implementation is proposed, and the index is subsequently used with mean standardized uptake value (SUV {sub mean}) to correct TLA estimates obtained from approximate lesion contours. The methodology is illustrated on fractal and synthetic objects contaminated by partial volume effects (PVEs), validated on realistic {sup 18}F-fluorodeoxyglucose PET simulations and tested for its robustness using a clinical {sup 18}F-fluorothymidine PET test–retest dataset. Results: TLA estimates were stable for a range of resolutions typical in PET oncology (4–6 mm). By contrast, the space-filling index and intensity estimates were resolution dependent. TLA was generally recovered within 15% of ground truth on postfiltered PET images affected by PVEs. Volumes were recovered within 15% variability in the repeatability study. Results indicated that TLA is a more robust index than other traditional metrics such as SUV {sub mean} or TV measurements across imaging protocols. Conclusions: The fractal procedure reported here is proposed as a simple and effective computational alternative to existing methodologies which require the incorporation of image preprocessing steps (i.e., partial volume correction and automatic segmentation) prior to quantification.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-J-135: An Investigation of Ultrasound Imaging for 3D Intra-Fraction Prostate Motion Estimation

Purpose: This study investigates the use of a mechanically swept 3D ultrasound (US) probe to estimate intra-fraction motion of the prostate during radiation therapy using an US phantom and simulated transperineal imaging. Methods: A 3D motion platform was used to translate an US speckle phantom while simulating transperineal US imaging. Motion patterns for five representative types of prostate motion, generated from patient data previously acquired with a Calypso system, were using to move the phantom in 3D. The phantom was also implanted with fiducial markers and subsequently tracked using the CyberKnife kV x-ray system for comparison. A normalised cross correlation block matching algorithm was used to track speckle patterns in 3D and 2D US data. Motion estimation results were compared with known phantom translations. Results: Transperineal 3D US could track superior-inferior (axial) and anterior-posterior (lateral) motion to better than 0.8 mm root-mean-square error (RMSE) at a volume rate of 1.7 Hz (comparable with kV x-ray tracking RMSE). Motion estimation accuracy was poorest along the US probe's swept axis (right-left; RL; RMSE < 4.2 mm) but simple regularisation methods could be used to improve RMSE (< 2 mm). 2D US was found to be feasible for slowly varying motion (RMSE < 0.5 mm). 3D US could also allow accurate radiation beam gating with displacement thresholds of 2 mm and 5 mm exhibiting a RMSE of less than 0.5 mm. Conclusion: 2D and 3D US speckle tracking is feasible for prostate motion estimation during radiation delivery. Since RL prostate motion is small in magnitude and frequency, 2D or a hybrid (2D/3D) US imaging approach which also accounts for potential prostate rotations could be used. Regularisation methods could be used to ensure the accuracy of tracking data, making US a feasible approach for gating or tracking in standard or hypo-fractionated prostate treatments.

60 APPLIED LIFE SCIENCES↗
NASA NTRSDate not supplied

AERACEPT (Aerosol Rapid Analysis Combined Entry Probe/sonde Technology): Enabling Technology for Missions to the Venus Clouds

AERACEPT (AErosol Rapid Analysis Combined Entry Probe/sonde Technology) is an early-stage technology allowing a single aeroshell body to act as both an entry vehicle and aerosol-sampling passive descent sonde. AERACEPT does not require heat shield separation, deployable parachutes, or descent control, thus reducing the mass, volume, and complexity of planetary aerosol sampling. AERACEPT is particularly well suited for a Venus mission, where the particles of greatest interest are within the subsonic descent regime. AERACEPT uses the aeroshell’s own velocity to drive aerosol capture and separation through a series of embedded inlets. It takes advantage of recently developed thermal protection materials (3D-CC and HEEET) in combination with heritage aerosol sampling technologies from both planetary and airborne science (high-speed inlets and particle separation). The trade space for a given descent trajectory includes the particle capture efficiency for a given size, the bias introduced in the sampled particle size and concentration distributions, and the thermal alteration experienced by the particles during their brief exposure to the internal flow environment. AERACEPT is included in the Nephele mission concept study for a small spacecraft targeting the Venus middle and lower cloud layers. Nephele complements larger missions targeting Venus atmospheric gas analysis, such as DAVINCI and Venera-D, by specifically targeting cloud and haze particles. Because of the short lifetime of the probe in the lower atmosphere, Nephele requires a fast cadence of analysis of the captured particles, and includes the VOLTR dual optical spectrometer (SERS/LIBS) as part of its notional payload. Preliminary modeling based on the Nephele trajectory at 63 km to 39 km indicates AERACEPT can limit sample heating to 30-60 K above ambient. A modified particle tracking model has been implemented to estimate capture efficiency of particles larger than 0.1 µm and total sample volume as part of an inlet and interal flow path geometry trade study. Further modeling and empirical testing is underway to improve these estimates.

AERACEPT↗
DOE OSTI2017

Modeling the Effect of Glass Microballoon (GMB) Volume Fraction on Behavior of Sylgard/GMB Composites.

This work was done to support customer questions about whether a Sylgard/Glass Microballoon (GMB) potting material in current use could be replaced with pure Sylgard and if this would significantly change stresses imparted to internal components under thermal cycling conditions. To address these questions, we provide micromechanics analysis of Sylgard/GMB materials using both analytic composite theory and finite element simulations to better understand the role of the GMB volume fraction in determining thermal expansion coefficient, elastic constants, and behavior in both confined and unconfined compression boundary value problems. A key finding is that damage accumulation in the material from breakage of GMBs significantly limits the global stress magnitude and results in a plateau stress behavior over large ranges of compressive strain. The magnitude of this plateau stress is reduced with higher volume fractions of GMBs. This effect is particularly pronounced in confined compression, which we estimate bears the most similarity to the application of interest. This stress-limiting damage mechanism is not present in pure Sylgard, however, and the result is much higher stresses under confined compression. Thus, we recommend that some volume fraction greater than 10% GMBs be used for confined deformation applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Using Sentinel-1 and GRACE satellite data to monitor the hydrological variations within the Tulare Basin, California

Abstract Subsidence induced by groundwater depletion is a grave problem in many regions around the world, leading to a permanent loss of groundwater storage within an aquifer and even producing structural damage at the Earth’s surface. California’s Tulare Basin is no exception, experiencing about a meter of subsidence between 2015 and 2020. However, understanding the relationship between changes in groundwater volumes and ground deformation has proven difficult. We employ surface displacement measurements from Interferometric Synthetic Aperture Radar (InSAR) and gravimetric estimates of terrestrial water storage from the Gravity Recovery and Climate Experiment (GRACE) satellite pair to characterize the hydrological dynamics within the Tulare basin. The removal of the long-term aquifer compaction from the InSAR time series reveals coherent short-term variations that correlate with hydrological features. For example, in the winter of 2018–2019 uplift is observed at the confluence of several rivers and streams that drain into the southeastern edge of the basin. These observations, combined with estimates of mass changes obtained from the orbiting GRACE satellites, form the basis for imaging the monthly spatial variations in water volumes. This approach facilitates the quick and effective synthesis of InSAR and gravimetric datasets and will aid efforts to improve our understanding and management of groundwater resources around the world.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2022

Using Sentinel-1 and GRACE Satellite Data to Monitor the Hydrological Variations Within the Tulare Basin, California

Subsidence induced by groundwater depletion is a grave problem in many regions around the world, leading to a permanent loss of groundwater storage within an aquifer and even producing structural damage at the Earth’s surface. California’s Tulare Basin is no exception, experiencing about a meter of subsidence between 2015 and 2020. However, understanding the relationship between changes in groundwater volumes and ground deformation has proven difficult. We employ surface displacement measurements from Interferometric Synthetic Aperture Radar (InSAR) and gravimetric estimates of terrestrial water storage from the Gravity Recovery and Climate Experiment (GRACE) satellite pair to characterize the hydrological dynamics within the Tulare basin. The removal of the long-term aquifer compaction from the InSAR time series reveals coherent short-term variations that correlate with hydrological features. For example, in the winter of 2018–2019 uplift is observed at the confluence of several rivers and streams that drain into the southeastern edge of the basin. These observations, combined with estimates of mass changes obtained from the orbiting GRACE satellites, form the basis for imaging the monthly spatial variations in water volumes. This approach facilitates the quick and effective synthesis of InSAR and gravimetric datasets and will aid efforts to improve our understanding and management of groundwater resources around the world.

Donald W. Vasco↗
DOE OSTI2014

A Prospective Study of Ventilation Rates and Illness Absence in California Office Buildings

Background – This study investigated the associations of ventilation rates (VRs), estimated from indoor CO2 concentrations, in offices with the amount of respiratory infections, illness absences, and building-related health symptoms in occupants. Methods – Office buildings were recruited from three California climate zones. In one or more study spaces within each building, real-time logging sensors measured carbon dioxide, temperature, and relative humidity for one year. Ventilation rates were estimated using daily peak CO2 levels, and also using an alternative metric. Data on occupants and health outcomes were collected through web-based surveys every three months. Multivariate models were used to assess relationships between metrics of ventilation rate or CO2 and occupant outcomes. For all outcomes, negative associations were hypothesized with VR metrics, and positive associations with CO2 metrics. Results – Difficulty recruiting buildings and low survey response limited sample size and study power. In 16 studied spaces within 9 office buildings, VRs were uniformly high over the year, from twice to over nine times the California office VR standard (7 L/s or 15 cfm per person). VR and CO2 metrics had no statistically significant relationships with occupant outcomes, except for a small significantly positive association of the alternative VR metric with respiratory illness-related absence, contrary to hypotheses. Conclusions– The very high time-averaged VRs in the California office buildings studied presumably resulted from “economizer cycles” bringing in large volumes of outdoor air; however, in almost all buildings even the estimated minimum VRs supplied (without the economizer) substantially exceeded the minimum required VR. These high VRs may explain the absence of hypothesized relationships with occupant outcomes. Among uniformly high VRs, little variation in contaminant concentration and occupant effects would be expected. These findings may provide initial evidence for an upper bound of the range of VRs within which increased VRs provide benefits in reducing illness absence.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2005

Development of Parametric Mass and Volume Models for an Aerospace SOFC/Gas Turbine Hybrid System

In aerospace power systems, mass and volume are key considerations to produce a viable design. The utilization of fuel cells is being studied for a commercial aircraft electrical power unit. Based on preliminary analyses, a SOFC/gas turbine system may be a potential solution. This paper describes the parametric mass and volume models that are used to assess an aerospace hybrid system design. The design tool utilizes input from the thermodynamic system model and produces component sizing, performance, and mass estimates. The software is designed such that the thermodynamic model is linked to the mass and volume model to provide immediate feedback during the design process. It allows for automating an optimization process that accounts for mass and volume in its figure of merit. Each component in the system is modeled with a combination of theoretical and empirical approaches. A description of the assumptions and design analyses is presented.

Tornabene, Robert↗
DOE OSTI2015

SU-E-T-594: Out-Of-Field Neutron and Gamma Dose Estimated Using TLD-600/700 Pairs in the Wobbling Proton Therapy System

Purpose: Secondary fast neutrons and gamma rays are mainly produced due to the interaction of the primary proton beam with the beam delivery nozzle. These secondary radiation dose to patients and radiation workers are unwanted. The purpose of this study is to estimate the neutron and gamma dose equivalent out of the treatment volume during the wobbling proton therapy system. Methods: Two types of thermoluminescent (TL) dosimeters, TLD-600 ({sup 6}LiF: Mg, Ti) and TLD-700 ({sup 7}LiF: Mg, Ti) were used in this study. They were calibrated in the standard neutron and gamma sources at National Standards Laboratory. Annealing procedure is 400°C for 1 hour, 100°C for 2 hours and spontaneously cooling down to the room temperature in a programmable oven. Two-peak method (a kind of glow curve analysis technique) was used to evaluate the TL response corresponding to the neutron and gamma dose. The TLD pairs were placed outside the treatment field at the neutron-gamma mixed field with 190-MeV proton beam produced by the wobbling system through the polyethylene plate phantom. The results of TLD measurement were compared to the Monte Carlo simulation. Results: The initial experiment results of calculated dose equivalents are 0.63, 0.38, 0.21 and 0.13 mSv per Gy outside the field at the distance of 50, 100, 150 and 200 cm. Conclusion: The TLD-600 and TLD-700 pairs are convenient to estimate neutron and gamma dosimetry during proton therapy. However, an accurate and suitable glow curve analysis technique is necessary. During the wobbling system proton therapy, our results showed that the neutron and gamma doses outside the treatment field are noticeable. This study was supported by the grants from the Chang Gung Memorial Hospital (CMRPD1C0682)

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRS2008

Reducing Spaceborne-Doppler-Radar Rainfall-Velocity Error

A combined frequency-time (CFT) spectral moment estimation technique has been devised for calculating rainfall velocity from measurement data acquired by a nadir-looking spaceborne Doppler weather radar system. Prior spectral moment estimation techniques used for this purpose are based partly on the assumption that the radar resolution volume is uniformly filled with rainfall. The assumption is unrealistic in general but introduces negligible error in application to airborne radar systems. However, for spaceborne systems, the combination of this assumption and inhomogeneities in rainfall [denoted non-uniform beam filling (NUBF)] can result in velocity measurement errors of several meters per second. The present CFT spectral moment estimation technique includes coherent processing of a series of Doppler spectra generated in a standard manner from data over measurement volumes that are partially overlapping in the along-track direction. Performance simulation of this technique using high-resolution data from an airborne rain-mapping radar shows that a spaceborne Ku-band Doppler radar operating at signal-to-noise ratios greater than 10 dB can achieve root-mean-square accuracy between 0.5 and 0.6 m/s in vertical-velocity estimates.

Tanelli, Simone↗
DOE OSTI2011

Variable Circular Collimator in Robotic Radiosurgery: A Time-Efficient Alternative to a Mini-Multileaf Collimator?

Purpose: Compared with many small circular beams used in CyberKnife treatments, beam's eye view-shaped fields are generally more time-efficient for dose delivery. However, beam's eye view-shaping devices, such as a mini-multileaf collimator (mMLC), are not presently available for CyberKnife, although a variable-aperture collimator (Iris, 12 field diameters; 5-60 mm) is available. We investigated whether the Iris can mimic noncoplanar mMLC treatments using a limited set of principal beam orientations (nodes) to produce time-efficient treatment plans. Methods and Materials: The data from 10 lung cancer patients and the beam-orientation optimization algorithm 'Cycle' were used to generate stereotactic treatment plans (3 x 20 Gy) for a CyberKnife virtually equipped with a mMLC. Typically, 10-16 favorable beam orientations were selected from 117 available robot node positions using beam's eye view-shaped fields with uniform fluence. Second, intensity-modulated Iris plans were generated by inverse optimization of nonisocentric circular candidate beams targeted from the same nodes selected in the mMLC plans. The plans were evaluated using the mean lung dose, lung volume receiving {>=}20 Gy, conformality index, number of nodes, beams, and monitor units, and estimated treatment time. Results: The mMLC plans contained an average of 12 nodes and 11,690 monitor units. For a comparable mean lung dose, the Iris plans contained 12 nodes, 64 beams, and 21,990 monitor units. The estimated fraction duration was 12.2 min (range, 10.8-13.5) for the mMLC plans and 18.4 min (range, 12.9-28.5) for the Iris plans. In contrast to the mMLC plans, the treatment time for the Iris plans increased with an increasing target volume. The Iris plans were, on average, 40% longer than the corresponding mMLC plans for small targets (<80 cm{sup 3}) and {<=}121% longer for larger targets. For a comparable conformality index, similar results were obtained. Conclusion: For stereotactic lung irradiation, time-efficient and high-quality plans were obtained for robotic-controlled noncoplanar treatments using a mMLC. Iris is a time-efficient alternative for small targets, with similar or better plan quality.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

A technique for estimating 4D-CBCT using prior knowledge and limited-angle projections

Purpose: To develop a technique to estimate onboard 4D-CBCT using prior information and limited-angle projections for potential 4D target verification of lung radiotherapy.Methods: Each phase of onboard 4D-CBCT is considered as a deformation from one selected phase (prior volume) of the planning 4D-CT. The deformation field maps (DFMs) are solved using a motion modeling and free-form deformation (MM-FD) technique. In the MM-FD technique, the DFMs are estimated using a motion model which is extracted from planning 4D-CT based on principal component analysis (PCA). The motion model parameters are optimized by matching the digitally reconstructed radiographs of the deformed volumes to the limited-angle onboard projections (data fidelity constraint). Afterward, the estimated DFMs are fine-tuned using a FD model based on data fidelity constraint and deformation energy minimization. The 4D digital extended-cardiac-torso phantom was used to evaluate the MM-FD technique. A lung patient with a 30 mm diameter lesion was simulated with various anatomical and respirational changes from planning 4D-CT to onboard volume, including changes of respiration amplitude, lesion size and lesion average-position, and phase shift between lesion and body respiratory cycle. The lesions were contoured in both the estimated and “ground-truth” onboard 4D-CBCT for comparison. 3D volume percentage-difference (VPD) and center-of-mass shift (COMS) were calculated to evaluate the estimation accuracy of three techniques: MM-FD, MM-only, and FD-only. Different onboard projection acquisition scenarios and projection noise levels were simulated to investigate their effects on the estimation accuracy.Results: For all simulated patient and projection acquisition scenarios, the mean VPD (±S.D.)/COMS (±S.D.) between lesions in prior images and “ground-truth” onboard images were 136.11% (±42.76%)/15.5 mm (±3.9 mm). Using orthogonal-view 15°-each scan angle, the mean VPD/COMS between the lesion in estimated and “ground-truth” onboard images for MM-only, FD-only, and MM-FD techniques were 60.10% (±27.17%)/4.9 mm (±3.0 mm), 96.07% (±31.48%)/12.1 mm (±3.9 mm) and 11.45% (±9.37%)/1.3 mm (±1.3 mm), respectively. For orthogonal-view 30°-each scan angle, the corresponding results were 59.16% (±26.66%)/4.9 mm (±3.0 mm), 75.98% (±27.21%)/9.9 mm (±4.0 mm), and 5.22% (±2.12%)/0.5 mm (±0.4 mm). For single-view scan angles of 3°, 30°, and 60°, the results for MM-FD technique were 32.77% (±17.87%)/3.2 mm (±2.2 mm), 24.57% (±18.18%)/2.9 mm (±2.0 mm), and 10.48% (±9.50%)/1.1 mm (±1.3 mm), respectively. For projection angular-sampling-intervals of 0.6°, 1.2°, and 2.5° with the orthogonal-view 30°-each scan angle, the MM-FD technique generated similar VPD (maximum deviation 2.91%) and COMS (maximum deviation 0.6 mm), while sparser sampling yielded larger VPD/COMS. With equal number of projections, the estimation results using scattered 360° scan angle were slightly better than those using orthogonal-view 30°-each scan angle. The estimation accuracy of MM-FD technique declined as noise level increased.Conclusions: The MM-FD technique substantially improves the estimation accuracy for onboard 4D-CBCT using prior planning 4D-CT and limited-angle projections, compared to the MM-only and FD-only techniques. It can potentially be used for the inter/intrafractional 4D-localization verification.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS1971

NONAME system description, volume 1

Description of NONAME program and system for determining and analyzing satellite orbits and estimating geodetic parameters

Williamson, R. G.↗