<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">

<head>
<meta http-equiv=Content-Type content="text/html; charset=us-ascii">
<meta name=Generator content="Microsoft Word 12 (filtered medium)">
<style>
<!--
 /* Font Definitions */
 @font-face
        {font-family:Wingdings;
        panose-1:5 0 0 0 0 0 0 0 0 0;}
@font-face
        {font-family:Wingdings;
        panose-1:5 0 0 0 0 0 0 0 0 0;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
        {font-family:Tahoma;
        panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
        {font-family:Consolas;
        panose-1:2 11 6 9 2 2 4 3 2 4;}
 /* Style Definitions */
 p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0cm;
        margin-bottom:.0001pt;
        font-size:12.0pt;
        font-family:"Times New Roman","serif";
        color:black;}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:purple;
        text-decoration:underline;}
pre
        {mso-style-priority:99;
        mso-style-link:"HTML Preformatted Char";
        margin:0cm;
        margin-bottom:.0001pt;
        font-size:10.0pt;
        font-family:"Courier New";
        color:black;}
span.HTMLPreformattedChar
        {mso-style-name:"HTML Preformatted Char";
        mso-style-priority:99;
        mso-style-link:"HTML Preformatted";
        font-family:Consolas;
        color:black;}
span.EmailStyle19
        {mso-style-type:personal;
        font-family:"Calibri","sans-serif";
        color:#1F497D;}
span.EmailStyle20
        {mso-style-type:personal;
        font-family:"Calibri","sans-serif";
        color:#1F497D;}
span.EmailStyle21
        {mso-style-type:personal;
        font-family:"Calibri","sans-serif";
        color:#1F497D;}
span.EmailStyle22
        {mso-style-type:personal;
        font-family:"Calibri","sans-serif";
        color:#1F497D;}
span.EmailStyle23
        {mso-style-type:personal-reply;
        font-family:"Calibri","sans-serif";
        color:#1F497D;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-size:10.0pt;}
@page WordSection1
        {size:612.0pt 792.0pt;
        margin:70.85pt 70.85pt 70.85pt 70.85pt;}
div.WordSection1
        {page:WordSection1;}
-->
</style>
<!--[if gte mso 9]><xml>
 <o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
 <o:shapelayout v:ext="edit">
  <o:idmap v:ext="edit" data="1" />
 </o:shapelayout></xml><![endif]-->
</head>

<body bgcolor=white lang=EN-US link=blue vlink=purple>

<div class=WordSection1>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>A table of tiles form a large coverage as one. Every process
might be done one tile at a time but SQL GROUP BY statements help regrouping
things together if needed.<o:p></o:p></span></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'><o:p> </o:p></span></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>Pierre<o:p></o:p></span></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'><o:p> </o:p></span></p>

<div style='border:none;border-left:solid blue 1.5pt;padding:0cm 0cm 0cm 4.0pt'>

<div>

<div style='border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm'>

<p class=MsoNormal><b><span style='font-size:10.0pt;font-family:"Tahoma","sans-serif";
color:windowtext'>From:</span></b><span style='font-size:10.0pt;font-family:
"Tahoma","sans-serif";color:windowtext'> Peter Baumann
[mailto:p.baumann@jacobs-university.de] <br>
<b>Sent:</b> 28 juillet 2010 04:27<br>
<b>To:</b> Pierre Racine<br>
<b>Cc:</b> PostGIS Users Discussion<br>
<b>Subject:</b> Re: [postgis-users] WKTRaster : gdal2wktraster.py cannot read
AIG/Arc/Info Binary Grid<o:p></o:p></span></p>

</div>

</div>

<p class=MsoNormal><o:p> </o:p></p>

<p class=MsoNormal>Pierre,<br>
<br>
admittedly this was a little inexact, and I certainly will correct it. A more
accurate phrasing I guess is "does not support retrieval on a raster array
consisting of a set of tiles, rather tiles have to be treated
individually."<br>
>From what I understand WKT Raster can tile an incoming image, but cannot treat
a set of tiles making up a larger coverage as one.<br>
<br>
Would you agree to that?<br>
<br>
-Peter<br>
<br>
<br>
<br>
On 07/27/2010 05:52 PM, Pierre Racine wrote: <o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>BTW Peter, </span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'> </span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>in your site (<a
href="http://kahlua.eecs.jacobs-university.de/trac/rasdaman/wiki/WhyRasdaman">http://kahlua.eecs.jacobs-university.de/trac/rasdaman/wiki/WhyRasdaman</a>)
you say: “First and foremost, it (PostGIS WKT Raster) does not support
tiling (partitioning) of raster arrays.” This statement is false since
WKT Raster actually does support tiling of large rasters. You can partition any
raster coverage by importing it as separate tiles or by using the –k
option of gdal2wktraster.py. Each tile become a separate row in the table and
the whole table is the coverage. The table is spatially indexed using the
PostgreSQL GiST index.</span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'> </span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>I think it would be just fair if you could rectify this
statement.</span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'> </span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>Thanks,</span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'> </span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>Pierre</span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'> </span><o:p></o:p></p>

<div style='border:none;border-left:solid windowtext 1.5pt;padding:0cm 0cm 0cm 4.0pt;
border-color:-moz-use-text-color -moz-use-text-color -moz-use-text-color blue'>

<div>

<div style='border:none;border-top:solid windowtext 1.0pt;padding:3.0pt 0cm 0cm 0cm;
border-color:-moz-use-text-color -moz-use-text-color'>

<p class=MsoNormal><b><span style='font-size:10.0pt;font-family:"Tahoma","sans-serif";
color:windowtext'>From:</span></b><span style='font-size:10.0pt;font-family:
"Tahoma","sans-serif";color:windowtext'> Peter Baumann [<a
href="mailto:p.baumann@jacobs-university.de">mailto:p.baumann@jacobs-university.de</a>]
<br>
<b>Sent:</b> 21 juillet 2010 17:28<br>
<b>To:</b> Pierre Racine<br>
<b>Cc:</b> PostGIS Users Discussion<br>
<b>Subject:</b> Re: [postgis-users] WKTRaster : gdal2wktraster.py cannot read
AIG/Arc/Info Binary Grid</span><o:p></o:p></p>

</div>

</div>

<p class=MsoNormal> <o:p></o:p></p>

<p class=MsoNormal>Pierre-<br>
<br>
rasdaman per se is domain agnostic - you can just as well maintain human brain
images (as we did). The petascope layer, which is part of the distro, adds geo
semantics by providing OGC interfaces.<br>
PostgreSQL integration is not yet done, we are working on this; same for GDAL /
MapServer.<br>
<br>
-Peter<br>
<br>
<br>
On 07/21/2010 08:29 PM, Pierre Racine wrote: <o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>Peter,</span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'> </span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>Are rasters in rasdaman georeferenced? How do I convert my
vectorial PostGIS points coordinates into rasdaman coordinates? Do I have to do
all the math myself?</span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'> </span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>Pierre</span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'> </span><o:p></o:p></p>

<div style='border:none;border-left:solid windowtext 1.5pt;padding:0cm 0cm 0cm 4.0pt;
border-color:-moz-use-text-color -moz-use-text-color -moz-use-text-color blue'>

<div>

<div style='border:none;border-top:solid windowtext 1.0pt;padding:3.0pt 0cm 0cm 0cm;
border-color:-moz-use-text-color'>

<p class=MsoNormal><b><span style='font-size:10.0pt;font-family:"Tahoma","sans-serif";
color:windowtext'>From:</span></b><span style='font-size:10.0pt;font-family:
"Tahoma","sans-serif";color:windowtext'> Peter Baumann [<a
href="mailto:p.baumann@jacobs-university.de">mailto:p.baumann@jacobs-university.de</a>]
<br>
<b>Sent:</b> 21 juillet 2010 04:41<br>
<b>To:</b> Pierre Racine<br>
<b>Subject:</b> Re: [postgis-users] WKTRaster : gdal2wktraster.py cannot read AIG/Arc/Info
Binary Grid</span><o:p></o:p></p>

</div>

</div>

<p class=MsoNormal> <o:p></o:p></p>

<p class=MsoNormal>Pierre-<br>
<br>
let's see. First, I would not store 100k raster points in the relational table
as its terribly inefficient. We do have input data formats, though, which are
similar to what a PG dump of such a table would deliver. We import this into an
array and then can query it. Import via rasql or through an import tool which
internally generates rasql calls.<br>
<br>
Then, for corner points (x0,y0) and (x1,y1) we retrieve as<br>
<br>
    select png( M[x0:y1,y0:y1] ) from MyMap as M<br>
<br>
...in PNG format via the rasql cmd line utility or the C++ or Java API.<br>
<br>
Alternatively, sending this this query over the API gets you the data into main
memory ready for nested for loops = in the main memory format of your client:<br>
<br>
    select M[x0:y1,y0:y1] from MyMap as M<br>
<br>
-Peter<br>
<br>
<br>
<br>
On 07/20/2010 10:03 PM, Pierre Racine wrote: <o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>So if I already have a table of,  say 100 000 points,
uploaded in PostGIS, how do I use rasdaman to get the raster values of the
pixels located under the points?</span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'> </span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>Pierre</span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'> </span><o:p></o:p></p>

<div style='border:none;border-left:solid windowtext 1.5pt;padding:0cm 0cm 0cm 4.0pt;
border-color:-moz-use-text-color -moz-use-text-color -moz-use-text-color blue'>

<div>

<div style='border:none;border-top:solid windowtext 1.0pt;padding:3.0pt 0cm 0cm 0cm;
border-color:-moz-use-text-color'>

<p class=MsoNormal><b><span style='font-size:10.0pt;font-family:"Tahoma","sans-serif";
color:windowtext'>From:</span></b><span style='font-size:10.0pt;font-family:
"Tahoma","sans-serif";color:windowtext'> Peter Baumann [<a
href="mailto:p.baumann@jacobs-university.de">mailto:p.baumann@jacobs-university.de</a>]
<br>
<b>Sent:</b> 20 juillet 2010 12:08<br>
<b>To:</b> PostGIS Users Discussion<br>
<b>Cc:</b> Pierre Racine<br>
<b>Subject:</b> Re: [postgis-users] WKTRaster : gdal2wktraster.py cannot read
AIG/Arc/Info Binary Grid</span><o:p></o:p></p>

</div>

</div>

<p class=MsoNormal> <o:p></o:p></p>

<p class=MsoNormal>On 07/20/2010 05:47 PM, Pierre Racine wrote: <o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>Since we have this bug with big rasters, </span><o:p></o:p></p>

<p class=MsoNormal><br>
proven, any-size, ... ;-)<br>
<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>
<br>
-Peter<br>
<br>
<br>
<br>
<br>
<br>
<br>
<o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>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><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'> </span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>As I said previously you can then query the raster in SQL like
this:</span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'> </span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>SELECT ST_Value(rast, ST_Geomfromtext('Point(-78.1 58.1)', 4326))</span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>FROM srtm_tiled_100x100</span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>WHERE ST_Intersects(rast::geometry, ST_Geomfromtext('Point(-78.1
58.1)', 4326)) AND whatever you want.</span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'> </span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>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:11.0pt;font-family:Wingdings;
color:#1F497D'>J</span><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'> Maybe not that prooven. But who is?</span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'> </span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>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><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'> </span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>Pierre</span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'> </span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'> </span><o:p></o:p></p>

<div style='border:none;border-left:solid windowtext 1.5pt;padding:0cm 0cm 0cm 4.0pt;
border-color:-moz-use-text-color -moz-use-text-color -moz-use-text-color blue'>

<div>

<div style='border:none;border-top:solid windowtext 1.0pt;padding:3.0pt 0cm 0cm 0cm;
border-color:-moz-use-text-color'>

<p class=MsoNormal><b><span style='font-size:10.0pt;font-family:"Tahoma","sans-serif"'>From:</span></b><span
style='font-size:10.0pt;font-family:"Tahoma","sans-serif"'> <a
href="mailto:postgis-users-bounces@postgis.refractions.net">postgis-users-bounces@postgis.refractions.net</a>
[<a href="mailto:postgis-users-bounces@postgis.refractions.net">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><o:p></o:p></p>

</div>

</div>

<p class=MsoNormal> <o:p></o:p></p>

<p class=MsoNormal style='margin-bottom:12.0pt'>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">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>
<br>
<br>
<br>
<br>
<o:p></o:p></p>

<div>

<p class=MsoNormal>On Tue, Jul 20, 2010 at 2:31 PM, Pierre Racine <<a
href="mailto:Pierre.Racine@sbf.ulaval.ca">Pierre.Racine@sbf.ulaval.ca</a>>
wrote:<o:p></o:p></p>

<div>

<div>

<p class=MsoNormal><span style='font-size:11.0pt;color:#1F497D'>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><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;color:#1F497D'> </span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;color:#1F497D'>Pierre</span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;color:#1F497D'> </span><o:p></o:p></p>

<div style='border:none;border-left:solid windowtext 1.5pt;padding:0cm 0cm 0cm 4.0pt;
border-color:-moz-use-text-color -moz-use-text-color -moz-use-text-color blue'>

<div>

<div style='border:none;border-top:solid windowtext 1.0pt;padding:3.0pt 0cm 0cm 0cm;
border-color:-moz-use-text-color'>

<p class=MsoNormal><b><span style='font-size:10.0pt'>From:</span></b><span
style='font-size:10.0pt'> <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><o:p></o:p></p>

</div>

</div>

<div>

<div>

<p class=MsoNormal> <o:p></o:p></p>

<p class=MsoNormal style='margin-bottom:12.0pt'>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>
<br>
<br>
<o:p></o:p></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:<o:p></o:p></p>

<div>

<div>

<p class=MsoNormal><span style='font-size:11.0pt;color:#1F497D'>Hi Sebastian,</span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;color:#1F497D'> </span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;color:#1F497D'>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><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;color:#1F497D'> </span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;color:#1F497D'>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><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;color:#1F497D'> </span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;color:#1F497D'>Could you
provide us with a file sample?</span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;color:#1F497D'> </span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;color:#1F497D'>Pierre</span><o:p></o:p></p>

<p class=MsoNormal><span style='font-size:11.0pt;color:#1F497D'> </span><o:p></o:p></p>

<div style='border:none;border-left:solid windowtext 1.5pt;padding:0cm 0cm 0cm 4.0pt;
border-color:-moz-use-text-color -moz-use-text-color -moz-use-text-color blue'>

<div>

<div style='border:none;border-top:solid windowtext 1.0pt;padding:3.0pt 0cm 0cm 0cm;
border-color:-moz-use-text-color'>

<p class=MsoNormal><b><span style='font-size:10.0pt'>From:</span></b><span
style='font-size:10.0pt'> <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><o:p></o:p></p>

</div>

</div>

<div>

<div>

<p class=MsoNormal> <o:p></o:p></p>

<p class=MsoNormal style='margin-bottom:12.0pt'>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<o:p></o:p></p>

</div>

</div>

</div>

</div>

</div>

<p class=MsoNormal style='margin-bottom:12.0pt'><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><o:p></o:p></p>

</div>

<p class=MsoNormal style='margin-bottom:12.0pt'><br>
<br clear=all>
<br>
-- <br>
Sebastian E. Ovide<br>
<br>
skype: sebastian.ovide<br>
<br>
+353 (0) 87 6340149<o:p></o:p></p>

</div>

</div>

</div>

</div>

</div>

<p class=MsoNormal style='margin-bottom:12.0pt'><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><o:p></o:p></p>

</div>

<p class=MsoNormal style='margin-bottom:12.0pt'><br>
<br clear=all>
<br>
-- <br>
Sebastian E. Ovide<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<o:p></o:p></p>

</div>

<p class=MsoNormal><br>
<br>
<br>
<br>
<br>
<o:p></o:p></p>

<pre>-- <o:p></o:p></pre><pre>Dr. Peter Baumann<o:p></o:p></pre><pre> - Professor of Computer Science, Jacobs University Bremen<o:p></o:p></pre><pre>   <a
href="http://www.faculty.jacobs-university.de/pbaumann">www.faculty.jacobs-university.de/pbaumann</a><o:p></o:p></pre><pre>   mail: <a
href="mailto:p.baumann@jacobs-university.de">p.baumann@jacobs-university.de</a><o:p></o:p></pre><pre>   tel: +49-421-200-3178, fax: +49-421-200-493178<o:p></o:p></pre><pre> - Executive Director, rasdaman GmbH Bremen (HRB 147737)<o:p></o:p></pre><pre>   <a
href="http://www.rasdaman.com">www.rasdaman.com</a>, mail: <a
href="mailto:baumann@rasdaman.com">baumann@rasdaman.com</a><o:p></o:p></pre><pre>   tel: 0800-rasdaman, fax: 0800-rasdafax, mobile: +49-173-5837882<o:p></o:p></pre><pre>"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)<o:p></o:p></pre><pre> <o:p></o:p></pre><pre> <o:p></o:p></pre></div>

<p class=MsoNormal><br>
<br>
<br>
<br>
<o:p></o:p></p>

<pre>-- <o:p></o:p></pre><pre>Dr. Peter Baumann<o:p></o:p></pre><pre> - Professor of Computer Science, Jacobs University Bremen<o:p></o:p></pre><pre>   <a
href="http://www.faculty.jacobs-university.de/pbaumann">www.faculty.jacobs-university.de/pbaumann</a><o:p></o:p></pre><pre>   mail: <a
href="mailto:p.baumann@jacobs-university.de">p.baumann@jacobs-university.de</a><o:p></o:p></pre><pre>   tel: +49-421-200-3178, fax: +49-421-200-493178<o:p></o:p></pre><pre> - Executive Director, rasdaman GmbH Bremen (HRB 147737)<o:p></o:p></pre><pre>   <a
href="http://www.rasdaman.com">www.rasdaman.com</a>, mail: <a
href="mailto:baumann@rasdaman.com">baumann@rasdaman.com</a><o:p></o:p></pre><pre>   tel: 0800-rasdaman, fax: 0800-rasdafax, mobile: +49-173-5837882<o:p></o:p></pre><pre>"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)<o:p></o:p></pre><pre> <o:p></o:p></pre><pre> <o:p></o:p></pre></div>

<p class=MsoNormal><br>
<br>
<br>
<o:p></o:p></p>

<pre>-- <o:p></o:p></pre><pre>Dr. Peter Baumann<o:p></o:p></pre><pre> - Professor of Computer Science, Jacobs University Bremen<o:p></o:p></pre><pre>   <a
href="http://www.faculty.jacobs-university.de/pbaumann">www.faculty.jacobs-university.de/pbaumann</a><o:p></o:p></pre><pre>   mail: <a
href="mailto:p.baumann@jacobs-university.de">p.baumann@jacobs-university.de</a><o:p></o:p></pre><pre>   tel: +49-421-200-3178, fax: +49-421-200-493178<o:p></o:p></pre><pre> - Executive Director, rasdaman GmbH Bremen (HRB 147737)<o:p></o:p></pre><pre>   <a
href="http://www.rasdaman.com">www.rasdaman.com</a>, mail: <a
href="mailto:baumann@rasdaman.com">baumann@rasdaman.com</a><o:p></o:p></pre><pre>   tel: 0800-rasdaman, fax: 0800-rasdafax, mobile: +49-173-5837882<o:p></o:p></pre><pre>"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)<o:p></o:p></pre><pre> <o:p></o:p></pre><pre> <o:p></o:p></pre></div>

<p class=MsoNormal><br>
<br>
<o:p></o:p></p>

<pre>-- <o:p></o:p></pre><pre>Dr. Peter Baumann<o:p></o:p></pre><pre> - Professor of Computer Science, Jacobs University Bremen<o:p></o:p></pre><pre>   <a
href="http://www.faculty.jacobs-university.de/pbaumann">www.faculty.jacobs-university.de/pbaumann</a><o:p></o:p></pre><pre>   mail: <a
href="mailto:p.baumann@jacobs-university.de">p.baumann@jacobs-university.de</a><o:p></o:p></pre><pre>   tel: +49-421-200-3178, fax: +49-421-200-493178<o:p></o:p></pre><pre> - Executive Director, rasdaman GmbH Bremen (HRB 147737)<o:p></o:p></pre><pre>   <a
href="http://www.rasdaman.com">www.rasdaman.com</a>, mail: <a
href="mailto:baumann@rasdaman.com">baumann@rasdaman.com</a><o:p></o:p></pre><pre>   tel: 0800-rasdaman, fax: 0800-rasdafax, mobile: +49-173-5837882<o:p></o:p></pre><pre>"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)<o:p></o:p></pre><pre><o:p> </o:p></pre><pre><o:p> </o:p></pre></div>

</div>

</body>

</html>