/*-------------------------------------------------------------------- * The MB-system: mbsys_hdcs.h 3/16/99 * $Id: mbsys_hdcs.h 1917 2012-01-10 19:25:33Z caress $ * * Copyright (c) 1999-2012 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. *--------------------------------------------------------------------*/ /* * mbsys_hdcs.h defines the data structure used by MBIO functions * to store swath sonar data in the UNB OMG HDCS format: * MBF_OMGHDCSJ : MBIO ID 151 * * Author: D. W. Caress * Date: March 16, 1999 * * $Log: mbsys_hdcs.h,v $ * Revision 5.6 2008/03/14 18:33:03 caress * Updated support for JHC format 151. * * Revision 5.5 2005/11/05 00:48:04 caress * Programs changed to register arrays through mb_register_array() rather than allocating the memory directly with mb_realloc() or mb_malloc(). * * Revision 5.4 2003/04/17 21:05:23 caress * Release 5.0.beta30 * * Revision 5.3 2002/09/18 23:32:59 caress * Release 5.0.beta23 * * Revision 5.2 2001/07/20 00:32:54 caress * Release 5.0.beta03 * * Revision 5.1 2001/01/22 07:43:34 caress * Version 5.0.beta01 * * Revision 5.0 2000/12/01 22:48:41 caress * First cut at Version 5.0. * * Revision 4.1 2000/09/30 06:31:19 caress * Snapshot for Dale. * * Revision 4.0 1999/03/31 18:29:20 caress * MB-System 4.6beta7 * * Revision 1.1 1999/03/31 18:11:35 caress * Initial revision * * */ /* * Notes on the MBSYS_HDCS data structure: * 1. The OMG-HDCS format is a collection of similar data * formats used by the seafloor mapping software * developed by Dr. John Hughes Clarke (Ocean Mapping * Group of the University of New Brunswick). A variety * of sonars are supported in OMG-HDCS. * 2. OMG-HDCS files all begin with a summary header that * specifies the format version, the data source * (type of sonar), the number of records, and the * minimum and maximum values of position and data * values. * 3. The summary header is followed by a set of uniformly * sized data records. The data record size is determined * by the format version and data source. Ecch record is * is divided into a profile (header) and an array of * beam structures. * 4. Sidescan imagery can be stored in parallel files in * the same directory as the primary bathymetry and * amplitude files. The sidescan files have the suffix * ".ss_data" added to the end of the primary file's name. * 5. Comment records are encoded in MBIO by setting the first * eight bytes of the data record header (profile) to '#' * values. The profile part of the comment record will have * the same size as expected for data from the tool type * listed in the summary. However, the data section for * comment records will always be 256 bytes long regardless * of the tool type. This mechanism is supported only by * MB-System. * */ /* defines sizes and maximums */ #define MBSYS_HDCS_SUMMARY_SIZE 96 #define MBSYS_HDCS_MAX_COMMENT 252 #define MBSYS_HDCS_MAX_BEAMS 1440 #define MBSYS_HDCS_MAX_PIXELS 1024 /* define tools (sonars) supported by OMG_HDCS */ #define MBSYS_HDCS_NUM_TOOLS 46 /* as of June 2010 */ #define MBSYS_HDCS_None -1 #define MBSYS_HDCS_SingleBeam 0 #define MBSYS_HDCS_ELAC_BottomChart 1 #define MBSYS_HDCS_EM12_dual 2 #define MBSYS_HDCS_EM100_depth 3 #define MBSYS_HDCS_FanSweep 4 #define MBSYS_HDCS_SeaBeam 5 #define MBSYS_HDCS_EM3000 6 #define MBSYS_HDCS_ROSS_Profiler 7 #define MBSYS_HDCS_EM12_single 8 #define MBSYS_HDCS_EM100_depth_ss 9 #define MBSYS_HDCS_EM1000 10 #define MBSYS_HDCS_LADS_2ndary 11 #define MBSYS_HDCS_EM3000D 12 #define MBSYS_HDCS_SB2100 13 #define MBSYS_HDCS_ISIS_Submetrix 14 #define MBSYS_HDCS_EM1000_ampl 15 #define MBSYS_HDCS_SB2K 16 #define MBSYS_HDCS_Seabat9001 17 #define MBSYS_HDCS_FanSweep_10A 18 #define MBSYS_HDCS_FanSweep_20 19 #define MBSYS_HDCS_ISIS_SWA 20 #define MBSYS_HDCS_SeaBeam_1180_MkII 21 #define MBSYS_HDCS_SeaBat_8101 22 #define MBSYS_HDCS_EM300 23 #define MBSYS_HDCS_EM121A 24 #define MBSYS_HDCS_SM2000 25 #define MBSYS_HDCS_HydroSweep_MD2 26 #define MBSYS_HDCS_EM1002 27 #define MBSYS_HDCS_Humminbird 28 #define MBSYS_HDCS_Knudsen_320 29 #define MBSYS_HDCS_EM120 30 #define MBSYS_HDCS_SeaBat_8125 31 #define MBSYS_HDCS_SeaBat_8111 32 #define MBSYS_HDCS_SeaBat_8150 33 #define MBSYS_HDCS_EM3002 34 #define MBSYS_HDCS_Optech_Laser 35 #define MBSYS_HDCS_EM710 36 #define MBSYS_HDCS_EM3002D 37 #define MBSYS_HDCS_SeaBat_8160 38 #define MBSYS_HDCS_SEA_SwathPlus 39 #define MBSYS_HDCS_EM122 40 #define MBSYS_HDCS_EM302 41 #define MBSYS_HDCS_SeaBat_7125 42 #define MBSYS_HDCS_R2Sonic_2024 43 #define MBSYS_HDCS_SeaBat_7150 44 #define MBSYS_HDCS_OMG_GLORIA 45 #define MBSYS_HDCS_COMMENT 999 static char *mbsys_hdcs_tool_names[MBSYS_HDCS_NUM_TOOLS] = { "Single Beam Echo-Sounder", "HoneyWell Elac BottomChart Mk I", "Simrad EM-12 (dual system)", "Simrad EM-100 (depths only)", "Krupp-Atlas FanSweep 10", "General Instruments SeaBeam (Classic)", "Simrad - EM3000S", "ROSS Sweep - MV Profiler DPW", "Simrad EM-12 (single system)", "Simrad EM-100 (depths and amplitudes)", "Simrad EM-1000", "RAN -- LADS (secondary format) ", "Simrad - EM3000 dual head", "SeaBeam 2100 series", "ISIS Submetrix 100/ 2000", "EM1000 with reflectivities", "Seabeam 2000 ", "Reson Seabat 9001 ", "STN-Atlas FanSweep 10A", "STN-Atlas FanSweep 20", "ISIS Submetrix 100SWA", "SeaBeam 1180 MkII", "Reson 8101", "Simrad EM300", "Simrad EM121A", "Simrad SM2000", "HydroSweep MD2", "Simrad EM1002", "Hummin'bird 3D", "Knudsen 320", "Kongsberg Simrad EM120", "Reson 8125", "Reson 8111", "Reson 8150", "Kongsberg - EM3002", "Optech - Laser", "Kongsberg - EM710", "Kongsberg - EM3002D", "Reson 8160", "SEA SwathPlus - Dual sided", "Kongsberg - EM122", "Kongsberg - EM302", "Reson 7125", "R2Sonic 2024", "Reson 7150", "OMG GLORIA" }; /* define beam record structure */ struct mbsys_hdcs_beam_struct { int status; /* status is either OK (0) or bad (other) */ mb_u_char scaling_factor; /* V4 */ int observedDepth; /* Depth (mm) */ int acrossTrack; /* Across track position of depth (mm) */ int alongTrack; /* Along track position of depth (mm) */ int latOffset; /* Latitude offset wrt. profile */ int longOffset; /* Longitude offset wrt. profile */ int processedDepth; /* Depth (mm) */ int timeOffset; int depthAccuracy; /* Depth accuracy (mm) */ mb_u_char reflectivity; char Q_factor; /* phase or amplitude detection */ char beam_no; char freq; /* 12.7, 13.0, 13.3, 95.0, Smii, GLORIA */ char calibratedBackscatter; /* effects of power/TVG and atten. removed*/ char mindB; char maxdB; mb_u_char pseudoAngleIndependentBackscatter; /* corrected for mean angular dependence for geological visualisation */ int range; /* other option on EM 12 */ int no_samples; int offset; int centre_no; char sample_unit; /* whether in time or distance */ char sample_interval; /* seconds or metres */ char dummy[2]; mb_u_char samp_win_length; short beam_depress_angle; unsigned short beam_heading_angle; /* V4 */ unsigned short other_range; signed short Tx_steer; signed short Rc_steer; mb_u_char TxSector; unsigned int timestampOffset; /* really is a 64 bit integer, trying to compress */ /* would'nt even need if didn't have to relate wavefile by this number */ unsigned short no_RAMAN; unsigned short no_IR; unsigned short no_GAPD; unsigned short no_PMT; mb_u_char prim_depth_conf; mb_u_char seco_depth_conf; signed short scan_azimuth; /* 100ths of degree */ unsigned short nadir_angle; /* 100ths of degree */ /* always dynamically compressed for v4 using scaling factor */ signed int secondaryDepth; /* Depth (mm) remember can be +ve or -ve */ signed short wave_height; signed int opaqueDepth_PMT; /* Depth (mm) */ signed int extinctionDepth_PMT;/* Depth (mm) */ signed int pimDepth_PMT; /* Depth (mm) */ signed int opaqueDepth_GAPD; /* Depth (mm) */ signed int extinctionDepth_GAPD;/* Depth (mm) */ signed int pimDepth_GAPD; /* Depth (mm) */ float twtt; unsigned int snippet_first_sample; unsigned int snippet_last_sample; float intensity; }; /* specific FOR ATLAS SAPI data */ struct mbsys_omghdcs_profile_subparams_struct { unsigned short txBeamIndex; unsigned short txLevel; short txBeamAngle; unsigned short txPulseLength; unsigned int ss_offset; unsigned short no_skipped_ss; unsigned short no_acquired_ss; unsigned short ss_sample_interval; unsigned short bscatClass; unsigned short nrActualGainSets; short rxGup; short rxGain; short ar; unsigned short rxtime[20]; short rxgain[20]; }; /* structure to hold everything */ struct mbsys_hdcs_struct { int kind; int read_summary; int profile_size; int num_beam; int beam_size; int data_size; int image_size; int sensorNumber; /* 1 = depth file */ int subFileID; /* 1 = data (as opposed to index) */ int fileVersion; /* 1 = original format */ /* 2 = packed format for EM1000 and others */ /* 3 = packed format for EM300 */ int toolType; /* Tool Type implies Profile Record Size and Depth Record Size and Image Record Size */ int numProfiles; /* # of profiles in the file */ int numDepths_sum; /* # of depths in the file */ int timeScale; /* time scale (# of uSec. units) */ int refTime; /* Reference time (100 sec. units) */ int minTime; /* Minimum time (offset wrt. ref.) */ int maxTime; /* Maximum time (offset wrt. ref.) */ int positionType; /* Geographic(1)/ UTM(2) */ int positionScale; /* Position scale (# of nRad. units) */ int refLat; /* Reference latitude (100 nRadians) */ int minLat; /* Minimum latitude (offset wrt. ref.) */ int maxLat; /* Maximum latitude (offset wrt. ref.) */ int refLong; /* Reference longitude (100 nRadians) */ int minLong; /* Minimum longitude (offset wrt. ref.) */ int maxLong; /* Maximum longitude (offset wrt. ref.) */ int minObsDepth; /* Minimum depth (mm) */ int maxObsDepth; /* Maximum depth (mm) */ int minProcDepth; /* Minimum depth (mm) */ int maxProcDepth; /* Maximum depth (mm) */ int status_sum; /* status not actually used at all .... */ /* V4 */ int totalProfileBytes; int Profile_BitsDefining[20]; int totalBeamBytes; int Beam_BitsDefining[20]; /* End V4 Summary */ int status_pro; /* status is either OK (0) or no nav (1) or unwanted for gridding (2) or comment record (999) (MB-System only) */ int numDepths_pro; /* Number of depths in profile */ int numSamples; /* Number of sidescan samples in parallel file */ int timeOffset; /* Time offset wrt. header */ int vesselLatOffset; /* Latitude offset wrt. header */ int vesselLongOffset;/* Longitude offset wrt. header */ int vesselHeading; /* Heading (100 nRadians) */ int vesselHeave; /* Heave (mm) */ int vesselPitch; /* Vessel pitch (100 nRadians) */ int vesselRoll; /* Vessel roll (100 nRadians) */ int tide; /* Tide (mm) */ int vesselVelocity; /* Vessel Velocity (mm/s) note - transducer pitch is generally tucked into the vel field */ char power; char TVG; char attenuation; char edflag; int soundVelocity; /* mm/s */ int lengthImageDataField; int pingNo; char mode; char Q_factor; char pulseLength; /* centisecs*/ mb_u_char unassigned; unsigned short td_sound_speed; unsigned short samp_rate; mb_u_char z_res_cm; mb_u_char xy_res_cm; mb_u_char ssp_source; mb_u_char filter_ID; unsigned short absorp_coeff; unsigned short tx_pulse_len; unsigned short tx_beam_width; unsigned short max_swath_width; mb_u_char tx_power_reduction; mb_u_char rx_beam_width; mb_u_char rx_bandwidth; mb_u_char rx_gain_reduction; mb_u_char tvg_crossover; mb_u_char beam_spacing; mb_u_char coverage_sector; mb_u_char yaw_stab_mode; /* V4 */ struct mbsys_omghdcs_profile_subparams_struct params[2]; int transducerDepth; /* transducer or towfish depth */ int transducerPitch; /* Transducer pitch (100 nRadians) */ int transducerRoll; /* Transducer roll (100 nRadians) */ /* enough for dyn. stab transducer */ int transducerHeading; /* Transducer pitch (100 nRadians) */ int transducerLatOffset; /* Latitude offset wrt. vessel */ int transducerLongOffset; /* Longitude offset wrt. vessel */ int transducerSlantRange; /*slantRange(mm) wrt. vessel (cable out) */ int transducerAcross; /* horizontal Range (mm) wrt. vessel */ int transducerAlong; /* horizontal Range (mm) wrt. vessel */ int transducerBearing; /* Bearing (100nRads) wrt. vessel */ short longperiod_heaveCorrection; short dynamic_draftCorrection; short deepdraftoffset_in_metres; short draft_at_Tx; short alternateRoll; short alternatePitch; short alternateHeave; short alternateHeading; short standaloneHeading; short RTK_at_RP; /* in cm units so that can support +/- 320m. */ short Lowpass_RTK_at_RP; /* in cm units so that can support +/- 320m. */ short WLZ; unsigned short samp_rate_SecondHead; signed int clock_drift_millis; unsigned int watercol_offset; unsigned int watercol_size; unsigned int watercol_offset_2nd; unsigned int watercol_size_2nd; unsigned short range_to_normal_incidence; unsigned int laser_timestampRef; unsigned int tx_sector_offset; unsigned short num_tx_sectors; unsigned int sonar_settings_offset; unsigned int ping_number; unsigned short multi_ping_sequence; unsigned int num_beams; /* Which is different than numDepths...*/ /* (which is usually the number of possible depths and not ACTUAL depths */ /* e.g. high-density vs. low-density mode in KM systems) */ unsigned char layer_compensation_flag; float bs_beam_position; unsigned int bs_control_flags; unsigned short bs_num_beams_per_side; unsigned short bs_current_beam_number; unsigned char bs_sample_descriptor; unsigned int snippet_sample_descriptor; /* End V4 Profile */ struct mbsys_hdcs_beam_struct *beams; mb_s_char *ss_raw; int pixel_size; /* processed sidescan pixel size in mm */ int pixels_ss; /* number of processed sidescan pixels stored */ short ss_proc[MBSYS_HDCS_MAX_PIXELS]; /* the processed sidescan ordered port to starboard */ int ssalongtrack[MBSYS_HDCS_MAX_PIXELS]; /* the processed sidescan alongtrack distances in mm */ char comment[MBSYS_HDCS_MAX_COMMENT]; }; /* system specific function prototypes */ int mbsys_hdcs_alloc(int verbose, void *mbio_ptr, void **store_ptr, int *error); int mbsys_hdcs_deall(int verbose, void *mbio_ptr, void **store_ptr, int *error); int mbsys_hdcs_dimensions(int verbose, void *mbio_ptr, void *store_ptr, int *kind, int *nbath, int *namp, int *nss, int *error); int mbsys_hdcs_extract(int verbose, void *mbio_ptr, void *store_ptr, int *kind, int time_i[7], double *time_d, double *navlon, double *navlat, double *speed, double *heading, int *nbath, int *namp, int *nss, char *beamflag, double *bath, double *amp, double *bathacrosstrack, double *bathalongtrack, double *ss, double *ssacrosstrack, double *ssalongtrack, char *comment, int *error); int mbsys_hdcs_insert(int verbose, void *mbio_ptr, void *store_ptr, int kind, int time_i[7], double time_d, double navlon, double navlat, double speed, double heading, int nbath, int namp, int nss, char *beamflag, double *bath, double *amp, double *bathacrosstrack, double *bathalongtrack, double *ss, double *ssacrosstrack, double *ssalongtrack, char *comment, int *error); int mbsys_hdcs_ttimes(int verbose, void *mbio_ptr, void *store_ptr, int *kind, int *nbeams, double *ttimes, double *angles, double *angles_forward, double *angles_null, double *heave, double *alongtrack_offset, double *draft, double *ssv, int *error); int mbsys_hdcs_detects(int verbose, void *mbio_ptr, void *store_ptr, int *kind, int *nbeams, int *detects, int *error); int mbsys_hdcs_extract_altitude(int verbose, void *mbio_ptr, void *store_ptr, int *kind, double *transducer_depth, double *altitude, int *error); int mbsys_hdcs_insert_altitude(int verbose, void *mbio_ptr, void *store_ptr, double transducer_depth, double altitude, int *error); int mbsys_hdcs_extract_nav(int verbose, void *mbio_ptr, void *store_ptr, int *kind, int time_i[7], double *time_d, double *navlon, double *navlat, double *speed, double *heading, double *draft, double *roll, double *pitch, double *heave, int *error); int mbsys_hdcs_insert_nav(int verbose, void *mbio_ptr, void *store_ptr, int time_i[7], double time_d, double navlon, double navlat, double speed, double heading, double draft, double roll, double pitch, double heave, int *error); int mbsys_hdcs_copy(int verbose, void *mbio_ptr, void *store_ptr, void *copy_ptr, int *error);