Трассировка лучей, OpenGl, C++ Опции

[legat]

Необходимо реализовать трассировку лучей на OpenGl и C++, задача вроде проста, есть базовая сфера( радиус =1, центр в начале координат) необходимо ее отрисовать. У меня почти получилось, только вот точки пересечения не получается найти, по моему что-то с заданием направления луча. Может кто-то сталкивался с этим? Буду благодарен за код=)

void Reshape(int width, int height)
{if (width>height)
	glViewport((width-height)/2,0,height,height);
	else glViewport(0,(height-width)/2,width,width);
 //glEnable(GL_DEPTH_TEST);
 glMatrixMode(GL_PROJECTION);
 glLoadIdentity();
 gluOrtho2D(-10,10,-10,10);
 glMatrixMode(GL_MODELVIEW);
}
void Display(void)
{
	glClear(GL_COLOR_BUFFER_BIT);
	Ray myray;	
	myray.start.x=0;
	myray.start.y=0;
	myray.start.z=5;
	for(float i=0;i<200;i++)
		for(float j=0;j<200;j++)
		{
			Point t;
			t.x=-10+i/20;
			t.y=10-j/20;
			t.z=0;
			//FILE * f=fopen("output.txt","r");
			myray.dir.x=t.x-myray.start.x; //задаю направление луча, думаю ошибка где-то здесь, хотя //может неправильно перевожу экранные координаты
			myray.dir.y=t.y-myray.start.y;
			myray.dir.z=t.z-myray.start.z;
			//fprintf(f,"x=%f, y=%f, z=%f\n",myray.dir.x,myray.dir.y,myray.dir.z);
			float A=sqr(myray.dir,myray.dir);
			float B=sqr(myray.dir,myray.start);
			float C=sqr(myray.start,myray.start)-1;
			float D=B*B-4*A*C;
			if (D>0)
			printf("i=%d j=%d %f\n",i,j,D);
			if (D<0)
			{
				glColor3f(0,0,0);
				glRecti(i,j,i+1,j+1);						
			}
			else
			{
				//float t1,t2;
				//t1=(-B-sqrt(D))/(2*A);
				//t1=(-B+sqrt(D))/(2*A);
				glColor3f(1,1,1);
				glRecti(i,j,i+1,j+1);				
			}
			//glFlush();
		}	
 glFlush();
}

Может кто знает, где можно найти книгу Лорен Хейни "Построение изображений метод слежения луча"?

Сообщение отредактировано: Lapp - 3.03.2010 5:22

[legat]

Пишу трассировку лучей по книге "Программирование компьютерной графики" Ф. Хилла
направление луча задается след. образом:

Код

float x=-W+(float)col*2*W/(float)nCols;
           float y=-H+(float)row*2*H/(float)nRows;
    //вычисляем направление луча
    theRay.setDir(-nearDist*n.x+x*u.x+y*v.x,
            -nearDist*n.y+x*u.y+y*v.y,
            -nearDist*n.z+x*u.z+y*v.z);

вектора n, u, v - поля класса Camera
а камеру задаю:

Код

Point3 eye(0,0,100);
Point3 look(0,0,0);
Vector3 up(0,1,0);
Camera cam;
cam.set(eye,look,up);

может кто задавал параметры для камеры? а то с моими получается, что луч никогда не пересекает даже хотя бы базовую сферу

 Raytracer.v.2.rar ( 3.69 килобайт ) Кол-во скачиваний: 590
(а вот весь код, который у меня есть)