I am not so good at the raster rendering chain, I usually put a break point in here and then step through the renderImage <https://github.com/geotools/geotools/blob/master/modules/library/render/src/main/java/org/geotools/renderer/lite/gridcoverage2d/GridCoverageRenderer.java#L770> method checking that each operation makes sense and is not being fed bad data. So in this case I would check with a world+image file that works, and compare against this geotiff file that does not work.
-- Jody Garnett On 27 October 2015 at 11:52, Neal Dillman <ndill...@allweatherinc.com> wrote: > Jody, > > > > The results from ImageLab are: > > > > PROJCS["WGS 84 / World Mercator", > > GEOGCS["WGS 84", > > DATUM["World Geodetic System 1984", > > SPHEROID["WGS 84", 6378137.0, 298.257223563, > AUTHORITY["EPSG","7030"]], > > AUTHORITY["EPSG","6326"]], > > PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], > > UNIT["degree", 0.017453292519943295], > > AXIS["Geodetic longitude", EAST], > > AXIS["Geodetic latitude", NORTH], > > AUTHORITY["EPSG","4326"]], > > PROJECTION["Mercator_1SP", AUTHORITY["EPSG","9804"]], > > PARAMETER["latitude_of_origin", 0.0], > > PARAMETER["central_meridian", 0.0], > > PARAMETER["scale_factor", 1.0], > > PARAMETER["false_easting", 0.0], > > PARAMETER["false_northing", 0.0], > > UNIT["m", 1.0], > > AXIS["Easting", EAST], > > AXIS["Northing", NORTH], > > AUTHORITY["EPSG","3395"]] > > > > Those results are the same whether or not I call GeoTiffReader directly, > and it appears that in either case the right thing is happening, just no > painting. > > > > Regards, > > Neal > > > > *From:* Jody Garnett [mailto:jody.garn...@gmail.com] > *Sent:* Tuesday, October 27, 2015 11:13 AM > *To:* Neal Dillman <ndill...@allweatherinc.com> > *Cc:* Devon Tucker <devonrtuc...@gmail.com>; > geotools-gt2-users@lists.sourceforge.net; ijtur...@gmail.com > *Subject:* Re: [Geotools-gt2-users] ImageLab and TIFF > > > > So can we stop in a debugger and see what projection information is being > extracted from the GeoTiff? Would like to confirm it was able to correctly > identify EPSG:3395 before digging deep into math transforms. > > > > Note we may also have some luck with matching a CRS > <http://docs.geotools.org/latest/userguide/library/referencing/crs.html#matching-a-coordinatereferencesystem> > - and using the official definition of EPSG:3395 (rather than what comes > out of the prj or geotiff metadata). > > > > -- > > Jody Garnett > > > > On 26 October 2015 at 18:28, Neal Dillman <ndill...@allweatherinc.com> > wrote: > > Hi Devon, > > > > Thanks for looking further into this. I agree that adding the hint works > for the bluemarble tif. Unfortunately that trick is not working for my > examples are EPSG:3395 (Mercator/WGS84). I still get the same silent (as > in blank) fail from those files. > > > > Regards, > > Neal > > > > *From:* Devon Tucker [mailto:devonrtuc...@gmail.com] > *Sent:* Monday, October 26, 2015 5:26 PM > *To:* Neal Dillman <ndill...@allweatherinc.com> > *Cc:* geotools-gt2-users@lists.sourceforge.net; ijtur...@gmail.com; > jody.garn...@gmail.com > *Subject:* Re: [Geotools-gt2-users] ImageLab and TIFF > > > > Hi all, > > I did some playing around with this. I'm able to get BlueMarble to render > by either: > > Giving the LONGITUDE_FIRST hint when creating the GeoTiff reader: > > new GeoTiffReader(rasterFile, new > Hints(Hints.FORCE_LONGITUDE_FIRST_AXIS_ORDER, Boolean.TRUE)); > > or using an EPSG:4326 .prj sidecar file along with it. > > Not quite sure what the issue is otherwise. > > >
------------------------------------------------------------------------------
_______________________________________________ GeoTools-GT2-Users mailing list GeoTools-GT2-Users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/geotools-gt2-users
