Changeset 993 for trunk/src/read_cine_phantom.m

Timestamp:

Feb 6, 2017, 11:29:01 AM (7 years ago)

Author:

sommeria

Message:

various updates

File:

• trunk/src/read_cine_phantom.m (modified) (3 diffs)

trunk/src/read_cine_phantom.m

@@ -3 +3 @@
 
 
-[CineFileHeader, BitmapInfoHeader, CameraSetup, imageLocations] = readCineHeader(cineFilePath);
+% [CineFileHeader, BitmapInfoHeader, CameraSetup, imageLocations] = readCineHeader(cineFilePath);       % old
+[CineFileHeader, BitmapInfoHeader, CameraSetup, TimeOnlyBlock, ExposureOnlyBlock, TimeCodeBlock, imageLocations] = readCineHeader(cineFilePath);
 fid = fopen(cineFilePath);
 lookupTable = lookupTablePackedFun;
@@ -22 +23 @@
     fseek(fid, imageLocations( frames(ii) ), 'bof');
     if ~BitmapInfoHeader.biCompression
-%         [A,count]=fread(fid, 10000)
         imTemp = fread(fid, [BitmapInfoHeader.biWidth BitmapInfoHeader.biHeight],'uint16');
         imMat(:,:,ii) = imTemp';
     else
-        imTemp = fread(fid, [BitmapInfoHeader.biWidth BitmapInfoHeader.biHeight], 'ubit10','b');
-        im = imTemp';
-        im( im < 1 ) = 1;
-        im =reshape( interp1( 1:1024, lookupTable, im(:) ), BitmapInfoHeader.biHeight, BitmapInfoHeader.biWidth);
-        im( im < CameraSetup.BlackLevel ) = CameraSetup.BlackLevel;
-        im( im > 4064 )=4064;
-        imMat(:,:,ii) = reshape(interp1( CameraSetup.BlackLevel:4064, linspace(0,4095,4064 - CameraSetup.BlackLevel+1), im(:)),...
-            BitmapInfoHeader.biHeight, BitmapInfoHeader.biWidth);
-        
-        if mod( round(0.01*Nf), ii)
-            if ~exist('dispStr', 'var')
-                dispStr = ' ';
-                disp( dispStr )
-            end
-            
-            lenDispStr = length( dispStr );
-            dispStr = ['Reading is ' num2str( round( 100*ii/Nf ) ) '% complete'];
-            disp( [char(8)*ones(1,lenDispStr+1) dispStr] )
-        end
+    imTemp = fread(fid, [BitmapInfoHeader.biWidth BitmapInfoHeader.biHeight], 'ubit10','b');
+    im = imTemp';
+    im( im < 1 ) = 1;
+    im =reshape( interp1( 1:1024, lookupTable, im(:) ), BitmapInfoHeader.biHeight, BitmapInfoHeader.biWidth);
+    im( im < CameraSetup.BlackLevel ) = CameraSetup.BlackLevel;
+    im( im > 4064 )=4064;
+    imMat(:,:,ii) = reshape(interp1( CameraSetup.BlackLevel:4064, linspace(0,4095,4064 - CameraSetup.BlackLevel+1), im(:)),...
+        BitmapInfoHeader.biHeight, BitmapInfoHeader.biWidth); 
+    
+    if mod( round(0.01*Nf), ii)
+       if ~exist('dispStr', 'var')
+           dispStr = ' ';
+           disp( dispStr )
+       end
+       
+     lenDispStr = length( dispStr );
+     dispStr = ['Reading is ' num2str( round( 100*ii/Nf ) ) '% complete'];
+     disp( [char(8)*ones(1,lenDispStr+1) dispStr] )
+    end
     end
 end
@@ -51 +51 @@
 end
 
-function lookupTable = lookupTablePackedFun()
+function lookupTable = lookupTablePackedFun() 
 % function to transform the compressed 10 bit images back, close to the 12 bit camera images
 lookupTable = [ 2,   5,   6,   7,   8,   9,  10,  11,  12,  13,  14,  15,  16,  17,  17,  18,...