博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
Oracle 空间查询, 数据类型为 sdo_geometry
阅读量:6991 次
发布时间:2019-06-27

本文共 4796 字,大约阅读时间需要 15 分钟。

因网上搜索到的相关资料大部分都是关于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查询结果一致。

转载地址:http://uibvl.baihongyu.com/

你可能感兴趣的文章
分享非常漂亮的WPF界面框架源码及其实现原理
查看>>
如何获取ResultSet的行数和列数(转)
查看>>
绑定列ORA-24816: 在实际的 LONG 或 LOB 列之后提供了扩展的非 LONG 绑定数据
查看>>
Mobile Web调试工具Weinre
查看>>
Android巴士转发
查看>>
未能进入中断模式,原因如下:源文件“XXXXXX”不属于正在调试的项目。
查看>>
linux编程掌握常用命令
查看>>
asp.net获取客户端IP,Mac
查看>>
Linux内核-内核线程
查看>>
zoj 1152 A Mathematical Curiosity
查看>>
UML时序图总结
查看>>
【2013Esri全球用户大会精彩案例】Horry Count GIS--南卡罗来那州霍里县企业级应用...
查看>>
c++虚函数表 Brew VTBL
查看>>
SQL Server 2008开启sa账户以及如何用JDBC进行连接
查看>>
读取同一文件夹下多个txt文件中的特定内容并做统计
查看>>
为sourceinsight添加makefile、kconfig、*.S文件支持
查看>>
sharepoint2010问卷调查(1)-实现问卷的图片调查(采用JS实现)
查看>>
linux下如何挂接(mount)光盘镜像文件、移动硬盘、U盘、Windows网络共享和NFS网络共享...
查看>>
python开发_函数的参数传递
查看>>
利用mysqldump 实现每天备份方案
查看>>