Re: [gdal-dev] Projecting points on a mercator projection

2010-04-29 Thread Brice Lambi
This actually works just as expected but I was being dense. The odd points with the negative coords are in Alaska and Hawaii which is outside of this domain. After removing these and accounting for the image size: gx = gx % 1073 gy = gy % 689 everything works as expected. la...@ucar.edu wrot

Re: [gdal-dev] Projecting points on a mercator projection

2010-04-29 Thread Brice Lambi
Also, this is the gdalinfo of the image I am trying to match: Driver: GTiff/GeoTIFF Files: /tmp/1725.tif Size is 1073, 689 Coordinate System is: PROJCS["unnamed", GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7

[gdal-dev] Projecting points on a mercator projection

2010-04-29 Thread lambi
Hello, I'm trying to figure out how to take a list of points and put them on a mercator projected image using the gdal python bindings. I'm trying to match it with a NDFD grid that I've been able to create and put on google maps from a GRIB2 file. This is what I have: geo = [-14483681.357462916