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/[email protected]:1521/orcl
SQL>@D:\shp2oracle\DJQ.sql
SQL>quit
(4) 导入数据
cmd下输入:sqlldr username/[email protected]:1521/orcl
control=D:\shp2oracle\DJQ.ctl
(5) 建立空间索引
cmd下输入:sqlplususername/[email protected]: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读取。