org.jgrasstools.gears.io.rasterreader.OmsRasterReader Maven / Gradle / Ivy
The newest version!
/*
* This file is part of JGrasstools (http://www.jgrasstools.org)
* (C) HydroloGIS - www.hydrologis.com
*
* JGrasstools is free software: you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation, either version 3 of the License, or
* (at your option) any later version.
*
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with this program. If not, see envelopeParams = getRegionParamsFromGridCoverage(outRaster);
double west = envelopeParams.get(WEST);
double south = envelopeParams.get(SOUTH);
double east = envelopeParams.get(EAST);
double north = envelopeParams.get(NORTH);
double xres = envelopeParams.get(XRES);
double yres = envelopeParams.get(YRES);
if (pBounds == null) {
pBounds = new double[]{north, south, west, east};
}
if (pRes == null) {
pRes = new double[]{xres, yres};
}
double n = pBounds[0];
double s = pBounds[1];
double w = pBounds[2];
double e = pBounds[3];
if (pRes != null || pRowcol != null) {
int newRows = 0;
int newCols = 0;
if (pRowcol != null) {
newRows = (int) pRowcol[0];
newCols = (int) pRowcol[1];
if (pRes == null) {
pRes = new double[2];
}
pRes[0] = (e - w) / (double) newCols;
pRes[1] = (n - s) / (double) newRows;
} else if (pRes != null) {
pRowcol = new int[2];
newRows = (int) Math.round((n - s) / pRes[1]);
newCols = (int) Math.round((e - w) / pRes[0]);
}
pRowcol[0] = newRows;
pRowcol[1] = newCols;
}
HashMap newParams = makeRegionParamsMap(n, s, w, e, pRes[0], pRes[1], (int) pRowcol[1], (int) pRowcol[0]);
CoordinateReferenceSystem crs = outRaster.getCoordinateReferenceSystem();
GridGeometry2D gg = gridGeometryFromRegionParams(newParams, crs);
outRaster = (GridCoverage2D) Operations.DEFAULT.resample(outRaster, crs, gg, null);
}
private void checkNovalues() {
// TODO make this nice, this can't be the way
if (fileNovalue == null || geodataNovalue == null) {
return;
}
if (isNovalue(internalFileNovalue) && isNovalue(internalGeodataNovalue)) {
return;
}
if (!NumericsUtilities.dEq(internalFileNovalue, internalGeodataNovalue)) {
HashMap params = getRegionParamsFromGridCoverage(outRaster);
int height = params.get(ROWS).intValue();
int width = params.get(COLS).intValue();
WritableRaster tmpWR = createDoubleWritableRaster(width, height, null, null, null);
WritableRandomIter tmpIter = RandomIterFactory.createWritable(tmpWR, null);
RenderedImage readRI = outRaster.getRenderedImage();
RandomIter readIter = RandomIterFactory.create(readRI, null);
int minX = readRI.getMinX();
int minY = readRI.getMinY();
for( int r = 0; r < height; r++ ) {
for( int c = 0; c < width; c++ ) {
double value = readIter.getSampleDouble(c + minX, r + minY, 0);
if (isNovalue(value) || value == internalFileNovalue || value == -Float.MAX_VALUE || value == Float.MAX_VALUE) {
tmpIter.setSample(c, r, 0, internalGeodataNovalue);
} else {
tmpIter.setSample(c, r, 0, value);
}
}
}
readIter.done();
tmpIter.done();
outRaster = buildCoverage(new File(file).getName(), tmpWR, params, outRaster.getCoordinateReferenceSystem());
}
}
private boolean hasBoundsRequest() {
return pNorth != null && pSouth != null && pWest != null && pEast != null;
}
private boolean hasRowColsRequest() {
return pRows != null && pCols != null;
}
private boolean hasResolutionRequest() {
return pXres != null && pYres != null;
}
/**
* Utility method to quickly read a grid in default mode.
*
* @param path the path to the file.
* @return the read coverage.
* @throws Exception
*/
public static GridCoverage2D readRaster( String path ) throws Exception {
OmsRasterReader reader = new OmsRasterReader();
reader.file = path;
reader.process();
GridCoverage2D geodata = reader.outRaster;
return geodata;
}
}