基于规则库和网格算法的土地利用现状图自动数字注记

图形自动注记可极大地提高制图的工作效率。针对土地利用现状图数字注记的特点,本文从地图认知角度出发,建立注记自动配置的规则库,提出面状要素自动注记的一种新的算法——网格法。该算法最大的优点是容易实现、效率高,可方便地解决注记压盖冲突等问题。通过开发一个程序实例,对图斑要素进行了自动数字注记,实验表明95％以上的图形注记效果都较为理想。     

Aspherical correction for free vibration of elastically anisotropic solid by means of theXYZ algorithm

By applying the perturbation theory to theXYZ algorithm (a kind of variational method), the difference f in free vibration frequencies between sphere and ellipsoid was approximated as , where i and _i (i = x,y andz) (i=x, y andz) are aspherical coefficients and asphericities of the ellipsoid, respectively. We developed an analytic method to compute the aspherical coefficients719-4 by using theXYZ algorithm. A numerical example was given for an ellipsoidal olivine, and an attempt was made to estimate the asphericities of the specimen by a least-squares method, based on the relationship between frequency shift and asphericity.     

Dynamics of the mid-ocean ridge plate boundary: Recent observations and theory

The global mid-ocean ridge system is one of the most active plate boundaries on the earth and understanding the dynamic processes at this plate boundary is one of the most important problems in geodynamics. In this paper I present recent results of several aspects of mid-ocean ridge studies concerning the dynamics of oceanic lithosphere at these diverging plate boundaries. I show that the observed rift valley to no-rift valley transition (globally due to the increase of spreading rate or locally due to the crustal thickness variations and/or thermal anomalies) can be explained by the strong temperature dependence of the power law rheology of the oceanic lithosphere, and most importantly, by the difference in the rheological behavior of the oceanic crust from the underlying mantle. The effect of this weaker lower crust on ridge dynamics is mainly influenced by spreading rate and crustal thickness variations. The accumulated strain pattern from a recently developed lens model, based on recent seismic observations, was proposed as an appealing mechanism for the observed gabbro layering sequence in the Oman Ophiolite. It is now known that the mid-ocean ridges at all spreading rates are offset into individual spreading segments by both transform and nontransform discontinuities. The tectonics of ridge segmentation are also spreading-rate dependent: the slow-spreading Mid-Atlantic Ridge is characterized by distinct bulls-eye shaped gravity lows, suggesting large along-axis variations in melt production and crustal thickness, whereas the fast-spreading East-Pacific Rise is associated with much smaller along-axis variations. These spreading-rate dependent changes have been attributed to a fundamental differences in ridge segmentation mechanisms and mantle upwelling at mid-ocean ridges: the mantle upwelling may be intrinsically plume-like (3-D) beneath a slow-spreading ridge but more sheet-like (2-D) beneath a fast-spreading ridge.     

Textural segmentation of digital rock images into bedding units using texture energy and cluster labels

The texture of digital rock images, as recorded, for instance, with borehole imaging devices, is shown to reflect different bedding types. Textural segmentation of borehole images, therefore, subdivides the recorded sequence into bedding units. We show that a textural segmentation algorithm based on the concept of texture energy achieves good results when compared with synthetic as well as real data in which petroleum geologists have performed zonations on cores. Texture energy involves filtering of the original image with a set of texture sensitive masks. The filtering is done as a finite convolution over the size of the masks. On the resulting images the variance is computed over a relatively large sliding window, which, in its practical implementation, covers the full width of the image. The resulting nine one-dimensional curves are then clustered hierarchically into a user-determined number of image texture or lithological bedding classes. Principal component analysis previous to clustering can be used to reduce redundancy in the data. A recurring and relatively ill-defined problem in this field are macro-textures, i.e., the cyclic interbedding of two or more bedding types. We show that sliding Fourier transforms and variable mask scale can successfully address the zonation of macro-textures. In general, the method gives best results with mask sizes equivalent to 2–4 centimeters, reflecting the length scale at which the investigated geological bedding seems to have its highest variation.     

Estimation for partially observed Markov processes

Many stochastic process models for environmental data sets assume a process of relatively simple structure which is in some sense partially observed. That is, there is an underlying process (X_n, n 0) or (X_t, t 0) for which the parameters are of interest and physically meaningful, and an observable process (Y_n, n 0) or (Y_t, t 0) which depends on the X process but not otherwise on those parameters. Examples are wide ranging: the Y process may be the X process with missing observations; the Y process may be the X process observed with a noise component; the X process might constitute a random environment for the Y process, as with hidden Markov models; the Y process might be a lower dimensional function or reduction of the X process. In principle, maximum likelihood estimation for the X process parameters can be carried out by some form of the EM algorithm applied to the Y process data. In the paper we review some current methods for exact and approximate maximum likelihood estimation. We illustrate some of the issues by considering how to estimate the parameters of a stochastic Nash cascade model for runoff. In the case of k reservoirs, the outputs of these reservoirs form a k dimensional vector Markov process, of which only the kth coordinate process is observed, usually at a discrete sample of time points.     

Reservoir description using dynamic parameterisation selection with a combined stochastic and gradient search

There is a correspondence between flow in a reservoir and large scale permeability trends. This correspondence can be derived by constraining reservoir models using observed production data. One of the challenges in deriving the permeability distribution of a field using production data involves determination of the scale of resolution of the permeability. The Adaptive Multiscale Estimation (AME) seeks to overcome the problems related to choosing the resolution of the permeability field by a dynamic parameterisation selection. The standard AME uses a gradient algorithm in solving several optimisation problems with increasing permeability resolution. This paper presents a hybrid algorithm which combines a gradient search and a stochastic algorithm to improve the robustness of the dynamic parameterisation selection. At low dimension, we use the stochastic algorithm to generate several optimised models. We use information from all these produced models to find new optimal refinements, and start out new optimisations with several unequally suggested parameterisations. At higher dimensions we change to a gradient-type optimiser, where the initial solution is chosen from the ensemble of models suggested by the stochastic algorithm. The selection is based on a predefined criterion. We demonstrate the robustness of the hybrid algorithm on sample synthetic cases, which most of them were considered insolvable using the standard AME algorithm.     

Temporal behavior of the background seismicity rate in central Japan, 1998 to mid-2003

We studied the temporal behavior of the background shallow seismicity rate in 700 circular areas across inland Japan. To search for and test the significance of the possible rate changes in background seismicity, we developed an efficient computational method that applies the space–time ETAS model proposed by Ogata in 1998 to the areas. Also, we conducted Monte Carlo tests using a simulated catalog to validate the model we applied. Our first finding was that the activation anomalies were found so frequently that the constant background seismicity hypothesis may not be appropriate and/or the triggered event model with constraints on the parameters may not adequately describe the observed seismicity. However, quiescence occasionally occurs merely by chance. Another outcome of our study was that we could automatically find several anomalous background seismicity rate changes associated with the occurrence of large earthquakes. Very significant seismic activation was found before the M6.1 Mt. Iwate earthquake of 1998. Also, possible seismic quiescence was found in an area 150 km southwest of the focal region of the M7.3 Western Tottori earthquake of 2000. The seismicity rate in the area recovered after the mainshock.     

角度格式转换过程中可靠算法研究

提出了一种可靠算法是将原始度.分秒加上一个与之同符号,且不影响计算精度的有效微小量ｅ,然后再进行格式转换,能有效减少工程返工和损失。     

基于可分解马尔科夫网的极端椒盐噪声图像滤波

提出了一种基于可分解马尔科夫网(decomposableMarkovnetworks,DMN)的极端椒盐噪声的均值滤波方法,指出其网络节点的阈值衰减特性和网络节点的连接特性具有很好的对椒盐噪声污染图像的噪声定位的作用,并提出一种在该网络控制下,只对与噪声相关的像素进行均值计算以替代噪声像素的亚均值滤波算法,实现了图像的较强自适应滤波。实验表明,本文方法具有良好的滤波性能。