sql server Geometry 类型操作 笔记

sqlGeometry 类型为sql server 2008之后的版本 新加的一种CLR扩展数据类型，为广大sql server开发人员存储几何类型及空间运算提供极大的便利，下面说明geometry类型的具体操作

示例SQL语句代码

IF OBJECT_ID ( 'dbo.SpatialTable', 'U' ) IS NOT NULL
DROP TABLE dbo.SpatialTable;
GO

CREATE TABLE SpatialTable
( id int IDENTITY (1,1),
GeomCol1 geometry,
GeomCol2 AS GeomCol1.STAsText() );
GO

INSERT INTO SpatialTable (GeomCol1)
VALUES (geometry::STGeomFromText('LINESTRING (100 100, 20 180, 180 180)', 0));

INSERT INTO SpatialTable (GeomCol1)
VALUES (geometry::STGeomFromText('POLYGON ((0 0, 150 0, 150 150, 0 150, 0 0))', 0));
GO

INSERT INTO SpatialTable (GeomCol1)
VALUES(geometry::STGeomFromText('LINESTRING (116.387112 39.920977,116.385243 39.913063,116.394226 39.917988,116.401772 39.921364,116.41248 39.927893,116.387112 39.920977)', 4326))

　　

几何类型操作

 Database db = DatabaseFactory.CreateDatabase();
 DbCommand dbCmd = db.GetSqlStringCommand("SELECT GeomCol1 FROM SpatialTable WHERE id=4");
 using (IDataReader reader = db.ExecuteReader(dbCmd))
 {
 if (reader.Read())
 {
 SqlGeometry o = reader[] as SqlGeometry;

 }
 }

 SqlGeometry geo = SqlGeometry.Parse("POLYGON ((0 0, 150 0, 150 150, 0 150, 0 0))");
 //DataSet ds = db.ExecuteDataSet(dbCmd);

 //SqlGeometry geo = ds.Tables[0].Rows[0][1] as SqlGeometry;

 //SqlDouble area = geo.STArea();

 dbCmd = db.GetSqlStringCommand(@"INSERT INTO SpatialTable (GeomCol1) values(@p)");

 //构建多边形
 SqlGeometryBuilder sb = new SqlGeometryBuilder();
 sb.SetSrid();
 sb.BeginGeometry(OpenGisGeometryType.Polygon);
 sb.BeginFigure(, );
 sb.AddLine(, );
 sb.AddLine(, );
 sb.AddLine(, );
 sb.AddLine(, );
 sb.EndFigure();
 sb.EndGeometry();

 dbCmd = db.GetSqlStringCommand(string.Format("INSERT INTO SpatialTable (GeomCol1) values(geometry::STGeomFromText('{0}', {1}))"
 , sb.ConstructedGeometry.ToString(), ));
 //db.AddInParameter(dbCmd, "@p", DbType.Binary, sb.ConstructedGeometry.STAsBinary().Buffer);

 int cnt = db.ExecuteNonQuery(dbCmd);

　　

数据库对数据进行过滤操作

 SELECT c_geomCol.MakeValid().STCentroid().STAsText() FROM t_green_point
 WHERE c_geomCol IS NOT NULL
 --AND c_geomCol2<>'POLYGON EMPTY'
 AND c_geomCol.STIsValid()=1

通过这个sql获得系统的坐标系（Sql server中）：Select * from sys.spatial_reference_systems

//MultiPolygon 多个多边形结合处理

 SqlGeometryBuilder sb = new SqlGeometryBuilder(); //构造多个多边形实例
 sb.SetSrid();
 sb.BeginGeometry(OpenGisGeometryType.MultiPolygon);

 sb.BeginGeometry(OpenGisGeometryType.Polygon); //构造第一个多边形
 sb.BeginFigure(, );
 sb.AddLine(, -);
 sb.AddLine(-,-);
 sb.AddLine(-,);
 sb.AddLine(,);
 sb.EndFigure();
 sb.EndGeometry();

 sb.BeginGeometry(OpenGisGeometryType.Polygon); //构造第二个多边形
 sb.BeginFigure(,);
 sb.AddLine(,);
 sb.AddLine(,);
 sb.AddLine(,);
 sb.AddLine(,);
 sb.EndFigure();
 sb.EndGeometry();

 sb.BeginGeometry(OpenGisGeometryType.Polygon); //构造第N个多边形
 sb.BeginFigure(,);
 sb.AddLine(,);
 sb.AddLine(,);
 sb.AddLine(,);
 sb.AddLine(,);
 sb.EndFigure();
 sb.EndGeometry();

 sb.EndGeometry();

 geo = sb.ConstructedGeometry;

 //下面是解析代码
 int numGeometries = geo.STNumGeometries().Value; //有多少个多边形

 SqlGeometry geoN = null;

 for (int i = ; i <= numGeometries; i++)
 {
 geoN = geo.STGeometryN(i);//第几个多边形
 for (int j = ; j <= geoN.STNumPoints(); j++) //转到百度地图多边形最后一个点可以去掉
 {
 point = geoN.STPointN(i);
 Console.WriteLine("第{0}个多边形，第{1}点,X={2},Y={3}", i, j, point.STX.Value, point.STY);
 }

 }

　　

输出结果：

参数化的SQL语句传值

   DbCommand dbCmd = db.GetSqlStringCommand(@"INSERT INTO SpatialTable (GeomCol1) values(@p)");

  db.AddInParameter(dbCmd, "@p", DbType.Binary, sb.ConstructedGeometry.Serialize());

  int cnt = db.ExecuteNonQuery(dbCmd);

　　

判断多边形是否相交

DECLARE @bigGeo geometry= 'POLYGON((0 0, 3 0, 3 3, 0 3,0 0))';
DECLARE @smallGeo geometry='POLYGON((1 1 ,2 1,2 2,1 2,1 1))';
DECLARE @midGeo geometry='POLYGON((0 0, 1.5 0, 1.5 1.5, 0 1.5,0 0))';
DECLARE @Geo3 geometry='POLYGON((2 2, 3 2,3 3,2 3,2 2))'; 

SELECT @bigGeo.STIntersection(@smallGeo).STAsText()

SELECT @midGeo.STIntersection(@smallGeo).STIsEmpty()

SELECT @midGeo.STIntersection(@Geo3).STIsEmpty()

如果相交则结果不为空

参考资料：

 https://msdn.microsoft.com/zh-cn/library/bb964711.aspx   

 https://msdn.microsoft.com/zh-cn/library/ee642119.aspx

 https://msdn.microsoft.com/zh-cn/library/microsoft.sqlserver.types.sqlgeometrybuilder(v=sql.120).aspx

 

       

空间数据类型相关Transact-SQL

         https://msdn.microsoft.com/zh-cn/library/cc280487(v=sql.120).aspx