[geotk] Reading netcdf GMT GRD as coverage

Brian Schlining bschlining at gmail.com
Thu Feb 17 18:16:00 EST 2011


Hi Martin,

I've attached the output below. Looks pretty good to me; let me know if
something is amiss. Also, I dumped a GRD file that you can play with at
http://www.mbari.org/staff/brian/pub/Monterey_UTM50.grd.zip . You can unzip
it and try it out. You'll have to change the extension to .nc so that
geotoolkit will read it.

Also, although the header on the ASC file didn't get written correctly, when
I edited it and added the missing info it rendered just fine. So the only
issue I'm having is with the header.

Its look like that the "gridToCRS" information is missing for some reason...
> Could you try the following please?
>
> System.out.println(gridGeometry);
>

*GridGeometry:*
GridGeometry2D[GeneralGridEnvelope[0..5082, 0..4066], PARAM_MT["Affine",
  PARAMETER["num_row", 3],
  PARAMETER["num_col", 3],
  PARAMETER["elt_0_0", 50.0],
  PARAMETER["elt_0_2", 422985.2775983825],
  PARAMETER["elt_1_1", -50.0],
  PARAMETER["elt_1_2", 4118864.9067777675]]]



>
> System.out.println(gridGeometry.gridToCRS);
>

*GridGeometry.gridToCRS:*
PARAM_MT["Affine",
  PARAMETER["num_row", 3],
  PARAMETER["num_col", 3],
  PARAMETER["elt_0_0", 50.0],
  PARAMETER["elt_0_2", 422985.2775983825],
  PARAMETER["elt_1_1", -50.0],
  PARAMETER["elt_1_2", 4118864.9067777675]]



>
> System.out.println(gridCoverage.sampleDimensions);
>

*GridCoverage.sampleDimensions*:
[RenderedSampleDimension["z"]]

Cheers
-- B
~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
Brian Schlining
bschlining at gmail.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/geotoolkit/attachments/20110217/a30d6bb8/attachment.html


More information about the Geotoolkit mailing list