基于泊松小波径向基函数融合多源数据的局部重力场建模

融合多源数据的高精度、高分辨率的局部重力场建模是物理大地测量学的前沿和热点问题.本文研究了基于径向基函数融合多源数据的局部重力场建模方法,利用Monte-Carlo方差分量估计实现了不同类型的观测数据的合理定权,引入了最小标准差法确定基函数的适宜网络,分析了地形因素对于基函数网络确定及局部重力场建模精度的影响.以泊松小波基函数为构造基函数,结合残差地形模型,融合实测的陆地重力异常、船载重力异常及航空重力扰动数据构建了局部区域陆海统一的似大地水准面模型.研究结果表明:引入残差地形模型平滑了地形质量引入的高频扰动信号,简化了基函数的网络设计;并提高了重力似大地水准面的精度,平原地区其精度提高了4mm,地形起伏较大的山区其精度提高了约5cm.总体而言,基于"三步法"构建的局部重力似大地水准面在荷兰、比利时及德国相关区域,其精度分别达到1.12cm、2.80cm以及2.92cm.     

闽东南沿海地区景观指数粒度效应的高分辨率遥感分析

以闽东南沿海地区3个流域SPOT 5影像解译所得的土地覆被类型分布图为原始矢量数据,通过优势规则聚合为一系列不同粒度大小的栅格类型图,统计各粒度下28种常用景观指数值;采用尺度检测图分析各种景观指数的粒度变化响应规律,对比中分辨率数据与高分辨率数据所建立尺度函数的差异,评价预测响应型指数尺度下推的相对误差。结果表明：闽东南沿海地区景观指数的粒度效应明显,可分为预测响应型、分段预测型、阶梯变化型和波动变化型4类,其中分段预测型指数由已有报道的预测响应型分化而来,其尺度转折点位于5 m 、7.5 m或10 m等粒度下,即高分辨率影像对景观指数粒度效应的分类有重要影响。在预测响应型指数中,高分辨率影像对ED、SHAPE_MN、CONTIG_MN和AI 4种指数的尺度函数影响较大,由中分辨率数据所建立的尺度函数进行尺度下推,会带来较大误差;这些指数的区域差异性也较大,对同一指数来说,不同区域可选用不同的尺度函数进行尺度下推。     

Development of design flood hydrographs using probability density functions

Probability density functions (PDFs) are used to fit the shape of hydrographs and have been popularly used for the development of synthetic unit hydrographs by many hydrologists. Nevertheless, modelling the shapes of continuous stream flow hydrographs, which are probabilistic in nature, is rare. In the present study, a novel approach was followed to model the shape of stream flow hydrographs using PDF and subsequently to develop design flood hydrographs for various return periods. Four continuous PDFs, namely, two parameter Beta, Weibull, Gamma and Lognormal, were employed to fit the shape of the hydrographs of 22 years at a site of Brahmani River in eastern India. The shapes of the observed and PDF fitted hydrographs were compared and root mean square errors, error of peak discharge (EQ_P) and error of time to peak (ET_P) were computed. The best‐fitted shape and scale parameters of all PDFs were subjected to frequency analysis and the quartiles corresponding to 20‐, 50‐, 100‐ and 200‐year were estimated. The estimated parameters of each return period were used to develop the flood hydrographs for 20‐, 50‐, 100‐ and 200‐year return periods. The peak discharges of the developed design flood hydrographs were compared with the design discharges estimated from the frequency analysis of 22 years of annual peak discharges at that site. Lognormal‐produced peak discharge was very close to the estimated design discharge in case of 20‐year flood hydrograph. On the other hand, peak discharge obtained using the Weibull PDF had close agreement with the estimated design discharge obtained from frequency analysis in case of 50‐, 100‐ and 200‐year return periods. The ranking of the PDFs based on estimation of peak of design flood hydrograph for 50‐, 100‐ and 200‐year return periods was found to have the following order: Weibull > Beta > Lognormal > Gamma. Copyright © 2009 John Wiley & Sons, Ltd.     

EXPLORATION ON METHOD OF AUTO-CLASSIFICATION FOR MAIN GROUND OBJECTS OF THREE GORGES RESERVOIR AREA

1IN TR O DU C TIO N With thedevelopmentof remotesensingtechnolog,y thestudyandapplicatioonftheautomaticclassifitciaon method usingcomputerhad achievedsatisfyi ndgevel- opment. From spectrumstatisticcallassificatitoonthe neural network classificatimoenthod…     

Modelling Thermohaline Properties in an Estuarine Upwelling Ecosystem (R?́a de Vigo: NW Spain) Using Box-Jenkins Transfer Function Models

Since 1987, twice weekly, hydrological variables have been monitored at a fixed station in the R?́a de Vigo (NW Spain), aiming to examine the time scales of variability and the relationships to meteorological conditions. The present paper analyses: (1) the advantage of Box-Jenkins transfer function (TF) models (single output–multiple input), a type of linear stochastic model, to describe the dynamic behaviour of the system; and (2) the coupling between the R?́a and meteorological events at the time scale of autonomy of this coastal inlet affected by the Iberian coastal upwelling, approximately a fortnightly period. In order to achieve these objectives, thermohaline properties have been used to characterize the estuarine ecosystem (output variables), while wind regime, runoff in the drainage basin and incoming solar radiation have been considered as the main forcing variables (input variables). The use of the amplitude time series, derived from principal component analysis (PCA) applied to the deseasonalized meteorological variables, is also explored as a different set of input variables.When compared with standard regression models, all TF models built to describe thermohaline behaviour had reduced residual variance. Similar TF models, as well as percentage of explained variance, were also obtained when meteorological variables or the amplitude time series were used as input variables. The fitted TF models provided an insight into the ‘ inertial ’ behaviour of the system and the time scales of coupling of the system with the forcing variables. The plausible physical mechanisms which link the response of the system with the observed meteorological variability are also discussed. As could be expected, bottom thermohaline properties show a stronger inertial behaviour than the surface ones, which is particularly marked for bottom temperature. Besides, the shelf domain, by means of upwelling-downwelling events, strongly influences surface and bottom temperature, as well as bottom salinity; by contrast, surface salinity is mainly influenced by the effect of wind along the main axis of the R?́a and runoff. In relation to the time scales of coupling between the system and the forcing variables, thermohaline properties show a dependance with the meteorological conditions in, at least, the immediately preceding fortnight period. It was concluded that: (1) TF models that incorporate meteorological information described the dynamic behaviour of the system adequately; and (2) this type of model can be useful as a first approximation to develop more sophisticated (deterministic) models, since, with the purpose of modelling any state variable of the system, both the coupling between different domains and the time scales of the interactions must be taken into account.     

BOLTZMANN方程与碰撞矩阵元的计算

作者用Ｂｕｒｎｅｔｔ函数作为展开基底表示分子速度分布函数和Ｂｏｌｔｚｍａｎｎ方程，并由此证明，Ｂｏｌｔｚｍａｎｎ碰撞矩阵元的计算归结为约化积分、Ｃｌｅｂｓｃｈ－Ｇｏｒｄａｎ系数和Ｂｒｏｄｙ－Ｍｏｓｈｉｎｓｋｙ系数的计算。为了提高计算效率，作者将这些系数用超几何函数予以表示，并由此提出了高效的计算方法。     

An elliptical basis function network for classification of remote sensing images

An elliptical basis function (EBF) network is employed in this study for the classification of remotely sensed images. Though similar in structure, the EBF network differs from the well-known radial basis function (RBF) network by incorporating full covariance matrices and employing the expectation-maximization (EM) algorithm to estimate the basis functions. Since remotely sensed data often take on mixture-density distributions in the feature space, the network not only possesses the advantage of the RBF mechanism, but also utilizes the EM algorithm to compute the maximum likelihood estimates of the mean vectors and covariance matrices of a Gaussian mixture distribution in the training phase. Experimental results show that the EM-based EBF network is more effective in training and simpler in structure than an RBF network constructed for the same task.The research was supported by grant 40101021 from the Natural Science Foundation of China, and grant 2002AA135230 from Hi-Tech research and development program of China. The authors would like to thank the reviewers for their valuable comments.     

Blending-Based Stochastic Simulator

This paper presents a new method of constructing random functions whose realizations can be evaluated efficiently. The basic idea is to blend, both stochastically and linearly, a limited set of independent initial realizations previously generated by any chosen simulation method. The blending stochastic coefficients are determined in such a way that the new random function so generated has the same mean and covariance functions as the random function used for generating the initial realizations.     

Point Sources of Elastic Deformation: Elementary Sources, Static Displacements

Summary. We construct a catalogue of all possible elementary point sources for static deformation in an elastic solid. the familiar double-couples, CLVD's, centres of compression and dilatation, etc., are all members of the complete catalogue. the sources are classified according to the rank of the seismic moment tensor, and according to the weight (or order) of the irreducible tensor representation of the 3-D rotation group. These sources can be classified as belonging to one of three general classes. the static excitation functions are calculated for an infinite, homogeneous, isotropic medium for all these sources. We show that, except for sources belonging to these three general classes, all other sources — which are numerous for the tensors of high rank — are null static sources. That is, sources that do not produce any static displacement outside of the source region. Due to the presence of null sources, an inversion of the static deformation data is non-unique. the expansion of the equivalent-force tensors and the stress glut tensors (or seismic moment tensors) into a set of the symmetric trace-free source tensors is proposed. the sources corresponding to seismic moment tensors of the second, third and fourth ranks are considered in more detail. We identify the third-rank sources with rotational dislocations or disclinations.