what is wrong with my code?

2 views (last 30 days)
yanjie qi
yanjie qi on 28 Feb 2016
Commented: Walter Roberson on 1 Jan 2018
function [lxw,XO,wl]= HSIFileOpen(lxwfilepath,HSIfilepath)
%HSIFileOpen函数用来打开高光谱影像数据
fp1 = fopen(lxwfilepath,'r');
fp2 = fopen(HSIfilepath,'r');
lxw=fread(fp1,'%f');
bands = lxw(1);%波段数
datatype = lxw(2);%字节数
samples = lxw(3);%列数
lines = lxw(4);%行数
columns = samples*lines;%像元个数
switch lxw(5)%数据格式
case 0
interleave = 'bsq';
case 1
interleave = 'bil';
case 2
interleave = 'bip';
case 3
interleave = 'mat';
end
%读取高光谱数据二进制文件
switch datatype
case 1
precision = 'uint8';
case 2
precision = 'uint16';
case 4
precision = 'float32';
end
if interleave == 'bsq'
XO = fread(fp2,[columns,bands],precision);
XO = XO';
elseif interleave == 'bil'
tmp = fread(fp2,[samples,bands*lines],precision);
XO = zeros(bands,columns);
for i=1:bands
for j=1:lines
XO(i,((j-1)*samples+1):j*samples) = tmp(:,(j-1)*bands+i);
end
end
elseif interleave == 'bip'
XO = fread(fp2,[bands,columns],precision);
elseif interleave == 'mat'
XO = fread(fp2,[bands,columns],precision);
end
b = lxw(6);%该参数判断有没有波长数据
if b==1
token = strtok(HSIfilepath,'.');
wavelengthpath = strcat(token,'.wl');
fp3 = fopen(wavelengthpath,'r');
wl=fscanf(fp3,'%f');
%关闭文件
fclose(fp3);
else
wl = 0;
end
fclose(fp1);
fclose(fp2);
[EDITED, Jan, copied from tags:]
but the error is
input argument "lxwfilepath" is undefined.
error in ==> hsifileopen at 6 fp1 = fopen(lxwfilepath">
  1 Comment
Jan
Jan on 28 Feb 2016
Please post the complete error message in your question, not a partial copy in the tags. Thanks.

Sign in to comment.

Accepted Answer

Jan
Jan on 28 Feb 2016
Edited: Jan on 28 Feb 2016
Guessing that the line 6 is this:
fp1 = fopen(lxwfilepath,'r');
I assume, that you call your function without defining input arguments. It has to be called like this:
[lxw, XO, wl] = HSIFileOpen(lxwfilepath, HSIfilepath)
Note: This will work in your case, because interleave has 3 characters in all cases:
if interleave == 'bsq'
But prefer strcmp for the comparison of strings:
if strcmp(interleave, 'bsq')
  3 Comments
yanjie qi
yanjie qi on 1 Mar 2016
sir,
This, I mean the following code. Should I change this code? There is something wrong with it. I change it as you suggest

Sign in to comment.

More Answers (3)

yanjie qi
yanjie qi on 28 Feb 2016
Edited: Walter Roberson on 1 Mar 2016
I run the M file,and the answer is this:
>> [lxw,XO,wl]= HSIFileOpen('F:\多光谱材料\pca\2.hdr','F:\多光谱材料\pca\2.raw')
??? Error using ==> fread
Invalid precision.
Error in ==> HSIFileOpen at 10
lxw=fread(fp1,'%f');
  3 Comments
Walter Roberson
Walter Roberson on 1 Mar 2016
That code does not protect against the possibility of unexpected content in the data file. That code only defines the variable named interleave if a particular location in the file contains 0, 1, 2, or 3. At the very least the code should have an "otherwise" on the switch statement that generates an error saying that the file is not in the expected format.

Sign in to comment.


Walter Roberson
Walter Roberson on 1 Mar 2016
You should be considering using http://www.mathworks.com/matlabcentral/fileexchange/29344-read-medical-data-3d from the File Exchange, as it reads .raw files with .hdr . It is tested code.
  2 Comments
yanjie qi
yanjie qi on 1 Mar 2016
SIR,
??? Error using ==> analyze75info>parseInputs at 457 F:\多光谱材料\pca\2.hdr is not a valid Analyze 7.5 format file.
It looks like something wrong with my .hdr file?

Sign in to comment.


Mehdi Mosafer
Mehdi Mosafer on 1 Jan 2018
I think this issue is because of using non-ASCII alphabets for naming a folder in which the data file is located.
  1 Comment
Walter Roberson
Walter Roberson on 1 Jan 2018
No, the user did not pass the file name in to the function.

Sign in to comment.

Tags

No tags entered yet.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!