使用SQL直接操纵FeatureClass(oracle),

zldessktin

摘要：使用SQL直接操纵FeatureClass(oracle),这里主要讲存储类型为ST_Geometry的要素类。对FeatureClass的操作主要包括下面几点：1数据的插入，删除，更新；数据的插入直接使用insert语句来进行，构造ST_Geometry的时候可以通过两种方法来完成:1）使用WKT编码 2）使用WKB编码 使用SQL直接操纵FeatureClass(oracle),这里主要讲存储类型为ST_Geometry的要素类。
对FeatureClass的操作主要包括下面几点：
1数据的插入，删除，更新
数据的插入直接使用insert语句来进行，构造ST_Geometry的时候可以通过两种方法来完成:
1）使用WKT编码
2）使用WKB编码
上面两种编码都是OGC规范的编码方式，分别通过ST_PolyFromText()和ST_PointFromWKB() 以及一系列类似函数还完成从WKT或WKB到ST_Geometry转换。
下面为所有相关ST函数：
ST_GeomFromText—Creates an ST_Geometry from a text representation of any geometry type
ST_PointFromText—Creates an ST_Point from a point text representation
ST_LineFromText—Creates an ST_LineString from a linestring text representation
ST_PolyFromText—Creates an ST_Polygon from a polygon text representation
ST_MPointFromText—Creates an ST_MultiPoint from a multipoint representation
ST_MLineFromText—Creates an ST_MultiLineString from a multilinestring representation
ST_MPolyFromText—Creates an ST_MultiPolygon from a multipolygon representation
ST_AsText—converts an existing geometry into a text representation
ST_GeomFromWKB—Creates an ST_Geometry from a WKB representation of any geometry type
ST_PointFromWKB—Creates an ST_Point from a point WKB representation
ST_LineFromWKB—Creates an ST_LineString from a linestring WKB representation
ST_PolyFromWKB—Creates an ST_Polygon from a polygon WKB representation
ST_MPointFromWKB—Creates an ST_MultiPoint from a multipoint WKB representation
ST_MLineFromWKB—Creates an ST_MultiLineString from a multilinestring WKB representation
ST_MPolyFromWKB—Creates an ST_MultiPolygon from a multipolygon WKB representation
ST_AsBinary—converts an existing geometry value into well-known binary representation.
关于wkt 和 wkb的编码方式会在下一节讲。
下面就使用sql语句来进行数据的插入，执行下面的sql语句：
INSERT INTO sensitive_areas  (area_id , name, area_size, type , zone,fid)
VALUES (1, 'Summerhill Elementary School', 67920.64, 'school', sde.ST_PolyFromText('polygon
((52 28,58 28,58 23,52 23,52 28))', 0));
在sensitive_areas表中插入一条记录。
   
打开ArcMap，加载这个featureClass，可以看到新增加的这条记录。
   
更新记录
UPDATEsensitive_areas
SET zone= sde.st_pointfromtext('polygon ((52 30,58 30,58 23,50 23,50 28))', 0))
WHERE area_id = 1;
删除记录
DELETE FROM sensitive_areas WHERE names
(SELECT sa.names
FROM sensitive_areas sa, hazardous_sites hs
WHERE sde.st_overlaps (sa.zone, sde.st_buffer (hs.location,.01)) = 1);
2查询
在进行空间查询的时候会用到索引，其中infomix使用R树索引，PostgreSQ使用Generalized Search Tree (GiST) R-tree索引， oracle和 DB2 使用网格索引。
在oracle中进行查询会执行下面的步骤：
1）首先比较grid和查询范围，找出在查询范围内的所有grid。
2）找出在这些grid内的所有要素。
3）将这些要素的外包络矩形和查询范围比较，找出所有在查询范围内以及和查询范围相交的要素。
4）使用ST函数进行最终过滤（一般使用ST_Intersects或ST_Contains），找到完全符合条件的要素。
下面是一个进行空间查询的例子：
SELECT sa.name "Sensitive Areas", hs.name "Hazardous Sites"
FROM sensitive_areas sa, hazardous_sites hs
WHERE sde.st_overlaps (sa.zone, sde.st_buffer(hs.location,.01)) = 1;
3创建索引
CREATE INDEX sa_idx ON sensitive_areas(zone)
INDEXTYPE IS sde.st_spatial_index
PARAMETERS('st_grids=1,3,0 st_srid=0');

4创建空间参考系
ESRI建议最好使用ArcGIS Desktop来进行空间参考系的创建，但也可以使用SQL语句进行空间参考系的创建。执行下面的sql语句：
INSERT INTO SDE.ST_SPATIAL_REFERENCES VALUES (
   GCS_North_American_1983,
   1,
   -400,
   -400,
   1000000000,
   -100000,
   10000,
   -100000,
   10000,
   9.999E35,
   -9.999E35,
   9.999E35,
   -9.999E35,
   9.999E35,
   -9.999E35,
   9.999E35,
   -9.999E35,
   4269,
   'GCS_North_American_1983',
   'PROJECTED',
   NULL,
'GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID
["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT
["Degree",0.0174532925199433]]',
   'ArcSDE SpRef'
   );