[gdal-dev] Re: Dateline gap

[Jerl Simpson]

Does anyone on the GDAL team have any suggestions on how to do this?


Thank you,

Jerl


On Tue, Apr 12, 2011 at 2:59 PM, Jerl Simpson <jsimpson at wxtrends.com> wrote:
> Hello:
>
> I stopped by the #gdal channel on IRC to ask this question.  Thank you
> for the help provided there.  I thought it might be better if I asked
> a question here so I can better explain what's going on in one spot
> and not flood the channel.
> I understand this is a pixel center vs pixel edge issue.
>
> My ultimate goal is to get from a netCDF file to a google map overlay contour.
> The general steps I've taken, which work all the way up to about
> 178.29 then it just cuts off.
>
> netCDF -> color-relief PNG using "gdaldem color-relief ...  "
> PNG -> Geocoded PNG using "gdal_translate -a_srs EPSG:4326 -outsize
> 4096 4096 -a_ullr -180.25 90.25 180.25 -90.25" on the PNG (Have also
> tried this with just "-180 90 180 -90"
> PNG -> tiles using "gdal2tiles.py -s EPSG:4326 -z 0 -n -w google" on
> the Geocoded PNG
> I guess first question, is my methodology sound?  Maybe there's a
> better way/tool to approach the problem.
>
> Following is the output from gdalinfo, other than that, what
> information might I be able to provide that could help?
>
> Thank you,
>
> Jerl
>
> gdalinfo for the source netcdf: (My apologies, I'm not certain what is
> relevant so I'm including most everything.)
>
> Driver: netCDF/Network Common Data Format
> Files: test.nc
> Size is 721, 361
> Coordinate System is `'
> Origin = (-180.250000000000000,90.250000000000000)
> Pixel Size = (0.500000000000000,-0.500000000000000)
> Metadata:
>  NC_GLOBAL#Conventions=CF-1.0
>  NC_GLOBAL#cdm_data_type=GRID
>  NC_GLOBAL#file_format=GEMPAK
>  NC_GLOBAL#location=test.grd
>  NC_GLOBAL#history=Direct read of GEMPAK into NetCDF-Java 4.0 API
>  NC_GLOBAL#_CoordinateModelRunDate=2011-04-12T00:00:00Z
>  UVIN_NONE#long_name=testdata
>  UVIN_NONE#units=
>  UVIN_NONE#missing_value=-9.999000e+03
>  UVIN_NONE#VectorComponentFlag=gridRelative
>  time#long_name=forecast time
>  time#units=minute since 2011-04-12T00:00:00Z
>  time#_CoordinateAxisType=Time
>  lat#units=degrees_north
>  lat#long_name=latitude coordinate
>  lat#standard_name=latitude
>  lat#grid_spacing=0.5 degrees_north
>  lat#_CoordinateAxisType=Lat
>  lon#units=degrees_east
>  lon#long_name=longitude coordinate
>  lon#standard_name=longitude
>  lon#grid_spacing=0.5 degrees_east
>  lon#_CoordinateAxisType=Lon
> Corner Coordinates:
> Upper Left  (-180.2500000,  90.2500000)
> Lower Left  (-180.2500000, -90.2500000)
> Upper Right ( 180.2500000,  90.2500000)
> Lower Right ( 180.2500000, -90.2500000)
> Center      (   0.0000000,   0.0000000)
>
> gdalinfo of geocoded PNG:
>
> Driver: PNG/Portable Network Graphics
> Files: testWorld.png
>       testWorld.png.aux.xml
> Size is 4096, 4096
> Coordinate System is:
> GEOGCS["WGS 84",
>    DATUM["WGS_1984",
>        SPHEROID["WGS 84",6378137,298.257223563,
>            AUTHORITY["EPSG","7030"]],
>        AUTHORITY["EPSG","6326"]],
>    PRIMEM["Greenwich",0,
>        AUTHORITY["EPSG","8901"]],
>    UNIT["degree",0.01745329251994328,
>        AUTHORITY["EPSG","9122"]],
>    AUTHORITY["EPSG","4326"]]
> Origin = (-180.000000000000000,90.000000000000000)
> Pixel Size = (0.087890625000000,-0.043945312500000)
> Image Structure Metadata:
>  INTERLEAVE=PIXEL
> Corner Coordinates:
> Upper Left  (-180.0000000,  90.0000000) (180d 0'0.00"W, 90d 0'0.00"N)
> Lower Left  (-180.0000000, -90.0000000) (180d 0'0.00"W, 90d 0'0.00"S)
> Upper Right ( 180.0000000,  90.0000000) (180d 0'0.00"E, 90d 0'0.00"N)
> Lower Right ( 180.0000000, -90.0000000) (180d 0'0.00"E, 90d 0'0.00"S)
> Center      (   0.0000000,   0.0000000) (  0d 0'0.01"E,  0d 0'0.01"N)
> Band 1 Block=4096x1 Type=Byte, ColorInterp=Red
>  Mask Flags: PER_DATASET ALPHA
> Band 2 Block=4096x1 Type=Byte, ColorInterp=Green
>  Mask Flags: PER_DATASET ALPHA
> Band 3 Block=4096x1 Type=Byte, ColorInterp=Blue
>  Mask Flags: PER_DATASET ALPHA
> Band 4 Block=4096x1 Type=Byte, ColorInterp=Alpha
>
>
>
> --
> This communication is privileged and may contain confidential
> information.  It's intended only for the use of the person or entity
> named above. If you are not the intended recipient, do not distribute
> or copy this communication. If you have received this communication in
> error, please notify the sender immediately and return the original to
> the email address above. © Copyright 2011 Weather Trends
> International, Inc.
>



-- 
This communication is privileged and may contain confidential
information.  It's intended only for the use of the person or entity
named above. If you are not the intended recipient, do not distribute
or copy this communication. If you have received this communication in
error, please notify the sender immediately and return the original to
the email address above. © Copyright 2011 Weather Trends
International, Inc.