X-Git-Url: http://git.tdb.fi/?a=blobdiff_plain;f=uemis.c;h=259d30c53758f9d8ad3e8f2868ef47fa30534efa;hb=f2566ba561ff8eb560e71dd89de5e66e46e7143d;hp=63541353b0b09ef2ccbc655cc8b55ec8a89cd5f2;hpb=518ec33ba3bf9c5e9966ef108c1673223337b0fa;p=ext%2Fsubsurface.git diff --git a/uemis.c b/uemis.c index 6354135..259d30c 100644 --- a/uemis.c +++ b/uemis.c @@ -16,10 +16,7 @@ #include #include -#include - #include "dive.h" -#include "display.h" #include "uemis.h" /* @@ -134,6 +131,24 @@ int matchit(FILE *infd, char *regex, char *typeregex, char **found) { return 0; } +/* + * pressure_to_depth: In centibar. And when converting to + * depth, I'm just going to always use saltwater, because I + * think "true depth" is just stupid. From a diving standpoint, + * "true depth" is pretty much completely pointless, unless + * you're doing some kind of underwater surveying work. + * + * So I give water depths in "pressure depth", always assuming + * salt water. So one atmosphere per 10m. + */ +static int pressure_to_depth(uint16_t value) +{ + double atm, cm; + + atm = (value / 100.0) / 1.01325; + cm = 100 * atm + 0.5; + return( (cm > 0) ? 10 * (long)cm : 0); +} /* * convert the base64 data blog @@ -155,7 +170,7 @@ int uemis_convert_base64(char *base64, uint8_t **data) { } decode(base64, *data, len); - if (memcmp(data,"Dive\01\00\00",7)) + if (memcmp(*data,"Dive\01\00\00",7)) fprintf(stderr,"Missing Dive100 header\n"); bail: @@ -165,32 +180,64 @@ bail: /* * parse uemis base64 data blob into struct dive */ -static void parse_divelog_binary(char *base64, struct dive **divep) { +void uemis_parse_divelog_binary(char *base64, struct dive **divep) { int datalen; int i; uint8_t *data; struct sample *sample; + struct dive *dive = *divep; + int template, gasoffset; datalen = uemis_convert_base64(base64, &data); + /* dive template in use: + 0 = air + 1 = nitrox (B) + 2 = nitrox (B+D) + 3 = nitrox (B+T+D) + uemis cylinder data is insane - it stores seven tank settings in a block + and the template tells us which of the four groups of tanks we need to look at + */ + gasoffset = template = *(uint8_t *)(data+115); + if (template == 3) + gasoffset = 4; + for (i = 0; i < template; i++) { + float volume = *(float *)(data+116+25*(gasoffset + i)) * 1000.0; + /* uemis always assumes a working pressure of 3000psi / 206bar - even though that's wrong */ + /* I also think that the unit that it stores (cuft for me) might change with SDA settings */ + // dive->cylinder[i].type.size.mliter = volume * 206.84 / 28.317; + dive->cylinder[i].type.size.mliter = volume * 200 / 28.317; + dive->cylinder[i].gasmix.o2.permille = *(uint8_t *)(data+120+25*(gasoffset + i)) * 10 + 0.5; + dive->cylinder[i].gasmix.he.permille = 0; + } /* first byte of divelog data is at offset 0x123 */ i = 0x123; while ((i < datalen) && (*(uint16_t *)(data+i))) { /* it seems that a dive_time of 0 indicates the end of the valid readings */ + /* the SDA usually records more samples after the end of the dive -- + * we want to discard those, but not cut the dive short; sadly the dive + * duration in the header is a) in minutes and b) up to 3 minutes short */ + if (*(uint16_t *)(data+i) > dive->duration.seconds + 180) + break; sample = prepare_sample(divep); + dive = *divep; /* prepare_sample might realloc the dive */ sample->time.seconds = *(uint16_t *)(data+i); - sample->depth.mm = (*(uint16_t *)(data+i+2) - 100) / 0.101428571 + 0.5; + sample->depth.mm = pressure_to_depth(*(uint16_t *)(data+i+2)); sample->temperature.mkelvin = (*(uint16_t *)(data+i+4) * 100) + 273150; sample->cylinderpressure.mbar= *(uint16_t *)(data+i+23) * 10; sample->cylinderindex = *(uint8_t *)(data+i+22); - finish_sample(*divep, sample); + finish_sample(dive, sample); i += 0x25; } + dive->duration.seconds = sample->time.seconds - 1; return; } +/* parse a single file + * TODO: we don't report any errors when the parse fails - we simply don't add them to the list + */ void -parse_uemis_file(char *divelogfilename,GError **error) { +parse_uemis_file(char *divelogfilename) { char *found=NULL; struct tm tm; struct dive *dive; @@ -217,14 +264,14 @@ parse_uemis_file(char *divelogfilename,GError **error) { /* some error handling */ goto bail; } - dive->maxdepth.mm = atof(found) / 0.10143 + 0.5; + dive->maxdepth.mm = pressure_to_depth(atoi(found)); if (! matchit(divelogfile,"", - ">\\([a-zA-Z0-9+/]*\\)<", &found)) { + ">\\([a-zA-Z0-9+/=]*\\)<", &found)) { /* some error handling */ goto bail; } - parse_divelog_binary(found,&dive); + uemis_parse_divelog_binary(found,&dive); record_dive(dive); bail: if (found) @@ -236,40 +283,7 @@ bail: */ void uemis_import() { - GtkWidget *dialog; - GtkFileFilter *filter = gtk_file_filter_new (); - gtk_file_filter_add_pattern (filter, "*.SDA"); - gtk_file_filter_set_name(filter, "uemis Zurich SDA files"); - dialog = gtk_file_chooser_dialog_new("Open File", - GTK_WINDOW(main_window), - GTK_FILE_CHOOSER_ACTION_OPEN, - GTK_STOCK_CANCEL, GTK_RESPONSE_CANCEL, - GTK_STOCK_OPEN, GTK_RESPONSE_ACCEPT, - NULL); - gtk_file_chooser_set_select_multiple(GTK_FILE_CHOOSER(dialog), TRUE); - gtk_file_chooser_add_filter(GTK_FILE_CHOOSER(dialog),filter); - - if (gtk_dialog_run(GTK_DIALOG(dialog)) == GTK_RESPONSE_ACCEPT) { - GSList *filenames; - char *filename; - filenames = gtk_file_chooser_get_filenames(GTK_FILE_CHOOSER(dialog)); - - GError *error = NULL; - while(filenames != NULL) { - filename = (char *)filenames->data; - parse_uemis_file(filename, &error); - if (error != NULL) - { - report_error(error); - g_error_free(error); - error = NULL; - } - - g_free(filename); - filenames = g_slist_next(filenames); - } - g_slist_free(filenames); + if (open_import_file_dialog("*.SDA","uemis Zurich SDA files", + &parse_uemis_file)) report_dives(); - } - gtk_widget_destroy(dialog); -} +}