栅格数据地图投影动态分块转换研究.docVIP

栅格数据地图投影动态分块转换研究 　　[摘 要]栅格数据是地理信息系统数据源的重要形式,而地图投影转换作为地理数据处理的重要环节,在数据入库、数据发布等工作中不可或缺。本文以地图投影转换原理为依据,设计了一种地图投影动态分块转换系统,并辅以实验说明。最后提出了借助分布式思想进一步提高转换速度的思路。 　　[关键词]栅格数据投影转换动态分块分布式 　　[中图分类号]TP17[文献标识码]A[文章编号]1007-9416(2010)03-0097-02 　　 　　引言 　　地球体的自然表面凹凸不平,形态复杂,不能直接将这样一个球体表面的数据映射到二维平面上,必须借助参考椭球与符合局部地理范围的基准面并按照一定的数学法则来完成。由于投影的目的与所要应用的地理范围不同,所采用的对地球体进行模拟表达的参考椭球与基准面也不尽相同,从而形成不同的空间参考系统。我们知道,完整的地理信息不仅应该包括数据体,而且应该包含数据所对应的空间参考信息。实际中,不同的地图比例尺往往采用不同的投影方式,如我国基本比例尺地形图中,大于等于1:50万的均采用高斯―克吕格投影(Gauss-Kruger),小于1:50万的地形图采用正轴等角割圆锥投影。 　　由于应用需求的经常变化,时常需要将同一种数据在不同的空间参考系统下进行表达以进行数据存储或显示,从而对地理信息数据投影转换提出了需求,而且动态投影显示对投影切换的速度要求更高。对于栅格数据而言,像素点的数据量极大,如果简单的进行一一坐标对应转换,无论是从转换的效率还是转换结果的连续性上分析,都是很难达到要求的,而简单的固定大小分块转换,则会造成数据相邻数据衔接的不连续性,不能较好的满足实际需求。所以对于栅格数据的投影转换,需要考虑更多的存储与显示相关的问题。 　　地图投影转换的本质是将A坐标参考系统下点坐标位置信息P(x,y)转换为B坐标参考系统下的P’(x,y)来表示。现在流行的商用软件如:arcgis,mapgis等都支持栅格数据的投影转换,但是存在一些列问题:灵活性差,无法根据机器物理配置动态调整转换速度;扩展性差,无法根据自定义的数据需求来进行新的转换,如google merc投影与经纬度表示的数据之间进行转换。 　　1 地图投影转换模型与数据结构分析 　　完成单幅地图投影转换分为以下步骤进行 　　1.1 根据空间参考生成转换器 　　转换器是能够独立完成一组坐标变换的算子,它是完成整个数据转换的可重用的模板。转换器由转换函数与参考系统等信息组成。以下是转换器的结构定义。其中poForwardTransform与poReverseTransform为封装了原始数据与结果数据空间参考坐标系统的正向转换与逆向转换,transinfoFunc是坐标转换函数,内含正向转换与逆向转换。 　　struct ProjectionTransformer 　　{ProjectionTransformerFunc transinfoFunc; //转换函数 　　Transformation *poForwardTransform; //正向转换 　　Transformation *poReverseTransform; 　　}//逆向转换 　　1.2 结果数据范围与仿射信息计算 　　创建空的结果数据集,首先需要获取结果数据集的范围与仿射系数信息。通过对原始数据边界进行高密度插值,可以计算出到结果数据的较为精确的范围,从而创建有效的结果数据集。 　　1.3 地图投影转换对象WarpOption的构建 　　typedef struct{double dfWarpMemoryLimit; 　　//内存限制 　　WarpResampleAlg eResampleAlg; 　　//抽样算法选择 　　WarpDataTypeeWorkingDataType; 　　//数据类型 　　WarpDatasetH 　　hSrcDS; 　　//原始数据集 　　WarpDatasetH 　　hDstDS; 　　//结果数据集 　　int 　　nBandCount; 　　//数据波段数 　　int 　　*panSrcBands; 　　//原始数据处理波段信息 　　int 　　*panDstBands; 　　//结果数据处理波段信息 　　WarpTransformerFunc pfnTransformer; //转换函数 　　void 　　*pTransformerArg; 　　//转换模板 　　} WarpOption; 　　1.4 动态分块 　　根据可用内存大小限制分块大小参数,设定分块Block大小。Block的大小决定了在内存中单次所能完成的数据读