int msProjectShape(projectionObj *in, projectionObj *out, shapeObj *shape)
{
#ifdef USE_PROJ
int i;
for( i = shape->numlines-1; i >= 0; i-- )
{
if( shape->type == MS_SHAPE_LINE || shape->type == MS_SHAPE_POLYGON )
{
if( msProjectShapeLine( in, out, shape, i ) == MS_FAILURE )
msShapeDeleteLine( shape, i );
}
else if( msProjectLine(in, out, shape->line+i ) == MS_FAILURE )
{
msShapeDeleteLine( shape, i );
}
}
if( shape->numlines == 0 ) {
msFreeShape( shape );
return MS_FAILURE;
} else {
msComputeBounds( shape ); /* fixes bug 1586 */
return(MS_SUCCESS);
}
#else
msSetError(MS_PROJERR, "Projection support is not available.", "msProjectShape()");
return(MS_FAILURE);
#endif
}
int msProjectShape(projectionObj *in, projectionObj *out, shapeObj *shape)
{
#ifdef USE_PROJ
int i;
#ifdef USE_PROJ_FASTPATHS
int j;
if(in->wellknownprojection == wkp_lonlat && out->wellknownprojection == wkp_gmerc) {
for( i = shape->numlines-1; i >= 0; i-- ) {
for( j = shape->line[i].numpoints-1; j >= 0; j-- ) {
#define p_x shape->line[i].point[j].x
#define p_y shape->line[i].point[j].y
p_x *= MAXEXTENTby180;
p_y = log(tan((90 + p_y) * M_PIby360)) * MS_RAD_TO_DEG;
p_y *= MAXEXTENTby180;
if (p_x > MAXEXTENT) p_x = MAXEXTENT;
if (p_x < -MAXEXTENT) p_x = -MAXEXTENT;
if (p_y > MAXEXTENT) p_y = MAXEXTENT;
if (p_y < -MAXEXTENT) p_y = -MAXEXTENT;
#undef p_x
#undef p_y
}
}
msComputeBounds( shape ); /* fixes bug 1586 */
return MS_SUCCESS;
}
#endif
for( i = shape->numlines-1; i >= 0; i-- ) {
if( shape->type == MS_SHAPE_LINE || shape->type == MS_SHAPE_POLYGON ) {
if( msProjectShapeLine( in, out, shape, i ) == MS_FAILURE )
msShapeDeleteLine( shape, i );
} else if( msProjectLine(in, out, shape->line+i ) == MS_FAILURE ) {
msShapeDeleteLine( shape, i );
}
}
if( shape->numlines == 0 ) {
msFreeShape( shape );
return MS_FAILURE;
} else {
msComputeBounds( shape ); /* fixes bug 1586 */
return(MS_SUCCESS);
}
#else
msSetError(MS_PROJERR, "Projection support is not available.", "msProjectShape()");
return(MS_FAILURE);
#endif
}
