untitled.m

来自「用于对天文图像进行去噪和探测」· M 代码 · 共 19 行

str{1}='G:\projects\matlab works\863-紫金山天文台星图处理070531\Data\Reduction\Image-793.fits';

Height=1024;
Width=1024;

file1 = fopen(str{1},'rb');    
fseek(file1,2880,'bof');

i44 = fread(file1,[Width,Height],'uint32','b');
% i44 = fread(file1,[Width,Height],'uint16');
fclose(file1);

i1 = uint32(i44');
i1=uint16(i1./2^16);

i1=imadjust(i1);
figure,imshow(i1,[]);