/**********************************************************************
 * $Id: Depth.cpp,v 1.2.2.1 2005/05/23 18:41:51 strk Exp $
 *
 * GEOS - Geometry Engine Open Source
 * http://geos.refractions.net
 *
 * Copyright (C) 2001-2002 Vivid Solutions Inc.
 *
 * This is free software; you can redistribute and/or modify it under
 * the terms of the GNU Lesser General Public Licence as published
 * by the Free Software Foundation. 
 * See the COPYING file for more information.
 *
 **********************************************************************
 * $Log: Depth.cpp,v $
 * Revision 1.2.2.1  2005/05/23 18:41:51  strk
 * Replaced sprintf uses with ostringstream
 *
 * Revision 1.2  2004/07/02 13:28:26  strk
 * Fixed all #include lines to reflect headers layout change.
 * Added client application build tips in README.
 *
 * Revision 1.1  2004/03/19 09:48:45  ybychkov
 * "geomgraph" and "geomgraph/indexl" upgraded to JTS 1.4
 *
 * Revision 1.10  2003/11/07 01:23:42  pramsey
 * Add standard CVS headers licence notices and copyrights to all cpp and h
 * files.
 *
 *
 **********************************************************************/


#include <sstream>
#include <geos/geomgraph.h>
#define DEPTHNULL -1

namespace geos {

int Depth::depthAtLocation(int location) {
	if (location == Location::EXTERIOR) return 0;
	if (location == Location::INTERIOR) return 1;
	return DEPTHNULL;
}

Depth::Depth() {
	// initialize depth array to a sentinel value
	for (int i=0; i<2; i++) {
		for (int j=0; j<3;j++) {
			depth[i][j]=DEPTHNULL;
		}
	}
}

Depth::~Depth() {
//	delete[] &depth;
}

int Depth::getDepth(int geomIndex,int posIndex){
	return depth[geomIndex][posIndex];
}

void Depth::setDepth(int geomIndex,int posIndex,int depthValue){
	depth[geomIndex][posIndex] = depthValue;
}

int Depth::getLocation(int geomIndex,int posIndex) {
	if (depth[geomIndex][posIndex] <= 0) return Location::EXTERIOR;
	return Location::INTERIOR;
}

void Depth::add(int geomIndex,int posIndex,int location){
	if (location == Location::INTERIOR)
		depth[geomIndex][posIndex]++;
}

/**
 * A Depth object is null (has never been initialized) if all depths are null.
 */
bool Depth::isNull() {
	for (int i=0; i<2; i++) {
		for (int j=0; j<3; j++) {
			if (depth[i][j] != DEPTHNULL)
				return false;
		}
	}
	return true;
}

bool Depth::isNull(int geomIndex) {
	return depth[geomIndex][1] == DEPTHNULL;
}

bool Depth::isNull(int geomIndex, int posIndex){
	return depth[geomIndex][posIndex] == DEPTHNULL;
}

int Depth::getDelta(int geomIndex) {
	return depth[geomIndex][Position::RIGHT]-depth[geomIndex][Position::LEFT];
}

/**
 * Normalize the depths for each geometry, if they are non-null.
 * A normalized depth
 * has depth values in the set { 0, 1 }.
 * Normalizing the depths
 * involves reducing the depths by the same amount so that at least
 * one of them is 0.  If the remaining value is > 0, it is set to 1.
 */
void Depth::normalize() {
	for (int i=0; i<2; i++) {
		if (!isNull(i)) {
			int minDepth=depth[i][1];
			if (depth[i][2]<minDepth)
				minDepth=depth[i][2];
			if (minDepth<0) minDepth = 0;
			for (int j=1; j<3; j++) {
				int newValue=0;
				if (depth[i][j]>minDepth)
					newValue = 1;
				depth[i][j] = newValue;
			}
		}
	}
}

void Depth::add(Label* lbl){
	for (int i=0; i<2; i++) {
		for (int j=1; j<3; j++) {
			int loc=lbl->getLocation(i,j);
			if (loc==Location::EXTERIOR || loc==Location::INTERIOR) {
			// initialize depth if it is null, otherwise add this location value
				if (isNull(i,j)) {
					depth[i][j]=depthAtLocation(loc);
				} else
					depth[i][j]+=depthAtLocation(loc);
			}
		}
	}
}

string Depth::toString() {
	ostringstream s;
	s<<"A:"<<depth[0][1]<<","<<depth[0][2]<<" ";
	s<<"B:"<<depth[1][1]<<","<<depth[1][2]<<"]";
	return s.str();
}


}