Exemplo n.º 1
0
  /**
   * Checking if NoData and ROI are defined correctly
   *
   * @param indexed
   * @param image
   * @param roi
   * @param nodata
   * @param destNoData
   */
  private void checkNoDataROI(
      RenderedImage indexed, RenderedImage image, ROI roi, Range nodata, int destNoData) {
    // Ensure the dimensions are the same
    assertEquals(indexed.getMinX(), image.getMinX());
    assertEquals(indexed.getMinY(), image.getMinY());
    assertEquals(indexed.getWidth(), image.getWidth());
    assertEquals(indexed.getHeight(), image.getHeight());

    boolean roiExists = roi != null;
    boolean nodataExists = nodata != null;
    // Simply ensure no exception is thrown
    if (!nodataExists && !roiExists) {
      PlanarImage.wrapRenderedImage(indexed).getTiles();
      return;
    }

    if (nodataExists) {
      nodata = RangeFactory.convertToDoubleRange(nodata);
    }
    RandomIter roiIter = null;
    Rectangle roiBounds = null;
    if (roiExists) {
      PlanarImage roiIMG = roi.getAsImage();
      roiIter = RandomIterFactory.create(roiIMG, null, true, true);
      roiBounds = roi.getBounds();
    }

    // Else check ROI and NoData
    RandomIter sourceIter = RandomIterFactory.create(image, null, true, true);
    RandomIter destIter = RandomIterFactory.create(indexed, null, true, true);
    // Start the iteration (we iterate only the first band)
    int w = image.getWidth();
    int h = image.getHeight();
    int minX = image.getMinX();
    int minY = image.getMinY();
    int maxX = minX + w;
    int maxY = minY + h;
    int limx = minX - image.getTileGridXOffset();
    int limy = minY - image.getTileGridYOffset();
    Rectangle translated = new Rectangle(limx, limy, w, h);

    for (int y = minY; y < maxY; y++) {
      for (int x = minX; x < maxX; x++) {

        double src = sourceIter.getSampleDouble(x, y, 0);
        double dest = destIter.getSampleDouble(x, y, 0);

        boolean valid = true;

        // ROI Check
        if (roiExists
            && !(roiBounds.contains(x, y) && roiIter.getSample(x, y, 0) > 0)
            && translated.contains(x, y)) {
          valid = false;
        }

        // NoData Check
        if (nodataExists && nodata.contains(src)) {
          valid = false;
        }
        if (!valid) {
          assertEquals(destNoData, dest, TOLERANCE);
        }
      }
    }
  }
Exemplo n.º 2
0
  /**
   * Returns a tile for reading.
   *
   * @param tileX The X index of the tile.
   * @param tileY The Y index of the tile.
   * @return The tile as a <code>Raster</code>.
   */
  public Raster computeTile(int tileX, int tileY) {
    // STATISTICAL ELABORATIONS
    // selection of the format tags
    RasterFormatTag[] formatTags = getFormatTags();
    // Selection of the RasterAccessor parameters
    Raster source = getSourceImage(0).getTile(tileX, tileY);
    // Control if the Period is bigger than the tile dimension, in that case, the
    // statistics are not updated
    if (xPeriod > getTileWidth() || yPeriod > getTileHeight()) {
      return source;
    }

    Rectangle srcRect = getSourceImage(0).getBounds().intersection(source.getBounds());
    // creation of the RasterAccessor
    RasterAccessor src =
        new RasterAccessor(source, srcRect, formatTags[0], getSourceImage(0).getColorModel());

    // ROI calculations if roiAccessor is used
    RasterAccessor roi = null;
    RandomIter roiIter = null;
    if (useROIAccessor) {
      // Note that the getExtendedData() method is not called because the input images are padded.
      // For each image there is a check if the rectangle is contained inside the source image;
      // if this not happen, the data is taken from the padded image.
      Raster roiRaster = null;
      if (srcROIImage.getBounds().contains(srcRect)) {
        roiRaster = srcROIImage.getData(srcRect);
      } else {
        roiRaster = srcROIImgExt.getData(srcRect);
      }

      // creation of the rasterAccessor
      roi =
          new RasterAccessor(
              roiRaster,
              srcRect,
              RasterAccessor.findCompatibleTags(new RenderedImage[] {srcROIImage}, srcROIImage)[0],
              srcROIImage.getColorModel());
    } else if (hasROI) {
      roiIter = RandomIterFactory.create(srcROIImage, srcROIImage.getBounds(), true, true);
    }

    // Computation of the statistics
    switch (src.getDataType()) {
      case DataBuffer.TYPE_BYTE:
        byteLoop(src, srcRect, roi, stats, roiIter);
        break;
      case DataBuffer.TYPE_USHORT:
        ushortLoop(src, srcRect, roi, stats, roiIter);
        break;
      case DataBuffer.TYPE_SHORT:
        shortLoop(src, srcRect, roi, stats, roiIter);
        break;
      case DataBuffer.TYPE_INT:
        intLoop(src, srcRect, roi, stats, roiIter);
        break;
      case DataBuffer.TYPE_FLOAT:
        floatLoop(src, srcRect, roi, stats, roiIter);
        break;
      case DataBuffer.TYPE_DOUBLE:
        doubleLoop(src, srcRect, roi, stats, roiIter);
        break;
    }

    return source;
  }