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 : [email protected] 
[[email protected]] de la part de 
[email protected] [[email protected]]
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
[email protected]<mailto:[email protected]>
http://www.inspq.qc.ca<http://www.inspq.qc.ca/>



_______________________________________________
postgis-users mailing list
[email protected]
http://postgis.refractions.net/mailman/listinfo/postgis-users

Reply via email to