Question

我创建了一个RDD，其中包含特定区域的DGM栅格并具有crs : EPSG 32632。我想将此数据覆盖在具有crs:EPSG 4326也是RDD的图像上以进行进一步分析。

尽管两个图像在QGIS中彼此重叠，但是投影系统和范围有所不同，这就是为什么稍后我尝试将两个RDD重新投影到webmercator时会给我一个黑色图像的原因。

val (_,inputrdd) : (Int, MultibandTileLayerRDD[SpatialKey]) = biggis.landuse.spark.s3.layerFromS3(bucket, name, Some(tileSize))

   val store = new S3AttributeStore(bucket,name){
    override def s3Client = s3ClientURL()
  }
  val writer = S3COGLayerWriter(store)
  val (_,rasterrdd) : (Int, MultibandTileLayerRDD[SpatialKey]) = biggis.landuse.spark.s3.layerFromS3(bucket1, name1, Some(tileSize))

  val raster_store = new S3AttributeStore(bucket1,name1){
    override def s3Client = s3ClientURL()
  }
  val raster_writer = S3COGLayerWriter(raster_store)

inputrdd与4326一起使用。

如何将rasterrdd重新投影到inputrdd投影系统？

Answer 1

由于rasterrdd的类型为MultibandTileLayerRDD[K]，因此您可以使用所有地理网格函数进行重新投影；它会像：

rasterrdd.reproject(LatLng, ZoomedLayoutScheme)

欢迎您在docs中查找更详细的信息，并加入我们的gitter channel，以解决任何其他问题和讨论。

如何将包含栅格的RDD重新投影到不同的crs（WGS84）？

1 个答案: