Changeset 993 for trunk/src/read_cine_phantom.m
- Timestamp:
- Feb 6, 2017, 11:29:01 AM (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/read_cine_phantom.m
r991 r993 3 3 4 4 5 [CineFileHeader, BitmapInfoHeader, CameraSetup, imageLocations] = readCineHeader(cineFilePath); 5 % [CineFileHeader, BitmapInfoHeader, CameraSetup, imageLocations] = readCineHeader(cineFilePath); % old 6 [CineFileHeader, BitmapInfoHeader, CameraSetup, TimeOnlyBlock, ExposureOnlyBlock, TimeCodeBlock, imageLocations] = readCineHeader(cineFilePath); 6 7 fid = fopen(cineFilePath); 7 8 lookupTable = lookupTablePackedFun; … … 22 23 fseek(fid, imageLocations( frames(ii) ), 'bof'); 23 24 if ~BitmapInfoHeader.biCompression 24 % [A,count]=fread(fid, 10000)25 25 imTemp = fread(fid, [BitmapInfoHeader.biWidth BitmapInfoHeader.biHeight],'uint16'); 26 26 imMat(:,:,ii) = imTemp'; 27 27 else 28 29 30 31 32 33 34 35 BitmapInfoHeader.biHeight, BitmapInfoHeader.biWidth);36 37 38 39 40 41 42 43 44 45 46 28 imTemp = fread(fid, [BitmapInfoHeader.biWidth BitmapInfoHeader.biHeight], 'ubit10','b'); 29 im = imTemp'; 30 im( im < 1 ) = 1; 31 im =reshape( interp1( 1:1024, lookupTable, im(:) ), BitmapInfoHeader.biHeight, BitmapInfoHeader.biWidth); 32 im( im < CameraSetup.BlackLevel ) = CameraSetup.BlackLevel; 33 im( im > 4064 )=4064; 34 imMat(:,:,ii) = reshape(interp1( CameraSetup.BlackLevel:4064, linspace(0,4095,4064 - CameraSetup.BlackLevel+1), im(:)),... 35 BitmapInfoHeader.biHeight, BitmapInfoHeader.biWidth); 36 37 if mod( round(0.01*Nf), ii) 38 if ~exist('dispStr', 'var') 39 dispStr = ' '; 40 disp( dispStr ) 41 end 42 43 lenDispStr = length( dispStr ); 44 dispStr = ['Reading is ' num2str( round( 100*ii/Nf ) ) '% complete']; 45 disp( [char(8)*ones(1,lenDispStr+1) dispStr] ) 46 end 47 47 end 48 48 end … … 51 51 end 52 52 53 function lookupTable = lookupTablePackedFun() 53 function lookupTable = lookupTablePackedFun() 54 54 % function to transform the compressed 10 bit images back, close to the 12 bit camera images 55 55 lookupTable = [ 2, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 17, 18,...
Note: See TracChangeset
for help on using the changeset viewer.