/*-------------------------------------------------------------------- * The MB-system: mbnavedit_prog.c 6/23/95 * $Id: mbnavedit_prog.c 1829 2010-02-05 02:53:39Z caress $ * * Copyright (c) 1995-2009 by * David W. Caress (caress@mbari.org) * Monterey Bay Aquarium Research Institute * Moss Landing, CA 95039 * and Dale N. Chayes (dale@ldeo.columbia.edu) * Lamont-Doherty Earth Observatory * Palisades, NY 10964 * * See README file for copying and redistribution conditions. *--------------------------------------------------------------------*/ /* * MBNAVEDIT is an interactive navigation editor for swath sonar data. * It can work with any data format supported by the MBIO library. * This file contains the code that does not directly depend on the * MOTIF interface - the companion files mbnavedit.c and * mbnavedit_callbacks.c contain the user interface related code. * * Author: D. W. Caress * Date: June 23, 1995 * Date: August 28, 2000 (New version - no buffered i/o) * * $Log: mbnavedit_prog.c,v $ * Revision 5.24 2009/03/10 05:11:22 caress * Added Gaussian mean smoothing to MBnavedit. * * Revision 5.23 2009/03/09 16:58:31 caress * Release 5.1.2beta01 * * Revision 5.22 2008/07/20 15:32:14 caress * Lengthened filename char arrays to prevent possible overflows. * * Revision 5.21 2008/05/16 23:05:05 caress * Release 5.1.1beta18. * * Revision 5.20 2007/10/08 16:17:00 caress * Repeated positions are now automatically flagged for modeling or smooth inversion. * * Revision 5.19 2006/02/16 21:15:07 caress * Made smooth inversion weights work as small as 0.01 in the interface. Redimensioned some strings too. * * Revision 5.18 2006/01/24 19:19:24 caress * Version 5.0.8 beta. * * Revision 5.17 2006/01/06 18:24:13 caress * Working towards 5.0.8 * * Revision 5.16 2005/11/05 00:58:10 caress * Programs changed to register arrays through mb_register_array() rather than allocating the memory directly with mb_realloc() or mb_malloc(). * * Revision 5.15 2005/06/04 04:45:50 caress * Added feature to apply longitude and latitude offsets to the navigation. * * Revision 5.14 2005/03/25 04:35:55 caress * Added capability to interpolate over repeated values. * * Revision 5.13 2004/12/18 01:36:40 caress * Working towards release 5.0.6. * * Revision 5.12 2004/05/21 23:33:03 caress * Moved to new version of BX GUI builder * * Revision 5.11 2003/07/02 18:13:13 caress * Release 5.0.0 * * Revision 5.10 2003/04/17 21:09:06 caress * Release 5.0.beta30 * * Revision 5.9 2002/10/02 23:54:07 caress * Release 5.0.beta24 * * Revision 5.8 2001/08/02 01:51:39 caress * Fixed call to mb_pr_ function. * * Revision 5.7 2001/07/20 00:34:00 caress * Release 5.0.beta03 * * Revision 5.6 2001/06/30 17:41:04 caress * Release 5.0.beta01 * * Revision 5.5 2001/06/03 07:06:25 caress * Release 5.0.beta01 * * Revision 5.4 2001/04/06 22:16:01 caress * Fixed unflagging. * * Revision 5.3 2001/03/22 21:10:37 caress * Trying to make release 5.0.beta0. * * Revision 5.2 2001/01/22 07:47:40 caress * Version 5.0.beta01 * * Revision 5.1 2000/12/10 20:30:08 caress * Version 5.0.alpha02 * * Revision 5.0 2000/12/01 22:56:08 caress * First cut at Version 5.0. * * Revision 4.22 2000/10/11 01:05:17 caress * Convert to ANSI C * * Revision 4.21 2000/09/30 07:03:14 caress * Snapshot for Dale. * * Revision 4.20 2000/09/30 07:02:34 caress * Snapshot for Dale. * * Revision 4.19 2000/08/28 22:45:11 caress * About to kick off new version. * * Revision 4.18 1999/12/11 04:42:03 caress * Moved xgraphics.c to src/xgraphics * * Revision 4.17 1999/11/12 21:25:21 caress * Fixed start date shown below plots * * Revision 4.16 1999/07/16 19:22:56 caress * Smaller window with new dialogs for Linux. * * Revision 4.15 1999/04/14 04:33:10 caress * Final (?) version 4.6 release * * Revision 4.14 1999/04/09 22:34:08 caress * Added time interval plot. * * Revision 4.13 1999/02/04 23:54:13 caress * MB-System version 4.6beta7 * * Revision 4.12 1998/12/18 01:40:25 caress * MB-System version 4.6beta5 * * Revision 4.11 1998/10/09 18:43:08 caress * MB-System version 4.6beta. * * Revision 4.10 1997/09/15 19:10:20 caress * Real Version 4.5 * * Revision 4.9 1997/04/22 19:25:57 caress * Fixed startup mode. * * Revision 4.8 1997/04/21 17:07:38 caress * MB-System 4.5 Beta Release. * * Revision 4.8 1997/04/16 21:40:29 caress * Version for MB-System 4.5. * * Revision 4.7 1996/08/26 17:33:15 caress * Release 4.4 revision. * * Revision 4.6 1996/04/22 13:22:24 caress * Now have DTR and MIN/MAX defines in mb_define.h * * Revision 4.5 1996/04/05 20:07:02 caress * Added GUI mode so done means quit for real. Also changed done and * quit handling in browse mode so that the program doesn't read the * entire data file before closing it. * * Revision 4.4 1995/11/02 19:22:45 caress * Fixed mb_error calls. * * Revision 4.3 1995/09/28 18:01:01 caress * Improved handling of .mbxxx file suffix convention. * * Revision 4.2 1995/09/18 22:40:44 caress * Fixed bug that caused "select all" function to miss last data * point on plots of entire data set. * * Revision 4.1 1995/08/17 14:58:12 caress * Revision for release 4.3. * * Revision 4.0 1995/08/07 18:33:22 caress * First cut. * * */ /*--------------------------------------------------------------------*/ /* standard include files */ #include #include #include #include #include #include #include #include #include /* MBIO include files */ #include "../../include/mb_format.h" #include "../../include/mb_status.h" #include "../../include/mb_define.h" #include "../../include/mb_io.h" #include "../../include/mb_process.h" #include "../../include/mb_aux.h" #include "../../include/mb_xgraphics.h" /* define global control parameters */ #include "mbnavedit.h" /* ping structure definition */ struct mbnavedit_ping_struct { int id; int record; int time_i[7]; double time_d; double file_time_d; double tint; double lon; double lat; double speed; double heading; double draft; double roll; double pitch; double heave; double time_d_org; double tint_org; double lon_org; double lat_org; int mean_ok; double lon_dr; double lat_dr; double speed_org; double heading_org; double draft_org; double speed_made_good; double course_made_good; int tint_x; int tint_y; int lon_x; int lon_y; int lat_x; int lat_y; int speed_x; int speed_y; int heading_x; int heading_y; int draft_x; int draft_y; int tint_select; int lon_select; int lat_select; int speed_select; int heading_select; int draft_select; int lonlat_flag; }; /* plot structure definition */ struct mbnavedit_plot_struct { int type; int ixmin; int ixmax; int iymin; int iymax; double xmin; double xmax; double ymin; double ymax; double xscale; double yscale; double xinterval; double yinterval; char xlabel[128]; char ylabel1[128]; char ylabel2[128]; }; /* id variables */ static char rcs_id[] = "$Id: mbnavedit_prog.c 1829 2010-02-05 02:53:39Z caress $"; static char program_name[] = "MBNAVEDIT"; static char help_message[] = "MBNAVEDIT is an interactive navigation editor for swath sonar data.\n\tIt can work with any data format supported by the MBIO library.\n"; static char usage_message[] = "mbnavedit [-Byr/mo/da/hr/mn/sc -D -Eyr/mo/da/hr/mn/sc \n\t-Fformat -Ifile -Ooutfile -X -V -H]"; /* status variables */ static int error = MB_ERROR_NO_ERROR; static int verbose = 0; static char *message = NULL; /* MBIO control parameters */ static int nav_source; static int heading_source; static int vru_source; static int svp_source; static int pings; static int lonflip; static double bounds[4]; static int btime_i[7]; static int etime_i[7]; static double btime_d; static double etime_d; static double speedmin; static double timegap; static int beams_bath; static int beams_amp; static int pixels_ss; static void *imbio_ptr = NULL; /* mbio read and write values */ static void *store_ptr = NULL; static int kind; static double distance; static double altitude; static double sonardepth; static int nbath; static int namp; static int nss; static char *beamflag = NULL; static double *bath = NULL; static double *bathacrosstrack = NULL; static double *bathalongtrack = NULL; static double *amp = NULL; static double *ss = NULL; static double *ssacrosstrack = NULL; static double *ssalongtrack = NULL; static char comment[MB_COMMENT_MAXLINE]; /* buffer control variables */ #define MBNAVEDIT_BUFFER_SIZE 1000000 static int file_open = MB_NO; static int nfile_open = MB_NO; static FILE *nfp; static int hold_size = 100; static int nload = 0; static int ndump = 0; static int nbuff = 0; static int current_id = 0; static int nload_total = 0; static int ndump_total = 0; static int first_read = MB_NO; /* plotting control variables */ #define NUMBER_PLOTS_MAX 9 #define DEFAULT_PLOT_WIDTH 767 #define DEFAULT_PLOT_HEIGHT 300 #define MBNAVEDIT_PICK_DISTANCE 50 #define MBNAVEDIT_ERASE_DISTANCE 10 static struct mbnavedit_ping_struct ping[MBNAVEDIT_BUFFER_SIZE]; static double plot_start_time; static double plot_end_time; static int nplot; static void *mbnavedit_xgid; static struct mbnavedit_plot_struct mbnavplot[NUMBER_PLOTS_MAX]; static int data_save; static double file_start_time_d; /* color control values */ #define WHITE 0 #define BLACK 1 #define RED 2 #define GREEN 3 #define BLUE 4 #define ORANGE 5 #define PURPLE 6 #define CORAL 7 #define LIGHTGREY 8 #define XG_SOLIDLINE 0 #define XG_DASHLINE 1 static int ncolors; static int pixel_values[256]; /*--------------------------------------------------------------------*/ int mbnavedit_init_globals() { /* local variables */ char *function_name = "mbnavedit_init_globals"; int status = MB_SUCCESS; /* set default global control parameters */ output_mode = OUTPUT_MODE_OUTPUT; run_mbprocess = MB_NO; gui_mode = MB_NO; data_show_max = 2000; data_show_size = 1000; data_step_max = 2000; data_step_size = 750; mode_pick = PICK_MODE_PICK; mode_set_interval = MB_NO; plot_tint = MB_YES; plot_tint_org = MB_YES; plot_lon = MB_YES; plot_lon_org = MB_YES; plot_lon_dr = MB_NO; plot_lat = MB_YES; plot_lat_org = MB_YES; plot_lat_dr = MB_NO; plot_speed = MB_YES; plot_speed_org = MB_YES; plot_smg = MB_YES; plot_heading = MB_YES; plot_heading_org = MB_YES; plot_cmg = MB_YES; plot_draft = MB_YES; plot_draft_org = MB_YES; plot_draft_dr = MB_NO; plot_roll = MB_NO; plot_pitch = MB_NO; plot_heave = MB_NO; mean_time_window = 100; drift_lon = 0; drift_lat = 0; strcpy(ifile,"\0"); plot_width = DEFAULT_PLOT_WIDTH; plot_height = DEFAULT_PLOT_HEIGHT; number_plots = 0; if (plot_tint == MB_YES) number_plots++; if (plot_lon == MB_YES) number_plots++; if (plot_lat == MB_YES) number_plots++; if (plot_speed == MB_YES) number_plots++; if (plot_heading == MB_YES) number_plots++; if (plot_draft == MB_YES) number_plots++; timestamp_problem = MB_NO; use_ping_data = MB_NO; strip_comments = MB_NO; model_mode = MODEL_MODE_OFF; weight_speed = 100.0; weight_acceleration = 100.0; scrollcount = 0; offset_lon = 0.0; offset_lat = 0.0; offset_lon_applied = 0.0; offset_lat_applied = 0.0; /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_init(int argc, char **argv) { /* local variables */ char *function_name = "mbnavedit_init"; int status = MB_SUCCESS; int fileflag = 0; struct stat file_status; int fstat; int i; /* parsing variables */ extern char *optarg; int errflg = 0; int c; int help = 0; int flag = 0; /* set default values */ status = mb_defaults(verbose,&format,&pings,&lonflip,bounds, btime_i,etime_i,&speedmin,&timegap); pings = 1; lonflip = 0; bounds[0] = -360.; bounds[1] = 360.; bounds[2] = -90.; bounds[3] = 90.; btime_i[0] = 1962; btime_i[1] = 2; btime_i[2] = 21; btime_i[3] = 10; btime_i[4] = 30; btime_i[5] = 0; btime_i[6] = 0; etime_i[0] = 2062; etime_i[1] = 2; etime_i[2] = 21; etime_i[3] = 10; etime_i[4] = 30; etime_i[5] = 0; etime_i[6] = 0; speedmin = 0.0; timegap = 1000000000.0; strcpy(ifile,"\0"); /* process argument list */ while ((c = getopt(argc, argv, "VvHhB:b:DdE:e:F:f:GgI:i:NnPpXx")) != -1) switch (c) { case 'H': case 'h': help++; break; case 'V': case 'v': verbose++; break; case 'B': case 'b': sscanf (optarg,"%d/%d/%d/%d/%d/%d", &btime_i[0],&btime_i[1],&btime_i[2], &btime_i[3],&btime_i[4],&btime_i[5]); btime_i[6] = 0; flag++; break; case 'D': case 'd': output_mode = OUTPUT_MODE_BROWSE; flag++; break; case 'E': case 'e': sscanf (optarg,"%d/%d/%d/%d/%d/%d", &etime_i[0],&etime_i[1],&etime_i[2], &etime_i[3],&etime_i[4],&etime_i[5]); etime_i[6] = 0; flag++; break; case 'F': case 'f': sscanf (optarg,"%d", &format); flag++; break; case 'G': case 'g': gui_mode = MB_YES; flag++; break; case 'I': case 'i': sscanf (optarg,"%s", ifile); flag++; fileflag++; break; case 'N': case 'n': strip_comments = MB_YES; flag++; break; case 'P': case 'p': use_ping_data = MB_YES; flag++; break; case 'X': case 'x': run_mbprocess = MB_YES; flag++; break; case '?': errflg++; } /* if error flagged then print it and exit */ if (errflg) { fprintf(stderr,"usage: %s\n", usage_message); fprintf(stderr,"\nProgram <%s> Terminated\n", program_name); error = MB_ERROR_BAD_USAGE; exit(error); } /* print starting message */ if (verbose == 1 || help) { fprintf(stderr,"\nProgram %s\n",program_name); fprintf(stderr,"Version %s\n",rcs_id); fprintf(stderr,"MB-system Version %s\n",MB_VERSION); } /* print starting debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 Program <%s>\n",program_name); fprintf(stderr,"dbg2 Version %s\n",rcs_id); fprintf(stderr,"dbg2 MB-system Version %s\n",MB_VERSION); fprintf(stderr,"dbg2 Control Parameters:\n"); fprintf(stderr,"dbg2 verbose: %d\n",verbose); fprintf(stderr,"dbg2 help: %d\n",help); fprintf(stderr,"dbg2 format: %d\n",format); fprintf(stderr,"dbg2 input file: %s\n",ifile); } /* if help desired then print it and exit */ if (help) { fprintf(stderr,"\n%s\n",help_message); fprintf(stderr,"\nusage: %s\n", usage_message); exit(error); } /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 argc: %d\n",argc); for (i=0;i 0) { /* get the output filename */ strcpy(nfile,ifile); strcat(nfile,".nve"); /* check if output file exists */ fstat = stat(nfile, &file_status); if (fstat != 0 || (file_status.st_mode & S_IFMT) == S_IFDIR) { /* open file directly */ do_open_file(MB_NO); } else { /* check if previous edits to be used - then open the file */ do_checkuseprevious(); } } /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_set_graphics(void *xgid, int ncol, unsigned int *pixels) { /* local variables */ char *function_name = "mbnavedit_set_graphics"; int status = MB_SUCCESS; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 xgid: %lu\n",(size_t)xgid); fprintf(stderr,"dbg2 ncolors: %d\n",ncol); for (i=0;i= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_open(int useprevious) { /* local variables */ char *function_name = "mbnavedit_action_open"; int status = MB_SUCCESS; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); } /* clear the screen */ status = mbnavedit_clear_screen(); /* open the file */ status = mbnavedit_open_file(useprevious); /* load the buffer */ if (status == MB_SUCCESS) status = mbnavedit_load_data(); /* set up plotting */ if (nbuff > 0) { /* set time span to zero so plotting resets it */ data_show_size = 0; /* turn file button off */ do_filebutton_off(); /* now plot it */ status = mbnavedit_plot_all(); } /* if no data read show error dialog */ else do_error_dialog("No data were read from the input", "file. You may have specified an", "incorrect MB-System format id!"); /* reset data_save */ data_save = MB_NO; /* print debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 File open attempted in MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Buffer values:\n"); fprintf(stderr,"dbg2 nload: %d\n",ndump); fprintf(stderr,"dbg2 nload: %d\n",nload); fprintf(stderr,"dbg2 nbuff: %d\n",nbuff); fprintf(stderr,"dbg2 current_id: %d\n",current_id); fprintf(stderr,"dbg2 error: %d\n",error); } /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_open_file(int useprevious) { /* local variables */ char *function_name = "mbnavedit_open_file"; int status = MB_SUCCESS; char ifile_use[MB_PATH_MAXLINE]; char command[MB_PATH_MAXLINE]; int format_use; int form; int format_error; struct stat file_status; int fstat; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 file: %s\n",ifile); fprintf(stderr,"dbg2 format: %d\n",format); fprintf(stderr,"dbg2 useprevious: %d\n",useprevious); } /* get format if required */ if (format == 0) { if (mb_get_format(verbose, ifile, NULL, &form, &format_error) == MB_SUCCESS) { format = form; } } /* if output on and using previously edited nav first copy old nav and then read it as input instead of specified input file */ if (useprevious == MB_YES && output_mode != OUTPUT_MODE_BROWSE) { /* get temporary file name */ sprintf(ifile_use, "%s.tmp", nfile); /* copy old edit save file to tmp file */ sprintf(command, "cp %s %s\n", nfile, ifile_use); format_use = MBF_MBPRONAV; system(command); fstat = stat(ifile_use, &file_status); if (fstat != 0 || (file_status.st_mode & S_IFMT) == S_IFDIR) { do_error_dialog("Unable to copy previously edited", "navigation. You may not have read", "permission in this directory!"); status = MB_FAILURE; return(status); } } /* if output off and using previously edited nav reset input names */ else if (useprevious == MB_YES) { sprintf(ifile_use, "%s", nfile); format_use = MBF_MBPRONAV; } /* else just read from previously edited nav */ else { strcpy(ifile_use, ifile); format_use = format; } /* initialize reading the input multibeam file */ status = mb_format_source(verbose, &format_use, &nav_source, &heading_source, &vru_source, &svp_source, &error); if ((status = mb_read_init( verbose,ifile_use,format_use,pings,lonflip,bounds, btime_i,etime_i,speedmin,timegap, &imbio_ptr,&btime_d,&etime_d, &beams_bath,&beams_amp,&pixels_ss,&error)) != MB_SUCCESS) { mb_error(verbose,error,&message); fprintf(stderr,"\nMBIO Error returned from function :\n%s\n",message); fprintf(stderr,"\nMultibeam File <%s> not initialized for reading\n",ifile); status = MB_FAILURE; do_error_dialog("Unable to open input file.", "You may not have read", "permission in this directory!"); return(status); } /* allocate memory for data arrays */ beamflag = NULL; bath = NULL; amp = NULL; bathacrosstrack = NULL; bathalongtrack = NULL; ss = NULL; ssacrosstrack = NULL; ssalongtrack = NULL; if (error == MB_ERROR_NO_ERROR) status = mb_register_array(verbose, imbio_ptr, MB_MEM_TYPE_BATHYMETRY, sizeof(char), (void **)&beamflag, &error); if (error == MB_ERROR_NO_ERROR) status = mb_register_array(verbose, imbio_ptr, MB_MEM_TYPE_BATHYMETRY, sizeof(double), (void **)&bath, &error); if (error == MB_ERROR_NO_ERROR) status = mb_register_array(verbose, imbio_ptr, MB_MEM_TYPE_AMPLITUDE, sizeof(double), (void **)&, &error); if (error == MB_ERROR_NO_ERROR) status = mb_register_array(verbose, imbio_ptr, MB_MEM_TYPE_BATHYMETRY, sizeof(double), (void **)&bathacrosstrack, &error); if (error == MB_ERROR_NO_ERROR) status = mb_register_array(verbose, imbio_ptr, MB_MEM_TYPE_BATHYMETRY, sizeof(double), (void **)&bathalongtrack, &error); if (error == MB_ERROR_NO_ERROR) status = mb_register_array(verbose, imbio_ptr, MB_MEM_TYPE_SIDESCAN, sizeof(double), (void **)&ss, &error); if (error == MB_ERROR_NO_ERROR) status = mb_register_array(verbose, imbio_ptr, MB_MEM_TYPE_SIDESCAN, sizeof(double), (void **)&ssacrosstrack, &error); if (error == MB_ERROR_NO_ERROR) status = mb_register_array(verbose, imbio_ptr, MB_MEM_TYPE_SIDESCAN, sizeof(double), (void **)&ssalongtrack, &error); /* if error initializing memory then quit */ if (error != MB_ERROR_NO_ERROR) { mb_error(verbose,error,&message); fprintf(stderr,"\nMBIO Error allocating data arrays:\n%s\n",message); fprintf(stderr,"\nProgram <%s> Terminated\n", program_name); exit(error); } /* initialize the buffer */ nbuff = 0; first_read = MB_NO; /* reset plotting time span */ plot_start_time = 0.0; plot_end_time = data_show_size; /* now deal with new nav save file */ nfile_open = MB_NO; if (status == MB_SUCCESS && output_mode != OUTPUT_MODE_BROWSE) { /* get nav edit save file */ sprintf(nfile, "%s.nve", ifile); /* open the nav edit save file */ if ((nfp = fopen(nfile,"w")) != NULL) { nfile_open = MB_YES; } else { nfile_open = MB_NO; fprintf(stderr, "\nUnable to open new nav save file %s\n", nfile); do_error_dialog("Unable to open new nav edit save file.", "You may not have write", "permission in this directory!"); } } /* if we got here we must have succeeded */ if (verbose >= 1) { if (useprevious == MB_YES) { fprintf(stderr,"\nSwath data file <%s> specified for input\n",ifile); fprintf(stderr,"MB-System Data Format ID: %d\n",format); fprintf(stderr,"Navigation data file <%s> initialized for reading\n",ifile_use); fprintf(stderr,"MB-System Data Format ID: %d\n",format_use); } else { fprintf(stderr,"\nSwath data file <%s> initialized for reading\n",ifile_use); fprintf(stderr,"MB-System Data Format ID: %d\n",format_use); } if (output_mode == OUTPUT_MODE_OUTPUT) fprintf(stderr,"Navigation File <%s> initialized for writing\n",nfile); } file_open = MB_YES; /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_close_file() { /* local variables */ char *function_name = "mbnavedit_close_file"; int status = MB_SUCCESS; char command[MB_PATH_MAXLINE]; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); } /* close the files */ status = mb_close(verbose,&imbio_ptr,&error); if (nfile_open == MB_YES) { /* close navigation file */ fclose(nfp); nfile_open = MB_NO; /* update mbprocess parameter file */ status = mb_pr_update_format(verbose, ifile, MB_YES, format, &error); status = mb_pr_update_nav(verbose, ifile, MBP_NAV_ON, nfile, 9, MBP_NAV_ON, MBP_NAV_ON, MBP_NAV_ON, MBP_NAV_ON, MBP_NAV_LINEAR, (double) 0.0, &error); /* run mbprocess if desired */ if (run_mbprocess == MB_YES) { /* turn message on */ do_message_on("Navigation edits being applied using mbprocess..."); /* run mbprocess */ if (strip_comments == MB_YES) sprintf(command, "mbprocess -I %s -N\n",ifile); else sprintf(command, "mbprocess -I %s\n",ifile); system(command); /* turn message off */ do_message_off(); } } /* check memory */ if (verbose >= 4) status = mb_memory_list(verbose,&error); /* if we got here we must have succeeded */ if (verbose >= 1) { fprintf(stderr,"\nMultibeam Input File <%s> closed\n",ifile); if (output_mode == OUTPUT_MODE_OUTPUT) fprintf(stderr,"Navigation Output File <%s> closed\n",nfile); fprintf(stderr,"%d data records loaded\n",nload_total); fprintf(stderr,"%d data records dumped\n",ndump_total); } file_open = MB_NO; nload_total = 0; ndump_total = 0; /* reset offsets */ offset_lon = 0.0; offset_lat = 0.0; offset_lon_applied = offset_lon; offset_lat_applied = offset_lat; /* turn file button on */ do_filebutton_on(); /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_dump_data(int hold) { /* local variables */ char *function_name = "mbnavedit_dump_data"; int status = MB_SUCCESS; int iping; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 hold: %d\n",hold); } /* write out edited data */ if (nfile_open == MB_YES) { for (iping=0;iping 0) { /* turn message on */ do_message_on("MBnavedit is clearing data..."); /* copy data to be held */ for (iping=0;iping 0) current_id = current_id - ndump; if (current_id < 0) current_id = 0; if (current_id > nbuff - 1) current_id = nbuff - 1; /* print out information */ if (verbose >= 1) { if (output_mode == OUTPUT_MODE_OUTPUT) fprintf(stderr,"\n%d data records dumped to output file <%s>\n", ndump,nfile); else fprintf(stderr,"\n%d data records dumped from buffer\n", ndump); fprintf(stderr,"%d data records remain in buffer\n",nbuff); } /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return (status); } /*--------------------------------------------------------------------*/ int mbnavedit_load_data() { /* local variables */ char *function_name = "mbnavedit_load_data"; int status = MB_SUCCESS; int i; char string[MB_PATH_MAXLINE]; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); } /* turn message on */ nload = 0; timestamp_problem = MB_NO; sprintf(string, "MBnavedit: %d records loaded so far...", nload); do_message_on(string); /* load data */ if (status == MB_SUCCESS) do { status = mb_get_all(verbose,imbio_ptr,&store_ptr,&kind, ping[nbuff].time_i, &ping[nbuff].time_d, &ping[nbuff].lon, &ping[nbuff].lat, &ping[nbuff].speed, &ping[nbuff].heading, &distance,&altitude,&sonardepth, &nbath,&namp,&nss, beamflag,bath,amp, bathacrosstrack,bathalongtrack, ss,ssacrosstrack,ssalongtrack, comment,&error); if (error <= MB_ERROR_NO_ERROR && (kind == nav_source || (kind == MB_DATA_DATA && use_ping_data == MB_YES)) && (error == MB_ERROR_NO_ERROR || error == MB_ERROR_TIME_GAP || error == MB_ERROR_OUT_BOUNDS || error == MB_ERROR_OUT_TIME || error == MB_ERROR_SPEED_TOO_SMALL)) { status = MB_SUCCESS; error = MB_ERROR_NO_ERROR; } else if (error <= MB_ERROR_NO_ERROR) { status = MB_FAILURE; error = MB_ERROR_OTHER; } if (error == MB_ERROR_NO_ERROR && (kind == nav_source || (kind == MB_DATA_DATA && use_ping_data == MB_YES))) { status = mb_extract_nav(verbose,imbio_ptr, store_ptr,&kind, ping[nbuff].time_i, &ping[nbuff].time_d, &ping[nbuff].lon, &ping[nbuff].lat, &ping[nbuff].speed, &ping[nbuff].heading, &ping[nbuff].draft, &ping[nbuff].roll, &ping[nbuff].pitch, &ping[nbuff].heave, &error); } if (status == MB_SUCCESS) { /* get first time value if first record */ if (first_read == MB_NO) { file_start_time_d = ping[nbuff].time_d; first_read = MB_YES; } /* get original values */ ping[nbuff].id = nload; ping[nbuff].record = ping[nbuff].id + ndump_total; ping[nbuff].lon_org = ping[nbuff].lon; ping[nbuff].lat_org = ping[nbuff].lat; ping[nbuff].speed_org = ping[nbuff].speed; ping[nbuff].heading_org = ping[nbuff].heading; ping[nbuff].draft_org = ping[nbuff].draft; ping[nbuff].file_time_d = ping[nbuff].time_d - file_start_time_d; /* apply offsets */ ping[nbuff].lon += offset_lon; ping[nbuff].lat += offset_lat; /* set starting dr */ ping[nbuff].mean_ok = MB_NO; ping[nbuff].lon_dr = ping[nbuff].lon; ping[nbuff].lat_dr = ping[nbuff].lat; /* set everything deselected */ ping[nbuff].tint_select = MB_NO; ping[nbuff].lon_select = MB_NO; ping[nbuff].lat_select = MB_NO; ping[nbuff].speed_select = MB_NO; ping[nbuff].heading_select = MB_NO; ping[nbuff].draft_select = MB_NO; ping[nbuff].lonlat_flag = MB_NO; /* select repeated data */ if (nbuff > 0 && ping[nbuff].lon == ping[nbuff-1].lon && ping[nbuff].lat == ping[nbuff-1].lat) { ping[nbuff].lonlat_flag = MB_YES; } /* print output debug statements */ if (verbose >= 5) { fprintf(stderr,"\ndbg5 Next good data found in function <%s>:\n", function_name); fprintf(stderr,"dbg5 %4d %4d %4d %d/%d/%d %2.2d:%2.2d:%2.2d.%6.6d %11.6f %11.6f %5.2f %5.1f %5.1f %5.2f %5.2f %5.2f\n", nbuff,ping[nbuff].id,ping[nbuff].record, ping[nbuff].time_i[1],ping[nbuff].time_i[2], ping[nbuff].time_i[0],ping[nbuff].time_i[3], ping[nbuff].time_i[4],ping[nbuff].time_i[5], ping[nbuff].time_i[6], ping[nbuff].lon, ping[nbuff].lat, ping[nbuff].speed, ping[nbuff].heading, ping[nbuff].draft, ping[nbuff].roll, ping[nbuff].pitch, ping[nbuff].heave); } /* increment counting variables */ nbuff++; nload++; /* update message every 250 records */ if (nload % 250 == 0) { sprintf(string, "MBnavedit: %d records loaded so far...", nload); do_message_on(string); } } } while (error <= MB_ERROR_NO_ERROR && nbuff < MBNAVEDIT_BUFFER_SIZE); nload_total += nload; /* define success */ if (nbuff > 0) { status = MB_SUCCESS; error = MB_ERROR_NO_ERROR; } /* check for time stamp repeats */ timestamp_problem = MB_NO; for (i=0;i 1) { for (i=1;i 0) { data_show_size = 0; plot_start_time = ping[0].file_time_d; plot_end_time = ping[nbuff-1].file_time_d; nplot = nbuff; } /* calculate speed-made-good and course-made-good */ for (i=0;i= 1) { fprintf(stderr,"\n%d data records loaded from input file <%s>\n", nload,ifile); fprintf(stderr,"%d data records now in buffer\n",nbuff); fprintf(stderr,"Current data record: %d\n", current_id); fprintf(stderr,"Current global data record: %d\n", current_id + ndump_total); } /* put up warning if timestamp problem detected */ if (timestamp_problem == MB_YES) { do_error_dialog("Duplicate or reverse order time", "stamps detected!! Time interpolation", "available under Controls menu."); } /* update controls */ do_set_controls(); /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return (status); } /*--------------------------------------------------------------------*/ int mbnavedit_clear_screen() { /* local variables */ char *function_name = "mbnavedit_clear_screen"; int status = MB_SUCCESS; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); } /* clear screen */ xg_fillrectangle(mbnavedit_xgid,0,0, plot_width,NUMBER_PLOTS_MAX*plot_height, pixel_values[WHITE],XG_SOLIDLINE); /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return (status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_next_buffer(int *quit) { /* local variables */ char *function_name = "mbnavedit_action_next_buffer"; int status = MB_SUCCESS; int save_dumped = 0; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); } /* clear the screen */ status = mbnavedit_clear_screen(); /* set quit off */ *quit = MB_NO; /* check if a file has been opened */ if (file_open == MB_YES) { /* dump the buffer */ status = mbnavedit_dump_data(hold_size); /* load the buffer */ status = mbnavedit_load_data(); /* if end of file reached then dump last buffer and close file */ if (nload <= 0) { save_dumped = ndump; status = mbnavedit_dump_data(0); status = mbnavedit_close_file(); ndump = ndump + save_dumped; /* if in normal mode last next_buffer does not mean quit, if in gui mode it does mean quit */ if (gui_mode == MB_YES) *quit = MB_YES; else *quit = MB_NO; /* if quitting let the world know... */ if (*quit == MB_YES && verbose >= 1) fprintf(stderr,"\nQuitting MBnavedit\nBye Bye...\n"); } /* else plot it */ else { status = mbnavedit_plot_all(); } } /* if no file open set failure status */ else { status = MB_FAILURE; ndump = 0; nload = 0; current_id = 0; } /* reset data_save */ data_save = MB_NO; /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 quit: %d\n",*quit); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_offset() { /* local variables */ char *function_name = "mbnavedit_action_offset"; int status = MB_SUCCESS; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); } /* check if a file has been opened */ if (file_open == MB_YES) { /* apply position offsets to the data */ for (i=0;i= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_close() { /* local variables */ char *function_name = "mbnavedit_action_close"; int status = MB_SUCCESS; int save_nloaded = 0; int save_ndumped = 0; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); } /* clear the screen */ status = mbnavedit_clear_screen(); /* if file has been opened and browse mode just dump the current buffer and close the file */ if (file_open == MB_YES && output_mode == OUTPUT_MODE_BROWSE) { /* dump the buffer */ status = mbnavedit_dump_data(0); /* now close the file */ status = mbnavedit_close_file(); } /* if file has been opened deal with it */ else if (file_open == MB_YES) { /* dump and load until the end of the file is reached */ save_ndumped = 0; save_nloaded = 0; do { /* dump the buffer */ status = mbnavedit_dump_data(0); save_ndumped += ndump; /* load the buffer */ status = mbnavedit_load_data(); save_nloaded += nload; } while (nload > 0); ndump = save_ndumped; nload = save_nloaded; /* now close the file */ status = mbnavedit_close_file(); } else { ndump = 0; nload = 0; nbuff = 0; current_id = 0; status = MB_FAILURE; } /* reset data_save */ data_save = MB_NO; /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_done(int *quit) { /* local variables */ char *function_name = "mbnavedit_action_done"; int status = MB_SUCCESS; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); } /* if in normal mode done does not mean quit, if in gui mode done does mean quit */ if (gui_mode == MB_YES) *quit = MB_YES; else *quit = MB_NO; /* if quitting let the world know... */ if (*quit == MB_YES && verbose >= 1) fprintf(stderr,"\nShutting MBnavedit down without further ado...\n"); /* call routine to deal with saving the current file, if any */ if (file_open == MB_YES) status = mbnavedit_action_close(); /* if quitting let the world know... */ if (*quit == MB_YES && verbose >= 1) fprintf(stderr,"\nQuitting MBnavedit\nBye Bye...\n"); /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 quit: %d\n",*quit); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_quit() { /* local variables */ char *function_name = "mbnavedit_action_quit"; int status = MB_SUCCESS; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); } /* let the world know... */ if (verbose >= 1) fprintf(stderr,"\nShutting MBnavedit down without further ado...\n"); /* call routine to deal with saving the current file, if any */ if (file_open == MB_YES) status = mbnavedit_action_close(); /* let the world know... */ if (verbose >= 1) fprintf(stderr,"\nQuitting MBnavedit\nBye Bye...\n"); /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_step(int step) { /* local variables */ char *function_name = "mbnavedit_action_step"; int status = MB_SUCCESS; int old_id, new_id; int set; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 step: %d\n",step); } /* check if a file has been opened */ if (file_open == MB_YES && nbuff > 0) { /* if current time span includes last data don't step */ if (step >= 0 && plot_end_time < ping[nbuff-1].file_time_d) { plot_start_time = plot_start_time + step; plot_end_time = plot_start_time + data_show_size; } else if (step < 0 && plot_start_time > ping[0].file_time_d) { plot_start_time = plot_start_time + step; plot_end_time = plot_start_time + data_show_size; } /* get current start of plotting data */ set = MB_NO; old_id = current_id; for (i=0;i= plot_start_time) { new_id = i; set = MB_YES; } } if (new_id < 0) new_id = 0; if (new_id >= nbuff) new_id = nbuff - 1; if (step < 0 && new_id > 0 && new_id == old_id) new_id--; if (step > 0 && new_id < nbuff - 1 && new_id == old_id) new_id++; current_id = new_id; /* replot */ if (nbuff > 0) { status = mbnavedit_plot_all(); } /* set failure flag if no step was made */ if (new_id == old_id) status = MB_FAILURE; } /* if no file open set failure status */ else { status = MB_FAILURE; current_id = 0; } /* print out information */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 Current buffer values:\n"); fprintf(stderr,"dbg2 nload: %d\n",nload); fprintf(stderr,"dbg2 nbuff: %d\n",nbuff); fprintf(stderr,"dbg2 nbuff: %d\n",nbuff); fprintf(stderr,"dbg2 nbuff: %d\n",nbuff); fprintf(stderr,"dbg2 current_id: %d\n",current_id); } /* reset data_save */ data_save = MB_NO; /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_mouse_pick(int xx, int yy) { /* local variables */ char *function_name = "mbnavedit_action_mouse_pick"; int status = MB_SUCCESS; int deselect; int iplot; int active_plot; int range, range_min; int iping; int ix, iy; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 xx: %d\n",xx); fprintf(stderr,"dbg2 yy: %d\n",yy); } /* don't try to do anything if no data */ active_plot = -1; if (nplot > 0) { /* figure out which plot the cursor is in */ for (iplot=0;iplot= mbnavplot[iplot].ixmin && xx <= mbnavplot[iplot].ixmax && yy <= mbnavplot[iplot].iymin && yy >= mbnavplot[iplot].iymax) active_plot = iplot; } } /* don't try to do anything if no data or not in plot */ if (nplot > 0 && active_plot > -1) { /* deselect everything in non-active plots */ deselect = MB_NO; for (iplot=0;iplot= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_mouse_select(int xx, int yy) { /* local variables */ char *function_name = "mbnavedit_action_mouse_select"; int status = MB_SUCCESS; int deselect; int iplot; int active_plot; int range; int ix, iy; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 xx: %d\n",xx); fprintf(stderr,"dbg2 yy: %d\n",yy); } /* don't try to do anything if no data */ active_plot = -1; if (nplot > 0) { /* figure out which plot the cursor is in */ for (iplot=0;iplot= mbnavplot[iplot].ixmin && xx <= mbnavplot[iplot].ixmax && yy <= mbnavplot[iplot].iymin && yy >= mbnavplot[iplot].iymax) active_plot = iplot; } } /* don't try to do anything if no data or not in plot */ if (nplot > 0 && active_plot > -1) { /* deselect everything in non-active plots */ deselect = MB_NO; for (iplot=0;iplot= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_mouse_deselect(int xx, int yy) { /* local variables */ char *function_name = "mbnavedit_action_mouse_deselect"; int status = MB_SUCCESS; int deselect; int iplot; int active_plot; int range; int ix, iy; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 xx: %d\n",xx); fprintf(stderr,"dbg2 yy: %d\n",yy); } /* don't try to do anything if no data */ active_plot = -1; if (nplot > 0) { /* figure out which plot the cursor is in */ for (iplot=0;iplot= mbnavplot[iplot].ixmin && xx <= mbnavplot[iplot].ixmax && yy <= mbnavplot[iplot].iymin && yy >= mbnavplot[iplot].iymax) active_plot = iplot; } } /* don't try to do anything if no data or not in plot */ if (nplot > 0 && active_plot > -1) { /* deselect everything in non-active plots */ deselect = MB_NO; for (iplot=0;iplot= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_mouse_selectall(int xx, int yy) { /* local variables */ char *function_name = "mbnavedit_action_mouse_selectall"; int status = MB_SUCCESS; int iplot; int active_plot; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 xx: %d\n",xx); fprintf(stderr,"dbg2 yy: %d\n",yy); } /* don't try to do anything if no data */ active_plot = -1; if (nplot > 0) { /* figure out which plot the cursor is in */ for (iplot=0;iplot= mbnavplot[iplot].ixmin && xx <= mbnavplot[iplot].ixmax && yy <= mbnavplot[iplot].iymin && yy >= mbnavplot[iplot].iymax) active_plot = iplot; } } /* don't try to do anything if no data or not in plot */ if (nplot > 0 && active_plot > -1) { /* deselect everything in non-active plots */ for (iplot=0;iplot= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_mouse_deselectall(int xx, int yy) { /* local variables */ char *function_name = "mbnavedit_action_mouse_deselectall"; int status = MB_SUCCESS; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 xx: %d\n",xx); fprintf(stderr,"dbg2 yy: %d\n",yy); } /* don't try to do anything if no data */ if (nplot > 0) { /* deselect all data points in all plots - this logic follows from deselecting all active plots plus all non-active plots */ for (i=current_id;i= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_deselect_all(int type) { /* local variables */ char *function_name = "mbnavedit_action_mouse_deselect_all"; int status = MB_SUCCESS; int ndeselect; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 type: %d\n",type); } /* don't try to do anything if no data */ if (nplot > 0) { /* deselect all data points in specified data type */ ndeselect = 0; for (i=0;i 0) status = MB_SUCCESS; else status = MB_FAILURE; } /* if no data then set failure flag */ else status = MB_FAILURE; /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_set_interval(int xx, int yy, int which) { /* local variables */ char *function_name = "mbnavedit_action_set_interval"; int status = MB_SUCCESS; static int interval_bound1; static int interval_bound2; static double interval_time1; static double interval_time2; static int interval_set1 = MB_NO; static int interval_set2 = MB_NO; int itmp; double dtmp; int set; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 xx: %d\n",xx); fprintf(stderr,"dbg2 yy: %d\n",yy); fprintf(stderr,"dbg2 which: %d\n",which); } /* don't try to do anything if no data */ if (nplot > 0 && number_plots > 0) { /* if which = 0 set first bound and draw dashed lines */ if (which == 0) { /* unplot old line on all plots */ if (interval_set1 == MB_YES) for (i=0;i mbnavplot[0].ixmax) xx = mbnavplot[0].ixmax; /* get lower bound time and location */ interval_bound1 = xx; interval_time1 = mbnavplot[0].xmin + (xx - mbnavplot[0].ixmin)/mbnavplot[0].xscale; interval_set1 = MB_YES; /* plot line on all plots */ for (i=0;i mbnavplot[0].ixmax) xx = mbnavplot[0].ixmax; /* get lower bound time and location */ interval_bound2 = xx; interval_time2 = mbnavplot[0].xmin + (xx - mbnavplot[0].ixmin)/mbnavplot[0].xscale; interval_set2 = MB_YES; /* plot line on all plots */ for (i=0;i interval_bound2) { itmp = interval_bound2; dtmp = interval_time2; interval_bound2 = interval_bound1; interval_time2 = interval_time1; interval_bound1 = itmp; interval_time1 = dtmp; } /* reset plotting parameters */ plot_start_time = interval_time1; plot_end_time = interval_time2; data_show_size = plot_end_time - plot_start_time; /* reset time stepping parameters */ data_step_size = data_show_size / 4; if (data_step_size > data_step_max) data_step_max = 2 * data_step_size; /* get current start of plotting data */ set = MB_NO; for (i=0;i= plot_start_time) { current_id = i; set = MB_YES; } } if (current_id < 0) current_id = 0; if (current_id >= nbuff) current_id = nbuff - 1; /* replot */ mbnavedit_plot_all(); } /* else if which = 3 unset bounds */ else if (which == 3) { interval_set1 = MB_NO; interval_set2 = MB_NO; } /* else failure */ else status = MB_FAILURE; } /* if no data then set failure flag */ else status = MB_FAILURE; /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_use_dr() { /* local variables */ char *function_name = "mbnavedit_action_use_dr"; int status = MB_SUCCESS; int iplot; int active_plot; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); } /* don't try to do anything if no data */ if (nplot > 0) { /* make sure either a lon or lat plot is active */ active_plot = -1; for (iplot=0;iplot -1) { for (i=current_id;i= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_use_smg() { /* local variables */ char *function_name = "mbnavedit_action_use_smg"; int status = MB_SUCCESS; int iplot; int active_plot; int speedheading_change; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); } /* don't try to do anything if no data */ if (nplot > 0) { /* figure out which plot is speed */ active_plot = -1; for (iplot=0;iplot -1) { speedheading_change = MB_NO; for (i=current_id;i= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_use_cmg() { /* local variables */ char *function_name = "mbnavedit_action_use_cmg"; int status = MB_SUCCESS; int iplot; int active_plot; int speedheading_change; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); } /* don't try to do anything if no data */ if (nplot > 0) { /* figure out which plot is heading */ active_plot = -1; for (iplot=0;iplot -1) { speedheading_change = MB_NO; for (i=current_id;i= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_interpolate() { /* local variables */ char *function_name = "mbnavedit_action_interpolate"; int status = MB_SUCCESS; int iping; int ibefore, iafter; int timelonlat_change; int speedheading_change; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); } /* don't try to do anything if no data */ if (nplot > 0) { /* look for position or time changes */ timelonlat_change = MB_NO; speedheading_change = MB_NO; /* do expected time */ for (iping=0;iping=0;i--) if (ping[i].tint_select == MB_NO && ibefore == iping) ibefore = i; iafter = iping; for (i=iping+1;i iping) { ping[iping].time_d = ping[ibefore].time_d + (ping[iafter].time_d - ping[ibefore].time_d) *((double)(iping - ibefore)) / ((double)(iafter - ibefore)); ping[iping].tint = ping[iping].time_d - ping[iping-1].time_d; timelonlat_change = MB_YES; } else if (ibefore < iping && ibefore > 0) { ping[iping].time_d = ping[ibefore].time_d + (ping[ibefore].time_d - ping[ibefore-1].time_d) * (iping - ibefore); ping[iping].tint = ping[iping].time_d - ping[iping-1].time_d; timelonlat_change = MB_YES; } else if (ibefore < iping) { ping[iping].time_d = ping[ibefore].time_d; ping[iping].tint = ping[iping].time_d - ping[iping-1].time_d; timelonlat_change = MB_YES; } else if (iafter > iping && iafter < nbuff - 1) { ping[iping].time_d = ping[iafter].time_d + (ping[iafter+1].time_d - ping[iafter].time_d) *(iping - iafter); ping[iping].tint = 0.0; timelonlat_change = MB_YES; } else if (iafter > iping) { ping[iping].time_d = ping[iafter].time_d; ping[iping].tint = ping[iping].time_d - ping[iping-1].time_d; timelonlat_change = MB_YES; } ping[iping].file_time_d = ping[iping].time_d - file_start_time_d; status = mb_get_date(verbose, ping[iping].time_d, ping[iping].time_i); if (iping < nbuff - 1) if (ping[iping+1].tint_select == MB_NO) ping[iping+1].tint = ping[iping+1].time_d - ping[iping].time_d; } } /* do longitude */ for (iping=0;iping=0;i--) if (ping[i].lon_select == MB_NO && ibefore == iping) ibefore = i; iafter = iping; for (i=iping+1;i iping) { ping[iping].lon = ping[ibefore].lon + (ping[iafter].lon - ping[ibefore].lon) *(ping[iping].time_d - ping[ibefore].time_d) /(ping[iafter].time_d - ping[ibefore].time_d); timelonlat_change = MB_YES; } else if (ibefore < iping && ibefore > 0) { ping[iping].lon = ping[ibefore].lon + (ping[ibefore].lon - ping[ibefore-1].lon) *(ping[iping].time_d - ping[ibefore].time_d) /(ping[ibefore].time_d - ping[ibefore-1].time_d); timelonlat_change = MB_YES; } else if (ibefore < iping) { ping[iping].lon = ping[ibefore].lon; timelonlat_change = MB_YES; } else if (iafter > iping && iafter < nbuff - 1) { ping[iping].lon = ping[iafter].lon + (ping[iafter+1].lon - ping[iafter].lon) *(ping[iping].time_d - ping[iafter].time_d) /(ping[iafter+1].time_d - ping[iafter].time_d); timelonlat_change = MB_YES; } else if (iafter > iping) { ping[iping].lon = ping[iafter].lon; timelonlat_change = MB_YES; } } } /* do latitude */ for (iping=0;iping=0;i--) if (ping[i].lat_select == MB_NO && ibefore == iping) ibefore = i; iafter = iping; for (i=iping+1;i iping) { ping[iping].lat = ping[ibefore].lat + (ping[iafter].lat - ping[ibefore].lat) *(ping[iping].time_d - ping[ibefore].time_d) /(ping[iafter].time_d - ping[ibefore].time_d); timelonlat_change = MB_YES; } else if (ibefore < iping && ibefore > 0) { ping[iping].lat = ping[ibefore].lat + (ping[ibefore].lat - ping[ibefore-1].lat) *(ping[iping].time_d - ping[ibefore].time_d) /(ping[ibefore].time_d - ping[ibefore-1].time_d); timelonlat_change = MB_YES; } else if (ibefore < iping) { ping[iping].lat = ping[ibefore].lat; timelonlat_change = MB_YES; } else if (iafter > iping && iafter < nbuff - 1) { ping[iping].lat = ping[iafter].lat + (ping[iafter+1].lat - ping[iafter].lat) *(ping[iping].time_d - ping[iafter].time_d) /(ping[iafter+1].time_d - ping[iafter].time_d); timelonlat_change = MB_YES; } else if (iafter > iping) { ping[iping].lat = ping[iafter].lat; timelonlat_change = MB_YES; } } } /* do speed */ for (iping=0;iping=0;i--) if (ping[i].speed_select == MB_NO && ibefore == iping) ibefore = i; iafter = iping; for (i=iping+1;i iping) { ping[iping].speed = ping[ibefore].speed + (ping[iafter].speed - ping[ibefore].speed) *(ping[iping].time_d - ping[ibefore].time_d) /(ping[iafter].time_d - ping[ibefore].time_d); speedheading_change = MB_YES; } else if (ibefore < iping) { ping[iping].speed = ping[ibefore].speed; speedheading_change = MB_YES; } else if (iafter > iping) { ping[iping].speed = ping[iafter].speed; speedheading_change = MB_YES; } } } /* do heading */ for (iping=0;iping=0;i--) if (ping[i].heading_select == MB_NO && ibefore == iping) ibefore = i; iafter = iping; for (i=iping+1;i iping) { ping[iping].heading = ping[ibefore].heading + (ping[iafter].heading - ping[ibefore].heading) *(ping[iping].time_d - ping[ibefore].time_d) /(ping[iafter].time_d - ping[ibefore].time_d); speedheading_change = MB_YES; } else if (ibefore < iping) { ping[iping].heading = ping[ibefore].heading; speedheading_change = MB_YES; } else if (iafter > iping) { ping[iping].heading = ping[iafter].heading; speedheading_change = MB_YES; } } } /* do draft */ for (iping=0;iping=0;i--) if (ping[i].draft_select == MB_NO && ibefore == iping) ibefore = i; iafter = iping; for (i=iping+1;i iping) { ping[iping].draft = ping[ibefore].draft + (ping[iafter].draft - ping[ibefore].draft) *(ping[iping].time_d - ping[ibefore].time_d) /(ping[iafter].time_d - ping[ibefore].time_d); timelonlat_change = MB_YES; } else if (ibefore < iping) { ping[iping].draft = ping[ibefore].draft; timelonlat_change = MB_YES; } else if (iafter > iping) { ping[iping].draft = ping[iafter].draft; } } } /* recalculate speed-made-good and course-made-good */ if (timelonlat_change == MB_YES) for (i=0;i= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_interpolaterepeats() { /* local variables */ char *function_name = "mbnavedit_action_interpolaterepeats"; int status = MB_SUCCESS; int iping; int ibefore, iafter; int timelonlat_change; int speedheading_change; int found; int i, j; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); } /* don't try to do anything if no data */ if (nplot > 0) { /* look for position or time changes */ timelonlat_change = MB_NO; speedheading_change = MB_NO; /* do expected time */ for (iping=1;iping= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_revert() { /* local variables */ char *function_name = "mbnavedit_action_revert"; int status = MB_SUCCESS; int iplot; int timelonlat_change; int speedheading_change; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); } /* don't try to do anything if no data */ if (nplot > 0) { /* look for position changes */ timelonlat_change = MB_NO; speedheading_change = MB_NO; /* loop over each of the plots */ for (iplot=0;iplot= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_flag() { /* local variables */ char *function_name = "mbnavedit_action_flag"; int status = MB_SUCCESS; int iplot; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); } /* don't try to do anything if no data */ if (nplot > 0) { /* loop over each of the plots */ for (iplot=0;iplot= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_unflag() { /* local variables */ char *function_name = "mbnavedit_action_unflag"; int status = MB_SUCCESS; int iplot; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); } /* don't try to do anything if no data */ if (nplot > 0) { /* loop over each of the plots */ for (iplot=0;iplot= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_fixtime() { /* local variables */ char *function_name = "mbnavedit_action_fixtime"; int status = MB_SUCCESS; int istart, iend; double start_time_d, end_time_d; int i, j; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); } /* loop over the data */ for (i=0;i start_time_d) { iend = i; end_time_d = ping[i].time_d; for (j=istart+1;j0) ping[j-1].tint = ping[j].time_d - ping[j-1].time_d; if (j= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_deletebadtime() { /* local variables */ char *function_name = "mbnavedit_action_deletebadtime"; int status = MB_SUCCESS; double lastgood_time_d; int i, j, nbuffnew; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); } /* loop over the data looking for bad times */ lastgood_time_d = ping[0].time_d; for (i=1;i 60.0) { if (i == nbuff - 1) ping[i].id = -1; else if (ping[i+1].time_d - ping[i].time_d <= 0.0) ping[i].id = -1; else lastgood_time_d = ping[i].time_d; } else if (ping[i].time_d > ping[nbuff-1].time_d) { ping[i].id = -1; } else { lastgood_time_d = ping[i].time_d; } } /* loop over the data in reverse deleting data with bad times */ nbuffnew = nbuff; for (i=nbuff-1;i>=0;i--) { if (ping[i].id == -1) { for (j=i;j 0) ping[i-1].tint = ping[i].time_d - ping[i-1].time_d; if (i == nbuffnew - 2 && i > 0) ping[i].tint = ping[i-1].tint; else if (i == nbuffnew - 2 && i == 0) ping[i].tint = 0.0; nbuffnew--; } } fprintf(stderr,"Data deleted: nbuff:%d nbuffnew:%d\n",nbuff,nbuffnew); nbuff = nbuffnew; /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_action_showall() { /* local variables */ char *function_name = "mbnavedit_action_showall"; int status = MB_SUCCESS; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); } /* reset plotting time span */ if (nbuff > 0) { plot_start_time = ping[0].file_time_d; plot_end_time = ping[nbuff-1].file_time_d; data_show_size = 0; current_id = 0; } /* replot */ status = mbnavedit_plot_all(); /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_get_smgcmg(int i) { /* local variables */ char *function_name = "mbnavedit_get_smgcmg"; int status = MB_SUCCESS; double time_d1, lon1, lat1; double time_d2, lon2, lat2; double mtodeglon, mtodeglat; double del_time, dx, dy, dist; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 i: %d\n",i); } /* calculate speed made good and course made for ping i */ if (i < nbuff) { if (i == 0) { time_d1 = ping[i].time_d; lon1 = ping[i].lon; lat1 = ping[i].lat; time_d2 = ping[i+1].time_d; lon2 = ping[i+1].lon; lat2 = ping[i+1].lat; } else if (i == nbuff - 1) { time_d1 = ping[i-1].time_d; lon1 = ping[i-1].lon; lat1 = ping[i-1].lat; time_d2 = ping[i].time_d; lon2 = ping[i].lon; lat2 = ping[i].lat; } else { time_d1 = ping[i-1].time_d; lon1 = ping[i-1].lon; lat1 = ping[i-1].lat; time_d2 = ping[i].time_d; lon2 = ping[i].lon; lat2 = ping[i].lat; } mb_coor_scale(verbose,lat1,&mtodeglon,&mtodeglat); del_time = time_d2 - time_d1; dx = (lon2 - lon1)/mtodeglon; dy = (lat2 - lat1)/mtodeglat; dist = sqrt(dx*dx + dy*dy); if (del_time > 0.0) ping[i].speed_made_good = 3.6*dist/del_time; else ping[i].speed_made_good = 0.0; if (dist > 0.0) ping[i].course_made_good = RTD*atan2(dx/dist,dy/dist); else ping[i].course_made_good = ping[i].heading; if (ping[i].course_made_good < 0.0) ping[i].course_made_good = ping[i].course_made_good + 360.0; status = MB_SUCCESS; } else status = MB_FAILURE; /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_get_model() { /* local variables */ char *function_name = "mbnavedit_get_model"; int status = MB_SUCCESS; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); } /* only model if data available */ if (nbuff > 0) { /* call correct modeling function */ if (model_mode == MODEL_MODE_MEAN) mbnavedit_get_gaussianmean(); else if (model_mode == MODEL_MODE_DR) mbnavedit_get_dr(); else if (model_mode == MODEL_MODE_INVERT) mbnavedit_get_inversion(); } /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_get_gaussianmean() { /* local variables */ char *function_name = "mbnavedit_get_gaussianmean"; int status = MB_SUCCESS; double timewindow; double sumlon, sumlat, w, weight; double dt, a; int jstart, nsum, npos, nneg; int jbefore, jafter; int i, j; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); } /* loop over navigation calculating gaussian mean positions */ timewindow = 0.1 * mean_time_window; a = -4.0 / (timewindow * timewindow); jstart = 0; for (i=0;i= 0.0) npos++; weight += w; sumlon += w * ping[j].lon; sumlat += w * ping[j].lat; if (nsum == 1) jstart = j; } } if (npos > 0 && nneg > 0) { ping[i].mean_ok = MB_YES; ping[i].lon_dr = sumlon / weight; ping[i].lat_dr = sumlat / weight; } else { ping[i].mean_ok = MB_NO; ping[i].lon_dr = ping[i].lon; ping[i].lat_dr = ping[i].lat; } } /* loop over navigation performing linear interpolation to fill gaps */ jbefore = -1; for (i=0;i i && ping[j].lonlat_flag == MB_NO) jafter = j; } if (jbefore >= 0 && jafter > i) { dt = (ping[i].time_d - ping[jbefore].time_d) / (ping[jafter].time_d - ping[jbefore].time_d); ping[i].lon_dr = ping[jbefore].lon + dt * (ping[jafter].lon - ping[jbefore].lon); ping[i].lat_dr = ping[jbefore].lat + dt * (ping[jafter].lat - ping[jbefore].lat); } else if (jbefore >= 0) { ping[i].lon_dr = ping[jbefore].lon; ping[i].lat_dr = ping[jbefore].lat; } else if (jafter > i) { ping[i].lon_dr = ping[jafter].lon; ping[i].lat_dr = ping[jafter].lat; } else { ping[i].lon_dr = ping[i].lon; ping[i].lat_dr = ping[i].lat; } } } /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_get_dr() { /* local variables */ char *function_name = "mbnavedit_get_dr"; int status = MB_SUCCESS; double mtodeglon, mtodeglat; double del_time, dx, dy; double driftlon, driftlat; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); } /* calculate dead reckoning */ driftlon = 0.00001 * drift_lon; driftlat = 0.00001 * drift_lat; for (i=0;i= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_get_inversion() { /* local variables */ char *function_name = "mbnavedit_get_inversion"; int status = MB_SUCCESS; int nnz; int nrows; int ncols; double *a; int *ia; int *nia; int *nx; double *x; double *dx; double *d; double *sigma; double *work; int ncyc, nsig; double smax, sup, err, supt, slo, errlsq; int ncycle; double bandwidth; int nr, nc; int nlon_avg, nlat_avg; double lon_avg, lat_avg; double mtodeglon, mtodeglat; double dtime_d, dtime_d_sq; char string[MB_PATH_MAXLINE]; int first, last; int i, ii, j, k; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); } /* set maximum dimensions of the inverse problem */ nrows = nplot + (nplot - 1) + (nplot - 2); ncols = nplot; nnz = 3; ncycle = 512; bandwidth = 10000.0; /* get average lon value */ lon_avg = 0.0; nlon_avg = 0; lat_avg = 0.0; nlat_avg = 0; first = current_id; last = current_id; for (i=current_id;i 0) lon_avg /= nlon_avg; if (nlat_avg > 0) lat_avg /= nlat_avg; mb_coor_scale(verbose, lat_avg, &mtodeglon, &mtodeglat); /* allocate space for the inverse problem */ status = mb_mallocd(verbose, __FILE__, __LINE__, nnz * nrows * sizeof(double), (void **)&a,&error); status = mb_mallocd(verbose, __FILE__, __LINE__, nnz * nrows * sizeof(int), (void **)&ia,&error); status = mb_mallocd(verbose, __FILE__, __LINE__, nrows * sizeof(int), (void **)&nia,&error); status = mb_mallocd(verbose, __FILE__, __LINE__, nrows * sizeof(double), (void **)&d,&error); status = mb_mallocd(verbose, __FILE__, __LINE__, ncols * sizeof(double), (void **)&x,&error); status = mb_mallocd(verbose, __FILE__, __LINE__, ncols * sizeof(int), (void **)&nx,&error); status = mb_mallocd(verbose, __FILE__, __LINE__, ncols * sizeof(double), (void **)&dx,&error); status = mb_mallocd(verbose, __FILE__, __LINE__, ncycle * sizeof(double), (void **)&sigma,&error); status = mb_mallocd(verbose, __FILE__, __LINE__, ncycle * sizeof(double), (void **)&work,&error); /* do inversion */ if (error == MB_ERROR_NO_ERROR) { /* set message */ sprintf(string,"Setting up inversion of %d longitude points", nplot); do_message_on(string); /* initialize arrays */ for (i=0;i 0.0 && ii > 0 && ping[i].time_d > ping[i-1].time_d) { /* get time difference */ dtime_d = ping[i].time_d - ping[i-1].time_d; /* constrain lon speed */ k = nnz * nr; d[nr] = 0.0; nia[nr] = 2; ia[k] = ii - 1; a[k] = -weight_speed / dtime_d; ia[k+1] = ii; a[k+1] = weight_speed / dtime_d; nr++; } /* constrain acceleration */ if (weight_acceleration > 0.0 && ii > 0 && ii < nplot - 1 && ping[i+1].time_d > ping[i-1].time_d) { /* get time difference */ dtime_d = ping[i+1].time_d - ping[i-1].time_d; dtime_d_sq = dtime_d * dtime_d; /* constrain lon acceleration */ k = nnz * nr; d[nr] = 0.0; nia[nr] = 3; ia[k] = ii - 1; a[k] = weight_acceleration / dtime_d_sq; ia[k+1] = ii; a[k+1] = -2.0 * weight_acceleration / dtime_d_sq; ia[k+2] = ii + 1; a[k+2] = weight_acceleration / dtime_d_sq; nr++; } } /* set message */ sprintf(string,"Inverting %dX%d for smooth longitude...", nc, nr); do_message_on(string); /* compute upper bound on maximum eigenvalue */ ncyc = 0; nsig = 0; lspeig(a, ia, nia, nnz, nc, nr, ncyc, &nsig, x, dx, sigma, work, &smax, &err, &sup); supt = smax + err; if (sup > supt) supt = sup; if (verbose > 1) fprintf(stderr, "Initial lspeig: %g %g %g %g\n", sup, smax, err, supt); ncyc = 16; for (i=0;i<4;i++) { lspeig(a, ia, nia, nnz, nc, nr, ncyc, &nsig, x, dx, sigma, work, &smax, &err, &sup); supt = smax + err; if (sup > supt) supt = sup; if (verbose > 1) fprintf(stderr, "lspeig[%d]: %g %g %g %g\n", i, sup, smax, err, supt); } /* calculate chebyshev factors (errlsq is the theoretical error) */ slo = supt / bandwidth; chebyu(sigma, ncycle, supt, slo, work); errlsq = errlim(sigma, ncycle, supt, slo); if (verbose > 1) fprintf(stderr, "Theoretical error: %f\n", errlsq); if (verbose > 1) for (i=0;i 0.0 && ii > 0 && ping[i].time_d > ping[i-1].time_d) { /* get time difference */ dtime_d = ping[i].time_d - ping[i-1].time_d; /* constrain lat speed */ k = nnz * nr; d[nr] = 0.0; nia[nr] = 2; ia[k] = ii - 1; a[k] = -weight_speed / dtime_d; ia[k+1] = ii; a[k+1] = weight_speed / dtime_d; nr++; } /* constrain acceleration */ if (weight_acceleration > 0.0 && ii > 0 && ii < nplot - 1 && ping[i+1].time_d > ping[i-1].time_d) { /* get time difference */ dtime_d = ping[i+1].time_d - ping[i-1].time_d; dtime_d_sq = dtime_d * dtime_d; /* constrain lat acceleration */ k = nnz * nr; d[nr] = 0.0; nia[nr] = 3; ia[k] = ii - 1; a[k] = weight_acceleration / dtime_d_sq; ia[k+1] = ii; a[k+1] = -2.0 * weight_acceleration / dtime_d_sq; ia[k+2] = ii + 1; a[k+2] = weight_acceleration / dtime_d_sq; nr++; } } /* set message */ sprintf(string,"Inverting %dX%d for smooth latitude...", nc, nr); do_message_on(string); /* compute upper bound on maximum eigenvalue */ ncyc = 0; nsig = 0; lspeig(a, ia, nia, nnz, nc, nr, ncyc, &nsig, x, dx, sigma, work, &smax, &err, &sup); supt = smax + err; if (sup > supt) supt = sup; if (verbose > 1) fprintf(stderr, "Initial lspeig: %g %g %g %g\n", sup, smax, err, supt); ncyc = 16; for (i=0;i<4;i++) { lspeig(a, ia, nia, nnz, nc, nr, ncyc, &nsig, x, dx, sigma, work, &smax, &err, &sup); supt = smax + err; if (sup > supt) supt = sup; if (verbose > 1) fprintf(stderr, "lspeig[%d]: %g %g %g %g\n", i, sup, smax, err, supt); } /* calculate chebyshev factors (errlsq is the theoretical error) */ slo = supt / bandwidth; chebyu(sigma, ncycle, supt, slo, work); errlsq = errlim(sigma, ncycle, supt, slo); if (verbose > 1) fprintf(stderr, "Theoretical error: %f\n", errlsq); if (verbose > 1) for (i=0;i= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return(status); } /*--------------------------------------------------------------------*/ int mbnavedit_plot_all() { /* local variables */ char *function_name = "mbnavedit_plot_all"; int status = MB_SUCCESS; double time_min; double time_max; double tint_min; double tint_max; double lon_min; double lon_max; double lat_min; double lat_max; double speed_min; double speed_max; double heading_min; double heading_max; double draft_min; double draft_max; double roll_min; double roll_max; double pitch_min; double pitch_max; double heave_min; double heave_max; double center, range; int margin_x, margin_y; int iyzero; int iplot; int center_x, center_y; double dx, x; int xtime_i[7]; int ix; int swidth, sascent, sdescent; char yformat[10]; int i; char string[MB_PATH_MAXLINE]; int fpx, fpdx, fpy, fpdy; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); } /* figure out which pings to plot */ nplot = 0; if (data_show_size > 0 && nbuff > 0) { plot_start_time = ping[current_id].file_time_d; plot_end_time = plot_start_time + data_show_size; for (i=current_id;i 0) { plot_start_time = ping[0].file_time_d; plot_end_time = ping[nbuff-1].file_time_d; data_show_size = plot_end_time - plot_start_time + 1; if (data_show_max < data_show_size) data_show_max = data_show_size; nplot = nbuff; } /* deselect data outside plots */ for (i=0;i 0) { /* find min max values */ time_min = plot_start_time; time_max = plot_end_time; tint_min = ping[current_id].tint; tint_max = ping[current_id].tint; lon_min = ping[current_id].lon; lon_max = ping[current_id].lon; lat_min = ping[current_id].lat; lat_max = ping[current_id].lat; speed_min = 0.0; speed_max = ping[current_id].speed; heading_min = ping[current_id].heading; heading_max = ping[current_id].heading; draft_min = ping[current_id].draft; draft_max = ping[current_id].draft; roll_min = ping[current_id].roll; roll_max = ping[current_id].roll; pitch_min = ping[current_id].pitch; pitch_max = ping[current_id].pitch; heave_min = ping[current_id].heave; heave_max = ping[current_id].heave; for (i=current_id+1;i (lat_max - lat_min)) { center = 0.5*(lat_min + lat_max); lat_min = center - 0.5*(lon_max - lon_min); lat_max = center + 0.5*(lon_max - lon_min); } else { center = 0.5*(lon_min + lon_max); lon_min = center - 0.5*(lat_max - lat_min); lon_max = center + 0.5*(lat_max - lat_min); } } /* make sure min max values aren't too small */ if ((tint_max - tint_min) < 0.01) { center = 0.5*(tint_min + tint_max); tint_min = center - 0.005; tint_max = center + 0.005; } if ((lon_max - lon_min) < 0.001) { center = 0.5*(lon_min + lon_max); lon_min = center - 0.0005; lon_max = center + 0.0005; } if ((lat_max - lat_min) < 0.001) { center = 0.5*(lat_min + lat_max); lat_min = center - 0.0005; lat_max = center + 0.0005; } if (speed_max < 10.0) speed_max = 10.0; if ((heading_max - heading_min) < 10.0) { center = 0.5*(heading_min + heading_max); heading_min = center - 5; heading_max = center + 5; } if ((draft_max - draft_min) < 10.0) { center = 0.5*(draft_min + draft_max); draft_min = center - 5; draft_max = center + 5; } if ((roll_max - roll_min) < 2.0) { center = 0.5*(roll_min + roll_max); roll_min = center - 1; roll_max = center + 1; } if ((pitch_max - pitch_min) < 2.0) { center = 0.5*(pitch_min + pitch_max); pitch_min = center - 1; pitch_max = center + 1; } if ((heave_max - heave_min) < 0.02) { center = 0.5*(heave_min + heave_max); heave_min = center - 0.01; heave_max = center + 0.01; } /* print out information */ if (verbose >= 2) { fprintf(stderr,"\n%d data records set for plotting (%d desired)\n", nplot,data_show_size); for (i=current_id;i 4 * margin_x) fpx = center_x + 2 * margin_x - fpdx; /* plot file position bar */ xg_drawrectangle(mbnavedit_xgid, center_x - 2 * margin_x, fpy, 4 * margin_x, fpdy, pixel_values[BLACK],XG_SOLIDLINE); xg_drawrectangle(mbnavedit_xgid, center_x - 2 * margin_x - 1, fpy - 1, 4 * margin_x + 2, fpdy + 2, pixel_values[BLACK],XG_SOLIDLINE); xg_fillrectangle(mbnavedit_xgid, fpx, fpy, fpdx, fpdy, pixel_values[LIGHTGREY],XG_SOLIDLINE); xg_drawrectangle(mbnavedit_xgid, fpx, fpy, fpdx, fpdy, pixel_values[BLACK],XG_SOLIDLINE); sprintf(string,"0 "); xg_justify(mbnavedit_xgid,string,&swidth, &sascent,&sdescent); xg_drawstring(mbnavedit_xgid, (int)(center_x - 2 * margin_x - swidth), fpy + sascent, string, pixel_values[BLACK],XG_SOLIDLINE); sprintf(string," %d", nbuff); xg_drawstring(mbnavedit_xgid, (int)(center_x + 2 * margin_x), fpy + sascent, string, pixel_values[BLACK],XG_SOLIDLINE); /* plot x label */ xg_justify(mbnavedit_xgid, mbnavplot[iplot].xlabel, &swidth,&sascent,&sdescent); xg_drawstring(mbnavedit_xgid, (int)(center_x - swidth/2), (int)(mbnavplot[iplot].iymin+0.75*margin_y), mbnavplot[iplot].xlabel, pixel_values[BLACK],XG_SOLIDLINE); /* plot y labels */ xg_justify(mbnavedit_xgid, mbnavplot[iplot].ylabel1, &swidth,&sascent,&sdescent); xg_drawstring(mbnavedit_xgid, (int)(mbnavplot[iplot].ixmin - swidth/2 - 0.75*margin_x), (int)(center_y-sascent), mbnavplot[iplot].ylabel1, pixel_values[BLACK],XG_SOLIDLINE); xg_justify(mbnavedit_xgid, mbnavplot[iplot].ylabel2, &swidth,&sascent,&sdescent); xg_drawstring(mbnavedit_xgid, (int)(mbnavplot[iplot].ixmin - swidth/2 - 0.75*margin_x), (int)(center_y+2*sascent), mbnavplot[iplot].ylabel2, pixel_values[BLACK],XG_SOLIDLINE); /* plot x axis time annotation */ dx = (plot_end_time - plot_start_time) / 5; for (i=0;i<6;i++) { /* get x position */ x = plot_start_time + i * dx; ix = mbnavplot[iplot].ixmin + mbnavplot[iplot].xscale * (x - mbnavplot[iplot].xmin); x += file_start_time_d; /* draw tickmarks */ xg_drawline(mbnavedit_xgid, ix, mbnavplot[iplot].iymin, ix, mbnavplot[iplot].iymin + 5, pixel_values[BLACK],XG_SOLIDLINE); /* draw annotations */ mb_get_date(verbose, x, xtime_i); sprintf(string, "%2.2d:%2.2d:%2.2d.%3.3d", xtime_i[3], xtime_i[4], xtime_i[5], (int)(0.001 * xtime_i[6])); xg_justify(mbnavedit_xgid, string, &swidth,&sascent,&sdescent); xg_drawstring(mbnavedit_xgid, (int)(ix - swidth/2), (int)(mbnavplot[iplot].iymin + 5 + 1.75*sascent), string, pixel_values[BLACK],XG_SOLIDLINE); } /* plot y min max values */ if (mbnavplot[iplot].type == PLOT_LONGITUDE || mbnavplot[iplot].type == PLOT_LATITUDE) strcpy(yformat, "%11.6f"); else strcpy(yformat, "%6.2f"); sprintf(string, yformat, mbnavplot[iplot].ymin); xg_justify(mbnavedit_xgid, string, &swidth,&sascent,&sdescent); xg_drawstring(mbnavedit_xgid, (int)(mbnavplot[iplot].ixmin - swidth - 0.03*margin_x), (int)(mbnavplot[iplot].iymin + 0.5*sascent), string, pixel_values[BLACK],XG_SOLIDLINE); sprintf(string, yformat, mbnavplot[iplot].ymax); xg_justify(mbnavedit_xgid, string, &swidth,&sascent,&sdescent); xg_drawstring(mbnavedit_xgid, (int)(mbnavplot[iplot].ixmin - swidth - 0.03*margin_x), (int)(mbnavplot[iplot].iymax + 0.5*sascent), string, pixel_values[BLACK],XG_SOLIDLINE); /* plot zero values */ if ((mbnavplot[iplot].ymax > 0.0 && mbnavplot[iplot].ymin < 0.0) || (mbnavplot[iplot].ymax < 0.0 && mbnavplot[iplot].ymin > 0.0)) { if (mbnavplot[iplot].type == PLOT_LONGITUDE || mbnavplot[iplot].type == PLOT_LATITUDE) strcpy(yformat, "%11.6f"); else strcpy(yformat, "%6.2f"); sprintf(string, yformat, 0.0); xg_justify(mbnavedit_xgid, string, &swidth,&sascent,&sdescent); iyzero = mbnavplot[iplot].iymin - mbnavplot[iplot].yscale*mbnavplot[iplot].ymin; xg_drawstring(mbnavedit_xgid, (int)(mbnavplot[iplot].ixmin - swidth - 0.03*margin_x), (int)(iyzero + 0.5*sascent), string, pixel_values[BLACK],XG_SOLIDLINE); xg_drawline(mbnavedit_xgid, mbnavplot[iplot].ixmin, iyzero, mbnavplot[iplot].ixmax, iyzero, pixel_values[BLACK],XG_DASHLINE); } /* plot bounding box */ xg_drawrectangle(mbnavedit_xgid, mbnavplot[iplot].ixmin, mbnavplot[iplot].iymax, mbnavplot[iplot].ixmax - mbnavplot[iplot].ixmin, mbnavplot[iplot].iymin - mbnavplot[iplot].iymax, pixel_values[BLACK],XG_SOLIDLINE); xg_drawrectangle(mbnavedit_xgid, mbnavplot[iplot].ixmin-1, mbnavplot[iplot].iymax-1, mbnavplot[iplot].ixmax - mbnavplot[iplot].ixmin + 2, mbnavplot[iplot].iymin - mbnavplot[iplot].iymax + 2, pixel_values[BLACK],XG_SOLIDLINE); /* now plot the data */ if (mbnavplot[iplot].type == PLOT_TINT) mbnavedit_plot_tint(iplot); else if (mbnavplot[iplot].type == PLOT_LONGITUDE) mbnavedit_plot_lon(iplot); else if (mbnavplot[iplot].type == PLOT_LATITUDE) mbnavedit_plot_lat(iplot); else if (mbnavplot[iplot].type == PLOT_SPEED) mbnavedit_plot_speed(iplot); else if (mbnavplot[iplot].type == PLOT_HEADING) mbnavedit_plot_heading(iplot); else if (mbnavplot[iplot].type == PLOT_DRAFT) mbnavedit_plot_draft(iplot); else if (mbnavplot[iplot].type == PLOT_ROLL) mbnavedit_plot_roll(iplot); else if (mbnavplot[iplot].type == PLOT_PITCH) mbnavedit_plot_pitch(iplot); else if (mbnavplot[iplot].type == PLOT_HEAVE) mbnavedit_plot_heave(iplot); } } /* set status */ if (nplot > 0) status = MB_SUCCESS; else status = MB_FAILURE; /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return (status); } /*--------------------------------------------------------------------*/ int mbnavedit_plot_tint(int iplot) { /* local variables */ char *function_name = "mbnavedit_plot_tint"; int status = MB_SUCCESS; int ixmin, ixmax, iymin, iymax; double xmin, xmax, ymin, ymax; double xscale, yscale; int tint_x1, tint_y1, tint_x2, tint_y2; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 iplot: %d\n",iplot); } /* get scaling values */ ixmin = mbnavplot[iplot].ixmin; ixmax = mbnavplot[iplot].ixmax; iymin = mbnavplot[iplot].iymin; iymax = mbnavplot[iplot].iymax; xmin = mbnavplot[iplot].xmin; xmax = mbnavplot[iplot].xmax; ymin = mbnavplot[iplot].ymin; ymax = mbnavplot[iplot].ymax; xscale = mbnavplot[iplot].xscale; yscale = mbnavplot[iplot].yscale; /* plot original expected time data */ if (plot_tint_org == MB_YES) { tint_x1 = ixmin + xscale*(ping[current_id].file_time_d - xmin); tint_y1 = iymin + yscale*(ping[current_id].tint_org - ymin); for (i=current_id+1;i= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return (status); } /*--------------------------------------------------------------------*/ int mbnavedit_plot_lon(int iplot) { /* local variables */ char *function_name = "mbnavedit_plot_lon"; int status = MB_SUCCESS; int ixmin, ixmax, iymin, iymax; double xmin, xmax, ymin, ymax; double xscale, yscale; int lon_x1, lon_y1, lon_x2, lon_y2; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 iplot: %d\n",iplot); } /* get scaling values */ ixmin = mbnavplot[iplot].ixmin; ixmax = mbnavplot[iplot].ixmax; iymin = mbnavplot[iplot].iymin; iymax = mbnavplot[iplot].iymax; xmin = mbnavplot[iplot].xmin; xmax = mbnavplot[iplot].xmax; ymin = mbnavplot[iplot].ymin; ymax = mbnavplot[iplot].ymax; xscale = mbnavplot[iplot].xscale; yscale = mbnavplot[iplot].yscale; /* plot original longitude data */ if (plot_lon_org == MB_YES) { lon_x1 = ixmin + xscale*(ping[current_id].file_time_d - xmin); lon_y1 = iymin + yscale*(ping[current_id].lon_org - ymin); for (i=current_id+1;i= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return (status); } /*--------------------------------------------------------------------*/ int mbnavedit_plot_lat(int iplot) { /* local variables */ char *function_name = "mbnavedit_plot_lat"; int status = MB_SUCCESS; int ixmin, ixmax, iymin, iymax; double xmin, xmax, ymin, ymax; double xscale, yscale; int lat_x1, lat_y1, lat_x2, lat_y2; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 iplot: %d\n",iplot); } /* get scaling values */ ixmin = mbnavplot[iplot].ixmin; ixmax = mbnavplot[iplot].ixmax; iymin = mbnavplot[iplot].iymin; iymax = mbnavplot[iplot].iymax; xmin = mbnavplot[iplot].xmin; xmax = mbnavplot[iplot].xmax; ymin = mbnavplot[iplot].ymin; ymax = mbnavplot[iplot].ymax; xscale = mbnavplot[iplot].xscale; yscale = mbnavplot[iplot].yscale; /* plot original latitude data */ if (plot_lat_org == MB_YES) { lat_x1 = ixmin + xscale*(ping[current_id].file_time_d - xmin); lat_y1 = iymin + yscale*(ping[current_id].lat_org - ymin); for (i=current_id+1;i= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return (status); } /*--------------------------------------------------------------------*/ int mbnavedit_plot_speed(int iplot) { /* local variables */ char *function_name = "mbnavedit_plot_speed"; int status = MB_SUCCESS; int ixmin, ixmax, iymin, iymax; double xmin, xmax, ymin, ymax; double xscale, yscale; int speed_x1, speed_y1, speed_x2, speed_y2; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 iplot: %d\n",iplot); } /* get scaling values */ ixmin = mbnavplot[iplot].ixmin; ixmax = mbnavplot[iplot].ixmax; iymin = mbnavplot[iplot].iymin; iymax = mbnavplot[iplot].iymax; xmin = mbnavplot[iplot].xmin; xmax = mbnavplot[iplot].xmax; ymin = mbnavplot[iplot].ymin; ymax = mbnavplot[iplot].ymax; xscale = mbnavplot[iplot].xscale; yscale = mbnavplot[iplot].yscale; /* plot original speed data */ if (plot_speed_org == MB_YES) { speed_x1 = ixmin + xscale*(ping[current_id].file_time_d - xmin); speed_y1 = iymin + yscale*(ping[current_id].speed - ymin); for (i=current_id+1;i= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return (status); } /*--------------------------------------------------------------------*/ int mbnavedit_plot_heading(int iplot) { /* local variables */ char *function_name = "mbnavedit_plot_heading"; int status = MB_SUCCESS; int ixmin, ixmax, iymin, iymax; double xmin, xmax, ymin, ymax; double xscale, yscale; int heading_x1, heading_y1, heading_x2, heading_y2; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 iplot: %d\n",iplot); } /* get scaling values */ ixmin = mbnavplot[iplot].ixmin; ixmax = mbnavplot[iplot].ixmax; iymin = mbnavplot[iplot].iymin; iymax = mbnavplot[iplot].iymax; xmin = mbnavplot[iplot].xmin; xmax = mbnavplot[iplot].xmax; ymin = mbnavplot[iplot].ymin; ymax = mbnavplot[iplot].ymax; xscale = mbnavplot[iplot].xscale; yscale = mbnavplot[iplot].yscale; /* plot original heading data */ if (plot_heading_org == MB_YES) { heading_x1 = ixmin + xscale*(ping[current_id].file_time_d - xmin); heading_y1 = iymin + yscale*(ping[current_id].heading - ymin); for (i=current_id+1;i= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return (status); } /*--------------------------------------------------------------------*/ int mbnavedit_plot_draft(int iplot) { /* local variables */ char *function_name = "mbnavedit_plot_draft"; int status = MB_SUCCESS; int ixmin, ixmax, iymin, iymax; double xmin, xmax, ymin, ymax; double xscale, yscale; int draft_x1, draft_y1, draft_x2, draft_y2; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 iplot: %d\n",iplot); } /* get scaling values */ ixmin = mbnavplot[iplot].ixmin; ixmax = mbnavplot[iplot].ixmax; iymin = mbnavplot[iplot].iymin; iymax = mbnavplot[iplot].iymax; xmin = mbnavplot[iplot].xmin; xmax = mbnavplot[iplot].xmax; ymin = mbnavplot[iplot].ymin; ymax = mbnavplot[iplot].ymax; xscale = mbnavplot[iplot].xscale; yscale = mbnavplot[iplot].yscale; /* plot original draft data */ if (plot_draft_org == MB_YES) { draft_x1 = ixmin + xscale*(ping[current_id].file_time_d - xmin); draft_y1 = iymin + yscale*(ping[current_id].draft - ymin); for (i=current_id+1;i= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return (status); } /*--------------------------------------------------------------------*/ int mbnavedit_plot_roll(int iplot) { /* local variables */ char *function_name = "mbnavedit_plot_roll"; int status = MB_SUCCESS; int ixmin, ixmax, iymin, iymax; double xmin, xmax, ymin, ymax; double xscale, yscale; int roll_x1, roll_y1, roll_x2, roll_y2; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 iplot: %d\n",iplot); } /* get scaling values */ ixmin = mbnavplot[iplot].ixmin; ixmax = mbnavplot[iplot].ixmax; iymin = mbnavplot[iplot].iymin; iymax = mbnavplot[iplot].iymax; xmin = mbnavplot[iplot].xmin; xmax = mbnavplot[iplot].xmax; ymin = mbnavplot[iplot].ymin; ymax = mbnavplot[iplot].ymax; xscale = mbnavplot[iplot].xscale; yscale = mbnavplot[iplot].yscale; /* plot roll data */ if (plot_roll == MB_YES) { roll_x1 = ixmin + xscale*(ping[current_id].file_time_d - xmin); roll_y1 = iymin + yscale*(ping[current_id].roll - ymin); for (i=current_id+1;i= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return (status); } /*--------------------------------------------------------------------*/ int mbnavedit_plot_pitch(int iplot) { /* local variables */ char *function_name = "mbnavedit_plot_pitch"; int status = MB_SUCCESS; int ixmin, ixmax, iymin, iymax; double xmin, xmax, ymin, ymax; double xscale, yscale; int pitch_x1, pitch_y1, pitch_x2, pitch_y2; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 iplot: %d\n",iplot); } /* get scaling values */ ixmin = mbnavplot[iplot].ixmin; ixmax = mbnavplot[iplot].ixmax; iymin = mbnavplot[iplot].iymin; iymax = mbnavplot[iplot].iymax; xmin = mbnavplot[iplot].xmin; xmax = mbnavplot[iplot].xmax; ymin = mbnavplot[iplot].ymin; ymax = mbnavplot[iplot].ymax; xscale = mbnavplot[iplot].xscale; yscale = mbnavplot[iplot].yscale; /* plot pitch data */ if (plot_pitch == MB_YES) { pitch_x1 = ixmin + xscale*(ping[current_id].file_time_d - xmin); pitch_y1 = iymin + yscale*(ping[current_id].pitch - ymin); for (i=current_id+1;i= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return (status); } /*--------------------------------------------------------------------*/ int mbnavedit_plot_heave(int iplot) { /* local variables */ char *function_name = "mbnavedit_plot_heave"; int status = MB_SUCCESS; int ixmin, ixmax, iymin, iymax; double xmin, xmax, ymin, ymax; double xscale, yscale; int heave_x1, heave_y1, heave_x2, heave_y2; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 iplot: %d\n",iplot); } /* get scaling values */ ixmin = mbnavplot[iplot].ixmin; ixmax = mbnavplot[iplot].ixmax; iymin = mbnavplot[iplot].iymin; iymax = mbnavplot[iplot].iymax; xmin = mbnavplot[iplot].xmin; xmax = mbnavplot[iplot].xmax; ymin = mbnavplot[iplot].ymin; ymax = mbnavplot[iplot].ymax; xscale = mbnavplot[iplot].xscale; yscale = mbnavplot[iplot].yscale; /* plot heave data */ if (plot_heave == MB_YES) { heave_x1 = ixmin + xscale*(ping[current_id].file_time_d - xmin); heave_y1 = iymin + yscale*(ping[current_id].heave - ymin); for (i=current_id+1;i= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return (status); } /*--------------------------------------------------------------------*/ int mbnavedit_plot_tint_value(int iplot, int iping) { /* local variables */ char *function_name = "mbnavedit_plot_tint_value"; int status = MB_SUCCESS; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 iplot: %d\n",iplot); fprintf(stderr,"dbg2 iping: %d\n",iping); } /* unplot basic expected time data value */ xg_drawrectangle(mbnavedit_xgid, ping[iping].tint_x-2, ping[iping].tint_y-2, 4, 4, pixel_values[WHITE],XG_SOLIDLINE); xg_fillrectangle(mbnavedit_xgid, ping[iping].tint_x-2, ping[iping].tint_y-2, 4, 4, pixel_values[WHITE],XG_SOLIDLINE); /* replot basic expected time data value */ if (ping[iping].tint_select == MB_YES) xg_drawrectangle(mbnavedit_xgid, ping[iping].tint_x-2, ping[iping].tint_y-2, 4, 4, pixel_values[RED],XG_SOLIDLINE); else if (ping[iping].tint != ping[iping].tint_org) xg_drawrectangle(mbnavedit_xgid, ping[iping].tint_x-2, ping[iping].tint_y-2, 4, 4, pixel_values[PURPLE],XG_SOLIDLINE); else xg_fillrectangle(mbnavedit_xgid, ping[iping].tint_x-2, ping[iping].tint_y-2, 4, 4, pixel_values[BLACK],XG_SOLIDLINE); /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return (status); } /*--------------------------------------------------------------------*/ int mbnavedit_plot_lon_value(int iplot, int iping) { /* local variables */ char *function_name = "mbnavedit_plot_lon_value"; int status = MB_SUCCESS; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 iplot: %d\n",iplot); fprintf(stderr,"dbg2 iping: %d\n",iping); } /* unplot basic lon data value */ xg_drawrectangle(mbnavedit_xgid, ping[iping].lon_x-2, ping[iping].lon_y-2, 4, 4, pixel_values[WHITE],XG_SOLIDLINE); xg_fillrectangle(mbnavedit_xgid, ping[iping].lon_x-2, ping[iping].lon_y-2, 4, 4, pixel_values[WHITE],XG_SOLIDLINE); /* replot basic lon data value */ if (ping[iping].lon_select == MB_YES) xg_drawrectangle(mbnavedit_xgid, ping[iping].lon_x-2, ping[iping].lon_y-2, 4, 4, pixel_values[RED],XG_SOLIDLINE); else if (ping[iping].lonlat_flag == MB_YES) xg_drawrectangle(mbnavedit_xgid, ping[iping].lon_x-2, ping[iping].lon_y-2, 4, 4, pixel_values[ORANGE],XG_SOLIDLINE); else if (ping[iping].lon != ping[iping].lon_org) xg_drawrectangle(mbnavedit_xgid, ping[iping].lon_x-2, ping[iping].lon_y-2, 4, 4, pixel_values[PURPLE],XG_SOLIDLINE); else xg_fillrectangle(mbnavedit_xgid, ping[iping].lon_x-2, ping[iping].lon_y-2, 4, 4, pixel_values[BLACK],XG_SOLIDLINE); /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return (status); } /*--------------------------------------------------------------------*/ int mbnavedit_plot_lat_value(int iplot, int iping) { /* local variables */ char *function_name = "mbnavedit_plot_lat_value"; int status = MB_SUCCESS; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 iplot: %d\n",iplot); fprintf(stderr,"dbg2 iping: %d\n",iping); } /* unplot basic lat data value */ xg_drawrectangle(mbnavedit_xgid, ping[iping].lat_x-2, ping[iping].lat_y-2, 4, 4, pixel_values[WHITE],XG_SOLIDLINE); xg_fillrectangle(mbnavedit_xgid, ping[iping].lat_x-2, ping[iping].lat_y-2, 4, 4, pixel_values[WHITE],XG_SOLIDLINE); /* replot basic lat data value */ if (ping[iping].lat_select == MB_YES) xg_drawrectangle(mbnavedit_xgid, ping[iping].lat_x-2, ping[iping].lat_y-2, 4, 4, pixel_values[RED],XG_SOLIDLINE); else if (ping[iping].lonlat_flag == MB_YES) xg_drawrectangle(mbnavedit_xgid, ping[iping].lat_x-2, ping[iping].lat_y-2, 4, 4, pixel_values[ORANGE],XG_SOLIDLINE); else if (ping[iping].lat != ping[iping].lat_org) xg_drawrectangle(mbnavedit_xgid, ping[iping].lat_x-2, ping[iping].lat_y-2, 4, 4, pixel_values[PURPLE],XG_SOLIDLINE); else xg_fillrectangle(mbnavedit_xgid, ping[iping].lat_x-2, ping[iping].lat_y-2, 4, 4, pixel_values[BLACK],XG_SOLIDLINE); /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return (status); } /*--------------------------------------------------------------------*/ int mbnavedit_plot_speed_value(int iplot, int iping) { /* local variables */ char *function_name = "mbnavedit_plot_speed_value"; int status = MB_SUCCESS; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 iplot: %d\n",iplot); fprintf(stderr,"dbg2 iping: %d\n",iping); } /* unplot basic speed data value */ xg_drawrectangle(mbnavedit_xgid, ping[iping].speed_x-2, ping[iping].speed_y-2, 4, 4, pixel_values[WHITE],XG_SOLIDLINE); xg_fillrectangle(mbnavedit_xgid, ping[iping].speed_x-2, ping[iping].speed_y-2, 4, 4, pixel_values[WHITE],XG_SOLIDLINE); /* replot basic speed data value */ if (ping[iping].speed_select == MB_YES) xg_drawrectangle(mbnavedit_xgid, ping[iping].speed_x-2, ping[iping].speed_y-2, 4, 4, pixel_values[RED],XG_SOLIDLINE); else if (ping[iping].speed != ping[iping].speed_org) xg_drawrectangle(mbnavedit_xgid, ping[iping].speed_x-2, ping[iping].speed_y-2, 4, 4, pixel_values[PURPLE],XG_SOLIDLINE); else xg_fillrectangle(mbnavedit_xgid, ping[iping].speed_x-2, ping[iping].speed_y-2, 4, 4, pixel_values[BLACK],XG_SOLIDLINE); /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return (status); } /*--------------------------------------------------------------------*/ int mbnavedit_plot_heading_value(int iplot, int iping) { /* local variables */ char *function_name = "mbnavedit_plot_heading_value"; int status = MB_SUCCESS; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 iplot: %d\n",iplot); fprintf(stderr,"dbg2 iping: %d\n",iping); } /* unplot basic heading data value */ xg_drawrectangle(mbnavedit_xgid, ping[iping].heading_x-2, ping[iping].heading_y-2, 4, 4, pixel_values[WHITE],XG_SOLIDLINE); xg_fillrectangle(mbnavedit_xgid, ping[iping].heading_x-2, ping[iping].heading_y-2, 4, 4, pixel_values[WHITE],XG_SOLIDLINE); /* replot basic heading data value */ if (ping[iping].heading_select == MB_YES) xg_drawrectangle(mbnavedit_xgid, ping[iping].heading_x-2, ping[iping].heading_y-2, 4, 4, pixel_values[RED],XG_SOLIDLINE); else if (ping[iping].heading != ping[iping].heading_org) xg_drawrectangle(mbnavedit_xgid, ping[iping].heading_x-2, ping[iping].heading_y-2, 4, 4, pixel_values[PURPLE],XG_SOLIDLINE); else xg_fillrectangle(mbnavedit_xgid, ping[iping].heading_x-2, ping[iping].heading_y-2, 4, 4, pixel_values[BLACK],XG_SOLIDLINE); /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return (status); } /*--------------------------------------------------------------------*/ int mbnavedit_plot_draft_value(int iplot, int iping) { /* local variables */ char *function_name = "mbnavedit_plot_draft_value"; int status = MB_SUCCESS; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n", function_name); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 iplot: %d\n",iplot); fprintf(stderr,"dbg2 iping: %d\n",iping); } /* unplot basic draft data value */ xg_drawrectangle(mbnavedit_xgid, ping[iping].draft_x-2, ping[iping].draft_y-2, 4, 4, pixel_values[WHITE],XG_SOLIDLINE); xg_fillrectangle(mbnavedit_xgid, ping[iping].draft_x-2, ping[iping].draft_y-2, 4, 4, pixel_values[WHITE],XG_SOLIDLINE); /* replot basic draft data value */ if (ping[iping].draft_select == MB_YES) xg_drawrectangle(mbnavedit_xgid, ping[iping].draft_x-2, ping[iping].draft_y-2, 4, 4, pixel_values[RED],XG_SOLIDLINE); else if (ping[iping].draft != ping[iping].draft_org) xg_drawrectangle(mbnavedit_xgid, ping[iping].draft_x-2, ping[iping].draft_y-2, 4, 4, pixel_values[PURPLE],XG_SOLIDLINE); else xg_fillrectangle(mbnavedit_xgid, ping[iping].draft_x-2, ping[iping].draft_y-2, 4, 4, pixel_values[BLACK],XG_SOLIDLINE); /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n", function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return */ return (status); } /*--------------------------------------------------------------------*/