gdal支持Hammer投影的做法

版权声明：本文为博主原创文章，未经博主允许不得转载。 https://blog.csdn.net/clever101/article/details/83510698

作者：朱金灿

来源：http://blog.csdn.net/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源码即可。