diff --git a/notebooks/Socet2ISIS.ipynb b/notebooks/Socet2ISIS.ipynb
index c5130953f65be68c481936fb3d3731abc33f58ba..fe3a81a97930daff220df4be6cdf336f5081aaf4 100644
--- a/notebooks/Socet2ISIS.ipynb
+++ b/notebooks/Socet2ISIS.ipynb
@@ -2,39 +2,28 @@
"cells": [
{
"cell_type": "code",
- "execution_count": 67,
+ "execution_count": 12,
"metadata": {},
- "outputs": [
- {
- "data": {
- "text/plain": [
- "{'ATF_FILE': 'CTX_Athabasca_Middle_step0.atf',\n",
- " 'GP_FILE': 'CTX_Athabasca_Middle.gpf',\n",
- " 'IMAGE_IPF': ['P19_008344_1894_XN_09N203W.ipf',\n",
- " 'P20_008845_1894_XN_09N203W.ipf',\n",
- " 'P03_002371_1888_XI_08N204W.ipf',\n",
- " 'P01_001540_1889_XI_08N204W.ipf',\n",
- " 'P01_001606_1897_XI_09N203W.ipf',\n",
- " 'P03_002226_1895_XI_09N203W.ipf'],\n",
- " 'IMAGE_SUP': ['P19_008344_1894_XN_09N203W.sup',\n",
- " 'P20_008845_1894_XN_09N203W.sup',\n",
- " 'P03_002371_1888_XI_08N204W.sup',\n",
- " 'P01_001540_1889_XI_08N204W.sup',\n",
- " 'P01_001606_1897_XI_09N203W.sup',\n",
- " 'P03_002226_1895_XI_09N203W.sup'],\n",
- " 'PATH': '/home/tthatcher/Desktop/Projects/Plio/plio/notebooks',\n",
- " 'PROJECT': 'D:\\\\data\\\\CTX_Athabasca_Middle.prj'}"
- ]
- },
- "execution_count": 67,
- "metadata": {},
- "output_type": "execute_result"
- }
- ],
+ "outputs": [],
"source": [
"import os\n",
+ "import sys\n",
+ "from functools import singledispatch\n",
+ "import warnings\n",
+ "\n",
+ "import pandas as pd\n",
"import numpy as np\n",
"\n",
+ "from plio.examples import get_path\n",
+ "from plio.io.io_bae import read_gpf"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 2,
+ "metadata": {},
+ "outputs": [],
+ "source": [
"# Reads a .atf file and outputs all of the \n",
"# .ipf, .gpf, .sup, .prj, and path to locate the \n",
"# .apf file (should be the same as all others) \n",
@@ -82,64 +71,1952 @@
" \n",
" return files_dict\n",
" \n",
- "read_atf('CTX_Athabasca_Middle_step0.atf')"
+ "@singledispatch\n",
+ "def read_ipf(arg):\n",
+ " return str(arg)\n",
+ "\n",
+ "@read_ipf.register(str)\n",
+ "def read_ipf_str(input_data):\n",
+ " \"\"\"\n",
+ " Read a socet ipf file into a pandas data frame\n",
+ "\n",
+ " Parameters\n",
+ " ----------\n",
+ " input_data : str\n",
+ " path to the an input data file\n",
+ "\n",
+ " Returns\n",
+ " -------\n",
+ " df : pd.DataFrame\n",
+ " containing the ipf data with appropriate column names and indices\n",
+ " \"\"\"\n",
+ "\n",
+ " # Check that the number of rows is matching the expected number\n",
+ " with open(input_data, 'r') as f:\n",
+ " for i, l in enumerate(f):\n",
+ " if i == 1:\n",
+ " cnt = int(l)\n",
+ " elif i == 2:\n",
+ " col = l\n",
+ " break\n",
+ " \n",
+ " columns = np.genfromtxt(input_data, skip_header=2, dtype='unicode',\n",
+ " max_rows = 1, delimiter = ',')\n",
+ "\n",
+ " # TODO: Add unicode conversion\n",
+ " d = [line.split() for line in open(input_data, 'r')]\n",
+ " d = np.hstack(np.array(d[3:]))\n",
+ " \n",
+ " d = d.reshape(-1, 12)\n",
+ " \n",
+ " df = pd.DataFrame(d, columns=columns)\n",
+ " df['ipf_file'] = pd.Series(np.full((len(df['pt_id'])), input_data), index = df.index)\n",
+ "\n",
+ " assert int(cnt) == len(df), 'Dataframe length {} does not match point length {}.'.format(int(cnt), len(df))\n",
+ " \n",
+ " # Soft conversion of numeric types to numerics, allows str in first col for point_id\n",
+ " df = df.apply(pd.to_numeric, errors='ignore')\n",
+ "\n",
+ " return df\n",
+ "\n",
+ "@read_ipf.register(list)\n",
+ "def read_ipf_list(input_data_list):\n",
+ " \"\"\"\n",
+ " Read a socet ipf file into a pandas data frame\n",
+ "\n",
+ " Parameters\n",
+ " ----------\n",
+ " input_data_list : list\n",
+ " list of paths to the a set of input data files\n",
+ "\n",
+ " Returns\n",
+ " -------\n",
+ " df : pd.DataFrame\n",
+ " containing the ipf data with appropriate column names and indices\n",
+ " \"\"\"\n",
+ " frames = []\n",
+ "\n",
+ " for input_file in input_data_list:\n",
+ " frames.append(read_ipf(input_file))\n",
+ "\n",
+ " df = pd.concat(frames)\n",
+ "\n",
+ " return df"
]
},
{
"cell_type": "code",
- "execution_count": null,
- "metadata": {},
- "outputs": [],
- "source": []
- },
- {
- "cell_type": "code",
- "execution_count": null,
+ "execution_count": 22,
"metadata": {},
- "outputs": [],
- "source": []
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "metadata": {},
- "outputs": [],
- "source": []
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "metadata": {},
- "outputs": [],
- "source": []
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "metadata": {},
- "outputs": [],
- "source": []
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "metadata": {},
- "outputs": [],
- "source": []
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "metadata": {},
- "outputs": [],
- "source": []
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "metadata": {},
- "outputs": [],
- "source": []
+ "outputs": [
+ {
+ "name": "stderr",
+ "output_type": "stream",
+ "text": [
+ "/Users/adampaquette/anaconda/envs/pysat/lib/python3.6/site-packages/ipykernel_launcher.py:12: UserWarning: The following points found in ipf files missing from gpf file: \n",
+ "\n",
+ "['P03_002226_1895_XI_09N203W_15', 'P03_002226_1895_XI_09N203W_16', 'P03_002226_1895_XI_09N203W_17', 'P03_002226_1895_XI_09N203W_18', 'P03_002226_1895_XI_09N203W_19', 'P03_002226_1895_XI_09N203W_20', 'P03_002226_1895_XI_09N203W_21', 'P03_002226_1895_XI_09N203W_22', 'P03_002226_1895_XI_09N203W_24', 'P03_002226_1895_XI_09N203W_26', 'P03_002226_1895_XI_09N203W_30', 'P03_002226_1895_XI_09N203W_31', 'P03_002226_1895_XI_09N203W_32', 'P03_002226_1895_XI_09N203W_34', 'P03_002226_1895_XI_09N203W_36', 'P03_002226_1895_XI_09N203W_37', 'P03_002226_1895_XI_09N203W_44', 'P03_002226_1895_XI_09N203W_48', 'P03_002226_1895_XI_09N203W_49', 'P03_002226_1895_XI_09N203W_56', 'P03_002226_1895_XI_09N203W_57', 'P03_002226_1895_XI_09N203W_61', 'P03_002226_1895_XI_09N203W_62', 'P03_002226_1895_XI_09N203W_63', 'P03_002226_1895_XI_09N203W_65', 'P19_008344_1894_XN_09N203W_4', 'P20_008845_1894_XN_09N203W_15']. \n",
+ "\n",
+ "Continuing, but these points will be missing from the control network\n",
+ " if sys.path[0] == '':\n"
+ ]
+ },
+ {
+ "data": {
+ "text/html": [
+ "<div>\n",
+ "<style scoped>\n",
+ " .dataframe tbody tr th:only-of-type {\n",
+ " vertical-align: middle;\n",
+ " }\n",
+ "\n",
+ " .dataframe tbody tr th {\n",
+ " vertical-align: top;\n",
+ " }\n",
+ "\n",
+ " .dataframe thead th {\n",
+ " text-align: right;\n",
+ " }\n",
+ "</style>\n",
+ "<table border=\"1\" class=\"dataframe\">\n",
+ " <thead>\n",
+ " <tr style=\"text-align: right;\">\n",
+ " <th></th>\n",
+ " <th>val</th>\n",
+ " <th>fid_val</th>\n",
+ " <th>no_obs</th>\n",
+ " <th>l.</th>\n",
+ " <th>s.</th>\n",
+ " <th>sig_l</th>\n",
+ " <th>sig_s</th>\n",
+ " <th>res_l</th>\n",
+ " <th>res_s</th>\n",
+ " <th>fid_x</th>\n",
+ " <th>...</th>\n",
+ " <th>known</th>\n",
+ " <th>lat_Y_North</th>\n",
+ " <th>long_X_East</th>\n",
+ " <th>ht</th>\n",
+ " <th>sig0</th>\n",
+ " <th>sig1</th>\n",
+ " <th>sig2</th>\n",
+ " <th>res0</th>\n",
+ " <th>res1</th>\n",
+ " <th>res2</th>\n",
+ " </tr>\n",
+ " </thead>\n",
+ " <tbody>\n",
+ " <tr>\n",
+ " <th>10_8344_8845_4r</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>-2919.380615</td>\n",
+ " <td>1052.729004</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.059914</td>\n",
+ " <td>0.930311</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.159378</td>\n",
+ " <td>2.724649</td>\n",
+ " <td>-2523.828227</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>25.000000</td>\n",
+ " <td>18.301328</td>\n",
+ " <td>44.206259</td>\n",
+ " <td>416.201741</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>10_8344_8845_4r</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>-4044.697510</td>\n",
+ " <td>1008.950928</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.063678</td>\n",
+ " <td>0.661294</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.159378</td>\n",
+ " <td>2.724649</td>\n",
+ " <td>-2523.828227</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>25.000000</td>\n",
+ " <td>18.301328</td>\n",
+ " <td>44.206259</td>\n",
+ " <td>416.201741</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>10_8344_8845_4r</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>1700.584473</td>\n",
+ " <td>-2390.001709</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-0.063695</td>\n",
+ " <td>-1.055619</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.159378</td>\n",
+ " <td>2.724649</td>\n",
+ " <td>-2523.828227</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>25.000000</td>\n",
+ " <td>18.301328</td>\n",
+ " <td>44.206259</td>\n",
+ " <td>416.201741</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>10_8344_8845_4r</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>2006.141113</td>\n",
+ " <td>-2234.915283</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.318779</td>\n",
+ " <td>-0.585138</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.159378</td>\n",
+ " <td>2.724649</td>\n",
+ " <td>-2523.828227</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>25.000000</td>\n",
+ " <td>18.301328</td>\n",
+ " <td>44.206259</td>\n",
+ " <td>416.201741</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>11_8344_8845_4r</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>370.151917</td>\n",
+ " <td>2340.653076</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-0.239226</td>\n",
+ " <td>-1.113693</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.164905</td>\n",
+ " <td>2.721815</td>\n",
+ " <td>-2445.237027</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>30.000000</td>\n",
+ " <td>-22.046575</td>\n",
+ " <td>103.403228</td>\n",
+ " <td>173.210013</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>11_8344_8845_4r</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>-761.216064</td>\n",
+ " <td>2303.787109</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-0.247307</td>\n",
+ " <td>-0.587299</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.164905</td>\n",
+ " <td>2.721815</td>\n",
+ " <td>-2445.237027</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>30.000000</td>\n",
+ " <td>-22.046575</td>\n",
+ " <td>103.403228</td>\n",
+ " <td>173.210013</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>11_8344_8845_4r</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>4985.437988</td>\n",
+ " <td>-1070.364990</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.239207</td>\n",
+ " <td>1.080801</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.164905</td>\n",
+ " <td>2.721815</td>\n",
+ " <td>-2445.237027</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>30.000000</td>\n",
+ " <td>-22.046575</td>\n",
+ " <td>103.403228</td>\n",
+ " <td>173.210013</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>11_8344_8845_4r</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>5293.700195</td>\n",
+ " <td>-993.390625</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.672877</td>\n",
+ " <td>0.626316</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.164905</td>\n",
+ " <td>2.721815</td>\n",
+ " <td>-2445.237027</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>30.000000</td>\n",
+ " <td>-22.046575</td>\n",
+ " <td>103.403228</td>\n",
+ " <td>173.210013</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>12_8344_8845_4r</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>236.395218</td>\n",
+ " <td>984.833252</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.073683</td>\n",
+ " <td>-0.237114</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.164949</td>\n",
+ " <td>2.724076</td>\n",
+ " <td>-2606.935163</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>100.000000</td>\n",
+ " <td>-7.549561</td>\n",
+ " <td>93.170584</td>\n",
+ " <td>180.058858</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>12_8344_8845_4r</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>-889.364441</td>\n",
+ " <td>966.533997</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-0.251753</td>\n",
+ " <td>0.008170</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.164949</td>\n",
+ " <td>2.724076</td>\n",
+ " <td>-2606.935163</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>100.000000</td>\n",
+ " <td>-7.549561</td>\n",
+ " <td>93.170584</td>\n",
+ " <td>180.058858</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>12_8344_8845_4r</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>4856.525391</td>\n",
+ " <td>-2439.154785</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.133406</td>\n",
+ " <td>0.068024</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.164949</td>\n",
+ " <td>2.724076</td>\n",
+ " <td>-2606.935163</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>100.000000</td>\n",
+ " <td>-7.549561</td>\n",
+ " <td>93.170584</td>\n",
+ " <td>180.058858</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>12_8344_8845_4r</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>5158.205566</td>\n",
+ " <td>-2295.737549</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.487209</td>\n",
+ " <td>0.036399</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.164949</td>\n",
+ " <td>2.724076</td>\n",
+ " <td>-2606.935163</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>100.000000</td>\n",
+ " <td>-7.549561</td>\n",
+ " <td>93.170584</td>\n",
+ " <td>180.058858</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>13_8344_8845_4r</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>-1431.023804</td>\n",
+ " <td>1817.888672</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-0.056901</td>\n",
+ " <td>-0.383259</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.161840</td>\n",
+ " <td>2.723059</td>\n",
+ " <td>-2551.901554</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>4.536068</td>\n",
+ " <td>-1.643694</td>\n",
+ " <td>74.244153</td>\n",
+ " <td>319.554535</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>13_8344_8845_4r</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>-2559.871338</td>\n",
+ " <td>1777.522827</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-0.131265</td>\n",
+ " <td>-0.305102</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.161840</td>\n",
+ " <td>2.723059</td>\n",
+ " <td>-2551.901554</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>4.536068</td>\n",
+ " <td>-1.643694</td>\n",
+ " <td>74.244153</td>\n",
+ " <td>319.554535</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>13_8344_8845_4r</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>3186.063232</td>\n",
+ " <td>-1598.743530</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.019201</td>\n",
+ " <td>0.485372</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.161840</td>\n",
+ " <td>2.723059</td>\n",
+ " <td>-2551.901554</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>4.536068</td>\n",
+ " <td>-1.643694</td>\n",
+ " <td>74.244153</td>\n",
+ " <td>319.554535</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>13_8344_8845_4r</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>3493.876221</td>\n",
+ " <td>-1493.029175</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.568393</td>\n",
+ " <td>0.146750</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.161840</td>\n",
+ " <td>2.723059</td>\n",
+ " <td>-2551.901554</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>4.536068</td>\n",
+ " <td>-1.643694</td>\n",
+ " <td>74.244153</td>\n",
+ " <td>319.554535</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>14_8344_8845_4r</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>3518.268066</td>\n",
+ " <td>2495.108643</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.268449</td>\n",
+ " <td>1.000455</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.170415</td>\n",
+ " <td>2.720880</td>\n",
+ " <td>-2505.953426</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>5.000000</td>\n",
+ " <td>-50.069808</td>\n",
+ " <td>144.316524</td>\n",
+ " <td>-164.540707</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>14_8344_8845_4r</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>2385.278320</td>\n",
+ " <td>2476.032227</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.389532</td>\n",
+ " <td>-0.299258</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.170415</td>\n",
+ " <td>2.720880</td>\n",
+ " <td>-2505.953426</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>5.000000</td>\n",
+ " <td>-50.069808</td>\n",
+ " <td>144.316524</td>\n",
+ " <td>-164.540707</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>14_8344_8845_4r</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>8131.361816</td>\n",
+ " <td>-902.645325</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-0.469784</td>\n",
+ " <td>0.023018</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.170415</td>\n",
+ " <td>2.720880</td>\n",
+ " <td>-2505.953426</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>5.000000</td>\n",
+ " <td>-50.069808</td>\n",
+ " <td>144.316524</td>\n",
+ " <td>-164.540707</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>14_8344_8845_4r</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>8437.115234</td>\n",
+ " <td>-842.309326</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.289877</td>\n",
+ " <td>-0.829150</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.170415</td>\n",
+ " <td>2.720880</td>\n",
+ " <td>-2505.953426</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>5.000000</td>\n",
+ " <td>-50.069808</td>\n",
+ " <td>144.316524</td>\n",
+ " <td>-164.540707</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>15_8344_8845_4r_mt_z</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>3522.998535</td>\n",
+ " <td>1049.406982</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.466921</td>\n",
+ " <td>1.168062</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>1</td>\n",
+ " <td>0.170723</td>\n",
+ " <td>2.723266</td>\n",
+ " <td>-2502.470000</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>5.000000</td>\n",
+ " <td>-34.817656</td>\n",
+ " <td>145.864550</td>\n",
+ " <td>0.214023</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>15_8344_8845_4r_mt_z</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>2395.869385</td>\n",
+ " <td>1038.165405</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-0.170850</td>\n",
+ " <td>0.372485</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>1</td>\n",
+ " <td>0.170723</td>\n",
+ " <td>2.723266</td>\n",
+ " <td>-2502.470000</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>5.000000</td>\n",
+ " <td>-34.817656</td>\n",
+ " <td>145.864550</td>\n",
+ " <td>0.214023</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>15_8344_8845_4r_mt_z</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>8142.054688</td>\n",
+ " <td>-2364.445557</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-0.072035</td>\n",
+ " <td>-0.832436</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>1</td>\n",
+ " <td>0.170723</td>\n",
+ " <td>2.723266</td>\n",
+ " <td>-2502.470000</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>5.000000</td>\n",
+ " <td>-34.817656</td>\n",
+ " <td>145.864550</td>\n",
+ " <td>0.214023</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>15_8344_8845_4r_mt_z</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>8440.085938</td>\n",
+ " <td>-2219.049805</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.274030</td>\n",
+ " <td>-0.850105</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>1</td>\n",
+ " <td>0.170723</td>\n",
+ " <td>2.723266</td>\n",
+ " <td>-2502.470000</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>5.000000</td>\n",
+ " <td>-34.817656</td>\n",
+ " <td>145.864550</td>\n",
+ " <td>0.214023</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>16_8344_8845_4r</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>1885.346313</td>\n",
+ " <td>1807.367065</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-0.028473</td>\n",
+ " <td>-0.238015</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.167682</td>\n",
+ " <td>2.722359</td>\n",
+ " <td>-2558.312931</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>1.707214</td>\n",
+ " <td>-29.465246</td>\n",
+ " <td>121.908506</td>\n",
+ " <td>19.874949</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>16_8344_8845_4r</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>756.099792</td>\n",
+ " <td>1785.494751</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.165338</td>\n",
+ " <td>-0.327408</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.167682</td>\n",
+ " <td>2.722359</td>\n",
+ " <td>-2558.312931</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>1.707214</td>\n",
+ " <td>-29.465246</td>\n",
+ " <td>121.908506</td>\n",
+ " <td>19.874949</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>16_8344_8845_4r</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>6502.070313</td>\n",
+ " <td>-1593.964233</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-0.007916</td>\n",
+ " <td>0.463232</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.167682</td>\n",
+ " <td>2.722359</td>\n",
+ " <td>-2558.312931</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>1.707214</td>\n",
+ " <td>-29.465246</td>\n",
+ " <td>121.908506</td>\n",
+ " <td>19.874949</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>16_8344_8845_4r</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>6805.896973</td>\n",
+ " <td>-1493.774048</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.328896</td>\n",
+ " <td>0.019526</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.167682</td>\n",
+ " <td>2.722359</td>\n",
+ " <td>-2558.312931</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>1.707214</td>\n",
+ " <td>-29.465246</td>\n",
+ " <td>121.908506</td>\n",
+ " <td>19.874949</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>17_8344_8845_2r_mt_z</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>-2015.364380</td>\n",
+ " <td>465.268188</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.010045</td>\n",
+ " <td>0.244017</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>1</td>\n",
+ " <td>0.161092</td>\n",
+ " <td>2.725426</td>\n",
+ " <td>-2590.130000</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>2.000000</td>\n",
+ " <td>-10.779830</td>\n",
+ " <td>38.761214</td>\n",
+ " <td>0.236077</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>17_8344_8845_2r_mt_z</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>-3138.269531</td>\n",
+ " <td>442.515503</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-0.023370</td>\n",
+ " <td>-0.248643</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>1</td>\n",
+ " <td>0.161092</td>\n",
+ " <td>2.725426</td>\n",
+ " <td>-2590.130000</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>2.000000</td>\n",
+ " <td>-10.779830</td>\n",
+ " <td>38.761214</td>\n",
+ " <td>0.236077</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>...</th>\n",
+ " <td>...</td>\n",
+ " <td>...</td>\n",
+ " <td>...</td>\n",
+ " <td>...</td>\n",
+ " <td>...</td>\n",
+ " <td>...</td>\n",
+ " <td>...</td>\n",
+ " <td>...</td>\n",
+ " <td>...</td>\n",
+ " <td>...</td>\n",
+ " <td>...</td>\n",
+ " <td>...</td>\n",
+ " <td>...</td>\n",
+ " <td>...</td>\n",
+ " <td>...</td>\n",
+ " <td>...</td>\n",
+ " <td>...</td>\n",
+ " <td>...</td>\n",
+ " <td>...</td>\n",
+ " <td>...</td>\n",
+ " <td>...</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_18</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>1</td>\n",
+ " <td>8918.806641</td>\n",
+ " <td>-1701.755249</td>\n",
+ " <td>0.022017</td>\n",
+ " <td>0.022017</td>\n",
+ " <td>0.734173</td>\n",
+ " <td>-0.074947</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.171455</td>\n",
+ " <td>2.722281</td>\n",
+ " <td>-2514.431453</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-38.388316</td>\n",
+ " <td>195.050373</td>\n",
+ " <td>-391.289238</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_2</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>1</td>\n",
+ " <td>-2277.691406</td>\n",
+ " <td>-2240.823486</td>\n",
+ " <td>0.144273</td>\n",
+ " <td>0.144273</td>\n",
+ " <td>-0.061422</td>\n",
+ " <td>-0.263554</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.161184</td>\n",
+ " <td>2.729925</td>\n",
+ " <td>-2535.712262</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>42.067289</td>\n",
+ " <td>-12.658160</td>\n",
+ " <td>627.742625</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_2</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>-3389.000000</td>\n",
+ " <td>-2155.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-0.073150</td>\n",
+ " <td>-0.271173</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.161184</td>\n",
+ " <td>2.729925</td>\n",
+ " <td>-2535.712262</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>42.067289</td>\n",
+ " <td>-12.658160</td>\n",
+ " <td>627.742625</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_2</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>1</td>\n",
+ " <td>-5474.326660</td>\n",
+ " <td>842.207031</td>\n",
+ " <td>0.275627</td>\n",
+ " <td>0.275627</td>\n",
+ " <td>0.164525</td>\n",
+ " <td>0.349547</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.161184</td>\n",
+ " <td>2.729925</td>\n",
+ " <td>-2535.712262</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>42.067289</td>\n",
+ " <td>-12.658160</td>\n",
+ " <td>627.742625</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_2</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>-4613.785645</td>\n",
+ " <td>643.217224</td>\n",
+ " <td>0.091086</td>\n",
+ " <td>0.091086</td>\n",
+ " <td>0.132736</td>\n",
+ " <td>0.239899</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.161184</td>\n",
+ " <td>2.729925</td>\n",
+ " <td>-2535.712262</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>42.067289</td>\n",
+ " <td>-12.658160</td>\n",
+ " <td>627.742625</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_3</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>1</td>\n",
+ " <td>-1801.096436</td>\n",
+ " <td>-2078.479004</td>\n",
+ " <td>0.062058</td>\n",
+ " <td>0.062058</td>\n",
+ " <td>-0.198696</td>\n",
+ " <td>-0.105860</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.161992</td>\n",
+ " <td>2.729563</td>\n",
+ " <td>-2525.215515</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>36.033823</td>\n",
+ " <td>-4.170818</td>\n",
+ " <td>616.480181</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_3</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>-2913.000000</td>\n",
+ " <td>-2000.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-0.088180</td>\n",
+ " <td>-0.049593</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.161992</td>\n",
+ " <td>2.729563</td>\n",
+ " <td>-2525.215515</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>36.033823</td>\n",
+ " <td>-4.170818</td>\n",
+ " <td>616.480181</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_3</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>1</td>\n",
+ " <td>-4997.791504</td>\n",
+ " <td>1000.375854</td>\n",
+ " <td>0.132941</td>\n",
+ " <td>0.132941</td>\n",
+ " <td>0.257260</td>\n",
+ " <td>0.117777</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.161992</td>\n",
+ " <td>2.729563</td>\n",
+ " <td>-2525.215515</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>36.033823</td>\n",
+ " <td>-4.170818</td>\n",
+ " <td>616.480181</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_3</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>-4137.536133</td>\n",
+ " <td>774.480347</td>\n",
+ " <td>0.185544</td>\n",
+ " <td>0.185544</td>\n",
+ " <td>0.199301</td>\n",
+ " <td>0.065777</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.161992</td>\n",
+ " <td>2.729563</td>\n",
+ " <td>-2525.215515</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>36.033823</td>\n",
+ " <td>-4.170818</td>\n",
+ " <td>616.480181</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_4</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>1</td>\n",
+ " <td>-2103.971680</td>\n",
+ " <td>184.334869</td>\n",
+ " <td>0.146920</td>\n",
+ " <td>0.146920</td>\n",
+ " <td>0.182557</td>\n",
+ " <td>-0.000683</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.160994</td>\n",
+ " <td>2.725912</td>\n",
+ " <td>-2562.446851</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-7.233876</td>\n",
+ " <td>33.324946</td>\n",
+ " <td>574.887719</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_4</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>1</td>\n",
+ " <td>-3226.000000</td>\n",
+ " <td>167.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-0.195774</td>\n",
+ " <td>-0.000066</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.160994</td>\n",
+ " <td>2.725912</td>\n",
+ " <td>-2562.446851</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-7.233876</td>\n",
+ " <td>33.324946</td>\n",
+ " <td>574.887719</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_5</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>1</td>\n",
+ " <td>-1813.285278</td>\n",
+ " <td>3.928788</td>\n",
+ " <td>0.019047</td>\n",
+ " <td>0.019047</td>\n",
+ " <td>-0.235171</td>\n",
+ " <td>-0.174079</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.161544</td>\n",
+ " <td>2.726149</td>\n",
+ " <td>-2560.812028</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>7.498790</td>\n",
+ " <td>33.109007</td>\n",
+ " <td>581.149470</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_5</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>1</td>\n",
+ " <td>-2934.000000</td>\n",
+ " <td>-6.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-0.049005</td>\n",
+ " <td>0.437424</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.161544</td>\n",
+ " <td>2.726149</td>\n",
+ " <td>-2560.812028</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>7.498790</td>\n",
+ " <td>33.109007</td>\n",
+ " <td>581.149470</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_5</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>-4165.228516</td>\n",
+ " <td>2507.076660</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.447100</td>\n",
+ " <td>-0.219764</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.161544</td>\n",
+ " <td>2.726149</td>\n",
+ " <td>-2560.812028</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>7.498790</td>\n",
+ " <td>33.109007</td>\n",
+ " <td>581.149470</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_6</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>1</td>\n",
+ " <td>-2250.549072</td>\n",
+ " <td>1572.788086</td>\n",
+ " <td>0.025015</td>\n",
+ " <td>0.025015</td>\n",
+ " <td>-0.146990</td>\n",
+ " <td>0.192274</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.160447</td>\n",
+ " <td>2.723639</td>\n",
+ " <td>-2605.266130</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>7.687475</td>\n",
+ " <td>59.859663</td>\n",
+ " <td>375.680663</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_6</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>3</td>\n",
+ " <td>-3378.000000</td>\n",
+ " <td>1532.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.026541</td>\n",
+ " <td>0.063346</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.160447</td>\n",
+ " <td>2.723639</td>\n",
+ " <td>-2605.266130</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>7.687475</td>\n",
+ " <td>59.859663</td>\n",
+ " <td>375.680663</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_6</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>1</td>\n",
+ " <td>2367.587646</td>\n",
+ " <td>-1849.272095</td>\n",
+ " <td>0.025110</td>\n",
+ " <td>0.025110</td>\n",
+ " <td>-0.040890</td>\n",
+ " <td>-0.136247</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.160447</td>\n",
+ " <td>2.723639</td>\n",
+ " <td>-2605.266130</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>7.687475</td>\n",
+ " <td>59.859663</td>\n",
+ " <td>375.680663</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_6</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>1</td>\n",
+ " <td>2675.137207</td>\n",
+ " <td>-1730.558105</td>\n",
+ " <td>0.064276</td>\n",
+ " <td>0.064276</td>\n",
+ " <td>0.546035</td>\n",
+ " <td>-0.185103</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.160447</td>\n",
+ " <td>2.723639</td>\n",
+ " <td>-2605.266130</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>7.687475</td>\n",
+ " <td>59.859663</td>\n",
+ " <td>375.680663</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_7</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>1</td>\n",
+ " <td>-2249.931396</td>\n",
+ " <td>1712.971436</td>\n",
+ " <td>0.010596</td>\n",
+ " <td>0.010596</td>\n",
+ " <td>-0.102909</td>\n",
+ " <td>0.079389</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.160419</td>\n",
+ " <td>2.723407</td>\n",
+ " <td>-2619.498291</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>6.358447</td>\n",
+ " <td>61.278242</td>\n",
+ " <td>375.840988</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_7</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>3</td>\n",
+ " <td>-3378.000000</td>\n",
+ " <td>1672.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.055623</td>\n",
+ " <td>0.013456</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.160419</td>\n",
+ " <td>2.723407</td>\n",
+ " <td>-2619.498291</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>6.358447</td>\n",
+ " <td>61.278242</td>\n",
+ " <td>375.840988</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_7</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>1</td>\n",
+ " <td>2367.507813</td>\n",
+ " <td>-1706.241821</td>\n",
+ " <td>0.047435</td>\n",
+ " <td>0.047435</td>\n",
+ " <td>-0.157931</td>\n",
+ " <td>-0.042074</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.160419</td>\n",
+ " <td>2.723407</td>\n",
+ " <td>-2619.498291</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>6.358447</td>\n",
+ " <td>61.278242</td>\n",
+ " <td>375.840988</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_7</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>1</td>\n",
+ " <td>2675.892578</td>\n",
+ " <td>-1596.527100</td>\n",
+ " <td>0.058218</td>\n",
+ " <td>0.058218</td>\n",
+ " <td>0.589561</td>\n",
+ " <td>-0.111421</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.160419</td>\n",
+ " <td>2.723407</td>\n",
+ " <td>-2619.498291</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>6.358447</td>\n",
+ " <td>61.278242</td>\n",
+ " <td>375.840988</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_8</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>1</td>\n",
+ " <td>-1664.131104</td>\n",
+ " <td>1875.622437</td>\n",
+ " <td>0.017266</td>\n",
+ " <td>0.017266</td>\n",
+ " <td>-0.239689</td>\n",
+ " <td>-0.402891</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.161417</td>\n",
+ " <td>2.723012</td>\n",
+ " <td>-2608.028730</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-0.219662</td>\n",
+ " <td>71.414267</td>\n",
+ " <td>336.409684</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_8</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>3</td>\n",
+ " <td>-2793.000000</td>\n",
+ " <td>1837.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-0.117483</td>\n",
+ " <td>-0.229337</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.161417</td>\n",
+ " <td>2.723012</td>\n",
+ " <td>-2608.028730</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-0.219662</td>\n",
+ " <td>71.414267</td>\n",
+ " <td>336.409684</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_8</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>1</td>\n",
+ " <td>2953.094727</td>\n",
+ " <td>-1538.956543</td>\n",
+ " <td>0.023364</td>\n",
+ " <td>0.023364</td>\n",
+ " <td>0.177406</td>\n",
+ " <td>0.403161</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.161417</td>\n",
+ " <td>2.723012</td>\n",
+ " <td>-2608.028730</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-0.219662</td>\n",
+ " <td>71.414267</td>\n",
+ " <td>336.409684</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_8</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>1</td>\n",
+ " <td>3261.290039</td>\n",
+ " <td>-1439.873169</td>\n",
+ " <td>0.023538</td>\n",
+ " <td>0.023538</td>\n",
+ " <td>0.575510</td>\n",
+ " <td>0.180705</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.161417</td>\n",
+ " <td>2.723012</td>\n",
+ " <td>-2608.028730</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-0.219662</td>\n",
+ " <td>71.414267</td>\n",
+ " <td>336.409684</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_9</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>1</td>\n",
+ " <td>332.072113</td>\n",
+ " <td>-2045.152222</td>\n",
+ " <td>0.232865</td>\n",
+ " <td>0.232865</td>\n",
+ " <td>0.095817</td>\n",
+ " <td>-0.483899</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.165741</td>\n",
+ " <td>2.729053</td>\n",
+ " <td>-2586.686862</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>15.191825</td>\n",
+ " <td>22.094037</td>\n",
+ " <td>544.874936</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_9</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>-780.000000</td>\n",
+ " <td>-1953.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.298743</td>\n",
+ " <td>-0.045963</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.165741</td>\n",
+ " <td>2.729053</td>\n",
+ " <td>-2586.686862</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>15.191825</td>\n",
+ " <td>22.094037</td>\n",
+ " <td>544.874936</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_9</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>-2866.689453</td>\n",
+ " <td>1031.006104</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>-0.038437</td>\n",
+ " <td>0.471309</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.165741</td>\n",
+ " <td>2.729053</td>\n",
+ " <td>-2586.686862</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>15.191825</td>\n",
+ " <td>22.094037</td>\n",
+ " <td>544.874936</td>\n",
+ " </tr>\n",
+ " <tr>\n",
+ " <th>P20_008845_1894_XN_09N203W_9</th>\n",
+ " <td>1</td>\n",
+ " <td>0</td>\n",
+ " <td>0</td>\n",
+ " <td>-2002.366211</td>\n",
+ " <td>824.604126</td>\n",
+ " <td>0.397616</td>\n",
+ " <td>0.397616</td>\n",
+ " <td>-0.141172</td>\n",
+ " <td>0.105223</td>\n",
+ " <td>0.0</td>\n",
+ " <td>...</td>\n",
+ " <td>0</td>\n",
+ " <td>0.165741</td>\n",
+ " <td>2.729053</td>\n",
+ " <td>-2586.686862</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.0</td>\n",
+ " <td>0.000000</td>\n",
+ " <td>15.191825</td>\n",
+ " <td>22.094037</td>\n",
+ " <td>544.874936</td>\n",
+ " </tr>\n",
+ " </tbody>\n",
+ "</table>\n",
+ "<p>919 rows × 23 columns</p>\n",
+ "</div>"
+ ],
+ "text/plain": [
+ " val fid_val no_obs l. s. \\\n",
+ "10_8344_8845_4r 1 0 0 -2919.380615 1052.729004 \n",
+ "10_8344_8845_4r 1 0 0 -4044.697510 1008.950928 \n",
+ "10_8344_8845_4r 1 0 0 1700.584473 -2390.001709 \n",
+ "10_8344_8845_4r 1 0 0 2006.141113 -2234.915283 \n",
+ "11_8344_8845_4r 1 0 0 370.151917 2340.653076 \n",
+ "11_8344_8845_4r 1 0 0 -761.216064 2303.787109 \n",
+ "11_8344_8845_4r 1 0 0 4985.437988 -1070.364990 \n",
+ "11_8344_8845_4r 1 0 0 5293.700195 -993.390625 \n",
+ "12_8344_8845_4r 1 0 0 236.395218 984.833252 \n",
+ "12_8344_8845_4r 1 0 0 -889.364441 966.533997 \n",
+ "12_8344_8845_4r 1 0 0 4856.525391 -2439.154785 \n",
+ "12_8344_8845_4r 1 0 0 5158.205566 -2295.737549 \n",
+ "13_8344_8845_4r 1 0 0 -1431.023804 1817.888672 \n",
+ "13_8344_8845_4r 1 0 0 -2559.871338 1777.522827 \n",
+ "13_8344_8845_4r 1 0 0 3186.063232 -1598.743530 \n",
+ "13_8344_8845_4r 1 0 0 3493.876221 -1493.029175 \n",
+ "14_8344_8845_4r 1 0 0 3518.268066 2495.108643 \n",
+ "14_8344_8845_4r 1 0 0 2385.278320 2476.032227 \n",
+ "14_8344_8845_4r 1 0 0 8131.361816 -902.645325 \n",
+ "14_8344_8845_4r 1 0 0 8437.115234 -842.309326 \n",
+ "15_8344_8845_4r_mt_z 1 0 0 3522.998535 1049.406982 \n",
+ "15_8344_8845_4r_mt_z 1 0 0 2395.869385 1038.165405 \n",
+ "15_8344_8845_4r_mt_z 1 0 0 8142.054688 -2364.445557 \n",
+ "15_8344_8845_4r_mt_z 1 0 0 8440.085938 -2219.049805 \n",
+ "16_8344_8845_4r 1 0 0 1885.346313 1807.367065 \n",
+ "16_8344_8845_4r 1 0 0 756.099792 1785.494751 \n",
+ "16_8344_8845_4r 1 0 0 6502.070313 -1593.964233 \n",
+ "16_8344_8845_4r 1 0 0 6805.896973 -1493.774048 \n",
+ "17_8344_8845_2r_mt_z 1 0 0 -2015.364380 465.268188 \n",
+ "17_8344_8845_2r_mt_z 1 0 0 -3138.269531 442.515503 \n",
+ "... ... ... ... ... ... \n",
+ "P20_008845_1894_XN_09N203W_18 1 0 1 8918.806641 -1701.755249 \n",
+ "P20_008845_1894_XN_09N203W_2 1 0 1 -2277.691406 -2240.823486 \n",
+ "P20_008845_1894_XN_09N203W_2 1 0 0 -3389.000000 -2155.000000 \n",
+ "P20_008845_1894_XN_09N203W_2 1 0 1 -5474.326660 842.207031 \n",
+ "P20_008845_1894_XN_09N203W_2 1 0 0 -4613.785645 643.217224 \n",
+ "P20_008845_1894_XN_09N203W_3 1 0 1 -1801.096436 -2078.479004 \n",
+ "P20_008845_1894_XN_09N203W_3 1 0 0 -2913.000000 -2000.000000 \n",
+ "P20_008845_1894_XN_09N203W_3 1 0 1 -4997.791504 1000.375854 \n",
+ "P20_008845_1894_XN_09N203W_3 1 0 0 -4137.536133 774.480347 \n",
+ "P20_008845_1894_XN_09N203W_4 1 0 1 -2103.971680 184.334869 \n",
+ "P20_008845_1894_XN_09N203W_4 1 0 1 -3226.000000 167.000000 \n",
+ "P20_008845_1894_XN_09N203W_5 1 0 1 -1813.285278 3.928788 \n",
+ "P20_008845_1894_XN_09N203W_5 1 0 1 -2934.000000 -6.000000 \n",
+ "P20_008845_1894_XN_09N203W_5 1 0 0 -4165.228516 2507.076660 \n",
+ "P20_008845_1894_XN_09N203W_6 1 0 1 -2250.549072 1572.788086 \n",
+ "P20_008845_1894_XN_09N203W_6 1 0 3 -3378.000000 1532.000000 \n",
+ "P20_008845_1894_XN_09N203W_6 1 0 1 2367.587646 -1849.272095 \n",
+ "P20_008845_1894_XN_09N203W_6 1 0 1 2675.137207 -1730.558105 \n",
+ "P20_008845_1894_XN_09N203W_7 1 0 1 -2249.931396 1712.971436 \n",
+ "P20_008845_1894_XN_09N203W_7 1 0 3 -3378.000000 1672.000000 \n",
+ "P20_008845_1894_XN_09N203W_7 1 0 1 2367.507813 -1706.241821 \n",
+ "P20_008845_1894_XN_09N203W_7 1 0 1 2675.892578 -1596.527100 \n",
+ "P20_008845_1894_XN_09N203W_8 1 0 1 -1664.131104 1875.622437 \n",
+ "P20_008845_1894_XN_09N203W_8 1 0 3 -2793.000000 1837.000000 \n",
+ "P20_008845_1894_XN_09N203W_8 1 0 1 2953.094727 -1538.956543 \n",
+ "P20_008845_1894_XN_09N203W_8 1 0 1 3261.290039 -1439.873169 \n",
+ "P20_008845_1894_XN_09N203W_9 1 0 1 332.072113 -2045.152222 \n",
+ "P20_008845_1894_XN_09N203W_9 1 0 0 -780.000000 -1953.000000 \n",
+ "P20_008845_1894_XN_09N203W_9 1 0 0 -2866.689453 1031.006104 \n",
+ "P20_008845_1894_XN_09N203W_9 1 0 0 -2002.366211 824.604126 \n",
+ "\n",
+ " sig_l sig_s res_l res_s fid_x \\\n",
+ "10_8344_8845_4r 0.000000 0.000000 0.059914 0.930311 0.0 \n",
+ "10_8344_8845_4r 0.000000 0.000000 0.063678 0.661294 0.0 \n",
+ "10_8344_8845_4r 0.000000 0.000000 -0.063695 -1.055619 0.0 \n",
+ "10_8344_8845_4r 0.000000 0.000000 0.318779 -0.585138 0.0 \n",
+ "11_8344_8845_4r 0.000000 0.000000 -0.239226 -1.113693 0.0 \n",
+ "11_8344_8845_4r 0.000000 0.000000 -0.247307 -0.587299 0.0 \n",
+ "11_8344_8845_4r 0.000000 0.000000 0.239207 1.080801 0.0 \n",
+ "11_8344_8845_4r 0.000000 0.000000 0.672877 0.626316 0.0 \n",
+ "12_8344_8845_4r 0.000000 0.000000 0.073683 -0.237114 0.0 \n",
+ "12_8344_8845_4r 0.000000 0.000000 -0.251753 0.008170 0.0 \n",
+ "12_8344_8845_4r 0.000000 0.000000 0.133406 0.068024 0.0 \n",
+ "12_8344_8845_4r 0.000000 0.000000 0.487209 0.036399 0.0 \n",
+ "13_8344_8845_4r 0.000000 0.000000 -0.056901 -0.383259 0.0 \n",
+ "13_8344_8845_4r 0.000000 0.000000 -0.131265 -0.305102 0.0 \n",
+ "13_8344_8845_4r 0.000000 0.000000 0.019201 0.485372 0.0 \n",
+ "13_8344_8845_4r 0.000000 0.000000 0.568393 0.146750 0.0 \n",
+ "14_8344_8845_4r 0.000000 0.000000 0.268449 1.000455 0.0 \n",
+ "14_8344_8845_4r 0.000000 0.000000 0.389532 -0.299258 0.0 \n",
+ "14_8344_8845_4r 0.000000 0.000000 -0.469784 0.023018 0.0 \n",
+ "14_8344_8845_4r 0.000000 0.000000 0.289877 -0.829150 0.0 \n",
+ "15_8344_8845_4r_mt_z 0.000000 0.000000 0.466921 1.168062 0.0 \n",
+ "15_8344_8845_4r_mt_z 0.000000 0.000000 -0.170850 0.372485 0.0 \n",
+ "15_8344_8845_4r_mt_z 0.000000 0.000000 -0.072035 -0.832436 0.0 \n",
+ "15_8344_8845_4r_mt_z 0.000000 0.000000 0.274030 -0.850105 0.0 \n",
+ "16_8344_8845_4r 0.000000 0.000000 -0.028473 -0.238015 0.0 \n",
+ "16_8344_8845_4r 0.000000 0.000000 0.165338 -0.327408 0.0 \n",
+ "16_8344_8845_4r 0.000000 0.000000 -0.007916 0.463232 0.0 \n",
+ "16_8344_8845_4r 0.000000 0.000000 0.328896 0.019526 0.0 \n",
+ "17_8344_8845_2r_mt_z 0.000000 0.000000 0.010045 0.244017 0.0 \n",
+ "17_8344_8845_2r_mt_z 0.000000 0.000000 -0.023370 -0.248643 0.0 \n",
+ "... ... ... ... ... ... \n",
+ "P20_008845_1894_XN_09N203W_18 0.022017 0.022017 0.734173 -0.074947 0.0 \n",
+ "P20_008845_1894_XN_09N203W_2 0.144273 0.144273 -0.061422 -0.263554 0.0 \n",
+ "P20_008845_1894_XN_09N203W_2 0.000000 0.000000 -0.073150 -0.271173 0.0 \n",
+ "P20_008845_1894_XN_09N203W_2 0.275627 0.275627 0.164525 0.349547 0.0 \n",
+ "P20_008845_1894_XN_09N203W_2 0.091086 0.091086 0.132736 0.239899 0.0 \n",
+ "P20_008845_1894_XN_09N203W_3 0.062058 0.062058 -0.198696 -0.105860 0.0 \n",
+ "P20_008845_1894_XN_09N203W_3 0.000000 0.000000 -0.088180 -0.049593 0.0 \n",
+ "P20_008845_1894_XN_09N203W_3 0.132941 0.132941 0.257260 0.117777 0.0 \n",
+ "P20_008845_1894_XN_09N203W_3 0.185544 0.185544 0.199301 0.065777 0.0 \n",
+ "P20_008845_1894_XN_09N203W_4 0.146920 0.146920 0.182557 -0.000683 0.0 \n",
+ "P20_008845_1894_XN_09N203W_4 0.000000 0.000000 -0.195774 -0.000066 0.0 \n",
+ "P20_008845_1894_XN_09N203W_5 0.019047 0.019047 -0.235171 -0.174079 0.0 \n",
+ "P20_008845_1894_XN_09N203W_5 0.000000 0.000000 -0.049005 0.437424 0.0 \n",
+ "P20_008845_1894_XN_09N203W_5 0.000000 0.000000 0.447100 -0.219764 0.0 \n",
+ "P20_008845_1894_XN_09N203W_6 0.025015 0.025015 -0.146990 0.192274 0.0 \n",
+ "P20_008845_1894_XN_09N203W_6 0.000000 0.000000 0.026541 0.063346 0.0 \n",
+ "P20_008845_1894_XN_09N203W_6 0.025110 0.025110 -0.040890 -0.136247 0.0 \n",
+ "P20_008845_1894_XN_09N203W_6 0.064276 0.064276 0.546035 -0.185103 0.0 \n",
+ "P20_008845_1894_XN_09N203W_7 0.010596 0.010596 -0.102909 0.079389 0.0 \n",
+ "P20_008845_1894_XN_09N203W_7 0.000000 0.000000 0.055623 0.013456 0.0 \n",
+ "P20_008845_1894_XN_09N203W_7 0.047435 0.047435 -0.157931 -0.042074 0.0 \n",
+ "P20_008845_1894_XN_09N203W_7 0.058218 0.058218 0.589561 -0.111421 0.0 \n",
+ "P20_008845_1894_XN_09N203W_8 0.017266 0.017266 -0.239689 -0.402891 0.0 \n",
+ "P20_008845_1894_XN_09N203W_8 0.000000 0.000000 -0.117483 -0.229337 0.0 \n",
+ "P20_008845_1894_XN_09N203W_8 0.023364 0.023364 0.177406 0.403161 0.0 \n",
+ "P20_008845_1894_XN_09N203W_8 0.023538 0.023538 0.575510 0.180705 0.0 \n",
+ "P20_008845_1894_XN_09N203W_9 0.232865 0.232865 0.095817 -0.483899 0.0 \n",
+ "P20_008845_1894_XN_09N203W_9 0.000000 0.000000 0.298743 -0.045963 0.0 \n",
+ "P20_008845_1894_XN_09N203W_9 0.000000 0.000000 -0.038437 0.471309 0.0 \n",
+ "P20_008845_1894_XN_09N203W_9 0.397616 0.397616 -0.141172 0.105223 0.0 \n",
+ "\n",
+ " ... known lat_Y_North long_X_East \\\n",
+ "10_8344_8845_4r ... 0 0.159378 2.724649 \n",
+ "10_8344_8845_4r ... 0 0.159378 2.724649 \n",
+ "10_8344_8845_4r ... 0 0.159378 2.724649 \n",
+ "10_8344_8845_4r ... 0 0.159378 2.724649 \n",
+ "11_8344_8845_4r ... 0 0.164905 2.721815 \n",
+ "11_8344_8845_4r ... 0 0.164905 2.721815 \n",
+ "11_8344_8845_4r ... 0 0.164905 2.721815 \n",
+ "11_8344_8845_4r ... 0 0.164905 2.721815 \n",
+ "12_8344_8845_4r ... 0 0.164949 2.724076 \n",
+ "12_8344_8845_4r ... 0 0.164949 2.724076 \n",
+ "12_8344_8845_4r ... 0 0.164949 2.724076 \n",
+ "12_8344_8845_4r ... 0 0.164949 2.724076 \n",
+ "13_8344_8845_4r ... 0 0.161840 2.723059 \n",
+ "13_8344_8845_4r ... 0 0.161840 2.723059 \n",
+ "13_8344_8845_4r ... 0 0.161840 2.723059 \n",
+ "13_8344_8845_4r ... 0 0.161840 2.723059 \n",
+ "14_8344_8845_4r ... 0 0.170415 2.720880 \n",
+ "14_8344_8845_4r ... 0 0.170415 2.720880 \n",
+ "14_8344_8845_4r ... 0 0.170415 2.720880 \n",
+ "14_8344_8845_4r ... 0 0.170415 2.720880 \n",
+ "15_8344_8845_4r_mt_z ... 1 0.170723 2.723266 \n",
+ "15_8344_8845_4r_mt_z ... 1 0.170723 2.723266 \n",
+ "15_8344_8845_4r_mt_z ... 1 0.170723 2.723266 \n",
+ "15_8344_8845_4r_mt_z ... 1 0.170723 2.723266 \n",
+ "16_8344_8845_4r ... 0 0.167682 2.722359 \n",
+ "16_8344_8845_4r ... 0 0.167682 2.722359 \n",
+ "16_8344_8845_4r ... 0 0.167682 2.722359 \n",
+ "16_8344_8845_4r ... 0 0.167682 2.722359 \n",
+ "17_8344_8845_2r_mt_z ... 1 0.161092 2.725426 \n",
+ "17_8344_8845_2r_mt_z ... 1 0.161092 2.725426 \n",
+ "... ... ... ... ... \n",
+ "P20_008845_1894_XN_09N203W_18 ... 0 0.171455 2.722281 \n",
+ "P20_008845_1894_XN_09N203W_2 ... 0 0.161184 2.729925 \n",
+ "P20_008845_1894_XN_09N203W_2 ... 0 0.161184 2.729925 \n",
+ "P20_008845_1894_XN_09N203W_2 ... 0 0.161184 2.729925 \n",
+ "P20_008845_1894_XN_09N203W_2 ... 0 0.161184 2.729925 \n",
+ "P20_008845_1894_XN_09N203W_3 ... 0 0.161992 2.729563 \n",
+ "P20_008845_1894_XN_09N203W_3 ... 0 0.161992 2.729563 \n",
+ "P20_008845_1894_XN_09N203W_3 ... 0 0.161992 2.729563 \n",
+ "P20_008845_1894_XN_09N203W_3 ... 0 0.161992 2.729563 \n",
+ "P20_008845_1894_XN_09N203W_4 ... 0 0.160994 2.725912 \n",
+ "P20_008845_1894_XN_09N203W_4 ... 0 0.160994 2.725912 \n",
+ "P20_008845_1894_XN_09N203W_5 ... 0 0.161544 2.726149 \n",
+ "P20_008845_1894_XN_09N203W_5 ... 0 0.161544 2.726149 \n",
+ "P20_008845_1894_XN_09N203W_5 ... 0 0.161544 2.726149 \n",
+ "P20_008845_1894_XN_09N203W_6 ... 0 0.160447 2.723639 \n",
+ "P20_008845_1894_XN_09N203W_6 ... 0 0.160447 2.723639 \n",
+ "P20_008845_1894_XN_09N203W_6 ... 0 0.160447 2.723639 \n",
+ "P20_008845_1894_XN_09N203W_6 ... 0 0.160447 2.723639 \n",
+ "P20_008845_1894_XN_09N203W_7 ... 0 0.160419 2.723407 \n",
+ "P20_008845_1894_XN_09N203W_7 ... 0 0.160419 2.723407 \n",
+ "P20_008845_1894_XN_09N203W_7 ... 0 0.160419 2.723407 \n",
+ "P20_008845_1894_XN_09N203W_7 ... 0 0.160419 2.723407 \n",
+ "P20_008845_1894_XN_09N203W_8 ... 0 0.161417 2.723012 \n",
+ "P20_008845_1894_XN_09N203W_8 ... 0 0.161417 2.723012 \n",
+ "P20_008845_1894_XN_09N203W_8 ... 0 0.161417 2.723012 \n",
+ "P20_008845_1894_XN_09N203W_8 ... 0 0.161417 2.723012 \n",
+ "P20_008845_1894_XN_09N203W_9 ... 0 0.165741 2.729053 \n",
+ "P20_008845_1894_XN_09N203W_9 ... 0 0.165741 2.729053 \n",
+ "P20_008845_1894_XN_09N203W_9 ... 0 0.165741 2.729053 \n",
+ "P20_008845_1894_XN_09N203W_9 ... 0 0.165741 2.729053 \n",
+ "\n",
+ " ht sig0 sig1 sig2 res0 \\\n",
+ "10_8344_8845_4r -2523.828227 0.0 0.0 25.000000 18.301328 \n",
+ "10_8344_8845_4r -2523.828227 0.0 0.0 25.000000 18.301328 \n",
+ "10_8344_8845_4r -2523.828227 0.0 0.0 25.000000 18.301328 \n",
+ "10_8344_8845_4r -2523.828227 0.0 0.0 25.000000 18.301328 \n",
+ "11_8344_8845_4r -2445.237027 0.0 0.0 30.000000 -22.046575 \n",
+ "11_8344_8845_4r -2445.237027 0.0 0.0 30.000000 -22.046575 \n",
+ "11_8344_8845_4r -2445.237027 0.0 0.0 30.000000 -22.046575 \n",
+ "11_8344_8845_4r -2445.237027 0.0 0.0 30.000000 -22.046575 \n",
+ "12_8344_8845_4r -2606.935163 0.0 0.0 100.000000 -7.549561 \n",
+ "12_8344_8845_4r -2606.935163 0.0 0.0 100.000000 -7.549561 \n",
+ "12_8344_8845_4r -2606.935163 0.0 0.0 100.000000 -7.549561 \n",
+ "12_8344_8845_4r -2606.935163 0.0 0.0 100.000000 -7.549561 \n",
+ "13_8344_8845_4r -2551.901554 0.0 0.0 4.536068 -1.643694 \n",
+ "13_8344_8845_4r -2551.901554 0.0 0.0 4.536068 -1.643694 \n",
+ "13_8344_8845_4r -2551.901554 0.0 0.0 4.536068 -1.643694 \n",
+ "13_8344_8845_4r -2551.901554 0.0 0.0 4.536068 -1.643694 \n",
+ "14_8344_8845_4r -2505.953426 0.0 0.0 5.000000 -50.069808 \n",
+ "14_8344_8845_4r -2505.953426 0.0 0.0 5.000000 -50.069808 \n",
+ "14_8344_8845_4r -2505.953426 0.0 0.0 5.000000 -50.069808 \n",
+ "14_8344_8845_4r -2505.953426 0.0 0.0 5.000000 -50.069808 \n",
+ "15_8344_8845_4r_mt_z -2502.470000 0.0 0.0 5.000000 -34.817656 \n",
+ "15_8344_8845_4r_mt_z -2502.470000 0.0 0.0 5.000000 -34.817656 \n",
+ "15_8344_8845_4r_mt_z -2502.470000 0.0 0.0 5.000000 -34.817656 \n",
+ "15_8344_8845_4r_mt_z -2502.470000 0.0 0.0 5.000000 -34.817656 \n",
+ "16_8344_8845_4r -2558.312931 0.0 0.0 1.707214 -29.465246 \n",
+ "16_8344_8845_4r -2558.312931 0.0 0.0 1.707214 -29.465246 \n",
+ "16_8344_8845_4r -2558.312931 0.0 0.0 1.707214 -29.465246 \n",
+ "16_8344_8845_4r -2558.312931 0.0 0.0 1.707214 -29.465246 \n",
+ "17_8344_8845_2r_mt_z -2590.130000 0.0 0.0 2.000000 -10.779830 \n",
+ "17_8344_8845_2r_mt_z -2590.130000 0.0 0.0 2.000000 -10.779830 \n",
+ "... ... ... ... ... ... \n",
+ "P20_008845_1894_XN_09N203W_18 -2514.431453 0.0 0.0 0.000000 -38.388316 \n",
+ "P20_008845_1894_XN_09N203W_2 -2535.712262 0.0 0.0 0.000000 42.067289 \n",
+ "P20_008845_1894_XN_09N203W_2 -2535.712262 0.0 0.0 0.000000 42.067289 \n",
+ "P20_008845_1894_XN_09N203W_2 -2535.712262 0.0 0.0 0.000000 42.067289 \n",
+ "P20_008845_1894_XN_09N203W_2 -2535.712262 0.0 0.0 0.000000 42.067289 \n",
+ "P20_008845_1894_XN_09N203W_3 -2525.215515 0.0 0.0 0.000000 36.033823 \n",
+ "P20_008845_1894_XN_09N203W_3 -2525.215515 0.0 0.0 0.000000 36.033823 \n",
+ "P20_008845_1894_XN_09N203W_3 -2525.215515 0.0 0.0 0.000000 36.033823 \n",
+ "P20_008845_1894_XN_09N203W_3 -2525.215515 0.0 0.0 0.000000 36.033823 \n",
+ "P20_008845_1894_XN_09N203W_4 -2562.446851 0.0 0.0 0.000000 -7.233876 \n",
+ "P20_008845_1894_XN_09N203W_4 -2562.446851 0.0 0.0 0.000000 -7.233876 \n",
+ "P20_008845_1894_XN_09N203W_5 -2560.812028 0.0 0.0 0.000000 7.498790 \n",
+ "P20_008845_1894_XN_09N203W_5 -2560.812028 0.0 0.0 0.000000 7.498790 \n",
+ "P20_008845_1894_XN_09N203W_5 -2560.812028 0.0 0.0 0.000000 7.498790 \n",
+ "P20_008845_1894_XN_09N203W_6 -2605.266130 0.0 0.0 0.000000 7.687475 \n",
+ "P20_008845_1894_XN_09N203W_6 -2605.266130 0.0 0.0 0.000000 7.687475 \n",
+ "P20_008845_1894_XN_09N203W_6 -2605.266130 0.0 0.0 0.000000 7.687475 \n",
+ "P20_008845_1894_XN_09N203W_6 -2605.266130 0.0 0.0 0.000000 7.687475 \n",
+ "P20_008845_1894_XN_09N203W_7 -2619.498291 0.0 0.0 0.000000 6.358447 \n",
+ "P20_008845_1894_XN_09N203W_7 -2619.498291 0.0 0.0 0.000000 6.358447 \n",
+ "P20_008845_1894_XN_09N203W_7 -2619.498291 0.0 0.0 0.000000 6.358447 \n",
+ "P20_008845_1894_XN_09N203W_7 -2619.498291 0.0 0.0 0.000000 6.358447 \n",
+ "P20_008845_1894_XN_09N203W_8 -2608.028730 0.0 0.0 0.000000 -0.219662 \n",
+ "P20_008845_1894_XN_09N203W_8 -2608.028730 0.0 0.0 0.000000 -0.219662 \n",
+ "P20_008845_1894_XN_09N203W_8 -2608.028730 0.0 0.0 0.000000 -0.219662 \n",
+ "P20_008845_1894_XN_09N203W_8 -2608.028730 0.0 0.0 0.000000 -0.219662 \n",
+ "P20_008845_1894_XN_09N203W_9 -2586.686862 0.0 0.0 0.000000 15.191825 \n",
+ "P20_008845_1894_XN_09N203W_9 -2586.686862 0.0 0.0 0.000000 15.191825 \n",
+ "P20_008845_1894_XN_09N203W_9 -2586.686862 0.0 0.0 0.000000 15.191825 \n",
+ "P20_008845_1894_XN_09N203W_9 -2586.686862 0.0 0.0 0.000000 15.191825 \n",
+ "\n",
+ " res1 res2 \n",
+ "10_8344_8845_4r 44.206259 416.201741 \n",
+ "10_8344_8845_4r 44.206259 416.201741 \n",
+ "10_8344_8845_4r 44.206259 416.201741 \n",
+ "10_8344_8845_4r 44.206259 416.201741 \n",
+ "11_8344_8845_4r 103.403228 173.210013 \n",
+ "11_8344_8845_4r 103.403228 173.210013 \n",
+ "11_8344_8845_4r 103.403228 173.210013 \n",
+ "11_8344_8845_4r 103.403228 173.210013 \n",
+ "12_8344_8845_4r 93.170584 180.058858 \n",
+ "12_8344_8845_4r 93.170584 180.058858 \n",
+ "12_8344_8845_4r 93.170584 180.058858 \n",
+ "12_8344_8845_4r 93.170584 180.058858 \n",
+ "13_8344_8845_4r 74.244153 319.554535 \n",
+ "13_8344_8845_4r 74.244153 319.554535 \n",
+ "13_8344_8845_4r 74.244153 319.554535 \n",
+ "13_8344_8845_4r 74.244153 319.554535 \n",
+ "14_8344_8845_4r 144.316524 -164.540707 \n",
+ "14_8344_8845_4r 144.316524 -164.540707 \n",
+ "14_8344_8845_4r 144.316524 -164.540707 \n",
+ "14_8344_8845_4r 144.316524 -164.540707 \n",
+ "15_8344_8845_4r_mt_z 145.864550 0.214023 \n",
+ "15_8344_8845_4r_mt_z 145.864550 0.214023 \n",
+ "15_8344_8845_4r_mt_z 145.864550 0.214023 \n",
+ "15_8344_8845_4r_mt_z 145.864550 0.214023 \n",
+ "16_8344_8845_4r 121.908506 19.874949 \n",
+ "16_8344_8845_4r 121.908506 19.874949 \n",
+ "16_8344_8845_4r 121.908506 19.874949 \n",
+ "16_8344_8845_4r 121.908506 19.874949 \n",
+ "17_8344_8845_2r_mt_z 38.761214 0.236077 \n",
+ "17_8344_8845_2r_mt_z 38.761214 0.236077 \n",
+ "... ... ... \n",
+ "P20_008845_1894_XN_09N203W_18 195.050373 -391.289238 \n",
+ "P20_008845_1894_XN_09N203W_2 -12.658160 627.742625 \n",
+ "P20_008845_1894_XN_09N203W_2 -12.658160 627.742625 \n",
+ "P20_008845_1894_XN_09N203W_2 -12.658160 627.742625 \n",
+ "P20_008845_1894_XN_09N203W_2 -12.658160 627.742625 \n",
+ "P20_008845_1894_XN_09N203W_3 -4.170818 616.480181 \n",
+ "P20_008845_1894_XN_09N203W_3 -4.170818 616.480181 \n",
+ "P20_008845_1894_XN_09N203W_3 -4.170818 616.480181 \n",
+ "P20_008845_1894_XN_09N203W_3 -4.170818 616.480181 \n",
+ "P20_008845_1894_XN_09N203W_4 33.324946 574.887719 \n",
+ "P20_008845_1894_XN_09N203W_4 33.324946 574.887719 \n",
+ "P20_008845_1894_XN_09N203W_5 33.109007 581.149470 \n",
+ "P20_008845_1894_XN_09N203W_5 33.109007 581.149470 \n",
+ "P20_008845_1894_XN_09N203W_5 33.109007 581.149470 \n",
+ "P20_008845_1894_XN_09N203W_6 59.859663 375.680663 \n",
+ "P20_008845_1894_XN_09N203W_6 59.859663 375.680663 \n",
+ "P20_008845_1894_XN_09N203W_6 59.859663 375.680663 \n",
+ "P20_008845_1894_XN_09N203W_6 59.859663 375.680663 \n",
+ "P20_008845_1894_XN_09N203W_7 61.278242 375.840988 \n",
+ "P20_008845_1894_XN_09N203W_7 61.278242 375.840988 \n",
+ "P20_008845_1894_XN_09N203W_7 61.278242 375.840988 \n",
+ "P20_008845_1894_XN_09N203W_7 61.278242 375.840988 \n",
+ "P20_008845_1894_XN_09N203W_8 71.414267 336.409684 \n",
+ "P20_008845_1894_XN_09N203W_8 71.414267 336.409684 \n",
+ "P20_008845_1894_XN_09N203W_8 71.414267 336.409684 \n",
+ "P20_008845_1894_XN_09N203W_8 71.414267 336.409684 \n",
+ "P20_008845_1894_XN_09N203W_9 22.094037 544.874936 \n",
+ "P20_008845_1894_XN_09N203W_9 22.094037 544.874936 \n",
+ "P20_008845_1894_XN_09N203W_9 22.094037 544.874936 \n",
+ "P20_008845_1894_XN_09N203W_9 22.094037 544.874936 \n",
+ "\n",
+ "[919 rows x 23 columns]"
+ ]
+ },
+ "execution_count": 22,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "atf_dict = read_atf(get_path('CTX_Athabasca_Middle_step0.atf'))\n",
+ "\n",
+ "gpf_file = os.path.join(atf_dict['PATH'], atf_dict['GP_FILE']);\n",
+ "ipf_list = [os.path.join(atf_dict['PATH'], i) for i in atf_dict['IMAGE_IPF']]\n",
+ "\n",
+ "gpf_df = read_gpf(gpf_file).set_index('point_id')\n",
+ "ipf_df = read_ipf(ipf_list).set_index('pt_id')\n",
+ "\n",
+ "point_diff = ipf_df.index.difference(gpf_df.index)\n",
+ "\n",
+ "if len(point_diff) != 0:\n",
+ " warnings.warn(\"The following points found in ipf files missing from gpf file: \\n\\n{}. \\\n",
+ " \\n\\nContinuing, but these points will be missing from the control network\".format(list(point_diff)))\n",
+ "\n",
+ "new_df = ipf_df.merge(gpf_df, left_index=True, right_index=True)\n",
+ "new_df"
+ ]
},
{
"cell_type": "code",
@@ -165,7 +2042,7 @@
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
- "version": "3.6.4"
+ "version": "3.6.3"
}
},
"nbformat": 4,
diff --git a/plio/examples/SocetGXP/InSightE08_XW.csv b/plio/examples/SocetGXP/InSightE08_XW.csv
index 187caf9d7778cf67c318bad830f95ea29b846eac..35f01d10c3484923367707aefd60081f7a895ba1 100644
--- a/plio/examples/SocetGXP/InSightE08_XW.csv
+++ b/plio/examples/SocetGXP/InSightE08_XW.csv
@@ -1,4 +1,4 @@
-point_id,stat,known,lat_Y_North,long_X_East,ht,sigma0,sigma1,sigma2,res0,res1,res2
+point_id,stat,known,lat_Y_North,long_X_East,ht,sig0,sig1,sig2,res0,res1,res2
1,1,1,0.08538133869186999,2.38013146162178,-2697.2374469465,0.0,0.0,1.707214,16.348041000000002,-13.917623,-0.151001
2,1,1,0.08508542650082002,2.3808135855714703,-2702.0,0.0,0.0,1.0,15.983285999999998,-16.869508,0.103341
3,1,1,0.08626480095809,2.3870838392611,-2714.16003756809,0.0,0.0,3.828854,12.831839,-34.360639,-2.2058459999999998
diff --git a/plio/examples/SocetSet/P20_008845_1894_XN_09N203W.csv b/plio/examples/SocetSet/P20_008845_1894_XN_09N203W.csv
new file mode 100644
index 0000000000000000000000000000000000000000..a137bf0f982373fb8d94c1a7d97a4e439e8027c2
--- /dev/null
+++ b/plio/examples/SocetSet/P20_008845_1894_XN_09N203W.csv
@@ -0,0 +1,122 @@
+pt_id,val,fid_val,no_obs,l.,s.,sig_l,sig_s,res_l,res_s,fid_x,fid_y
+1_8344_8845_4r,1,0,0,-4058.982422,-2318.0107420000004,0.0,0.0,-0.062556,-0.21471300000000001,0.0,0.0
+2_8344_8845_4r,1,0,0,-3969.065186,-606.849243,0.0,0.0,0.22866,0.105249,0.0,0.0
+3_8344_8845_4r_mt_z,1,0,0,-1019.7390140000001,-2300.877197,0.0,0.0,-0.025129,-0.002447,0.0,0.0
+4_8344_8845_4r_mt_z,1,0,0,-1037.0991210000002,-548.180237,0.0,0.0,-0.0007559999999999999,0.227752,0.0,0.0
+5_8344_8845_4r,1,0,0,2438.9841309999997,-2304.8435059999997,0.0,0.0,0.062022,-0.10997699999999999,0.0,0.0
+6_8344_8845_4r,1,0,0,2397.826904,-562.432861,0.0,0.0,-0.07267799999999999,-0.056103999999999994,0.0,0.0
+7_8344_8845_4r_mt_z,1,0,0,-2510.927734,-1267.364868,0.0,0.0,0.21525,-0.040427,0.0,0.0
+8_8344-8845_4r_mt_z,1,0,0,840.825317,-1028.345337,0.0,0.0,-0.24192199999999997,0.382924,0.0,0.0
+9_8344_8845_4r,1,0,0,-4070.962158,2465.817139,0.0,0.0,-0.140837,0.00598,0.0,0.0
+10_8344_8845_4r,1,0,0,-4044.69751,1008.950928,0.0,0.0,0.063678,0.6612939999999999,0.0,0.0
+11_8344_8845_4r,1,0,0,-761.216064,2303.787109,0.0,0.0,-0.24730700000000003,-0.587299,0.0,0.0
+12_8344_8845_4r,1,0,0,-889.364441,966.5339970000001,0.0,0.0,-0.251753,0.00817,0.0,0.0
+13_8344_8845_4r,1,0,0,-2559.871338,1777.522827,0.0,0.0,-0.13126500000000002,-0.305102,0.0,0.0
+14_8344_8845_4r,1,0,0,2385.27832,2476.032227,0.0,0.0,0.389532,-0.29925799999999997,0.0,0.0
+15_8344_8845_4r_mt_z,1,0,0,2395.869385,1038.165405,0.0,0.0,-0.17085,0.372485,0.0,0.0
+16_8344_8845_4r,1,0,0,756.099792,1785.4947510000002,0.0,0.0,0.165338,-0.327408,0.0,0.0
+17_8344_8845_2r_mt_z,1,0,0,-3138.269531,442.515503,0.0,0.0,-0.02337,-0.24864299999999998,0.0,0.0
+18_8344_8845_2r,1,0,0,-1773.8649899999998,354.084259,0.0,0.0,-0.002223,-0.000684,0.0,0.0
+19_8344_8845_2r_mt_z,1,0,0,-92.482826,723.305237,0.0,0.0,-0.214498,0.14519200000000002,0.0,0.0
+20_8344_8845_2r_mt_z,1,0,0,1434.079712,742.064026,0.0,0.0,-0.058132,-0.171782,0.0,0.0
+21_8344_8845_4r_xyz,1,0,0,-838.9910279999999,-614.524109,0.0,0.0,-0.131313,0.9971200000000001,0.0,0.0
+22_8344_8845_4r_mt_z,1,0,0,-3994.073975,-2160.606445,0.0,0.0,0.21777399999999997,-0.141372,0.0,0.0
+23_8344_8845_4r_mt_z,1,0,0,-3851.86499,2088.046875,0.0,0.0,-0.07844,-0.22165,0.0,0.0
+24_8344_8845_4r_mt_z,1,0,0,2343.634033,-2073.817871,0.0,0.0,-0.27434400000000003,0.06779099999999999,0.0,0.0
+25_8344_8845_4r_mt_z,1,0,0,2350.053711,2168.792236,0.0,0.0,-0.311081,-0.152226,0.0,0.0
+26_8344_8845_4r_mt_z,1,0,0,-1844.550781,-996.0164789999999,0.0,0.0,-0.108472,0.122981,0.0,0.0
+27_8344_8845_2r_mt_z,1,0,0,-1618.0053710000002,700.699158,0.0,0.0,-0.020846,0.073456,0.0,0.0
+28_8344_8845_4r_mt_z,1,0,0,521.921509,-1602.4414060000001,0.0,0.0,-0.009861,0.363227,0.0,0.0
+29_8344_8845_4r_mt_z,1,0,0,820.691284,1469.345459,0.0,0.0,-0.664959,-0.06624400000000001,0.0,0.0
+30_8344_8845_3r_mt_z,1,0,0,-326.405365,-283.346985,0.0,0.0,-0.20681100000000002,0.180573,0.0,0.0
+61_8344_8845_4r_mt_z,1,0,0,-4067.466064,-2193.302002,0.0,0.0,-0.017458,-0.229608,0.0,0.0
+62_8344_8845_4r_mt_z,1,0,0,2040.637939,2052.019043,0.0,0.0,0.187054,-0.19953900000000002,0.0,0.0
+63_8344_8845_4r_mt_z,1,0,0,2267.361328,-1402.6448970000001,0.0,0.0,-0.25967399999999996,-0.350489,0.0,0.0
+65_8344_8845_4r_mt_z,1,0,0,-675.820618,1659.5936279999999,0.0,0.0,0.285702,-0.673177,0.0,0.0
+P19_008344_1894_XN_09N203W_3,1,0,1,-3543.180176,-2052.497559,0.09999,0.09999,-0.271723,0.046268000000000004,0.0,0.0
+P19_008344_1894_XN_09N203W_5,1,0,1,-3693.414551,118.131378,0.191225,0.191225,-0.007534999999999999,3e-06,0.0,0.0
+P19_008344_1894_XN_09N203W_8,1,0,1,-3529.115479,1634.99353,0.010654,0.010654,-0.11299000000000001,0.190108,0.0,0.0
+P19_008344_1894_XN_09N203W_9,1,0,1,-1530.085083,-2021.120728,0.025998,0.025998,-0.28780100000000003,-0.05082,0.0,0.0
+P19_008344_1894_XN_09N203W_10,1,0,1,-1163.659546,-2110.94873,0.054246,0.054246,-0.209774,-0.028914999999999996,0.0,0.0
+P19_008344_1894_XN_09N203W_11,1,0,1,-649.603027,-291.606293,0.144968,0.144968,0.179422,0.48035,0.0,0.0
+P19_008344_1894_XN_09N203W_12,1,0,1,-1879.1154789999998,1650.314575,0.044342,0.044342,0.194544,-0.551285,0.0,0.0
+P19_008344_1894_XN_09N203W_13,1,0,1,1195.076782,-2043.631714,0.081316,0.081316,-0.184309,0.153021,0.0,0.0
+P19_008344_1894_XN_09N203W_14,1,0,1,1826.979736,-2180.584961,0.185789,0.185789,0.208304,-0.026378,0.0,0.0
+P19_008344_1894_XN_09N203W_15,1,0,1,1833.592407,-1740.516968,0.23812199999999997,0.23812199999999997,0.652892,-0.509659,0.0,0.0
+P19_008344_1894_XN_09N203W_16,1,0,1,1370.180664,-273.589996,0.11324100000000001,0.11324100000000001,0.296132,-0.060485000000000004,0.0,0.0
+P19_008344_1894_XN_09N203W_17,1,0,1,1198.482544,-429.786133,0.055307,0.055307,-0.091967,0.10811400000000002,0.0,0.0
+P19_008344_1894_XN_09N203W_18,1,0,1,1516.067993,1830.024048,0.009269,0.009269,-0.002093,-0.264173,0.0,0.0
+P19_008344_1894_XN_09N203W_19,1,0,1,1338.967651,1679.986084,0.004901,0.004901,-0.030115,-0.306342,0.0,0.0
+P19_008344_1894_XN_09N203W_20,1,0,1,1350.758423,1952.026978,0.019969999999999998,0.019969999999999998,-0.011408,-0.409686,0.0,0.0
+P20_008845_1894_XN_09N203W_1,1,0,0,-3224.0,-2003.0,0.0,0.0,0.055124,0.018517,0.0,0.0
+P20_008845_1894_XN_09N203W_2,1,0,0,-3389.0,-2155.0,0.0,0.0,-0.07315,-0.271173,0.0,0.0
+P20_008845_1894_XN_09N203W_3,1,0,0,-2913.0,-2000.0,0.0,0.0,-0.08818,-0.049593,0.0,0.0
+P20_008845_1894_XN_09N203W_4,1,0,1,-3226.0,167.0,0.0,0.0,-0.195774,-6.6e-05,0.0,0.0
+P20_008845_1894_XN_09N203W_5,1,0,1,-2934.0,-6.0,0.0,0.0,-0.049005,0.437424,0.0,0.0
+P20_008845_1894_XN_09N203W_6,1,0,3,-3378.0,1532.0,0.0,0.0,0.026541000000000002,0.063346,0.0,0.0
+P20_008845_1894_XN_09N203W_7,1,0,3,-3378.0,1672.0,0.0,0.0,0.055623,0.013456000000000001,0.0,0.0
+P20_008845_1894_XN_09N203W_8,1,0,3,-2793.0,1837.0,0.0,0.0,-0.11748299999999999,-0.22933699999999999,0.0,0.0
+P20_008845_1894_XN_09N203W_9,1,0,0,-780.0,-1953.0,0.0,0.0,0.298743,-0.045963,0.0,0.0
+P20_008845_1894_XN_09N203W_10,1,0,2,423.0,-1995.0,0.0,0.0,0.067585,0.401377,0.0,0.0
+P20_008845_1894_XN_09N203W_11,1,0,1,-21.0,196.0,0.0,0.0,0.06955,-0.000458,0.0,0.0
+P20_008845_1894_XN_09N203W_12,1,0,1,2885.0,-2124.0,0.0,0.0,0.25413600000000003,-0.273798,0.0,0.0
+P20_008845_1894_XN_09N203W_13,1,0,1,2746.0,-2253.0,0.0,0.0,-0.407167,-0.176964,0.0,0.0
+P20_008845_1894_XN_09N203W_14,1,0,0,2728.0,-1980.0,0.0,0.0,-0.885848,0.137668,0.0,0.0
+P20_008845_1894_XN_09N203W_15,0,0,0,2942.0,58.0,0.0,0.0,0.0,0.0,0.0,0.0
+P20_008845_1894_XN_09N203W_16,1,0,2,2865.0,2024.0,0.0,0.0,-0.053859000000000004,0.035981,0.0,0.0
+P20_008845_1894_XN_09N203W_17,1,0,2,3175.0,1974.0,0.0,0.0,-0.064058,0.044374000000000004,0.0,0.0
+P20_008845_1894_XN_09N203W_18,1,0,2,2872.0,1565.0,0.0,0.0,0.256228,0.21886,0.0,0.0
+P03_002371_1888_XI_08N204W_49,1,0,1,-2629.182861,1085.034424,0.012761,0.012761,0.20840300000000003,-0.004282,0.0,0.0
+P03_002371_1888_XI_08N204W_54,1,0,1,-1095.4067380000001,1127.088013,0.274421,0.274421,0.019106,-0.331897,0.0,0.0
+P03_002371_1888_XI_08N204W_55,1,0,1,-941.496277,1247.825806,0.05563099999999999,0.05563099999999999,0.152919,-0.306029,0.0,0.0
+P03_002371_1888_XI_08N204W_56,1,0,1,-1098.9589839999999,1441.816895,0.11484200000000001,0.11484200000000001,-0.297223,-0.308728,0.0,0.0
+P03_002371_1888_XI_08N204W_63,1,0,1,119.94445800000001,1116.177612,0.023727,0.023727,-0.104057,-0.034587,0.0,0.0
+P03_002371_1888_XI_08N204W_64,1,0,1,-19.255381,1270.651611,0.007022,0.007022,-0.29496999999999995,-0.13664600000000002,0.0,0.0
+P03_002371_1888_XI_08N204W_65,1,0,1,-271.636139,1030.2692869999998,0.14131400000000002,0.14131400000000002,-0.007378,0.127183,0.0,0.0
+P03_002371_1888_XI_08N204W_71,1,0,1,1793.756836,1419.46814,0.070811,0.070811,0.008126999999999999,-0.059650999999999996,0.0,0.0
+P03_002371_1888_XI_08N204W_77,1,0,1,3770.185547,1056.4335939999999,0.051015,0.051015,-0.344149,0.7918810000000001,0.0,0.0
+P01_001540_1889_XI_08N204W_3,1,0,0,-819.994568,-476.940125,0.0,0.0,-0.010086,0.238933,0.0,0.0
+P01_001540_1889_XI_08N204W_4,1,0,0,-1038.7459720000002,-225.617737,0.0,0.0,-0.10251500000000001,0.293944,0.0,0.0
+P01_001540_1889_XI_08N204W_5,1,0,0,-1043.086548,-229.29058799999999,0.0,0.0,-0.103618,0.22217399999999998,0.0,0.0
+P01_001540_1889_XI_08N204W_6,1,0,0,-1016.6671140000001,-212.669418,0.0,0.0,0.098336,0.182504,0.0,0.0
+P01_001540_1889_XI_08N204W_7,1,0,0,-1089.852539,-179.61308300000002,0.0,0.0,-0.09901499999999999,0.24748699999999998,0.0,0.0
+P01_001540_1889_XI_08N204W_43,1,0,1,231.430023,1286.035034,0.179721,0.179721,0.043216000000000004,0.007601999999999999,0.0,0.0
+P01_001540_1889_XI_08N204W_44,1,0,1,53.3703,1130.200195,0.035373,0.035373,-0.239021,0.05197,0.0,0.0
+P01_001606_1897_XI_09N203W_20,1,0,1,-3490.634277,-1071.889893,0.074986,0.074986,-0.10951300000000001,0.027899,0.0,0.0
+P01_001606_1897_XI_09N203W_21,1,0,1,-3489.029297,-1344.598267,0.117221,0.117221,0.11984000000000002,-0.07396599999999999,0.0,0.0
+P01_001606_1897_XI_09N203W_26,1,0,1,-2843.921875,-1365.05249,0.066368,0.066368,0.053526,-0.129321,0.0,0.0
+P01_001606_1897_XI_09N203W_27,1,0,1,-2872.278076,-1055.845581,0.044944,0.044944,0.014887000000000001,-0.175307,0.0,0.0
+P01_001606_1897_XI_09N203W_32,1,0,1,-702.726318,-1309.760742,0.097152,0.097152,-0.12127,0.083388,0.0,0.0
+P01_001606_1897_XI_09N203W_33,1,0,1,-425.8797,-1332.264404,0.141272,0.141272,0.27812,0.018378,0.0,0.0
+P01_001606_1897_XI_09N203W_39,1,0,1,957.507324,-1314.205933,0.214219,0.214219,-0.08834600000000001,0.135983,0.0,0.0
+P01_001606_1897_XI_09N203W_40,1,0,1,1286.863525,-1325.946045,0.022701,0.022701,0.025672000000000004,-0.015907,0.0,0.0
+P01_001606_1897_XI_09N203W_41,1,0,1,1415.3436279999999,-1585.812134,0.062944,0.062944,0.17824700000000002,0.100219,0.0,0.0
+P01_001606_1897_XI_09N203W_42,1,0,1,-713.122681,-1899.2719730000001,0.10971099999999999,0.10971099999999999,-0.24755700000000003,0.240463,0.0,0.0
+P01_001606_1897_XI_09N203W_48,1,0,1,2357.769287,-1605.21167,0.142259,0.142259,-0.046652,-0.341859,0.0,0.0
+P01_001606_1897_XI_09N203W_49,1,0,1,2368.67749,-1310.2536619999998,0.084816,0.084816,0.16844800000000001,-0.178884,0.0,0.0
+P01_001606_1897_XI_09N203W_50,1,0,1,2938.904541,-1286.650635,0.268181,0.268181,0.22215,-0.28928200000000004,0.0,0.0
+P03_002226_1895_XI_09N203W_41,1,0,1,-2964.822021,-752.72998,0.130511,0.130511,0.19989500000000002,0.18193199999999998,0.0,0.0
+P03_002226_1895_XI_09N203W_42,1,0,1,-2676.941162,-549.997925,0.015522999999999999,0.015522999999999999,-0.189075,0.282047,0.0,0.0
+P03_002226_1895_XI_09N203W_43,1,0,1,-2503.68457,-1243.333618,0.091619,0.091619,-0.183131,-0.034143,0.0,0.0
+P03_002226_1895_XI_09N203W_46,1,0,1,-984.080139,-1301.115845,0.02843,0.02843,-0.10387,-0.014166,0.0,0.0
+P03_002226_1895_XI_09N203W_47,1,0,1,-1119.626953,-926.1320800000001,0.195907,0.195907,-0.217611,0.08741499999999999,0.0,0.0
+P03_002226_1895_XI_09N203W_53,1,0,1,3304.79834,-528.821411,0.090225,0.090225,-0.12703599999999998,-0.24331999999999998,0.0,0.0
+P03_002226_1895_XI_09N203W_54,1,0,1,3597.587646,-608.987488,0.18746300000000002,0.18746300000000002,0.13215,-0.402642,0.0,0.0
+P03_002226_1895_XI_09N203W_55,1,0,1,3277.619873,-1106.364136,0.141123,0.141123,-0.258301,-0.48746300000000004,0.0,0.0
+P03_002226_1895_XI_09N203W_64,1,0,1,-1488.9279789999998,306.37027,0.141123,0.141123,0.11354000000000002,0.000408,0.0,0.0
+P03_002226_1895_XI_09N203W_67,1,0,1,-1501.843384,294.863434,0.141123,0.141123,-0.038216,-0.001066,0.0,0.0
+P03_002226_1895_XI_09N203W_68,1,0,1,-1487.291992,301.690338,0.141123,0.141123,-0.16136,0.000601,0.0,0.0
+P03_002226_1895_XI_09N203W_71,1,0,1,-1476.8244630000002,421.30285599999996,0.141123,0.141123,0.045458,0.000436,0.0,0.0
+P03_002226_1895_XI_09N203W_66,1,0,1,-1473.926636,445.16494800000004,0.141123,0.141123,-0.039173,-6.1e-05,0.0,0.0
+P03_002226_1895_XI_09N203W_72,1,0,1,-1478.19104,398.059052,0.141123,0.141123,0.03951,-0.0005690000000000001,0.0,0.0
+P03_002226_1895_XI_09N203W_74,1,0,1,-1463.276489,270.004364,0.141123,0.141123,0.097608,0.000572,0.0,0.0
+P03_002226_1895_XI_09N203W_76,1,0,1,-1399.385376,210.815338,0.141123,0.141123,-0.021339,-0.000335,0.0,0.0
+P03_002226_1895_XI_09N203W_82,1,0,1,-1110.612061,-112.804382,0.141123,0.141123,-0.179249,0.431591,0.0,0.0
+P19_008344_1894_XN_09N203W_2,1,0,1,-3987.0656740000004,-1756.185791,0.141123,0.141123,0.018394,-0.141215,0.0,0.0
+P19_008344_1894_XN_09N203W_7,1,0,1,-3991.17334,1967.117554,0.141123,0.141123,0.242488,-0.04689,0.0,0.0
+P03_002371_1888_XI_08N204W_76,1,0,1,4066.5737299999996,1038.447998,0.141123,0.141123,-0.272976,0.9673510000000001,0.0,0.0
+P03_002371_1888_XI_08N204W_78,1,0,1,4083.8732909999994,1277.689209,0.141123,0.141123,-0.328486,0.8362860000000001,0.0,0.0
+P01_001606_1897_XI_09N203W_19,1,0,1,-3917.203125,-1175.195679,0.141123,0.141123,-0.022262999999999998,-0.15115599999999998,0.0,0.0
+P01_001606_1897_XI_09N203W_65,1,0,1,4070.48584,-2315.9921879999997,0.141123,0.141123,-0.12137,-0.47907,0.0,0.0
+P03_002226_1895_XI_09N203W_1,1,0,0,-1197.647949,-569.925598,0.0,0.0,-0.070627,0.194853,0.0,0.0
+P03_002226_1895_XI_09N203W_2,1,0,0,-1241.685791,-499.006836,0.0,0.0,-0.066482,0.270447,0.0,0.0
diff --git a/plio/io/io_bae.py b/plio/io/io_bae.py
index e0c3d4fb9e87d84c6c28bfa9894a35c88ab9c81d..9f1dbcc99b4399f805d7f58bed893471e752d250 100644
--- a/plio/io/io_bae.py
+++ b/plio/io/io_bae.py
@@ -1,5 +1,6 @@
import json
import re
+from functools import singledispatch
import numpy as np
import pandas as pd
@@ -46,12 +47,83 @@ def socetset_keywords_to_json(keywords, ell=None):
# Case where the values are on a newline after the key
nums = numeric_matcher.findall(l)
stream[key] += map(float, nums)
-
+
parse(keywords)
if ell:
parse(ell)
return json.dumps(stream)
-
+
+@singledispatch
+def read_ipf(arg): # pragma: no cover
+ return str(arg)
+
+@read_ipf.register(str)
+def read_ipf_str(input_data):
+ """
+ Read a socet ipf file into a pandas data frame
+
+ Parameters
+ ----------
+ input_data : str
+ path to the an input data file
+
+ Returns
+ -------
+ df : pd.DataFrame
+ containing the ipf data with appropriate column names and indices
+ """
+
+ # Check that the number of rows is matching the expected number
+ with open(input_data, 'r') as f:
+ for i, l in enumerate(f):
+ if i == 1:
+ cnt = int(l)
+ elif i == 2:
+ col = l
+ break
+
+ columns = np.genfromtxt(input_data, skip_header=2, dtype='unicode',
+ max_rows = 1, delimiter = ',')
+
+ # TODO: Add unicode conversion
+ d = [line.split() for line in open(input_data, 'r')]
+ d = np.hstack(np.array(d[3:]))
+
+ d = d.reshape(-1, 12)
+
+ df = pd.DataFrame(d, columns=columns)
+
+ assert int(cnt) == len(df), 'Dataframe length {} does not match point length {}.'.format(int(cnt), len(df))
+
+ # Soft conversion of numeric types to numerics, allows str in first col for point_id
+ df = df.apply(pd.to_numeric, errors='ignore')
+
+ return df
+
+@read_ipf.register(list)
+def read_ipf_list(input_data_list):
+ """
+ Read a socet ipf file into a pandas data frame
+
+ Parameters
+ ----------
+ input_data_list : list
+ list of paths to the a set of input data files
+
+ Returns
+ -------
+ df : pd.DataFrame
+ containing the ipf data with appropriate column names and indices
+ """
+ frames = []
+
+ for input_file in input_data_list:
+ frames.append(read_ipf(input_file))
+
+ df = pd.concat(frames)
+
+ return df
+
def read_gpf(input_data):
"""
Read a socet gpf file into a pandas data frame
@@ -76,22 +148,35 @@ def read_gpf(input_data):
col = l
break
+ default_columns = np.genfromtxt(input_data, skip_header=2, dtype='unicode',
+ max_rows = 1, delimiter = ',')
+
+ columns = []
+
+ for column in default_columns:
+
+ if '(' in column and ')' in column:
+ column_name ,suffix = column.split('(')
+ num = int(suffix.split(')')[0])
+
+ for column_num in range(int(num)):
+ new_column = '{}{}'.format(column_name, column_num)
+ columns.append(new_column);
+
+ else:
+ columns.append(column)
+
# Mixed types requires read as unicode - let pandas soft convert
d = np.genfromtxt(input_data, skip_header=3, dtype='unicode')
d = d.reshape(-1, 12)
- #TODO: cols should be used to dynamically generate the column names
-
- df = pd.DataFrame(d, columns=['point_id', 'stat', 'known',
- 'lat_Y_North', 'long_X_East','ht',
- 'sigma0', 'sigma1', 'sigma2',
- 'res0', 'res1', 'res2'])
+ df = pd.DataFrame(d, columns=columns)
# Soft conversion of numeric types to numerics, allows str in first col for point_id
df = df.apply(pd.to_numeric, errors='ignore')
# Validate the read data with the header point count
- assert int(cnt) == len(df)
+ assert int(cnt) == len(df), 'Dataframe length {} does not match point length {}.'.format(int(cnt), len(df))
return df
@@ -117,7 +202,7 @@ def save_gpf(df, output_file):
try:
outGPF = open(output_file, 'w', newline='\r\n')
except:
- print ('Unable to open output gpf file: {0}'.format(output_file))
+ print('Unable to open output gpf file: {0}'.format(output_file))
return 1
#grab number of rows in pandas dataframe
@@ -132,8 +217,8 @@ def save_gpf(df, output_file):
#Output coordinates to gpf file
outGPF.write('{0} {1} {2}\n'.format(row['point_id'], row['stat'], row['known']))
outGPF.write('{0} {1} {2}\n'.format(row['lat_Y_North'], row['long_X_East'], row['ht']))
- outGPF.write('{0} {1} {2}\n'.format(row['sigma0'], row['sigma1'], row['sigma2']))
+ outGPF.write('{0} {1} {2}\n'.format(row['sig0'], row['sig1'], row['sig2']))
outGPF.write('{0} {1} {2}\n\n'.format(row['res0'], row['res1'], row['res2']))
outGPF.close()
- return
\ No newline at end of file
+ return
diff --git a/plio/io/tests/test_io_bae.py b/plio/io/tests/test_io_bae.py
index 330361160638e0c3459aad7b17c476b59ec9048d..105680bf6116b8923955777f113617fb87f498b3 100644
--- a/plio/io/tests/test_io_bae.py
+++ b/plio/io/tests/test_io_bae.py
@@ -1,9 +1,10 @@
import json
+
import numpy as np
import pandas as pd
from pandas.util.testing import assert_frame_equal
-from plio.io.io_bae import socetset_keywords_to_json, read_gpf, save_gpf
+from plio.io.io_bae import socetset_keywords_to_json, read_gpf, save_gpf, read_ipf
from plio.examples import get_path
import pytest
@@ -13,27 +14,40 @@ def insight_gpf():
return get_path('InSightE08_XW.gpf')
@pytest.fixture()
-def insight_expected():
+def insight_expected_gpf():
return pd.read_csv(get_path('InSightE08_XW.csv'))
-@pytest.mark.parametrize('gpf, expected', [(insight_gpf(),insight_expected())])
+@pytest.fixture
+def insight_ipf():
+ return get_path('P20_008845_1894_XN_09N203W.ipf')
+
+@pytest.fixture()
+def insight_expected_ipf():
+ return pd.read_csv(get_path('P20_008845_1894_XN_09N203W.csv'))
+
+@pytest.mark.parametrize('gpf, expected', [(insight_gpf(),insight_expected_gpf())])
def test_read_gfp(gpf, expected):
df = read_gpf(gpf)
assert_frame_equal(df, expected)
-@pytest.mark.parametrize('gpf', [(insight_gpf())])
-def test_write_gpf(gpf):
+@pytest.mark.parametrize('ipf, expected', [([insight_ipf()],insight_expected_ipf())])
+def test_read_ifp(ipf, expected):
+ df = read_ipf(ipf)
+ assert_frame_equal(df, expected)
+
+@pytest.mark.parametrize('gpf, file', [(insight_gpf(), 'out.gpf')])
+def test_write_gpf(gpf, file):
"""
We test by manually comparing files and not using filecmp so that we
are not testing float point precision differences, e.g. 0.0 == 0.00000000.
"""
df = read_gpf(gpf)
- save_gpf(df, 'out.gpf')
+ save_gpf(df, file)
with open(gpf) as f:
fl = f.readlines()
- with open('out.gpf') as f:
+ with open(file) as f:
fs = f.readlines()
# Check that the header is the same
@@ -41,7 +55,7 @@ def test_write_gpf(gpf):
assert fl[i] == fs[i]
truth_arr = np.genfromtxt(gpf, skip_header=3)
- test_arr = np.genfromtxt('out.gpf', skip_header=3)
+ test_arr = np.genfromtxt(file, skip_header=3)
np.testing.assert_array_almost_equal(truth_arr, test_arr)
def test_create_from_socet_lis():
@@ -49,8 +63,8 @@ def test_create_from_socet_lis():
socetell = get_path('ellipsoid.ell')
js = json.loads(socetset_keywords_to_json(socetlis))
assert isinstance(js, dict) # This is essentially a JSON linter
- # Manually validated
+ # Manually validated
assert 'RECTIFICATION_TERMS' in js.keys()
assert 'SEMI_MAJOR_AXIS' in js.keys() # From ellipsoid file
assert 'NUMBER_OF_EPHEM' in js.keys()
- assert len(js['EPHEM_PTS']) / 3 == js['NUMBER_OF_EPHEM']
\ No newline at end of file
+ assert len(js['EPHEM_PTS']) / 3 == js['NUMBER_OF_EPHEM']