我有以下查询:
select distance(GeomFromText('POINT(8 49)',4326),GeomFromText('LINESTRING(8.329969 49.919323,8.330181 49.919468)',4326))
这给我0.97度.但我需要它在米,不知道哪个SRID要转换.
有人可以给我一个例子,如何得到结果米为空间的?
职位全部在欧洲.
解决方法
只需将度数乘以111195 – 该值为(地球平均半径)* PI / 180 – 即“地球表面上以米为单位的一个大圆度的平均长度”.
使用此方法的最大误差为〜0.1%
编辑
好的,我上面的答案仍然是一个问题:“如何将弧度以度数转换成米长”,但是这不是你问的问题(应该问的).
我没有专业地使用Spatialite,所以我假设你的样品查询确实返回了“以度为单位的长度”.这不是真的.
不幸的是,Spatialite似乎无法计算“地理意义上的距离”.尽管使用SRID 4326定义了几何形状,但它们将它们视为在飞机上.
这是一个简单的证明:
select Distance(GeomFromText('POINT(0 0)',GeomFromText('POINT(3 4)',4326));
返回5.0.
这是一个耻辱 …
让我们来看看你原来的查询:
select Distance( GeomFromText('POINT(8 49)',4326) )
SELECT (geography::STGeomFromText('POINT(8 49)',4326)).STDistance(geography::STGeomFromText('LINESTRING(8.329969 49.919323,4326));
立即获得正确的结果:105006.59673084648,以米为单位,没有任何额外的brouhaha.
那么你有什么选择Spatialite?
实际上,正如您在评论中所说,一个选项是投影您的几何,并计算这些.对于欧洲,使用SRID 3035也是有道理的(如果你的位置主要在德国,我会考虑SRID 25832).
select Distance( Transform(GeomFromText('POINT(8 49)',25832),Transform(GeomFromText('LINESTRING(8.329969 49.919323,25832) )
返回104969.401605453.
至于您的其他样本(在评论中):
select distance( Transform(GeomFromText('POINT(8.328957 49.920900)',3035),Transform(GeomFromText('POINT(8.339665 49.918000)',3035) )
有一个更简单的方法(如果您有两个POINT,而不是一个POINT和一个LINESTRING):使用您的POINT创建一个LINESTRING并使用GeodesicLength函数,如下所示:
select GeodesicLength(GeomFromText('LINESTRING(8.328957 49.920900,8.339665 49.918000)',4326))
按预期返回833.910006698673.