• Home
• com.azavea.geotrellis
• geotrellis_2.10
• 0.9.2
• source code
• Histogram.scala

×

Project download

Many resources are needed to download a project. Please understand that we have to compensate our server costs. Thank you in advance.
Project price only 1 $

You can buy this project and download/modify it how often you want.

geotrellis.raster.op.zonal.summary.Histogram.scala Maven / Gradle / Ivy

package geotrellis.raster.op.zonal.summary

import geotrellis._
import geotrellis.source._
import geotrellis.feature._
import geotrellis.feature.rasterize._
import geotrellis.statistics._

object Histogram extends TileSummary[Histogram,Histogram,ValueSource[Histogram]] {
  def handlePartialTile[D](pt:PartialTileIntersection[D]):Histogram = {
    val PartialTileIntersection(r,polygons) = pt
    val histogram = FastMapHistogram()
    for(p <- polygons.asInstanceOf[List[Polygon[D]]]) {
      Rasterizer.foreachCellByFeature(p, r.rasterExtent)(
        new Callback[Geometry,D] {
          def apply (col:Int, row:Int, g:Geometry[D]) {
            val z = r.get(col,row)
            if (isData(z)) histogram.countItem(z, 1)
          }
        }
      )
    }

    histogram
  }

  def handleFullTile(ft:FullTileIntersection):Histogram = {
    val histogram = FastMapHistogram()
    ft.tile.foreach((z:Int) => if (isData(z)) histogram.countItem(z, 1))
    histogram
  }

  def converge(ds:DataSource[Histogram,_]) =
    ds.map(x=>x).converge // Map to kick in the CanBuildFrom for HistogramDS
}