Efficient GOCE satellite gravity field recovery based on least-squares using QR decomposition

We develop and apply an efficient strategy for Earth gravity field recovery from satellite gravity gradiometry data. Our approach is based upon the Paige-Saunders iterative least-squares method using QR decomposition (LSQR). We modify the original algorithm for space-geodetic applications: firstly, we investigate how convergence can be accelerated by means of both subspace and block-diagonal preconditioning. The efficiency of the latter dominates if the design matrix exhibits block-dominant structure. Secondly, we address Tikhonov-Phillips regularization in general. Thirdly, we demonstrate an effective implementation of the algorithm in a high-performance computing environment. In this context, an important issue is to avoid the twofold computation of the design matrix in each iteration. The computational platform is a 64-processor shared-memory supercomputer. The runtime results prove the successful parallelization of the LSQR solver. The numerical examples are chosen in view of the forthcoming satellite mission GOCE (Gravity field and steady-state Ocean Circulation Explorer). The closed-loop scenario covers 1 month of simulated data with 5 s sampling. We focus exclusively on the analysis of radial components of satellite accelerations and gravity gradients. Our extensions to the basic algorithm enable the method to be competitive with well-established inversion strategies in satellite geodesy, such as conjugate gradient methods or the brute-force approach. In its current development stage, the LSQR method appears ready to deal with real-data applications.     

多尺度空间线状实体形状相似关系的表达与度量

本文对空间数据多尺度表达中线状实体综合前后的形状相似性进行了描述、分析和研究。通过调查统计,得出了衡量线状实体形状相似程度的两个重要参量SIMⅠ、SIMⅡ,基于这两个参量给出了空间线状实体形状相似性的定义以及在不同尺度下线状实体形状相似性的度量方法,并用实例进行验证。为空间数据多尺度表达时线状实体的形状表达提供了一种分析和评价手段。     

经验模式分解的尺度阈值滤波在变形分析中的应用

本文将信号处理领域的经验模式分解算法应用于变形信息的提取中,通过引入阈值函数,建立了基于经验模式分解的尺度阈值滤波模型,采用优化模型确定了经验模式分解的次数。分别通过模拟试验和实测数据与小波阈值法和多项式拟合法进行了比对,分析表明:在低噪声情况下,三种方法都有一定的滤波效果;在高噪声情况下,经验模式分解的尺度阈值滤波法具有与小波阈值去噪法相等的精度,而且瞬时强噪声识别能力更好,优于多项式拟合法。     

空间索引与多尺度表达的一体化模型研究

在分析各种常用索引的基础上,依据各自的适用范围,提出了对点对象建立格网索引,对线、面对象建立改进四叉树索引的组合优化策略。设计的改进四叉树索引避免了传统四叉树的重复索引,提高了数据访问的效率。     

Multi-Scale Texture Modeling

Patterns of crystallographic preferred orientation are referred to as texture. The specific subject of texture analysis is the experimental determination and interpretation of the statistical distribution of orientations of crystals within a specimen of polycrystalline material, which could be metals or rocks. The objective is to relate an observed pattern of preferred orientation to its generating processes and vice versa. In geosciences, texture of minerals in rocks is used to infer constraints on their tectono-metamorphic history. Since most physical properties of crystals, such as elastic moduli, the coefficients of thermal expansion, or chemical resistance to etching depends on crystal symmetry and orientation, the presence of texture imparts directional properties to the polycrystalline material. A major issue of mathematical texture analysis is the resolution of the inverse problem to determine a reasonable orientation density function on SO(3) from measured pole intensities on , which relates to the inverse of the totally geodesic Radon transform. This communication introduces a wavelet approach into mathematical texture analysis. Wavelets on the two-dimensional sphere and on the rotational group SO(3) are discussed, and an algorithms for a wavelet decomposition on both domains following the ideas of Ta-Hsin Li is given. The relationship of these wavelets on both domains with respect to the totally geodesic Radon transform is investigated. In particular, it is shown that the Radon transform of these wavelets on SO(3) are again wavelets on . A novel algorithm for the inversion of experimental pole intensities to an orientation density function based on this relationship is developed.     

电子地图中多尺度地图数据显示的研究

将计算机三维模型简化中的LOD思想引入到电子地图的组织和显示中,提出了同一比例尺数字地图不同详细程度的分层显示方法和不同比例尺数字地图嵌套显示方法,实现了电子地图的多级缩放。     

单频GPS整周模糊度动态快速求解的研究

提出了一种动态快速求解整周模糊度的方法,即先对系数阵进行QR分解,然后通过矩阵变换使模糊度参数和位置参数分离,从而降低矩阵的维数,满足实时动态求解的要求,最后应用LAMBDA方法搜索模糊度。为验证该算法,用单频GPS接收机进行了静态定位和动态定位两种试验。结果表明,静态定位误差小于1cm,动态定位误差小于3cm。由此可见,该方法能够为动态用户快速而精确地实施GPS单频动态定位。     

东亚地区春冬季与夏季大气环流异常相互关系的研究

利用1951~2000年共50年的北半球500 hPa月平均高度距平场资料和奇异值分解技术（SVD），重点对东亚地区季节间大气环流异常的相互关系进行了初步探讨。结果表明，东亚地区季节间大气环流异常存在着较为密切的关联，并且这种明显的非同步联系具有时空相关显著的特点，尤其是夏季大气环流异常与其前冬和前春大气环流异常的联系更为密切。当前冬和前春北半球东亚大槽和北美大槽及蒙古高压偏强（或偏弱），极涡偏弱（或偏强），中高纬度盛行经向环流（或纬向环流），以及低纬和热带地区高度正距平（或负距平）明显时，则对应于夏季东亚地区西太平洋副高和鄂霍次克海阻高强度偏强（或偏弱），位置偏南（或偏北），贝加尔湖阻高强度也偏强（或偏弱），但位置偏西（或偏东）的大尺度环流形势出现。当春季北半球大气环流具有上述特点以及夏季鄂霍次克海阻高和西太平洋副高强度偏强（或偏弱），位置偏南（或偏北），且极涡较弱（较强）时，则东亚地区秋季大气环流对应于蒙古高压加强（或较弱），西太平洋副高减弱（或加强），并向南和向东移动（或移动较慢），极涡向南扩散（或扩散减弱），大气环流向冬季过渡加快（或减慢）。另外，大气环流异常还具有一定的持续性特征。     

A recursive decomposition algorithm for network seismic reliability evaluation

A new probabilistic analytical approach to evaluate seismic system reliability of large lifeline systems is presented in this paper. The algorithm takes the shortest path from the source to the terminal of a node weight or edge weight network as decomposition policy, using the Boolean laws of set operation and probabilistic operation principal, a recursive decomposition process then could be constructed. For a general weight network, the modified Torrieri method (NTR/T method) is introduced to combine with the suggested algorithm. Therefore, the recursive decomposition algorithm may be applied to evaluate the seismic reliability of general lifeline systems. A series of case studies, including a practical district electric power network system and a large urban water supply system, show that the suggested algorithm supplies a useful probabilistic analysis means for the seismic reliability evaluation of large lifeline systems. Copyright © 2002 John Wiley & Sons, Ltd.     

A multi-scale study on land use and land cover quality change: The case of the Yellow River Delta in China

This paper presents a case study of the Yellow River Delta in China, to trace land use and land cover changes during the past 20 years, with an emphasis on land quality changes. Three sets of data are used in this case study: remote sensing data derived from satellite images; crop yield data from statistics; and soil data collected by the researchers in the field. Our study reveals that at the regional scale, LUCC has taken place in a positive direction: vegetation cover has been expanding and crop yields per hectare have been on rise. However, while the overall eco-environment has improved, the improvement is uneven across the Delta region. At local levels, some areas show signs of increased salinization and declining organic content. Both natural forces and human activities are responsible for the LUCC, but human activities play a more important role. While some impacts of human activities are positive, the damages are often long-lasting and irreversible. We also conclude that it is necessary to use both macro data (such as remote sensing data) and micro data (data collected in the field) to study land quality change. The former are efficient in examining land quality changes at the regional scale, the latter can serve to verify ground patterns revealed from macro data and help to identify local variations, so as to get a comprehensive understanding of LUCC and promote sustainable land use and land management. This revised version was published online in July 2006 with corrections to the Cover Date.