Я хотел бы реализовать быстрый PLDA (вероятностный линейный дискриминантный анализ) в OpenCV. в этом LINK быстрая PLDA реализована в Matlab
и Python
. Одной из частей PLDA является LDA. Я написал следующий код для реализации LDA в OpenCV:
int LDA_dim = 120;
// Load data
FileStorage fs("newStorageFile.yml", FileStorage::READ);
// Read data
Mat train_data, train_labels;
fs["train_data"] >> train_data;
fs["train_labels"] >> train_labels;
// LDA
if (LDA_dim > 0)
{
LDA lda(LDA_dim);
lda.compute(train_data, train_labels); // compute eigenvectors
Mat eigenvectors = lda.eigenvectors();
}
Я преобразовал базу данных, которая была введена в ссылке выше от .mat
до .yml
. В результате newStorageFile.yml
я загрузил здесь. train_data
имеют 650 строк и 600 столбцов, а train_labels - 650 строк и 1 столбик. Я не знаю, почему собственные векторы и собственное значение становятся равными нулю?? PLZ поможет мне исправить этот код.
Лучше привести код, который преобразует данные с .mat
в .yml
:
function matlab2opencv( variable, fileName, flag)
[rows cols] = size(variable);
% Beware of Matlab linear indexing
variable = variable';
% Write mode as default
if ( ~exist('flag','var') )
flag = 'w';
end
if ( ~exist(fileName,'file') || flag == 'w' )
% New file or write mode specified
file = fopen( fileName, 'w');
fprintf( file, '%%YAML:1.0\n');
else
% Append mode
file = fopen( fileName, 'a');
end
% Write variable header
fprintf( file, ' %s: !!opencv-matrix\n', inputname(1));
fprintf( file, ' rows: %d\n', rows);
fprintf( file, ' cols: %d\n', cols);
fprintf( file, ' dt: f\n');
fprintf( file, ' data: [ ');
% Write variable data
for i=1:rows*cols
fprintf( file, '%.6f', variable(i));
if (i == rows*cols), break, end
fprintf( file, ', ');
if mod(i+1,4) == 0
fprintf( file, '\n ');
end
end
fprintf( file, ']\n');
fclose(file);
Изменить 1) Я попробовал LDA с некоторым образцом, который сам генерировал:
Mat train_data = (Mat_<double>(3, 3) << 25, 45, 44, 403, 607, 494, 2900, 5900, 2200);
Mat train_labels = (Mat_<int>(3, 1) << 1, 2, 3 );
LDA lda(LDA_dim);
lda.compute(train_data, train_labels); // compute eigenvectors
Mat_<double> eigenvectors = lda.eigenvectors();
Mat_<double> eigenvalues = lda.eigenvalues();
cout << eigenvectors << endl << eigenvalues;
но я должен получить тот же результат: собственное значение и собственный вектор становятся равными нулю: