/*!
* \brief MainWindow::do_EPI_peaks
*/
void MainWindow::do_EPI_peaks(){
if (DEBUG) cout << Q_FUNC_INFO << endl;
if (EPI_coord_ifft->size() != 0){
//Peaks detection
detect_peaks(*EPI_coord_ifft, EPI_coord_peaks, EPI_threshold);
//Plot
plot_EPI_peaks();
//Status
status_info("");
}
else{
status_error("");
}
}
int main(int argc, char * argv[]){
FILE * signal, * peaks;
int i, j;
time_series * ankle_joint;
motion m;
unidimentional_series ankle, ankle_peaks;
bvh_load_data( argv[1], &m );
unidimentional_series_init( &ankle, 1000);
unidimentional_series_init( &ankle_peaks, 100);
ankle_joint = m.data.get( m.data_ptr, "RightFoot");
get_component( ankle_joint, &ankle, 'y');
detect_peaks( &ankle, &ankle_peaks, 5);
signal = fopen("signal.txt", "w");
peaks = fopen("peaks.txt", "w");
for(i = 0; i < ankle_joint -> length; ++i)
fprintf(signal, "%d %lf\n", i, ankle_joint -> begin[i].y);
for(i = 0; i < ankle_peaks.length; ++i){
j = ankle_peaks.begin[i];
fprintf(peaks, "%d %lf\n", j, ankle_joint -> begin[j].y);
}
fclose(signal);
fclose(peaks);
bvh_unload_data( &m );
ankle.destroy( &ankle );
ankle_peaks.destroy( &ankle_peaks );
return 0;
}
