vb_get_vivid_proj

 Project SPM-R coordinate to 2D screen
  [V,F,C,X,Y,Z] = vb_get_vivid_proj(Coord, LR)
 --- Input
 Coord : Left / Right /Front  head coordinate
  'MGF'- File
        .LV   / .RV   / .FV   : 3D-Coordinate   : [Npoint 3]  [m]
        .LF   / .RF   / .FF   : Patch index     : [Npatch 4]
        .Limg / .Rimg / .Fimg : 3D RGB image    : [Npoint 3]
 LR    = 'L' or 'R' or 'F' : Left/Right/Front
 --- Output
 V : 3D coordinate corresponding to view point 'Left/Right/Front'
 X : X-coordinate in projected 2D screen-space
 Y : Y-coordinate in projected 2D screen-space
 Z : Z-coordinate = 0 
 F : patch index
 C : Color of vertex

 Copyright (C) 2011, ATR All Rights Reserved.
 License : New BSD License(see VBMEG_LICENSE.txt)

0001 function [V,F,C,X,Y,Z] = vb_get_vivid_proj(Coord, LR)
0002 % Project SPM-R coordinate to 2D screen
0003 %  [V,F,C,X,Y,Z] = vb_get_vivid_proj(Coord, LR)
0004 % --- Input
0005 % Coord : Left / Right /Front  head coordinate
0006 %  'MGF'- File
0007 %        .LV   / .RV   / .FV   : 3D-Coordinate   : [Npoint 3]  [m]
0008 %        .LF   / .RF   / .FF   : Patch index     : [Npatch 4]
0009 %        .Limg / .Rimg / .Fimg : 3D RGB image    : [Npoint 3]
0010 % LR    = 'L' or 'R' or 'F' : Left/Right/Front
0011 % --- Output
0012 % V : 3D coordinate corresponding to view point 'Left/Right/Front'
0013 % X : X-coordinate in projected 2D screen-space
0014 % Y : Y-coordinate in projected 2D screen-space
0015 % Z : Z-coordinate = 0
0016 % F : patch index
0017 % C : Color of vertex
0018 %
0019 % Copyright (C) 2011, ATR All Rights Reserved.
0020 % License : New BSD License(see VBMEG_LICENSE.txt)
0021 
0022 % SPM-R
0023 %   X: º¸->±¦   Left(-191/2)   -> Right(191/2)
0024 %   Y: ¸å->Á°   Back(-256/2)   -> Front(256/2)
0025 %   Z: ²¼->¾å   Bottom(-256/2) -> Top(256/2)
0026 %
0027 %view([-90, 0]);    % Left view  (from -X axis)
0028 %view([ 90, 0]);    % Right view (from +X axis)
0029 %view([180, 0]);    % Front view (from +Y axis)
0030 %view([  0, 90]);   % Top view   (from +Z axis)
0031 
0032 angleL = [-90, 0];
0033 angleR = [ 90, 0];
0034 angleF = [180, 0];
0035 
0036 switch    LR
0037 case    'L'
0038     V = Coord.LV;
0039     F = Coord.LF;
0040     C = Coord.Limg;
0041     T = viewmtx(angleL(1),angleL(2));
0042 case    'R'
0043     V = Coord.RV;
0044     F = Coord.RF;
0045     C = Coord.Rimg;
0046     T = viewmtx(angleR(1),angleR(2));
0047 case    'F'
0048     V = Coord.FV;
0049     F = Coord.FF;
0050     C = Coord.Fimg;
0051     T = viewmtx(angleF(1),angleF(2));
0052 end
0053 
0054 NV = size(V,1);
0055 X  = [V ones(NV,1)]*T(1,:)';
0056 Y  = [V ones(NV,1)]*T(2,:)';
0057 Z  = zeros(NV,1);
0058 
0059 return
0060 %%%%%%%%%%%%%%%%%%% END %%%%%%%%%%%%%%%%%%%%
0061 
0062 h = figure;
0063 
0064 subplot(1,2,1)
0065 % Plot face on projected 2D screen
0066 patch('Faces',F,'Vertices',[X Y Z],'FaceColor','interp',...
0067     'FaceVertexCData',C,'EdgeColor','none',...
0068     'FaceLighting','none');
0069 hold on
0070 view([ 0, 90]);
0071 axis equal
0072 
0073 subplot(1,2,2)
0074 % Plot 3D face
0075 patch('Faces',F,'Vertices',V,'FaceColor','interp',...
0076     'FaceVertexCData',C,'EdgeColor','none',...
0077     'FaceLighting','none');
0078 hold on
0079 view(T);
0080 axis equal
0081 
0082 return