Hi,
I'm trying to load rasters with the code below and I get negative y axis values.
I've got the same problem with the geotiff downloaded at
http://www.csc.noaa.gov/crs/tcm/sample_geotiff.html
So it seems to come from my code and the use of geotools.
I've used both
GeoTiffReader
=============
GeoTiffReader reader = (GeoTiffReader) ((new GeoTiffFormat()).getReader(new File(NOM_FICHIER_GEOTIFF)));
GridCoverage gc = reader.read(null);
GridCoverageReader
==================
StreamGridCoverageExchange gce = new StreamGridCoverageExchange();
GridCoverageReader reader = gce.getReader(new File(NOM_FICHIER_GEOTIFF));
Format format = reader.getFormat();
// get the parameters and set them
/* ParameterDescriptorGroup paramDescriptor = format.getReadParameters();
ParameterGroup params = (ParameterGroup) paramDescriptor.createValue();
params.parameter( "Compressed" ).setValue( true );
params.parameter( "GRASS" ).setValue( true );
*/
// read the grid
GridCoverage gc = reader.read( null );
Any ideas where it does come from ?
Jérémy
***************
code :
------
MapContext mapContext = new DefaultMapContext();
StyleBuilder styleBuilder = new StyleBuilder();
Envelope envelope = new Envelope();
// On charge le GridCoverage
GeoTiffReader reader = (GeoTiffReader) ((new
GeoTiffFormat()).getReader(new File(NOM_FICHIER_GEOTIFF)));
GridCoverage gc = reader.read(null);
debug("gc.getCoordinateReferenceSystem() >" +
gc.getCoordinateReferenceSystem());
debug("gc.getEnvelope() >" + gc.getEnvelope());
output :
--------
Could not find mediaLib accelerator wrapper classes. Continuing in pure
Java mode.
[DEBUG_GeotiffTest] gc.getCoordinateReferenceSystem() >PROJCS["NAD83 / UTM
zone 17N",
GEOGCS["NAD83",
DATUM["North_American_Datum_1983",
SPHEROID["GRS 1980", 6378137.0, 298.257222101,
AUTHORITY["EPSG","7019"]],
TOWGS84[0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0],
AUTHORITY["EPSG","6269"]],
PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],
UNIT["degree", 0.017453292519943295],
AXIS["Lon", EAST],
AXIS["Lat", NORTH],
AUTHORITY["EPSG","4269"]],
PROJECTION["Transverse_Mercator", AUTHORITY["OGC","Transverse_Mercator"]],
PARAMETER["central_meridian", -81.0],
PARAMETER["latitude_of_origin", 0.0],
PARAMETER["scale_factor", 0.9996],
PARAMETER["false_easting", 500000.0],
PARAMETER["false_northing", 0.0],
UNIT["m", 1.0],
AXIS["x", EAST],
AXIS["y", NORTH],
AUTHORITY["EPSG","26917"]]
[DEBUG_GeotiffTest] gc.getEnvelope() >GeneralEnvelope[703983.695935,
-3738421.777099, 704893.695935, -3737739.777099]
-------------------------------------------------------
This SF.net email is sponsored by: Splunk Inc. Do you grep through log files
for problems? Stop! Download the new AJAX search engine that makes
searching your log files as easy as surfing the web. DOWNLOAD SPLUNK!
http://ads.osdn.com/?ad_id=7637&alloc_id=16865&op=click
_______________________________________________
Geotools-gt2-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/geotools-gt2-users
