Parallel computing solutions for Markov chain spatial sequential simulation of categorical fields     

Weixing Zhang  Weidong Li  Tian Zhao 《International Journal of Digital Earth》2019,12(5):566-582

The Markov chain random field (MCRF) model is a spatial statistical approach for modeling categorical spatial variables in multiple dimensions. However, this approach tends to be computationally costly when dealing with large data sets because of its sequential simulation processes. Therefore, improving its computational efficiency is necessary in order to run this model on larger sizes of spatial data. In this study, we suggested four parallel computing solutions by using both central processing unit (CPU) and graphics processing unit (GPU) for executing the sequential simulation algorithm of the MCRF model, and compared them with the nonparallel computing solution on computation time spent for a land cover post-classification. The four parallel computing solutions are: (1) multicore processor parallel computing (MP), (2) parallel computing by GPU-accelerated nearest neighbor searching (GNNS), (3) MP with GPU-accelerated nearest neighbor searching (MP-GNNS), and (4) parallel computing by GPU-accelerated approximation and GPU-accelerated nearest neighbor searching (GA-GNNS). Experimental results indicated that all of the four parallel computing solutions are at least 1.8× faster than the nonparallel solution. Particularly, the GA-GNNS solution with 512 threads per block is around 83× faster than the nonparallel solution when conducting a land cover post-classification with a remotely sensed image of 1000?×?1000 pixels.  相似文献   

Algorithms for Stellar Perturbation Computations on Oort Cloud Comets     

Hans Rickman  Marc Fouchard  Giovanni B. Valsecchi  Christiane Froeschlé 《Earth, Moon, and Planets》2005,97(3-4):411-434

We investigate different approximate methods of computing the perturbations on the orbits of Oort cloud comets caused by passing stars, by checking them against an accurate numerical integration using Everhart’s RA15 code. The scenario under study is the one relevant for long-term simulations of the cloud’s response to a predefined set of stellar passages. Our sample of stellar encounters simulates those experienced by the Solar System currently, but extrapolated over a time of 10¹⁰ years. We measure the errors of perihelion distance perturbations for high-eccentricity orbits introduced by several estimators – including the classical impulse approximation and Dybczyński’s (1994, Celest. Mech. Dynam. Astron. 58, 1330–1338) method – and we study how they depend on the encounter parameters (approach distance and relative velocity). We introduce a sequential variant of Dybczyński’s approach, cutting the encounter into several steps whereby the heliocentric motion of the comet is taken into account. For the scenario at hand this is found to offer an efficient means to obtain accurate results for practically any domain of the parameter space.  相似文献   

On the local multiscale determination of the Earth’s disturbing potential from discrete deflections of the vertical     

T. Fehlinger  W. Freeden  C. Mayer  M. Schreiner 《Computational Geosciences》2008,12(4):473-490

As a first approximation, the Earth is a sphere; as a second approximation, it may be considered an ellipsoid of revolution. The deviations of the actual Earth’s gravity field from the ellipsoidal “normal” field are so small that they can be understood to be linear. The splitting of the Earth’s gravity field into a “normal” and a remaining small “disturbing” field considerably simplifies the problem of its determination. Under the assumption of an ellipsoidal Earth model, high observational accuracy is achievable only if the deviation (deflection of the vertical) of the physical plumb line, to which measurements refer, from the ellipsoidal normal is not ignored. Hence, the determination of the disturbing potential from known deflections of the vertical is a central problem of physical geodesy. In this paper, we propose a new, well-promising method for modelling the disturbing potential locally from the deflections of the vertical. Essential tools are integral formulae on the sphere based on Green’s function with respect to the Beltrami operator. The determination of the disturbing potential from deflections of the vertical is formulated as a multiscale procedure involving scale-dependent regularized versions of the surface gradient of the Green function. The modelling process is based on a multiscale framework by use of locally supported surface curl-free vector wavelets.   相似文献   

Using global node-based velocity in random walk particle tracking in variably saturated porous media: application to contaminant leaching from road constructions   总被引：2，自引：0，他引：2  

Chan-Hee Park  Christof Beyer  Sebastian Bauer  Olaf Kolditz 《Environmental Geology》2008,55(8):1755-1766

Precise and efficient numerical simulation of transport processes in subsurface systems is a prerequisite for many site investigation or remediation studies. Random walk particle tracking (RWPT) methods have been introduced in the past to overcome numerical difficulties when simulating propagation processes in porous media such as advection-dominated mass transport. Crucial for the precision of RWPT methods is the accuracy of the numerically calculated ground water velocity field. In this paper, a global node-based method for velocity calculation is used, which was originally proposed by Yeh (Water Resour Res 7:1216–1225, 1981). This method is improved in three ways: (1) extension to unstructured grids, (2) significant enhancement of computational efficiency, and (3) extension to saturated (groundwater) as well as unsaturated systems (soil water). The novel RWPT method is tested with numerical benchmark examples from the literature and used in two field scale applications of contaminant transport in saturated and unsaturated ground water. To evaluate advective transport of the model, the accuracy of the velocity field is demonstrated by comparing several published results of particle pathlines or streamlines. Given the chosen test problem, the global node-based velocity estimation is found to be as accurate as the CK method (Cordes and Kinzelbach in Water Resour Res 28(11):2903–2911, 1992) but less accurate than the mixed or mixed-hybrid finite element methods for flow in highly heterogeneous media. To evaluate advective–diffusive transport, a transport problem studied by Hassan and Mohamed (J Hydrol 275(3–4):242–260, 2003) is investigated here and evaluated using different numbers of particles. The results indicate that the number of particles required for the given problem is decreased using the proposed method by about two orders of magnitude without losing accuracy of the concentration contours as compared to the published numbers.  相似文献   

点状地名信息的加权泰森多边形检索法   总被引：1，自引：0，他引：1  

张宇  王琦  吴文周  苏奋振 《测绘学报》2017,46(11):1919-1926

鉴于地名库中多数地名仅以中心点坐标形式记录其空间位置信息,缺乏其空间范围的具体描述,造成地名检索中的片面性和局限性,本文在深入分析地名及其空间属性、空间关系的基础上,考虑同类型地名的不同性质在检索中的重要作用,利用泰森多边形在地名边界近似中的优势,提出了针对点状地名信息的加权泰森多边形检索法,依据同类型地名的面积属性作为权重指标构建泰森多边形,以近似表达地名的空间范围,进而描述地名间的各类空间关系;并从拓扑关系、方位关系与距离关系3个方面给出了具体公式以计算地名间的空间相似性;最后以行政区划边界近似与检索为例验证该方法。试验结果表明,该方法可较好地近似表达地名空间范围及空间关系,增强了检索词与地理信息资源在空间位置上的相似性度量,检索结果优于传统方法。  相似文献   

A New Parameterization of Canopy Radiative Transfer for Land Surface Radiation Models     

Feng ZHANG  Yadong LEI  Jia-Ren YAN  Jian-Qi ZHAO  Jiangnan LI  Qiudan DAI 《大气科学进展》2017,34(5):613-622

A new parameterization of canopy asymmetry factor on phase function,which is dependent on the leaf normal distribution and leaf reflection/transmission,is derived. This new parameterization is much more accurate than the existing scheme. In addition,the new solutions for both the diffuse and direct radiation can be obtained using the Eddington approximation. It is found that the direct radiation can be described as a function of the diffuse radiation. This new approach offers a substantial improvement in accuracy,as compared with the hemispheric constant method,for both isotropic and anisotropic cases. Given the analytical nature of the solution and its high accuracy,we recommend the new parameterization for application in land surface radiation modeling.  相似文献   

基于拟线性Born近似的叠前深度偏移方法   总被引：14，自引：7，他引：14       下载免费PDF全文

陈生昌  曹景忠  马在田 《地球物理学报》2001,44(5):704-710

在频率波数域和频率空间域实现了一种基于拟线性Born近似的叠前深度偏移方法，并在2-D空间进行了Marmousi模型炮集数据的处理.通过与Split-StepFourier、Phase-Screen和稳定的Born近似叠前深度偏移等方法比较，认为基于拟线性Born近似的叠前深度偏移方法不仅在效果上要优于前三者，而且还能更好地处理速度横向变化.在散射波场计算中，使用了一个更稳定的散射波场计算公式，扩大了拟线性Born近似的应用范围，使基于拟线性Born近似的叠前深度偏移方法能够适应更强的横向速度变化.  相似文献   

线性和非线性地形罗斯贝波   总被引：1，自引：0，他引：1  

赵平  孙淑清 《大气科学》1991,15(2):33-40

本文在半地转概念下,讨论了线性和非线性地形罗斯贝波的稳定性及其解。指出:线性和非线性稳定性判据形式一致。在线性时,地形东西向及南北向坡度对稳定波动的周期和传播速度有明显影响,其解为周期函数;在不稳定时,解为非周期函数。二级近似时其解为孤立波形式——不稳定时,在地形不同位置可形成东移或西移的孤立波槽或孤立波脊;而在稳定时只形成孤立波槽。三级近似时解出现间断点。  相似文献   

太平洋大尺度环流数值模拟 I:数学模式及其性能   总被引：2，自引：4，他引：2       下载免费PDF全文

张荣华  曾庆存  张学洪 《大气科学》1991,15(4):1-16

本工作基于曾庆存的海洋环流模式设计思想和大气物理研究所海洋环流模式(IAPOGCM)基本动力框架,设计并实现了四层太平洋环流模式;同时引进稳定有效的时间积分方案和海表热通量参数化方法,对太平洋大尺度环流进行了较为系统的数值模拟. 本文首先导出消除“刚盖”近似、扣除海洋标准层结、保持总有效能量守恒的OGCM基本方程组;给出空间差分格式和时间积分方案,引入海表热通量计算方法;然后,在太平洋区域对所设计的四层OGCM的一些模式性能进行了分析、表明,时间积分方案采用基于“灵活性系数”技术的加速收敛方法和适应过程中流场正压斜压进一步分解算法,同时采取一些技术处理,可确保模式稳定而有效地进行长时间数值积分、有效能量分析表明,与海面起伏相联系的有效表面位能是总有效能量的重要组成部分,它与动能同量级甚至更大些,因此在OGCM中须考虑有效表面位能与其它能量间的转换机制,这在物理上更加合理,而采用“刚盖”近似则是不合理的;同时,扣除海洋标准层结,缩小了位能与动能之间的量级差别,有利于长期数值积分的稳定进行.此外,这个模式能直接计算海面起伏和分析有效能量,使得表面流场和上层海洋热储率的模拟比已有文献的结果有明显的改进.除模式设计外,本部分只给出模式的基本性能,系统的模拟结果将在第Ⅱ,Ⅲ两部分发表.  相似文献   

Fresnel zone inversion for lateral heterogeneities in the earth   总被引：2，自引：0，他引：2  

Kiyoshi Yomogida 《Pure and Applied Geophysics》1992,138(3):391-406

We propose a different kind of seismic inversion from travel-time or waveform inversion for lateral heterogeneities in the earth: Fresnel zone inversion. Amplitude and phase delay of data in several frequency ranges are inverted for model space around ray paths with a width corresponding to the considered frequency so that primary effect of finiteness of wavelength be included. For vertically heterogeneous media, Fréchet derivatives for inversion are obtained very efficiently using the paraxial ray approximation, with nearly similar amounts of computation compared to travel-time inversion. As an example, Fréchet derivatives are computed for a teleseismic observation system for a three-dimensional structure in the lithosphere beneath an array of seismic stations. Even if the used frequency is around 2 Hz, the width of Fréchet derivatives cannot be neglected, particularly near the bottom of the lithosphere. Sensitivity of model parameters to observations is, moreover, different in our approach from conventional travel-time inversion: it is zero along ray paths but large slightly away from them. Some model calculations show that travel-time inversion, particularly with models divided into very fine meshes or blocks, might give misleading results. An example of inversion for a simple Camembert model, in the event that travel-time inversion gives no reliable results, shows how this technique works with much smaller data sets and computation than waveform inversions.  相似文献