Hi All,<br><br>trying to create a SQL with gdal2wktraster.py. It works on PNG but it doesn't on AIG files... <br><br>Note: As Gdal works fine.<br><br>C:\Program Files\PostgreSQL\8.4\bin>gdalinfo c:\tmp\raster<br>Driver: AIG/Arc/Info Binary Grid<br>
Files: c:\tmp\raster<br> c:\tmp\raster\dblbnd.adf<br> c:\tmp\raster\hdr.adf<br> c:\tmp\raster\metadata.xml<br> c:\tmp\raster\prj.adf<br> c:\tmp\raster\sta.adf<br> c:\tmp\raster\vat.adf<br>
c:\tmp\raster\w001000.adf<br> c:\tmp\raster\w001000x.adf<br> c:\tmp\raster\w001001.adf<br> c:\tmp\raster\w001001x.adf<br> c:\tmp\raster\z001001.adf<br> c:\tmp\raster\z001001x.adf<br> c:\tmp\raster\z001002.adf<br>
c:\tmp\raster\z001002x.adf<br> c:\tmp\raster\z001003.adf<br> c:\tmp\raster\z001003x.adf<br> c:\tmp\raster\z001004.adf<br> c:\tmp\raster\z001004x.adf<br> c:\tmp\raster\z001005.adf<br> c:\tmp\raster\z001005x.adf<br>
c:\tmp\raster\z001006.adf<br> c:\tmp\raster\z001006x.adf<br> c:\tmp\raster\z001007.adf<br> c:\tmp\raster\z001007x.adf<br> c:\tmp\raster\z001008.adf<br> c:\tmp\raster\z001008x.adf<br> c:\tmp\raster\z001009.adf<br>
c:\tmp\raster\z001009x.adf<br> c:\tmp\raster\z001010.adf<br> c:\tmp\raster\z001010x.adf<br> c:\tmp\raster\z001011.adf<br> c:\tmp\raster\z001011x.adf<br> c:\tmp\raster\z001012.adf<br> c:\tmp\raster\z001012x.adf<br>
c:\tmp\raster\z001013.adf<br> c:\tmp\raster\z001013x.adf<br> c:\tmp\raster\z001014.adf<br> c:\tmp\raster\z001014x.adf<br> c:\tmp\raster\z001015.adf<br> c:\tmp\raster\z001015x.adf<br>Size is 107759, 168633<br>
Coordinate System is:<br>PROJCS["unnamed",<br> GEOGCS["Unknown datum based upon the Airy 1830 ellipsoid",<br> DATUM["Not_specified_based_on_Airy_1830_ellipsoid",<br> SPHEROID["Airy 1830",6377563.396,299.3249646,<br>
AUTHORITY["EPSG","7001"]],<br> AUTHORITY["EPSG","6001"]],<br> PRIMEM["Greenwich",0,<br> AUTHORITY["EPSG","8901"]],<br>
UNIT["degree",0.01745329251994328,<br> AUTHORITY["EPSG","9122"]],<br> AUTHORITY["EPSG","4001"]],<br> PROJECTION["Transverse_Mercator"],<br>
PARAMETER["latitude_of_origin",49],<br> PARAMETER["central_meridian",-2],<br> PARAMETER["scale_factor",0.9996012717],<br> PARAMETER["false_easting",400000],<br> PARAMETER["false_northing",-100000],<br>
UNIT["METERS",1]]<br>Origin = (128110.000000000000000,813270.000000000000000)<br>Pixel Size = (5.000000000000000,-5.000000000000000)<br>Corner Coordinates:<br>Upper Left ( 128110.000, 813270.000) ( 6d29'37.32"W, 57d 7'47.53"N)<br>
Lower Left ( 128110.000, -29895.000) ( 5d45'40.00"W, 49d34'10.24"N)<br>Upper Right ( 666905.000, 813270.000) ( 2d24'41.72"E, 57d 7'58.04"N)<br>Lower Right ( 666905.000, -29895.000) ( 1d41'32.29"E, 49d34'18.23"N)<br>
Center ( 397507.500, 391687.500) ( 2d 2'15.04"W, 53d25'18.19"N)<br>Band 1 Block=256x4 Type=Byte, ColorInterp=Undefined<br> Min=1.000 Max=4.000<br> NoData Value=255<br clear="all"><br><br>C:\Program Files\PostgreSQL\8.4\bin>python gdal2wktraster.py -r c:\tmp\raster -t sebastable -o c:\tmp\sebas.sql<br>
gdal2wktraster.py:644: DeprecationWarning: 'H' format requires 0 <= number <= 65535<br> hexstr = binascii.hexlify(struct.pack(fmt_little, data)).upper()<br>gdal2wktraster.py:644: DeprecationWarning: integer argument expected, got float<br>
hexstr = binascii.hexlify(struct.pack(fmt_little, data)).upper()<br>ERROR 2: Multiplication overflow : 107759 * 168633 * 1<br>Traceback (most recent call last):<br> File "gdal2wktraster.py", line 1013, in <module><br>
main()<br> File "gdal2wktraster.py", line 976, in main<br> wkblify_raster(opts, filename, i)<br> File "gdal2wktraster.py", line 921, in wkblify_raster<br> summary = wkblify_raster_level(options, ds, options.overview_level, band_range, infile, i)<br>
File "gdal2wktraster.py", line 888, in wkblify_raster_level<br> hexwkb += wkblify_band(options, band, level, xoff, yoff, read_block_size, block_size, infile, b)<br> File "gdal2wktraster.py", line 777, in wkblify_band<br>
target_block_size[0], target_block_size[1])<br> File "C:\OSGeo4W\apps\gdal-16\pymod\osgeo\gdal.py", line 835, in ReadAsArray<br> buf_xsize, buf_ysize, buf_obj )<br> File "C:\OSGeo4W\apps\gdal-16\pymod\osgeo\gdal_array.py", line 140, in BandReadAsArray<br>
ar = numpy.reshape(ar, [buf_ysize,buf_xsize])<br> File "C:\OSGeo4W\apps\Python25\lib\site-packages\numpy\core\fromnumeric.py", line 116, in reshape<br> return reshape(newshape, order=order)<br>ValueError: total size of new array must be unchanged<br>
<br>Any ideas ?<br>-- <br>Sebastian E. Ovide<br><br><br><br>