Question

matlab 三维矩阵降维

对一个高光谱图像做波段聚类，用kmeans函数做。首先我想把它降维，请问该怎么做。
比如这个图像时500x300x200的，其中200是波段数
本来我的想法是这样：
for m = 1:204
DR[m]=[salinas_corrected(:,:,m)];
end
但是报错说Error: Unbalanced or unexpected parenthesis or bracket.
而如果用DR=[salinas_corrected(:,:)];
DR就变成了一个500X60000的矩阵了，波段貌似就混了……
所以请问到底该怎么降维呢……

Answer 1

如果每次只取一个波段进行计算的话，直接

for m = 1:200
    DR=[salinas_corrected(:,:,m)];
end

就可以了

每次使用DR进行计算，而DR的维数就是500×300。

追问

我是要对波段进行聚类，所以算聚类的时候用的是全部的波段0.0

追答

就算聚类，也是每次取一个波段的数据，然后看这个波段的数据属于哪一类吧？？
聚类的具体实践我不是很清楚。

追问

如果我自己编函数的话，是你说的那样，但是我想用matlab自带的kmeans函数来算，然后他输入参数是个二维矩阵，所以我想降维

追答

那如果这样的话，按照你的方法，把矩阵变成500X60000，也没有问题啊，反正函数会把相似的点放到一块。

追问

我的意思是我希望它变成200x150000这样，因为200才是波段数

追答

哦哦，这样的话，

那就好办了。

for m = 1:200
    DR(m,:)=reshape(salinas_corrected(:,:,m),1,[]);
end

DR就是200x150000的了。