[Gdal-dev] Envisat Meris GCPs bad georeferencement...

Tue May 30 09:59:26 EDT 2006

Hi,

I’d like to display Meris products in a Mapserver map but i can’t get the
image to be well projected and scaled.

The image appears on the map but its extents are far out bounds and its
lower left corner is always docked to the map center.

The destination map is in projection epsg:4326

Here’s the info on an image I’m using:

----------------------------------------------------------------------------
-------------------------------------------------------------------

Driver: GTiff/GeoTIFF

Size is 2241, 2241

Coordinate System is `'

GCP Projection = GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS
84",6378137,298.

2572235630016,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Gr
eenw

ich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]]

GCP[  0]: Id=1, Info=

          (0,0) -> (1.39256,58.1008,0)

GCP[  1]: Id=2, Info=

          (560.25,0) -> (3.76845,57.7639,0)

GCP[  2]: Id=3, Info=

          (1120.5,0) -> (6.09777,57.3834,0)

GCP[  3]: Id=4, Info=

          (1680.75,0) -> (8.37677,56.9608,0)

GCP[  4]: Id=5, Info=

          (2241,0) -> (10.6024,56.4975,0)

GCP[  5]: Id=6, Info=

          (0,448.2) -> (0.869172,56.9542,0)

GCP[  6]: Id=7, Info=

          (560.25,448.2) -> (3.17538,56.6237,0)

GCP[  7]: Id=8, Info=

          (1120.5,448.2) -> (5.43918,56.2513,0)

GCP[  8]: Id=9, Info=

          (1680.75,448.2) -> (7.65714,55.8383,0)

GCP[  9]: Id=10, Info=

          (2241,448.2) -> (9.82645,55.3862,0)

GCP[ 10]: Id=11, Info=

          (0,896.4) -> (0.370076,55.8061,0)

GCP[ 11]: Id=12, Info=

          (560.25,896.4) -> (2.61121,55.4814,0)

GCP[ 12]: Id=13, Info=

          (1120.5,896.4) -> (4.81357,55.1164,0)

GCP[ 13]: Id=14, Info=

          (1680.75,896.4) -> (6.97408,54.7123,0)

GCP[ 14]: Id=15, Info=

          (2241,896.4) -> (9.09013,54.2705,0)

GCP[ 15]: Id=16, Info=

          (0,1344.6) -> (-0.107114,54.6566,0)

GCP[ 16]: Id=17, Info=

          (560.25,1344.6) -> (2.07314,54.3372,0)

GCP[ 17]: Id=18, Info=

          (1120.5,1344.6) -> (4.21788,53.979,0)

GCP[ 18]: Id=19, Info=

          (1680.75,1344.6) -> (6.32427,53.5832,0)

GCP[ 19]: Id=20, Info=

          (2241,1344.6) -> (8.38995,53.1511,0)

GCP[ 20]: Id=21, Info=

          (0,1792.8) -> (-0.564501,53.5058,0)

GCP[ 21]: Id=22, Info=

          (560.25,1792.8) -> (1.55875,53.1912,0)

GCP[ 22]: Id=23, Info=

          (1120.5,1792.8) -> (3.64935,52.8393,0)

GCP[ 23]: Id=24, Info=

          (1680.75,1792.8) -> (5.70476,52.4512,0)

GCP[ 24]: Id=25, Info=

          (2241,1792.8) -> (7.72279,52.028,0)

GCP[ 25]: Id=26, Info=

          (0,2241) -> (-1.00394,52.3539,0)

GCP[ 26]: Id=27, Info=

          (560.25,2241) -> (1.06585,52.0436,0)

GCP[ 27]: Id=28, Info=

          (1120.5,2241) -> (3.10558,51.6975,0)

GCP[ 28]: Id=29, Info=

          (1680.75,2241) -> (5.11291,51.3165,0)

GCP[ 29]: Id=30, Info=

          (2241,2241) -> (7.08585,50.9018,0)

Metadata:

  AREA_OR_POINT=Area

  TIFFTAG_XRESOLUTION=1

  TIFFTAG_YRESOLUTION=1

  TIFFTAG_RESOLUTIONUNIT=1 (unitless)

Corner Coordinates:

Upper Left  (    0.0,    0.0)

Lower Left  (    0.0, 2241.0)

Upper Right ( 2241.0,    0.0)

Lower Right ( 2241.0, 2241.0)

Center      ( 1120.5, 1120.5)

Band 1 Block=2241x8 Type=Byte, ColorInterp=Red

Band 2 Block=2241x8 Type=Byte, ColorInterp=Green

Band 3 Block=2241x8 Type=Byte, ColorInterp=Blue

----------------------------------------------------------------------------
-------------------------------------------------------------------

I just wonder if there’s an option I’m missing in Mapserver or if the image
info is simply flawed

Thank for help !

Stéphane COLZY

stephane.colzy at magellium.fr
05 61 28 81 22

MAGELLIUM
10 avenue de l'Europe
31525 - Ramonville Saint Agne

http://www.magellium.fr <http://www.magellium.fr/> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20060530/ab660823/attachment.html