Hi Pierre,<br><br>Does gdal2wktraster.py have any limitation on the maximum number of columnsxrows ?<br><br>in my case, my raster is 107759 x 168633...<br><br>gdal works well:<br><br>sebas@SeansPC:~/rasters$ gdal_translate -of PNG raster/ test.png<br>
Input file size is 107759, 168633<br>0...10...20...30...40...50...60...70...80...90...100 - done.<br><br><br>sebas@SeansPC:~/rasters$ python gdal2wktraster.py -r raster/ -t sebastable -o ok.sql<br>gdal2wktraster.py:695: DeprecationWarning: 'H' format requires 0 <= number <= 65535<br>
hexwkb += wkblify('H', xsize)<br>gdal2wktraster.py:696: DeprecationWarning: 'H' format requires 0 <= number <= 65535<br> hexwkb += wkblify('H', ysize)<br>gdal2wktraster.py:727: DeprecationWarning: integer argument expected, got float<br>
hexwkb += wkblify(pt2fmt(pixtype), nodata)<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 "/usr/lib/python2.6/dist-packages/osgeo/gdal.py", line 895, in ReadAsArray<br> buf_xsize, buf_ysize, buf_obj )<br> File "/usr/lib/python2.6/dist-packages/osgeo/gdal_array.py", line 228, in BandReadAsArray<br>
ar = numpy.empty([buf_ysize,buf_xsize], dtype = typecode)<br>MemoryError<br><br><br>sebas@SeansPC:~/rasters$ python gdal2wktraster.py -r test.png -t sebastable -o ok.sql<br>gdal2wktraster.py:695: DeprecationWarning: 'H' format requires 0 <= number <= 65535<br>
hexwkb += wkblify('H', xsize)<br>gdal2wktraster.py:696: DeprecationWarning: 'H' format requires 0 <= number <= 65535<br> hexwkb += wkblify('H', ysize)<br>gdal2wktraster.py:727: DeprecationWarning: integer argument expected, got float<br>
hexwkb += wkblify(pt2fmt(pixtype), nodata)<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 "/usr/lib/python2.6/dist-packages/osgeo/gdal.py", line 895, in ReadAsArray<br> buf_xsize, buf_ysize, buf_obj )<br> File "/usr/lib/python2.6/dist-packages/osgeo/gdal_array.py", line 228, in BandReadAsArray<br>
ar = numpy.empty([buf_ysize,buf_xsize], dtype = typecode)<br>MemoryError<br><br><br><br><br><br><div class="gmail_quote">On Mon, Jul 19, 2010 at 5:56 PM, Pierre Racine <span dir="ltr"><<a href="mailto:Pierre.Racine@sbf.ulaval.ca">Pierre.Racine@sbf.ulaval.ca</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div link="blue" vlink="purple" lang="EN-US">
<div>
<p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);">Hi Sebastian,</span></p>
<p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);"> </span></p>
<p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);">I can convert ESRI Grid file to .sql without problem using
gdal2wktraster.py and the same parameters as you. I can do both integer and
floating point rasters.</span></p>
<p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);"> </span></p>
<p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);">Maybe this is a GDAL problem. Try to convert it using
gdal_translate (to tiff for example). This would be a better test than just
gdalinfo.</span></p>
<p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);"> </span></p>
<p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);">Could you provide us with a file sample?</span></p>
<p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);"> </span></p>
<p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);">Pierre</span></p>
<p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);"> </span></p>
<div style="border-width: medium medium medium 1.5pt; border-style: none none none solid; border-color: -moz-use-text-color -moz-use-text-color -moz-use-text-color blue; padding: 0cm 0cm 0cm 4pt;">
<div>
<div style="border-width: 1pt medium medium; border-style: solid none none; border-color: rgb(181, 196, 223) -moz-use-text-color -moz-use-text-color; padding: 3pt 0cm 0cm;">
<p class="MsoNormal"><b><span style="font-size: 10pt;">From:</span></b><span style="font-size: 10pt;"> <a href="mailto:postgis-users-bounces@postgis.refractions.net" target="_blank">postgis-users-bounces@postgis.refractions.net</a>
[mailto:<a href="mailto:postgis-users-bounces@postgis.refractions.net" target="_blank">postgis-users-bounces@postgis.refractions.net</a>] <b>On Behalf Of </b>Sebastian
E. Ovide<br>
<b>Sent:</b> 19 juillet 2010 12:28<br>
<b>To:</b> PostGIS Users Discussion<br>
<b>Subject:</b> [postgis-users] WKTRaster : gdal2wktraster.py cannot read
AIG/Arc/Info Binary Grid</span></p>
</div>
</div><div><div></div><div class="h5">
<p class="MsoNormal"> </p>
<p class="MsoNormal" style="margin-bottom: 12pt;">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>
</p>
</div></div></div>
</div>
</div>
<br>_______________________________________________<br>
postgis-users mailing list<br>
<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
<a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br>Sebastian E. Ovide<br><br>skype: sebastian.ovide<br><br>+353 (0) 87 6340149<br><br><br>