<p>Thanks for the GMT suggestions. I've never used it, I'll try it out.</p>
<p>Thanks!!</p>
<p>Typed with thumbs...</p>
<div class="gmail_quote">On Apr 12, 2011 7:24 PM, "Joaquim Luis" <<a href="mailto:jluis@ualg.pt">jluis@ualg.pt</a>> wrote:<br type="attribution">> Jerl,<br>> <br>> I would do that with GMT. By using ps2raster you can convert any GMT <br>
> figure (in geogs of course) into KML, and GMT has excellent care of date <br>> line transition issues.<br>> <br>> Joaquim<br>> <br>>> Hello:<br>>><br>>> I stopped by the #gdal channel on IRC to ask this question. Thank you<br>
>> for the help provided there. I thought it might be better if I asked<br>>> a question here so I can better explain what's going on in one spot<br>>> and not flood the channel.<br>>> I understand this is a pixel center vs pixel edge issue.<br>
>><br>>> My ultimate goal is to get from a netCDF file to a google map overlay contour.<br>>> The general steps I've taken, which work all the way up to about<br>>> 178.29 then it just cuts off.<br>
>><br>>> netCDF -> color-relief PNG using "gdaldem color-relief ... "<br>>> PNG -> Geocoded PNG using "gdal_translate -a_srs EPSG:4326 -outsize<br>>> 4096 4096 -a_ullr -180.25 90.25 180.25 -90.25" on the PNG (Have also<br>
>> tried this with just "-180 90 180 -90"<br>>> PNG -> tiles using "gdal2tiles.py -s EPSG:4326 -z 0 -n -w google" on<br>>> the Geocoded PNG<br>>> I guess first question, is my methodology sound? Maybe there's a<br>
>> better way/tool to approach the problem.<br>>><br>>> Following is the output from gdalinfo, other than that, what<br>>> information might I be able to provide that could help?<br>>><br>>> Thank you,<br>
>><br>>> Jerl<br>>><br>>> gdalinfo for the source netcdf: (My apologies, I'm not certain what is<br>>> relevant so I'm including most everything.)<br>>><br>>> Driver: netCDF/Network Common Data Format<br>
>> Files: <a href="http://test.nc">test.nc</a><br>>> Size is 721, 361<br>>> Coordinate System is `'<br>>> Origin = (-180.250000000000000,90.250000000000000)<br>>> Pixel Size = (0.500000000000000,-0.500000000000000)<br>
>> Metadata:<br>>> NC_GLOBAL#Conventions=CF-1.0<br>>> NC_GLOBAL#cdm_data_type=GRID<br>>> NC_GLOBAL#file_format=GEMPAK<br>>> NC_GLOBAL#location=test.grd<br>>> NC_GLOBAL#history=Direct read of GEMPAK into NetCDF-Java 4.0 API<br>
>> NC_GLOBAL#_CoordinateModelRunDate=2011-04-12T00:00:00Z<br>>> UVIN_NONE#long_name=testdata<br>>> UVIN_NONE#units=<br>>> UVIN_NONE#missing_value=-9.999000e+03<br>>> UVIN_NONE#VectorComponentFlag=gridRelative<br>
>> time#long_name=forecast time<br>>> time#units=minute since 2011-04-12T00:00:00Z<br>>> time#_CoordinateAxisType=Time<br>>> lat#units=degrees_north<br>>> lat#long_name=latitude coordinate<br>
>> lat#standard_name=latitude<br>>> lat#grid_spacing=0.5 degrees_north<br>>> lat#_CoordinateAxisType=Lat<br>>> lon#units=degrees_east<br>>> lon#long_name=longitude coordinate<br>
>> lon#standard_name=longitude<br>>> lon#grid_spacing=0.5 degrees_east<br>>> lon#_CoordinateAxisType=Lon<br>>> Corner Coordinates:<br>>> Upper Left (-180.2500000, 90.2500000)<br>>> Lower Left (-180.2500000, -90.2500000)<br>
>> Upper Right ( 180.2500000, 90.2500000)<br>>> Lower Right ( 180.2500000, -90.2500000)<br>>> Center ( 0.0000000, 0.0000000)<br>>><br>>> gdalinfo of geocoded PNG:<br>>><br>>> Driver: PNG/Portable Network Graphics<br>
>> Files: testWorld.png<br>>> testWorld.png.aux.xml<br>>> Size is 4096, 4096<br>>> Coordinate System is:<br>>> GEOGCS["WGS 84",<br>>> DATUM["WGS_1984",<br>
>> SPHEROID["WGS 84",6378137,298.257223563,<br>>> AUTHORITY["EPSG","7030"]],<br>>> AUTHORITY["EPSG","6326"]],<br>>> PRIMEM["Greenwich",0,<br>
>> AUTHORITY["EPSG","8901"]],<br>>> UNIT["degree",0.01745329251994328,<br>>> AUTHORITY["EPSG","9122"]],<br>>> AUTHORITY["EPSG","4326"]]<br>
>> Origin = (-180.000000000000000,90.000000000000000)<br>>> Pixel Size = (0.087890625000000,-0.043945312500000)<br>>> Image Structure Metadata:<br>>> INTERLEAVE=PIXEL<br>>> Corner Coordinates:<br>
>> Upper Left (-180.0000000, 90.0000000) (180d 0'0.00"W, 90d 0'0.00"N)<br>>> Lower Left (-180.0000000, -90.0000000) (180d 0'0.00"W, 90d 0'0.00"S)<br>>> Upper Right ( 180.0000000, 90.0000000) (180d 0'0.00"E, 90d 0'0.00"N)<br>
>> Lower Right ( 180.0000000, -90.0000000) (180d 0'0.00"E, 90d 0'0.00"S)<br>>> Center ( 0.0000000, 0.0000000) ( 0d 0'0.01"E, 0d 0'0.01"N)<br>>> Band 1 Block=4096x1 Type=Byte, ColorInterp=Red<br>
>> Mask Flags: PER_DATASET ALPHA<br>>> Band 2 Block=4096x1 Type=Byte, ColorInterp=Green<br>>> Mask Flags: PER_DATASET ALPHA<br>>> Band 3 Block=4096x1 Type=Byte, ColorInterp=Blue<br>>> Mask Flags: PER_DATASET ALPHA<br>
>> Band 4 Block=4096x1 Type=Byte, ColorInterp=Alpha<br>>><br>>><br>>><br>> <br></div>