1sphinx.addnodesdocument)}( rawsourcechildren](docutils.nodessubstitution_definition)}(h.. |rel| replace:: 10.71h]h Text10.71}parenth sba attributes}(ids]classes]names]reladupnames]backrefs]utagnameh source lineKhh _documenthubh )}(h.. |date| replace:: |today| h]h2024 136}(hh)h(hh%Nh'Nubah}(h]h]h]dateah ]h"]uh$h h%h&h'Khhh(hubh target)}(h.. _intro_prod_editing:h]h}(h]h]h]h ]h"]refidintro-prod-editinguh$h8h'Khhh(hh%F/chandler/home/floyd/public_html/gg/docs/source/intro/prod/editing.rstubh section)}(hhh](h title)}(h Data Editingh]h Data Editing}(hhNh(hh%Nh'Nubah}(h]h]h]h ]h"]uh$hLhhIh(hh%hFh'Kubh paragraph)}(hXYour solutions will be corrupted if you include data that are statistical outliers. For a time series from independent h-files, an error in one component of one site's coordinates will not affect any other site (and may have minimal effect on the other components), but as soon as you combine h-files with common parameters, whether multiple networks on a single day or the results from more than one day, the different estimates of the common parameter will clash, raising chi-square and distorting the solution. There are two ways you can remove outliers: using the :content:`rename` command of the :content:`eq_file` and using the :content:`sig_neu` command directly in the :program:`globk` command file. In the :content:`rename` command, if you designate the extent as :content:`_XCL` the site will always be excluded; with :content:`_XPS` the site will be excluded from any solution of more than one day but will still appear in daily time series; e.g.h](hX:Your solutions will be corrupted if you include data that are statistical outliers. For a time series from independent h-files, an error in one component of one site’s coordinates will not affect any other site (and may have minimal effect on the other components), but as soon as you combine h-files with common parameters, whether multiple networks on a single day or the results from more than one day, the different estimates of the common parameter will clash, raising chi-square and distorting the solution. There are two ways you can remove outliers: using the }(hh^h(hh%Nh'Nubh literal)}(h:content:`rename`h]hrename}(hhhh(hh%Nh'Nubah}(h]h](code highlightcontenttexteh]h ]h"]languagehvuh$hfhh^ubh command of the }(hh^h(hh%Nh'Nubhg)}(h:content:`eq_file`h]heq_file}(hhh(hh%Nh'Nubah}(h]h](hshthuhveh]h ]h"]languagehvuh$hfhh^ubh and using the }(hh^h(hh%Nh'Nubhg)}(h:content:`sig_neu`h]hsig_neu}(hhh(hh%Nh'Nubah}(h]h](hshthuhveh]h ]h"]languagehvuh$hfhh^ubh command directly in the }(hh^h(hh%Nh'Nubhliteral_strong)}(h:program:`globk`h]hglobk}(hhh(hh%Nh'Nubah}(h]h]programah]h ]h"]uh$hhh^ubh command file. In the }(hh^h(hh%Nh'Nubhg)}(h:content:`rename`h]hrename}(hhh(hh%Nh'Nubah}(h]h](hshthuhveh]h ]h"]languagehvuh$hfhh^ubh) command, if you designate the extent as }(hh^h(hh%Nh'Nubhg)}(h:content:`_XCL`h]h_XCL}(hhh(hh%Nh'Nubah}(h]h](hshthuhveh]h ]h"]languagehvuh$hfhh^ubh( the site will always be excluded; with }(hh^h(hh%Nh'Nubhg)}(h:content:`_XPS`h]h_XPS}(hhh(hh%Nh'Nubah}(h]h](hshthuhveh]h ]h"]languagehvuh$hfhh^ubhr the site will be excluded from any solution of more than one day but will still appear in daily time series; e.g.}(hh^h(hh%Nh'Nubeh}(h]h]h]h ]h"]uh$h\h%hFh'KhhIh(hubh literal_block)}(hCrename MIT algo_gps algo_xcl 1997 5 14 0 0 1997 5 17 24 0h]hCrename MIT algo_gps algo_xcl 1997 5 14 0 0 1997 5 17 24 0}hhsbah}(h]h]h]h ]h"] xml:spacepreserveforcelanguagetexthighlight_args}uh$hh%hFh'KhhIh(hubh])}(hXwill remove from the solution the observations of ALGO in any h-files containing the characters :content:`MIT` between 14 and 17 May 1997. Note that the date span specified must completely encompass the start time of the first h-file and the end time of the last h-file. The h-file selector entry and/or the dates may be omitted for a more expansive exclusion. With the :content:`sig_neu` command you can effectively remove the effect of an error in a particular component by adding sufficient noise that it has negligible weight, e.g.h](h`will remove from the solution the observations of ALGO in any h-files containing the characters }(hjh(hh%Nh'Nubhg)}(h:content:`MIT`h]hMIT}(hjh(hh%Nh'Nubah}(h]h](hshthuhveh]h ]h"]languagehvuh$hfhjubhX between 14 and 17 May 1997. Note that the date span specified must completely encompass the start time of the first h-file and the end time of the last h-file. The h-file selector entry and/or the dates may be omitted for a more expansive exclusion. With the }(hjh(hh%Nh'Nubhg)}(h:content:`sig_neu`h]hsig_neu}(hj+h(hh%Nh'Nubah}(h]h](hshthuhveh]h ]h"]languagehvuh$hfhjubh command you can effectively remove the effect of an error in a particular component by adding sufficient noise that it has negligible weight, e.g.}(hjh(hh%Nh'Nubeh}(h]h]h]h ]h"]uh$h\h%hFh'KhhIh(hubh)}(hDsig_neu MIT algo_gps 0 0 .5 1997 5 14 0 0 1997 5 18 0 0h]hDsig_neu MIT algo_gps 0 0 .5 1997 5 14 0 0 1997 5 18 0 0}hjDsbah}(h]h]h]h ]h"]j j j j textj}uh$hh%hFh'KhhIh(hubh])}(hwould reduce the effect of an error in height by adding :math:`0.25 \text{m}^2` in quadrature to the variance of the height component as estimated from the h-files specified.h](h8would reduce the effect of an error in height by adding }(hjTh(hh%Nh'Nubh math)}(h:math:`0.25 \text{m}^2`h]h0.25 \text{m}^2}(hj^h(hh%Nh'Nubah}(h]h]h]h ]h"]uh$j\hjTubh_ in quadrature to the variance of the height component as estimated from the h-files specified.}(hjTh(hh%Nh'Nubeh}(h]h]h]h ]h"]uh$h\h%hFh'KhhIh(hubh])}(hXWhich of these approaches you choose is a matter of taste and could depend on the tools you use for editing. For continuous data, when you can afford to remove many outliers without weakening the solution, using the :content:`nsigma` option of the :program:`tsfit` program (see Chapter 4 of the `GLOBK Reference Manual `_) to automatically generate :content:`_XPS` commands, or the MATLAB-based interactive program `tsview `_ to generate :content:`_XPS` commands with the click of a mouse, are quite convenient. For survey-mode date, when every observation counts and you often want to down-weight heights but not horizontal components or only one horizontal component, :content:`sig_neu` commands are more convenient. The GAMIT utility :program:`grw` can be used to generate :content:`sig_neu` commands with only a few key-strokes.h](hWhich of these approaches you choose is a matter of taste and could depend on the tools you use for editing. For continuous data, when you can afford to remove many outliers without weakening the solution, using the }(hjvh(hh%Nh'Nubhg)}(h:content:`nsigma`h]hnsigma}(hj~h(hh%Nh'Nubah}(h]h](hshthuhveh]h ]h"]languagehvuh$hfhjvubh option of the }(hjvh(hh%Nh'Nubh)}(h:program:`tsfit`h]htsfit}(hjh(hh%Nh'Nubah}(h]h]hah]h ]h"]uh$hhjvubh program (see Chapter 4 of the }(hjvh(hh%Nh'Nubh reference)}(hG`GLOBK Reference Manual `_h]hGLOBK Reference Manual}(hjh(hh%Nh'Nubah}(h]h]h]h ]h"]nameGLOBK Reference Manualrefuri+http://geoweb.mit.edu/gg/docs/GLOBK_Ref.pdfuh$jhjvubh9)}(h. h]h}(h]globk-reference-manualah]h]globk reference manualah ]h"]refurijuh$h8 referencedKhjvubh) to automatically generate }(hjvh(hh%Nh'Nubhg)}(h:content:`_XPS`h]h_XPS}(hjh(hh%Nh'Nubah}(h]h](hshthuhveh]h ]h"]languagehvuh$hfhjvubh3 commands, or the MATLAB-based interactive program }(hjvh(hh%Nh'Nubj)}(h0`tsview `_h]htsview}(hjh(hh%Nh'Nubah}(h]h]h]h ]h"]nametsviewj$http://geoweb.mit.edu/~tah/GGMatlab/uh$jhjvubh9)}(h' h]h}(h]tsviewah]h]tsviewah ]h"]refurijuh$h8jKhjvubh to generate }(hjvh(hh%Nh'Nubhg)}(h:content:`_XPS`h]h_XPS}(hjh(hh%Nh'Nubah}(h]h](hshthuhveh]h ]h"]languagehvuh$hfhjvubh commands with the click of a mouse, are quite convenient. For survey-mode date, when every observation counts and you often want to down-weight heights but not horizontal components or only one horizontal component, }(hjvh(hh%Nh'Nubhg)}(h:content:`sig_neu`h]hsig_neu}(hjh(hh%Nh'Nubah}(h]h](hshthuhveh]h ]h"]languagehvuh$hfhjvubh1 commands are more convenient. The GAMIT utility }(hjvh(hh%Nh'Nubh)}(h:program:`grw`h]hgrw}(hj$h(hh%Nh'Nubah}(h]h]hah]h ]h"]uh$hhjvubh can be used to generate }(hjvh(hh%Nh'Nubhg)}(h:content:`sig_neu`h]hsig_neu}(hj6h(hh%Nh'Nubah}(h]h](hshthuhveh]h ]h"]languagehvuh$hfhjvubh& commands with only a few key-strokes.}(hjvh(hh%Nh'Nubeh}(h]h]h]h ]h"]uh$h\h%hFh'KhhIh(hubeh}(h]( data-editinghEeh]h]( data editingintro_prod_editingeh ]h"]uh$hGhhh(hh%hFh'Kexpect_referenced_by_name}jUh:sexpect_referenced_by_id}hEh:subeh}(h]h]h]h ]h"]sourcehFtranslation_progress}(totalK translatedKuuh$hcurrent_sourceN current_lineNsettingsdocutils.frontendValues)}(outputNhLN generatorN datestampN source_linkN source_urlN toc_backlinksentryfootnote_backlinksK sectnum_xformKstrip_commentsNstrip_elements_with_classesN strip_classesN report_levelK halt_levelKexit_status_levelKdebugNwarning_streamN tracebackinput_encoding utf-8-siginput_encoding_error_handlerstrictoutput_encodingutf-8output_encoding_error_handlerjerror_encodingutf-8error_encoding_error_handlerbackslashreplace language_codeenrecord_dependenciesNconfigN id_prefixhauto_id_prefixid dump_settingsNdump_internalsNdump_transformsNdump_pseudo_xmlNexpose_internalsNstrict_visitorN_disable_configN_sourcehF _destinationN _config_files]file_insertion_enabled raw_enabledKline_length_limitM'pep_referencesN pep_base_urlhttps://peps.python.org/pep_file_url_templatepep-%04drfc_referencesN rfc_base_url&https://datatracker.ietf.org/doc/html/ tab_widthKtrim_footnote_reference_spacesyntax_highlightlong smart_quotessmartquotes_locales]character_level_inline_markupdoctitle_xform docinfo_xformKsectsubtitle_xform image_loadinglinkembed_stylesheetcloak_email_addressessection_self_linkenvNubreporterNindirect_targets]substitution_defs}(hh h5h)usubstitution_names}(relhdateh5urefnames}refids}hE]h:asnameids}(jUhEjTjQjjjju nametypes}(jUjTjjuh}(hEhIjQhIjjjju footnote_refs} citation_refs} autofootnotes]autofootnote_refs]symbol_footnotes]symbol_footnote_refs] footnotes] citations]autofootnote_startKsymbol_footnote_startK id_counter collectionsCounter}Rparse_messages]transform_messages]h system_message)}(hhh]h])}(hhh]h8Hyperlink target "intro-prod-editing" is not referenced.}hjsbah}(h]h]h]h ]h"]uh$h\hjubah}(h]h]h]h ]h"]levelKtypeINFOsourcehFlineKuh$juba transformerN include_log] decorationNh(hub.