GeoLocationData缩放因子返回错误的值

neyton · 发表于 2015-9-26 08:22:10

嗨，
我是一个民用3D用户，我正在尝试实现此代码以将本地坐标（x， y）转换为UTM网格坐标（北，东）：
Dim gd As GeoLocationData=DB.GeoDataObject.GetObject（OpenMode.ForRead）
Dim As Matrix3d=Matrix3d.Displacement（gd.DesignPoint.GetVectorTo（gd.ReferencePoint））
m=m*Matrix3d.Rotation（gd.NorthDirection，gd.UpDirection，gd.DesignPoint）
m=m*Matrix3d.Scaling（gd.ScaleFactor，gd.DesignPoint）
Dim ptUTM as Point3d=ptLOCAL.TransformBy（m）
但是，gd.ScaleFactor始终返回1，即使gd.ScaleEstimationMethod=ScaleEstimationMethod.ScaleEstMethodReferencePoint
注意：民用3D返回正确的值：
Autodesk.Civil.ApplicationServices.CivilApplication.ActiveDocument.Settings.DrawingSettings.TransformationSettings.GridScaleFactor
怎么了？我错过了什么吗？

本帖以下内容被隐藏保护；需要你回复后，才能看到！

游客，如果您要查看本帖隐藏内容请回复

Sheldon1874 · 发表于 2015-10-8 20:25:06

所以我对你用的autoCAD的味道一无所知。我所知道的是地理信息学。首先想到的是:你是否使用了足够的进动。特定投影的比例因子通常非常接近于1，0.997是我本地情况的起点。第二:你不能用一个点来缩放和旋转，你需要向量。就像我说的不确定“Matrix3d。旋转"和" Matrix3d。缩放“工作，但他们必须需要矢量而不是点？

neyton · 发表于 2015-10-28 14:37:59

您好，
计算UTM坐标所需的精度等于传统程序中出现的精度。实际上，我正在编写一个程序，该程序创建包裹的描述性纪念馆，以及谁将在AutoCAD中工作。在许多情况下，用户希望纪念馆出现在批次顶点的坐标所在的位置，但设计通常是在局部坐标中。然后程序必须具有转换能力，但应该为描述显示正确的结果。
参见此程序：http://tbn2net.appspot.com/download?pacid=C3DMEMO&lang=en

Gasty · 发表于 2015-10-31 19:23:52

嗨，
你说的“局部坐标”是什么意思，局部坐标系统可以是任何东西，x，y(单位？)没有地理意义上的全球参考就没有任何意义。
至于转换，假设您在Civil 3D(一种包含Map的产品)中，您可以使用查询从一个系统转换到另一个系统，前提是您有一个定义良好的本地系统。转换到UTM的程序*必须*知道(或者能够从本地系统定义中推断出)至少UTM区号。加斯顿·努内斯

neyton · 发表于 2015-11-4 09:10:59

局部坐标被仲裁。为了计算一个属性的面积，我们使用局部坐标，因为没有失真。地理参考以UTM坐标或甚至地理坐标给出。但是该地产的面积不能用UTM坐标来确定，因为存在失真，所以您需要转换
通常，它会转换为地心系统，获取顶点的地心坐标，然后计算面积。
geolocation data类将地理坐标转换为本地坐标(autocad wcs ),反之亦然，但不会转换为UTM坐标
为此，我实现了一个代码:

copy


    Function TransformToUTM(gd As GeoLocationData, ByVal local As Point3d) As Point3d
        Dim geo As Point3d = gd.TransformToLonLatAlt(local)
        Dim latd = geo.Y
        Dim lngd = geo.X
        Dim drad = Math.PI / 180.0
        Dim xDoc As New XmlDocument
        xDoc.LoadXml(gd.CoordinateSystem)
        Dim SemiMajorAxis As Double = CDbl(xDoc.GetElementsByTagName("SemiMajorAxis").Item(0).InnerText)
        Dim b As Double = CDbl(xDoc.GetElementsByTagName("SemiMinorAxis").Item(0).InnerText)
        Dim k0 = 0.9996
        Dim e = Math.Sqrt(1 - (b / SemiMajorAxis) * (b / SemiMajorAxis))
        Dim phi = latd * drad
        Dim lng = lngd * drad
        Dim utmz = 1 + Math.Floor((lngd + 180) / 6)
        Dim zcm = 3 + 6 * (utmz - 1) - 180
        Dim e0 = e / Math.Sqrt(1 - e * e) 'Called e prime In reference
        Dim esq = (1 - (b / SemiMajorAxis) * (b / SemiMajorAxis)) ' e squared For use In expansions
        Dim e0sq = e * e / (1 - e * e) ' e0 squared - always even powers
        Dim N = SemiMajorAxis / Math.Sqrt(1 - Math.Pow(e * Math.Sin(phi), 2))
        Dim T = Math.Pow(Math.Tan(phi), 2)
        Dim C = e0sq * Math.Pow(Math.Cos(phi), 2)
        Dim A = (lngd - zcm) * drad * Math.Cos(phi)
        Dim M0 = 0
        Dim M = phi * (1 - esq * (1 / 4 + esq * (3 / 64 + 5 * esq / 256)))
        M = M - Math.Sin(2 * phi) * (esq * (3 / 8 + esq * (3 / 32 + 45 * esq / 1024)))
        M = M + Math.Sin(4 * phi) * (esq * esq * (15 / 256 + esq * 45 / 1024))
        M = M - Math.Sin(6 * phi) * (esq * esq * esq * (35 / 3072))
        M = M * SemiMajorAxis 'Arc length along standard meridian
        Dim x = k0 * N * A * (1 + A * A * ((1 - T + C) / 6 + A * A * (5 - 18 * T + T * T + 72 * C - 58 * e0sq) / 120)) 'Easting relative To CM
        x = x + 500000 'Easting standard 
        Dim y = k0 * (M - M0 + N * Math.Tan(phi) * (A * A * (1 / 2 + A * A * ((5 - T + 9 * C + 4 * C * C) / 24 + A * A * (61 - 58 * T + T * T + 600 * C - 330 * e0sq) / 720)))) 'Northing from equator
        If (y < 0) Then y = 10000000 + y
        Return New Point3d(x, y, local.Z)
    End Function

neyton · 发表于 2015-11-4 09:17:58

该函数接收一个3dpoint作为局部坐标（autocad wcs）并转换为UTM

Sheldon1874 · 发表于 2015-12-8 20:27:15

这是一次有趣的讨论
当我听到“局部网格坐标”时，我首先想到的是基于任意原点的规划网格，类似于10002000100，选择了北向以简化设计（在建筑广场的前面）。可能来自真实世界的测量。然后，这将为网格移动到新网格生成向量。

GeoLocationData缩放因子返回错误的值

发表回复

精选推荐