Question

我正在用C＃编写WinForms应用程序。我需要确保Datarow中的两个Datatable相距不超过100公里。每行在单独的DataColumn中有一个UTM区域，东区和北区。所有坐标都使用相同的数据，但有些具有不同的区域（否则，我会使用pythag数学，因为UTM以米为单位）。到目前为止，我使用以下代码来执行此操作，但似乎我的DbGeography.PointFromText方法工作正常（请参阅代码中的*），就像代码到达时一样使用＆＃39; **＆＃39;的代码行在它的开头，我得到一个错误说＆＃34; 24201：纬度值必须在-90到90度之间＆＃34;。我也试过了：

dtTrap.Rows[i]["TrapGeog"] = DbGeography.PointFromText(pointWellKnownText: "POINT M(" + dtTrap.Rows[i][intEastingIndex].ToString() + " " + dtTrap.Rows[i][intNorthingIndex].ToString() + " " + dtTrap.Rows[i][Zone].ToString() + ")", coordinateSystemId: SRID);

只是让它抱怨＆＃34;没有列＃17＆＃34; （17是我的UTM专区）。

我发现使用这些东西的文档很少......据我所知，我的SRID是正确的（我从this site中取出它们）。我之前使用Lat + Longs做过这种事情，并且效果非常好。我无法找到适合UTM的语法。

using System.Data.Entity.Spatial;
...
DataColumn dcGeog = new DataColumn("TrapGeog", typeof(DbGeography));
dtTrap.Columns.Add(dcGeog);
byte Zone;
Int16 SRID;
for (int i = 0; i < dtTrap.Rows.Count; ++i)
{
    if (dtTrap.Rows[i][intZoneIndex] != null
    && dtTrap.Rows[i][intNorthingIndex] != null
    && dtTrap.Rows[i][intEastingIndex] != null
    && byte.TryParse(dtTrap.Rows[i][intZoneIndex].ToString(), out Zone) == true)
    {
        if (Zone == 15) { SRID = 26915; }
        else if (Zone == 16) { SRID = 26916; }
        else if (Zone == 17) { SRID = 26917; }
        else { SRID = 26918; }
        // shove it in:
        try
        {
            *dtTrap.Rows[i]["TrapGeog"] = DbGeography.PointFromText(pointWellKnownText: "POINT(" + dtTrap.Rows[i][intEastingIndex].ToString() + " " + dtTrap.Rows[i][intNorthingIndex].ToString() + ")", coordinateSystemId: SRID);
        }
        catch (Exception ex)
        {
            if (ex.InnerException != null)
            {
                **MessageBox.Show(ex.InnerException.Message);
            }
            else
            {
                MessageBox.Show(ex.Message);
            }
        }

    }
}
for (int i = 0; i < dtTrap.Rows.Count - 1; ++i)
{
    for (int k = i + 1; k < dtTrap.Rows.Count; ++i)
    {
        DbGeography iTrap = (DbGeography)dtTrap.Rows[i]["TrapGeog"];
        DbGeography kTrap = (DbGeography)dtTrap.Rows[k]["TrapGeog"];
        if (iTrap.Distance(kTrap) > 100000)
        {
            sbErrorsAndWarningsLog.Append(@"Warning:  Line number " + (i + 2).ToString() + " on the Trap spreadsheet has coordinates that are at least 100 km away from row " + (k + 2).ToString() + "'s point.  Please check that these coordinates are correct.").AppendLine();
            boolWarningsFound = true;
            break;
        }
    }
}

Answer 1

为了解决这个问题，我在姊妹网站here上偶然发现了这篇文章。我认为DbGeography基于SQL Server的地理数据类型，DbGeometry同样基于几何数据类型。

一些相关的花絮：

地理数据类型EPSG：4326，以度为单位，这就是您的POINT（257306 142708）失败的原因，因为它超出了（-180,180），（ - 90,90）范围。

...

解决方案......使用几何数据类型，即投影而非地理数据。

所以我将代码更改为：

DataColumn dcGeom = new DataColumn("TrapGeom", typeof(DbGeometry));
...
dtTrap.Rows[i]["TrapGeom"] = DbGeometry.PointFromText(pointWellKnownText: stPoint, coordinateSystemId: SRID);

似乎解析得很好。虽然，这篇文章的最后部分让我担心：

SQL Server坚持认为存在SRID，即使您实际上无法对其进行任何有用的操作，例如将一个坐标系转换为另一个坐标系。如果您不小心将几何图形与不同的SRID混合在一起会让您感到非常困惑，因为您将无法获得任何结果（没有警告，但这是暂时的）。

所以，当我最终做if (Trap1.Distance(Trap2) > 100000)时...当我在不同的SRID中有两个不同的点时，我有一半期望它出现故障。我会对我发现的内容进行测试和评论。

Answer 2

为了更好地调试此错误，我建议您先检查输入数据。

string point = "";
try
{
    string point = "POINT(" + dtTrap.Rows[i][intEastingIndex].ToString() + " " + dtTrap.Rows[i][intNorthingIndex].ToString() + ")";

    dtTrap.Rows[i]["TrapGeog"] = DbGeography.PointFromText(pointWellKnownText: point, coordinateSystemId: SRID);
}
catch (Exception ex)
{
    System.Diagnostics.Debug.WriteLine("correct point? " + point);
    if (ex.InnerException != null)
    {
        MessageBox.Show(ex.InnerException.Message + " at " + point);
    }
    else
    {
        MessageBox.Show(ex.Message);
    }
}

我的猜测是，你的字符串中有一些逗号而不是一个点或一些空格。

在C＃中将UTM坐标解析为DBGeography

2 个答案: