網(wǎng)上有個工具可以做到，可惜導入的中文是亂碼，搜索了一圈好像沒什么解決辦法，于是想自己搞一個，哎，一搞才知道麻煩重重啊~

首先查資料，了解了一下SQLServer中的空間數(shù)據(jù)分兩種，Geometry和Geography，一種是平面一種球面，比如在做距離計算時結(jié)果是不一樣的啊

所以我必須先Geography球面坐標系，但這玩意兒有限制啊，比如多邊形外環(huán)方向必須是逆時針，內(nèi)環(huán)方向必須是順時針，不能超過一個半球等等

不管了先做，加載GDAL和ORG，打開MapInfo圖層，遍歷之，取字段正常，取幾何個數(shù)正確，下一步將幾何數(shù)據(jù)導出，方法很多，支持OpenGIS的系統(tǒng)

都支持WKT和WKB，為了方便選ExportToWkt，然后生成 INSERT語句，差不多了，執(zhí)行之，問題來了：Sqlserver報錯,說什么

消息 6522，級別 16，狀態(tài) 1，第 1 行

在執(zhí)行用戶定義例程或聚合 "geography" 期間出現(xiàn) .NET Framework 錯誤:

Microsoft.SqlServer.Types.GLArgumentException: 24205: 因為指定的輸入超過了一個半球，所以它不代表有效的地理實例。每個地理實例必須能容納在一個半球內(nèi)。此錯誤的常見原因是某個多邊形的環(huán)方向錯誤。

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，沒問題可以正常生成幾何，哎，說明真的環(huán)方向有問題，改吧，于是有了以后代碼

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

{

///

/// 判斷環(huán)的是否為順時針方向

///

///

/// True順時針，F(xiàn)alse逆時針

private static bool IsRingClockwise(IList 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;

}

//判斷是否是環(huán)

public static bool IsRing(IList 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 pts, bool resverse)

{

var sl = new List();

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規(guī)則，多邊形外環(huán)應該逆時針，內(nèi)環(huán)順時針

//POLYGON(())中第一個()中的是外形，其它的是內(nèi)環(huán)

//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();

foreach (var p in pa)

{

var pt = Point.Parse(p);

if (pt != null)

ls.Add(pt);

}

//TODO 多邊形有內(nèi)環(huán)的不能逆轉(zhuǎn)，要順時針

if (IsRing(ls))

{

var rs = string.Join(",", Resverse(ls, IsRingClockwise(ls)));//如果是順時針則逆轉(zhuǎn)

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，再導面縣界，結(jié)果一些搞不定，導出WKT用SELECT在SQLER里試了一下，真不行啊，查了一下用網(wǎng)上工具導的記錄，WKT是GEOMETRYCOLLECTION，而我用GDAL/OGR查出來說是POLYGON類型

利用GDAL/OGR導空間數(shù)據(jù)到SQLServer2008 三聯(lián)

這明明是一條線和一個多邊形組成的啊，GDAL/OGR居然說是POLYGON，試了GDAL1.6、GDAL1.7，應該是OGR的問題，不搞了，累覺不愛啊

總之導ESRI Shape用Geometry模式還是可以的，中文問題解決了，配和SQLServer的Update語句基本能解決問題

更多信息請查看IT技術(shù)專欄