[gdal-dev] Question on writing a GTiff using python

[Stephen Woodbridge]

Hi all,

I'm writing a GTiff using GDAL/Python and it is mostly working. The 
source data is coming from a NETCDF file that is georeferenced in 
Longitude from 74.1600037 to 434.160 degrees

If I use this I get:

Size is 4500, 1782
Coordinate System is:
GEOGCS["WGS 84",
     DATUM["WGS_1984",
         SPHEROID["WGS 84",6378137,298.257223563,
             AUTHORITY["EPSG","7030"]],
         AUTHORITY["EPSG","6326"]],
     PRIMEM["Greenwich",0],
     UNIT["degree",0.0174532925199433],
     AUTHORITY["EPSG","4326"]]
Origin = (74.160003662109375,46.987300872802734)
Pixel Size = (0.079999959309896,0.079999959309896)

Corner Coordinates:
Upper Left  (  74.1600037,  46.9873009) ( 74d 9'36.01"E, 46d59'14.28"N)
Lower Left  (      74.160,     189.547) ( 74d 9'36.01"E,189d32'50.02"N)
Upper Right (     434.160,      46.987) (Invalid angle, 46d59'14.28"N)
Lower Right (     434.160,     189.547) (Invalid angle,189d32'50.02"N)
Center      (     254.160,     118.267) (254d 9'35.68"E,118d16' 2.15"N)

Looks like I also need to flip the Y axis, but ignore that for the moment.

Do I need to reorganize the pixels to be in the range of -180 to 180?
What is the best way t do this?
I can probably slice and splice the numpy array for each row as I copy 
it. Other options?

Thanks,
   -Steve

---
This email has been checked for viruses by Avast antivirus software.
https://www.avast.com/antivirus