基于postgis的空间数据计算.docxVIP

利用 PostGis 的空间数据库自带的分析计算功能，采用 WebService 技术 在 GIS 服务器端发布 GIS 计算服务。GIS 计算服务提供多种参考系和多种精度 的几何计算和地理计算功能， 客服端可采用任意开辟语言连接 GIS 计算服务， 调用对应的计算接口来完成空间计算功能。 一、服务端设计 1 、 计算点点之间的距离 1.1 计算点点距离的方法 1.1.1 在点与点之间计算时， 采用 ST_Distance(geometry, geometry)，计算 两个点在特定参考系下的距离， sql 示例语句如下。 1.1.2 也可以将两个点传入， 组成在特定参考系下的线， 然后调用 postgis 中 ST_Length(geometry)函数计算这条线的长度。 sql 示例语句如下 1.2 需要向外开放的参数有。 1) 计算距离的两个点坐标 A (x1,y1)，B (x2,y2) ; 2) 指定的参考系代号，如 4326 ，900913 等; 2 、 计算多点线的长度 2.1 在计算多点组成的线时， 需要调用 postgis 中的 ST_Length(geometry) 函数， 需要注意的是， 在组成这条线的点要故意义， 比如 LINESTRING (0 0,0 0)这个线就是无意义的。具体的 sql 示例如下： 2.2 需要向外开放的参数有： 1) 组成一条线的若干个点 A (x1,y1)，B (x2,y2)?? 2) 指定的参考系代号，如 4326 ，900913 等; 3 、 计算多面要素的面积 3.1 在计算多面要素的面积时，需要做个限定，即，在组成多面的几个 面的关系属于包含于第一面内，否则无任何几何意义。示意图如下： 从上图可知，形如 multipolygon (((面 1)，(面 2)，(面 3)，(面 4))) 的一个多面图形中， “面 1”代表一个整体的多边行范围， 而其余几个面为其 内部的一些多边形，其结果面积为上图阴影部份所示。因此计算多面要素面 积的思路可以理解为，首先计算出“面 1”的面积 S1，然后在根据其内部几 个多边形的几何关系求出内部多边形的面积之和 S2 (如果两面相交， 则求出 其相交后的面积(即面相并面积) ；如果两面相离，则求出其各自的面积) 。 则图中阴影部份的面积为 S1-S2 。S2 的面积可以参考小一节的多个面相并后 的面积计算方法。 4 、 计算多个面相合并后的面积 在给定多个面后，求其相并后的面积，这就需要先将这些面进行合并， 形成一个整体的图形，这个图形是去掉重叠部份的图形，示意图如下所示： 图中面 1 的区域与面 2 的区域相交形成面 4，并且两个都与面 3 区域相 离，则这三个面所组成的区域的面积为 S1+S2-S4+S3。在 postgis 中，函数 ST_Union(geometry,geomtry)可以获得两个几何图形的合并结果。而 ST_Union(geometry set) 是获取一系列几何图形的合并结果，但是经过测验， 后者并没有实现预期的合并结果。因此，需要在程序端实现时，先使用 ST_Union(geometry,geomtry)将需要合并的多个面状几何图形进行合并， 然后 在对合并后的几何图形进行面积的计算。具体的 sql 如下 首先初始化第一个面： 然后循环合并面： 最后计算合并好的几何面的面积： 5 、 计算多个面相交部份面积 求多个面相交部份的面积的方式与求多个面相并后的面积方式类似， 首 先需要使用 postgis 中的 ST_Intersection (geometry ，geometry)函数，求出 两两的相交部份，然后在利用 ST_AREA 求出这个相交部份的面积。 6 、 计算简单多边行面积 关于计算简单多边形的面积时，需要注意 postgis 中关于面的定义形式， 如 polygon((0 0,5 0,5 5,0 5,0 0))，其中的坐标点要首尾相接。 示例 sql 如下图 所示。 需要向外提供的参数有 1) 围成一个多边形的的若干个点 A (x1,y1)，B (x2,y2)?? 2) 指定的参考系代号，如 4326 ，900913 等; 二、客户端调用服务 在客户端调用时，目前分为两种方式。 1、在 java 中调用 GetWSByAxis2 类中的响应方法 (根据返回值调用相 关函数)，然后将需要调用的 web service 服务中的函数，以及函数需要 的参数传入。此处需要注意的是，对于 web service 中的函数的参数要给 予正确的格式。如多个面求交求并，则参数的格式为： 在代表面的字符串中间，需要插入