Carving a Dinosaur: spacecarving.project( camera, world_X, world_Y, world_Z ) - File Exchange

Code covered by the BSD License

Highlights from
Carving a Dinosaur

makeFullAxes( axh ) makeFullAxes: make some axes fill the available space
spacecarving.carve( voxels, c... CARVE: remove voxels that are not in the silhouette
spacecarving.carveall( voxels... CARVEALL carve away voxels using all cameras
spacecarving.colorsurface( pt... COLORSURFACE: color in a surface based on image data
spacecarving.decomposeP(P) DECOMPOSEP decompose a projection matrix into internal and external params
spacecarving.findmodel( camer... FINDMODEL: locate the model to be carved relative to the cameras
spacecarving.getcameradirecti... GETCAMERADIRECTION: return the view direction of a camera
spacecarving.getsilhouette( i... GETSILHOUETTE - find the silhouette of an object centered in the image
spacecarving.loadcameradata(d... LOADCAMERADATA: Load the dino data
spacecarving.makevoxels(xlim,... MAKEVOXELS create a basic grid of voxels ready for carving
spacecarving.project( camera,... PROJECT: project a 3D point into an image
spacecarving.showcamera(ax,ca... SHOWCAMERA: draw a schematic of a camera
spacecarving.showscene(camera... SHOWSCENE: show a carve scene, including cameras, images and the model
spacecarving.showsurface( vox... SHOWSURFACE: draw a surface based on some voxels
tmwMultiWaitbar( label, value... TMWMULTIWAITBAR: add, remove or update an entry on the multi waitbar
space_carving_demo.m
View all files

from Carving a Dinosaur by Ben Tordoff
This is a demo of reconstructing a 3D shape from multiple images using space-carving.

spacecarving.project( camera, world_X, world_Y, world_Z )

function [im_x, im_y] = project( camera, world_X, world_Y, world_Z )
%PROJECT: project a 3D point into an image
%
%   [IM_X,IM_Y] = PROJECT(CAMERA,WORLD_X,WORLD_Y,WORLD_Z) projects one or
%   more 3D point in the world coordinate frame into image coordinates

%   Copyright 2005-2009 The MathWorks, Inc.
%   $Revision: 1.0 $    $Date: 2006/06/30 00:00:00 $

z = camera.rawP(3,1) * world_X ...
    + camera.rawP(3,2) * world_Y ...
    + camera.rawP(3,3) * world_Z ...
    + camera.rawP(3,4);
im_y = round( (camera.rawP(2,1) * world_X ...
    + camera.rawP(2,2) * world_Y ...
    + camera.rawP(2,3) * world_Z ...
    + camera.rawP(2,4)) ./ z);
im_x = round( (camera.rawP(1,1) * world_X ...
    + camera.rawP(1,2) * world_Y ...
    + camera.rawP(1,3) * world_Z ...
    + camera.rawP(1,4)) ./ z);

Highlights from Carving a Dinosaur

Highlights from
Carving a Dinosaur