Stroke构造、移位一体化的道路网示意化方法

针对当前常用示意性地图生成方法中往往存在简化程度不够、未考虑长度信息以及时间效率低等问题,提出了stroke构造、移位一体化的道路网示意化方法。该方法在对道路网进行stroke划分构造的同时,直接对其进行渐进式的移位示意化及拓扑检查。同时本文还提出利用相似分形维数来定量比较并验证不同示意化方法的有效性。试验表明,本算法在考虑原始线形的基础上,示意化过程简单,时间效率较高,减少了拓扑冲突问题,保证了拓扑一致性及路网的均衡分布,具有较好的清晰度和认知度。     

基于规则格网的DEM插值实验

通过地形建模,将6个地形隶属函数按照101×101,112×112,126×126,144×144,168×168,201×201,257×257 7种格网密度生成6种局部地形单元的规则格网DEM;使用反距离加权(IDW)等7种插值算法,将前6种格网密度下的DEM插值成257×257规格;从原始257×257DEM中随机抽取检查点计算残差,并对残差中误差进行分析。通过分组插值实验,运用控制变量法、方差分析等方法研究地貌类型、采样密度和插值算法对DEM插值精度的影响。     

平行轮廓线自适应插值的三维表面重构算法

针对稀疏点轮廓线在三维表面重构中存在的不够光滑、匹配易出错和难以交互等问题,提出了基于轮廓线插值的三维表面重构方法,并通过实验进行验证。文章分析了三维表面重构的形式化表达方法和准则,提出了基于虚跨距解决稀疏点轮廓线三维重构的基本思想;基于同步前进算法和虚跨距,设计了自适应等比例插值的同步前进三角面镶嵌算法;针对轮廓线插值违背三角表面重构拓扑准则的规律进行了分析,提出了相应的修正算法,保证了三角网的拓扑一致性。试验结果表明:本算法所构建的三角网表面更为光滑,各个三角面的大小更为均匀,具有更好的鲁棒性和交互性,对稀疏和稠密坐标点的轮廓线均能构建出更为理想的三维表面模型。     

抗差估计的多波束测深数据内插方法

针对传统多波束测深数据网格化内插方法抗差性不足的问题,该文将抗差估计理论应用于网格化内插方法中。基于不确定度指标建立了距离不确定度联合加权内插模型,在该模型的权函数中引入了等价权,同时给出了等价权设计方案及迭代初值确定原则,以增强模型的抗差性;介绍了一种节点邻域参考点合理选取方法,给出了改进内插法实现的基本步骤,最后通过实测数据对反距离加权法、联合加权法及改进方法内插水深质量进行比较。结果表明,改进方法内插结果要优于另外两种方法,可应用于多波束测深数据网格化处理中。     

华北克拉通基底花岗质片麻岩变形和流变学研究——以辽西寺儿堡地区为例

辽西寺儿堡镇新太古代花岗质片麻岩内发育的宏观、微观构造变形特征表明该地区曾遭受了强烈的韧性变形改造。花岗质岩石变形程度在初糜棱岩–糜棱岩之间,岩石经历了SWW向左行剪切作用改造。岩石中石英有限应变测量判别结果表明,构造岩类型为L-S型,为平面应变。岩石的剪应变平均值为1.43,运动学涡度值为0.788~0.829,指示岩石形成于以简单剪切为主的一般剪切变形中。此外,石英颗粒以亚颗粒旋转重结晶和颗粒边界迁移重结晶作用为主,长石颗粒塑性拉长,部分发生膨凸式重结晶作用;石英组构特征(EBSD)揭示石英以中–高温柱面滑移为主;石英颗粒边界具有明显的分形特征,分形维数值为1.151~1.201,指示了中高温变形条件。综合石英、长石的变形行为、石英组构特征以及分形法Kruhl温度计的判别结果,推断辽西寺儿堡镇新太古代花岗质片麻岩经历过480~600℃的中高温变形,其同构造变质相为高绿片岩相-低角闪岩相。花岗质岩石的古差异应力为10.62~12.21 MPa,估算的应变速率为10~(–11.67)~10~(–13.34) s~(–1),即缓慢的变形,可能记录早期中高温、低应变速率的韧性变形过程,反映华北克拉通基底中下部地壳变形特征。     

分形分维理论在地质灾害发育及空间分布规律中的应用——以长安区滑坡、崩塌地质灾害为例

长安区区内地层岩性、地质构造多样,加之人类经济工程活动、降雨、地震等诱发影响,秦岭北坡及黄土台塬边缘滑坡崩塌地质灾害发育,是陕西省地质灾害较严重的县(区)之一。本文以长安区为研究区域,在对长安区自然地理及区域地质环境条件等资料充分收集的基础上,结合长安区滑坡崩塌灾害点的实际调查资料,详细分析了长安区的滑坡崩塌灾害发育特征及空间分布规律,并采用网络覆盖法,对长安区滑坡崩塌地质灾害的空间分布进行了分形分维计算。求得地质灾害的空间分布分维值为0.876 6,反应出长安区滑坡崩塌地质灾害空间分布特性的复杂性。研究结果表明,分形分维理论能合理评价地质灾害空间分布特征,对类似地区的地质灾害空间分布特征分析具有借鉴意义。     

Dasymetric Spatiotemporal Interpolation

This research applies the principles of dasymetric mapping to spatiotemporal interpolation by extending the spatial concepts of zone and area to their temporal analogs of interval and duration, respectively. An example application of dasymetric spatiotemporal interpolation using crime event data is presented. Results indicate that dasymetric spatiotemporal interpolation significantly improves the accuracy of estimates over areal or duration weighting. In addition, even when dasymetric interpolation in either the spatial or temporal dimension is relatively weak, combining dasymetric estimation in both space and time dimensions simultaneously has the potential to amplify the accuracy of the overall dasymetric estimation.     

MORPHOMETRIC ANALYSIS OF KARST DEPRESSIONS ON A MEDITERRANEAN KARST MASSIF

Dolines are the most typical landform in karst landscapes. They play an important role in karst hydrogeology; they are effective sediment traps that provide information about past climates and karst evolution, and can also be indicators of tectonic activity. The morphometic analysis of those depressions can help in their study and interpretation. However, detailed identification and delineation of dolines becomes increasingly difficult when the study area is very large, of limited access or covered by vegetation. Here, we have used a method developed by the authors for automatic identification and delineation of karst depressions based on the use of a pit removal algorithm. Next a complete morphometric analysis of the karst depressions has been carried out. The main advantage of the methodology is that by delineating and identifying the depth of each doline from its rim, many morphometric statistics can be obtained, including the relation of doline depth and area with respect to altitude, the main doline orientations, doline density, the fractal dimension of doline area and frequency, and point field analysis, among others. The results are interpreted from a geomorphological point of view and the mapped dolines can be easily integrated into geomorphological mapping. A case study is presented from the large Sierra Gorda karst massif in southern Spain, where 3100 depressions have been identified. A strong structural control of doline shape and orientation is demonstrated; the main areas of doline density are highlighted and doline size is shown to follow a fractal law.     

Hybrid N-order Lagrangian Interpolation Eulerian?Lagrangian Method for Salinity Calculation

The Eulerian?Lagrangian method (ELM) has been used by many ocean models as the solution of the advection equation, but the numerical error caused by interpolation imposes restriction on its accuracy. In the present study, hybrid N-order Lagrangian interpolation ELM (LiELM) is put forward in which the N-order Lagrangian interpolation is used at first, then the lower order Lagrangian interpolation is applied in the points where the interpolation results are abnormally higher or lower. The calculation results of a step-shaped salinity advection model are analyzed, which show that higher order (N=3?8) LiELM can reduce the mean numerical error of salinity calculation, but the numerical oscillation error is still significant. Even number order LiELM makes larger numerical oscillation error than its adjacent odd number order LiELM. Hybrid N-order LiELM can remove numerical oscillation, and it significantly reduces the mean numerical error when N is even and the current is in fixed direction, while it makes less effect on mean numerical error when N is odd or the current direction changes periodically. Hybrid odd number order LiELM makes less mean numerical error than its adjacent even number order LiELM when the current is in the fixed direction, while the mean numerical error decreases as N increases when the current direction changes periodically, so odd number of N may be better for application. Among various types of Hybrid N-order LiELM, the scheme reducing N-order directly to 1st-order may be the optimal for synthetic selection of accuracy and computational efficiency.     

基于非均质地层参数表征的颗粒流模型构建

非均质地层参数的获取与表征方法是浅层地质体改造的重要依据。基于非均质地层的地质建模,表征具有随机不确定性和模糊不确定性的土体表观参数和物理力学参数。利用 Weibull 分布统计描述局部地层信息,并应用Diamond-Square分形插值方法,结合差分盒子维计算非均质地层分形维数,进一步合理演算整体地层参数。采用PFC3D数值模拟软件,结合马尔科夫链蒙特卡洛方法采样数据,对非均质地层块石颗粒及土颗粒进行建模。通过非均质地层模拟结果,发现其破坏过程与应力应变曲线要比均质土体复杂得多,由于块石大颗粒的含量、分布特性对土样力学性质影响显著。本文建模方法有助于非均质地层参数表征研究。