field ii 肾脏成像

field II的肾脏成像 源程序，包括field 的全套工具包，需要的可以下载下的。其中包括肾脏模型文件。

function [Xmapout3]=bmpread（filename）;
%BMPREAD Read a BMP （Microsoft Windows Bitmap） file from disk.
%	[XMAP]=BMPREAD（‘filename‘） reads the file ‘filename‘ and
%	returns the indexed image X and associated colormap
%	MAP. If no extension is given for the filename the
%	extension ‘.bmp‘ is assumed.
%	[RGB]=BMPREAD（‘filename‘） reads the 24-bit BMP file
%	from the file ‘filename‘.
%	BPP=BMPREAD（‘filename‘） returns the number of bits per
%	pixel in the BMP file.
%
%	BMPREAD does not read 1-bit or compressed BMP files.
%
%	See also: BMPWRITE GIFREAD HDFREAD PCXREAD TIFFREAD
%	          XWDREAD.

%	Mounil Patel 3/10/94
%	Revised Steve Eddins February 9 1995
%	Copyright （c） 1994 by The MathWorks Inc.
%	$Revision: 1.13 $  $Date: 1995/02/09 15:17:43 $


if （nargin~=1）
	error（‘Requires a filename as an argument.‘）;
end;

if （isstr（filename）~=1）
	error（‘Requires a string filename as an argument.‘）;
end;

if （isempty（findstr（filename‘.‘））==1）
	filename=[filename‘.bmp‘];
end;

fid=fopen（filename‘rb‘‘l‘）;
if （fid==-1）
	error（[‘Error opening ‘filename‘ for input.‘]）;
end;

bfType=fread（fid2‘uchar‘）;
if （bfType~=[66;77]）
	fclose（fid）;
	error（‘Not a BMP file.‘）;
end;

bfSize=fread（fid1‘uint‘）;
bfReserved1=fread（fid1‘ushort‘）;
bfReserved2=fread（fid1‘ushort‘）;
bfOffBytes=fread（fid1‘uint‘）;

biSize=fread（fid1‘uint‘）;
if （biSize~=40）
	fclose（fid）;
	error（‘Not a MS Windows Device Independent Bitmap BMP file.‘）;
end;

biWidth=fread（fid1‘uint‘）;
biHeight=fread（fid1‘uint‘）;
biPlanes=fread（fid1‘ushort‘）;
biBitCount=fread（fid1‘ushort‘）;
if （biBitCount == 1）
  error（‘BMPREAD does not read 1-bit BMP files.‘）;
end
biCompression=fread（fid1‘uint‘）;
if biCompression~=0 then
	error（‘Can‘‘t load compressed format bitmaps.‘）;
end;
biSizeImage=fread（fid1‘uint‘）;
biXPels=fread（fid1‘uint‘）;
biYPels=fread（fid1‘uint‘）;
biClrUsed=fread（fid1‘uint‘）;
biClrImportant=fread（fid1‘uint‘）;

if （nargout <= 1）
  X = biBitCount;
  map = [];
  fclose（fid）;
  return;
end

if （（nargout == 2） & （biBitCount == 24））
  error（‘Three output arguments required for 24-bit BMP file.‘）;
end

if （（nargout == 3） & （biBitCount < 24））
  error（‘Only two arguments needed for 4- and 8-bit BMP files.‘）;
end

if （biBitCount == 24）
  if （rem（biWidth4）~=0）
    XSize=biWidth+（4-rem（biWidth4））;
  else
    XSize=biWidth;
  end
  YSize=biHeight;
  Size=XSize*YSize;

  fseek（fid bfOffBytes ‘bof‘）;
  X=fread（fidSize*3‘uchar‘）;
  out3 = rot90（reshape（X（1:3:length（X）） XSize YSize）） * （1/255）;
  map = rot90（reshape（X（2:3:length（X）） XSize YSize）） * （1/255）;
  X = rot90（reshape（X（3:3:length（X）） XSize YSize）） * （1/255）;
  if （rem（biWidth4）~=0）
    X=X（:1:biWidth）;
    map = map（:1:biWidth）;
    out3 = out3（:1:biWidth）;
  end;
  
  fclose（fid）;
  return;
end  

if （biClrUsed==0）
	nColors=2.^biBitCount;
else
	nColors=biClrUsed;
end;

% load color map now
if （biClrUsed>256）
	map=[];	% 24-bit images have no colormap
else
	map=fread（fid4*nColors‘uchar‘）;
	map=reshape（map

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2014-05-13 14:33  kidney\
     文件      115200  2010-11-19 23:32  kidney\Mat_field.mexw32
     文件      157184  2010-11-19 23:39  kidney\Mat_field.mexw64
     文件        3981  2014-03-25 12:47  kidney\bmpread.m
     文件         792  1997-05-21 11:53  kidney\calc_h.m
     文件         846  1997-05-21 11:53  kidney\calc_hhp.m
     文件         764  1998-03-19 13:22  kidney\calc_hp.m
     文件        1197  1997-05-21 11:53  kidney\calc_scat.m
     文件        2504  2001-08-17 17:07  kidney\calc_scat_all.m
     文件        1386  1999-05-21 17:08  kidney\calc_scat_multi.m
     文件        1152  1998-07-01 08:22  kidney\ele_apodization.m
     文件        1137  1998-07-01 08:22  kidney\ele_delay.m
     文件        1143  1998-07-01 08:22  kidney\ele_waveform.m
     文件       15700  2014-03-25 14:16  kidney\fast_int.c
     文件       31232  2014-03-27 16:50  kidney\fast_int.dll
     文件       10502  2014-03-24 21:38  kidney\fast_int.mexglx
     文件       10332  2014-03-24 21:42  kidney\fast_int.mexlx
     文件        9284  2014-03-25 12:32  kidney\fast_int.mexmac
     文件        9308  2014-03-24 21:43  kidney\fast_int.mexsol
     文件       20480  2014-03-25 15:07  kidney\fast_int.mexw32
     文件         417  1997-05-21 11:53  kidney\field_debug.m
     文件         295  1997-05-21 11:53  kidney\field_end.m
     文件         883  2011-05-06 22:34  kidney\field_guide.m
     文件         370  2001-12-17 17:24  kidney\field_info.m
     文件         826  1999-01-21 10:08  kidney\field_init.m
     文件         393  2007-08-10 21:17  kidney\field_logo.m
     文件        1700  2014-03-24 21:44  kidney\human_kidney_phantom.m
     文件      246782  2014-03-28 15:19  kidney\kidney_cut.bmp
     文件       76373  2014-03-26 17:13  kidney\kidney_phantom_example.tar.gz
     文件      160979  2014-03-29 11:02  kidney\kidney_scan.jpg
     文件      285392  2001-08-21 13:44  kidney\logo_field.mat
............此处省略170个文件信息