untitled.m

本程序是针对sar的原始数据到成像后幅度图像的成像过程。。

clear all

%    disp(' ');
%    chip_name = input('Name of  target chips: ','s');

%    fid = 1;
t1='HB0';
t2='.015';
% cyc=0;
% for ii=3333:5684
ii=3429;
    t3=sprintf('%d',ii);
    files=[t1 t3 t2];
    for i = 1:100
        header(i,:) = blanks(100);
    end

    i = 0;
    tp=[];
    fid = fopen(files,'r');
%     if fid==-1
%         continue
%     end
%     %fid
%     cyc=1+cyc;
    while (strcmp(tp,'[EndofPhoenixHeader]') == 0)

        z1 = [];

        z1 = fgets(fid);

        tp = sscanf(z1,'%s');


        if(isempty(tp) == 0)
            i = i+1;
            header(i,:) = zeros(1,100);
            header(i,1:(size(tp,2))) = tp;
        end

    end

    hdr_size_field = 'PhoenixHeaderLength=';
    hdr_size_flag = 0;
    i = 0;

    while(hdr_size_flag == 0)
        i = i+1;
        hdr_size_flag = strcmp(header(i,1:size(hdr_size_field,2)), hdr_size_field);
    end

    hdrsize = str2num(header(i,size(hdr_size_field,2)+1:size(header,2)));



    numcol_field = 'NumberOfColumns=';
    numcol_flag = 0;
    i = 0;

    while(numcol_flag == 0)
        i = i+1;
        numcol_flag = strcmp(header(i,1:size(numcol_field,2)),numcol_field);
    end

    numcol = str2num(header(i,size(numcol_field,2)+1:size(header,2)));

    numrow_field = 'TargetAz=';
    numrow_flag = 0;
    i = 0;

    while(numrow_flag == 0)
        i = i+1;
        numrow_flag = strcmp(header(i,1:size(numrow_field,2)),numrow_field);
    end

    TargetAz = str2num(header(i,size(numrow_field,2)+1:size(header,2)));
    %     TargetAz

    numrow_field = 'NumberOfRows=';
    numrow_flag = 0;
    i = 0;

    while(numrow_flag == 0)
        i = i+1;
        numrow_flag = strcmp(header(i,1:size(numrow_field,2)),numrow_field);
    end

    numrow = str2num(header(i,size(numrow_field,2)+1:size(header,2)));

    %i=1;
    fclose(fid);



    %     disp(['Processing chip image: ', chip_name]);

    fid = fopen(files,'r','b');
    fseek(fid,hdrsize,'bof');

    [mag_vector,count] = fread(fid,[numrow*numcol],'float32');
    imgdata(:,1:numrow) = reshape(mag_vector,numcol,numrow);%读出的幅度图像
    fclose(fid);
    imshow(imgdata);