<div>Sorry to bring this up again. I've been looking for why epsg_tr.py -proj4 is not including the +gamma (RectifiedToSkew angle) parameter for the +proj=omerc projection as used in EPSG:3168, EPSG:3375, etc. It seems that it is not set inĀ <a href="https://svn.osgeo.org/fdo/trunk/Thirdparty/gdal/ogr/ogr_srs_proj4.cpp">https://svn.osgeo.org/fdo/trunk/Thirdparty/gdal/ogr/ogr_srs_proj4.cpp</a>:</div>
<div><span class="Apple-style-span" style="font-family: monospace; white-space: pre-wrap; font-size: medium; "><br></span></div><div><span class="Apple-style-span" style="font-family: monospace; white-space: pre-wrap; font-size: medium; ">else if( EQUAL(pszProj,"omerc") )</span></div>
<div><span class="Apple-style-span" style="font-family: 'Times New Roman'; font-size: medium; "><pre style="word-wrap: break-word; white-space: pre-wrap; "> {
SetHOM( OSR_GDV( papszNV, "lat_0", 0.0 ),
OSR_GDV( papszNV, "lonc", 0.0 ),
OSR_GDV( papszNV, "alpha", 0.0 ),
0.0, /* ??? */
OSR_GDV( papszNV, "k", 1.0 ),
OSR_GDV( papszNV, "x_0", 0.0 ),
OSR_GDV( papszNV, "y_0", 0.0 ) );
}</pre><pre style="word-wrap: break-word; white-space: pre-wrap; ">...</pre></span><span class="Apple-style-span" style="font-family: 'Times New Roman'; font-size: medium; "><pre style="word-wrap: break-word; white-space: pre-wrap; ">
else
{
sprintf( szProj4+strlen(szProj4),
"+proj=omerc +lat_0=%.16g +lonc=%.16g +alpha=%.16g"
" +k=%.16g +x_0=%.16g +y_0=%.16g ",
GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN,0.0),
GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN,0.0),
GetNormProjParm(SRS_PP_AZIMUTH,0.0),
GetNormProjParm(SRS_PP_SCALE_FACTOR,1.0),
GetNormProjParm(SRS_PP_FALSE_EASTING,0.0),
GetNormProjParm(SRS_PP_FALSE_NORTHING,0.0) );
}</pre></span></div>There seems to be some uncertainty about this parameter. How can this be resolved and how can I help?<div><br></div><div>Thanks.</div><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>
</div>