作者:朱金灿
来源:http://blog.csdn/clever101
用记事本打开gdal源码ogr目录下的ogr_srs_proj4.cpp,在OGRErr OGRSpatialReference::exportToProj4( char ** ppszProj4 ) const
函数添加下面一段代码:
else if( EQUAL(pszProjection,"Hammer_Aitoff") )
{
CPLsprintf( szProj4+strlen(szProj4),
"+proj=hammer +lon_0=%.16g +x_0=%.16g +y_0=%.16g ",
GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN,0.0),
GetNormProjParm(SRS_PP_FALSE_EASTING,0.0),
GetNormProjParm(SRS_PP_FALSE_NORTHING,0.0) );
}
然后重新编译gdal源码即可。
更多推荐
gdal支持Hammer投影的做法
发布评论