[postgis-users] PostGis Raster and GDAL 1.9

Tue May 22 10:01:03 PDT 2012

Hi,

I did tests and all gave the result expected. It seem to be a problem
with the origin of the result file (see attached file). The application
has taken the origin of the file 056n05_0100_deme.tif instead of the
origin of the file 056n05_0100_demw.tif. What do you think ?

For your information the files are in geographic coordinate system.

David

> From: Pierre.Racine at sbf.ulaval.ca
> To: postgis-users at postgis.refractions.net
> Date: Tue, 22 May 2012 11:37:03 -0400
> Subject: Re: [postgis-users] PostGis Raster and GDAL 1.9
> 
> I would maybe do some tests:
> 
> 1) Remove the -a_nodata "-32767" parameter. Why did you have to set it? Is the nodata set properly in PostGIS for all rows?
> 
> 2) Extract one file at a time with mode=2 and then with mode=1
> 
> Pierre
> 
> > -----Original Message-----
> > From: postgis-users-bounces at postgis.refractions.net [mailto:postgis-users-
> > bounces at postgis.refractions.net] On Behalf Of David B?langer
> > Sent: Tuesday, May 22, 2012 10:25 AM
> > To: postgis-users at postgis.refractions.net
> > Subject: Re: [postgis-users] PostGis Raster and GDAL 1.9
> > 
> > Hi Pierre,
> > Yes I have a unique rid (Integer) with sequential number. I used the tool
> > raster2pgsql.exe to upload the data in the database. I uploaded 2 files
> > 
> > 056n05_0100_deme.tif
> > 056n05_0100_demw.tif
> > 
> > "C:\Program Files\PostgreSQL\9.1\bin\raster2pgsql.exe" -I -e -Y -F -s 4269
> > e:\dir\*.tif dbelange.elevationTest | "psql.exe" -U myuser -d myDB -h myhost -p
> > 14070
> > 
> > I used this command to extract data:
> > "C:\Program Files\GDAL\gdal_translate" -a_nodata "-32767" -of GTIFF
> > "PG:host=myhost port=14070 dbname='myDB' user='myuser'
> > password='mypassword' schema='dbelange' table='elevationTest' mode=2
> > where=\"filename LIKE \'%056n05%\'\" " 056N05_testGDAL1.9.tif
> > 
> > Like I said it’s working very well with gdal 1.8
> > 
> > I can send you the 2 files If you want to try on your side.
> > 
> > Do you have another idea?
> > 
> > David
> > 
> > 
> > > From: Pierre.Racine at sbf.ulaval.ca
> > > To: postgis-users at postgis.refractions.net
> > > Date: Sat, 19 May 2012 13:55:47 -0400
> > > Subject: Re: [postgis-users] PostGis Raster and GDAL 1.9
> > >
> > > Hi David,
> > >
> > > Do you have a "rid" column in your table with unique values?
> > >
> > > Pierre
> > >
> > > > -----Original Message-----
> > > > From: postgis-users-bounces at postgis.refractions.net [mailto:postgis-users-
> > > > bounces at postgis.refractions.net] On Behalf Of David B?langer
> > > > Sent: Thursday, May 17, 2012 8:13 PM
> > > > To: postgis-users at postgis.refractions.net
> > > > Subject: [postgis-users] PostGis Raster and GDAL 1.9
> > > >
> > > > Hi,
> > > > I stored DEM raster data in PostGIS 2.0 and I used gdal_translate (with
> > mode=2)
> > > > to extract a mosaic of tiles. When I’m using GDAL 1.8 is working but when
> > I’m
> > > > using GDAL 1.9 I can see the values of pixel for just one tile. For example I
> > > > uploaded in the database a tile with an attribute filename=056n05_w and an
> > > > other tile with an filename=056n05_e and I executed the command with
> > GDAL
> > > > 1.8 and GDAL 1.9:
> > > >
> > > > gdal_translate" -a_nodata "-32767" -of GTIFF "PG:host=myhost port=14070
> > > > dbname='mybd' user='myuser' password='mypassword'
> > > > schema='myschema'table='Elevation' mode=2 where=\"filename LIKE
> > > > \'%056n05%\'\" " 056N05.tif
> > > >
> > > > With GDAL 1.8 I obtained a mosaic of 2 tiles (the result expected) but with
> > GDAL
> > > > 1.9, I obtained a raster with the same extent but only the values of the tile
> > > > 056n05_w are present. The pixels of the tile 056n05_w have a value 0.
> > > >
> > > > I tried the stable version 1.9 MSVC2008 (Win32) and development version
> > > > available today at http://www.gisinternals.com/sdk/
> > > >
> > <https://email.nrcan.gc.ca/owa/redir.aspx?C=OzXx4wlk90m3TV7W8fzG7_FxObj
> > > >
> > PB88I1bAwXOtZColSmDyhIvj_1yhbCVJUJ0vrnCWWTuivB9Y.&URL=http://www.g
> > > > isinternals.com/sdk/> and I have still the same problem. For your
> > information,
> > > > I'm running GDAL on windows7 and PostGIS 2.0 is installed on a Linux server.
> > I
> > > > tried many version of GDAL 1.9 and I always have the same problem. The
> > version
> > > > 1.8 is working very well but I would like to take advantage of the
> > improvements
> > > > in the last version because I'm working on a project where processing time is
> > > > very important.
> > > >
> > > > Does anybody have any idea about this problem ?
> > > >
> > > > Thanks
> > > >
> > > > David Bélanger
> > >
> > > _______________________________________________
> > > postgis-users mailing list
> > > postgis-users at postgis.refractions.net
> > > http://postgis.refractions.net/mailman/listinfo/postgis-users
> 
> _______________________________________________
> postgis-users mailing list
> postgis-users at postgis.refractions.net
> http://postgis.refractions.net/mailman/listinfo/postgis-users

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/postgis-users/attachments/20120522/43aef4ed/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: GDAL19_Mosaicing_problem.doc
Type: application/msword
Size: 172032 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/postgis-users/attachments/20120522/43aef4ed/attachment.doc>