4

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

有些人建议使用诸如 ogr2ogr 之类的外部库，而这对我来说不是一个好选择，因为我需要即时进行。我在这里找到了一篇文章 Transform/Project a geometry from one SRID to another，但是这个转换为 lat 和 long 并且具有一些有限的能力。

我正在寻找一种解决方案，允许我在任何 srid 中进行重新投影，这与ST_Transformpostgresql (postgis) 中的非常相似

4

1 回答 1

6

好吧，如果您正在寻找一个简单的选择，恐怕这是不可能的。但是有一个解决这个问题的方法。前段时间我有一个类似的问题，我必须在.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# 库并将以下代码粘贴到那里（此代码获取每一行数据，遍历每个元素，将其分解，提取每个经纬度，重新投影并替换它）

 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 复制到一个新文件夹中，如果您愿意，可以将它们留在调试文件夹中（其中一些 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，因此您需要从某处获取等效的 srid）。
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

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

在此处输入图像描述

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

于 2018-10-03T01:33:57.600 回答