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读取。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值