Oracle spatial 空间数据SQL查询操作相关实例

Oracle spatial 空间数据SQL查询操作相关实例

--select dlbm,dlmc,trim(zldwdm) as zldwdm_1 from gzdt where nvl(zldwdm,'t')='t' or zldwdm='';
update gzdt set ZLDWDM='510113106' where  nvl(zldwdm,'t')='t' or zldwdm='';
commit;
delete from gzdt where ZLDWDM like '510113106%'
commit;
//获取空间面的面积
update gzdt set mj=mdsys.sdo_geom.sdo_area(geometry,0.0000000005);
//获取空间线的长度
update xzdw set cd=mdsys.sdo_geom.sdo_length(geometry,0.0000000005);

//删除空间数据 用SQL语句
//sql insert oracle spatial object  耕地
delete from spatial;
insert into spatial(dlbm,geometry)
select dlbm,geometry from v_dltb where dlbm in('011','012','013') ;
commit;

//插入空间数据 用SQL语句
insert into spatial(dlbm,geometry)
select dlbm,geometry from v_dltb where dlbm in('011','012','013') ;
commit;

//创建空间字段索引 oracle spatial table
//======================================
drop index index_spatial_v_gb_gdbhdk_h;
drop index index_spatial_v_jj_xzq_h;
drop index index_spatial_v_tdlygh_ytfq_xz_e;
drop index index_spatial_v_tdlyxz_dltb_h;
drop index index_spatial_v_tdly_nydfddj_k;


create index v_gb_gdbhdk_h_spatial_index      on v_gb_gdbhdk_h(geometry)        indextype mdsys.spatial_index;
create index v_jj_xzq_h_spatial_index         on v_jj_xzq_h(geometry)           indextype mdsys.spatial_index;
create index v_tdlygh_ytfq_xz_e_spatial_index on v_tdlygh_ytfq_xz_e(geometry)   indextype mdsys.spatial_index;
create index v_tdlyxz_dltb_h_spatial_index    on v_tdlyxz_dltb_h(geometry)      indextype mdsys.spatial_index;
create index v_tdly_nydfddj_k_spatial_index   on v_tdly_nydfddj_k(geometry)     indextype mdsys.spatial_index;
//======================================


//创建字段索引
//=======================================
drop   index index_fd_v_gb_gdbhdk_h_xzqdm;
drop   index index_fd_v_jj_xzq_h_xzqdm;
drop   index index_fd_v_tdlygh_ytfq_xz_e_xzqdm;
drop   index index_fd_v_tdlyxz_dltb_h_zldwdm;
drop   index index_fd_v_tdly_nydfddj_k_xzdm;


create index index_fd_v_gb_gdbhdk_h_xzqdm         on v_gb_gdbhdk_h(xzqdm);
create index index_fd_v_jj_xzq_h_xzqdm            on v_jj_xzq_h(xzqdm);
create index index_fd_v_tdlygh_ytfq_xz_e_xzqdm    on v_tdlygh_ytfq_xz_e(xzqdm);
create index index_fd_v_tdlyxz_dltb_h_zldwdm      on v_tdlyxz_dltb_h(zldwdm);
create index index_fd_v_tdly_nydfddj_k_xzdm       on v_tdly_nydfddj_k(xzdm);
//=======================================


//读取空间数据字段sql geometry 
select DLBM,dlmc,
       mdsys.sdo_geom.sdo_area(geometry,0.0000000005) as geo_mj,
       sdo_util.getnumelem(geometry) as num_elem,
       sdo_util.getVertices(geometry) as Vertices,
       sdo_util.GetNumRings(geometry) as Num_Rings,
       sdo_util.to_gmlgeometry(geometry) as gmlgeo,
       geometry
from v_dltb


//两空间图层相交运算
//任意相交运算mask=anyinteract
delete from gzdt;
insert into gzdt(dlbm,geometry)
select a.dlbm,
     SDO_GEOM.SDO_INTERSECTION(a.geometry, b.geometry, 0.0001) as geometry
from v_dltb as a
     v_ytfq as b
where sdo_relate(a.geometry,b.geometry,'mask=ANYINTERACT')='TRUE'


//在内部运算mask=inside
delete from gzdt;
insert into gzdt(dlbm,geometry)
select a.dlbm,
     SDO_GEOM.SDO_INTERSECTION(a.geometry, b.geometry, 0.0001) as geometry
from v_dltb as a
     v_ytfq as b
where sdo_relate(a.geometry,b.geometry,'mask=INSIDE')='TRUE'




//dltb_jbnt叠加分析
select * from v_dltb
where dlbm in('011','012','013') and dldwdm like '510112106%';
//
select d.dlbm,d.dlmc,
       d.tbmj,d.tbdlmj,d.xzdwmj,d.lxdwmj,d.tkmj,
       mdsys.sdo_geom.sdo_area(d.geometry,0.0000000005) as geo_mj,
       sdo_util.getnumelem(d.geometry) as num_elem,
       sdo_util.getVertices(d.geometry) as Vertices,
       sdo_util.GetNumRings(d.geometry) as Num_Rings,
       sdo_util.to_gmlgeometry(d.geometry) as gmlgeo,
       SDO_GEOM.SDO_INTERSECTION(d.geometry, y.geometry, 0.0001) as geometry 
from v_dltb d,
     v_ytfq y
where  d.dldwdm like '510112106%'  and  (d.dlbm in('021') or d.dlbz in('k','K')) and
       y.xzqdm like '510112%' and
       mdsys.sdo_geom.relate(d.geometry,'INSIDE',y.geometry,0.0001)='INSIDE';
//提取v_gbjj图层有效几何图形数据
select * from v_gbjj
where sdo_geom.validate_geometry_with_context(GEOMETRY,0.0001)='TRUE'
//提取v_gbjj图层无效几何图形数据
select * from v_gbjj
where sdo_geom.validate_geometry_with_context(GEOMETRY,0.0001)<>'TRUE'
//==the==end==