Programming Memoirs

3D (depth) composition of CUDA ray traced images with OpenGL rasterized images using CUDA Driver API

Depth composition of CUDA ray traced image with OpenGL rasterised object transformation gizmos

Ray tracing is a great method of generating synthetic images. It has many benefits over traditionally used (e.g. in computer game) rasterization. Ray tracing stays great up to the moment when you need to render e.g. a line segment placed in your 3D space (which is potential occluded by other 3D objects).

Why would you want to ray trace a line or a line segment? Say you want to create some kind of transformation gizmo for you 3D objects which blends nicely into the scene, or a bounding box depicting the boundaries of an object, or include wire frame of meshes in your ray traced scene, or …  There are many potential uses.

You cannot just mathematically test for camera ray vs. line segment collision and expect the line segment to appear on the ray traced rendering. Chance of a camera ray colliding a line segment are too slim for the line segment to be visible.

You can try simulating a line by drawing a thin cylinder, or an ‘x’ made of 2 quads. But this is far from an elegant solution. Additionally, your ‘line’ stops being a mathematical line, as it suddenly has a width. Thus, it becomes thicker closer to the camera, and thinner further away from it (considering you use perspective projection).

Fortunately there is a solution for drawing lines in ray traced content. The solution involves using some some kind of rasterization based renderer, such as OpenGL to draw the lines (or any objects) separately from the ray tracing pass and then performing a 3D composition of the two images.  The problem becomes even more interesting if you obtain the ray traced image using CUDA and CUDA Driver API.

CUDA and OpenGL interoperability using CUDA Driver API

Before we compose CUDA ray traced content with OpenGL rasterised content we need to obtain the former. If you’re a CUDA programmer you know that CUDA does not posses any native graphics output and You have to manually take care of displaying the generated data. One of the solutions for this is to write pixel data into a OpenGL Pixel Buffer Object (PBO), copying the content of PBO to an OpenGL texture, and rendering a quad textured using this texture into the viewport. Fortunately for us CUDA and OpenGL interoperability is quite robust and this can be done quite easily, even when using lower-level CUDA Driver API.

// general initialisation
init_CUDA_and_OpenGL();
set_up_camera(camera);

GLuint pbo; // Pixel buffer object
GLuint tex; // Texture
unsigned int width = 512; // horizontal render resolution
unsigned int height = 512; // vertical render resolution

// Init PBO
glGenBuffersARB(1, &pbo);
glBindBufferARB(GL_PIXEL_UNPACK_BUFFER_ARB, pbo);
glBufferDataARB(GL_PIXEL_UNPACK_BUFFER_ARB, width*height*sizeof(GLubyte)*4, 0, GL_STREAM_DRAW_ARB);
glBindBufferARB(GL_PIXEL_UNPACK_BUFFER_ARB, 0);

// Init Texture
glGenTextures(1, &tex);
glBindTexture(GL_TEXTURE_2D, tex);
glTexImage2D(GL_TEXTURE_2D, 0, GL_RGBA8, width, height, 0, GL_RGBA, GL_UNSIGNED_BYTE, NULL);
glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_NEAREST);
glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_NEAREST);
glBindTexture(GL_TEXTURE_2D, 0);

// Allocate CUDA graphic Resource (Notice that I'm using Driver API function calls)
CUgraphicsResource * cuda_pbo_resource = (CUgraphicsResource*) malloc(sizeof(CUgraphicsResource));
cuGraphicsGLRegisterBuffer(cuda_pbo_resource, // returns pointer to the CUDA resource
 pbo, // we pass here the OpenGL resource
 CU_GRAPHICS_MAP_RESOURCE_FLAGS_WRITE_DISCARD); // Flag to aid optimisation

Rendering pass #1 — Ray tracing in CUDA

We render our scene (only solid objects, i.e. no lines etc.) using ray tracing. We store the raw pixel data in a PBO mapped as CUDA resource.
Note that PBO resource has to be unmapped from CUDA for OpenGL to be able to reliably access it. Thus we map the PBO resource to CUDA only for the time of ray tracing.

// map PBO resource
cuGraphicsMapResources(1, cuda_pbo_resource, 0); //mapping is only temporary

// get pointer in CUDA memory space for to the resource
CUdeviceptr d_pbo; // CUDA pointer through which the mapped graphics resources may be accessed
size_t num_bytes; //size of memory which may be accessed from that pointer
cuGraphicsResourceGetMappedPointer(&d_pbo, &num_bytes, *cuda_pbo_resource); // the d_pbo will be passed to rendering kernel

// Start ray tracing kernel
start_CUDA_ray_tracing_kernel(d_pbo, num_bytes); // render to PBO. Pass the pointer to the memory space where the data is to be saved

// unmap PBO resource
cuGraphicsUnmapResources(1, cuda_pbo_resource, 0);

After the PBO is filled with data by CUDA kernel we copy the content of the PBO to an OpenGL texture. On current hardware this step takes less than a milisecond for a 24 bits per pixel with a 512×512 resolution, thus it’s quite fast. Alternatively. You could attempt to render from CUDA directly to OpenGL texture memory space.


// copy from PBO to texture
glPixelStorei(GL_UNPACK_ALIGNMENT, 1);
glBindBufferARB(GL_PIXEL_UNPACK_BUFFER_ARB, pbo);
glBindTexture(GL_TEXTURE_2D, tex);  // bind our texture!
glTexSubImage2D(GL_TEXTURE_2D, 0, 0, 0, width, height, GL_RGBA, GL_UNSIGNED_BYTE, 0); // copy data
glBindBufferARB(GL_PIXEL_UNPACK_BUFFER_ARB, 0); // unbind PBO

Displaying CUDA rendered content with OpenGL

We can now render the texture onto a quad which is filling the viewport (using orthogonal projection).


// select matrix mode
glMatrixMode(GL_PROJECTION); // we select to  work with projection matrix stack
glLoadIdentity(); // make sure the there is nothing on the stack

glOrtho(0,1, 0,1,
0,             // zNear
999.0f);    // zFar

// clear color and depth buffer
glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
glDisable(GL_DEPTH_TEST); // we dont need it right now

// draw binded textured to a quad filling the viewport
glEnable(GL_TEXTURE_2D);
glBegin(GL_QUADS);
	glTexCoord2f(0, 0); glVertex2f(0, 0);
	glTexCoord2f(1, 0); glVertex2f(1, 0);
	glTexCoord2f(1, 1); glVertex2f(1, 1);
	glTexCoord2f(0, 1); glVertex2f(0, 1);
glEnd();

glDisable(GL_TEXTURE_2D);
glBindTexture(GL_TEXTURE_2D, 0);  // unbind texture

Rendering pass #2 — Rasterisation with OpenGL

We render the whole 3D scene (solid object only) once more, this time using OpenGL. You would probably want to switch the view mode from orthogonal to perspective projection with the same parameters as you used in your ray tracer (i.e. camera position, FOV, camera ‘look at’ point). I recommend setting the perspective manually using glFrustum() as it allows to set the camera parameters in a similar manner as you have probably used them in the CUDA kernel for ray tracing.

Note, that now we want to render the scene only to the OpenGL depth buffer — set you color buffer channel’s masks to false, and have depth testing enabled. Also, we will not need any OpenGL lights, textures, fog etc. — have them disabled.  As we are rasterizing only depth information this additional rendering pass is quite fast even for complex scenes.


// set perspective to the same as in the ray tracing pass
GLdouble fl_left = camera->plane_height;    // width of left
GLdouble fl_right = -camera->plane_height; // width of right (note the MINUS)
GLdouble fl_top = camera->plane_width;	// width of top
GLdouble fl_bottom = -camera->plane_width; // width of bottom (note the MINUS)
GLdouble zNear = camera->plane_distance;  // near clipping distance
GLdouble zFar = camera->plane_distance + 100; // far clipping distance

glFrustum(fl_left, fl_right, fl_bottom, fl_top, zNear, zFar);

gluLookAt(
 // camera position
 camera->origin.x,
 camera->origin.y,
 camera->origin.z,
 // look_at is stored as a normal vector, thus we add it to the camera position to determine a point in space where to look
 camera->origin.x + camera->look_at.x,
 camera->origin.y + camera->look_at.y,
 camera->origin.z + camera->look_at.z,
 // camera up direction
 camera->up.x,
 -camera->up.y,
 camera->up.z);

// render only to depth buffer
glColorMask(GL_FALSE,GL_FALSE,GL_FALSE,GL_FALSE); // we dont want to render to color buffer
glEnable(GL_DEPTH_TEST); // it is important to be enabled now

// Render the meshes in the scene to depth buffer only (do not render lines yet)
render_scene_objects_with_OpenGL();

Now, we want to enable color mask, and while having the depth testing still on, we rasterize remaining objects (lines, points etc.) to the color and depth buffers of OpenGL. However we can choose to also rasterize any additional solid objects at this point, basically any OpenGL scene!

glColorMask(GL_TRUE, GL_TRUE, GL_TRUE, GL_TRUE); // we now want to render colors
render_lines_and_misc_stuff();

As the depth buffer already contains scene depth information, rasterized object and ray traced objects will seem to occlude each other in a manner which you would expect. :-) You can also achieve a partial occlution by CUDA object (transparency effect), but first rasterizing the misc. OpenGL objects partly transparent with GL_DEPTH_TEST disabled, and later rendering them fully opaque with depth testing enabled.

Oh, don’t forget to swap buffers if you were using double buffering!

Limitations

The major limitation of this method is the obvious fact that the lines (or any misc OpenGL rendered content) present behind translucent object (which in ray tracing refract light) won’t appear refracted. Same goes for reflections. Fortunately, for transformation gismos, bounding boxes etc. this is not really a problem. It might be even desirable.

Additional resources

You might also want to check this resource for more on depth composition of images in OpenGL.

More on PBOs can be found here and here.

Triers CUDA ray tracing tutorial — a great tutorial with source code available.

Hope this helps.

1,047 ResponsesLeave one →

  1. [url=http://sildenafilcitrate100.online/]sildenafil 100 mg lowest price[/url]

    Reply
  2. [url=https://sildenafilgeneric.online/]can you buy sildenafil over the counter[/url]

    Reply
  3. [url=https://ivertabs.online/]stromectol[/url]

    Reply
  4. [url=http://iverm.online/]buy oral ivermectin[/url]

    Reply
  5. [url=https://buyivermectin.monster/]stromectol ivermectin[/url]

    Reply
  6. [url=https://ekb-sharik.ru/]гелевые шары в ульяновске[/url]

    Reply
  7. SyaeiertBic

     /  2021/11/26

    Насосы: тепловые насосы – циркуляционные насосы для вашей системы отопления.
    Это душа всякий системы отопления: циркуляционный насос, нередко также называемый тепловым насосом. Единовластно от того, работает ли учение отопления на газе, масле, тепловом насосе или солнечной системе – без циркуляционного насоса нуль не работает. Он отвечает следовать циркуляцию отопительной воды в системе отопления и за то, для тепло доходило накануне потребителей тепла, таких будто радиаторы. Посмотри цена циркуляционный насос для отопления как выбрать смотри на сайте [url=http://www.nasoscirkulyacionnyi.ru/]nasoscirkulyacionnyi[/url] циркуляционные насосы какой выбрать. Поскольку насосы должны трудолюбивый прежде 5000 часов в год, они представляют собой одного из крупнейших потребителей энергии в домашнем хозяйстве, следовательно важно, дабы они соответствовали последнему слову техники, поскольку современные высокоэффективные насосы потребляют прибл. Для 80% меньше электроэнергии. Впрочем большинство старых насосов являются устаревшими и крупногабаритными, поэтому стоит взглянуть на котельную!
    Циркуляционные насосы чтобы отопления используются ради циркуляции отопительной воды в замкнутой системе отопления. Они надежно гарантируют, что теплая вода через теплогенератора направляется в радиаторы сиречь в систему теплого пола и возвращается через обратку. Циркуляционный насос вынужден исключительно преодолеть гидравлическое оппозиция арматуры и трубопроводов из-за замкнутой системы.
    Экономьте электроэнергию с помощью современных циркуляционных насосов.
    Современные циркуляционные насосы для систем отопления отличаются высокой эффективностью и низким энергопотреблением – всего порядочно ватт во пора работы. По сравнению с устаревшими насосами прошлых десятилетий существует огромный потенциал экономии энергии. Именно следовательно стоит заменить извечный насос отопления для современный высокоэффективный насос, что часто оживленно окупается после счет небольшого энергопотребления. Циркуляционные насосы для систем отопления доступны в различных исполнениях, которые могут разниться сообразно максимальному напору и по резьбе подключения.
    Циркуляционный насос также часто называют тепловым насосом или насосом чтобы теплого пола. По сути, это одни и те же водяные насосы, однако названы по-разному из-за их применения. В зависимости от области применения и, в частности, присоединительных размеров, у нас вы легко найдете счастливый насос.
    Циркуляционный насос – это насос, обеспечивающий снабжение систем отопления (теплой) водой. Циркуляционный насос можно найти в каждом доме либо здании, где ужинать радиаторы и / сиречь полы с подогревом. Когда ужинать и теплый пол, и радиаторы, то циркуляционных насосов несколько. В больших домах / зданиях в большинстве случаев также устанавливаются скольконибудь циркуляционных насосов ради обеспечения полной циркуляции. Действительно вы можете говорить, который циркуляционный насос – это сердце вашей системы центрального отопления или системы теплых полов. Ежели центральное отопление или полы с подогревом не работают сиречь работают плохо, это многократно связано с циркуляционным насосом, что необходимо заменить. В зависимости от системы отопления вы можете легко встречать циркуляционный насос за передней панелью вашего центрального отопления. Примерно во всех случаях размеры насоса, которым должен идти новый циркуляционный насос из нашего ассортимента, указаны на вашем старом насосе. Если это не так, вы можете свободно найти марку и серийный комната в таблице обмена . Наши инструкции, приведенные ниже, помогут вам заменить насос отопления.
    Мы выбрали три насоса (центрального) отопления иначе циркуляционные насосы чтобы нашего ассортимента. Обладая этими тремя размерами, мы можем предложить замену теплового насоса чтобы каждой установки, с через которой вы также напрямую сэкономите на расходах на электроэнергию. Подробнее об этом ниже! Исключая того, заменить циркуляционный насос, насос отопления либо насос теплого пола шабаш просто. Это экономит затраты на установку. Ознакомьтесь с инструкциями и / иначе видео ниже. Разумеется, вы также можете заказать установку теплового насоса у специалиста. Однако чтобы экономии средств рекомендуем приказывать циркуляционный насос у нас по конкурентоспособной цене. Тутто вы платите только следовать установку у специалиста.

    Reply
  8. [url=https://malegra.online/]malegra on line[/url]

    Reply
  9. [url=http://cleocinbuy.online/]clindamycin tablet coupon[/url]

    Reply
  10. [url=https://tadalafil.men/]tadalafil online us[/url]

    Reply
  11. [url=http://ivermectin.media/]ivermectin 4000 mcg[/url]

    Reply
  12. [url=https://triamterene.quest/]triamterene hctz 37.5 25 mg[/url]

    Reply
  13. [url=http://tadalafilgmed.com/]tadalafil generic india[/url]

    Reply
  14. [url=http://budesonide.quest/]budecort 1mg[/url]

    Reply
  15. [url=http://ivermectin.red/]stromectol 3mg[/url]

    Reply
  16. [url=https://buyhydroxychloroquine.quest/]hydroxychloroquine plaquenil cost[/url]

    Reply
  17. [url=https://sildenafil.today/]sildenafil 50mg[/url]

    Reply
  18. [url=https://ivermectin.gen.in/]ivermectin 2mg[/url]

    Reply
  19. [url=http://buyamoxicillin.online/]amoxicillin best price[/url]

    Reply
  20. [url=http://ivermectin.party/]stromectol price[/url]

    Reply
  21. [url=https://ivermecfin.com/]ivermectin 5 mg[/url]

    Reply
  22. [url=http://iveromectinforhumans.com/]ivermectin uk buy[/url]

    Reply
  23. [url=https://tadalafilwithoutdoctor.com/]buy cheap tadalafil[/url]

    Reply
  24. [url=https://buyiverpill.com/]stromectol for humans[/url]

    Reply
  25. [url=https://sildenafilcitrat.com/]sildenafil 100mg price[/url]

    Reply
  26. [url=http://ivermectinforsale.online/]buy ivermectin cream for humans[/url]

    Reply
  27. [url=https://tadalafilwithoutdoctor.com/]buy tadalafil online canada[/url]

    Reply
  28. [url=http://tizanidine.quest/]tizanidine 2mg capsules[/url]

    Reply
  29. [url=https://axsildenafil.com/]where can i buy over the counter sildenafil[/url]

    Reply
  30. [url=https://plavix.quest/]buy clopidogrel 75 mg[/url]

    Reply
  31. [url=http://buyantibiotics.online/]ceftin price[/url]

    Reply
  32. [url=http://tadalafilng.com/]buy generic tadalafil uk[/url]

    Reply
  33. [url=https://disulfiram.online/]disulfiram 250 mg price[/url]

    Reply
  34. [url=http://singulair.quest/]where to buy singulair[/url]

    Reply
  35. [url=https://yandex.ru/maps/org/49926867232]шары в ульяновске воздушные[/url]

    Reply
  36. Thank you, A good amount of stuff.
    do my essay cheap [url=https://englishessayhelp.com/]argumentative essay[/url] websites that write essays for you

    Reply
  37. Whoa tons of great data!
    nyu college essay [url=https://paperwritingservicestops.com/]https://paperwritingservicestops.com/[/url] website content writing

    Reply
  38. You’ve made your stand pretty nicely.!
    online essay writers wanted [url=https://quality-essays.com/]what is college essay[/url] online essay writing service review

    Reply
  39. Nicely put. Cheers.
    how to write a bibliography for an essay [url=https://theessayswriters.com/]https://theessayswriters.com/[/url] medical residency personal statement writing services

    Reply
  40. [url=https://maps.google.com/maps?cid=7613754105589455939]гелиевые шары ульяновск дешево адреса[/url]

    Reply
  41. [url=http://brandandgeneric.com/]stromectol how much it cost[/url]

    Reply

Leave a Reply