球面退化四叉树格网的剖分及变形分析

为有效解决经纬度格网与四元三角网(QTM)在全球空间数据管理与操作中存在的不足,提出一种新的全球离散格网剖分方法——球面退化四叉树格网(DQG),给出该格网的剖分原理与编码规则,并对其进行几何变形计算与分析。结果表明:DQG既具有经纬度格网结构简单的特征,又具有QTM几何变形稳定等优点。     

基于球面DQG的数字高程建模算法与可视化表达

为有效解决经纬度格网与四元三角网(Quaternary Triangular Mesh,QTM)在全球地形建模方面存在的不足,根据球面DQG(Degenerate Quadtree Grid)的几何结构特点,选择双线性多项式内插方法进行格网点高程内插,给出基于球面DQG的DEM建模算法与效率分析,并应用美国地质调查局提供的GTOPO30全球地形数据进行相关实验。结果表明:全球DEM建模时,DQG所需的格网数是经纬度格网或QTM的2/3左右,且可视化操作(图形放大、缩小及漫游等)时画面平滑、流畅,没有抖动。     

一种基于球面QTM格网的面要素边界跟踪填充算法

为实现四元三角网(Quaternary Triangular Mesh,QTM)对球面面状要素矢量数据的离散化及可视化,该文将基于栅格单元的"边界代数法"扩展到球面三角格网单元,提出了一种基于QTM格网的边界跟踪填充算法,主要内容包括:通过ETP投影实现球面QTM格网与平面三角格网的相互转换,建立矢量线角度与三角格元邻近搜索的对应关系;"边界跟踪"并记录边界格元进行分类处理;最终给出"边界跟踪"下两种不同的填充模式。实验表明,该算法实现了球面矢量多边形在QTM格网中的填充。     

一种改进的DQG格网编码与经纬度高效转换算法

格网编码与经纬度间的转换效率是影响格网应用的关键因素。现有转换算法难以满足海量空间大数据的实时计算与分析需求,效率有待提升。为此,该文提出一种改进的DQG格网编码与经纬度高效转换算法:首先,引入二进制DQG格网编码代替四进制编码;然后,根据DQG格元的分布特征推导出格元行号和经度差的关系,以此改进格元列号的计算过程;最后,使用查找表进行Morton码的编码和解码,以提高转换效率。实验结果表明,该文提出的DQG格网编码与经纬度转换算法的平均效率分别是四进制DQG、单层二维SDZ算法的20.15倍和4.58倍,基本满足海量格网数据的实时转换与计算需求,为海量空间大数据的高效计算与分析奠定了基础。     

基于Laplace方程的球面正交曲线格网生成方法

为克服基于传统经纬度格网的全球海洋及大气数值模式产生的极点问题,同时顾及全球格网编码便捷及信息查询,该文将立方体球心投影的格网经球极投影转换到二维平面区域,再应用Laplace方程数值求解生成正交曲线格网;然后采用球极逆投影,得到球面上一种新的正交曲线格网。通过与等角球心投影、保角变换两种球面格网生成方法的结果相比较,从格网两边的长度、单元的面积、正交性等指标进行统计分析,对所生成格网的质量进行评价。结果表明,等角球心投影法生成的格网单元大小较均匀,但正交性不理想;Laplace方程法生成的正交曲线格网稍优于保角变换法,可作为应用有限差分方法和有限体积方法的全球海洋及大气大尺度数值模拟的基础格网。     

基于QTM的海平面上升分析与模拟

针对海平面上升影响范围分析与模拟这一国际前沿问题,通过研究基于球面四元三角网(QTM)的关键技术问题,包括复杂拓扑关系计算、LOD剖分、球面水淹分析、基于QTM的多分辨率的DEM数据组织方法和分析精度的相关性评定等,以.Net和Direct3D为开发工具,设计开发了基于QTM的海平面上升影响范围评估模型。该研究结果可为全球海平面上升影响的防灾减灾决策提供有效支撑,并推动了球面数据模型和球面格网拓扑分析的理论成果在全球变化预测相关领域的应用进展。     

平面影像到QTM像元的转换算法及精度分析

平面影像数据到四元三角网(Quaternary Triangular Mesh,QTM)像元的转换是实现全球QTM格网影像无缝建模的首要任务。提出基于面积权的数据转换方法,给出不同类型的权重及格网值的计算,并利用ArcGIS中256级灰度的Wsiearth.tif图像数据分析了转换精度。结果表明:将4 km分辨率的影像数据转换到12层的QTM像元时,转换误差在0~2灰度内占96.27%,对一般应用的影响可以忽略不计。     

全球离散格网的建模方法及基本问题

全球离散格网(Discrete Global Grid,DGG)模型是数字地球及空间信息格网的基础,不同的建模方法不但影响空间数据的存储和管理效率,而且影响全球GIS的操作功能。该文介绍了DGG的评价标准,将DGG的建模方法归纳为3种类型:经纬度格网模型、自适应格网模型和正多面体格网模型,重点分析了不同类型球面离散格网模型的几何结构、单元特征和应用模式。最后,提出了DGG在Global GIS中亟待解决的基本问题,包括编码、精度、应用、误差、整合和定位问题。     

球面似均匀格网的剖分方法及特征分析

给出球面"似均匀格网"概念的内涵,并介绍3种常用的球面似均匀格网(三角形格网、菱形格网和六边形格网)的剖分方法.从几何上分析3种格网的单元特征及不同格网在不同领域的应用优势.通过对比得出3种格网的特征对照表及相互关系,即:三角形格网是最基本的球面剖分格网,球面六边形格网和菱形格网可以通过球面三角格网来生成,且这3种剖分格网各有优缺点并在一定条件下可以相互转换.     

菱形三十面体六边形全球离散格网系统构建算法

全球离散格网系统是递归剖分整个地球空间形成的多分辨率离散地球参考模型,有助于建立以空间位置为主键的多源数据关联融合新模式。正二十面体的三角面不利于建立格网单元行列索引,且生成的球面格网单元存在较大形状差异,对此,该文采用面数更多、结构更优的菱形三十面体构建六边形全球离散格网系统。首先,在单个菱形面上建立离散斜坐标系描述多种孔径六边形剖分产生的多分辨率格网;其次,根据菱形三十面体的几何属性确定跨面单元归属,据此拼接各面上的多分辨率单元,构成封闭三十面体上的格网系统;最后,通过Snyder等积多面体投影将封闭菱形三十面体上的格网系统映射至球面。实验结果表明,相比二十面体的同类成果,菱形三十面体六边形全球离散格网系统的几何属性和应用潜力更好,有助于建立地球空间信息数据的高精度承载模型。     

SRG剖分编码与地理坐标的转换算法研究

提出一种基于SRG剖分方法的剖分编码与经纬度坐标之间的转换算法,不涉及任何投影变换,计算过程只需加、减、乘、除简单算术运算,计算速度快,且SRG剖分编码本身就具有固定的方向性,利于邻近搜索。算法中采用坐标系辅助区分一些难以区分的菱形块,降低了转化误码率,进而提高了转换的精度。     

World population in a grid of spherical quadrilaterals

"We report on a project that converted subnational population data to a raster of cells on the earth. We note that studies using satellites as collection devices yield results indexed by latitude and longitude. Thus it makes sense to assemble the terrestrial arrangement of people in a compatible manner. This alternative is explored here, using latitude/longitude quadrilaterals as bins for population information.... The results to date of putting world boundary coordinates together with estimates of the number of people are described. The estimated 1994 population of 219 countries, subdivided into 19,032 polygons, has been assigned to over six million five minute by five minute quadrilaterals covering the world."     

Computation of complete synthetic seismograms for laterally heterogeneous models using the Direct Solution Method

We use the Direct Solution Method (DSM) together with the modified operators derived by Geller & Takeuchi (1995) and Takeuchi, Geller & Cummins (1996) to compute complete synthetic seismograms and their partial derivatives for laterally heterogeneous models in spherical coordinates. The methods presented in this paper are well suited to conducting waveform inversion for 3-D Earth structure. No assumptions of weak perturbation are necessary, although such approximations greatly improve computational efficiency when their use is appropriate.
An example calculation is presented in which the toroidal wavefield is calculated for an axisymmetric model for which velocity is dependent on depth and latitude but not longitude. The wavefield calculated using the DSM agrees well with wavefronts calculated by tracing rays. To demonstrate that our algorithm is not limited to weak, aspherical perturbations to a spherically symmetric structure, we consider a model for which the latitude-dependent part of the velocity structure is very strong.     

川渝地区气温随地形、经度和纬度的变化(英文)

Using the daily temperature data of 95 meteorological stations from Si-chuan-Chongqing Region and its surrounding areas, this paper adopted these methods (e.g., linear regression, trend coefficient, geographical statistics, gray relational analysis and spatial analysis functions of GIS) to analyze the relations of temperature variability with topography, latitude and longitude. Moreover, the rank of gray correlation between temperature variability and elevation, longitude, latitude, topographic position and surface roughness also was measured. These results indicated: (1) The elevation affected temperature variability most obviously, followed by latitude, and longitude. The slope of the linear regression between temperature change rate and elevation, latitude and longitude was 0.4142, 0.0293 and -0.3270, respectively. (2) The rank of gray correlation between temperature change rate and geographic factors was elevation > latitude > surface roughness > topographic position > longitude. The gray correla-tion degree between temperature change rate and elevation was 0.865, followed by latitude with 0.796, and longitude with 0.671. (3) The rate of temperature change enhanced with the increase of elevation. Especially, the warming trend was significant in the plateau and mountain areas of western Sichuan, and mountain and valley areas of southwestern Sichuan (with the warming rate of 0.74℃/10a during the 1990s). However, there was a weak warming trend in Sichuan Basin and its surrounding low mountain and hilly areas. (4) The effects of latitude on temperature change rate presented the specific regulation, which the warming rate of low-latitude areas was more significant than that of high-latitude areas. However, they were consistent with the regulation that the increasing of low temperature controlled most of the warming trend, due to the effects of terrain and elevation on annual mean temperature. (5) Ba-sically, temperature variability along longitude direction resulted from the regular change of elevation along longitude. It was suggested that, in Sichuan-Chongqing Region, special features of temperature variability largely depended on the terrain complexity (e.g., undulations, mutations and roughness). The elevation level controlled only high or low annual mean temperature and the range of temperature change rate in the macro sense.     

The Sinusoidal Projection: A New Importance in Relation to Global Image Data

The reprojection of image data causes the loss or duplication of original pixel values. This research investigated the feasibility of using the sinusoidal projection for global image database construction. Specifically, reprojection accuracies were tested with geographic latitude and longitude coordinates, and the Hammer‐Aitoff, Eckert IV, Mollweide, and sinusoidal projections. Reprojections between these five global projections and the Universal Transverse Mercator (UTM) projection and referencing system were performed using fifty‐four sample datasets. A statistical analysis of categorical accuracy, a measure describing the omission of pixel values during reprojection, was conducted. Geographic coordinates and the sinusoidal projection both showed very high accuracy rates (100.0 percent and 99.5 percent respectively) when sample datasets were reprojected from UTM. The geographic coordinates, however, showed very low accuracy (65.3 percent) when sample datasets were reprojected to the UTM projection, while the sinusoidal projection showed the highest accuracy (98.4 percent). The results strongly suggest that the sinusoidal projection is very accurate and efficient for building global image databases.     

基于GIS的定位观测数据空间化

属性数据空间化是当前GIS领域的前沿问题之一。在对中国624个气象站多年平均气温数据空间化过程中,通过使用30秒分辨率的数字高程数据,把气温分解为受经纬度、海拔高度影响的规律性成分和受其它因素影响的非规律性成分两部分,并分别用多元回归和反距离权重内插方法对二者实施空间化,最后将空间化结果进行合成得到基于栅格的中国多年平均气温分布数据。该数据既能反映气温在空间上的宏观变化,又能反映气温在局部地区的微观变化。该方法可供其它类型观测数据空间化、特别是在观测站点稀疏的情况下参考和借鉴。     

1947年中国南海断续线精准划定的地形依据

基于1947年中国政府内政部印制的1:400万《南海诸岛位置图》,采用地理信息系统及数字化仿射变换精准测量出11段断续线的地理坐标,同时利用4.4万多个实测水深点数据建立南海海底三维地形模型,通过叠加分析获得了1947年南海断续线所在位置的海底地形特征。结果显示：断续线东长西短,间距东西大、南北小,其走向、形态与下伏地形特征相符合,主要分布在大陆坡和大陆架上;并基于断续线下伏海底地形分析,推断出当时断续线划定时可能遵循如下依据：① 在浅海盆地、海峡处采用“等距离中间线”原则进行划分;② 在具有明显地形变化处沿海槽坡麓或海槽槽沟中线进行划分;③ 在西南部陆坡陆架区,结合岸线走向及下伏地形特征进行划分。