[postgis-commits] svn - r3856 - spike/wktraster/scripts

[postgis-commits at postgis.refractions.net]

Author: mloskot
Date: 2009-03-11 12:46:07 -0700 (Wed, 11 Mar 2009)
New Revision: 3856

Modified:
   spike/wktraster/scripts/gdal2wktraster.py
Log:
gdal2wktraster.py: Testing if loader perhaps breaks raster data, so there is a bands shift in the output. To be removed after tests.

Modified: spike/wktraster/scripts/gdal2wktraster.py
===================================================================
--- spike/wktraster/scripts/gdal2wktraster.py	2009-03-11 17:35:21 UTC (rev 3855)
+++ spike/wktraster/scripts/gdal2wktraster.py	2009-03-11 19:46:07 UTC (rev 3856)
@@ -33,6 +33,10 @@
 import binascii
 import sys, glob
 
+# TODO: *** DUMPER *** Testing if loading breaks raster data.
+# This is to dump source raster to text file for testing raw pixel values
+#fdumper = open('dump.txt', 'w')
+
 ###############################################################################
 # CONSTANTS - DO NOT CHANGE
 
@@ -198,6 +202,7 @@
 
 def wkblify_band(band):
     assert band is not None, "Error: No raster band"
+    import struct # TODO: Remove
 
     hexwkb = ""
     pixtype = gdt2pt(band.DataType)['id']
@@ -212,6 +217,12 @@
     # Rows of pixels of raster 
     for row in range(0, band.YSize):
         scanline = band.ReadRaster(0, row, band.XSize, 1, band.XSize, 1)
+
+        # TODO: *** DUMPER *** Dump scanline to file
+        #tuple = struct.unpack('B' * band.XSize, scanline)
+        #s = '\t'.join([str(v) for v in tuple]) + '\n'
+        #fdumper.write(s)
+
         hexwkb += binascii.hexlify(scanline)
         logit('.')
     logit('\n')
@@ -245,6 +256,9 @@
             logit("Cannot fetch raster band of index %d" % bn)
             continue
         
+        # TODO: *** DUMPER - Mark band being dumped
+        #fdumper.write('--- BAND %d ------------------------------\n' % bn)
+
         logit('MSG: Band #%d\n' % bn)
         hexwkb += wkblify_band(band)
         band = None