[postgis-users] RE : importing Raster with wkt raster

[Pierre Racine]

Hi Steve,

Refer to http://trac.osgeo.org/postgis/wiki/WKTRaster/SpecificationWorking01 for the most complete documentation about gdal2wktraster.py The official doc is still in development. Look in the PostgreSQL doc for a complete understanding of the -m parameter.

You have probably produced a .sql file. Now you have to load and execute it in PostgreSQL the same way you did with rtpostgis.sql. You should then have a table containing all your tiles...

Pierre


________________________________________
De : postgis-users-bounces at postgis.refractions.net [postgis-users-bounces at postgis.refractions.net] de la part de Steve.Toutant at inspq.qc.ca [Steve.Toutant at inspq.qc.ca]
Date d'envoi : 18 janvier 2010 16:32
À : PostGIS Users Discussion
Objet : [postgis-users] importing Raster with wkt raster

Hi,
I tried importing a tif using wktraster. I used this as an example
http://mateusz.loskot.net/2009/03/30/crunching-wkt-rasters/

1-
I don't understand the use of -m parameters in the example. this option doesn't exist in gdal2wktraster.py -h

2-
I need more information on how to use properly the -k option. Blocksize is it the Size of the image or should I enter the pixel size

Here is the syntax I used and below the gdalinfo of my .TIF
Can you please tell if I'm importing the raster properly?


gdal2wktraster.py -r estMtl.tif -t estMtl_tif_32198 -s 32198 -o wkt.sql -k 30x30

Here is the gdalinfo of the .tif
Driver: GTiff/GeoTIFF
Files: EstMtl.tif
Size is 7771, 7410
Coordinate System is:
PROJCS["NAD_1983_Lambert_Conformal_Conic",
    GEOGCS["NAD83",
        DATUM["North_American_Datum_1983",
            SPHEROID["GRS 1980",6378137,298.2572221010002,
                AUTHORITY["EPSG","7019"]],
            AUTHORITY["EPSG","6269"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4269"]],
    PROJECTION["Lambert_Conformal_Conic_2SP"],
    PARAMETER["standard_parallel_1",46],
    PARAMETER["standard_parallel_2",60],
    PARAMETER["latitude_of_origin",44],
    PARAMETER["central_meridian",-68.5],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (-464859.388848508130000,348513.395163321870000)
Pixel Size = (30.000000000000000,-30.000000000000000)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  ( -464859.389,  348513.395) ( 74d37'31.17"W, 46d57'11.19"N)
Lower Left  ( -464859.389,  126213.395) ( 74d23'5.99"W, 44d57'37.51"N)
Upper Right ( -231729.389,  348513.395) ( 71d33'32.58"W, 47d 5'16.37"N)
Lower Right ( -231729.389,  126213.395) ( 71d26'19.00"W, 45d 5'21.92"N)
Center      ( -348294.389,  237363.395) ( 73d 0'7.85"W, 46d 1'58.13"N)
Band 1 Block=256x256 Type=Byte, ColorInterp=Gray
  NoData Value=-3.4028234663852886e+038

thanks
Steve

Steve Toutant, M. Sc.
Analyste en géomatique
Secteur environnement
Direction de la santé environnementale et de la toxicologie
Institut national de santé publique du Québec
945, avenue Wolfe
Québec, Qc G1V 5B3

Tél.: (418) 650-5115 #5281
Fax.: (418) 654-3144
steve.toutant at inspq.qc.ca<mailto:steve.toutant at inspq.qc.ca>
http://www.inspq.qc.ca<http://www.inspq.qc.ca/>