转自:http://blog.csdn.net/wanglipo/article/details/5894249
首先介绍一下shapefile:
Shapefile文件是美国环境系统研究所(ESRI)所研制的GIS文件系统格式文件,是工业标准的矢量数据文件。 Shapefile将空间特征表中的非拓扑几何对象和属性信息存储在数据集中,特征表中的几何对象存为以坐标点集表示的图形文件—SHP文 件,Shapefile文件并不含拓扑(Topological)数据结构。一个Shape文件包括三个文件:一个主文件(*.shp),一个索引文件 (*.shx),和一个dBASE(*.dbf)表。
主文件是一个直接存取,变长度记录的文件,其中每个记录描述构成一个地理特征(Feature)的所有vertices坐标值。在索引文件中,每条记录包 含对应主文件记录距离主文件头开始的偏移量,dBASE表包含SHP文件中每一个Feature的特征属性,表中几何记录和属性数据之间的一一对应关系是 基于记录数目的ID。在dBASE文件中的属性记录必须和主文件中的记录顺序是相同的。图形数据和属性数据通过索引号建立一一对应的关系。
导入shapefile文件是通过ORACLE提供的一个工具shp2sdo.exe,网上导出可以找到。具体操作步骤如下:
1 下载shp2sdo.exe,把此文件复制到PATH变量包含的目录下,
如我的oracle客户端安装后自动注册的环境变量是 path :D:/oracle/product/10.2.0/client_1/BIN;
2 在dos(WINDOWS运行CMD)下定位到shp文件的目录,使用命令:cd /d 文件目录 进入需要操作的文件目录。(注意:文件目录不能含有汉字);
3 在文件目录下操作。输入命令为:shp2sdo.exe BOUNT_poly t_areainfo -g geom -d -x (-180,180) -y (-90,90)
-s 8307 -t 0.5 -v。 命令结束后会在当前目录生成t_areainfo.ctl 和t_areainfo.sql两个文件
4 连接数据库C:/Documents and Settings/Administrator>sqlplus user/password@orcl,登陆ORACLE成功后,执行脚本SQL> @t_areainfo.sql; 然后退出 SQL> exit;
5 导入控制文件t_areainfo.ctl : C:/Documents and Settings/Administrator>sqlldr user/password@orcl 回车,然根据提示control = 输入t_areainfo.ctl 回车,导入结束。
最后就是验证成果的时候,首先要创建MDSYS.SPATIAL_INDEX序列 语句如下:
CREATE INDEX areainfo_idx ON t_areainfo(geom) INDEXTYPE IS MDSYS.SPATIAL_INDEX;
接着进行查询,语句如下:
select * from t_areainfo t
where SDO_RELATE( t.geom,
MDSYS.Sdo_Geometry(2001, 8307, MDSYS.Sdo_Point_Type(113.119526, 34.84593, null), null, null ),
'mask=ANYINTERACT querytype=WINDOW' ) = 'TRUE';
结果如下,大功告成。
附注:各参数的含义
shp2sdo [-o] <shapefile> <tablename> -g <geometry column>
-i <id column> -n <start_id> -p -d
-x (xmin,xmax) -y (ymin,ymax) -s <srid>
or
shp2sdo -r <shapefile> <outlayer> -c <ordcount> -n <start_gid> -a -d
-x (xmin,xmax) -y (ymin,ymax)
shapefile - name of input shape file
(Do not include suffix .shp .dbf or .shx)
tablename - spatial table name
if not specified: same as input file name
Generic options:
-o - Convert to object/relational format (default)
-r - Convert to the relational format
-d - store data in the control file
if not specified: keep data in separate files
-x - bounds for the X dimension
-y - bounds for the Y dimension
-v - verbose output
-h or -? - print this message
Options valid for the object model only:
-g geometry column - Name of the column used for the SDO_GEOMETRY object
if not specified: GEOM
-i id_column - Name of the column used for numbering the geometries
if not specified, no key column will be generated
if specified without name, use ID
-n start_id - Start number for IDs
if not specified, start at 1
-p - Store points in the SDO_ORDINATES array
if not specified, store in SDO_POINT
-s - Load SRID field in geometry and metadata
if not specified, SRID field is NULL
-t - Load tolerance fields (x and y) in metadata
if not specified, tolerance fields are 0.00000005
-8 - Write control file in 8i format
if not specified, file written in 9i format
-f - Write geometry data with 10 digits of precision
if not specified, 6 digits of precision is used
Options valid for the relational model only:
-c ordcount - Number of ordinates in _SDOGOEM table
if not specified: 16 ordinates
-n start_gid - Start number for GIDs
if not specified, start at 1
-a - attributes go in _SDOGEOM table
if not specified, attributes are in separate table
本文转自http://www.blogjava.net/reeve/archive/2007/08/01/133759.html
最 近因为工作的关系,开始接触GIS,采用开源的Geoserver、Geotools、Mapbuilder来开发,数据库使用oracle spataral,刚开始是一直在看别人已经做好的一个半成品系统,自己模仿着来改,GIS采用的数据格式比较特殊,不能使用一般的导表方式来导,而且地 图文件的格式也比较特殊,一般是shp、或者是tab,这里主要是用的shp文件,因此需要将shp文件进行处理,转换成oracle数据格式,其实也没 什么难的,只是如果不熟悉的话,可能也需要浪费一点时间来看一下。
首先我们要把shp文件转换分割成oracle可以导入的数据库文件, 这个转换可以使用oracle提供的一个现成的工具shp2sdo(很小),可以到官方网站去找(不过可能不是很好找,需要的可以留下邮箱),下载后把此 文件复制到PATH变量包含的目录下,如我的oracle客户端安装后自动注册的环境变量是 path C:\Oracle\product\10.1.0\Client_1\bin;,我们可以把此文件拷贝到该目录下,然后在dos下运行该工具,定位到我 们的shp文件的位置,例如我们的shp文件名称是state.shp(shapefile包括至少三个文件state.shp state.dbf state.idx)在D:\data\目录下,我们在dos命令窗口下就应该这样操作
shp2sdo state statearea -i gid -s 8307 -g -d
其 中state即为本地的shp文件名(不加shp后缀),statearea为生成的文件名(同时也是将来导进数据库的表名),-d代表含义是将分解后的 ctl文件(控制文件)和data文件(数据存储文件)分别生成,如果没有该选项,则不会有单独的data文件生成,数据存储和控制都在ctl一个文件 中, 经常用到的还有选项 -i id_colum指定id序列列,默认是id,即作为生成数据的唯一性标志,此列是不可能重复的,一般作为索引列,-s **指定生成srid,默认是null,目前一般是8307吧,-g geometry column指定sdo_geometry,默认是GEOM,此选项一般使用默认即可,命令执行后,会生成三个文件,statearea.sql、 statearea.ctl、stateare.data。
分解完成后就是导入,这里仍然使用命令行的方式
D:\data\>sqlplus pgg/pgg@orcl
SQL>@statearea.sql
SQL>quit
导入ctl文件
D:\data\>sqlldr pgg/pgg@orcl statearea
建立空间索引
D:\data\>sqlplus pgg/pgg@orcl
SQL>CREATE Index STATEAREA_idx ON STATEAREA (GEOM) INDEXTYPE is MDSYS.SPATIAL_INDEX;
至此空间数据导入完毕。
另外还需要建立一个视图,具体作用说不大清楚,但在启动程序发布地图的时候会用到,如果没有此视图文件会出现图层无法map的错误,建立过程如下:
SQL> CREATE OR REPLACE VIEW CS_SRS AS
2 SELECT "CS_NAME","SRID","AUTH_SRID","AUTH_NAME","WKTEXT","CS_BOUNDS"
3 FROM MDSYS.CS_SRS
4 /
另外如果自己手工新建图层的话,过程如下:
DROP TABLE BZXX;
CREATE TABLE BZXX(
GID VARCHAR2(32) NOT NULL,
NAME VARCHAR2(20),
BZNR VARCHAR2(400),
BEIZ VARCHAR2(600),
COOR MDSYS.SDO_GEOMETRY
);
ALTER TABLE BZXX
ADD CONSTRAINT PK_BZXX PRIMARY KEY (GID);
DELETE FROM USER_SDO_GEOM_METADATA
WHERE TABLE_NAME = 'BZXX' AND COLUMN_NAME = 'COOR' ;
INSERT INTO USER_SDO_GEOM_METADATA (TABLE_NAME, COLUMN_NAME, DIMINFO, SRID)
VALUES ('BZXX', 'COOR',
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);
COMMIT;
CREATE INDEX BZXX_GEOM_INDEX ON BZXX(COOR)
INDEXTYPE IS MDSYS.SPATIAL_INDEX
PARAMETERS('SDO_INDX_DIMS=2, LAYER_GTYPE=POINT');