thanks for that Peter<br><br><div class="gmail_quote">On Thu, Jul 22, 2010 at 9:35 AM, Peter Baumann <span dir="ltr"><<a href="mailto:p.baumann@jacobs-university.de">p.baumann@jacobs-university.de</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 text="#000000" bgcolor="#ffffff">
Sebastian,<br>
<br>
as it will take some time until we complete our performance evaluation,
here a quick response taken from the log file on my notebook. I have
run a few rasql queries (while all other demon stuff plus my desktop
apps were running as well):<br>
<br>
cold access to 3 images of 256x211 8bit greyscale:<br>
--- snip -------------<br>
[2010-07-22 10:11:35] request from 127.0.0.1<br>
Request: executeQquery 'select mr from mr'...parsing...checking
semantics...optimizing (level 3)...No optimization library found<br>
evaluating...ok, result type 'set <marray <char,
[0:255,0:210]>>', 3 element(s).<br>
[2010-07-22 10:11:35] request completed; <b>request time=0.049s</b><br>
--- snip -------------<br>
<br>
same, hot access:<br>
--- snip -------------<br>
[2010-07-22 10:16:27] request completed; <b>request time=0.007s</b><br>
--- snip -------------<br>
<br>
accessing a 2D slice 512x512 (pixel type: unsigned short) from a <b>3D
weather data cube</b> 16384x144x73:<br>
--- snip -------------<br>
[2010-07-22 10:25:20] request from 127.0.0.1<br>
Request: executeQquery 'select m[0:511,0:511,0] from
ten_metre_U_wind_component as m'...parsing...checking
semantics...optimizing (level 3)...No optimization library found<br>
evaluating...ok, result type 'set <marray <ushort,
[0:511,0:143]>>', 1 element(s).<br>
[2010-07-22 10:25:21] request completed; <b>request time=0.290s</b><br>
--- snip -------------<br>
<br>
accessing a 3D subcube 100x100x73 from same cube:<br>
--- snip -------------<br>
[2010-07-22 10:26:49] request from 127.0.0.1<br>
Request: executeQquery 'select m[0:99,0:99,*:*] from
ten_metre_U_wind_component as m'...parsing...checking
semantics...optimizing (level 3)...No optimi<br>
zation library found<br>
evaluating...ok, result type 'set <marray <ushort,
[0:99,0:99,0:72]>>', 1 element(s).<br>
[2010-07-22 10:26:49] request completed; <b>request time=0.293s</b><br>
--- snip -------------<br>
<br>
extracting a 1024x73 x/t slice from same cube:<br>
--- snip -------------<br>
[2010-07-22 10:28:27] request from 127.0.0.1<br>
Request: executeQquery 'select m[0,0:1023,*:*] from
ten_metre_U_wind_component as m'...parsing...checking
semantics...optimizing (level 3)...No optimization library found<br>
evaluating...ok, result type 'set <marray <ushort,
[0:143,0:72]>>', 1 element(s).<br>
[2010-07-22 10:28:28] request completed; <b>request time=0.289s</b><br>
--- snip -------------<br>
<br>
As you see, slicing in all 3 dimensions delivers roughly the same times
(although tiling has not been optimized).<br>
<br>
laptop specs:<br>
CPU: Intel(R) Core(TM)2 Duo CPU     T8300  @ 2.40GHz<br>
Mem:   2GB (swap off)<br>
Linux floridita 2.6.32-23-generic #37-Ubuntu SMP Fri Jun 11 07:54:58
UTC 2010 i686 GNU/Linux<br>
rasdaman running without optimization package (cf log).<br>
<br>
Of course not a scientific test (that will follow), but it may give an
idea.<br>
<br>
cheers,<br>
Peter<br>
<br>
<br>
<br>
On 07/21/2010 11:00 AM, Peter Baumann wrote:
<blockquote type="cite">
  
Sebastian,<br>
  <br>
it's a XEON. Your question comes at the right time, we are in the
process of compiling performance data this summer. I will soon post
detail info here.<br>
  <br>
Loading 100k x 150k pixels is no problem. You load piecemeal, and the
system will align with its internal tiling + maintain the image
pyramid. <br>
  <br>
10 years back we have loaded the whole ortho image mosaic of Bavaria
into the database (server was a desktop PC of about 1000 EUR), that was
600 GB compressed (losslessly). At that time, import of an 230 MB image
took about 2min (array placement, re-tiling of input, index
maintenance, compression, and this for all pyramid levels). Read
access, OTOH, was about 150ms. A few years later the Thuringia Mapping
Agency, again with a desktop PC as a server, verified that access times
are stable under a load of about 15 concurrent users challenging the
system with WMS type interactive browsing. <br>
  <br>
For another mapping agency we have set up a database containing ortho
image, 30+ thematic raster layers, DEM. Access was via WMS, the
rasdaman system delivered a single image for display (Javascript was
not what it is today...). This included on-the-fly classification of
the DEM and styling the thematic layers.<br>
  <br>
So we have already done complete raster servers for several German
States.<br>
  <br>
BTW, standards: We have established the OGC Web Coverage Processing
Service (WCPS) standard which lifts the rasdaman concept of a raster
querfy language to an OGC standard. Being editor of currently 9 geo
raster service specs in OGC we also pursue reference implementations on
the basis of rasdaman. See <a href="http://www.earthlook.org" target="_blank">www.earthlook.org</a>
for demos.<br>
  <br>
-Peter<br>
  <br>
PS: Should data really exceed disk capacity (hard to believe today)
then we have a PhD thesis work which has implemented transparent tape
cabinetaccess. <br>
  <br>
  <br>
  <br>
On 07/21/2010 09:37 AM, Sebastian E. Ovide wrote:
  <blockquote type="cite"><br>
    <br>
    <div class="gmail_quote">On Tue, Jul 20, 2010 at 5:08 PM, Peter
Baumann <span dir="ltr"><<a href="mailto:p.baumann@jacobs-university.de" target="_blank">p.baumann@jacobs-university.de</a>></span>
wrote:<br>
    <blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
      <div text="#000000" bgcolor="#ffffff">
      <div>On 07/20/2010 05:47 PM, Pierre Racine wrote:
      <blockquote type="cite">
        <div>
        <p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);">Since
we have this bug with big rasters, </span></p>
        </div>
      </blockquote>
      <br>
      </div>
proven, any-size, ... ;-)<br>
      </div>
    </blockquote>
    <div><br>
what do you mean ? were you able to load a raster of 100kx150k ?<br>
 </div>
    <blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
      <div text="#000000" bgcolor="#ffffff"><br>
Who is proven? Well, running a dozen-TB seamless mosaic on PostgreSQL
for many years, having online-demos available since years, etc... you
decide, in comparison.<br>
      <font color="#888888"><br>
      </font></div>
    </blockquote>
    <div><br>
Peter, how big is your cluster to support dozen of TB ? and how fast
are the queries ? (for example a simple lookup)<br>
 </div>
    <blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
      <div text="#000000" bgcolor="#ffffff"><font color="#888888">-Peter</font>
      <div>
      <div><br>
      <br>
      <br>
      <blockquote type="cite">
        <div>
        <p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);">I
would suggest you
split your raster with something like GDAL gdal_retile.py and then
import your
tiles into a single WKT Raster table following the instruction provided
in our
tutorial. Each tile will be stored in a column cell of type RASTER
similar to
the PostGIS GEOMETRY type.</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);">As
I said previously you can then query the raster in SQL like
this:</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);">SELECT
ST_Value(rast, ST_Geomfromtext('Point(-78.1 58.1)',
4326))</span></p>
        <p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);">FROM
srtm_tiled_100x100</span></p>
        <p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);">WHERE
ST_Intersects(rast::geometry, ST_Geomfromtext('Point(-78.1
58.1)', 4326)) AND whatever you want.</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);">Basically
WKT Raster is the first true SQL interface with which is
is simple to do such things… It's proven, any-size, cloud-scalable, and
open source </span><span style="font-size: 11pt; font-family: Wingdings; color: rgb(31, 73, 125);">J</span><span style="font-size: 11pt; color: rgb(31, 73, 125);">
Maybe not that prooven. But who is?</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);">There
is also rasdaman but I don’t think you can use SQL.
It would be nice if you could compare both… I can’t find time for
this. Jorge has started comparing PostGIS WKT Raster with Oracle
Georaster
though.</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>
        <p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);"> </span></p>
        <div style="border-style: none none none solid; border-color: -moz-use-text-color -moz-use-text-color -moz-use-text-color blue; border-width: medium medium medium 1.5pt; padding: 0cm 0cm 0cm 4pt;">
        <div>
        <div style="border-style: solid none none; border-color: rgb(181, 196, 223) -moz-use-text-color -moz-use-text-color; border-width: 1pt medium medium; 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>
[<a href="mailto:postgis-users-bounces@postgis.refractions.net" target="_blank">mailto:postgis-users-bounces@postgis.refractions.net</a>]

        <b>On Behalf Of </b>Sebastian
E. Ovide<br>
        <b>Sent:</b> 20 juillet 2010 10:46<br>
        <b>To:</b> PostGIS Users Discussion<br>
        <b>Subject:</b> Re: [postgis-users] WKTRaster :
gdal2wktraster.py
cannot read
AIG/Arc/Info Binary Grid</span></p>
        </div>
        </div>
        <p class="MsoNormal"> </p>
        <p class="MsoNormal" style="margin-bottom: 12pt;">yes... with
png
worked... but
it was a different png (a smaller one)...<br>
        <br>
so I've converted the ESRI into a png and tried to import it... and it
didn't
work neither... <br>
        <br>
so this is the situation:<br>
A have huge rusters (from 150kx150k).. In Oracle I would just load it
(the huge
raster) in a single row of a GeoRaster table and then Oracle GeoRaster
would
split it in small tiles and store one tile per line of another table
(Raster
Data Table)... then I can run a query similar to this one: SELECT
getcellvalue(t.rastercolumn,x,y) from GeoRasterTable t where <a href="http://t.id" target="_blank">t.id</a>=1;
and GeoRaster
will query automatically the
spatial indexes and the Raster Data Table and it will find the right
tile
etc...<br>
        <br>
I do not know how WKTRaster works.... If I cannot import a such big
image, of
course I can split it in smaller georeferenced tiles... (how?)...
but.... My
main question is: after that, how will the table look like ? how can I
do the
same query (where id=1 or where name="UK" etc...) ?<br>
        <br>
        </p>
        <div>
        <p class="MsoNormal">On Tue, Jul 20, 2010 at 2:31 PM, Pierre
Racine
<<a href="mailto:Pierre.Racine@sbf.ulaval.ca" target="_blank">Pierre.Racine@sbf.ulaval.ca</a>>
wrote:</p>
        <div>
        <div>
        <p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);">Wait. You first said
the png was
working. Now it’s not? Did you try gdal_translate with the ESRI grid?
For
sure I haven’t test yet with such big rasters. Is this the result of a
merge or all your original raster are all this size? The point is that
with WKT
Raster you don’t have to merge your raster first into a gigantic raster
in order to get it store in a unique table like with Oracle Spatial.</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-style: none none none solid; border-color: -moz-use-text-color -moz-use-text-color -moz-use-text-color blue; border-width: medium medium medium 1.5pt; padding: 0cm 0cm 0cm 4pt;">
        <div>
        <div style="border-style: solid none none; border-color: -moz-use-text-color; border-width: 1pt medium medium; 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> 20 juillet 2010 05:51<br>
        <b>To:</b> PostGIS Users Discussion<br>
        <b>Subject:</b> Re: [postgis-users] WKTRaster :
gdal2wktraster.py
cannot read
AIG/Arc/Info Binary Grid</span></p>
        </div>
        </div>
        <div>
        <div>
        <p class="MsoNormal"> </p>
        <p class="MsoNormal" style="margin-bottom: 12pt;">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>
        </p>
        <div>
        <p class="MsoNormal">On
Mon, Jul 19, 2010 at 5:56 PM, Pierre Racine <<a href="mailto:Pierre.Racine@sbf.ulaval.ca" target="_blank">Pierre.Racine@sbf.ulaval.ca</a>>
wrote:</p>
        <div>
        <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-style: none none none solid; border-color: -moz-use-text-color -moz-use-text-color -moz-use-text-color blue; border-width: medium medium medium 1.5pt; padding: 0cm 0cm 0cm 4pt;">
        <div>
        <div style="border-style: solid none none; border-color: -moz-use-text-color; border-width: 1pt medium medium; 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>
        <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</p>
        </div>
        </div>
        </div>
        </div>
        </div>
        <p class="MsoNormal" style="margin-bottom: 12pt;"><br>
_______________________________________________<br>
postgis-users mailing list<br>
        <a href="mailto:postgis-users@postgis.refractions.net" target="_blank">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></p>
        </div>
        <p class="MsoNormal" style="margin-bottom: 12pt;"><br>
        <br clear="all">
        <br>
-- <br>
Sebastian E. Ovide<br>
        <br>
skype: sebastian.ovide<br>
        <br>
+353 (0) 87 6340149</p>
        </div>
        </div>
        </div>
        </div>
        </div>
        <p class="MsoNormal" style="margin-bottom: 12pt;"><br>
_______________________________________________<br>
postgis-users mailing list<br>
        <a href="mailto:postgis-users@postgis.refractions.net" target="_blank">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></p>
        </div>
        <p class="MsoNormal" style="margin-bottom: 12pt;"><br>
        <br clear="all">
        <br>
-- <br>
Sebastian E. Ovide<br>
        <br>
        <br>
        <br>
        </p>
        </div>
        </div>
      </blockquote>
      <br>
      </div>
      </div>
      <div>
      <pre cols="80">-- 
Dr. Peter Baumann
 - Professor of Computer Science, Jacobs University Bremen
   <a href="http://www.faculty.jacobs-university.de/pbaumann" target="_blank">www.faculty.jacobs-university.de/pbaumann</a>
   mail: <a href="mailto:p.baumann@jacobs-university.de" target="_blank">p.baumann@jacobs-university.de</a>
   tel: +49-421-200-3178, fax: +49-421-200-493178
 - Executive Director, rasdaman GmbH Bremen (HRB 147737)
   <a href="http://www.rasdaman.com" target="_blank">www.rasdaman.com</a>, mail: <a href="mailto:baumann@rasdaman.com" target="_blank">baumann@rasdaman.com</a>
   tel: 0800-rasdaman, fax: 0800-rasdafax, mobile: +49-173-5837882
"Si forte in alienas manus oberraverit hec peregrina epistola incertis ventis dimissa, sed Deo commendata, precamur ut ei reddatur cui soli destinata, nec preripiat quisquam non sibi parata." (mail disclaimer, AD 10xx)


    </pre>
      </div>
      </div>
      <br>
_______________________________________________<br>
postgis-users mailing list<br>
      <a href="mailto:postgis-users@postgis.refractions.net" target="_blank">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>
    <br>
    <br>
  </blockquote>
  <br>
  <pre cols="80">-- 
Dr. Peter Baumann
 - Professor of Computer Science, Jacobs University Bremen
   <a href="http://www.faculty.jacobs-university.de/pbaumann" target="_blank">www.faculty.jacobs-university.de/pbaumann</a>
   mail: <a href="mailto:p.baumann@jacobs-university.de" target="_blank">p.baumann@jacobs-university.de</a>
   tel: +49-421-200-3178, fax: +49-421-200-493178
 - Executive Director, rasdaman GmbH Bremen (HRB 147737)
   <a href="http://www.rasdaman.com" target="_blank">www.rasdaman.com</a>, mail: <a href="mailto:baumann@rasdaman.com" target="_blank">baumann@rasdaman.com</a>
   tel: 0800-rasdaman, fax: 0800-rasdafax, mobile: +49-173-5837882
"Si forte in alienas manus oberraverit hec peregrina epistola incertis ventis dimissa, sed Deo commendata, precamur ut ei reddatur cui soli destinata, nec preripiat quisquam non sibi parata." (mail disclaimer, AD 10xx)


  </pre>
</blockquote>
<br>
<pre cols="80">-- 
Dr. Peter Baumann
 - Professor of Computer Science, Jacobs University Bremen
   <a href="http://www.faculty.jacobs-university.de/pbaumann" target="_blank">www.faculty.jacobs-university.de/pbaumann</a>
   mail: <a href="mailto:p.baumann@jacobs-university.de" target="_blank">p.baumann@jacobs-university.de</a>
   tel: +49-421-200-3178, fax: +49-421-200-493178
 - Executive Director, rasdaman GmbH Bremen (HRB 147737)
   <a href="http://www.rasdaman.com" target="_blank">www.rasdaman.com</a>, mail: <a href="mailto:baumann@rasdaman.com" target="_blank">baumann@rasdaman.com</a>
   tel: 0800-rasdaman, fax: 0800-rasdafax, mobile: +49-173-5837882
"Si forte in alienas manus oberraverit hec peregrina epistola incertis ventis dimissa, sed Deo commendata, precamur ut ei reddatur cui soli destinata, nec preripiat quisquam non sibi parata." (mail disclaimer, AD 10xx)


</pre>
</div>

</blockquote></div><br><br clear="all"><br>-- <br>Sebastian E. Ovide<br><br><br>