Question

我正在尝试在SQL Server中进行一些重新投影，但是不幸的是我找不到任何解决方案。我知道这不是在SQL Server中实现的，所以想知道是否有任何快捷方式。

有人建议使用外部库（例如ogr2ogr），但这对我来说不是一个好选择，因为我需要即时进行操作。我在这里找到了一个帖子 Transform/ Project a geometry from one SRID to another，但此转换为lat和long且功能有限。

我正在寻找一种解决方案，以使我可以在任何狭缝中进行重新投影，这与postgresql（postgis）中的ST_Transform非常相似

Answer 1

好吧，如果您正在寻找一个简单的选择，恐怕是不可能的。但是有解决此问题的方法。不久前，我遇到了类似的问题，我必须在.NET中实现CLR，然后将其作为程序集导入到MSSQL中，虽然有点慢，但可以正常工作。输入为：源投影，目标投影和几何体为文本，输出为重新投影的几何体为文本。请从.Net和PostgreSQL彼此相邻（几乎相同）查看此程序集的结果。

输入的是：

     POLYGON ((1755828.2027002387 5944302.7072003055, 1755826.3549002428 5944302.8313003061, 1755825.1724002406 5944285.1574003045, 1755809.5710002393 5944286.4007003047, ......))";// geomobj.line; //"POINT (1736946.0983 5923253.9175)";

因此，解决方案如下：

1-打开.NET并创建ac＃库，并将以下代码粘贴到此处（此代码获取每一行数据，遍历每个元素，将其分解，提取每个long和lat，重新投影，然后替换）< / p>

 using Microsoft.SqlServer.Server;
 using System;
 using System.Collections;
 using System.Collections.Generic;
 using System.Linq;
 using System.Text;
 using System.Threading.Tasks;

public class classval
{
    public string line;
    public string src;
    public string dst;

    public classval(string line, string src, string dst)
    {
        this.line = line;
        this.src = src;
        this.dst = dst;
    }



}
public class RPG_Transform
{
    private static IEnumerable<classval> ConvertedEnumerable(string line, string src, string dst)
    {
        return new List<classval> { new classval(line, src, dst) };
    }

    [SqlFunction(FillRowMethodName = "FillRow")]
    public static IEnumerable STRPG_Transform(string Geometry, string src, string dst)
    {
        return ConvertedEnumerable(Geometry, src, dst);
    }

    private static void FillRow(Object classvalobj, out string Geometry, out string srcprj, out string dstprj)
    {

        classval geomobj = (classval)classvalobj;
        string _geometry = geomobj.line;  
        string proj4_src = geomobj.src;  
        string proj4_dst = geomobj.dst; 
        string _originalgeom = _geometry;
        _geometry = _geometry.Remove(0, _geometry.IndexOf('('));
        // _geometry = _geometry.Replace("(", "[ ");
        // _geometry = _geometry.Replace(")", " ]");
        string[] splitbycomma = _geometry.Split(',');

        foreach (var itembycomma in splitbycomma)
        {

            string tmpitem = itembycomma;
            tmpitem = tmpitem.Replace('(', ' ');
            tmpitem = tmpitem.Replace(')', ' ');
            tmpitem = tmpitem.Trim();
            string[] splitbyspace = tmpitem.Split(' ');
            for (int ibs = 0; ibs < splitbyspace.Length - 1; ibs++)
            {
                string originallonglat = splitbyspace[ibs] + " " + splitbyspace[ibs + 1];

                double[] yxval = new double[] { double.Parse(splitbyspace[ibs]), double.Parse(splitbyspace[ibs + 1]) };

                double[] z = new double[1] { 0 };




                DotSpatial.Projections.ProjectionInfo src =
                    DotSpatial.Projections.ProjectionInfo.FromProj4String(proj4_src);
                DotSpatial.Projections.ProjectionInfo trg =
                    DotSpatial.Projections.ProjectionInfo.FromProj4String(proj4_dst);

                DotSpatial.Projections.Reproject.ReprojectPoints(yxval, z, src, trg, 0, 1);

                string longlat = yxval[0] + " " + yxval[1];

                _originalgeom = _originalgeom.Replace(originallonglat, longlat);
            }
        }
        srcprj = proj4_src;
        dstprj = proj4_dst;
        Geometry = _originalgeom;
    }

}

2-运行此程序之前，您需要安装一些库，最简单的是使用Nugget Manager，如下所示。

3-接下来，您需要编译代码并获取dll，然后将这些dll复制到新文件夹中（如果您希望将其保留在debug文件夹中）（其中一些dll已添加到项目中取决于依赖性），如下所示。

4-在MSSQL中启用CLR，然后将程序集添加到MSSQL服务器（通过创建如下所示的函数）

  sp_configure 'show advanced options', 1;  
 GO  
 RECONFIGURE;  
 GO  
 sp_configure 'clr enabled', 1;  
 GO  
 RECONFIGURE;  
 GO  

 ALTER DATABASE Prod SET trustworthy ON
  CREATE ASSEMBLY CLRFunctionAssem
  FROM N'E:\CLR\RPG_Transform\RPG_Transform\bin\Debug\RPG_Transform.dll'
  WITH PERMISSION_SET = UNSAFE
  GO

--DROP ASSEMBLY CLRFunctionAssem ---if you need to drop it in the future

   CREATE FUNCTION dbo.RPG_STTransform(@Geometry nvarchar(max), @src 
  nvarchar(max),@dst nvarchar(max))
    RETURNS TABLE
   ( _geom  nvarchar(max) ,srcprj  nvarchar(max) ,dstprj  nvarchar(max) 
  ) with execute as caller
  AS
   EXTERNAL NAME CLRFunctionAssem.[RPG_Transform].STRPG_Transform

5-获取源和目标投影的公式。您需要从某个地方获取这两个公式。对我来说，要获取源投影和目标投影，请使用以下网站。（此解决方案使用Proj4j，因此您需要从某处获取等价物。）
https://epsg.io/ 如下所示：

与目的地相同

6-运行sql函数，如下所示：

         SELECT   top 1000 pk   , geom.STAsText() as input,conv._geom as ouput, 
    geometry::STGeomFromText(conv._geom,4326)
  FROM [Prod].[dbo].DPO_Geographic_Units as a  CROSS APPLY dbo.RPG_STTransform (
  a.Geom.STAsText()--geometry as text
 ,'+proj=tmerc +lat_0=0 +lon_0=173 +k=0.9996 +x_0=1600000 +y_0=10000000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs ',--source =2193
 '+proj=longlat +datum=WGS84 +no_defs'--destination =4326
 ) as conv
 where [Geom] is not null

以及输入和输出彼此相邻：

如果愿意，可以将这些投影公式保存在代码中，并使用切换盒为每个投影获取对应的公式，然后可以传递一个数字（例如2193），然后代码查找等效的公式

。

如何在SQL Server中进行空间重新投影

1 个答案: