ApacheSedona的ST_Buffer和ST_Distance之间为什么会出现意外差异?
创始人
2024-09-06 10:31:16
0

这可能是因为两者的EPSG(欧洲石油地理信息标准)编号不同导致的。ST_Buffer默认使用WGS84(EPSG:4326)坐标系,而ST_Distance默认使用平面坐标系(EPSG:3857)。如果两个空间对象在不同的坐标系中,那么它们的距离和缓冲结果将会有所不同。

以下是使用ST_Buffer和ST_Distance的示例代码:

ST_Buffer:

from pyspark.sql.functions import col
from sedona.register import SedonaRegistrator
from sedona.utils import SedonaKryoRegistrator
from sedona.core.enums import GridType
from sedona.core.SpatialRDD import PointRDD 
from sedona.core.formatMapper.shapefileParser import ShapefileReader
from sedona.core.formatMapper import GeoJsonReader, WktReader, WkbReader
from sedona.sql.utils import SedonaSQLRegistrator
from sedona.core.spatialOperator import JoinQuery
from sedona.core.enums import JoinBuildSide, DistanceFunction

SedonaRegistrator.registerAll(spark)
SedonaSQLRegistrator.registerAll(spark)

input_location = "data/shapefiles/parks.shp"
shapefile_rdd = ShapefileReader.readToGeometryRDD(sc, input_location)
shapefile_rdd.analyze()

query_point = WktReader.read("POINT (-120.139, 38.788)")

result = JoinQuery.SpatialBufferQuery(
        shapefile_rdd,
        query_point,
        False,
        True
    )

result.show()

ST_Distance:

from pyspark.sql.functions import col
from sedona.register import SedonaRegistrator
from sedona.utils import SedonaKryoRegistrator
from sedona.core.enums import GridType
from sedona.core.SpatialRDD import PointRDD 
from sedona.core.formatMapper.shapefileParser import ShapefileReader
from sedona.core.formatMapper import GeoJsonReader, WktReader, WkbReader
from sedona.sql.utils import SedonaSQLRegistrator
from sedona.core.spatialOperator import JoinQuery
from sedona.core.enums import JoinBuildSide, DistanceFunction

SedonaRegistrator.registerAll(spark)
SedonaSQLRegistrator.registerAll(spark)

input_location = "data/shapefiles/parks.shp"
shapefile_rdd = ShapefileReader.readToGeometryRDD(sc, input_location)
shapefile_rdd.analyze()

query_point = WktReader.read("POINT (-120.139, 38.788)")

result = JoinQuery

相关内容

热门资讯

前端-session、jwt 目录:   (1)session (2&#x...
linux入门---制作进度条 了解缓冲区 我们首先来看看下面的操作: 我们首先创建了一个文件并在这个文件里面添加了...
关于测试,我发现了哪些新大陆 关于测试 平常也只是听说过一些关于测试的术语,但并没有使用过测试工具。偶然看到编程老师...
前缀和与对数器与二分法 1. 前缀和 假设有一个数组,我们想大量频繁的去访问L到R这个区间的和,...
nodejs:本地安装nvm实... 一、背景-使用不同版本node的原因 vue3+ts、nuxt3版本,node...
JAVA集合知识整理 Java集合知识整理 HashMap相关 HashMap的底层数据结构:jdk1.8之...
无刷直流电机介绍及单片机控制实... 无刷直流电机介绍及单片机控制实例前言基本概念优势与劣势使用寿命基本结构使用单片机控制实例电子调速器&...
fwdiary(2) dp2 1.传纸条  AcWing 275. 传纸条 - AcWing 走两条路,走一条最大的...
常用的DOS命令 常用的DOS命令 DOS(Disk Operating System,磁...
<C++> 类和对象(下) 1.const成员函数将const修饰的“成员函数”称之为const成员函数,cons...