Supervised SAR Image MPM Segmentation Based on Region-Based Hierarchical Model

This letter presents a novel method of supervised multiresolution segmentation for synthetic aperture radar images. The method uses a region-based half-tree hierarchical Markov random field model for multiresolution segmentation. To form the region-based multilayer model, the watershed algorithm is employed at each resolution level independently. The nodes of a quadtree in the proposed model are defined as regions instead of pixels. The relationship over scale is studied, and the region-based upward and downward maximization of posterior marginal estimations are deduced. The experimental results for the segmentation of homogeneous areas prove the region-based model much better in terms of robustness to speckle and preservation of edges compared to the pixel-based hierarchical model and the Gibbs sampler with the single-resolution model     

Synergetic efficiency of Lidar and WorldView-2 for 3D urban cartography in Northeast Mexico

Three-dimensional urban cartography is needed for city changes’ assessment. The variety of studies using 3D calculations of urban elements grows each year. Building and vegetation volumes are necessary to assess and understand spatio-temporal urban changeable environments. However, there are technical questions as to which method can improve 3D urban cartographic accuracy. The innovative part of this current study is the creation of a six-band hybrid obtained from LIDAR and WorldView2 synergy. Two different enhancement algorithms demonstrated the most important spectral features for the urban development and vegetation classes. Results indicated an improvement in accuracy by up to 21.3%, according to the Kappa coefficient. Both infra-red band and intensity band were the most significant, according to the principal components analysis. The synergy delimited classes and polygons, as well as the direct display of information regarding heights of elements and improving the extraction of roads, buildings and vegetation classes.     

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.     

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.     

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.     

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.     

Attenuation differences in layer 2A in intermediate- and slow-spreading oceanic crust

In situ seismic attenuationQ⁻¹logs are derived from borehole velocity profiles and reveal sharp boundaries between morphologies of the extrusive volcanic layers in intermediate- and slow-spreading oceanic crust.Q⁻¹logs are calculated from the scattering attenuation associated with vertical velocity heterogeneity in Ocean Drilling Program Holes 504B and 896A and in Hole 395A, located in 5.9–7.3 Ma crust on the Pacific and Atlantic plates, respectively. Our results strongly tie crustal properties to seismic measurables and observed geological structures: we find that the scattering attenuation can be used to identify the extrusive volcanic sequence because it is closely related to changes in the degree of vertical heterogeneity. We interpret a distinct decrease in the Q⁻¹log at the transition below the extrusive volcanic layer to correspond with the seismic layer 2A/2B boundary. The boundary is located at 465 m depth below the sea floor in both Hole 395A and 504B, although this is likely to be a coincidence of the sediment thickness at these sites. Layer 2A is estimated to be approximately 150 m thick in Hole 504B and > 300 m thick in Hole 395A. Cyclic sequences of high-porosity pillows and low-porosity massive units in the uppermost 100 m of volcanics in Hole 395A result in large velocity heterogeneities which cause > 5 times more attenuation in this layer than in Hole 504B. In Hole 896A, by contrast, fewer pillows, more massive flows, and a greater volume of carbonate veins decrease the velocity heterogeneity and attenuation significantly over only 1 km distance from Hole 504B. We conclude that the attenuation in the extrusive volcanics of the ocean crust is largely controlled by variation in local heterogeneity and morphology as well as by subsequent hydrothermal alteration. The observed differences inQ⁻¹profiles and layer 2A thickness at these sites may be attributed to variations in the volume and duration of volcanic activity at mid-ocean spreading centers for these Pacific and Atlantic ridge segments.     

Regionalization of unit hydrograph parameters: 1. Comparison of regression analysis techniques

Hydrologic regionalization is a useful tool that allows for the transfer of hydrological information from gaged sites to ungaged sites. This study developed regional regression equations that relate the two parameters in Nash's IUH model to the basin characteristics for 42 major watersheds in Taiwan. In the process of developing the regional equations, different regression procedures including the conventional univariate regression, multivariate regression, and seemingly unrelated regression were used. Multivariate regression and seeming unrelated regression were applied because there exists a rather strong correlation between the Nash's IUH parameters. Furthermore, a validation study was conducted to examine the predictability of regional equations derived by different regression procedures. The study indicates that hydrologic regionalization involving several dependent variables should consider their correlations in the process of establishing the regional equations. The consideration of such correlation will enhance the predictability of resulting regional equations as compared with the ones from the conventional univariate regression procedure.     

Afternoon mid-latitude current system and low-latitude geomagnetic field asymmetry during geomagnetic storms

For four geomagnetic storms of middle intensity the relationship between the low-latitude magnetic field asymmetry using ASY indices and the intensity of the auroral eastward and westward electro-jet was considered. It was asked whether there exists a connection between ASY and the eastward electrojet. To answer this question equivalent current systems were estimated in mid-latitudes. It was found that the observations obviously show no correlative relationship between the low-latitude magnetic-field asymmetry and the eastward electrojet, whereas one exists between ASY and the westward electrojet. To explain the generally accepted common three-dimensional current system between the partial ring current and the eastward electrojet, a condensor model of the three-dimensional current system was developed. It could be shown that the short periodic variations of the partial ring current are shielded by the condensor and cannot influence the eastward-electrojet current.