STILTS-TMATCHN(1)

πŸ’– STILTS-TMATCHN πŸ’–

okay bestie, welcome to the manual page for stilts-tmatchn β€” part of the Stilts commands family. sit down, grab a lil snack πŸ“, we're crossmatching MULTIPLE tables and it's gonna be cute.


✨ NAME ✨

stilts-tmatchn β€” Crossmatches multiple tables using flexible criteria πŸ¦‹


✨ SYNOPSIS ✨

here's the whole vibe in one go, lowkey iconic:

stilts tmatchn [nin=<count>] [ifmtN=<in-format>] [inN=<tableN>]
               [icmdN=<cmds>] [ocmd=<cmds>]
               [omode=out|meta|stats|count|checksum|cgi|discard|topcat|samp|tosql|gui]
               [out=<out-table>] [ofmt=<out-format>]
               [multimode=pairs|group] [iref=<table-index>]
               [matcher=<matcher-name>] [params=<match-params>]
               [tuning=<tuning-params>] [valuesN=<expr-list>]
               [joinN=default|match|nomatch|always]
               [fixcols=none|dups|all] [suffixN=<label>]
               [progress=none|log|time|profile]
               [runner=parallel|parallel<n>|parallel-
               all|sequential|classic|partest]

✨ DESCRIPTION ✨

so tmatchn does efficient and flexible crossmatching between multiple tables. it can match rows based on where they sit relative to each other in the sky πŸ›°οΈ, or honestly using a whole bunch of other criteria β€” like separation in some isotropic or anisotropic Cartesian space, identity of a key value, or literally some combo of all of those. the full range of match criteria is discussed in SUN/256 (iykyk, go read it).

now, since the match criteria are what define a "match" between two objects, it's not immediately obvious what a multi-table match even means. and that's real. so the command can work in one of two distinct modes, controlled by the multimode parameter:

  • pairs mode β€” one table (by default the first input table) is designated the reference table, and pair matches between each of the other tables and that one are identified.
  • group mode β€” groups of objects from all the input tables are identified, as discussed in SUN/256.

currently, in both cases an output matched row cannot contain more than one object from each input table. 😭 options for output of multiple rows per input table per match may be forthcoming in future releases if there is demand β€” so like, speak up bestie.

one more thing: tmatchn is intended for use with more than two input tables. for 1- and 2-table crossmatching respectively, go see tmatch1 and tmatch2 instead. πŸ’…


✨ OPTIONS ✨

okay here's the whole menu, no skips, no skimming:

  • nin=<count> β€” The number of input tables for this task. For each of the input tables N there will be associated parameters ifmtN, inN and icmdN. ✨
  • ifmtN=<in-format> β€” Specifies the format of input table #N as specified by parameter inN. The known formats are listed in SUN/256. this flag is your bestie when you already know what format your table is in. If it has the special value (auto) β€” which is the default β€” then an attempt will be made to detect the format of the table automatically. This cannot always be done correctly though 😬, in which case the program will exit with an error explaining which formats were attempted. This parameter is ignored for scheme-specified tables.
  • inN=<tableN> β€” The location of input table #N. This may take one of the following forms:
    • A filename.
    • A URL. 🌐
    • The special value "-", meaning standard input. In this case the input format must be given explicitly using the ifmtN parameter. Note that not all formats can be streamed in this way. 🫠
    • A scheme specification of the form :<scheme-name>:<scheme-args>.
    • A system command line with either a "<" character at the start, or a "|" character at the end ("<syscmd" or "syscmd|"). This executes the given pipeline and reads from its standard output. This will probably only work on unix-like systems.

In any case, compressed data in one of the supported compression formats (gzip, Unix compress or bzip2) will be decompressed transparently. πŸŽ€

  • icmdN=<cmds> β€” Specifies processing to be performed on input table #N as specified by parameter inN, before any other processing has taken place. The value of this parameter is one or more of the filter commands described in SUN/256. If more than one is given, they must be separated by semicolon characters (";"). This parameter can be repeated multiple times on the same command line to build up a list of processing steps. The sequence of commands given in this way defines the processing pipeline which is performed on the table. πŸ’«

Commands may alternatively be supplied in an external file, by using the indirection character '@'. Thus a value of "@filename" causes the file filename to be read for a list of filter commands to execute. The commands in the file may be separated by newline characters and/or semicolons, and lines which are blank or which start with a '#' character are ignored. A backslash character \ at the end of a line joins it with the following line.

  • ocmd=<cmds> β€” Specifies processing to be performed on the output table, after all other processing has taken place. The value of this parameter is one or more of the filter commands described in SUN/256. If more than one is given, they must be separated by semicolon characters (";"). This parameter can be repeated multiple times on the same command line to build up a list of processing steps. The sequence of commands given in this way defines the processing pipeline which is performed on the table. 🌈

Commands may alternatively be supplied in an external file, by using the indirection character '@'. Thus a value of "@filename" causes the file filename to be read for a list of filter commands to execute. The commands in the file may be separated by newline characters and/or semicolons, and lines which are blank or which start with a '#' character are ignored. A backslash character \ at the end of a line joins it with the following line.

  • omode=out|meta|stats|count|checksum|cgi|discard|topcat|samp|tosql|gui β€” The mode in which the result table will be output. The default mode is out, which means the result will be written as a new table to disk or elsewhere, as determined by the out and ofmt parameters. But there are other possibilities πŸ‘€, which correspond to uses to which a table can be put other than outputting it β€” like displaying metadata, calculating statistics, or populating a table in an SQL database. For some values of this parameter, additional parameters (<mode-args>) are required to determine the exact behaviour.

Possible values are:

  • out
  • meta
  • stats
  • count
  • checksum
  • cgi
  • discard
  • topcat
  • samp
  • tosql
  • gui

Use the help=omode flag or see SUN/256 for more information. πŸ’…

  • out=<out-table> β€” The location of the output table. This is usually a filename to write to. If it is equal to the special value "-" (the default) the output table will be written to standard output.
This parameter must only be given if omode has its default value of "out". ⚠️
  • ofmt=<out-format> β€” Specifies the format in which the output table will be written (one of the ones in SUN/256 β€” matching is case-insensitive and you can use just the first few letters, cute right?). If it has the special value "(auto)" (the default), then the output filename will be examined to try to guess what sort of file is required, usually by looking at the extension. If it's not obvious from the filename what output format is intended, an error will result. 😭
This parameter must only be given if omode has its default value of "out". ⚠️
  • multimode=pairs|group β€” Defines what is meant by a multi-table match. There are two possibilities:
    • pairs: Each output row corresponds to a single row of the reference table (see parameter iref) and contains entries from other tables which are pair matches to that. If a reference table row matches multiple rows from one of the other tables, only the best one is included.
    • group: Each output row corresponds to a group of entries from the input tables which are mutually linked by pair matches between them. This means that although you can get from any entry to any other entry via one or more pair matches, there is no guarantee that any entry is a pair match with any other entry. No table has privileged status in this case. If there are multiple entries from a given table in the match group, an arbitrary one is chosen for inclusion (there is no unique way to select the best). See SUN/256 for more discussion. 🧚

Note that which rows actually appear in the output is also influenced by the joinN parameter.

  • iref=<table-index> β€” If multimode=pairs this parameter gives the index of the table in the input table list which is to serve as the reference table (the one which must be matched by other tables). Ignored in other modes.

Row ordering in the output table is usually tidiest if the default setting of 1 is used (i.e. if the first input table is used as the reference table). ✨

  • matcher=<matcher-name> β€” Defines the nature of the matching that will be performed. Depending on the name supplied, this may be positional matching using celestial or Cartesian coordinates, exact matching on the value of a string column, or other things. A list and explanation of the available matching algorithms is given in SUN/256. The value supplied for this parameter determines the meanings of the values required by the params, values* and tuning parameter(s). πŸ”₯
  • params=<match-params> β€” Determines the parameters of this match. This is typically one or more tolerances such as error radii. It may contain zero or more values; the values that are required depend on the match type selected by the matcher parameter. If it contains multiple values, they must be separated by spaces; values which contain a space can be 'quoted' or "quoted". 🫢
  • tuning=<tuning-params> β€” Tuning values for the matching process, if appropriate. It can hold zero or more values; which values are permitted depends on the match type selected by the matcher parameter. If it contains multiple values, they must be separated by spaces; values that contain a space can be 'quoted' or "quoted". If you don't supply this optional parameter, sensible defaults will be chosen. Slay, defaults are your bestie. ✨
  • valuesN=<expr-list> β€” Defines the values from table N that are used to decide whether a match has occurred. These are typically coordinate values like RA and Dec, and maybe some per-row error values too, though exactly which values you need is determined by the kind of match from matcher. Depending on the kind of match, the number and type of required values will be different. Separate multiple values with whitespace; if whitespace happens inside a single value it must be 'quoted' or "quoted". Elements of the expression list are commonly just column names, but they can be algebraic expressions calculated from zero or more columns, as explained in SUN/256. πŸ’…
  • joinN=default|match|nomatch|always β€” Decides which rows from input table N make it into the output table. The matching algorithm figures out which rows in each input table correspond to which rows in the other input tables, and this parameter says what to do with that info.

The default vibe: a row appears in the output table if it represents a match of rows from two or more of the input tables. You can change that per input table with the non-default choices below:

  • match β€” Rows are included only if they contain an entry from input table N.
  • nomatch β€” Rows are included only if they do not contain an entry from input table N.
  • always β€” Rows are included if they contain an entry from input table N (overrides any match and nomatch settings of other tables).
  • default β€” Input table N has no special effect on whether rows are included.
  • fixcols=none|dups|all β€” Controls how input columns are renamed before use in the output table. The choices are:
    • none β€” columns are not renamed
    • dups β€” columns that would otherwise have duplicate names in the output will be renamed to show which table they came from
    • all β€” all columns will be renamed to show which table they came from

If columns are renamed, the new names are determined by the suffix* parameters. πŸ‘‘

  • suffixN=<label> β€” If fixcols is set so input columns are renamed for insertion into the output table, this parameter decides how the renaming happens. It gives a suffix that is appended to all renamed columns from table N. Cute. πŸŽ€
  • progress=none|log|time|profile β€” Decides whether info about the match progress gets output to the standard error stream while it runs. For lengthy matches this is a useful reassurance and can give guidance on how much longer it'll take. It's also useful as a performance diagnostic.

The options:

  • none β€” no progress is shown
  • log β€” progress information is shown
  • time β€” progress information and some time profiling information is shown
  • profile β€” progress information and limited time/memory profiling information are shown
  • runner=parallel|parallel<n>|parallel-all|sequential|classic|partest β€” Selects the threading implementation. Current options:
    • parallel β€” uses multithreaded implementation for large tables, with default parallelism, which is the smaller of 6 and the number of available processors
    • parallel<n> β€” uses multithreaded implementation for large tables, with parallelism given by the supplied value <n>
    • parallel-all β€” uses multithreaded implementation for large tables, with a parallelism given by the number of available processors
    • sequential β€” uses multithreaded implementation but with only a single thread
    • classic β€” uses legacy sequential implementation
    • partest β€” uses multithreaded implementation even when tables are small

The parallel* options should normally run faster than sequential or classic (which are provided mainly for testing purposes), at least for large matches and when multiple processing cores are available.

The default value "parallel" is currently limited to a parallelism of 6 since larger values yield diminishing returns given that some parts of the matching algorithms run sequentially (Amdahl's Law), and using too many threads can sometimes end up doing more work or impacting on other operations on the same machine. But you can experiment with other concurrencies, e.g. "parallel16" to run on 16 cores (if available) or "parallel-all" to run on all available cores. 🌈

The value of this parameter should make no difference to the matching results. If you notice any discrepancies please report them. No cap, report the weirdness. πŸ‘€

πŸ‘€ SEE ALSO πŸ‘€

stilts(1)

If the package stilts-doc is installed, the full documentation SUN/256 is available in HTML format: file:///usr/share/doc/stilts/sun256/index.html

πŸ’– VERSION πŸ’–

STILTS version 3.5.4-debian

This is the Debian version of Stilts, which lack the support of some file formats and network protocols. For differences see file:///usr/share/doc/stilts/README.Debian

✨ AUTHOR ✨

Mark Taylor (Bristol University)


And that's a wrap, bestie β€” this page is STILTS-TMATCHN(1), dated Mar 2017, and it's for STILTS version 3.5.4-debian. Go match those tables and stay iconic. πŸ’…πŸ¦‹πŸŒΈ

πŸ“¦ where this page comes from