-
Notifications
You must be signed in to change notification settings - Fork 4
/
sequential_read.c
46 lines (35 loc) · 1.3 KB
/
sequential_read.c
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
// an example programme that uses slow5lib to read a SLOW5 file sequentially
#include <stdio.h>
#include <stdlib.h>
#include <slow5/slow5.h>
#define FILE_PATH "examples/example.slow5"
#define TO_PICOAMPS(RAW_VAL,DIGITISATION,OFFSET,RANGE) (((RAW_VAL)+(OFFSET))*((RANGE)/(DIGITISATION)))
int main(){
//open the SLOW5 file
slow5_file_t *sp = slow5_open(FILE_PATH,"r");
if(sp==NULL){
fprintf(stderr,"Error in opening file\n");
exit(EXIT_FAILURE);
}
slow5_rec_t *rec = NULL; //slow5 record to be read
int ret=0; //for return value
//iterate through the file until end
while((ret = slow5_get_next(&rec,sp)) >= 0){
printf("%s\t",rec->read_id);
uint64_t len_raw_signal = rec->len_raw_signal;
for(uint64_t i=0;i<len_raw_signal;i++){ //iterate through the raw signal and print in picoamperes
double pA = TO_PICOAMPS(rec->raw_signal[i],rec->digitisation,rec->offset,rec->range);
printf("%f ",pA);
}
printf("\n");
}
if(ret != SLOW5_ERR_EOF){ //check if proper end of file has been reached
fprintf(stderr,"Error in slow5_get_next. Error code %d\n",ret);
exit(EXIT_FAILURE);
}
//free the SLOW5 record
slow5_rec_free(rec);
//close the SLOW5 file
slow5_close(sp);
return 0;
}