I am reading SRMT+ data using the GTopo30Reader. The 9.5 version works
great, but 10-RC2 gives incorrect values for coordinates in the middle of
the grid. This looks like a problem with how tiles are set up in the
underlying GridCoverage, but I'm not sure what initialization parameters
are set wrong, and whether this is a problem with how the GTopo30Reader
sets up the GridCoverage2D or whether the error spans all GridCoverage2D
instances.
Specific steps to reproduce:
I am loading the SMTP+ elevation data from *
ftp://topex.ucsd.edu/pub/srtm30_plus/srtm30/erm/**
w140n40.Bathymetry.srtm.dem* (with associated ers,hdr,prj, and stx files)
with the code
GTopo30Reader reader = new GTopo30Reader( demFile );
GridCoverage2D coverage = reader.read(null);
This is a large 4800x6000 cell dataset. The returned coverage is broken up
into 480x500 pixel tiles internally (see
GridCoverage2D#evaluate(finalGridCoordinates2D coord,
final int[] dest)
In 10-RC2 these tiles are not drawing data from their 480x500 pixel subset
of the file, but rather from a 4800x50 pixel subset. Thus, only the first
and last 480 pixels return the correct elevation.
For this specific file, the following test passes in 9.5 and fails in 10-RC*
pos = new GridCoordinates2D(480,0);
assertEquals(-4337,coverage.evaluate(pos,(int[])null)[0]);
I'll be happy to provide the auxiliary files I'm working with if they're
needed to reproduce the bug.
-Spencer
------------------------------------------------------------------------------
Learn the latest--Visual Studio 2012, SharePoint 2013, SQL 2012, more!
Discover the easy way to master current and previous Microsoft technologies
and advance your career. Get an incredible 1,500+ hours of step-by-step
tutorial videos with LearnDevNow. Subscribe today and save!
http://pubads.g.doubleclick.net/gampad/clk?id=58041391&iu=/4140/ostg.clktrk
_______________________________________________
GeoTools-GT2-Users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/geotools-gt2-users
