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

Question

我有一群没有SRS的地质学家。可以使用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。