Hi,
I'm trying to operate several .asc coverages and write de result to new
GridCoverage2D but when I finish my operations and write to disk, if I
do gdalinfo over .asc result file the statistics data there aren't in
the file.
For example:
gdalinfo over one source file (prados_clip.asc):
|Driver: AAIGrid/Arc/Info ASCII Grid
Files: prado_clip.asc
prado_clip.asc.aux.xml
Size is 540, 492
Coordinate System is `'
Origin = (617151.000000000000000,4814861.000435999594629)
Pixel Size = (20.000000000000000,-20.000000000000000)
Corner Coordinates:
Upper Left ( 617151.000, 4814861.000)
Lower Left ( 617151.000, 4805021.000)
Upper Right ( 627951.000, 4814861.000)
Lower Right ( 627951.000, 4805021.000)
Center ( 622551.000, 4809941.000)
Band 1 Block=540x1 Type=Float32, ColorInterp=Undefined
Min=0.022 Max=0.816
Minimum=0.022, Maximum=0.816, Mean=0.391, StdDev=0.195
NoData Value=-99999
Metadata:
STATISTICS_MAXIMUM=0.81599998474121
STATISTICS_MEAN=0.39072748197934
STATISTICS_MINIMUM=0.021999999880791
STATISTICS_STDDEV=0.19466570119278
|
gdalinfo over result coverage created:
|Files: salidaWS.asc
salidaWS.prj
Size is 540, 492
Coordinate System is:
LOCAL_CS["Wildcard 2D cartesian plane in metric unit",
LOCAL_DATUM["Unknown",0],
UNIT["m",1.0],
AXIS["x",EAST],
AXIS["y",NORTH],
AUTHORITY["EPSG","404000"]]
Origin = (617151.000000000000000,4814861.000435999594629)
Pixel Size = (20.000000000000000,-20.000000000000000)
Corner Coordinates:
Upper Left ( 617151.000, 4814861.000)
Lower Left ( 617151.000, 4805021.000)
Upper Right ( 627951.000, 4814861.000)
Lower Right ( 627951.000, 4805021.000)
Center ( 622551.000, 4809941.000)
Band 1 Block=540x1 Type=Float32, ColorInterp=Undefined
NoData Value=-99999
|
Code to create the coverage:
| /**
*
* @param name Future name of coverage
* @param coverage One of source coverage to copy metadata
* @param data future data content of coverage
* @return new GridCoverage2D
*/
private static GridCoverage2D createCoverage(String name,
GridCoverage2D coverage, double[] data){
Raster raster = coverage.getRenderedImage().getData();
WritableRaster wr2 = raster.createCompatibleWritableRaster();
wr2.setPixels(raster.getMinX(), raster.getMinY(),
raster.getWidth(), raster.getHeight(), data);
TiledImage ti = new TiledImage(coverage.getRenderedImage(), false);
ti.setData(wr2);
RenderingHints rh =
ImageUtilities.getRenderingHints(coverage.getRenderedImage());
ImageLayout imgLayout =
ImageUtilities.getImageLayout(coverage.getRenderedImage());
imgLayout.setColorModel(ti.getColorModel());
ti.setProperty("GC_NODATA", getNoDataValue(coverage));
RenderedImage ri = FormatDescriptor.create(ti,
DataBuffer.TYPE_DOUBLE, rh);
GridSampleDimension[] gsds = new GridSampleDimension[1];
Color[] colorramp = null; //create gray scale color
Category[] categories = {
coverage.getSampleDimensions()[0].getCategories().get(0),
// To get NoDataValue Category
new Category(name, colorramp,
NumberRange.create(0.0,false,1.0,false), 1,0) // Category range (0.0 , 1.0)
};
gsds[0] = new GridSampleDimension(name, categories, Unit.ONE);
GridCoverageFactory fac =
CoverageFactoryFinder.getGridCoverageFactory(null);
GridCoverage2D coverageTemp = fac.create(name, ri,
coverage.getGridGeometry(), gsds,null, null);
return coverageTemp;
}
|
And my code to write result operations to GridCoverage2D:
| public static int writeCoverateToDisk(GridCoverage2D coverage,
String name){
try {
File imageFile = new File(name);
ArcGridWriter gridWriter = new ArcGridWriter(imageFile);
gridWriter.write((GridCoverage)coverage, null);
return 0;
} catch (IllegalArgumentException ex) {
logger.log(Level.SEVERE, null, ex);
return -1;
} catch (IOException ex) {
logger.log(Level.SEVERE, null, ex);
return -1;
}
}
|
Somebody knows how insert statistics data in result file?
Thanks in Advance
--
Pablo Díaz Redondo
LaboraTe - USC
http://laborate.usc.es
[email protected] / [email protected]
Telf.: 982823301
------------------------------------------------------------------------------
Everyone hates slow websites. So do we.
Make your web apps faster with AppDynamics
Download AppDynamics Lite for free today:
http://p.sf.net/sfu/appdyn_d2d_mar
_______________________________________________
GeoTools-GT2-Users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/geotools-gt2-users
