MSTAR 数据转格式

MSTAR是陆地目标的雷达图像数据集，但给出的文件是二进制格式的文件，在下载数据集是说明中给了一个地址可以下载到转换工具（https://www.sdms.afrl.af.mil，需翻墙），我用的不是很明白。

README中对MSTAR数据格式说明如下：



FILE FORMAT


Each file is constructed with a prepended, variable-length, 

Phoenix formatted (ASCII) header which contains detailed ground 

truth and sensor information for the specific chip.  Following 

the Phoenix header is the data block.  The data block is written 

in Sun floating point format and is divided into two blocks, a 

magnitude block followed by a phase block.



用记事本任意打开一个文件可以看到文件的PhoenixHeader，中间我们需要的信息是NumberOfColumns和NumberOfRows这两个变量，PhoenixHeader以单独一行“[EndofPhoenixHeader]”结束，后面跟的是模值块与相位块（SAR图像是一个复矩阵），我们需要提取模值块并保存为图片。

下面是代码：

  
  
clear
ReadPath = 'D:\Document\Graduation Project\MSTAR\DAA7B02AA\TARGETS\TEST\15_DEG\T72\SN_S7\';
SavePath = 'D:\Document\Graduation Project\MSTAR\DAA7B02AA\TARGETSJPG\TEST\15_DEG\T72\SN_S7\';
FileType = '*.017';
Files = dir([ReadPath FileType]);
NumberOfFiles = length(Files);
for i = 1 : NumberOfFiles
FileName = Files(i).name;
NameLength = length(FileName);
FID = fopen([ReadPath FileName], 'rb','ieee-be');
ImgColumns = 0;
ImgRows = 0;
while ~feof(FID) % 在PhoenixHeader找到图片尺寸大小
Text = fgetl(FID);
if ~isempty(strfind( Text, 'NumberOfColumns'))
ImgColumns = str2double( Text( 18: end));
Text = fgetl(FID);
ImgRows = str2double( Text( 15: end));
break;
end
end
while ~feof(FID) % 跳过PhoenixHeader
Text = fgetl(FID);
if ~isempty(strfind( Text, '[EndofPhoenixHeader]'))
break
end
end
Mag = fread(FID,ImgColumns*ImgRows, 'float32','ieee-be');
Img = reshape(Mag,[ImgColumns ImgRows]);
imwrite(uint8(imadjust(Img)* 255),[SavePath FileName( 1:NameLength -3) 'jpg']); % 调整对比度后保存
fclose(FID);
end
  
  

更新一下

最近又用到了这个小脚本，不过要自己建文件夹挺麻烦的，把它改了一下，写成了函数形式。当目标文件夹不存在时可以自动建立目标文件夹，并且支持递归建立子文件夹并将生成文件保存至对应的子文件夹。

以下是代码：

  
  
function MSTAR2JPG(sourcePath, targetPath)
if ~ exist (targetPath,'dir')
mkdir (targetPath) ;
end
Files = dir(sourcePath);
for i = 1:length(Files)
if Files(i).isdir == 0
FID = fopen([sourcePath '\' Files(i). name], 'rb', 'ieee-be');
while ~feof(FID) % 在PhoenixHeader找到图片尺寸大小
Text = fgetl(FID);
if ~isempty(strfind(Text, 'NumberOfColumns'))
ImgColumns = str2double(Text( 18: end));
Text = fgetl(FID);
ImgRows = str2double(Text( 15: end));
break;
end
end
while ~feof(FID) % 跳过PhoenixHeader
Text = fgetl(FID);
if ~isempty(strfind(Text, '[EndofPhoenixHeader]'))
break
end
end
Mag = fread(FID,ImgColumns*ImgRows, 'float32', 'ieee-be');
Img = reshape(Mag,[ImgColumns ImgRows]);
imwrite(uint8(imadjust(Img)* 255),[targetPath '\' Files(i). name( 1: end- 3) 'JPG']); % 调整对比度后保存
clear ImgColumns ImgRows
fclose(FID);
else
if strcmp(Files(i). name, '.') ~= 1 && strcmp(Files(i). name, '..') ~= 1
if ~exist([targetPath '\' Files(i). name], 'dir')
mkdir([targetPath '\' Files(i). name]);
end
MSTAR2JPG([sourcePath '\' Files(i). name],[targetPath '\' Files(i). name]);
end
end
end
end