Rのkmeansによる教師なし分類
衛星画像(5バンド)の時系列があり、Rのkmeansでそれらを分類したいと思います。私のスクリプトは正常に動作しています(画像をループして、画像をdata.frameに変換し、クラスター化して、元の画像に戻します。ラスター): for (n in files) { image <- stack(n) image <- clip(image,subset) ###classify raster image.df <- as.data.frame(image) cluster.image <- kmeans(na.omit(image.df), 10, iter.max = 10, nstart = 25) ### kmeans, with 10 clusters #add back NAs using the NAs in band 1 (identic NA positions in all bands), see http://stackoverflow.com/questions/12006366/add-back-nas-after-removing-them/12006502#12006502 image.df.factor …