Как исправить недостающую SRS в геотифе с помощью геотреллиса

Question

У меня есть куча геотифов без SRS. Исправить SRS можно с помощью команды GDAL:

gdal_translate -a_srs EPSG:25832 srcfile targetfile

Как заменить команду 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 если вы работаете с несколькими группами.