matlab读取raw格式文件

function [imagedata] = dpcrawprocessing（filename）;
%fnum = str2num（filename（1））;
dot = max（find（filename == ‘.‘））;

suffix = filename（dot+1:dot+3）;% 提取图像文件后缀
name = filename（1:dot-1）;% 提取图像文件名

if strcmp（suffix‘raw‘）   % 要求是raw格式

    disp（sprintf（‘\nOpens %s file\n‘filename））;

    fp = fopen（filename）;  % “Big-endian“ byte order.

    if （fp<0）

        error（[‘Cannot open ‘ filename ‘.‘]）;

    end
	
    % Creates a gray palette and scale it to [01].
    %disp（‘Creating gray palette...‘）;
    % 灰度映射表
   % for k=1:512
     %   map（k[1:2]）=[k/512k/512];
   % end;
    
    %  读取raw文件，并以矩阵显示
	image = fread（fp‘uint16‘）;
	imagesize = sqrt（size（image（:1）））;
	for i=1:imagesize
		for j= 1:imagesize