如何将ESRI的Shapefile导入到Oracle数据库中(转)

本文转自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');