У меня есть группа геотипов без SRS. можно зафиксировать с помощью SRS GDAL команды:Как исправить недостающие SRS в геотифе с использованием geotrellis
gdal_translate -a_srs EPSG:25832 srcfile targetfile
Как я могу заменить команду GDAL с помощью geotrellis коды?
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, если вы работаете с несколькими полосами.
но это фактически имитирует gdal_warp, правильно? –
@ViliamSimko загляните в https://github.com/geotrellis/geotrellis/pull/1697 (я дважды проверил, что тиифы без crs будут загружены как tiff с помощью LatLng), однако вы можете указать дополнительные crs во время процесса проглатывания, можете ли вы проверить эту новую функцию? Если нет, любой пример плитки будет отличным! – DaunnC
@ViliamSimko выглядит как да, это _reprojects_, не присваивает crs; – DaunnC