利用GDAL/OGR导空间数据到SQLServer2008
2014-04-12 20:47
316 查看
网上有个工具可以做到,可惜导入的中文是乱码,搜索了一圈好像没什么解决办法,于是想自己搞一个,哎,一搞才知道麻烦重重啊~
首先查资料,了解了一下SQLServer中的空间数据分两种,Geometry和Geography,一种是平面一种球面,比如在做距离计算时结果是不一样的啊
所以我必须先Geography球面坐标系,但这玩意儿有限制啊,比如多边形外环方向必须是逆时针,内环方向必须是顺时针,不能超过一个半球等等
不管了先做,加载GDAL和ORG,打开MapInfo图层,遍历之,取字段正常,取几何个数正确,下一步将几何数据导出,方法很多,支持OpenGIS的系统
都支持WKT和WKB,为了方便选ExportToWkt,然后生成 INSERT语句,差不多了,执行之,问题来了:Sqlserver报错,说什么
将geography::STGeomFromText换成geometry::STGeomFromText,没问题可以正常生成几何,哎,说明真的环方向有问题,改吧,于是有了以后代码
以上用的.net 2.0所以string[].ToArray这种都不支持,虽然MutiPolygon和GeomertyCollection还是搞不定, 不过一般Polygon是没问题的,再次测试了一个MapInfo Table
(中国面省界)基本OK,再导面县界,结果一些搞不定,导出WKT用SELECT在SQLER里试了一下,真不行啊,查了一下用网上工具导的记录,WKT是GEOMETRYCOLLECTION,而我用GDAL/OGR查出来说是POLYGON类型
这明明是一条线和一个多边形组成的啊,GDAL/OGR居然说是POLYGON,试了GDAL1.6、GDAL1.7,应该是OGR的问题,不搞了,累觉不爱啊
总之导ESRI Shape用Geometry模式还是可以的,中文问题解决了,配和SQLServer的Update语句基本能解决问题
以下给出源码链接:点击下载
首先查资料,了解了一下SQLServer中的空间数据分两种,Geometry和Geography,一种是平面一种球面,比如在做距离计算时结果是不一样的啊
所以我必须先Geography球面坐标系,但这玩意儿有限制啊,比如多边形外环方向必须是逆时针,内环方向必须是顺时针,不能超过一个半球等等
不管了先做,加载GDAL和ORG,打开MapInfo图层,遍历之,取字段正常,取几何个数正确,下一步将几何数据导出,方法很多,支持OpenGIS的系统
都支持WKT和WKB,为了方便选ExportToWkt,然后生成 INSERT语句,差不多了,执行之,问题来了:Sqlserver报错,说什么
消息 6522,级别 16,状态 1,第 1 行 在执行用户定义例程或聚合 "geography" 期间出现 .NET Framework 错误: Microsoft.SqlServer.Types.GLArgumentException: 24205: 因为指定的输入超过了一个半球,所以它不代表有效的地理实例。每个地理实例必须能容纳在一个半球内。此错误的常见原因是某个多边形的环方向错误。 Microsoft.SqlServer.Types.GLArgumentException: 在 Microsoft.SqlServer.Types.GLNativeMethods.ThrowExceptionForHr(GL_HResult errorCode) 在 Microsoft.SqlServer.Types.GLNativeMethods.GeodeticIsValid(GeoData g) 在 Microsoft.SqlServer.Types.SqlGeography.IsValidExpensive() 在 Microsoft.SqlServer.Types.SqlGeography.ConstructGeographyFromUserInput(GeoData g, Int32 srid) 在 Microsoft.SqlServer.Types.SqlGeography.GeographyFromText(OpenGisType type, SqlChars taggedText, Int32 srid) 在 Microsoft.SqlServer.Types.SqlGeography.STGeomFromText(SqlChars geometryTaggedText, Int32 srid) 。
将geography::STGeomFromText换成geometry::STGeomFromText,没问题可以正常生成几何,哎,说明真的环方向有问题,改吧,于是有了以后代码
using System; using System.Collections.Generic; using System.Text.RegularExpressions; namespace SqlSpatialTools { public class Point { public double X { get; set; } public double Y { get; set; } public Point(double x, double y) { X = x; Y = y; } public static Point Parse(string pt) { var p = pt.Split(new[] { " " }, StringSplitOptions.RemoveEmptyEntries); return p.Length == 2 ? new Point(Convert.ToDouble(p[0]), Convert.ToDouble(p[1])) : null; } public override string ToString() { return string.Format("{0:F8} {1:F8}", X, Y); } } public static class Wkt { /// <summary> /// 判断环的是否为顺时针方向 /// </summary> /// <param name="pp"></param> /// <returns>True顺时针,False逆时针</returns> private static bool IsRingClockwise(IList<Point> pp) { double t = 0; for (var i = 1; i < pp.Count; i++) t += pp[i - 1].X * pp[i].Y - pp[i].X * pp[i - 1].Y; t += pp[pp.Count - 1].X * pp[0].Y - pp[0].X * pp[pp.Count - 1].Y; return t < 0; } //判断是否是环 public static bool IsRing(IList<Point> pp) { var t = pp.Count - 1; return (t > 2 && pp[0].X == pp[t].X && pp[0].Y == pp[t].Y); } public static string[] Resverse(List<Point> pts, bool resverse) { var sl = new List<string>(); foreach (var pt in pts) { sl.Add(pt.ToString()); } if (resverse) sl.Reverse(); return sl.ToArray(); } public static string ReverseWKT(string wkt) { if (!Regex.IsMatch(wkt, "POLYGON")) return wkt; var retwkt = wkt; var ms = Regex.Matches(retwkt, @"\(\d+.*?\d+\)"); //按geography规则,多边形外环应该逆时针,内环顺时针 //POLYGON(())中第一个()中的是外形,其它的是内环 //LINESTRING()直接判断是否为逆时针即可(判断一下需要需要逆时针) if (ms.Count > 0) { for (var i = ms.Count - 1; i >= 0; i--) { var m = ms[i]; var s = m.ToString().TrimStart(new[] { '(' }).TrimEnd(new[] { ')' }); var pa = s.Split(new[] { "," }, StringSplitOptions.RemoveEmptyEntries); var ls = new List<Point>(); foreach (var p in pa) { var pt = Point.Parse(p); if (pt != null) ls.Add(pt); } //TODO 多边形有内环的不能逆转,要顺时针 if (IsRing(ls)) { var rs = string.Join(",", Resverse(ls, IsRingClockwise(ls)));//如果是顺时针则逆转 retwkt = retwkt.Remove(m.Index, m.Length).Insert(m.Index, string.Format("({0})", rs)); } } return retwkt; } return wkt; } } }
以上用的.net 2.0所以string[].ToArray这种都不支持,虽然MutiPolygon和GeomertyCollection还是搞不定, 不过一般Polygon是没问题的,再次测试了一个MapInfo Table
(中国面省界)基本OK,再导面县界,结果一些搞不定,导出WKT用SELECT在SQLER里试了一下,真不行啊,查了一下用网上工具导的记录,WKT是GEOMETRYCOLLECTION,而我用GDAL/OGR查出来说是POLYGON类型
这明明是一条线和一个多边形组成的啊,GDAL/OGR居然说是POLYGON,试了GDAL1.6、GDAL1.7,应该是OGR的问题,不搞了,累觉不爱啊
总之导ESRI Shape用Geometry模式还是可以的,中文问题解决了,配和SQLServer的Update语句基本能解决问题
以下给出源码链接:点击下载
相关文章推荐
- Python空间数据处理1: GDAL读写遥感图像
- 报表统计(十) 访问数据库 利用Ajax与数据控件空间交互
- GDAL\OGR读取数据示例 C#版本
- [转]开源栅格/矢量空间数据转换库-GDAL/OGR
- 使用GDAL/OGR读取多边形数据
- GDAL python教程(1)——用OGR读写矢量数据
- 利用EXP/IMP进行数据迁移,如何转换表空间操作(完整版)
- SQL Server重新利用数据表的空间
- 利用可传输表空间技术实现数据的高效迁移
- Python空间数据处理3: HDF格式转GDAL支持格式
- Python空间数据处理1: GDAL读写遥感图像
- Oracle利用传输表空间导出导入数据的步骤
- 利用ST_GEOMETRY 从一个图层更新另一图层空间坐标数据
- 利用Python实现分布于各工作空间中的分幅数据合并
- 利用GDAL读取图片,将数据以BIP的方式存在内存中。
- Sqlserver2008数据引擎特性之空间数据预先体验
- Python地理数据处理之GDAL/OGR使用
- 利用游标批量清空数据同时批量查询数据空间和条数
- 操作系统以外的硬盘空间-利用DOS/WINDOWS9X的隐含扇区保护数据 - 企鹅人生 - CSDNBlog
- 改個名字再發一次--本例運用了javascript,html控件,服务器控件,最重要的就是javascript对空间的操作--GridView套GridViewtao--- 利用GridView显示主细表并一次编辑明细表所有数据的例子