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
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
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