Changeset 5577 for spike/wktraster/scripts/gdal2wktraster.py

Timestamp:

04/27/10 16:06:05 (2 years ago)

Author:

jorgearevalo

Message:

Not necessary to change the geotransform array returned by GDAL driver GetGeoTransform? function for non-referenced rasters. Each driver returns a correct geotransform array. Related ticket: #498.

Files:

• spike/wktraster/scripts/gdal2wktraster.py (modified) (3 diffs)

spike/wktraster/scripts/gdal2wktraster.py

@@ -191 +191 @@
     """If verbose mode requested, sends extra progress information to stderr"""
     if VERBOSE is True:
-        sys.stderr.write(msg)
+        sys.stderr.write(msg)        
 
 
@@ -569 +569 @@
     return (xpad, ypad)
 
-def is_georeferenced(gt):
-    if gt[0] == 0.0 and gt[1] == 1.0 and gt[3] == 0.0 and gt[5] == -1.0:
-        return False
-    else:
-        return True
-
 def get_gdal_geotransform(ds):
     assert ds is not None
     gt = list(ds.GetGeoTransform())
-    # GDAL returns a positive Y pixelsize making ungeoreferenced images 
-    # to be badly georeferenced. Here we fix this by negating the Y  
-    # pixelsize when the image is not georeferenced.
-    if gt[0] == 0.0 and gt[1] == 1.0 and gt[3] == 0.0 and gt[5] == 1.0:
-        gt[5] = -1.0  
     return tuple(gt)
 
@@ -591 +580 @@
     assert len(xy) == 2
 
-    if is_georeferenced(gt):
-        xgeo = gt[0] + gt[1] * xy[0] + gt[2] * xy[1];
-        ygeo = gt[3] + gt[4] * xy[0] + gt[5] * xy[1];
-    else:
-        xgeo = xy[0]
-        ygeo = xy[1]
-    
+    xgeo = gt[0] + gt[1] * xy[0] + gt[2] * xy[1];
+    ygeo = gt[3] + gt[4] * xy[0] + gt[5] * xy[1];
+
     return (xgeo, ygeo)