[postgis-devel] [WKT Raster] Regular blocking and null extent
Pierre Racine
Pierre.Racine at sbf.ulaval.ca
Mon May 11 06:19:19 PDT 2009
>>>>> Shall we allow generation of regular blocking tables
>>>> >from non-georeferenced datasets reported by GDAL as follows:
>>>>> $ gdalinfo chessboard-64x64.tif
>>>>> Driver: GTiff/GeoTIFF
>>>>> Files: chessboard-64x64.tif
>>>>> Size is 64, 64
>>>>> Coordinate System is `'
>>>>> Metadata:
>>>>> TIFFTAG_DOCUMENTNAME=/home/mloskot/data/test/sample/chessboard-64x64.tif
>>>>> TIFFTAG_IMAGEDESCRIPTION=Sample 4-fields greyscale chessboard
>>>>> TIFFTAG_XRESOLUTION=72
>>>>> TIFFTAG_YRESOLUTION=72
>>>>> TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
>>>>> Image Structure Metadata:
>>>>> INTERLEAVE=BAND
>>>>> Corner Coordinates:
>>>>> Upper Left ( 0.0, 0.0)
>>>>> Lower Left ( 0.0, 64.0)
>>>>> Upper Right ( 64.0, 0.0)
>>>>> Lower Right ( 64.0, 64.0)
>>>>> Center ( 32.0, 32.0)
>>>>> Band 1 Block=64x64 Type=Byte, ColorInterp=Gray
>>>>>
>>>>> and assume pixelsize_x = 1 and pixelsize_y = -1 ?
>>>> That's what I would do also. There is an extent but no (geographical)
>>>> coordinate system. That's what ArcGIS does when importing an ungeoreferenced
>>>> raster.
>>> Givne what GDAL reports for chessboard-64x64.tif, could you explain
>>> what values you put into RASTER header in each row (tile) and
>>> what goes to RASTER_COLUMNS, etc.?
>>
>> In the raster header (there should be only one tile, right?):
>>
>> scaleX = 1
>> scaleY = -1
>> ipX = 0.0
>> ipY = 0.0
>> skewX = 0.0
>> skewY = 0.0
>> srid = -1 (or 0?)
>> width = 64
>> height = 64
>>
>> In RASTER_COLUMNS:
>>
>> srid = -1 (or 0?)
>> pixel_types = 8BUI
>> regular_blocking = true
>> pixelsize_x = 1.0
>> pixelsize_y = -1.0
>> blocksize_x = 64
>> blocksize_y = 64
>>
>> No?
>
>What's in the extent field of RASTER_COLUMNS?
It would be: ST_GeomFromEWKT('SRID=-1;POLYGON(0.0 0.0, 0.0 64.0, 64.0 64.0, 64.0 0.0, 0.0 0.0)')
Pierre
More information about the postgis-devel
mailing list