On the evaluation of the gravity effects of polyhedral bodies and a consistent treatment of related singularities

We show that the singularities which can affect the computation of the gravity effects (potential, gravity and tensor gradient fields) can be systematically addressed by invoking distribution theory and suitable formulas of differential calculus. Thus, differently from previous contributions on the subject, the use of a-posteriori corrections of the formulas derived in absence of singularities can be ruled out. The general approach presented in the paper is further specialized to the case of polyhedral bodies and detailed for a rectangular prism having a constant mass density. With reference to this last case, we derive novel expressions for the related gravitational field, as well as for its first and second derivative, at an observation point coincident with a prism vertex and show that they turn out to be more compact than the ones reported in the specialized literature.     

Crop Residue Discrimination Using Ground-Based Hyperspectral Data

Crop residue has become an increasingly important factor in agriculture management. It assists in the reduction of soil erosion and is an important source of soil organic carbon (soil carbon sequestration). In recent past, remote sensing, especially narrowband, data have been explored for crop residue assessment. In this context, a study was carried out to identify different narrow-bands and evaluate the performance of SWIR region based spectral indices for crop residue discrimination. Ground based hyperspectral data collected for wheat crop residue was analyzed using Stepwise Discriminant Analysis (SDA) technique to select significant bands for discrimination. Out of the seven best bands selected to discriminate between matured crop, straw heap, combine-harvested field with stubbles and soil, four bands were from SWIR (1980, 2030, 2200, 2440 nm) region. Six spectral indices were computed, namely CAI, LCA, SINDRI, NDSVI, NDI5 and hSINDRI for crop residue discrimination. LCA and CAI showed to be best (F?>?115) in discriminating above classes, while LCA and SINDRI were best (F?>?100) among all indices in discriminating crop residue under different harvesting methods. Comparison of different spectral resolution (from 1 nm to 150 nm) showed that for crop residue discrimination a resolution of 100 nm at 2100–2300 m region would be sufficient to discriminate crop residue from other co-existing classes.     

A Composite Semisupervised SVM for Classification of Hyperspectral Images

This letter presents a novel composite semisupervised support vector machine (SVM) for the spectral-spatial classification of hyperspectral images. In particular, the proposed technique exploits the following: 1) unlabeled data for increasing the reliability of the training phase when few training samples are available and 2) composite kernel functions for simultaneously taking into account spectral and spatial information included in the considered image. Experiments carried out on a hyperspectral image pointed out the effectiveness of the presented technique, which resulted in a significant increase of the classification accuracy with respect to both supervised SVMs and progressive semisupervised SVMs with single kernels, as well as supervised SVMs with composite kernels.     

The determination of potential difference by the joint application of measured and synthetical gravity data: a case study in Hungary

In an elementary approach every geometrical height difference between the staff points of a levelling line should have a corresponding average g value for the determination of potential difference in the Earth’s gravity field. In practice this condition requires as many gravity data as the number of staff points if linear variation of g is assumed between them. Because of the expensive fieldwork, the necessary data should be supplied from different sources. This study proposes an alternative solution, which is proved at a test bed located in the Mecsek Mountains, Southwest Hungary, where a detailed gravity survey, as dense as the staff point density (~1 point/34 m), is available along a 4.3-km-long levelling line. In the first part of the paper the effect of point density of gravity data on the accuracy of potential difference is investigated. The average g value is simply derived from two neighbouring g measurements along the levelling line, which are incrementally decimated in the consecutive turns of processing. The results show that the error of the potential difference between the endpoints of the line exceeds 0.1 mm in terms of length unit if the sampling distance is greater than 2 km. Thereafter, a suitable method for the densification of the decimated g measurements is provided. It is based on forward gravity modelling utilising a high-resolution digital terrain model, the normal gravity and the complete Bouguer anomalies. The test shows that the error is only in the order of 10⁻³mm even if the sampling distance of g measurements is 4 km. As a component of the error sources of levelling, the ambiguity of the levelled height difference which is the Euclidean distance between the inclined equipotential surfaces is also investigated. Although its effect accumulated along the test line is almost zero, it reaches 0.15 mm in a 1-km-long intermediate section of the line.     

On the use of iterative re-weighting least-squares and outlier detection for empirically modelling rates of vertical displacement

The proper identification and removal of outliers in the combination of rates of vertical displacements derived from GPS, tide gauges/satellite altimetry, and GRACE observations is presented. Outlier detection is a necessary pre-screening procedure in order to ensure reliable estimates of stochastic properties of the observations in the combined least-squares adjustment (via rescaling of covariance matrices) and to ensure that the final vertical motion model is not corrupted and/or distorted by erroneous data. Results from this study indicate that typical data snooping methods are inadequate in dealing with these heterogeneous data sets and their stochastic properties. Using simulated vertical displacement rates, it is demonstrated that a large variety of outliers (random scattered and adjacent, as well as jointly influential) can be dealt with if an iterative re-weighting least-squares adjustment is combined with a robust median estimator. Moreover, robust estimators are efficient in areas weakly constrained by the data, where even high quality observations may appear to be erroneous if their estimates are largely influenced by outliers. Four combined models for the vertical motion in the region of the Great Lakes are presented. The computed vertical displacements vary between  − 2 mm/year (subsidence) along the southern shores and 3 mm/year (uplift) along the northern shores. The derived models provide reliable empirical constraints and error bounds for postglacial rebound models in the region.     

Phase center modeling for LEO GPS receiver antennas and its impact on precise orbit determination

Most satellites in a low-Earth orbit (LEO) with demanding requirements on precise orbit determination (POD) are equipped with on-board receivers to collect the observations from Global Navigation Satellite systems (GNSS), such as the Global Positioning System (GPS). Limiting factors for LEO POD are nowadays mainly encountered with the modeling of the carrier phase observations, where a precise knowledge of the phase center location of the GNSS antennas is a prerequisite for high-precision orbit analyses. Since 5 November 2006 (GPS week 1400), absolute instead of relative values for the phase center location of GNSS receiver and transmitter antennas are adopted in the processing standards of the International GNSS Service (IGS). The absolute phase center modeling is based on robot calibrations for a number of terrestrial receiver antennas, whereas compatible antenna models were subsequently derived for the remaining terrestrial receiver antennas by conversion (from relative corrections), and for the GNSS transmitter antennas by estimation. However, consistent receiver antenna models for space missions such as GRACE and TerraSAR-X, which are equipped with non-geodetic receiver antennas, are only available since a short time from robot calibrations. We use GPS data of the aforementioned LEOs of the year 2007 together with the absolute antenna modeling to assess the presently achieved accuracy from state-of-the-art reduced-dynamic LEO POD strategies for absolute and relative navigation. Near-field multipath and cross-talk with active GPS occultation antennas turn out to be important and significant sources for systematic carrier phase measurement errors that are encountered in the actual spacecraft environments. We assess different methodologies for the in-flight determination of empirical phase pattern corrections for LEO receiver antennas and discuss their impact on POD. By means of independent K-band measurements, we show that zero-difference GRACE orbits can be significantly improved from about 10 to 6 mm K-band standard deviation when taking empirical phase corrections into account, and assess the impact of the corrections on precise baseline estimates and further applications such as gravity field recovery from kinematic LEO positions.     

Changes in timber haul emissions in the context of shifting forest management and infrastructure

Background  

Although significant amounts of carbon may be stored in harvested wood products, the extraction of that carbon from the forest generally entails combustion of fossil fuels. The transport of timber from the forest to primary milling facilities may in particular create emissions that reduce the net sequestration value of product carbon storage. However, attempts to quantify the effects of transport on the net effects of forest management typically use relatively sparse survey data to determine transportation emission factors. We developed an approach for systematically determining transport emissions using: 1) -remotely sensed maps to estimate the spatial distribution of harvests, and 2) - industry data to determine landscape-level harvest volumes as well as the location and processing totals of individual mills. These data support spatial network analysis that can produce estimates of fossil carbon released in timber transport.     

Mapping of aggregated floodplain plant communities using image fusion of CASI and LiDAR data

Combined optical and laser altimeter data offer the potential to map and monitor plant communities based on their spectral and structural characteristics. A problem unresolved is, however, that narrowly defined plant communities, i.e. plant communities at a low hierarchical level of classification in the Braun-Blanquet system, often cannot be linked directly to remote sensing data for vegetation mapping. We studied whether and how a floristic dataset can be aggregated into a few major discrete, mappable classes without substantial loss of ecological meaning. Multi-source airborne data (CASI and LiDAR) and floristic field data were collected for a floodplain along the river Waal in the Netherlands. Mapping results based on floristic similarity alone did not achieve highest levels of accuracy. Ordination of floristic data showed that terrain elevation and soil moisture were the main underlying environmental drivers shaping the floodplain vegetation, but grouping of plant communities based on their position in the ordination space is not always obvious. Combined ordination-based grouping with floristic similarity clustering led to syntaxonomically relevant aggregated plant assemblages and yielded highest mapping accuracies.     

Using spectral information from the NIR water absorption features for the retrieval of canopy water content

Canopy water content (CWC) is important for mapping and monitoring the condition of the terrestrial ecosystem. Spectral information related to the water absorption features at 970 nm and 1200 nm offers possibilities for deriving information on CWC. In this study, we compare the use of derivative spectra, spectral indices and continuum removal techniques for these regions. Hyperspectral reflectance data representing a range of canopies were simulated using the combined PROSPECT + SAILH model. Best results in estimating CWC were obtained by using spectral derivatives at the slopes of the 970 nm and 1200 nm water absorption features. Real data from two different test sites were analysed. Spectral information at both test sites was obtained with an ASD FieldSpec spectrometer, whereas at the second site HyMap airborne imaging spectrometer data were also acquired. Best results were obtained for the derivative spectra. In order to avoid the potential influence of atmospheric water vapour absorption bands the derivative of the reflectance on the right slope of the canopy water absorption feature at 970 nm can best be used for estimating CWC.