如何将含有栅格的RDD重新投影到不同的CRS(WGS84)

Question

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

尽管这两个图像在QGIS中彼此都位于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]，您可以使用所有geotrellis函数执行重新投入；它会像：

一样

rasterrdd.reproject(LatLng, ZoomedLayoutScheme)

欢迎您在文档中寻找更详细的信息，并加入我们的吉特频道进行任何其他问题和讨论。