An improved building detection approach using L-band POLSAR two-dimensional time-frequency decomposition over oriented built-up areas

The objective of this study is to efficiently extract detailed information about various man-made targets in oriented built-up areas using polarimetric synthetic aperture radar (POLSAR) images. This paper develops an improved approach for building detection by utilizing Two-Dimensional Time-Frequency (2-D TF) decomposition. This method performs outstandingly in distinguishing between man-made and natural targets based on the isotropic behaviors, frequency-sensitive responses, and scattering mechanisms of objects. The proposed method can preserve the spatial resolution and exploit the advantages of TF decomposition; specifically, the exact outlines of buildings can be effectively located, and more types of features (e.g., flat roofs, roads, and walls that are oblique to the radar illumination) can be distinguished from forests in complex built-up areas by 2-D TF decomposition. The coarser-resolution subaperture images that are produced in the azimuth direction, which correspond to different looking angles, are beneficial for detecting man-made structures with main scattering centers oriented at oblique angles with respect to the radar illumination. In the range direction, the obtained subaperture images, which correspond to various observation frequencies, can be helpful in distinguishing flat roofs and roads from forests. This method was successfully implemented to analyze both NASA/JPL L-band AIRSAR and L-band EMISAR data sets. The building detection results of the proposed method exhibit a significant improvement over those of other methods and reach an overall accuracy over 80%, with approximately 20% higher than the accuracies of K-means clustering and the entropy/alpha-Wishart classifier and approximately 10% higher than the accuracy of the support vector machine method. Moreover, building details can be precisely detected, obliquely oriented buildings can be identified, and the distinction between buildings and forests is significantly improved, as both visually and statistically indicated. This method is highly adaptable and has substantial application value.     

改进分数阶达尔文粒子群的多Renyi熵图像分割算法

针对智能优化图像分割算法易陷入局部最优、分割精度不高等问题,本文融合改进的分数阶达尔文粒子群算法和二维Renyi熵多阈值,提出了一种新的多阈值遥感图像分割算法。算法利用粒子自身进化信息来定义进化因子,结合进化因子并利用高斯图函数调整分数阶次a系数以实现精确计算和快速收敛;根据局部最优概率因子对局部最优位置进行Levy飞行随机扰动以提高算法跳出局部最优的能力;同时将二维Renyi熵单阈值扩展到多阈值分割上,并结合改进的分数阶达尔文粒子群算法,将二维Renyi熵多阈值应用于遥感图像分割中仿真结果表明,与其他2种智能优化分割算法相比,本文分割算法在细节处理和分割精度上均有明显优势,在PRI上至少提升7.27%、VOI至少降低6.5%、GCE至少降低10.4%.     

准噶尔盆地中部永进地区走滑断裂发育特征及成因物理模拟

永进地区位于准噶尔盆地中部,最近发现了多个与走滑断层相关的含油气构造,但关于走滑断层的发育特征及成因机制研究程度不够深入。本文通过三维地震资料精细解释,在研究区三叠系—侏罗系内识别出近东西向、北西西向以及北东东向的三组走滑断裂体系,平面上呈“网格状”展布,剖面上具有不同深度几何学形态差异展布特征。在此基础上基于相似性原理设计四组砂箱模拟对比实验,重现研究区构造演化过程。模拟结果表明,这类走滑断裂的形成与基底先存断层的发育位置有关,是受先存构造和地层属性双重控制的广布式走滑断裂系统,从而建立了研究区的断裂系统成因模式。研究成果对具有相似地质背景地区的走滑断裂成因解释具有借鉴意义。     

Water and salt balance modelling to predict the effects of land-use changes in forested catchments. 1. Small catchment water balance model

A long-term water balance model has been developed to predict the hydrological effects of land-use change (especially forest clearing) in small experimental catchments in the south-west of Western Australia. This small catchment model has been used as the building block for the development of a large catchment-scale model, and has also formed the basis for a coupled water and salt balance model, developed to predict the changes in stream salinity resulting from land-use and climate change. The application of the coupled salt and water balance model to predict stream salinities in two small experimental catchments, and the application of the large catchment-scale model to predict changes in water yield in a medium-sized catchment that is being mined for bauxite, are presented in Parts 2 and 3, respectively, of this series of papers. The small catchment model has been designed as a simple, robust, conceptually based model of the basic daily water balance fluxes in forested catchments. The responses of the catchment to rainfall and pan evaporation are conceptualized in terms of three interdependent subsurface stores A, B and F. Store A depicts a near-stream perched aquifer system; B represents a deeper, permanent groundwater system; and F is an intermediate, unsaturated infiltration store. The responses of these stores are characterized by a set of constitutive relations which involves a number of conceptual parameters. These parameters are estimated by calibration by comparing observed and predicted runoff. The model has performed very well in simulations carried out on Salmon and Wights, two small experimental catchments in the Collie River basin in south-west Western Australia. The results from the application of the model to these small catchments are presented in this paper.     

APPLICATION OF A GROUPED RESPONSE UNIT HYDROLOGICAL MODEL TO A NORTHERN WETLAND REGION

Applications of hydrological models to northern wetland-dominated regions have been limited in the past to a few case studies on small basins employing ‘lumped’ models. Only recently have there been attempts to apply the grouped response unit (GRU) distributed modelling approach using terrain classifications to these same basins. This study summarizes recent efforts in applying such a model. For the purposes of implementing the GRU approach, terrain types that are hydrologically significant and characteristic to the wetland-dominated regime were successfully discriminated using a principal component analysis and a hybrid unsupervised/supervised classification technique on Landsat–Thematic Mapper imagery. The terrain classifications were then used as input into a distributed hydrological model for calibration and validation using recorded spring runoff events. Preliminary model applications and results are described. Calibration to a historic spring runoff event yielded an r² value of 0.86. Model validation, however, yielded much poorer results. The problems of model applicability to this region and limitations of sparse data networks are highlighted. The need for more field research in this type of hydrological regime, and associated improvements to the model parameter set are also identified.     

MEASURING AND MODELLING OF SOIL WATER DYNAMICS AND RUNOFF GENERATION IN AN AGRICULTURAL LOESSIAL HILLSLOPE

Surface runoff may be generated when the rainfall intensity exceeds the infiltration capacity, or when the soil profile is saturated with water. Indications exist that both types of overland flow may occur in hilly agricultural loess regions. Here, for a loessial hillslope under maize in the southern part of The Netherlands, it was shown, with pressure head and runoff measurements, that Hortonian overland flow occurs during typical summer rain events. Surface runoff was initiated after saturation of the top 5–10 cm of the soil. Deeper in the soil, unsaturated conditions prevailed while runoff took place. Peak runoff discharges at the outlet of the subcatchment occurred a few minutes after peak rainfall intensities were measured. It appeared that SWMS_2D, a two-dimensional water flow model, was capable in simulating observed pressure head changes and runoff. Simulated potential runoff for the transect studied was higher by a magnitude of three than the measured areal average. This indicates effects of surface ponding, and the probable location of this particular transect in a region with high runoff production.     

FLOOD MODELLING WITH GIS-DERIVED DISTRIBUTED UNIT HYDROGRAPHS

The concept of a spatially distributed unit hydrograph is based on the fact that the unit hydrograph can be derived from the time–area curve of a watershed by the S-curve method. The time–area diagram is a graph of cumulative drainage area contributing to discharge at the watershed outlet within a specified time of travel. Accurate determination of the time–area diagram is made possible by using a GIS. The GIS is used to describe the connectivity of the links in the watershed flow network and to calculate distances and travel times to the watershed outlet for various points within the watershed. Overland flow travel times are calculated by the kinematic wave equation for time to equilibrium; channel flow times are based on the Manning and continuity equations. To account for channel storage, travel times for channel reaches are increased by a percentage depending on the channel reach length and geometry. With GIS capability for rainfall mapping, the assumption of a uniform spatial rainfall distribution is no longer necessary; hence the term, spatially distributed unit hydrograph. An example of the application for the Waiparous Creek in the Alberta Foothills is given. IDRISI is used to develop a simple digital elevation model of the 229 km² watershed, using 1 km × 1 km grid cells. A grid of flow directions is developed and used to create an equivalent channel network. Excess rainfall for each 1 km × 1 km cell is individually computed by the Soil Conservation Service (SCS) runoff curve method and routed through the equivalent channel network to obtain the time–area curve. The derived unit hydrograph gave excellent results in simulating an observed flood hydrograph. The distributed unit hydrograph is no longer a lumped model, since it accounts for internal distribution of rainfall and runoff. It is derived for a watershed without the need for observed rainfall and discharge data, because it is essentially a geomorphoclimatic approach. As such, it allows the derivation of watershed responses (hydrographs) to inputs of various magnitudes, thus eliminating the assumption of proportionality of input and output if needed. The superposition of outputs is retained in simulating flood hydrographs by convolution, since it has been shown that some non-linear systems satisfy the principle of superposition. The distributed unit hydrograph appears to be a very promising rainfall runoff model based on GIS technology.     

起伏海面的雷电电磁场传播特征及其对闪电定位精度的影响

为了研究起伏海面对雷电电磁传播的影响,本文利用Barrick表面阻抗理论和Wait近似算法,采用改进二维分形海面模型模拟起伏海面,利用数值模式,分析起伏海面的雷电电磁传播特征,并进一步讨论了起伏海面对时差法闪电定位系统定位精度的影响。结果表明:起伏海面对垂直电场和磁场的峰值的影响不显著,但会引起波形的上升期时间的延长,浪高越大,影响越明显;随着观测距离的增加,雷电垂直电场波形的上升时间逐渐变长;风速的变化与雷电垂直电场波形的上升时间成正比;由于海面起伏引起雷电电磁场波形在传播中的变化会影响基于时差法闪电定位系统的定位精度,定位误差可达几至十几公里。     

基于纵横波保幅分离的粘弹介质弹性波正演模拟

地下介质往往表现为粘弹特性,研究基于粘弹假设的数值模拟方法对于正确认识地震波的传播规律和提高地震勘探精度具有重要意义。波动方程正演过程中纵横波的保幅解耦是研究准确的粘弹介质中地震波传播机理的前提,基于散度和旋度算子的纵横波解耦方法会使波场的相位和振幅产生畸变,且解耦后的波场在极性反转位置上无法与分离前混合波场各分量对应。在散度和旋度算子上再做一次梯度和旋度运算的波场分离方法虽然能够克服上述缺陷,但存在保幅性差等问题。本文从粘弹介质中的一阶速度-应力方程出发,推导了矢量纵横波分离的波数域表达式,结合有限差分思路给出了其在空间域的求解方法。本文方法利用纵横波的传播速度对现有的矢量波场分离方法进行振幅校正,并将校正结果分别作为纵波与横波对时间的二阶偏导,实现了粘弹介质中的纵横波分离。模型试算结果表明,本文方法能够克服现有方法的缺陷,获得更具保幅性的波场分离结果。     

Quantifying sand delivery to deep water during changing sea-level: Numerical models from the Quaternary Brazos Icehouse continental margin

Sequence stratigraphy for clastic continental margins predicts the development of sand-rich turbidite deposits during specific times in relation to base-level cycles. It is now widely understood that deltas can extend to the shelf-edge forced by high sediment flux and/or base level, providing a direct connection to transfer sediment and sand to the slope and basin floor even during high base level periods. Herein, we build a stratigraphic forward model for the last 120 kyr of the fluvio-deltaic to deep-water Brazos system (USA) where sediment partitioning along an Icehouse continental margin can be evaluated. The reduced-complexity stratigraphic forward model employs geologically constrained input parameters and mass balance. The modelled architecture is consistent with the location of depositional units previously mapped in the shelf. Sand bypasses the shelf and upper slope between 35 to 15 kyr before present and only about 20%–30% of all the sediment and sand supplied to the system is transferred to deep water. Several scenarios based on the initial Brazos model investigate the relationships between base level and deep-water sand ratio (DWSR). DWSR is defined as the relative amount of sand transferred to the deep-water portions of the system subdivided by the total sand input to the model. Linear correlations between DWSR and base level change rates or base level are very poor. Short-term variability due to local processes (for example avulsions) is superimposed to the long-term trends and mask the base level signal. DWSR for an entire base-level cycle is mainly controlled by the proportion of time the delta stays docked at the shelf-edge. Stratigraphic forward models are useful to complement field observations and quantify how different processes control stratigraphy, which is important for making predictions in areas with limited information.