.\" Man page generated from reStructuredText.
.
.TH "X2SYS_SOLVE" "1gmt" "Nov 05, 2016" "5.3.1" "GMT"
.SH NAME
x2sys_solve \- Determine least-squares systematic correction from crossovers
.
.nr rst2man-indent-level 0
.
.de1 rstReportMargin
\\$1 \\n[an-margin]
level \\n[rst2man-indent-level]
level margin: \\n[rst2man-indent\\n[rst2man-indent-level]]
-
\\n[rst2man-indent0]
\\n[rst2man-indent1]
\\n[rst2man-indent2]
..
.de1 INDENT
.\" .rstReportMargin pre:
. RS \\$1
. nr rst2man-indent\\n[rst2man-indent-level] \\n[an-margin]
. nr rst2man-indent-level +1
.\" .rstReportMargin post:
..
.de UNINDENT
. RE
.\" indent \\n[an-margin]
.\" old: \\n[rst2man-indent\\n[rst2man-indent-level]]
.nr rst2man-indent-level -1
.\" new: \\n[rst2man-indent\\n[rst2man-indent-level]]
.in \\n[rst2man-indent\\n[rst2man-indent-level]]u
..
.SH SYNOPSIS
.sp
\fBx2sys_solve\fP  \fB\-C\fP\fIcolumn\fP  \fB\-T\fP\fITAG\fP  \fB\-E\fP\fImode\fP
[ \fICOE_list.d\fP ]
[  \fB\-V\fP[\fIlevel\fP] ]
[  \fB\-W\fP\fB[u]\fP ]
[ \fB\-bi\fPbinary ]
[ \fB\-di\fPnodata ]
[ \fB\-x\fP[[\-]\fIn\fP] ]
.sp
\fBNote:\fP No space is allowed between the option flag and the associated arguments.
.SH DESCRIPTION
.sp
\fBx2sys_solve\fP will use the supplied crossover information to solve
for systematic corrections that can then be applied per track to improve
data quality. Several systematic corrections can be solved for using a
least\-squares approach. Note: Only one data column can be processed at
the time.
.SH REQUIRED ARGUMENTS
.INDENT 0.0
.TP
.B \fICOE_list.d\fP
Name of file with the required crossover columns as produced by
\fBx2sys_list\fP\&. NOTE: If \fB\-bi\fP is used
then the first two columns are expected to hold the integer track
IDs; otherwise we expect those columns to hold the text string names
of the two tracks. If no file is given we will read from \fIstdin\fP\&.
.UNINDENT
.INDENT 0.0
.TP
\fB\-T\fP\fITAG\fP
Specify the x2sys \fITAG\fP which tracks the attributes of this data type.
.UNINDENT
.INDENT 0.0
.TP
\fB\-C\fP\fIcolumn\fP
Specify which data column you want to process. Needed for proper
formatting of the output correction table and must match the same
option used in \fBx2sys_list\fP when preparing the input data.
.UNINDENT
.INDENT 0.0
.TP
\fB\-E\fP\fImode\fP
The correction type you wish to model. Choose among the following
functions f(\fIp\fP) , where \fIp\fP are the \fIm\fP
parameters per track that we will fit simultaneously using a least
squares approach:
.sp
\fBc\fP will fit f(\fIp\fP) = \fIa\fP (a constant offset);
records must contain track ID1, ID2, COE.
.sp
\fBd\fP will fit f(\fIp\fP) = \fIa\fP + \fIb\fP * \fId\fP (linear
drift; \fId\fP is distance; records must contain track ID1, ID2, d1,
d2, COE.
.sp
\fBg\fP will fit f(\fIp\fP) = \fIa\fP + \fIb\fP sin(y)^2
(1980\-1930 gravity correction); records must contain track ID1,
ID2, latitude y, COE.
.sp
\fBh\fP will fit f(\fIp\fP) = \fIa\fP + \fIb\fP cos(H) + \fIc\fP
cos(2H) + \fId\fP sin(H) + \fIe\fP sin(2H)
(magnetic heading correction); records must contain track ID1, ID2,
heading H, COE.
.sp
\fBs\fP will fit f(\fIp\fP) = \fIa\fP * z (a unit scale
correction); records must contain track ID1, ID2, z1, z2.
.sp
\fBt\fP will fit f(\fIp\fP) = \fIa\fP + \fIb\fP * (\fIt \- t0\fP)
(linear drift; \fIt0\fP is the start time of the track); records must
contain track ID1, ID2, t1\-t0, t2\-t0, COE.
.UNINDENT
.SH OPTIONAL ARGUMENTS
.INDENT 0.0
.TP
\fB\-V\fP[\fIlevel\fP] (more ...)
Select verbosity level [c].  
.UNINDENT
.INDENT 0.0
.TP
\fB\-W\fP
Means that each input records has an extra column with the composite
weight for each crossover record. These are used to obtain a
weighted least squares solution [no weights]. Append \fBu\fP to report
unweighted mean/std [Default, report weighted stats].
.UNINDENT
.INDENT 0.0
.TP
\fB\-bi\fP[\fIncols\fP][\fBt\fP] (more ...)
Select native binary input.  
.UNINDENT
.INDENT 0.0
.TP
\fB\-di\fP\fInodata\fP (more ...)
Replace input columns that equal \fInodata\fP with NaN.  
.UNINDENT
.INDENT 0.0
.TP
\fB\-x\fP[[\-]\fIn\fP] (more ...)
Limit number of cores used in multi\-threaded algorithms (OpenMP required).
.UNINDENT
.INDENT 0.0
.TP
\fB\-^\fP or just \fB\-\fP
Print a short message about the syntax of the command, then exits (NOTE: on Windows use just \fB\-\fP).
.TP
\fB\-+\fP or just \fB+\fP
Print an extensive usage (help) message, including the explanation of
any module\-specific option (but not the GMT common options), then exits.
.TP
\fB\-?\fP or no arguments
Print a complete usage (help) message, including the explanation of options, then exits.
.UNINDENT
.SH NOTES
.sp
Most of the model corrections in  \fB\-E\fP involve a constant offset.
Because crossovers are differences between values, any absolute level
will cancel out and hence the constant offsets we obtain are relative
to an undetermined absolute level.  To obtain a solvable solution we
add the constraint that the sum of all constant offsets equal zero.
If the tracks form clusters in which no tracks from one cluster cross
any track from another cluster then these are two independent data
sets and require they own constraint equation for their offsets.  We
determine the number of clusters and automatically add the required
constraint equations.  If you need a particular reference track to have
a particular offset (e.g., 0) then you can subtract the offset you
found from every track correction and add in the desired offset.
.SH EXAMPLES
.sp
To fit a simple bias offset to faa for all tracks under the MGD77 tag, try
.INDENT 0.0
.INDENT 3.5
.INDENT 0.0
.INDENT 3.5
.sp
.nf
.ft C
gmt x2sys_list COE_data.txt \-V \-TMGD77 \-Cfaa \-Fnc > faa_coe.txt
gmt x2sys_solve faa_coe.txt \-V \-TMGD77 \-Cfaa \-Ec > coe_table.txt
.ft P
.fi
.UNINDENT
.UNINDENT
.UNINDENT
.UNINDENT
.sp
To fit a faa linear drift with time instead, try
.INDENT 0.0
.INDENT 3.5
.INDENT 0.0
.INDENT 3.5
.sp
.nf
.ft C
gmt x2sys_list COE_data.txt \-V \-TMGD77 \-Cfaa \-FnTc > faa_coe.txt
gmt x2sys_solve faa_coe.txt \-V \-TMGD77 \-Cfaa \-Et > coe_table.txt
.ft P
.fi
.UNINDENT
.UNINDENT
.UNINDENT
.UNINDENT
.sp
To estimate heading corrections based on magnetic crossovers associated
with the tag MGD77 from the file COE_data.txt, try
.INDENT 0.0
.INDENT 3.5
.INDENT 0.0
.INDENT 3.5
.sp
.nf
.ft C
gmt x2sys_list COE_data.txt \-V \-TMGD77 \-Cmag \-Fnhc > mag_coe.txt
gmt x2sys_solve mag_coe.txt \-V \-TMGD77 \-Cmag \-Eh > coe_table.txt
.ft P
.fi
.UNINDENT
.UNINDENT
.UNINDENT
.UNINDENT
.sp
To estimate unit scale corrections based on bathymetry crossovers, try
.INDENT 0.0
.INDENT 3.5
.INDENT 0.0
.INDENT 3.5
.sp
.nf
.ft C
gmt x2sys_list COE_data.txt \-V \-TMGD77 \-Cdepth \-Fnz > depth_coe.txt
gmt x2sys_solve depth_coe.txt \-V \-TMGD77 \-Cdepth \-Es > coe_table.txt
.ft P
.fi
.UNINDENT
.UNINDENT
.UNINDENT
.UNINDENT
.SH SEE ALSO
.sp
x2sys_binlist,
x2sys_cross,
x2sys_datalist,
x2sys_get,
x2sys_init,
x2sys_list,
x2sys_put,
x2sys_report
.SH COPYRIGHT
2016, P. Wessel, W. H. F. Smith, R. Scharroo, J. Luis, and F. Wobbe
.\" Generated by docutils manpage writer.
.