Non-tidal variations in the deflection of the vertical at Beijing Observatory

At Beijing Observatory both astrometric and gravimetric observations are available to study the non-tidal variations in the deflections of the vertical (or non-tidal plumbline variations, PLVs). From repeated gravimetric observations performed in a network around the observatory, the PLVs at Beijing Observatory during the period 1987.0–1996.0 have been calculated. After comparison with the observational residuals (which also contain the PLV components) of the photoelectric astrolabe located at the observatory, the accuracy of the obtained PLV results has been examined. It is shown that, due to the asymmetry of the gravimetric network, the qualities of the two different PLV components are not equal. The longitude component of the PLV at Beijing has been determined successfully, to be of the order of 0.05, with a contribution of about 0.02 in the inter-annual time scale. The result for the latitude PLV component is not good enough to draw a conclusion. Although both techniques are able to measure the PLV, the result of the determination depends very much on the availability of observational data.     

The hierarchical tessellation model and its use in spatial analysis

The hierarchid tessellation model belongs to a class of spatial data models based on the recursive decomposition of space. The quadtree is one such tessellation and is characterized by square cells and a 1:4 decomposition ratio. To relax these constraints in the tessellation, a generalized hierarchical tessellation data model, called Adaptive Recursive Tessellations (ART), has been proposed. ART increases flexibility in the tessellation by the use of rectangular cells and variable decomposition ratios. In ART, users can specify cell sizes which are intuitively meaningful to their applications, or which can reflect the scales of data. ART is implemented in a data structure called Adaptive Recursive Run-Encoding (ARRE), which is a variant of two-dimensional run-encoding whose running path can vary with the different tessellation structures incorporated in an ART model. Given the recognition of the benefits of implementing statistical spatial analysis in GIS, the use of hierarchical tessellation models such as ART in spatial analysis is discussed. Three examples are introduced to show how ART can: (1) be applied to solve the quadrat size problem in quadrat analysis of point patterns; (2) act as the data model in the variable resolution block kriging technique for geostatistical data to reduce variation in kriging error; and (3) facilitate the evaluation of spatial autocorrelation for area data at multiple map resolutions via the construction of a connectivity matrix for calculating spatial autocorrelation indices based on ARRE.     

Delineating agricultural field boundaries from TM imagery using dyadic wavelet transforms

Agricultural field boundary information is important in the study of agricultural crops by remote sensing. In this paper, the dyadic wavelet transform is adopted to extract field boundaries directly from TM imagery. A field boundary is a set of locally connected pixels characterized by their abrupt spectral intensity variation (singularity) on the image. Preliminary extraction is made on the near-infrared Band 4 (0.76–0.90 μm) image. It is found that most of the boundaries are well characterized by the wavelet modulus maxima curves, particularly for fully vegetated fields, while better results are obtained from the NDVI image for cropped field boundaries. To extract additional boundaries while suppressing unwanted noisy edges, the wavelet transform is then applied to the Wetness image calculated by the Tasselled Cap Transform. A technique is implemented to integrate the edges extracted from each source. The combined results show that the majority of the field boundaries are delineated. The use of a single date image fails to delineate some of the fallow field boundaries and boundaries between two cropped fields with similar spectral properties. Although techniques have been applied to reduce noise, still great difficulties are encountered when trying to delineate edges that are field boundaries from all the detected singularities. Manual editing is still required to trace the remaining boundaries to obtain a complete boundary map. It is recommended that a multi-temporal data set be used to obtain a complete delineation with a fully automated algorithm.     

Response of GPS occultation signals to atmospheric gravity waves and retrieval of gravity wave parameters

We show that the amplitude of the Global Positioning System (GPS) signals in the radio occultation (RO) experiments is an indicator of the activity of the gravity waves (GW) in the atmosphere. The amplitude of the GPS RO signals is more sensitive to the atmospheric wave structures than is the phase. Early investigations used only the phase of the GPS occultation signals for statistical investigation of the GW activity in the height interval 10–40 km on a global scale. In this study, we use the polarization equations and Hilbert transform to find the 1-D GW radio image in the atmosphere by analyzing the amplitude of the RO signal. The radio image, also called the GW portrait, consists of the phase and amplitude of the GW as functions of height. We demonstrate the potential of this method using the amplitude data from GPS/Meteorology (GPS/MET) and satellite mission Challenge Mini-satellite Payload (CHAMP) RO events. The GW activity is nonuniformly distributed with the main contribution associated with the tropopause and the secondary maximums related to the GW breaking regions. Using our method we find the vertical profiles of the horizontal wind perturbations and its vertical gradient associated with the GW influence. The estimated values of the horizontal wind perturbations are in fairly good agreement with radiosonde data. The horizontal wind perturbations v(h) are ±1 to ±5 m s with vertical gradients dv/dh ±0.5 to ±15 m s km at height 10–40 km. The height dependence of the GW vertical wavelength was inferred through the differentiation of the GW phase. Analysis of this dependence using the dispersion relationship for the GW gives the estimation of the projection of the horizontal background wind velocity on the direction of the GW propagation. For the event considered, the magnitude of this projection changes between 1.5 and 10 m s at heights of 10–40 km. We conclude that the amplitude of the GPS occultation signals contain important information about the wave processes in the atmosphere on a global scale.     

Urban change detection related to earthquakes using an adaptive nonlinear mapping of high-resolution images

We propose a nonlinear mapping method for detecting geographical changes. In our system, two images photographed at different points of time are checked based on an iterative nonlinear mapping. The system automatically extracts feature changes of aerial imageries by computing distribution of the image matching score. We evaluate the change detection ability by using the receiver operating characteristic under misregistration cases, including rotational misalignments.     

Evaluation of Information Loss in Digital Elevation Models With Digital Photogrammetric Systems

Information loss is caused when a surface is sampled with a finite interval, such as in the production of a digital elevation model (DEM). This information loss can become the dominant part of the error in a DEM. The ability to quantify information loss enables guidance to be provided for an appropriate choice of grid interval and better accuracy assessment for the DEM. With the use of digital photogrammetric systems, evaluation of information loss has become much easier. This paper describes three methods of evaluating information loss. An example is given of the method which is most appropriate for use with a digital photogrammetric system, based on rock cliff surface data and the VirtuoZo system.     

Maximum-likelihood ambiguity resolution based on Bayesian principle

 Based on the Bayesian principle and the fact that GPS carrier-phase ambiguities are integers, the posterior distribution of the ambiguities and the position parameters is derived. This is then used to derive the maximum posterior likelihood solution of the ambiguities. The accuracy of the integer ambiguity solution and the position parameters is also studied according to the posterior distribution. It is found that the accuracy of the integer solution depends not only on the variance of the corresponding float ambiguity solution but also on its values. Received: 27 July 1999 / Accepted: 22 November 2000     

A statistical analysis applied to the Rf/φ method

A statistical analysis for the R_f/φ method is developed in the paper. The technique provides estimates of R_s, R_i and θ according to R_f and φ. Moreover, we can restore the pattern of initial orientations of marker objects by means of their frequency histogram, provided the strain is homogeneous. It is unnecessary to postulate random initial orientations if the principal directions of strain ellipsoid are known.