Thanks very much Even. Good to know where the error is.<div><br clear="all"><span style="font-size:large"><i><font face="garamond, serif"><b><font color="#FF6600">Hilmy</font></b></font></i></span><br>
<br><br><div class="gmail_quote">On Sat, Apr 9, 2011 at 4:50 PM, Even Rouault <span dir="ltr"><<a href="mailto:even.rouault@mines-paris.org">even.rouault@mines-paris.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Selon Hilmy Hashim <<a href="mailto:hilmyh@gmail.com">hilmyh@gmail.com</a>>:<br>
<br>
I can confirm the bug. It is in fact a bug in libgeotiff. I have reported it as<br>
<a href="http://trac.osgeo.org/geotiff/ticket/38" target="_blank">http://trac.osgeo.org/geotiff/ticket/38</a> with an attached patch to fix it.<br>
<div><div></div><div class="h5"><br>
> Hi list,<br>
><br>
> I assigned geotiff raster an srs using gdal_translate with the EPSG number<br>
> as follows:<br>
><br>
> gdal_translate -a_srs EPSG:3168 sourcefile destinationfile<br>
><br>
> Using gdalinfo, I notice that the false_easting parameter is 0 where it<br>
> should be 804670.24<br>
><br>
> Coordinate System is:<br>
> PROJCS["Kertau (RSO) / RSO Malaya (m)",<br>
> GEOGCS["Kertau (RSO)",<br>
> DATUM["Kertau_RSO",<br>
> SPHEROID["Everest 1830 (RSO<br>
> 1969)",6377295.664,300.8017000000015,<br>
> AUTHORITY["EPSG","7056"]],<br>
> AUTHORITY["EPSG","6751"]],<br>
> PRIMEM["Greenwich",0],<br>
> UNIT["degree",0.0174532925199433],<br>
> AUTHORITY["EPSG","4751"]],<br>
> PROJECTION["Hotine_Oblique_Mercator"],<br>
> PARAMETER["latitude_of_center",4],<br>
> PARAMETER["longitude_of_center",102.25],<br>
> PARAMETER["azimuth",323.0257905],<br>
> PARAMETER["rectified_grid_angle",323.1301023611111],<br>
> PARAMETER["scale_factor",0.99984],<br>
> PARAMETER["false_easting",0],<br>
> PARAMETER["false_northing",0],<br>
> UNIT["metre",1,<br>
> AUTHORITY["EPSG","9001"]],<br>
> AUTHORITY["EPSG","3168"]]<br>
><br>
> I am using OSGeo4W on Windows 7 64 bit (Gdal 1.4.5.0) shell. It is the same<br>
> with gdal 1.6, 1.7, gdal-dev, gdal 1.7.2 on Ubuntu 10.04 and QGis 1.7 (gdal<br>
> 1.8).<br>
><br>
> I checked the OSGeo4W\share\proj\epsg file, but the values there are<br>
> correct. The anomaly also affects EPSG:3375 which uses the same projection.<br>
> QGis or Geoserver doesn't have a problem displaying the raster with the<br>
> assigned srs, so it looks like its just a problem with gdalinfo.<br>
><br>
> regards<br>
><br>
</div></div>> *Hilmy*<br>
><br>
<br>
<br>
</blockquote></div><br></div>