gis - 如何使用 geotrellis 修复 geotiff 中丢失的 SRS

Question

我有一堆没有 SRS 的 geotiffs。可以使用 GDAL 命令修复 SRS：

gdal_translate -a_srs EPSG:25832 srcfile targetfile

如何使用 geotrellis 代码替换 GDAL 命令？

Answer 1

import geotrellis.proj4._
import geotrellis.raster.io.geotiff._

val fileName: String = ???
val tif = SinglebandGeoTiff(fileName)
val targetCrs = CRS.fromName("EPSG:25832”)
val reprojected = tif.projectedRaster.reproject(targetCrs)
GeoTiff(reprojected.raster,  reprojected.crs).write("/Users/eugene/reprojected.tif")

MultibandGeoTiff如果您正在使用多个乐队，则可以使用。