#include "ellipsoidTechnique.h"


ellipsoidTechnique::ellipsoidTechnique()
{
	_height = 0.0;
}

ellipsoidTechnique::~ellipsoidTechnique()
{
}

void ellipsoidTechnique::modifyHeightfield(region& modificationROI, osg::HeightField* h, region tileExtends)
{
	//OSG_NOTIFY( osg::ALWAYS ) << "ellipsoidTechnique::modifyHeightfield()" << std::endl;
	//OSG_NOTIFY( osg::ALWAYS ) << "LAT: " << tileExtends._lat_min << " | " << tileExtends._lat_max << std::endl;
	//OSG_NOTIFY( osg::ALWAYS ) << "LON: " << tileExtends._lon_min << " | " << tileExtends._lon_max << std::endl;

	// calculate colum start/end and row start/end of affected vertices
	int startX=0, startY=0, endX=h->getNumColumns(), endY=h->getNumRows();

	
	// Determine loop variables..
	if(tileExtends._lat_min<modificationROI._lat_min)
		startY = (modificationROI._lat_min - tileExtends._lat_min) / h->getYInterval();
	if(tileExtends._lat_max>modificationROI._lat_max)
		endY = (modificationROI._lat_max - tileExtends._lat_min) / h->getYInterval();

	if(tileExtends._lon_min<modificationROI._lon_min)
		startX = (modificationROI._lon_min - tileExtends._lon_min) / h->getXInterval();
	if(tileExtends._lon_max>modificationROI._lon_max)
		endX = (modificationROI._lon_max - tileExtends._lon_min) / h->getXInterval();


	// Modify height value of affected vertices
	for(int x=startX;x<endX;x++)
	{
		for(int y=startY;y<endY;y++)
		{
			h->setHeight( x, y, _height);
		}
	}
}