Interpolation by regularized spline with tension: I. Theory and implementation

Bivariate and trivariate functions for interpolation from scattered data are derived. They are constructed by explicit minimization of a general smoothness functional, and they include a tension parameter that controls the character of the interpolation function (e.g., for bivariate case the surface can be tuned from a membrane to a thin steel plate), Tension can be applied also in a chosen direction, for modeling of phenomena with a simple type of anisotropy. The functions have regular derivatives of all orders everywhere. This makes them suitable for analysis of surface geometry and for direct application in models where derivatives are necessary. For processing of large datasets (thousands of data points), which are now common in geosciences, a segmentation algorithm with a flexible size of overlapping neighborhood is presented. Simple examples demonstrating flexibility and accuracy of the functions are presented.On leave from the Department of Physical Geography and Cartography, Comenius University, Mlynská dolina, Bratislava, Czechoslovakia.On leave from the Institute of Physics, Dúbravská cesta 9, Bratislava, Czechoslovakia.     

利用主成分分析改善土地利用变化的遥感监测精度——以珠江三角洲城市用地扩张为例

近年来,珠江三角洲由于经济的快速发展,城市用地急剧增加,利用多时相的遥感图,可以定量地监测这种城市化的现象。但理,由一般的遥感动态监测方法所得的结果往往夸大变化的程度,以及获得一些不合理的结论．该文提出主成分分析的方法来改善遥感动态监测的精度。将该方法应用应用于珠江三角洲发展最快的东莞市,获得了较满意的结果。     

“4D”技术及其应用

介绍“４Ｄ”（即数字地面模型、数字正射影像、数字栅格地图和数字线划地图）系列产品的生产流程、关键技术、及其在耕地保护、防灾减灾救灾和城市规划等方面的应用方法。     

区域地下水信息管理系统研究

本文通过《河南省地下水动态信息管理系统》的建立,讨论了区域地下水信息管理系统的设计目标,数据结构和软件实现,着重介绍了系统的软件设计特点,最后指出了这种管理系统在水文地质部门中具有广泛的应用前景。     

A realistic model for point-sources imaged on array detectors: The model and initial results

We have constructed a computer model for simulation of point-sources imaged on two-dimensional detectors. An attempt has been made to ensure that the model produces data that mimic real data taken with 2-D detectors. To be realistic, such simulations must include randomly generated noise of the appropriate type from all sources (e.g. source, background, and detector). The model is generic and accepts input values for parameters such as pixel size, read noise, source magnitude, and sky brightness. Point-source profiles are then generated with noise and detector characteristics added via our model. The synthetic data are output as simple integrations (onedimensional), as radial slices (two-dimensional), and as intensity-contour plots (three-dimensional). Each noise source can be turned on or off so that they can be studied separately as well as in combination to yield a realistic view of an image. This paper presents the basic properties of the model and some examples of how it can be used to simulate the effects of changing image position, image scale, signal strength, noise characteristics, and data reduction procedures.Use of the model has allowed us to confirm and quantify three points: 1) The use of traditionalsize apertures for photometry of faint point-sources adds substantial noise to the measurement which can significantly degrade the quality of the observation; 2) The number of pixels used to estimate the background is important and must be considered when estimating errors; and 3) The CCD equation normally used by the astronomical community consistently overestimates the signal-to-noise obtainable by a measurement while a revised equation, discussed here, provides a better estimator.     

A combined study of macroseismic data and focal mechanisms applied to the West-Bohemian earthquake,Czechoslovakia

The main shock of the West-Bohemian earthquake swarm, Czechoslovakia, (magnitudem=4.5, depthh=10 km) exhibits an irregular areal distribution of macroseismic intensities 6° to 7° MSK-64. Four lobes of the 6° isoseismal are found and the maximum observed intensity is located at a distance of 8 km from the instrumentally determined epicentre. This distribution can be explained by the energy flux of the directS wave generated by a circular source, the hypocentral location and focal mechanism of which are taken from independent instrumental studies. The theoretical intensity, which is assumed to be logarithmically proportional to the integrated squared ground-motion velocity (i.e.,I=const+log v ² (t)dt), fits the observed intensity with an overall root-mean-square error less than 0.5°. It is important that the present intensity data can also be equally well explained by the isotropic source. The fit was attained by means of a horizontally layered model though large fault zones and an extended sedimentary basin suggest a significant lateral heterogeneity of the epicentral region. The results encourage a broader application of the simple modelling technique used.     

Structural and remote sensing studies of the southern Betsimisaraka Suture, Madagascar

An assessment of the southern Betsimisaraka Suture (B.S.) of southeastern Madagascar using remote sensing and field investigation reveals a complex deformation history. Image processing of Landsat ETM+data and JERS-I Synthetic Aperture Radar (SAR) imagery was integrated with field observations of structural geology and field petrography. The southern B.S. divides the Precambrian basement rocks of Madagascar in two parts. The western part includes Proterozoic rocks whereas the eastern part is an Archean block, named the Masora block. The southern part of the B.S. includes high-grade metamorphic rocks, recording strong deformation and has mineral deposits including chromite, nickel, and emerald, characteristic of oceanic material that is compatible with a suture zone.Large-scale structural features indicate ductile deformation including three generations of folding (F1, F2, and F3) associated with dextral shearing. The first folding event (F1) shows a succession of folds with NE striking axial planes. The second folding event (F2) mainly has north–south striking axial planes and the last event (F3) is represented by mega folds that have ENE–WSW axial plane directions and have NNW and SSE contractional strain patterns. Closure of the Mozambique Ocean between two components of Gondwana sandwiched rocks of the B.S. and formed upright folds and shortening zones which produced N–S trending lineaments. Later dextral movements followed the contraction and formed NW–SE trending lineaments and N–S trending normal faults associated with dextral strike slip faults and fractures.     

GIS技术支持下的招远市金矿区崩塌遥感调查

在总结前人研究成果的基础上,全面系统地分析了诱发招远市金矿区崩塌的自然和人为因数,然后运用遥感技术对金矿区遥感图像进行处理,提取诱发崩塌的条件因子,聘请有经验的专家对各项诱发因子进行诊断分析,以此作为判断条件,在G IS技术的支持下对诱发崩塌的条件因子进行空间分析,预测出招远金矿区发生崩塌的危害程度,为防治崩塌提供科学依据。预测结果表明,招远金矿区崩塌有进一步发展的可能,其中,高危险区有36820 m2,中易发区有50 610 m2,低易发区有67 200 m2,需要采取有效措施加以防治。     

用于湿地气候效应模拟的三江平原下垫面数据获取

三江平原是我国沼泽湿地的重要分布区。在20世纪的后50 a中,三江平原在经过几次大规模的开荒后,其下垫面发生了明显的阶段性变化,湿地大面积消失,农田面积不断增加。为了认识20世纪下半叶不同时期三江平原湿地的气候效应,拟采用第三代区域气候模式(RegCM3)来开展研究。而目前被广泛用于RegCM3的GL-CC数据,所反映的研究区下垫面几乎不存在沼泽湿地,这与事实严重不符,所以必须获取三江平原有沼泽湿地存在的真实的下垫面数据,才能得到可靠的湿地气候效应的模拟结果。分别选择了4期(1954年、1985年、1996年和2000年)有阶段代表性的航片、TM影像数据为主要数据源,利用ERDAS8.7和Arcview3.2等软件,首先得到4期研究区的土地利用数据,然后通过土地利用数据类型向模式数据类型的转换,从4期研究区土地利用数据中提取出了三江平原下垫面的模式数据。经过误差分析,认为获取的4期三江平原下垫面模式数据准确的反映出研究区几个时期的湿地下垫面实况,可以作为三江平原湿地气候效应模拟的基础数据。