Oracle 空间查询, 数据类型为 sdo_geometry

Posted 2020-10-24 疯子加天才

因网上搜索到的相关资料大部分都是关于sdo_geometry的介绍和以及通过sql语句添加要素，查询要素等等。没有找到存储过程相关的例子，所以只好自己动手啦。 

准备

环境：windowsxp系统，安装Arcgis 10.1，oracle11g32位，plsql。

数据：通过arcgis直连数据库，导入测试数据，导入时数据类型选择SDO_GEOMETRY。

测试数据包含点线面，点数据MAP_USER_POINT, 线数据MAP_USER_LINE,面数据MAP_USER_POLY。在arcmap中显示如下：

 

 

 

创建存储过程 添加几何要素

创建函数  创建 几何要素

create or replace FUNCTION FUN_CreateGeometry

( wkt IN VARCHAR2

  ,srid IN NUMBER

) RETURN sdo_geometry AS

BEGIN

  RETURN SDO_GEOMETRY(wkt,srid);

END FUN_CreateGeometry;

 

创建存储过程 向矢量表插入要素

create or replace

PROCEDURE AddNEWPOINT(

  OBJECTID IN NUMBER

  WKT IN VARCHAR2

  ,SRID IN NUMBER

) AS

 v_Return SDO_GEOMETRY;

BEGIN

  v_Return := FUN_CREATEGEOMETRY(

    WKT => WKT,

    SRID => SRID

  );

 

  INSERT INTO MAP_USER_8332(OBJECTID, SHAPE)

  VALUES

    (OBJECTID, v_Return);

END AddNEWPOINT;

 

调用存储过程创建几何，注意srid必须一致

 

Oracle SQL Developter中调用：

DECLARE

  OBJECTID NUMBER;

  WKT VARCHAR2(200);

  SRID NUMBER;

BEGIN

  OBJECTID := 100;

  WKT := \'POINT (118.731963311458 32.0661417793197) \';

  SRID := 4326;

 

  ADDNEWPOINT(OBJECTID, WKT, SRID);

END;

Sql plus中调用：

   EXECUTE ADDNEWPOINT(1,’POINT(118 32)’, 4326)

   Commit

 

空间查询

Sql语句给定要素500缓冲区之内的要素，可变部分为输入要素wkt（点，线，面），查询图层，缓存距离，示例如下：

 

创建包，定义一个cursor

CREATE OR REPLACE

PACKAGE PACKAGE_GEOMETRY AS

  TYPE GEOMETRY_CURSOR IS ref cursor;

  /* TODO enter package declarations (types, exceptions, methods etc) here */

END PACKAGE_GEOMETRY;

 

创建存储过程，返回查询结果cursor

create or replace

PROCEDURE SEARCHGEOMETRY(

 WKT VARCHAR2

,SRID NUMBER

,buflength VARCHAR2

,p_cursor in out PACKAGE_GEOMETRY.CURSOR_GEOMETRY) AS

 srarch_geom SDO_GEOMETRY;

 bufferCondition VARCHAR2(35);

BEGIN

    srarch_geom := FUN_CREATEGEOMETRY(

    WKT => WKT,

    SRID => SRID

  );

  bufferCondition := \'distance=\' || buflength;

  open p_cursor for SELECT c_c.objectid,c_c.NAME OBJNAME  FROM  MAP_USER_0000008332 c_c

    WHERE SDO_WITHIN_DISTANCE(c_c.SHAPE, srarch_geom, bufferCondition) = \'TRUE\';

   

END SEARCHGEOMETRY;

 

调用存储过程查询要素

Oracle SQL Developter或plSql中调用：

DECLARE

  BUFLENGTH VARCHAR2(25);

  WKT VARCHAR2(200);

  SRID NUMBER;

  RESULT_CUR PACKAGE_GEOMETRY.CURSOR_GEOMETRY;

  OBJECTID NUMBER;

  OBJNAME VARCHAR2(200);

BEGIN

 

  WKT := \'POINT (118.731963311458 32.0661417793197) \';

  SRID := 4326;

  BUFLENGTH := \'500\';

  SEARCHGEOMETRY(WKT, SRID,BUFLENGTH, RESULT_CUR);

   loop

        fetch RESULT_CUR

            into OBJECTID,OBJNAME;

        exit when RESULT_CUR%notfound;

        dbms_output.put_line(OBJECTID);

dbms_output.put_line(OBJNAME);

   end loop;

CLOSE RESULT_CUR;

END;

 

 

空间查询测试

查询点周围500米内的点数据

选取点数据内某一个点获取其wkt坐标：

SELECT CC.SHAPE.GET_WKT() FROM MAP_USER_POINT CC WHERE CC.NAME=\'湖心花园\'

POINT (118.748698730662 32.035395102502)   

修改SEARCHGEOMETRY内查询的表名为MAP_USER_POINT，

查询结果：

8

车站小区

9

文体西村

7

蓓蕾小区

11

迎宾村

10

茶亭村

12

湖心花园

19

园中园公寓

Arcmap查询结果：

 

Plsql查询结果与arcmap查询结果一致。

 

查询点周围500米内的线数据

选取点数据内某一个点获取其wkt坐标：

SELECT CC.SHAPE.GET_WKT() FROM MAP_USER_POINT CC WHERE CC.NAME=\'万科金色家园\'

POINT (118.759385095746 32.0420345318682)

  

修改查询参数wkt为’ POINT (118.759385095746 32.0420345318682)’, 修改存储过程SEARCHGEOMETRY内查询的表名为MAP_USER_LINE。

查询结果：

线段2被选中，arcmap中做该点的缓冲区查询，线段2被选中。修改查询缓冲区大小为700米，线段2、3被选中。

Plsql查询结果与arcmap查询结果一致。

 

查询点周围500米内的面数据

选取点数据内某一个点获取其wkt坐标：

SELECT CC.SHAPE.GET_WKT() FROM MAP_USER_POINT CC WHERE CC.NAME=\'名湖雅居\'

POINT (118.759724313726 32.0354580244683)

 

修改查询参数wkt为’ POINT (118.759724313726 32.0354580244683)’, 修改存储过程SEARCHGEOMETRY内查询的表名为MAP_USER_POLY。

查询结果：

2
face2
6
face6
7
face7
面2、6、7被选中，arcmap中做该点的500米缓冲区查询，面2、6、7被选中。

 

Plsql查询结果与arcmap查询结果一致。

 

查询线周围500米内的点数据

选取线数据内某一条线段获取其wkt坐标：

SELECT CC.SHAPE.GET_WKT() FROM MAP_USER_LINE CC WHERE CC.NAME=\'line3\'

查询所得wkt：

LINESTRING (118.763605759798 32.0391445982447, 118.763755695669 32.0386259340397, 118.763874074329 32.0381089704527, 118.763949962721 32.0377700573423, 118.764039433573 32.0374028695468, 118.764145034666 32.0369562671181, 118.764232941597 32.0365690568166, 118.76430597554 32.0362041074337, 118.764344298351 32.0359930356501, 118.764388237427 32.0357066978059, 118.764487386784 32.0350466349943, 118.764565439843 32.0345753407794, 118.76460084975 32.0343352406788, 118.764633191169 32.0341192388119, 118.764646064065 32.0339927329783, 118.764660872302 32.0337829885941, 118.764663602643 32.0336185295715, 118.764666619869 32.0334691063143, 118.764654318942 32.0330732652191, 118.764633596763 32.032851987929, 118.764589327635 32.0325688795502, 118.764551264729 32.0323638844857, 118.764504484694 32.0321339440263, 118.764435108294 32.0317689748584, 118.764351598148 32.0313420450993, 118.764308733761 32.0311321028643, 118.764242330519 32.0307380145478, 118.76422386474 32.0305728504567)

 

修改查询参数wkt为上一步查询到的wkt，修改存储过程SEARCHGEOMETRY内查询的表名为MAP_USER_POINT。

查询结果：

5

碧虹苑

6

草根居

22

君园

23

莫愁东寓小区

24

名湖雅居

25

万科金色家园

36

罗廊巷2号小区

37

环宇轩小区

38

中兴新村

39

止马营社区

44

陶李王巷新苑

45

汉中苑物管小区

46

汉中苑物管小区

47

止马村

48

通宇花园

49

韩家苑小区

50

韩家苑小区

         在arcmap中通过line3的500米缓冲区查询点要素，查询结果如下图：

 

Plsql查询结果与arcmap查询结果一致。

查询线周围500米内的线数据

Plsql查询结果与arcmap查询结果一致。

查询线周围500米内的面数据

PLSQL 查询结果：

1
face1
3
face3
6
face6
7
face7

arcmap查询结果如下图：

 

 

Plsql查询结果与arcmap查询结果一致。