Question

我有以下查询：

select distance(GeomFromText('POINT(8 49)',4326),GeomFromText('LINESTRING(8.329969 49.919323,8.330181 49.919468)',4326))

这给了我0.97度。但我需要它以米为单位，并且不知道要转换为哪个SRID。

有人能举例说明如何以米为单位获得结果吗？

这些职位都在欧洲。

Answer 1

只需将度数乘以111195 - 此值为(Earth mean radius)*PI/180 - 即“地球表面上一个大圆度的平均长度”。

使用此方法获得的结果在WGS84椭圆体的测地距离的1％范围内。

修改

好吧，我上面的回答仍然代表一个问题：“如何将弧度转换为以米为单位的长度”，但是，这不是你问的问题（应该问）。

我没有专业地使用Spatialite，所以我假设您的示例查询确实返回'以度为单位的长度'。这不是真的。

不幸的是，似乎Spatialite无法计算“地理意义上的距离”。尽管您的几何形状是使用SRID 4326定义的，但它将它们视为在平面上。

这是一个简单的证明：

select Distance(GeomFromText('POINT(0 0)',4326),GeomFromText('POINT(3 4)',4326));

返回5.0。

这很遗憾......

让我们看看您的原始查询：

select Distance(
  GeomFromText('POINT(8 49)',4326),
  GeomFromText('LINESTRING(8.329969 49.919323,8.330181 49.919468)',4326)
)

MS SQL Server中的等效查询：

SELECT (geography::STGeomFromText('POINT(8 49)', 4326)).STDistance(geography::STGeomFromText('LINESTRING(8.329969 49.919323,8.330181 49.919468)', 4326));

立即为您提供正确的结果：105006.59673084648，以米为单位，无任何额外的骚动。

那么您对Spatialite的选择是什么？

事实上，正如您在评论中所说，一种选择是投射几何图形，并计算它们。在欧洲使用SRID 3035也是有意义的（如果您的位置主要在德国，我会考虑SRID 25832）。

select Distance(
  Transform(GeomFromText('POINT(8 49)',4326),25832),
  Transform(GeomFromText('LINESTRING(8.329969 49.919323,8.330181 49.919468)',4326),25832)
)

返回104969.401605453。

关于你的其他样本（在评论中）：

select distance(
  Transform(GeomFromText('POINT(8.328957 49.920900)',4326),3035),
  Transform(GeomFromText('POINT(8.339665 49.918000)',4326),3035)
)

有一种更简单的方法（如果你有两个POINT，而不是POINT和LINESTRING）：使用你的POINT创建一个LINESTRING并使用GeodesicLength函数，如下所示：

select GeodesicLength(GeomFromText('LINESTRING(8.328957 49.920900, 8.339665 49.918000)',4326))

按预期返回833.910006698673。

Answer 2

在SpatiaLite的functions reference guide中，您可以看到Distance()函数有两个版本。一个只接受两个参数并以CRS单位返回距离，另一个接受3个参数并以米为单位返回距离。

要以米为单位获取距离，只需将第三个参数传递给距离：

sqlite> select Distance(MakePoint(0, 0), MakePoint(3, 4));
5.0
sqlite> select Distance(MakePoint(0, 0), MakePoint(3, 4), 1);
554058.923752633

在Spatialite中获取米的距离而不是度数

2 个答案: