如何使用Spark/Geomesa将几何图形从一个EPSG投影到另一个EPSG？
Spark是一个开源的大数据处理框架，而Geomesa是一个基于Spark的地理空间数据处理库。它们可以用于将几何图形从一个EPSG投影到另一个EPSG。
要使用Spark/Geomesa进行几何图形的投影转换，可以按照以下步骤进行：
- 导入必要的库和依赖项：
import org.apache.spark.sql.SparkSession
import org.locationtech.geomesa.spark.jts._- 创建SparkSession：
val spark = SparkSession.builder()
.appName("Geometry Projection")
.master("local[*]") // 根据实际情况选择合适的master
.config("spark.sql.extensions", "org.locationtech.jts.spark.sql.CatalystSimpleFeatureEncoders")
.getOrCreate()- 读取包含几何图形的数据集：
val inputDF = spark.read
.format("geomesa")
.option("geomesa.feature", "your_feature_name")
.option("geomesa.converter", "your_converter_name")
.option("geomesa.query", "your_query")
.load()其中，your_feature_name是数据集的特征名称，your_converter_name是数据集的转换器名称，your_query是数据集的查询条件。
- 定义源EPSG和目标EPSG：
val sourceEPSG = "EPSG:4326" // 源EPSG
val targetEPSG = "EPSG:3857" // 目标EPSG- 进行几何图形的投影转换：
val projectedDF = inputDF.withColumn("projected_geometry", st_transform($"geometry", targetEPSG))其中，geometry是包含几何图形的列名，projected_geometry是转换后的几何图形列名。
- 查看转换结果：
projectedDF.show()以上步骤中，st_transform函数用于执行几何图形的投影转换。
关于Spark和Geomesa的更多详细信息和使用方法，可以参考腾讯云的相关产品和文档：
请注意，以上答案仅供参考，实际操作可能会因环境和数据集的不同而有所差异。
相关·内容
我正在“翻译”一些Postgis代码到Geomesa，我有一些Postgis代码，如下所示： select ST_Transform(ST_SetSRID(ST_Point(longitude, latitude), 4326), 27700) 例如，它将点几何图形从4326转换为27700。在Geomesa-Spark-sql文档https://www.geomesa.org/documentation/user/spark</
浏览 55提问于2020-06-24得票数 0
在Python3.7中，我尝试沿着一个公共属性(ZCTA Zip Code number)连接两个文件，并输出一个shapefile，其中包含沿着ZCTA的连接，并添加属性Manufa_EMP (制造雇佣)，然后更改投影的坐标系(PCS或CRS)。但是，由以下代码输出的shapefile会干扰输出文件的地理投影的比例和大小。当我在ArcMap中将shapfile作为地图上的图层打开，并将其与第一个导入的shapefile进行比较时，输出的shapfile的投影要小得多。
浏览 57提问于2020-06-14得票数 1
我将一个具有lat/lon的数据集转换为geopandas dataframe (让我们称之为gpd)。gpd没有CRS，我试图使用以下方法将这个gpd投影到"EPSG:3857“：
我说“无法转换天真的几何图形。这是否意味着对于任何没有CRS的gpd，我必须首先为其分配一个CRS，然后将其重新投影</
浏览 4提问于2020-04-05得票数 2
回答已采纳
1回答
这些地图(A)基于EPSG26912投影。现在我想使用这些作为我的openlayer web应用程序的底图；此外，我还想使用其他基础地图源(B)，例如OpenStreetmap，谷歌地图，ESRI等。由于这些地图中的大多数都是EPSG900913格式的，据我所知这是sphericalMercator投影，如果我只是将这些基础地图图层添加到openlayer地图组件中，A组和B组将不会正确叠加，这意味着它们将显示为单独的地图我知道这是与投影</e
浏览 0提问于2011-08-09得票数 2
1回答
我使用以下代码检索多边形的WKT (知名文本)格式: getGeometry().toString()；"POLYGON((2214109.1963538 -3576110.9131445,2150513.5888293我试过改变地图的投影，但没有结果。
浏览 3提问于2016-09-30得票数 0
回答已采纳
我正在使用python bokeh在世界地图上投影数据。我有一个带有crs="epsg:4326“的geopandas geopdwd init，我想用bokeh tile把它投影到墨卡托的世界地图上。然后我填充一些颜色选择的国家。基本上，我需要将我的geopandas几何图形(一个epsg:4326)转换为epsg=3857。为此，我这样做： geopdwd = geopdwd.to_crs(ep
浏览 61提问于2020-12-30得票数 0
回答已采纳
我被我认为是不知道如何为我的州平面(阿拉斯加州立飞机第4区NAD83英尺)找到正确的投影输入来使用function sp::CRS的问题所困扰。我查了一下?我正在使用Hadley Wickham的教程绘制shapefile，网址是：我的代码：require(maptools) aes(long,lat,group=group) +
geom_path(co
浏览 3提问于2014-02-21得票数 0
1回答
:4326"})#gdf.to_crs(10000) 4.将几何图形列更改为新的buffer_radius表列 gdf = gdf.set_geometry('buffer_radius')
gdf.geometry.name Out: 'buffer_radius' 上面的过
浏览 48提问于2020-10-05得票数 3
回答已采纳
1回答
我有一个形状文件(.shp)，其中包含EPSG:2180中的数据。我提取了像这样的位置：与我在等效.gml文件中看到的值相同。，但是我如何把它转换成纬度和经度？(来自注释GPS是WGS 84 (a.k.a )。第4326条)。我看到了关于EPSG:2180的参数central_meridian=19false_easting=500000
浏览 4提问于2020-06-05得票数 1
1回答
PARAMETER["Latitude_Of_Origin",23.0], ]你能解释一下所有这些参数的含义吗？谢谢!
浏览 1提问于2013-08-18得票数 0
回答已采纳
[-73.979681,40.7033127],
center: ol.proj.transform([-73.979681,40.7033127], 'EPSG:4326', 'EPSG:3857'),
});
我在上找到了如何将投影从一个转换到另一个。然而，在另一个教程中，这一点并不十分清楚。我不习惯那
浏览 4提问于2015-04-27得票数 0
回答已采纳
我正在构建一个使用Angular和mapbox-gl-js作为地图库的web应用程序。地图瓦片使用EPSG:4326plate-carrée投影，数据集为GeoJson (WGS84)。由于mapbox仅支持web墨卡托投影，因此数据不会显示在地图上的正确位置。有什么建议如何重新投影数据集，以便mapbox能够在正确的位置显示数据？谢谢。
浏览 3提问于2019-02-03得票数 0
我目前使用的是srsName，但是这些层的位置是错误的，如果我是对的，那么问题来自于，但是如何更改呢？现在我有了同样的问题，它给了我以下结果：
几何图形应该在这个国家的底部。我认为VectorSource有一个默认的srsname，但我不知道如何更改它。我使用的是"EPSG:2154“，但是地图和这个TileWMS完全匹配，我不知道怎么做到的。My VectorSource有格式-> dataProjection ->代码："<em
浏览 1提问于2020-02-13得票数 0
1回答
我有一个名为eq的包含地震数据的dataframe，其中包含列出纬度和经度的列。使用geopandas，我创建了一个包含以下内容的点列：from shapely.geometry importPoint(x,y) for x, y in zip(df['longitude'], df['latitude'])])eq.crs = {'ini
浏览 11提问于2017-03-13得票数 3
如何得到以米为单位的半径可能是存在的问题，但我没有得到正确的结果。我试图在postgis中创建具有相同半径和圆心的多边形，这些半径和圆心来自openlayers圆。var radiusInMeters = circleRadius * ol.proj.METERS_PER_UNIT['m'];
在得到中心，半径(以米为单位)后，我尝试使用postgis (服务器作业基本上，我对Circle的输入/输出将仅以米为单位。
浏览 2提问于2015-08-25得票数 5
回答已采纳
我正在编程一个应用程序使用Proj4和OpenLayers来显示来自荷兰的地图从外部(开放)的世界管理系统服务器。我已经让它工作了，但重新投影仍然不太正确。现在，通过在线阅读，我发现问题在于EPSG:28992 Proj4字符串将坐标重定向到与我们的视图不同的系统。我们的视图设置为EPSG:3857 (我们的大多数层都是)，而另一层显然重定向到EPSG:4326。有没有办法在层的投影级‘插入’(当定义来自源的
浏览 25提问于2017-02-23得票数 0
我已经创建了一个RDD，它包含特定区域的DGM光栅，并具有crs : EPSG 32632。我想将这些数据覆盖在一张图像上，以便进一步分析，该图像具有crs:EPSG 4326，这也是一种RDD。虽然这两个图像在QGIS中彼此重叠，但投影系统和范围是不同的，这就是为什么后来当我尝试将两个RDD重新投影到webmercator时，它给我一个黑色的图像。val (_,inputrdd) : (Int, Multiband
浏览 29提问于2019-05-06得票数 0
回答已采纳
我试图用默认的数据集制作一个与Geopandas交互的地图。countries.to_crs(epsg=3395)现在，我得到以下错误： 1 #Popultion mapping- Interactive
2 countries.to_crs(epsg
浏览 9提问于2022-02-13得票数 0
回答已采纳
1回答
我正在尝试将shapefile绘制为点，这是一张Folium地图，但我费了很大力气才能得到正确的投影。数据在EPSG:28992系统中，RD新坐标系。这是一个投影坐标系，因此它提供了以米为单位的坐标。要在folium地图上绘制，我需要将其转换为EPSG:3857，如果我正确的话，它可以与度数一起使用。它还具有正确的CRS，如print(WTG.crs)的输出所示 epsg:28992 现在，为了将其转换为其他CRS，我使用以下简单的代
浏览 60提问于2020-06-18得票数 0
回答已采纳
如何使用geotools或其他java库将我的经纬度转换为EPSG 3857投影？我很难找到合适的方法来使用。lon, lat: -71 1.25820, 42 22.0932final CoordinateReferenceSystem source = CRS.decode( "EPSG:4236" );
final CoordinateReferenceSystem dest = CRS.decode("EPSG<
浏览 10提问于2013-03-25得票数 4
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM活动推荐
腾讯云开发者
