Derivation of planetary topography using multi-image shape-from-shading

In many cases, the derivation of high-resolution digital terrain models (DTMs) from planetary surfaces using conventional digital image matching is a problem. The matching methods need at least one stereo pair of images with sufficient texture. However, many space missions provide only a few stereo images and planetary surfaces often possess insufficient texture.This paper describes a method for the generation of high-resolution DTMs from planetary surfaces, which has the potential to overcome the described problem. The suggested method, developed by our group, is based on shape-from-shading using an arbitrary number of digital optical images, and is termed “multi-image shape-from-shading” (MI-SFS). The paper contains an explanation of the theory of MI-SFS, followed by a presentation of current results, which were obtained using images from NASA's lunar mission Clementine, and constitute the first practical application with our method using extraterrestrial imagery. The lunar surface is reconstructed under the assumption of different kinds of reflectance models (e.g. Lommel-Seeliger and Lambert). The represented results show that the derivation of a high-resolution DTM of real digital planetary images by means of MI-SFS is feasible.     

Homoclinic and heteroclinic orbits in the photogravitational restricted three-body problem

We study numerically the asymptotic homoclinic and heteroclinic orbits around the hyperbolic Lyapunov periodic orbits which emanate from Euler's critical points L ₁ and L ₂, in the photogravitational restricted plane circular three-body problem. The invariant stable-unstable manifolds associated to these Lyapunov orbits, are also presented. Poincaré surface of sections of these manifolds on appropriate planes and several homoclinic and heteroclinic orbits for the gravitational case as well as for varying radiation factor q ₁, are displayed. Homoclinic-homoclinic and homoclinic-heteroclinic-homoclinic chains which link the interior with the exterior Hill's regions, are illustrated. We adopt the Sun-Jupiter system and assume that only the larger primary radiates. It is found that for small deviations of its value from the gravitational case (q ₁ = 1), the radiation pressure exerts a significant impact on the Hill's regions and on these asymptotic orbits.     

论地理科学体系建设

本文在论述重建地理科学体系必要性的基础上,提出和重点阐述一个新的地理科学体系的框架和内容,并由此分析我国未来地理科学学科建设的战略重点。     

A simple method using climatic variables to estimate canopy temperature,sensible and latent heat fluxes in a winter wheat field on the North China Plain

Estimation of evapotranspiration from a crop field is of great importance for detecting crop water status and proper irrigation scheduling. The Penman–Monteith equation is widely viewed as the best method to estimate evapotranspiration but it requires canopy resistance, which is very difficult to determine in practice. This paper presents a simple method simplified from the Penman–Monteith equation for estimating canopy temperature (T_c). The proposed method is a biophysically‐sound extended version of that proposed by Todorovic. The estimated canopy temperature is used to calculate sensible heat flux, and then latent heat flux is calculated as the residual of the surface energy balance. An eddy covariance (EC) system and an infrared thermometer (IRT) were installed in an irrigated winter wheat field on the North China Plain in 2004 and 2005, to measure T_c, and sensible and latent heat fluxes were used to test the modified Todorovic model (MTD). The results indicate that the original Todorovic model (TD) severely underestimates T_c and sensible heat flux, and hence severely overestimates the latent heat flux. However, the MTD model has good capability for estimating T_c, and gives acceptable results for latent heat flux at both half‐hourly and daily scales. The MTD model results also agreed well with the evapotranspiration calculated from the measured T_c. Copyright © 2008 John Wiley & Sons, Ltd.     

基于二次曲面的Shepard拟合法在GPS水准中的应用

针对GPS水准中常用拟合方法的不足,给出一种基于二次曲面的Shepard拟合法。通过某测区GPS水准实测数据,采用二次曲面拟合法、Shepard拟合法、多面函数法和基于二次曲面的Shepard拟合法对测区高程异常进行拟合,结果表明,基于二次曲面的Shepard拟合法计算得到的高程异常的精度和可靠性最高。     

无人机遥感真正射影像高精度制图

近年来无人机低空遥感技术不断发展,利用无人机影像生成的真正射影像（TDOM）在成图精度、制作流程方面仍有提高的空间。本文采用固定翼无人机和专业摄影相机采集影像,布设地面控制点,提出了利用运动恢复结构（SfM）和多视立体视觉（MVS）工作流来生成高精度数字表面模型（DSM）和数字正射影像（DOM）的方法;对遮蔽倾斜部分进行阴影检测、DSM修编和多视影像纹理补偿生成TDOM;最后用TDOM上随机分布的检查点进行精度检查,水平精度为3.3 cm,垂直精度为7.5 cm;消除了DOM中倾斜和阴影部分,使建筑物保持垂直视角,生成的满足1：500比例尺高精度并消除倾斜阴影的TDOM可用于农村宅基地确权、国土规划设计等领域。     

车载激光扫描数据分类支持下的路面数据提取

车载激光扫描系统可以快速采集道路及两旁的建筑物、植被、电杆等地物的点云数据,而点云数据的分类提取是车载激光扫描系统应用的关键。本文选用全景激光移动测量系统获取的激光点云数据,分析了路面点云数据的特征,采用渐进格网法进行了路面点云数据的提取研究;通过试验区的实例验证,取得了较好的分类效果。     

北京地区Landsat 8 OLI高空间分辨率气溶胶光学厚度反演

卫星气溶胶光学厚度(AOD)反演中,传统暗目标方法在反射率较低的水体、浓密植被覆盖区域取得了较好效果,在反射率较高且结构复杂的高反射地表上空目前多采用深蓝算法,但存在空间分辨率较低,对细节分布描述性较差等问题。为解决这一问题,本文首先以5年(2008年—2012年)长时间序列MODIS地表反射率产品为基础,采用最小值合成法建立500 m分辨率逐月地表反射率产品数据集,然后利用地物波谱库中典型地物波谱数据,分析建立MODIS与Landsat 8 OLI传感器蓝光波段反射率转换模型,最后北京地区AERONET地基观测数据确定了气溶胶光学物理参数,并反演获取了北京地区上空500 m分辨率的AOD分布。为验证反演算法的精度,分别将反演结果同AERONET及MODIS/Terra气溶胶产品(MOD04)进行交叉对比,同时利用相关系数R,均方根误差RMSE,平均绝对误差MAE以及MODIS AOD产品预期误差EE共4个指标进行衡量。结果表明:算法反演获取的AOD与AERONET观测值具有较高的一致性,各指标分别为R=0.963,RMSE=0.156,MAE=0.097,EE=85.3%,稍优于MOD04产品(R=0.962,RMSE=0.158,MAE=0.101,EE=75.8%),并且有效的对比点数也高于MOD04。通过与地基观测相比,卫星遥感获取的高分辨率城市地区AOD精度可作为定量评估城市空气质量的有效依据。     

渤海湾周缘高温地热异常特征及深部热结构分析

渤海湾周缘高温地热异常区是华北重要的地热资源区,研究其浅层高温地热异常及深层高温地热结构,是解决中国东部高温地热勘查开发瓶颈的重要基础,对加快构建东部现代能源体系、推进“双碳”目标进程具有重要的科学意义。利用最新的中国大陆岩石圈速度模型USTClitho2.0和航磁资料,在大地热流资料约束下,本文计算了渤海湾周缘浅层与深层高温热结构,结果表明:该区居里面平均深度20.8km,居里温度420℃,中、上地壳高温热结构与居里面起伏形态密切相关。上地壳生热率随深度增大而递减,反映了放射性生热元素向上迁移、在浅层富集的垂向分异特征。Moho面温度在600～800℃之间,岩石圈底界面1300℃等温面深度在66.3～97.5 km之间,平均值为76.8 km。依据USTClitho2.0模型Vp、Vs速度结构,本文修正了前人的地壳分层结构及上地壳生热率,计算了研究区内全部91个大地热流点的壳幔热流比,结果表明:Qc/Q=56.4%、Qm/Q=43.6%, Qc/Q>Qm/Q,显示该区为“热壳”结构。前人认为此区是“冷壳热幔”型热结构,本文的计算结果修正了前人的认识,指出上地壳层热流在地表总热流...     

近地表隐伏断层地震波场响应特征与初至波成像应用研究

近地表隐伏断层直接影响工程安全施工。通过构建浅覆盖层下方陡立界面断层数值模型,分析近地表隐伏断层地震波场的传播与转换规律。通过工程实测数据验证了数值模拟波场响应特征的可靠性,并利用初至波走时层析成像技术对地下介质波速进行成像。研究结果表明：(1)初至波转换形成的断点绕射波具有1频率升高、速度突变的特点,可清晰指示测区内断裂带的存在,可作为现场数据采集时断层的直观反映;(2)工程实测中通过原始波形到时、频率信息可初步判断断层的存在,基于射线追踪的初至波层析成像结果与后期钻孔揭露情况基本吻合,说明初至波层析成像在浅层地震勘探中的实用性;(3)接收测线排列内布设震源的方式,两侧震源延拓距离短,无需额外为震源激发提供场地,在工程施工中具有经济、快速、高效的特点,可作为野外复杂地表条件地下浅层结构探测的有效观测系统使用。