]> git.tdb.fi Git - ext/subsurface.git/blobdiff - uemis.c
Integrate loading of uemis SDA files into the regular xml parsing
[ext/subsurface.git] / uemis.c
diff --git a/uemis.c b/uemis.c
index 63541353b0b09ef2ccbc655cc8b55ec8a89cd5f2..769752915edb1f93e645645a8097af56d79b71c8 100644 (file)
--- a/uemis.c
+++ b/uemis.c
 #include <time.h>
 #include <regex.h>
 
-#include <gtk/gtk.h>
-
 #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,65 @@ 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, void *datap) {
        int datalen;
        int i;
        uint8_t *data;
        struct sample *sample;
+       struct dive **divep = datap;
+       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 +265,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,"<val key=\"file_content\">",
-                       ">\\([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 +284,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);
-}
+}