如何使用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的更多详细信息和使用方法，可以参考腾讯云的相关产品和文档：

请注意，以上答案仅供参考，实际操作可能会因环境和数据集的不同而有所差异。

相关·内容

Google Earth Engine（GEE）——Error: Projection: The CRS of a map projection could not be parsed.（坐标转换错）

世界地图 EPSG:4326 在世界地图方面，EPSG:4326是比较著名的一个，因为由美国主导的GPS系统就是在用它，它还有一个名气更大的别名叫作WGS84，WGS(World Geodetic System...EPSG:3857 另一个比较知名的编码是EPSG:3857，这也是一张世界地图，目前主要是各大互联网地图公司以它为基准，例如Google地图，Microsoft地图都在用它。...- 知乎解析： Earth Engine 可以将图像重新投影到非 EPSG 坐标参考系统。...Earth Engine 使用 Java GeoTools库进行坐标系转换，该库支持坐标系众所周知的文本 (WKT)格式和一组基础投影。...大家可以点击上面的WKT链接，将代码作为一个变量放入到自己的参数当中，然后进行投影变化 // 加载一个简单的Landsat8影像 var image = ee.Image('LANDSAT/LC08/C01

2351 0

如何利用现代化数据栈高效处理地理信息数据

但是当数据需要展示在地图上时则需要将其原始的空间数据投射到大地坐标系上（这个过程称为投影）才可以得到这个几何图形具体的地理坐标。...空间引用识别号 (SRID) 要将几何图形投影到坐标系，必须需要使用 SRID。SRID 可以理解为唯一标识了将某个几何体空间数据映射成某个具体坐标系中的方式。...从而形成最终的地图数据供我们使用。基于定位的点和地图的基准点的之间的偏差就可以完成整个定位到地图的转换这个过程就是坐标个定位到地图的转换这个过程就是坐标系转换。...地理信息数据应用的问题我们会从一个具体案例来和大家探讨地理信息数据应用中会遇到的实际问题。我们这个地理数据应用案例如下：如何知道地球上一块土地在一段时间内的使用情况？...然后将这个矩形区域存储在一个新的字段中求 GIS 数据几何图形的中心点提前裁剪数据，将清洗好、裁剪好的规整数据写入对端新型数据库长周期的实时地理信息数据同步 CloudCanal 不仅支持历史数据的迁移同时还支持异构数据源之间的实时数据同步

1.1K2 0

Python pyproj 实现地理坐标转换

本文记录使用 Python 库 pyproj 实现地理坐标转换的流程。简介 pyproj是一个Python库，用于执行坐标转换和投影变换。...它基于Proj库，后者是一个C库，用于处理地图投影和坐标变换。pyproj提供了Python语言的接口，使得用户可以方便地使用这些功能。...由于各种原因（如地图制作、技术限制等），这些坐标可能需要从一个系统转换到另一个系统。pyproj库正是用来执行这些转换的。...易用性：pyproj的API设计简单直观，使得用户可以轻松地进行坐标变换和投影操作。社区支持：pyproj有一个活跃的社区，不断更新和改进库的功能。...://pyproj4.github.io/pyproj/stable/ 安装方法 pip install pyproj 地理坐标转换 pyproj 开发了很多 API ，这里简单介绍常用的地理坐标转换使用方法

2911 0

GIS坐标系测绘原理:大地水准面基准面参考椭球体EPSGSRIWKT

因为该射线是垂直于该平面(perpendicular)，因此称为正交(orthogonal)投影。现实生活中，从一大堆统计点中拟合出一条有规律的线，就需要用最小二乘法，其实就是正交投影的思路。...假设一个椭圆柱面与地球椭球体面横切于某一条经线上，按照等角条件将中央经线东、西各3°或1.5°经线范围内的经纬线投影到椭圆柱面上，然后将椭圆柱面展开成平面而成的。...兰勃特投影有两种：等角圆锥投影设想用一个正圆锥切于或割于球面，应用等角条件将地球面投影到圆锥面上，然后沿一母线展开成平面。投影后纬线为同心圆圆弧，经线为同心圆半径。...Web墨卡托投影有两个相关的投影标准，经常搞混： EPSG4326：Web墨卡托投影后的平面地图，但仍然使用WGS84的经度、纬度表示坐标； EPSG3857：Web墨卡托投影后的平面地图，坐标单位为米...EPSG:3857 (Pseudo-Mercator) 伪墨卡托投影，也被称为球体墨卡托，Web Mercator。它是基于墨卡托投影的，把 WGS84坐标系投影到正方形。

4.1K1 1

（数据科学学习手札75）基于geopandas的空间数据分析——坐标参考系篇

2 坐标参考系基础 2.1 CRS 　　在一个二维的平面中，我们可以使用如图1所示的坐标系统，通过坐标\((x_{0},y_{0})\)唯一确定点的位置：图1 　　现实世界中的地球作为一个球体...国内常用的基准面有：BEIJING1954，XIAN1980，WGS84等) units=m：声明坐标系单位设置为米 ellps=WGS84：声明椭球面（如何计算地球的圆度）使用WGS84 　　...　　与投影坐标系相比，没有单位units的信息，因为地理坐标系通常单位为十进制度数；而上述两个示例中都带有towgs84=0,0,0，这是一个转换因子，在需要进行数据转换时使用。...:4326，在当前的CRS下将其绘制出来：图18 　　利用to_crs()将其再投影到EPSG:2381并进行绘制：图19 　　通过比较可以发现，再投影之后的中国形变失真情况得到缓解...='red', alpha=0.8) # 先再投影到EPSG:2381 cq.to_crs(crs='EPSG:2381').plot(ax=ax, color='orange', markersize

1.6K3 0

Python地信专题 | 基于geopandas的空间数据分析-坐标参考系篇

，创建了包裹全球的网格，譬如WGS84，将本初子午线与赤道的交点作为原点（图4）：图4 WGS84地理坐标系及其经纬网格 2.1.2 投影坐标系地理坐标系虽然解决了我们在地球球面上定位的问题，但纬度和经度位置没有使用统一的测量单位...国内常用的基准面有：BEIJING1954，XIAN1980，WGS84等) units=m：声明坐标系单位设置为米 ellps=WGS84：声明椭球面（如何计算地球的圆度）使用WGS84 上述例子记录了投影坐标系的...towgs84=0,0,0，这是一个转换因子，在需要进行数据转换时使用。...:4326，在当前的CRS下将其绘制出来：图19 利用to_crs()将其再投影到EPSG:2381并进行绘制：图20 通过比较可以发现，再投影之后的中国形变失真情况得到缓解，且坐标系单位范围也发生了变化...red', alpha=0.8) # 先再投影到EPSG:2381 cq.to_crs(crs='EPSG:2381').plot(ax=ax, color='orange', markersize=100

1.8K2 1

2019-2023年全球固定宽带和移动（蜂窝）网络性能（更新）

数据以 Shapefile 格式和 Apache Parquet 格式提供，几何图形以 EPSG:4326 中的已知文本 (WKT) 表示。...为了创建一个易于管理的数据集，我们将原始数据汇总到磁贴中。磁贴数据的大小定义为 "缩放级别"（或 "z"）的函数。在 z=0 时，数据块的大小就是整个世界的大小。...根据这一定义，磁贴尺寸实际上是根据 Web 墨卡托投影法（EPSG:3857）计算的地球宽度/高度的一部分。因此，瓦片大小会因纬度不同而略有差异，但可以米为单位估算瓦片大小。...每个平铺图的几何图形在平铺图字段中以 WGS 84（EPSG:4326）表示。...时间周期和更新频率图层根据一个季度（三个月）的数据生成，文件将按季度更新和添加。

591 0

GEE数据——全球固定宽带和移动（蜂窝）网络性能（网速）（2019-2024）

数据以 Shapefile 格式和 Apache Parquet 格式提供，几何图形以 EPSG:4326 中的已知文本 (WKT) 表示。...为了创建一个易于管理的数据集，我们将原始数据汇总到磁贴中。数据块的大小定义为 "缩放级别"（或 "z"）的函数。在 z=0 时，磁贴的大小是整个世界的大小。...根据这一定义，磁贴尺寸实际上是根据 Web 墨卡托投影法（EPSG:3857）计算的地球宽度/高度的一部分。因此，瓦片尺寸会因纬度不同而略有差异，但瓦片尺寸可以米为单位进行估算。...每个平块的几何图形在平块字段中以 WGS 84（EPSG:4326）表示。...这对于在空间上连接来自多个时期（季度）的数据、在不使用地理空间函数的情况下创建更粗略的空间聚合、空间索引、分区以及存储和导出平铺几何图形都很有用。

771 0

Basemap系列：管理投影

所有的地图都拥有一个投影。投影方式及其特征在Basemap对象创建时确定。这种方式和其他库（比如 GDAL）的做法具有很大不同。...使用 epsg 设置投影 EPSG代码是一种使用数值代码命名投影的标准。Basemap允许使用这一准则创建地图，但仅适用于特定情况。...使用这种方式设置投影时可传递 epsg 参数给 Basemap 构造器。...这就是为什么有些投影会失败的原因，因为有些以经纬度设置的方形框在投影单位中不会给出一个适合的边界框。在此例中，使用的是 UTM (Transverse Mercator)投影。...（译注：仅当你知道设置的参数将满足你的要求时使用）上例展示了如何使用 plot 函数使用几个点显示坐标由 0 到 width 和 height的范围。

1.5K2 0

geotrellis使用（十四）导出定制的GeoTiff

二、需求说明很多时候我们需要从一块（或者很多块）大的Tiff中根据需要截取一部分数据，并且需要采用某种采样方式转成特定的投影，并转成需要的数据类型。...2.数据导入要想处理大数据或者处理多块栅格数据就不能直接处理栅格数据，可以先将栅格数据导入到Accumulo中，当然导入之后是一块块的瓦片，这部分在geotrellis使用（三）geotrellis...4.3 获取dstCRS、resampleMethod、cellType 首先可以肯定的是前端传来的这三个参数都是字符串型的，这就要求后台将字符串转成相应的类型。 ...投影方式我这里投了个懒，传递的是EPSG CODE,EPSG是投影方式的一种数字编码，具体请见：http://www.epsg.org/。...总之每一个EPSG编码对应了一种投影方式，像常见的经纬度投影的编码是4326，WebMercator的编码是3857。

9446 0

Web墨卡托公开的小秘密

墨卡托投影是一种“正轴等角切圆柱投影”,由荷兰地图学家墨卡托(G.Mercator)于1569年创立:假设地球被围在一个中空的圆柱里，其赤道与圆柱相接触，然后再假想地球中心有一盏灯，把球面上的图形投影到圆柱体上...伪墨卡托投影 Web Mercator 坐标系使用的投影方法不是严格意义的墨卡托投影，而是一个被 EPSG（European Petroleum Survey Group）称为伪墨卡托的投影方法，这个伪墨卡托投影方法的大名是...在投影过程中，将表示地球的参考椭球体近似的作为正球体处理（正球体半径 R = 椭球体半长轴 a）。...Auxiliary Sphere 就是在告知你，这个坐标在投影过程中，将椭球体近似为正球体做投影变换，虽然基准面是WGS 1984 椭球面。...ESRI:102113 Esri内部使用ID，与 EPSG:3785 相应。已被弃用。 ESRI:102100 Esri内部使用ID，与 EPSG:3857 相应。

2.6K1 0

地图基础

，能够标示地球上的任何一个位置对于地球上的某个位置，使用不同的测量基准得到的坐标是不一样的使用最广泛的WGS 84北美NAD 83、欧洲ETRS 89和我国CGCS 2000同WGS 84差异都非常小，...但是并不一致投影坐标系将三维的地理坐标转化为二维，叫做投影投影方式有很多，比如圆柱形、圆锥形、方形等，但是都会导致变形和失真常用的投影有灯具矩形投影(platte carre) 和墨卡托投影(mercator...，所以并不适合航海地图，但是因为坐标与像素之间的对应关系十分简单，比较适合gis软件用于栅格图的展示EPSG在国际上，每个坐标系统都会被分配一个EPSG代码EPSG: 4326EPSG: 4326是WPGS...84的代码WGS 84是目前最流行的地理坐标系统GPS是基于WGS 84的一般在存储数据时，按照WGS 84存储EPSG: 3857伪墨卡投影，也被称为球体墨卡托，web mercator基于墨卡托投影...，把WGS 84坐标系投影到正方形，切掉了南北极附近地区由于墨卡托投影等正形性的特点，在不同层级的图层上物体的形状保持不变，一个正方形可以不断地被划分为更多更小的正方形以显示更清晰的细节伪墨卡托坐标系非常适合显示数据

1.3K4 0

GIS世界里坐标系“黑话”

向北为正，向南为负；大地高（H）：指某点沿法线方向到参考椭球面的距离 [d3a1a13a152de5ce7bec7971295f12a3.png] （四）什么是PCS [872ad4e2690ddd80177caa464cf20d8a.png...Beijing 54---EPSG:4214 Xian 80---EPSG:4610 CGCS 2000---EPSG:4490 WGS 84---EPSG:4326 （二）投影坐标系 WGS 84.../ Pseudo-Mercator---EPSG:3857 WGS_1984_Web_Mercator_Auxiliary_Sphere---ESRI:102100（Esri内部使用ID，与 EPSG:...内部使用ID，与 EPSG:3785 相应）三、特殊的坐标系 1、GCJ02（国测局02 或火星坐标）这是一种偏移算法，而不是单独的一个坐标系。...我们习惯将加了偏移算法的坐标称为gcj02坐标。

2.2K5 0

WKT_WK是什么品牌

下面是一个有效创建和插入一个OGC空间对象的语句：　INSERT INTO SPATIALTabLE ( THE_GEOM,　THE_NAME ) 　　　VALUES ( GeomFromText...PostGIS的EWKB, EWKT 　　　OGC的格式只支持２D的地理要素，而且其相关联的SRID不能嵌入到输入输出表达中。...下面是一些EWKT的例子：　　　POINT（0 0 0）– XYZ 　　　SRID=32632;POINT(0 0) – XY with SRID 　　　由于SRID可以嵌入到EWKT...ESRI亦在其shape文件格式(*.prj)中使用WKT。...degree”, 0.0174532925199433, AUTHORITY[“EPSG”, “9122”]], AUTHORITY[“EPSG”, “4326”]] 投影坐标系的格式：<projected

6083 0

大地测量学：EPSG:4326、EPSG:3857

显而易见的是，从三维到二维的转化，必然会导致变形和失真，失真是不可避免的，但是不同投影下会有不同的失真，这让我们可以有得选择。...需要注意的是，对于墨卡托投影来说，越到高纬度，大小扭曲越严重，到两极会被放到无限大，所以，墨卡托投影无法显示极地地区。下图来自维基百科，可以看到墨卡托投影下每个国家的大小和实际大小的差异。...3.1 EPSG:4326 (WGS84) 前面说了 WGS84 是目前最流行的地理坐标系统。在国际上，每个坐标系统都会被分配一个 EPSG 代码，EPSG:4326 就是 WGS84 的代码。...3.2 EPSG:3857 (Pseudo-Mercator) 伪墨卡托投影，也被称为球体墨卡托，Web Mercator。它是基于墨卡托投影的，把 WGS84坐标系投影到正方形。...我们前面已经知道 WGS84 是基于椭球体的，但是伪墨卡托投影把坐标投影到球体上，这导致两极的失真变大，但是却更容易计算。这也许是为什么被称为”伪“墨卡托吧。

2.7K2 1

Google Earth Engine（GEE）——全球固定宽带和移动（蜂窝）网络性能数据集

数据以Shapefile格式和Apache Parquet格式提供，几何图形以WKT格式表示，以EPSG:4326投影。...为了创建一个可管理的数据集，我们将原始数据汇总成瓦片。数据瓦片的大小被定义为 "缩放级别"（或 "Z"）的函数。在Z=0的时候，一个瓦片的大小是整个世界的大小。...在Z=1的时候，瓦片在垂直和水平方向上被分成两半，形成4个覆盖全球的瓦片。这种瓦片分割随着缩放级别的增加而继续进行，当我们放大到一个特定的区域时，瓦片会以指数形式变小。...根据这个定义，瓦片的大小实际上是根据Web墨卡托投影（EPSG:3857）的地球宽度/高度的一些分数。因此，瓦片大小因纬度不同而略有不同，但瓦片大小可以用米来估计。...时间周期和更新频率图层是根据一个季度的数据（三个月）生成的，文件将按季度更新和添加。

1411 0

打开栅格数据的正确方式

以一个简单例子说明如何打开栅格影像下面的例子打开一副GeoTIFF影像，输出了影像的一些信息，然后遍历了所有波段，输出波段的一些信息 import gdal # 打开栅格数据集 ds = gdal.Open...('example.tif') # 获得栅格数据的一些重要信息 print(f'投影信息：{ds.GetProjection()}') print(f'栅格波段数：{ds.RasterCount}')...NoData值：-28672.0 统计值（最大值最小值）：(-468.0, 6265.0) 数据类型：Int16 NoData值：-28672.0 统计值（最大值最小值）：(21.0, 7267.0) 如何将...下面介绍几种转换的方法：在Dataset级别进行转换，转换结果是一个三维数组，第一个维度是波段数在Band级别进行转换，转换的结果是一个二维数据使用gdal_array模块中的LoadFile()...image = ds.ReadAsArray() print(f'数据的尺寸：{image.shape}') # 输出结果为：数据的尺寸：(3, 4800, 4800) # 这说明ReadAsArray方法将每个波段都转换为了一个二维数组

1.1K2 0

【C++】开源：地图投影和坐标转换proj库配置使用

项目介绍官网：https://proj.org/en/9.4/ 项目Github地址：https://github.com/OSGeo/PROJ proj 库是一个开源的库，主要用于地理坐标系之间的转换和地图投影...2.地图投影：将地球的三维表面映射到二维平面上（如地图）的方法。不同的投影方法适用于不同的应用场景。 3.坐标转换：将一种坐标系中的点转换到另一种坐标系中，通常涉及复杂的数学计算。...主要功能： 1.坐标转换：支持多种地理坐标系之间的转换，如从WGS84到UTM，或者从地理坐标系到地心坐标系。 2.地图投影：支持多种地图投影方式，如墨卡托投影、兰伯特正形投影等。...:32633 - UTM zone 33N) 和 WGS84坐标系 (EPSG:4326) PJ *utm_crs = proj_create(ctx, "EPSG:32633"); // UTM...transformation object." << std::endl; proj_context_destroy(ctx); return 1; } // 将转换对象规范化为视觉化使用

2241 0

使用GDAL进行影像投影坐标、地理坐标、图上坐标的转换

我使用GDAL库写了四个函数分别进行投影坐标与地理坐标（经纬度）之间的转换，投影坐标和图上坐标（行列号）之间的转换。有需要的朋友可以参考。...geosrs = prosrs.CloneGeogCS() return prosrs, geosrs def geo2lonlat(dataset, x, y): ''' 将投影坐标转为经纬度坐标...coords = ct.TransformPoint(x, y) return coords[:2] def lonlat2geo(dataset, lon, lat): ''' 将经纬度坐标转为投影坐标...) 注：关于投影坐标和图上坐标转换的六参数模型可以参考我的另外一篇博文：经纬度坐标和投影坐标的转换，其实质就是一个仿射变换。...具体参数可以使用gdallocationinfo –help查看。

8.4K2 0

使用GDAL读取Sentinel数据

数据的读取，我这里使用Sentinel-2光学卫星数据给出使用GDAL工具对其进行读取的方法。...GDAL将Sentinel数据看做一个数据集（概念上类似HDF格式的数据集），里面包含了很多子数据文件。所以，对于Sentinel数据的读取就和对于HDF数据的读取是相同的啦。...下图显示的数据子集中包含四个波段的数据（红，绿，蓝，近红外） image.png 使用GDAL命令行工具将Sentinel数据转为GeoTIFF格式转换是针对具体的子数据集而言的，所以使用gdal_translate...() # 将数据集中的数据转为ndarray del visual_arr # 获得栅格数据的一些重要信息 print(f'打开数据为：{ds_list[0][1]}') print(f'投影信息：...'栅格行数（高度）：{visual_ds.RasterYSize}') 程序输出如下：打开数据为：Bands B2, B3, B4, B8 with 10m resolution, UTM 13N 投影信息

1.7K0 0

点击加载更多

扫码

添加站长进交流群

领取专属 10元无门槛券

手把手带您无忧上云