如何使用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。