改进在 Landsat 上组合 RGB 通道并转换为 PNG 时的输出

Question

我正在使用 Scala 和 Geotrellis 库对 Landsat-8 图像进行处理和分析，但是当我尝试组合波段 2、3 和 4（红色、绿色和蓝色）时，下面是由此生成的多波段 tiff（右）和从中生成的 png（左）：

https://i.stack.imgur.com/uKI58.jpg

组合单频 tiffs 的代码是：

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（）.write（bwPath）时，我得到的只是左图的彩色版本。

抱歉，如果这是一个菜鸟问题，请提前感谢您的任何帮助。

Answer 1

如果要

将 landsat 像素值视为字节，则需要将其归一化为 0 - 255。带有掩码0xFF渲染代码将值视为字节。

下面是一个 PR，显示了使代码正常工作的更改：https://github.com/ardilgulez/gt-1/pull/1

下面是一些其他代码的示例，这些代码对 landsat 图像进行规范化和轻微的颜色校正，以便将其渲染为 web 地图的 PNG：https://github.com/geotrellis/geotrellis-landsat-emr-demo/blob/57353688471805b133a04a97ff86b26021e8d7b5/server/src/main/scala/demo/Render.scala#L22