//
//  RMProjection.m
//
// Copyright (c) 2008-2009, Route-Me Contributors
// All rights reserved.
//
// Redistribution and use in source and binary forms, with or without
// modification, are permitted provided that the following conditions are met:
//
// * Redistributions of source code must retain the above copyright notice, this
//   list of conditions and the following disclaimer.
// * Redistributions in binary form must reproduce the above copyright notice,
//   this list of conditions and the following disclaimer in the documentation
//   and/or other materials provided with the distribution.
//
// THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
// AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
// IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
// ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
// LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
// CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
// SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
// INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
// CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
// ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
// POSSIBILITY OF SUCH DAMAGE.

#import "RMGlobalConstants.h"
#import "proj_api.h"
#import "RMProjection.h"

@implementation RMProjection

@synthesize internalProjection;
@synthesize planetBounds;
@synthesize projectionWrapsHorizontally;

- (id)initWithString:(NSString *)params inBounds:(RMProjectedRect)projectedBounds
{
	if (!(self = [super init]))
		return nil;

	internalProjection = pj_init_plus([params UTF8String]);
	if (internalProjection == NULL)
	{
		RMLog(@"Unhandled error creating projection. String is %@", params);
		[self release];
		return nil;
	}

	planetBounds = projectedBounds;
	projectionWrapsHorizontally = YES;

	return self;
}

- (id)initWithString:(NSString *)params
{
	RMProjectedRect theBounds;
	theBounds = RMMakeProjectedRect(0, 0, 0, 0);
	return [self initWithString:params inBounds:theBounds];
}

- (id)init
{
	return [self initWithString:@"+proj=latlong +ellps=WGS84"];
}

- (void)dealloc
{
	if (internalProjection)
		pj_free(internalProjection);

	[super dealloc];
}

- (RMProjectedPoint)wrapPointHorizontally:(RMProjectedPoint)aPoint
{
	if (!projectionWrapsHorizontally || planetBounds.size.width == 0.0f || planetBounds.size.height == 0.0f)
		return aPoint;

	while (aPoint.easting < planetBounds.origin.easting)
		aPoint.easting += planetBounds.size.width;

	while (aPoint.easting > (planetBounds.origin.easting + planetBounds.size.width))
		aPoint.easting -= planetBounds.size.width;

	return aPoint;
}

- (RMProjectedPoint)constrainPointToBounds:(RMProjectedPoint)aPoint
{
	if (planetBounds.size.width == 0.0f || planetBounds.size.height == 0.0f)
		return aPoint;

	[self wrapPointHorizontally:aPoint];

	if (aPoint.northing < planetBounds.origin.northing)
		aPoint.northing = planetBounds.origin.northing;
	else if (aPoint.northing > (planetBounds.origin.northing + planetBounds.size.height))
		aPoint.northing = planetBounds.origin.northing + planetBounds.size.height;

	return aPoint;
}

- (RMProjectedPoint)coordinateToProjectedPoint:(CLLocationCoordinate2D)aLatLong
{
	projUV uv = {
		aLatLong.longitude * DEG_TO_RAD,
		aLatLong.latitude * DEG_TO_RAD
	};

	projUV result = pj_fwd(uv, internalProjection);

	RMProjectedPoint result_point = {
		result.u,
		result.v,
	};

	return result_point;
}

- (CLLocationCoordinate2D)projectedPointToCoordinate:(RMProjectedPoint)aPoint
{
	projUV uv = {
		aPoint.easting,
		aPoint.northing,
	};

	projUV result = pj_inv(uv, internalProjection);

	CLLocationCoordinate2D result_coordinate = {
		result.v * RAD_TO_DEG,
		result.u * RAD_TO_DEG,
	};

	return result_coordinate;
}

static RMProjection *_google = nil;
static RMProjection *_latlong = nil;
static RMProjection *_osgb = nil;

+ (RMProjection *)googleProjection
{
	if (_google) {
		return _google;
	}
	else {
		RMProjectedRect theBounds = RMMakeProjectedRect(-20037508.34, -20037508.34, 20037508.34 * 2, 20037508.34 * 2);

		_google = [[RMProjection alloc] initWithString:@"+title= Google Mercator EPSG:900913 +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs"
											  inBounds: theBounds];
		return _google;
	}
}

+ (RMProjection *)EPSGLatLong
{
	if (_latlong) {
		return _latlong;
	}
	else {
		RMProjectedRect theBounds = RMMakeProjectedRect(-kMaxLong, -kMaxLat, 360, kMaxLong);

		_latlong = [[RMProjection alloc] initWithString:@"+proj=latlong +ellps=WGS84" inBounds:theBounds];
		return _latlong;
	}
}

@end