Oracle Spatial的简单应用
05年的时候就用过这个东西，不过当时完成了项目就了事，没把一些操作记录下来，这两天又用了一下，记录一下一个简单应用的操作步骤吧。
- --建立测试表
- create table test_geo (id number primary key, name varchar2(32), shape mdsys.sdo_geometry);
- --插入metadata
- insert into user_sdo_geom_metadata
- (table_name, COLUMN_NAME, DIMINFO, SRID)
- values
- ('test_geo',
- 'SHAPE',
- MDSYS.SDO_DIM_ARRAY(MDSYS.SDO_DIM_ELEMENT('X',
- -180.000000000,
- 180.000000000,
- 0.500000000),
- MDSYS.SDO_DIM_ELEMENT('Y',
- -90.000000000,
- 90.000000000,
- 0.500000000)),
- 8307);
- --建立索引
- CREATE INDEX test_geo_shape_idx ON test_geo(shape)
- INDEXTYPE IS MDSYS.SPATIAL_INDEX;
- --插入数据
- --这是一个矩形
- insert into test_geo
- (id, name, shape)
- values
- (2,
- 'test_area_2',
- mdsys.sdo_geometry(2003,
- 8307,
- null,
- mdsys.sdo_elem_info_array(1, 1003, 1),
- mdsys.SDO_ORDINATE_array(118, 32, 119, 32, 119,33, 118,33, 118,32)));
- --这是一个不规则多边形
- insert into test_geo
- (id, name, shape)
- values
- (3,
- 'test_area_3',
- mdsys.sdo_geometry(2003,
- 8307,
- null,
- mdsys.sdo_elem_info_array(1, 1003, 1),
- mdsys.SDO_ORDINATE_array(1,3, 2,2, 3,3, 4,2, 3,1, 2,1, 0,2, 1,3)));
- --执行查询
- --这个能查询到test_area_2
- select *
- from test_geo t
- where SDO_RELATE(t.shape,
- MDSYS.Sdo_Geometry(2001,
- null,
- MDSYS.Sdo_Point_Type(118.5, 32.5, null),
- null,
- null),
- 'mask=ANYINTERACT querytype=WINDOW') = 'TRUE';
- --这个能查询到test_area_3，说明落在地图边界上的点也是可以查询的
- --替换查询条件为2，2可查到，2，1.5可查到，2，2.5查不到
- select *
- from test_geo t
- where SDO_RELATE(t.shape,
- MDSYS.Sdo_Geometry(2001,
- 8307,
- MDSYS.Sdo_Point_Type(1.5, 2.5, null),
- null,
- null),
- 'mask=ANYINTERACT querytype=WINDOW') = 'TRUE';
--建立测试表
create table test_geo (id number primary key, name varchar2(32), shape mdsys.sdo_geometry);
--插入metadata
insert into user_sdo_geom_metadata
(table_name, COLUMN_NAME, DIMINFO, SRID)
values
('test_geo',
'SHAPE',
MDSYS.SDO_DIM_ARRAY(MDSYS.SDO_DIM_ELEMENT('X',
-180.000000000,
180.000000000,
0.500000000),
MDSYS.SDO_DIM_ELEMENT('Y',
-90.000000000,
90.000000000,
0.500000000)),
8307);
--建立索引
CREATE INDEX test_geo_shape_idx ON test_geo(shape)
INDEXTYPE IS MDSYS.SPATIAL_INDEX;
--插入数据
--这是一个矩形
insert into test_geo
(id, name, shape)
values
(2,
'test_area_2',
mdsys.sdo_geometry(2003,
8307,
null,
mdsys.sdo_elem_info_array(1, 1003, 1),
mdsys.SDO_ORDINATE_array(118, 32, 119, 32, 119,33, 118,33, 118,32)));
--这是一个不规则多边形
insert into test_geo
(id, name, shape)
values
(3,
'test_area_3',
mdsys.sdo_geometry(2003,
8307,
null,
mdsys.sdo_elem_info_array(1, 1003, 1),
mdsys.SDO_ORDINATE_array(1,3, 2,2, 3,3, 4,2, 3,1, 2,1, 0,2, 1,3)));
--执行查询
--这个能查询到test_area_2
select *
from test_geo t
where SDO_RELATE(t.shape,
MDSYS.Sdo_Geometry(2001,
null,
MDSYS.Sdo_Point_Type(118.5, 32.5, null),
null,
null),
'mask=ANYINTERACT querytype=WINDOW') = 'TRUE';
--这个能查询到test_area_3，说明落在地图边界上的点也是可以查询的
--替换查询条件为2，2可查到，2，1.5可查到，2，2.5查不到
select *
from test_geo t
where SDO_RELATE(t.shape,
MDSYS.Sdo_Geometry(2001,
8307,
MDSYS.Sdo_Point_Type(1.5, 2.5, null),
null,
null),
'mask=ANYINTERACT querytype=WINDOW') = 'TRUE';
还有一点疑惑写在下面
我最早接触到spatial的时候，地图数据是通过oracle的shp2sdo这个工具，从shp格式的地图里导入的，那时候做查询的时候，查询条件可以用MDSYS.Sdo_Geometry(2003, 8307,
MDSYS.Sdo_Point_Type(1.5, 2.5, null), null, null)， 但是现在自己插入的数据，如果第一个参数还用2003就会报一个ora-13031异常，必须改成2001才能正常查询，不知道是为什么。