0
Я использую библиотеку Scala и Geotrellis для обработки и анализа изображений Landsat-8, но когда я попытался объединить полосы 2, 3 и 4 (красный , зеленый и синий) и вот в результате многополосный TIFF (справа) и PNG генерируется из него (слева направо):Объединение каналов RGB на ландшафт и преобразование в png дает ужасный результат
код для объединения однозонных размолвки является:
class Tiler (outputPath : String, bandIds : Array[String]){
def bandPath(b: String) = s"data/landsat/LC81750342016185LGN00_B${b}.TIF"
def obtainBands() : (ArrayMultibandTile, SinglebandGeoTiff) = {
var tileArrayBuffer = ArrayBuffer[Tile]()
var geoTiff : SinglebandGeoTiff = null
bandIds.foreach((bandId) => {
geoTiff = SinglebandGeoTiff(bandPath(bandId))
tileArrayBuffer += geoTiff.tile
})
(ArrayMultibandTile(tileArrayBuffer), geoTiff)
}
def writeResultToFS(multiBand : MultibandTile, extent : Extent, crs : CRS) : Unit = {
MultibandGeoTiff(multiBand, extent, crs).write(outputPath)
}
def parseTiles() : Unit = {
val (multiBandTileArray, geotiff) = obtainBands()
writeResultToFS(multiBandTileArray, geotiff.extent, geotiff.crs)
}
}
Кроме , это мой код конвертера png:
class PngTest() {
val tiffPath = "data/output/output.tif"
val bwPath = "data/output/bw.png"
def makepng(): Unit = {
println("Rendering PNG and saving to disk...")
val conversionMap = {
val multibandTile = MultibandGeoTiff(tiffPath).tile.convert(IntConstantNoDataCellType)
multibandTile.combine(0, 1, 2) { (rBand, gBand, bBand) =>
val r = if (isData(rBand)) { rBand } else 0
val g = if (isData(gBand)) { gBand } else 0
val b = if (isData(bBand)) { bBand } else 0
if(r + g + b == 0) 0
else {
val color = (r + g + b)/3
((color & 0xFF) << 24) | ((color & 0xFF) << 16) | ((color & 0xFF) << 8) | 0xFF
}
}
}
conversionMap.renderPng().write(bwPath)
}
}
Когда я удалить карту преобразования и просто сделать MultibandGeoTiff (tiffPath) .tile.convert (IntConstantNoDataCellType) .renderPng(). Написать (bwPath), все это я получаю это красочная версия на картинке слева.
Извините, если это вопрос о нобе и спасибо за любую помощь заранее.