1 files changed, 109 insertions, 0 deletions

diff --git a/SD-VBS/common/toolbox/toolbox_basic/TOOLBOX_calib/loadppm.m b/SD-VBS/common/toolbox/toolbox_basic/TOOLBOX_calib/loadppm.m
new file mode 100755
index 0000000..0c004fc
--- /dev/null
+++ b/SD-VBS/common/toolbox/toolbox_basic/TOOLBOX_calib/loadppm.m
@@ -0,0 +1,109 @@
+%LOADPPM        Load a PPM image
+%
+%       I = loadppm(filename)
+%
+%       Returns a matrix containing the image loaded from the PPM format
+%       file filename.  Handles ASCII (P3) and binary (P6) PPM file formats.
+%
+%       If the filename has no extension, and open fails, a '.ppm' and
+%       '.pnm' extension will be tried.
+%
+% SEE ALSO:     saveppm loadpgm
+%
+%       Copyright (c) Peter Corke, 1999  Machine Vision Toolbox for Matlab
+
+
+% Peter Corke 1994
+
+function I = loadppm(file)
+        white = [' ' 9 10 13];  % space, tab, lf, cr
+        white = setstr(white);
+
+        fid = fopen(file, 'r');
+        if fid < 0,
+                fid = fopen([file '.ppm'], 'r');
+        end
+        if fid < 0,
+                fid = fopen([file '.pnm'], 'r');
+        end
+        if fid < 0,
+                error('Couldn''t open file');
+        end
+                
+        magic = fread(fid, 2, 'char');
+        while 1
+                c = fread(fid,1,'char');
+                if c == '#',
+                        fgetl(fid);
+                elseif ~any(c == white)
+                        fseek(fid, -1, 'cof');  % unputc()
+                        break;
+                end
+        end
+        cols = fscanf(fid, '%d', 1);
+        while 1
+                c = fread(fid,1,'char');
+                if c == '#',
+                        fgetl(fid);
+                elseif ~any(c == white)
+                        fseek(fid, -1, 'cof');  % unputc()
+                        break;
+                end
+        end
+        rows = fscanf(fid, '%d', 1);
+        while 1
+                c = fread(fid,1,'char');
+                if c == '#',
+                        fgetl(fid);
+                elseif ~any(c == white)
+                        fseek(fid, -1, 'cof');  % unputc()
+                        break;
+                end
+        end
+        maxval = fscanf(fid, '%d', 1);
+        while 1
+                c = fread(fid,1,'char');
+                if c == '#',
+                        fgetl(fid);
+                elseif ~any(c == white)
+                        fseek(fid, -1, 'cof');  % unputc()
+                        break;
+                end
+        end
+        if magic(1) == 'P',
+                if magic(2) == '3',
+                        %disp(['ASCII PPM file ' num2str(rows) ' x ' num2str(cols)])
+                        I = fscanf(fid, '%d', [cols*3 rows]);
+                elseif magic(2) == '6',
+                        %disp(['Binary PPM file ' num2str(rows) ' x ' num2str(cols)])
+                        if maxval == 1,
+                                fmt = 'unint1';
+                        elseif maxval == 15,
+                                fmt = 'uint4';
+                        elseif maxval == 255,
+                                fmt = 'uint8';
+                        elseif maxval == 2^32-1,
+                                fmt = 'uint32';
+                        end
+                        I = fread(fid, [cols*3 rows], fmt);
+                else
+                        disp('Not a PPM file');
+                end
+        end
+        %
+        % now the matrix has interleaved columns of R, G, B
+        %
+        I = I';
+        size(I);
+        R = I(:,1:3:(cols*3));
+        G = I(:,2:3:(cols*3));
+        B = I(:,3:3:(cols*3));
+        fclose(fid);
+   
+   
+   I = zeros(rows,cols,3);
+   I(:,:,1) = R;
+   I(:,:,2) = G;
+   I(:,:,3) = B;
+   I = uint8(I);
+   
\ No newline at end of file