Changeset 991 for trunk/src/read_cine_phantom.m

Timestamp:

Feb 2, 2017, 1:25:44 PM (7 years ago)

Author:

sommeria

Message:

various

File:

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

trunk/src/read_cine_phantom.m

@@ -22 +22 @@
     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
@@ -50 +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,...