便携Gis数据库spatialite使用及与Arcgis数据互转

最近看到有些程序员研究在手机上或者平板上研究基于spatialite的矢量数据处理，或者称之为离线数据绘制或者编辑的一种实现。

由于我也在学习研究android的开发，所以就对基于sqlite数据库的这个spatialite学习了一番。

 

spatialite，如果需要对其详细了解请进下面网站

 

http://www.gaia-gis.it/gaia-sins/  (就是spatialite项目的站点了)

 

spatialite是sqlite的一个扩展，使用相对来说并不复杂，按照网站上的Cook book，如下操作建立支持空间数据的sqlite库

1、创建表，这个没什么说的。就是sqlite创表了。

CREATE TABLE test_geom (  id INTEGER NOT NULL    PRIMARY KEY AUTOINCREMENT,  name TEXT NOT NULL,  measured_value DOUBLE NOT NULL); 

2、给表增加Geometry列，就是一个BLOB类型的图形数据列。

SELECT AddGeometryColumn('test_geom', 'the_geom',  4326, 'POINT', 'XY', 1);

参数简单解释一下AddGeometrt(

'表名'

'要加的图形列的列名'

坐标SRID对应的数字

'图形类型',

'XY'  这一列是标示二维三维等 如果是三维就是写'XYZ'

数字1标示该列不为null相反0则可为null

 

其中图形类型可以分为POINT  LINESTRING  POLYGON MULTIPOINT  MULTILINESTRING  MULTIPOLYGON  GEOMETRYCOLLECTION  GEOMETRY

图形类型常用的就是前三种了。

3、删除图形列，删除后可以用RecoverGeometryColumn恢复

SELECT DiscardGeometryColumn('test_geom', 'the_geom');SELECT RecoverGeometryColumn('test_geom', 'the_geom',  4326, 'POINT', 'XY');

 

4、空间索引,spatialite默认是RTree空间索引，建完表以后使用下面语句建立空间索引。

SELECT CreateSpatialIndex('表名','图形列名');    

以上简单介绍了Spatialite库的创建。下面可以使用.net 开发语言看如果调用这些sql 或者进行Spatialite二次开发

首先需要下载Spatialite的dll库,Spatialite需要在Linux环境c++编译，有其他人编译好的dll或者也可以自己按照说明编译

代码下载地址  https://www.gaia-gis.it/fossil/libspatialite/index  网站也提供了已经编译后的库 http://www.gaia-gis.it/spatialite-3.0.0-BETA/binaries.html

 http://download.csdn.net/detail/rex0y/7904053

 

以c#为例调用过程如下

 

 

//注意引用 sqlite的数据库驱动System.Data.SQLite.dllstatic bool ExecuteSql(string sql)        {            SQLiteCommand cmdinsert = new SQLiteCommand(sql, conn);            int r = cmdinsert.ExecuteNonQuery();            if (r > 0)                return true;            else                return false;        }private void InitDb()        {            SQLiteConnection conn = null;/*  打开连接 如果没有数据库文件自动创建*/            string dbpath = "Data Source=" + Application.StartupPath + "\\db.sqlite";            conn = new SQLiteConnection(dbpath);            conn.Open();/*  加载sqlite扩展库 spatialite */            string sql = "SELECT load_extension('libspatialite-4.dll');";            bool bRet = ExecuteSql(sql, conn);/* 初始化空间数据库描述的一些表 */            sql = "SELECT InitSpatialMetaData()";            bRet = ExecuteSql(sql, conn);/*  */            sql = "CREATE TABLEtest_geom (id INTEGER NOT NULL PRIMARY KEY AUTOINCREMENT,name TEXT NOT NULL,measured_value DOUBLE NOT NULL)";            bRet = ExecuteSql(sql, conn);/*  添加图形列 */            sql = "SELECT AddGeometryColumn('test_geom', 'the_geom',4326, 'POINT', 'XY')";            bRet = ExecuteSql(sql, conn);/* 添加空间索引 */            sql = "SELECT CreateSpatialIndex('test_geom', 'the_geom')";            bRet = ExecuteSql(sql, conn);/*  通过OGC标准的WTK 文件描述格式插入一个点记录*/            sql = "INSERT INTO test_geom(id, name, measured_value, the_geom) VALUES (NULL, 'first point', 1.23456,GeomFromText('POINT(1.01 2.02)', 4326))";            bRet = ExecuteSql(sql, conn);/*  通过OGC标准的WTK 文件描述格式插入一个点记录*/            sql = "INSERT INTO test_geom(the_geom, measured_value, name, id)VALUES (GeomFromText('POINT(11.01 11.02)', 4326),11.123456789, 'eleventh point', NULL)";            bRet = ExecuteSql(sql, conn);            conn.Close();        }

 

对于建立好的这个spatialite库可以使用工具 spatialite_gui.exe  打开查询一下

可以看到一些metadata的表，是由spatialite自动建立的，而Spatial Index分组下是空间索引。当然也可以用一般的sqlite工具打开。但可能会有一些警告。因为空间索引是虚表。

 

5、空间数据如果导入spatialite呢，以我们常见的arcgis数据转移到spatial为例，可以将arcgis的图形读取并转换为WKT (文本) 或者WKB (二进制)

过程简单实现如下：

AppendGeometryTaggedText(IGeometry geometry, StringWriter writer)        {            if (geometry == null)            {                throw new NullReferenceException("Geometry为空，无法转写WKT");            }            if (geometry is IPoint)            {                IPoint point = geometry as IPoint;                AppendPointTaggedText(point, writer);            }            if (geometry is IPolyline)            {                IPolyline polyline=geometry as IPolyline;                AppendPolylineTaggedText(polyline, writer);            }            if (geometry is IPolygon&&!ableMULTI)            {                IPolygon polygon = geometry as IPolygon;                AppendPolygonTaggedText(polygon, writer);            }            if (geometry is IPolygon && ableMULTI)            {                IPolygon polygon = geometry as IPolygon;                AppendMultiPolygonTaggedText(polygon, writer);            }                    }