Context Navigation

DistortionSetupStrategyProjectSyntropy.cpp @ 418

Last change on this file since 418 was 410, checked in by Torben Dannhauer, 12 years ago
New Distortion setup strategy added: distortionNG can now handle configuration files of "Projection Designer"
File size: 13.0 KB

Rev	Line
[372]	1	/* osgVisual test. distortionNG, experimental.
	2	*
	3	* Permission is hereby granted, free of charge, to any person obtaining a copy
	4	* of this software and associated documentation files (the "Software"), to deal
	5	* in the Software without restriction, including without limitation the rights
	6	* to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
	7	* copies of the Software, and to permit persons to whom the Software is
	8	* furnished to do so, subject to the following conditions:
	9	*
	10	* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
	11	* IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
	12	* FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
	13	* AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
	14	* LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
	15	* OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
	16	* THE SOFTWARE.
	17	*/
	18
[378]	19	#include <osg/Vec3d>
	20	#include <osg/ImageUtils>
[410]	21	#include <osgDB/FileUtils>
[376]	22	#include <osgDB/ReadFile>
	23	#include <osgDB/fstream>
	24
[372]	25	#include "DistortionSetupStrategyProjectSyntropy.h"
	26
	27	DistortionSetupStrategyProjectSyntropy::DistortionSetupStrategyProjectSyntropy()
	28	{
[377]	29	_blendmapFilename = "";
	30	_frustumFilename = "";
[394]	31	_distortionFilename = "";
[376]	32
	33	_distortionInitialized=false;
[372]	34	}
	35
	36	DistortionSetupStrategyProjectSyntropy::~DistortionSetupStrategyProjectSyntropy()
	37	{
[377]	38
[372]	39	}
	40
[394]	41	void DistortionSetupStrategyProjectSyntropy::setDistortionInputFiles(std::string distortionFile, std::string blendmapFile, std::string frustumFile)
[376]	42	{
[377]	43	_blendmapFilename = blendmapFile;
	44	_frustumFilename = frustumFile;
[394]	45	_distortionFilename = distortionFile;
[376]	46	}
	47
	48	int DistortionSetupStrategyProjectSyntropy::readFrustumFromCSVFile(std::string filePath, int numValues, float* frustumValues)
	49	{
	50	int numSeparator=0;
	51	std::string line, valueString;
	52	char separator=';';
	53	size_t valueStart, valueEnd;
	54
	55	osgDB::ifstream inputFile (const_cast<char*>(filePath.c_str()));
	56
	57	for (int i=0;i<numValues; i++) frustumValues[i]=0.0;
	58
	59	if (inputFile.is_open())
	60	{
	61	inputFile.clear();
	62	inputFile.seekg (0, std::ios::beg);
	63
	64	while (!inputFile.eof())
	65	{
	66	getline (inputFile,line);
	67
	68	if((line.size()>0) && (line.find('x')==std::string::npos))
	69	{new osg::Vec2Array;
	70	numSeparator=0;
	71	for (size_t numLineChar=0; numLineChar < line.length(); numLineChar++)
	72	{
	73	if(line.at(numLineChar)==separator) numSeparator++;
	74	}
	75	if(numSeparator!=numValues-1) //Number of Entries per Line is not correct
	76	{
	77	OSG_ALWAYS<< "Inconsistency in file " << filePath << " detected. Please check this file." <<std::endl;
	78	return 1;
	79	}
	80	valueStart=0;
	81	valueEnd=0;
	82	for (int j=0;j<numValues; j++)
	83	{
	84	if(j!=numValues-1)
	85	{
	86	if((valueEnd=line.find(separator,valueStart))!=std::string::npos)
	87	{
	88	valueString=line.substr(valueStart,valueEnd-valueStart);
	89	valueStart=valueEnd+1;
	90	}
	91	}
	92	else valueString=line.substr(valueStart,line.size());
	93
	94	frustumValues[j]=atof(const_cast<char*>(valueString.c_str()));
	95	}
	96	}
	97	}
	98	inputFile.close();
	99	}
	100	else //File couldn't be loaded
	101	{
	102	OSG_ALWAYS<< "Unable to open " << filePath << std::endl;
	103	return 1;
	104	}
	105	return 0;
	106	}
	107
	108	int DistortionSetupStrategyProjectSyntropy::readMeshDimensionsFromCSVFile(std::string filePath, int* meshRows, int* meshColumns)
	109	{
	110	int tempRows=0, tempCols=0, numSeparator=0, numValues=6;
	111	std::string line, valueString;
	112	char separator=';';
	113	size_t valueStart, valueEnd;
	114
	115	osgDB::ifstream inputFile (const_cast<char*>(filePath.c_str()));
	116
	117	if (inputFile.is_open())
	118	{
	119	inputFile.clear();
	120	inputFile.seekg (0, std::ios::beg);
	121
	122	while (!inputFile.eof())
	123	{
	124	getline (inputFile,line);
	125
	126	if((line.size()>0) && (line.find('x')==std::string::npos))
	127	{
	128	numSeparator=0;
	129	for (size_t numLineChar=0; numLineChar < line.length(); numLineChar++)
	130	{
	131	if(line.at(numLineChar)==separator) numSeparator++;
	132	}
	133	if(numSeparator!=numValues-1) //Number of Entries per Line is not correct
	134	{
	135	OSG_ALWAYS<< "Inconsistency in file " << filePath << " detected. Please check this file." <<std::endl;
	136	return 1;
	137	}
	138	//Get number of rows of the distortion grid, last valid line of input file contains actual value
	139	valueStart=line.rfind(separator,line.size());
	140	valueString=line.substr(valueStart+1,line.size()-valueStart-1);
	141	tempRows=atoi(const_cast<char*>(valueString.c_str()));
	142	if(tempRows>meshRows) meshRows=tempRows;
	143
	144	//Get number of columns of the distortion grid, last valid line of input file contains actual value
	145	valueEnd=valueStart-1;
	146	valueStart=line.rfind(separator,valueEnd);
	147	valueString=line.substr(valueStart+1,valueEnd-valueStart);
	148	tempCols=atoi(const_cast<char*>(valueString.c_str()));
	149	if(tempCols>meshColumns) meshColumns=tempCols;
	150	}
	151
	152	}
	153	inputFile.close();
	154	}
	155	else //File couldn't be loaded
	156	{
	157	OSG_ALWAYS<< "Unable to open " << filePath << std::endl;
	158	return 1;
	159	}
	160
	161	*meshColumns+=1;
	162	*meshRows+=1;
	163
	164	return 0;
	165	}
	166
[394]	167	int DistortionSetupStrategyProjectSyntropy::readMeshPointsFromCSVFile(std::string filePath, osg::Vec4Array* tmpMesh)
[376]	168	{
	169	int numValues=6;
	170	std::string line, valueString;
	171	char separator=';';
	172	float* values = new float[numValues];
	173	size_t valueStart, valueEnd;
	174	osg::Vec2 tempVec2;
	175	osg::Vec3 tempVec3;
	176
	177	osgDB::ifstream inputFile (const_cast<char*>(filePath.c_str()));
	178
	179	if (inputFile.is_open())
	180	{
	181	inputFile.clear();
	182	inputFile.seekg (0, std::ios::beg);
	183
	184	while (!inputFile.eof())
	185	{
	186	getline (inputFile,line);
	187
	188	if((line.size()>0) && (line.find('x')==std::string::npos))
	189	{
	190	for (int i=0;i<numValues; i++) values[i]=0.0;
	191
	192	valueStart=0;
	193	valueEnd=0;
	194	for (int j=0;j<numValues; j++)
	195	{
	196	if(j!=numValues-1)
	197	{
	198	if((valueEnd=line.find(separator,valueStart))!=std::string::npos)
	199	{
	200	valueString=line.substr(valueStart,valueEnd-valueStart);
	201	valueStart=valueEnd+1;
	202	}
	203	}
	204	else valueString=line.substr(valueStart,line.size());
	205
	206	values[j]=atof(const_cast<char*>(valueString.c_str()));
	207	}
	208
[382]	209	// Store values into array
[400]	210	tmpMesh->push_back(osg::Vec4(values[0], 1.0-values[1], values[2], 1.0-values[3]));
[376]	211	}
	212	}
	213	inputFile.close();
	214	}
	215	else
	216	{
	217	OSG_ALWAYS<< "Unable to open " << filePath << std::endl;
	218	return 1;
	219	}
	220
	221	return 0;
	222	}
	223
[372]	224	void DistortionSetupStrategyProjectSyntropy::delegateDistortionSetup(osgViewer::DistortionSet* distortionSet)
	225	{
	226	if(distortionSet == NULL)
	227	{
	228	OSG_ALWAYS<<"DistortionSetupStrategyProjectSyntropy::delegateDistortionSetup : Invalid DistortionSet"<<std::endl;
	229	return;
	230	}
	231
[394]	232	if( _blendmapFilename.empty() \|\| _frustumFilename.empty() \|\| _distortionFilename.empty() )
[377]	233	{
[394]	234	OSG_ALWAYS<<"DistortionSetupStrategyProjectSyntropy::delegateDistortionSetup : You have not specified the imput filenames correctly!"<<std::endl;
[377]	235	return;
	236	}
	237
[410]	238	if( !osgDB::fileExists(_blendmapFilename) \|\| !osgDB::fileExists(_frustumFilename) \|\| !osgDB::fileExists(_distortionFilename) )
	239	{
	240	OSG_ALWAYS<<"DistortionSetupStrategyProjectionDesigner::delegateDistortionSetup : One of the specified input files does not exist!"<<std::endl;
	241	return;
	242	}
	243
[376]	244	//OSG_ALWAYS<<"drin"<<std::endl;
[372]	245
[377]	246	if (_distortionInitialized==false )
[376]	247	{
[379]	248	_distortionInitialized=true;
[383]	249	distortionSet->dirtyMesh();
	250	distortionSet->dirtyMatrix();
[376]	251
[383]	252
[376]	253	// ******************** //
	254	// * IntensityMap * //
	255	// ******************** //
	256
	257	osg::Image* intensityMap=NULL;
[377]	258	intensityMap=osgDB::readImageFile(_blendmapFilename);
[376]	259	osg::GraphicsContext::WindowingSystemInterface* wsi = osg::GraphicsContext::getWindowingSystemInterface();
	260	if (!wsi)
	261	{
	262	OSG_NOTICE<<"Error, no WindowSystemInterface available, cannot create windows."<<std::endl;
	263	return;
	264	}
	265
	266	osg::GraphicsContext::ScreenIdentifier si;
	267	si.readDISPLAY();
	268
	269	if (si.displayNum<0) si.displayNum = 0;
	270
	271	si.screenNum = 0; //Fixed value, may be replaced later
	272
	273	unsigned int width=0, height=0;
	274	wsi->getScreenResolution(si, width, height);
	275
	276	if(intensityMap->s()!=width \|\| intensityMap->t()!=height)
	277	intensityMap->scaleImage(width, height, intensityMap->r());
	278
[378]	279	// Copy image data and flag as dirty to update texture.
	280	osg::copyImage(intensityMap, 0, 0, 0, intensityMap->s(), intensityMap->t(), intensityMap->r(), distortionSet->getIntensityMap(), 0, 0, 0, true);
	281	distortionSet->getIntensityMap()->dirty();
[376]	282
	283	// ******************************* //
	284	// * Frustum and View Offset * //
	285	// ******************************* //
	286
	287	//....Read Files and set parameters
	288	int numFrustumValues=16;
	289	float* frustumValues = new float[numFrustumValues];
[400]	290	double zNear=1.0, zFar=1000.0;
[376]	291
[377]	292	readFrustumFromCSVFile(_frustumFilename, numFrustumValues, frustumValues);
	293
[391]	294	osg::Matrixd viewOffset=osg::Matrix::identity();
[376]	295	osg::Matrixd projectionOffset=osg::Matrixd();
	296
	297	// View Offset Parameters: x, y, z, heading, pitch, bank
	298	// Translational Parameters are defines in millimeters (->Correction required?!)
	299	// Coordinate System / Directions may be looked up in "Koordinatensystem.png", trans/rot order confirmed bei PS
[391]	300	osg::Matrixd rotViewOffset=osg::Matrixd();
	301	rotViewOffset.makeRotate(
	302	osg::DegreesToRadians(frustumValues[3]), osg::Vec3(0,1,0), // heading
	303	osg::DegreesToRadians(-frustumValues[4]), osg::Vec3(1,0,0), // pitch
	304	osg::DegreesToRadians(frustumValues[5]), osg::Vec3(0,0,1)); // roll
[376]	305
[391]	306	osg::Matrixd transViewOffset=osg::Matrixd();
	307	transViewOffset.makeTranslate( frustumValues[0], frustumValues[1], frustumValues[2]);
	308
	309	viewOffset = viewOffset * rotViewOffset; // * transViewOffset;
	310
[399]	311	//Frustum Parameters: Left, Right, Bottom, Top, zNear, zFar.
[410]	312	// Important: By definition of OpenGL, Left, Right, Bottom and Top must be set as tangens, multiplied by zNear!
	313	projectionOffset.makeFrustum(frustumValues[10] * zNear, frustumValues[11] * zNear, frustumValues[12] * zNear, frustumValues[13] * zNear, zNear, zFar);
[376]	314
	315	distortionSet->setViewOffset(viewOffset);
	316	distortionSet->setProjectionOffset(projectionOffset);
[382]	317
[376]	318	// ******************************* //
	319	// ***** Distortion Mesh ***** //
	320	// ******************************* //
	321
	322	osg::Geometry* distortionMeshGeometry = distortionSet->getDistortionInternals()->getChild(osgViewer::DistortionSet::MESH)->asGeode()->getDrawable(0)->asGeometry();
	323	if(distortionMeshGeometry)
	324	{
	325	//*******************************************************************************************//
	326	//* PRINT VERTEXARRAY/TEXCOORDARRAY for Debug Purposes *//
	327	//* both Arrays are listed in "HorizontalLine-after-HorizontalLine" order, *//
	328	//* sorted from topleft (0,0,0 vertices; 0,0 texcoords) to bottomright *//
[382]	329	//* (1920,1080,0 vertices; 1,1 texcoords) *//
[376]	330	//*******************************************************************************************//
	331
[383]	332	#if 0
[382]	333	osg::Vec3Array* vertices = (osg::Vec3Array*) distortionMeshGeometry->getVertexArray();
	334	osg::Vec2Array* texcoords0 = (osg::Vec2Array*) distortionMeshGeometry->getTexCoordArray(0);
	335
	336	for (osg::Vec3Array::iterator it = vertices->begin(); it != vertices->end(); ++it)
[376]	337	{
[382]	338	OSG_ALWAYS<<"Vertice "<<(it).x()<<" "<<(it).y()<<" "<<(*it).z()<<std::endl;
	339	}
[376]	340
[382]	341	for (osg::Vec2Array::iterator it = texcoords0->begin(); it != texcoords0->end(); ++it)
	342	{
	343	OSG_ALWAYS<<"Texcoord "<<(it).x()<<" "<<(it).y()<<std::endl;
[376]	344	}
[382]	345	#endif
[376]	346
	347	//*******************************************************************************************//
[382]	348	//* MESH WARPING *//
[394]	349	//*** Both Vertex Warping and Texcoord Warping is supported.
	350	//* The input file contains . Although there might *//
[376]	351	//* be a preferred version for osgvisual, both options should be implemented in osg *//
	352	//*******************************************************************************************//
	353
[383]	354	// Exctract mesh dimensions from file and pass them to the container.
[394]	355
[376]	356	int distortionMeshRows=0, distortionMeshColumns=0;
[394]	357	readMeshDimensionsFromCSVFile(_distortionFilename, &distortionMeshRows, &distortionMeshColumns);
[377]	358
[376]	359	distortionSet->setDistortionMeshRows(distortionMeshRows);
	360	distortionSet->setDistortionMeshColumns(distortionMeshColumns);
	361
[383]	362	// Read in the mesh arrays
[394]	363	osg::Vec4Array* distortionMesh = new osg::Vec4Array;
	364	readMeshPointsFromCSVFile(_distortionFilename, distortionMesh);
[382]	365
[399]	366	#if 0
[387]	367	for (osg::Vec4Array::iterator it = distortionMesh->begin(); it != distortionMesh->end(); ++it)
	368	{
	369	OSG_ALWAYS<<"final vector: "<<(it).x()<<" "<<(it).y()<<" "<<(it).z()<<" "<<(it).w()<<std::endl;
	370	}
[388]	371	#endif
[387]	372
	373
[383]	374	distortionSet->setDistortionMesh(distortionMesh);
[393]	375	distortionSet->setDistortionMeshDimensions(distortionMeshRows, distortionMeshColumns);
[376]	376	}
	377	}
	378	}

Note: See TracBrowser for help on using the repository browser.

Context Navigation

source: experimental/distortionNG/DistortionSetupStrategyProjectSyntropy.cpp @ 418

Download in other formats: