No subject
Wed Jul 15 09:27:57 PDT 2009
eat
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=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'>BTW Peter, </span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'> </span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'>in your site (<a
href=3D"http://kahlua.eecs.jacobs-university.de/trac/rasdaman/wiki/WhyRasda=
man">http://kahlua.eecs.jacobs-university.de/trac/rasdaman/wiki/WhyRasdaman=
</a>)
you say: “First and foremost, it (PostGIS WKT Raster) does not suppor=
t
tiling (partitioning) of raster arrays.” This statement is false sinc=
e
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 a=
nd
the whole table is the coverage. The table is spatially indexed using the
PostgreSQL GiST index.</span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'> </span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'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=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'> </span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'>Thanks,</span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'> </span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'>Pierre</span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'> </span><o:p></o:p></p>
<div style=3D'border:none;border-left:solid windowtext 1.5pt;padding:0cm 0c=
m 0cm 4.0pt;
border-color:-moz-use-text-color -moz-use-text-color -moz-use-text-color bl=
ue'>
<div>
<div style=3D'border:none;border-top:solid windowtext 1.0pt;padding:3.0pt 0=
cm 0cm 0cm;
border-color:-moz-use-text-color -moz-use-text-color'>
<p class=3DMsoNormal><b><span style=3D'font-size:10.0pt;font-family:"Tahoma=
","sans-serif";
color:windowtext'>From:</span></b><span style=3D'font-size:10.0pt;font-fami=
ly:
"Tahoma","sans-serif";color:windowtext'> Peter Baumann [<a
href=3D"mailto:p.baumann at jacobs-university.de">mailto:p.baumann at jacobs-univ=
ersity.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 re=
ad
AIG/Arc/Info Binary Grid</span><o:p></o:p></p>
</div>
</div>
<p class=3DMsoNormal> <o:p></o:p></p>
<p class=3DMsoNormal>Pierre-<br>
<br>
rasdaman per se is domain agnostic - you can just as well maintain human br=
ain
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 GD=
AL /
MapServer.<br>
<br>
-Peter<br>
<br>
<br>
On 07/21/2010 08:29 PM, Pierre Racine wrote: <o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'>Peter,</span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'> </span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'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 t=
o do
all the math myself?</span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'> </span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'>Pierre</span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'> </span><o:p></o:p></p>
<div style=3D'border:none;border-left:solid windowtext 1.5pt;padding:0cm 0c=
m 0cm 4.0pt;
border-color:-moz-use-text-color -moz-use-text-color -moz-use-text-color bl=
ue'>
<div>
<div style=3D'border:none;border-top:solid windowtext 1.0pt;padding:3.0pt 0=
cm 0cm 0cm;
border-color:-moz-use-text-color'>
<p class=3DMsoNormal><b><span style=3D'font-size:10.0pt;font-family:"Tahoma=
","sans-serif";
color:windowtext'>From:</span></b><span style=3D'font-size:10.0pt;font-fami=
ly:
"Tahoma","sans-serif";color:windowtext'> Peter Baumann [<a
href=3D"mailto:p.baumann at jacobs-university.de">mailto:p.baumann at jacobs-univ=
ersity.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 re=
ad AIG/Arc/Info
Binary Grid</span><o:p></o:p></p>
</div>
</div>
<p class=3DMsoNormal> <o:p></o:p></p>
<p class=3DMsoNormal>Pierre-<br>
<br>
let's see. First, I would not store 100k raster points in the relational ta=
ble
as its terribly inefficient. We do have input data formats, though, which a=
re
similar to what a PG dump of such a table would deliver. We import this int=
o an
array and then can query it. Import via rasql or through an import tool whi=
ch
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 =3D in the main memory format of your cli=
ent:<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=3DMsoNormal><span style=3D'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=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'> </span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'>Pierre</span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'> </span><o:p></o:p></p>
<div style=3D'border:none;border-left:solid windowtext 1.5pt;padding:0cm 0c=
m 0cm 4.0pt;
border-color:-moz-use-text-color -moz-use-text-color -moz-use-text-color bl=
ue'>
<div>
<div style=3D'border:none;border-top:solid windowtext 1.0pt;padding:3.0pt 0=
cm 0cm 0cm;
border-color:-moz-use-text-color'>
<p class=3DMsoNormal><b><span style=3D'font-size:10.0pt;font-family:"Tahoma=
","sans-serif";
color:windowtext'>From:</span></b><span style=3D'font-size:10.0pt;font-fami=
ly:
"Tahoma","sans-serif";color:windowtext'> Peter Baumann [<a
href=3D"mailto:p.baumann at jacobs-university.de">mailto:p.baumann at jacobs-univ=
ersity.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 re=
ad
AIG/Arc/Info Binary Grid</span><o:p></o:p></p>
</div>
</div>
<p class=3DMsoNormal> <o:p></o:p></p>
<p class=3DMsoNormal>On 07/20/2010 05:47 PM, Pierre Racine wrote: <o:p></o:=
p></p>
<p class=3DMsoNormal><span style=3D'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=3DMsoNormal><br>
proven, any-size, ... ;-)<br>
<br>
Who is proven? Well, running a dozen-TB seamless mosaic on PostgreSQL for m=
any
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=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'>I would suggest you split your raster with something like GD=
AL
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 store=
d in
a column cell of type RASTER similar to the PostGIS GEOMETRY type.</span><o=
:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'> </span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'>As I said previously you can then query the raster in SQL li=
ke
this:</span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'> </span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'>SELECT ST_Value(rast, ST_Geomfromtext('Point(-78.1 58.1)', 4=
326))</span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'>FROM srtm_tiled_100x100</span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'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=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'> </span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'>Basically WKT Raster is the first true SQL interface with wh=
ich
is is simple to do such things… It's proven, any-size, cloud-scalable=
,
and open source </span><span style=3D'font-size:11.0pt;font-family:Wingding=
s;
color:#1F497D'>J</span><span style=3D'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=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'> </span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'>There is also rasdaman but I don’t think you can use S=
QL.
It would be nice if you could compare both… I can’t find time f=
or
this. Jorge has started comparing PostGIS WKT Raster with Oracle Georaster
though.</span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'> </span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'>Pierre</span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'> </span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;font-family:"Calibri",=
"sans-serif";
color:#1F497D'> </span><o:p></o:p></p>
<div style=3D'border:none;border-left:solid windowtext 1.5pt;padding:0cm 0c=
m 0cm 4.0pt;
border-color:-moz-use-text-color -moz-use-text-color -moz-use-text-color bl=
ue'>
<div>
<div style=3D'border:none;border-top:solid windowtext 1.0pt;padding:3.0pt 0=
cm 0cm 0cm;
border-color:-moz-use-text-color'>
<p class=3DMsoNormal><b><span style=3D'font-size:10.0pt;font-family:"Tahoma=
","sans-serif"'>From:</span></b><span
style=3D'font-size:10.0pt;font-family:"Tahoma","sans-serif"'> <a
href=3D"mailto:postgis-users-bounces at postgis.refractions.net">postgis-users=
-bounces at postgis.refractions.net</a>
[<a href=3D"mailto:postgis-users-bounces at postgis.refractions.net">mailto:po=
stgis-users-bounces at 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 re=
ad
AIG/Arc/Info Binary Grid</span><o:p></o:p></p>
</div>
</div>
<p class=3DMsoNormal> <o:p></o:p></p>
<p class=3DMsoNormal style=3D'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 woul=
d
split it in small tiles and store one tile per line of another table (Raste=
r
Data Table)... then I can run a query similar to this one: SELECT
getcellvalue(t.rastercolumn,x,y) from GeoRasterTable t where <a
href=3D"http://t.id">t.id</a>=3D1; and GeoRaster will query automatically t=
he
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.... M=
y
main question is: after that, how will the table look like ? how can I do t=
he
same query (where id=3D1 or where name=3D"UK" etc...) ?<br>
<br>
<br>
<br>
<br>
<br>
<o:p></o:p></p>
<div>
<p class=3DMsoNormal>On Tue, Jul 20, 2010 at 2:31 PM, Pierre Racine <<a
href=3D"mailto:Pierre.Racine at sbf.ulaval.ca">Pierre.Racine at sbf.ulaval.ca</a>=
>
wrote:<o:p></o:p></p>
<div>
<div>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;color:#1F497D'>Wait. Y=
ou first
said the png was working. Now it’s not? Did you try gdal_translate wi=
th
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 po=
int
is that with WKT Raster you don’t have to merge your raster first int=
o 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=3DMsoNormal><span style=3D'font-size:11.0pt;color:#1F497D'> <=
/span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;color:#1F497D'>Pierre<=
/span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;color:#1F497D'> <=
/span><o:p></o:p></p>
<div style=3D'border:none;border-left:solid windowtext 1.5pt;padding:0cm 0c=
m 0cm 4.0pt;
border-color:-moz-use-text-color -moz-use-text-color -moz-use-text-color bl=
ue'>
<div>
<div style=3D'border:none;border-top:solid windowtext 1.0pt;padding:3.0pt 0=
cm 0cm 0cm;
border-color:-moz-use-text-color'>
<p class=3DMsoNormal><b><span style=3D'font-size:10.0pt'>From:</span></b><s=
pan
style=3D'font-size:10.0pt'> <a
href=3D"mailto:postgis-users-bounces at postgis.refractions.net" target=3D"_bl=
ank">postgis-users-bounces at postgis.refractions.net</a>
[mailto:<a href=3D"mailto:postgis-users-bounces at postgis.refractions.net"
target=3D"_blank">postgis-users-bounces at 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 re=
ad
AIG/Arc/Info Binary Grid</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class=3DMsoNormal> <o:p></o:p></p>
<p class=3DMsoNormal style=3D'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 at 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 at SeansPC:~/rasters$ python gdal2wktraster.py -r raster/ -t sebastable =
-o
ok.sql<br>
gdal2wktraster.py:695: DeprecationWarning: 'H' format requires 0 <=3D nu=
mber
<=3D 65535<br>
hexwkb +=3D wkblify('H', xsize)<br>
gdal2wktraster.py:696: DeprecationWarning: 'H' format requires 0 <=3D nu=
mber
<=3D 65535<br>
hexwkb +=3D wkblify('H', ysize)<br>
gdal2wktraster.py:727: DeprecationWarning: integer argument expected, got f=
loat<br>
hexwkb +=3D 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 =3D wkblify_raster_level(options, ds,
options.overview_level, band_range, infile, i)<br>
File "gdal2wktraster.py", line 888, in wkblify_raster_leve=
l<br>
hexwkb +=3D wkblify_band(options, band, level, xoff, yof=
f,
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", lin=
e
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 =3D numpy.empty([buf_ysize,buf_xsize], dtype =3D type=
code)<br>
MemoryError<br>
<br>
<br>
sebas at SeansPC:~/rasters$ python gdal2wktraster.py -r test.png -t
sebastable -o ok.sql<br>
gdal2wktraster.py:695: DeprecationWarning: 'H' format requires 0 <=3D nu=
mber
<=3D 65535<br>
hexwkb +=3D wkblify('H', xsize)<br>
gdal2wktraster.py:696: DeprecationWarning: 'H' format requires 0 <=3D nu=
mber
<=3D 65535<br>
hexwkb +=3D wkblify('H', ysize)<br>
gdal2wktraster.py:727: DeprecationWarning: integer argument expected, got f=
loat<br>
hexwkb +=3D 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 =3D wkblify_raster_level(options, ds,
options.overview_level, band_range, infile, i)<br>
File "gdal2wktraster.py", line 888, in wkblify_raster_leve=
l<br>
hexwkb +=3D wkblify_band(options, band, level, xoff, yof=
f,
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", lin=
e
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 =3D numpy.empty([buf_ysize,buf_xsize], dtype =3D type=
code)<br>
MemoryError<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<o:p></o:p></p>
<div>
<p class=3DMsoNormal>On Mon, Jul 19, 2010 at 5:56 PM, Pierre Racine <<a
href=3D"mailto:Pierre.Racine at sbf.ulaval.ca" target=3D"_blank">Pierre.Racine=
@sbf.ulaval.ca</a>>
wrote:<o:p></o:p></p>
<div>
<div>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;color:#1F497D'>Hi Seba=
stian,</span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;color:#1F497D'> <=
/span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;color:#1F497D'>I can c=
onvert
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=3DMsoNormal><span style=3D'font-size:11.0pt;color:#1F497D'> <=
/span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;color:#1F497D'>Maybe t=
his 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=3DMsoNormal><span style=3D'font-size:11.0pt;color:#1F497D'> <=
/span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;color:#1F497D'>Could y=
ou
provide us with a file sample?</span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;color:#1F497D'> <=
/span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;color:#1F497D'>Pierre<=
/span><o:p></o:p></p>
<p class=3DMsoNormal><span style=3D'font-size:11.0pt;color:#1F497D'> <=
/span><o:p></o:p></p>
<div style=3D'border:none;border-left:solid windowtext 1.5pt;padding:0cm 0c=
m 0cm 4.0pt;
border-color:-moz-use-text-color -moz-use-text-color -moz-use-text-color bl=
ue'>
<div>
<div style=3D'border:none;border-top:solid windowtext 1.0pt;padding:3.0pt 0=
cm 0cm 0cm;
border-color:-moz-use-text-color'>
<p class=3DMsoNormal><b><span style=3D'font-size:10.0pt'>From:</span></b><s=
pan
style=3D'font-size:10.0pt'> <a
href=3D"mailto:postgis-users-bounces at postgis.refractions.net" target=3D"_bl=
ank">postgis-users-bounces at postgis.refractions.net</a>
[mailto:<a href=3D"mailto:postgis-users-bounces at postgis.refractions.net"
target=3D"_blank">postgis-users-bounces at 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=3DMsoNormal> <o:p></o:p></p>
<p class=3DMsoNormal style=3D'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>
&nb=
sp;
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 =3D (128110.000000000000000,813270.000000000000000)<br>
Pixel Size =3D (5.000000000000000,-5.000000000000000)<br>
Corner Coordinates:<br>
Upper Left ( 128110.000, 813270.000) ( 6d29'37.32&q=
uot;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=3D256x4 Type=3DByte, ColorInterp=3DUndefined<br>
Min=3D1.000 Max=3D4.000<br>
NoData Value=3D255<br clear=3Dall>
<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 <=3D nu=
mber
<=3D 65535<br>
hexstr =3D binascii.hexlify(struct.pack(fmt_little, data)).upper()<b=
r>
gdal2wktraster.py:644: DeprecationWarning: integer argument expected, got f=
loat<br>
hexstr =3D binascii.hexlify(struct.pack(fmt_little, data)).upper()<b=
r>
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 =3D wkblify_raster_level(options, ds,
options.overview_level, band_range, infile, i)<br>
File "gdal2wktraster.py", line 888, in wkblify_raster_leve=
l<br>
hexwkb +=3D wkblify_band(options, band, level, xoff, yof=
f,
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 8=
35,
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 =3D numpy.reshape(ar, [buf_ysize,buf_xsize])<br>
File "C:\OSGeo4W\apps\Python25\lib\site-packages\numpy\core\fro=
mnumeric.py",
line 116, in reshape<br>
return reshape(newshape, order=3Dorder)<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=3DMsoNormal style=3D'margin-bottom:12.0pt'><br>
_______________________________________________<br>
postgis-users mailing list<br>
<a href=3D"mailto:postgis-users at postgis.refractions.net" target=3D"_blank">=
postgis-users at postgis.refractions.net</a><br>
<a href=3D"http://postgis.refractions.net/mailman/listinfo/postgis-users"
target=3D"_blank">http://postgis.refractions.net/mailman/listinfo/postgis-u=
sers</a><o:p></o:p></p>
</div>
<p class=3DMsoNormal style=3D'margin-bottom:12.0pt'><br>
<br clear=3Dall>
<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=3DMsoNormal style=3D'margin-bottom:12.0pt'><br>
_______________________________________________<br>
postgis-users mailing list<br>
<a href=3D"mailto:postgis-users at postgis.refractions.net">postgis-users at post=
gis.refractions.net</a><br>
<a href=3D"http://postgis.refractions.net/mailman/listinfo/postgis-users"
target=3D"_blank">http://postgis.refractions.net/mailman/listinfo/postgis-u=
sers</a><o:p></o:p></p>
</div>
<p class=3DMsoNormal style=3D'margin-bottom:12.0pt'><br>
<br clear=3Dall>
<br>
-- <br>
Sebastian E. Ovide<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<o:p></o:p></p>
</div>
<p class=3DMsoNormal><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> - Pro=
fessor of Computer Science, Jacobs University Bremen<o:p></o:p></pre><pre>&=
nbsp; <a
href=3D"http://www.faculty.jacobs-university.de/pbaumann">www.faculty.jacob=
s-university.de/pbaumann</a><o:p></o:p></pre><pre> mail: <a
href=3D"mailto:p.baumann at jacobs-university.de">p.baumann at jacobs-university.=
de</a><o:p></o:p></pre><pre> tel: +49-421-200-3178, fax: +49-42=
1-200-493178<o:p></o:p></pre><pre> - Executive Director, rasdaman GmbH Brem=
en (HRB 147737)<o:p></o:p></pre><pre> <a
href=3D"http://www.rasdaman.com">www.rasdaman.com</a>, mail: <a
href=3D"mailto:baumann at rasdaman.com">baumann at rasdaman.com</a><o:p></o:p></p=
re><pre> tel: 0800-rasdaman, fax: 0800-rasdafax, mobile: +49-17=
3-5837882<o:p></o:p></pre><pre>"Si forte in alienas manus oberraverit =
hec peregrina epistola incertis ventis dimissa, sed Deo commendata, precamu=
r 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=3DMsoNormal><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> - Pro=
fessor of Computer Science, Jacobs University Bremen<o:p></o:p></pre><pre>&=
nbsp; <a
href=3D"http://www.faculty.jacobs-university.de/pbaumann">www.faculty.jacob=
s-university.de/pbaumann</a><o:p></o:p></pre><pre> mail: <a
href=3D"mailto:p.baumann at jacobs-university.de">p.baumann at jacobs-university.=
de</a><o:p></o:p></pre><pre> tel: +49-421-200-3178, fax: +49-42=
1-200-493178<o:p></o:p></pre><pre> - Executive Director, rasdaman GmbH Brem=
en (HRB 147737)<o:p></o:p></pre><pre> <a
href=3D"http://www.rasdaman.com">www.rasdaman.com</a>, mail: <a
href=3D"mailto:baumann at rasdaman.com">baumann at rasdaman.com</a><o:p></o:p></p=
re><pre> tel: 0800-rasdaman, fax: 0800-rasdafax, mobile: +49-17=
3-5837882<o:p></o:p></pre><pre>"Si forte in alienas manus oberraverit =
hec peregrina epistola incertis ventis dimissa, sed Deo commendata, precamu=
r 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=3DMsoNormal><br>
<br>
<br>
<o:p></o:p></p>
<pre>-- <o:p></o:p></pre><pre>Dr. Peter Baumann<o:p></o:p></pre><pre> - Pro=
fessor of Computer Science, Jacobs University Bremen<o:p></o:p></pre><pre>&=
nbsp; <a
href=3D"http://www.faculty.jacobs-university.de/pbaumann">www.faculty.jacob=
s-university.de/pbaumann</a><o:p></o:p></pre><pre> mail: <a
href=3D"mailto:p.baumann at jacobs-university.de">p.baumann at jacobs-university.=
de</a><o:p></o:p></pre><pre> tel: +49-421-200-3178, fax: +49-42=
1-200-493178<o:p></o:p></pre><pre> - Executive Director, rasdaman GmbH Brem=
en (HRB 147737)<o:p></o:p></pre><pre> <a
href=3D"http://www.rasdaman.com">www.rasdaman.com</a>, mail: <a
href=3D"mailto:baumann at rasdaman.com">baumann at rasdaman.com</a><o:p></o:p></p=
re><pre> tel: 0800-rasdaman, fax: 0800-rasdafax, mobile: +49-17=
3-5837882<o:p></o:p></pre><pre>"Si forte in alienas manus oberraverit =
hec peregrina epistola incertis ventis dimissa, sed Deo commendata, precamu=
r 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=3DMsoNormal><br>
<br>
<o:p></o:p></p>
<pre>-- <o:p></o:p></pre><pre>Dr. Peter Baumann<o:p></o:p></pre><pre> - Pro=
fessor of Computer Science, Jacobs University Bremen<o:p></o:p></pre><pre>&=
nbsp; <a
href=3D"http://www.faculty.jacobs-university.de/pbaumann">www.faculty.jacob=
s-university.de/pbaumann</a><o:p></o:p></pre><pre> mail: <a
href=3D"mailto:p.baumann at jacobs-university.de">p.baumann at jacobs-university.=
de</a><o:p></o:p></pre><pre> tel: +49-421-200-3178, fax: +49-42=
1-200-493178<o:p></o:p></pre><pre> - Executive Director, rasdaman GmbH Brem=
en (HRB 147737)<o:p></o:p></pre><pre> <a
href=3D"http://www.rasdaman.com">www.rasdaman.com</a>, mail: <a
href=3D"mailto:baumann at rasdaman.com">baumann at rasdaman.com</a><o:p></o:p></p=
re><pre> tel: 0800-rasdaman, fax: 0800-rasdafax, mobile: +49-17=
3-5837882<o:p></o:p></pre><pre>"Si forte in alienas manus oberraverit =
hec peregrina epistola incertis ventis dimissa, sed Deo commendata, precamu=
r 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>
--_000_87A96661E65C5541AB4D20721C2DD7F880A3C48D08EXCHMBXAulava_--