基于正交贴体网格的VTI介质地震波数值模拟

为了解决规则网格有限差分法在处理起伏地表自由边界条件时需进行复杂坐标转换和插值运算的问题,以及对起伏地表进行阶梯状离散近似所产生的虚假散射波问题,将计算流体力学中的正交贴体网格生成技术引入到起伏地表下VTI介质的网格剖分中,采用基于同位网格的DRP/opt MacCormack有限差分法模拟起伏地表下VTI介质中的地震波场。网格的正交性使实施自由边界条件时无需做复杂的坐标转换和插值运算。数值算例表明,本文提出的网格剖分法能够有效消除阶梯状网格产生的虚假散射波,从而提高起伏地表下VTI介质地震波数值模拟精度。     

基于Landsat TM影像不同地表温度反演算法的比较

地表温度是反映地表环境的一个重要参数,精确获取地表温度的方法对研究城市热岛效应、进行生态环境监测评价是必不可少的。利用遥感手段进行地表温度的反演是一种较新颖的方法,相对成熟的反演算法主要有单通道算法、辐射传输方程法及单窗算法。本文以山东省威海市Landsat TM数据为例,分别利用辐射传输方程法和单窗算法进行地表温度反演,通过遥感目视解译的方法提取研究区各类典型地物对应的地表温度,进行统计分析,结果表明,两者反演得到的地表温度具有一致的变化趋势,其中单窗算法中不同典型地物的地表温度整体标准差较小,如建设用地为0.530,起伏波动小,算法精度略优于辐射传输方程法。     

有限元法与伪谱法混合求解弹性波动方程

在地震波场数值模拟中,有限差分法、有限元法和伪谱法都是常用的基本方法,但它们各有不同的适应性和优缺点,如有限差分法、有限元法都存在减弱网格频散和提高计算效率的矛盾,而伪谱法的网格频散小且计算效率高.有限差分法和伪谱法在处理地表结构复杂或地表剧烈起伏以及地下结构复杂的情况时存在较大的难度,而有限元法可较为理想地拟合起伏地表和任意弯曲界面,且可方便地处理自由边界条件和界面边界条件.尝试将有限元法和伪谱法相结合,形成地震波场数值模拟的一种混合方法,利用二者的优点,克服二者的缺点,达到既减弱网格频散又提高计算精度和效率的目的.并采用所谓的‘过度区域‘技术解决两种不同算法的衔接问题.模拟实例表明,给出的混合模拟方法不失为弹性波场数值模拟的一种有效方法.     

基于卷积神经网络的遥感影像地表覆盖分类

遥感影像地表覆盖分类是地理国情监测和地理信息资源建设中至关重要的环节,利用卷积神经网络对遥感影像进行特征提取和分类,具有十分重要的科研和应用价值。为提高遥感影像的地表覆盖分类精度,在深度卷积神经网络VGGNet的基础上,采用SeLU函数作为激活函数,并将激活函数中的λ、α作为训练参数,得到改进的VGGNet,用逐层贪婪算法对网络参数初始化,并选择适当的学习次数利用迁移学习的方法对网络参数调整,以提高网络的泛化能力来提取遥感影像各类别的深层特征,从而有效进行地表覆盖分类。通过GF-1卫星影像的实验表明本文方法在地表覆盖分类精度方面的优越性。     

基于支持向量机理论的地下水动态遥感监测模型与应用

地下水是我国内陆干旱地区水资源的重要组成部分,也是极为敏感的生态环境因素之一。地下水动态变化影响着绿洲和湿地的演化,以及土地资源的开发。西北地区地下水监测网尚不完善,动态资料相对缺乏。遥感技术可以弥补传统地下水位监测手段的不足。由于降水极少,西北干旱区地表反射率与地下水水位埋深关系极其密切。选用归一化植被指数(NDVI)、地表温度(LST)数据,应用支持向量机回归方法,建立西北干旱地区地下水位遥感监测模型。提取MODIS影像中的NDVI和LST产品上的地表温度和植被指数信息,作为模型的输入,通过合理选择核函数进行支持向量机的回归分析,从而建立地表植被指数、地表温度与地下水位的相关数学模型,并分析了不同核函数所拟合结果。在河西走廊疏勒河流域的研究成果表明,运用MODIS数据开发地下水动态模型反演水位变化是可行的,模型拟合的结果比较符合实际情况,尤其是对于细土平原地下水浅埋地区模型应用效果更为理想。一次多项式核函数适合模拟埋深小于3m浅埋地下水,径向基函数RBF核函数和三次多项式核函数法则更适合模拟较大埋深情况。开发的地下水位遥感监测模型可为西北干旱区水循环研究和流域水资源管理提供技术手段。     

遥感影像空间格局变异函数分析研究进展

随着多光谱传感器的广泛运用,利用地物光谱响应特征提取地表信息的技术日益成熟,但是由于地表状况的复杂性和光谱响应的局限性,光谱方法在指示平均大小、空间异向性、空间分布、空间异质性等格局信息方面存在不足,因此挖掘遥感影像的空间格局特征日益受到研究者的重视。已有研究发现,变异参数与地表场景参数存在一定的对应关系,通过变异参数可以实现地表场景参数的提取,因此变异函数分析方法被广泛应用于真实遥感影像格局分析中,具体包括平均尺度提取、周期性格局探测、空间异质性表征与空间异向性描述等地表格局参数量化方面、最佳尺度选择与影像纹理分析等遥感影像信息提取方面。尽管变异函数分析方法在上述应用领域中都发挥了重要作用,但是当前利用变异函数进行遥感影像空间格局分析大多局限于定性描述层面,缺乏精确化的量化描述与分析,限制了变异函数分析方法应用的进一步拓展,究其原因在于对遥感影像格局变异函数分析的内在机制缺乏深入了解。本文回顾了近20年来变异函数分析方法在遥感格局分析领域的主要应用,并对该方法本身的优势和存在的不足进行了总结,可为变异函数这一工具在遥感影像格局分析方面的有效应用提供参考。     

多极化SAR数据反演额济纳冲积扇地表参数

本文采用遗传算法和后向散射模型相结合的方法,探讨了多极化SAR(Synthetic Aperture Radar)数据反演冲积扇地表参数(地表粗糙度和土壤湿度)的可行性.通过理论模拟和实地地表测量对比,表明该方法在反演冲积扇地表参数方面是切实可行的.该方法可以根据获取数据情况不同,灵活调整反演的代价函数式,并且用于反演的SAR图像必须大于(或等于)两景,数据越多,反演结果越精确.在此基础上,利用ENVISAT ASAR数据和ALOS PALSAR数据,对内蒙古额济纳冲积扇的地表参数进行了反演计算.结果表明,额济纳冲积扇地面相对比较平坦,粗糙度参数变化较小,大部分地区均方根高度小于1.0cm.黑河沿岸地区粗糙度参数较大,而远离黑河的戈壁滩地表粗糙度较小.该地区的土壤湿度反演结果表明,该地区属于极端干旱区域,大部分地区土壤体积含水量低于10%.     

基于BEMD分解的地貌分类研究

地貌是指地势高低起伏的变化,即地表的形态。地貌划分对气温、降水、太阳辐照等诸多应用领域都有着重要作用。本文选择空间分辨率为90 m的福建省的数字高程模型(DEM)数据作为地理信号,运用二维经验模态分解(BEMD)进行分解处理,得到多个具有不同尺度、不同物理意义的本征模函数(BIMF1-BIMF3)以及对应余量ORIG。这些BIMF分量对应不同尺度的微观地形,ORIG余量表现为该研究区的地貌分布趋势,体现了平原、丘陵与山地的大致分布区域。运用变点分析法确定最佳计算单元,利用地形起伏度对各个信号区域进行一级分类,依据绝对高度进行二级分类,最后将一级分类与二级分类相结合,实现对地形的分类,这一分类过程体现了研究地区地貌组合复杂的特征。结果表明:①叠加BIMF分量,提取出分量和大于74 m的区域为高频信号区域。该区域以小起伏度的低山为主,并伴随有丘陵和小起伏度中山。②ORIG余量中余量高度小于等于340 m的区域,去除其中包含的高频信号区域后为低频信号区域,该区域以平原、丘陵为主。③剩余区域定义为中频信号区域,该区域的地貌以平丘陵和小起伏的山地为主。研究成果表明福建地貌可分为7种主要类型:低频平原,低频丘陵,中频丘陵,高频丘陵,中频小起伏低山,高频小起伏低山,高频小起伏中山。     

面向无测风区域的复杂地形风场模拟研究

由于气象站点的布设受经济、地形、技术等因素限制,使得某风能资源丰富的无测风区域缺少气象站点数据,无法进行风场模拟。针对这一问题,尝试用NCEP再分析风速数据代替传统的气象站点数据,根据NCEP再分析风速数据的栅格数据特点,提出了一种面向无测风区域的复杂地形风场模拟方法。该风场模拟方法以SRTM DEM数据和NCEP再分析风速数据为数据源,首先通过均值变点分析法获取反映无测风区域地形起伏特点的最佳统计单元,得到无测风区域复杂地形的地形起伏度;然后对Cressman插值方法进行改进,将地形起伏度引入到权重函数中,并根据交叉检验中均方根误差最小原则求解权重函数中参数的最优解;最后将最优解代入权重函数中对风场进行加密。实验结果表明:风场模拟结果与气象站点风速之间的相关系数为0.778 2,模拟风速的变化趋势与实测风速基本一致;风电场选址区主要建设在模拟风速较大的地区,十分符合风电场一般建设在风能资源丰富的地方的特点,间接验证面向无测风区域的复杂地形风场模拟方法是有效的。     

基于数字高程模型的城市地表开阔度研究——以南京老城区为例

 城市地表开阔度是对地表形态开阔程度的定量描述,间接反映了通风、日照等城市物理性能。本文基于城市数字高程模型,以自然地形较为起伏的丘陵城市——南京为例,采用以光线追踪算法为核心的栅格计算模型,对南京老城区的地表开阔度进行了计算,并以典型样区分析和街区统计为方法,综合分析了建筑与自然地形对城市地表开阔度的影响,结果表明:南京老城区公共绿地、广场以及文物古迹用地的地表开阔度相对较高,均大于0.8,社区中心和居住用地的地表开阔度相对较低,均小于0.6,突出特点是居住用地地表开阔度小于金融商业用地。地表开阔度相对较低的街区主要临近主干道,形成较为明显的西北——东南方向分布的条带状特征,整体呈现出中心条带低四周高的宏观特征。