shp导入Oracle

ESRI Shapefile格式的矢量数据导入到OracleSpatial空间数据库，有四种方法：

1.        使用ArcGIS（版本10.0以上）将shp数据导入Oracle（需要安装Oracle客户端或者服务端）；

2.        使用小工具shp2sdo或者ogr2ogr将数据导入Oracle（需要安装Oracle客户端或者服务端）；

3.        代码中使用GDAL中的ogr将shp矢量数据导入Oracle；

4.        SQL语句导入Oracle。

 

第一种方法：使用ArcGIS将shp数据导入Oracle。

导入步骤如下：

如图1所示，在Catalog中添加一个空间数据库的连接。

图1 添加Oracle数据库连接

参数如下所示：

Database Platform：Oracle

Instance：192.168.202.1:1521/orcl

Authentication Type：Databaseauthentication

User Name：username

Password：password

其中，192.168.202.1是连接Oracle的IP地址，1521是连接端口号，orcl是Oracle的实例名。

连接上以后就可以对Oracle中的数据进行导入导出了。

数据导入导出：如图2所示，在连接名上点击右键，选择Import或者Export就可以批量进行矢量数据和表格数据的导入导出了，也可以在数据库中的某一个表上右键，选择Export导出单个数据。

图2 数据导入导出

 

第二种方法：使用小工具shp2sdo.exe或者ogr2ogr.exe将数据导入Oracle。

shp2sdo.exe（只有一个exe文件）下载地址：http://download.csdn.net/detail/zhzhx0318/9619438

ogr2ogr.exe（需添加OCI类库编译后才可用，暂不介绍）。

导入步骤：

(1)    打开cmd，进入shp2sdo.exe的目录（例如D:\shp2oracle\shp2sdo.exe），如图3所示。

图3 cmd操作

(2)    然后输入：shp2sdo D:\DJQ1 DJQ -i id -g geom -s 4610 -d

如图4所示。

图4 生成ctl和sql

其中：

D:\DJQ1是shp文件路径及文件名（此处注意不带后缀名.shp）；

DJQ是导入到Oracle中的表名；

-i id指定id序列的列名，默认是id，即生成数据的唯一性标志；

-s 4610是指生成的SRID（即shp文件的坐标系，可从ArcGIS中查询如图4所示）。

图5 坐标系SRID

-g geom是生成的空间坐标列的列名，默认geom；

-d表示将生成的ctl文件和data文件合并到一个文件。如果没有这个则会生成单独的ctl文件和data文件。

输完命令，回车以后发现，shp2sdo.exe目录下多出来两个文件，DJQ.ctl和DJQ.sql。

(3)    sqlplus连接数据库，执行sql脚本命令

cmd下输入：sqlplususername/password@192.168.202.1:1521/orcl

SQL>@D:\shp2oracle\DJQ.sql

SQL>quit

(4)    导入数据

cmd下输入：sqlldr username/password@192.168.202.1:1521/orcl

control=D:\shp2oracle\DJQ.ctl

(5)    建立空间索引

cmd下输入：sqlplususername/password@192.168.202.1:1521/orcl

SQL>create index DJQ_IDX on DJQ(geom)indextype is mdsys.spatial_index;

(6)    对Oracle地理空间数据进行升级

SQL>EXECUTESDO_MIGRATE.TO_CURRENT(‘DJQ’,’GEOM’);

如果没有这一步，geoserver发布图层后，Tomcat后台会报错，而且图层也不能正常显示。

 

第三种方法：使用GDAL中的ogr将shp矢量数据导入Oracle

此方法需要编译好的GDAL（编译的过程中需要把OCI相关的类库也编译进去）。

源码（包括整个解决方案）下载地址：http://download.csdn.net/detail/zhzhx0318/9619447

备注：不用安装Oracle客户端或者服务器端，源码中的项目包括GDAL相关的dll，源码中包括第三种方法和如下的第四种方法，都有注释说明。

 

第四种方法：SQL语句导入Oracle

导入步骤：

(1)      删除数据库中已存在的表：

DROP TABLE tablename;

(2)      删除表对应的元数据记录

DELETE FROM MDSYS.SDO_GEOM_METADATA_TABLEWHERE SDO_OWNER='username' AND SDO_TABLE_NAME='tablename';

(3)      新建表

CREATE TABLE tablename(

ID NUMBER(38) PRIMARY KEY,

SHAPE MDSYS.SDO_GEOMETRY);

(4)      插入元数据记录

INSERT INTO USER_SDO_GEOM_METADATA(TABLE_NAME, COLUMN_NAME, DIMINFO, SRID)

VALUES ('tablename', 'SHAPE',

MDSYS.SDO_DIM_ARRAY

(MDSYS.SDO_DIM_ELEMENT('X', 122, 124,0.000000050),

MDSYS.SDO_DIM_ELEMENT('Y', 47, 49,0.000000050)

),

4610);

(5)      插入数据

insert into tablename values(‘1’,sdo_geometry(Polygon((100.00 40.00, 150.00 40.00, 150.00 70.00, 100.00 70.00,100.00 40.00)),null));

备注：示例只插入了一条记录，可以插入多条记录

(6)      更新数据的SRID

update tablename t sett.SHAPE.SDO_SRID=4610;

(7)      创建索引

create index tablename on tablename (SHAPE)indextype is MDSYS.SPATIAL_INDEX;

 

总结：

第一种ArcGIS导入的方法最方便，而且导入的效率也比较高，缺点是导入的过程中容易出错，而且需要安装Oracle客户端或者服务器。

第二种使用工具的方法步骤比较麻烦，但是可以写到代码里边，在代码里调用cmd上传shp文件到Oracle。

第三种使用GDAL的方法效率也比较高，但有个最大的缺点，导入数据后投影不一样，导致发布会有问题。但是如果只是数据库导入导出的话，是没有问题的。

第四种使用SQL语句的方法效率较低，但是都可以在代码里直接实现，坐标信息可以使用GDAL读取。