shapefile文件导入 转载

转自：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

 

 

 

如何将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');