赞
踩
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
Copyright © 2003-2013 www.wpsshop.cn 版权所有,并保留所有权利。