Skip to content
Snippets Groups Projects
Commit e1a28769 authored by Mattia Mancini's avatar Mattia Mancini
Browse files

SSB-42: setting up directory's tree and the testing infrastructure.

parent 2cb4f3e5
No related branches found
No related tags found
1 merge request!44Merge back holography to master
Showing
with 1663 additions and 4 deletions
......@@ -7,13 +7,131 @@ CAL/CalibrationCommon/doc/CalibrationCommon.md -text
CAL/CalibrationCommon/doc/Holography_Data_Set.txt -text
CAL/CalibrationCommon/doc/package.dox -text
CAL/CalibrationCommon/lib/CMakeLists.txt -text
CAL/CalibrationCommon/lib/__init__.py -text
CAL/CalibrationCommon/lib/classes.py -text
CAL/CalibrationCommon/lib/coordinates.py -text
CAL/CalibrationCommon/lib/mshologextract.py -text
CAL/CalibrationCommon/lib/utils.py -text
CAL/CalibrationCommon/test/CMakeLists.txt -text
CAL/CalibrationCommon/test/t_holography_beam_specification_class.py -text
CAL/CalibrationCommon/test/t_holography_beam_specification_class.run -text
CAL/CalibrationCommon/test/t_holography_beam_specification_class.sh -text
CAL/CalibrationCommon/test/t_holography_dataset_class.py -text
CAL/CalibrationCommon/test/t_holography_dataset_class.run -text
CAL/CalibrationCommon/test/t_holography_dataset_class.sh -text
CAL/CalibrationCommon/test/t_holography_ms_class.py -text
CAL/CalibrationCommon/test/t_holography_ms_class.run -text
CAL/CalibrationCommon/test/t_holography_ms_class.sh -text
CAL/CalibrationCommon/var/CMakeLists.txt -text
CAL/CalibrationCommon/var/antenna-fields/CS001-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS001-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS002-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS002-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS003-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS003-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS004-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS004-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS005-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS005-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS006-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS006-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS007-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS007-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS011-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS011-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS013-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS013-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS017-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS017-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS021-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS021-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS024-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS024-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS026-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS026-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS028-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS028-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS030-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS030-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS031-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS031-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS032-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS032-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS101-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS101-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS103-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS103-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS201-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS201-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS301-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS301-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS302-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS302-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS401-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS401-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS501-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/CS501-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/DE601-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/DE601-iHBADeltas.JMA.conf -text
CAL/CalibrationCommon/var/antenna-fields/DE601-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/DE602-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/DE602-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/DE603-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/DE603-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/DE604-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/DE604-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/DE605-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/DE605-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/DE609-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/DE609-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/FI609-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/FI609-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/FR606-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/FR606-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/IE613-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/IE613-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/PL610-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/PL610-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/PL611-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/PL611-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/PL612-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/PL612-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS106-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS106-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS205-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS205-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS208-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS208-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS210-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS210-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS305-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS305-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS306-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS306-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS307-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS307-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS310-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS310-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS406-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS406-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS407-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS407-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS409-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS409-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS503-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS503-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS508-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS508-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS509-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS509-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS511-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/RS511-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/SE607-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/SE607-CableDelays.conf -text
CAL/CalibrationCommon/var/antenna-fields/SE607-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/UK608-AntennaField.conf -text
CAL/CalibrationCommon/var/antenna-fields/UK608-iHBADeltas.conf -text
CAL/CalibrationCommon/var/antenna-fields/iHBADeltas_default.conf -text
CEP/Calibration/BBSControl/include/BBSControl/CommandHandlerEstimator.h -text
CEP/Calibration/BBSControl/include/BBSControl/CommandHandlerReducer.h -text
CEP/Calibration/BBSControl/include/BBSControl/OptionParser.h -text
......
......@@ -17,10 +17,11 @@
# $Id$
lofar_package(CalibrationCommon 1.0)
lofar_package(CalibrationCommon 1.0 DEPENDS PyCommon)
include(PythonInstall)
add_subdirectory(var)
add_subdirectory(lib)
add_subdirectory(bin)
add_subdirectory(test)
......@@ -19,9 +19,8 @@
python_install(
__init__.py
holography_ms.py
holography_beam_specification.py
classes.py
coordinates.py
holography.py
utils.py
DESTINATION lofar/calibration/common)
import re
import os
from collections import defaultdict
import datetime
from casacore.tables import table as MS_Table
class HolographySpecification(object):
hs_name_pattern = r'Holog-(?P<date>\d{8})-(?P<comment>.*)-(?P<id>\d{3}).txt'
def __init__(self, name, path):
self.path = os.path.join(path, name)
self.name = name
self.id, self.date, self.comment = HolographySpecification.\
extract_id_date_comment_from_name(name)
self.station_specification_map = defaultdict(list)
self.start_datetime = None
self.end_datetime = None
self.rcu_mode = None
self.beam_set_interval = None
def __repr__(self):
return 'HolographySpecification(%s, %s, %s, %s, %s)' % (
self.id,
self.date,
self.comment,
self.name,
self.path,
)
def _read_lines(self):
with open(self.path, 'r') as fstream_in:
return fstream_in.read().splitlines()
@staticmethod
def _split_header(line):
date_regex = '\d{4}-\d{2}-\d{2}\s\d{2}:\d{2}:\d{2}'
format = r'(?P<start_date>{date_regex})\s*' \
'(?P<end_date>{date_regex})\s*' \
'(?P<rcu_mode>\d*)\s*' \
'(?P<beam_switch_delay>\d*.\d*)'.format(date_regex=date_regex)
match = re.match(format, line)
return match.groupdict()
def _parse_header(self, header):
splitted_header = HolographySpecification._split_header(header)
self.start_datetime = splitted_header['start_date']
self.end_datetime = splitted_header['end_date']
self.rcu_mode = splitted_header['rcu_mode']
self.beam_set_interval = splitted_header['beam_switch_delay']
@staticmethod
def _split_line(line):
range_regex = '(\d*\:\d*)|(\d*)'
ra_dec_regex = '\d*\.\d*,-?\d*\.\d*,\w*'
regex = r'^(?P<station_name>\w*)\s*' \
r'(?P<mode_description>\w*)\s*' \
r'(?P<sub_band>[\d,]*)\s*' \
r'(?P<beamlets>{range_regex})\s*' \
r'(?P<rcus>{range_regex})\s*' \
r'(?P<rcus_mode>(\d*))\s*' \
r'(?P<virtual_pointing>{ra_dec_regex})\s*' \
r'(?P<station_pointing>{ra_dec_regex})'.format(range_regex=range_regex,
ra_dec_regex=ra_dec_regex)
match = re.match(regex, line)
if match is None:
raise ValueError('Cannot parse line {}'.format(line))
return match.groupdict()
@staticmethod
def _split_lines(lines):
return [HolographySpecification._split_line(line)
for line in lines]
@staticmethod
def _parse_pointing(pointing_string):
ra, dec, coordinate_system = pointing_string.split(',')
ra = float(ra)
dec = float(dec)
return dict(ra=ra, dec=dec, coordinate_system=coordinate_system)
@staticmethod
def _parse_line(splitted_line):
rcus_mode = int(splitted_line['rcus_mode'])
sub_band_ids = [int(sub_band) for sub_band in splitted_line['sub_band'].split(',')]
mode_description = splitted_line['mode_description']
rcus_involved = splitted_line['rcus']
beamlets = splitted_line['beamlets']
station_pointing = HolographySpecification._parse_pointing(
splitted_line['station_pointing'])
virtual_pointing = HolographySpecification._parse_pointing(
splitted_line['virtual_pointing'])
if len(sub_band_ids) == 1:
station_type = 'target'
else:
station_type = 'reference'
station_name = splitted_line['station_name']
return dict(station_name=station_name,
rcus_mode=rcus_mode,
sub_band_ids=sub_band_ids,
mode_description=mode_description,
rcus_involved=rcus_involved,
beamlets=beamlets,
station_pointing=station_pointing,
virtual_pointing=virtual_pointing,
station_type=station_type)
def _parse_lines(self, lines):
splitted_lines = HolographySpecification._split_lines(lines)
for line in splitted_lines:
parsed_line = HolographySpecification._parse_line(line)
self.station_specification_map[parsed_line['station_name']] += [parsed_line]
def _update_class_attributes(self):
self.station_names = self.station_specification_map.keys()
self.reference_station_names = [station_name for station_name in
self.station_specification_map
if len(self.station_specification_map[station_name]) == 1]
self.target_station_names = [station_name for station_name in
self.station_specification_map
if len(self.station_specification_map[station_name]) > 1]
print(self.target_station_names, self.reference_station_names)
def read_file(self):
lines = self._read_lines()
self._parse_header(lines[0])
self._parse_lines(lines[1:])
self._update_class_attributes()
@staticmethod
def create_hs_list_from_name_list_and_path(name_list, path):
return [HolographySpecification(name, path) for name in name_list]
@staticmethod
def is_holography_specification_file_name(name):
return re.match(HolographySpecification.hs_name_pattern, name) is not None
@staticmethod
def extract_id_date_comment_from_name(name):
match = re.match(HolographySpecification.hs_name_pattern, name)
print(name, match)
date = match.group('date')
hs_id = int(match.group('id'))
comment = match.group('comment')
date = datetime.datetime.strptime(date, '%Y%m%d')
return hs_id, date, comment
class MeasurementSet(object):
ms_name_pattern = r'L(?P<sas_id>\d{6})_SB(?P<sub_band_id>\d{3})_uv\.MS'
def __init__(self, ms_name, ms_path):
self.path = os.path.join(ms_path, ms_name)
if not MeasurementSet.is_a_valid_ms_name(ms_name):
raise ValueError('The measurement set located in %s has not a valid name' % self.path,)
self.name = ms_name
self.sas_id, self.beamlet = \
MeasurementSet.parse_sas_id_and_sub_band_from_ms_name(self.name)
def get_data_table(self):
data_table = MS_Table(self.path)
return data_table
def get_antenna_table(self):
antenna_table = MS_Table(self.path + '/ANTENNA')
return antenna_table
def read_cross_correlation_per_station_names(self, reference, target):
data_table = self.get_data_table()
antennas_table = self.get_antenna_table()
antenna_name_id_map = {name:i for i, name in enumerate(antennas_table.getcol('NAME'))}
antenna1_list = data_table.getcol('ANTENNA1')
antenna2_list = data_table.getcol('ANTENNA2')
timestamp = data_table.getcol('TIME')
cross_correlation = data_table.getcol('DATA')
reference_antenna_id = antenna_name_id_map[reference]
target_antenna_id = antenna_name_id_map[target]
selected_data = [index for index, (a_i, a_j) in enumerate(zip(antenna1_list, antenna2_list))
if a_i == reference_antenna_id and a_j == target_antenna_id]
return (timestamp[selected_data], cross_correlation[selected_data])
def __repr__(self):
return 'MeasurementSet(%d) located in %s for sas_id %d and sub_band_id %d' % (id(self),
self.name,
self.sas_id,
self.beamlet)
@staticmethod
def create_ms_dict_from_ms_name_list_and_path(list_of_ms_names, path):
filtered_list_of_ms_names = MeasurementSet.filter_valid_ms_names(list_of_ms_names)
ms_list = [MeasurementSet(ms_name, path) for ms_name in filtered_list_of_ms_names]
return {ms.beamlet:ms for ms in ms_list}
@staticmethod
def parse_sas_id_and_sub_band_from_ms_name(ms_name):
if not MeasurementSet.is_a_valid_ms_name(ms_name):
raise ValueError('The measurement set %s has not a valid name' % ms_name,)
match = re.match(MeasurementSet.ms_name_pattern, ms_name)
return (int(match.group('sas_id')), int(match.group('sub_band_id')))
@staticmethod
def is_a_valid_ms_name(ms_name):
pattern = MeasurementSet.ms_name_pattern
return re.match(pattern, ms_name.strip()) # is not None
@staticmethod
def filter_valid_ms_names(list_of_ms_names):
return list(filter(MeasurementSet.is_a_valid_ms_name, list_of_ms_names))
class HolographyDataset():
def __init__(self):
self.rcu_list = () # array of ints
self.mode = None # integer
self.sas_ids = () # array of strings
self.target_station_name = None # string
self.target_station_position = None # list of 3 floats
self.source_name = None # string
self.source_position = None # list of 3 floats
self.start_time = None # date time when the observation started in MJD (float)
self.end_time = None # date time when the observation started in MJD (float)
self.rotation_matrix = None # array(3,3), translation matrix for
# (RA, DEC) <-> (l, m) conversion
self.antenna_field_position = () # coordinates of the antenna position in the target
# station
self.reference_stations = () # list of reference station names
self.frequencies = () # list of frequencies
self.ra_dec = () # array(Nfrequency, Nbeamlets, 2) contains the ra_dec of which a beam
# points at given a frequency and a beamlet number
self.data = () # array(NreferenceStations, Nfrequencies, Nbeamlets) that contains the
# 4 polarization crosscorrelation for the 4 polarizations, the l and m coordinates, and
# the timestamp in mjd of the sample
def load_from_specification_and_ms(self, station_name, hb_specifications, h_measurement_set):
"""
Loads the dataset from the specification files and the measurements for the given station
name
:param station_name: target station name
:param hb_specifications: list of holography beam specification files
:param h_measurement_set: map of the measurement set indexed by beamlet number
"""
@staticmethod
def load_from_file(path):
"""
It reads the dataset from file and returns a HolographyDataset class
:param path: path to file
:return: the read dataset
"""
result = HolographyDataset()
raise NotImplementedError
def store_to_file(self, path):
"""
Stores the holography dataset at the given path
:param path: path to file
"""
raise NotImplementedError
\ No newline at end of file
from numpy import cos, sin, arcsin, arctan2, sqrt, array, dot, float64
from astropy.coordinates import Angle
from astropy.coordinates import SkyCoord
import astropy.units as u
# Always try to download the most recent IERS tables.
from astropy.utils.data import download_file
from astropy.utils import iers
iers.IERS.iers_table = iers.IERS_A.open(download_file(iers.IERS_A_URL, cache=True))
def lm_from_radec(ra, dec, ra0, dec0):
r'''
Calculate direction cosines given a right ascension and
declination and the right ascension and declination of the
projection centre.
**Parameters**
ra : astropy.coordinates.Angle
Right ascension to convert.
dec : astropy.coordinates.Angle
Declination to convert.
ra0 : astropy.coordinates.Angle
Right ascension of the projection centre.
dec0 : astropy.coordinates.Angle
Declination of the projection centre.
**Returns**
A tuple of floats (l_rad, m_rad).
**Examples**
>>> import astropy.units as u
>>> lm_from_radec(ra=Angle(2.0, u.rad), dec=Angle(1.1, u.rad),
... ra0=Angle(2.0, u.rad), dec0=Angle(1.0, u.rad))
(0.0, 0.099833416646828321)
>>> lm_from_radec(ra=Angle(2.2, u.rad), dec=Angle(0.0, u.rad),
... ra0=Angle(2.0, u.rad), dec0=Angle(0.0, u.rad))
(0.19866933079506138, 0.0)
>>> lm_from_radec(ra=Angle(2.0, u.rad), dec=Angle(1.0, u.rad),
... ra0=Angle(2.0, u.rad), dec0=Angle(1.0, u.rad))
(0.0, 0.0)
>>> lm_from_radec(ra=Angle(1.8, u.rad), dec=Angle(0.9, u.rad),
... ra0=Angle(2.0, u.rad), dec0=Angle(1.0, u.rad))
(-0.1234948364118721, -0.089406906258670149)
'''
cos_dec = cos(dec.rad)
sin_dec = sin(dec.rad)
cos_dec0 = cos(dec0.rad)
sin_dec0 = sin(dec0.rad)
sin_dra = sin(float(ra.rad - ra0.rad))
cos_dra = cos(float(ra.rad - ra0.rad))
l_rad = cos_dec*sin_dra
m_rad = sin_dec*cos_dec0 - cos_dec*sin_dec0*cos_dra
return (l_rad, m_rad)
def radec_from_lm(l_rad, m_rad, ra0, dec0):
r'''
Calculate right ascension and declination given direction cosines
l and m and the RA and Dec of the projection centres.
**Parameters**
l_rad : float
Direction cosine parallel to right ascension at the projection
centre (in rad).
m_rad : float
Direction cosine parallel to declination at the projection
centre (in rad).
ra0 : astropy.coordinates.Angle
Right ascension of the projection centre.
dec0 : astropy.coordinates.Angle
Declination of the projection centre.
**Returns**
A tuple of astropy.coordinates.Angles (ra, dec)
**Examples**
>>> import astropy.units as u
>>> radec_from_lm(l_rad=0.0, m_rad=0.099833416646828321,
... ra0=Angle(2.0, u.rad), dec0=Angle(1.0, u.rad))
(<Angle 2.0 rad>, <Angle 1.1 rad>)
>>> radec_from_lm(l_rad=0.19866933079506138, m_rad=0.0,
... ra0=Angle(2.0, u.rad), dec0=Angle(0.0, u.rad))
(<Angle 2.2 rad>, <Angle 0.0 rad>)
>>> radec_from_lm(l_rad=-0.1234948364118721, m_rad=-0.089406906258670149,
... ra0=Angle(2.0, u.rad), dec0=Angle(1.0, u.rad))
(<Angle 1.8 rad>, <Angle 0.9 rad>)
'''
n_rad = sqrt(1.0 - l_rad*l_rad - m_rad*m_rad)
cos_dec0 = cos(dec0.rad)
sin_dec0 = sin(dec0.rad)
ra_rad = ra0.rad + arctan2(l_rad, cos_dec0*n_rad - m_rad*sin_dec0)
dec_rad = arcsin(m_rad*cos_dec0 + sin_dec0*n_rad)
return (Angle(ra_rad, u.rad), Angle(dec_rad, u.rad))
def icrs_from_itrs(unit_vector_itrs, obstime):
r'''
Convert a geocentric cartesian unit vector in the ITRS system into
an astropy.coordinates.SkyCoord in the ICRS system (equinox
J2000), given an observing date/time. This routine is used to find
out the RA and DEC of the direction to which a station's normal
vector points.
**Parameters**
unit_vector_itrs : numpy array of 3 floats
The direction to convert it is assumed to be in the ITRS
system.
obstime : astropy.time.Time or string
When a string is provided, it is assumed to be readable by an
astropy.time.Time instance.
**Returns**
A SkyCoord containing the ICRS position.
**Examples**
>>> icrs_from_itrs([0.0, 0.0, 1.0], '2000-01-01 00:00:00')
<SkyCoord (ICRS): (ra, dec) in deg
(349.95496154, 89.99588844)>
>>> icrs_from_itrs([0.0, 1.0, 0.0], '2000-03-21 12:00:00')
<SkyCoord (ICRS): (ra, dec) in deg
(89.30969017, 0.00329331)>
>>> icrs_from_itrs([0.0, 0.0, 1.0], '2015-01-01 00:00:00')
<SkyCoord (ICRS): (ra, dec) in deg
(358.7928571, 89.91033405)>
'''
x, y, z = array(unit_vector_itrs).T
c_itrs = SkyCoord(x, y, z, representation='cartesian',frame='itrs',
obstime=obstime, equinox='J2000')
return c_itrs.icrs
def itrs_from_icrs(icrs_position_rad, obstime):
r'''
Takes an array RA/Dec ICRS positions in radians and converts those
to geocentric ITRF unit vectors.
**Parameters**
icrs_position_rad : numpy array of floats (ra_rad, dec_rad)
The ICRS position to convert. May also be an array of RA/DEC
pairs.
obstime : astropy.time.Time or string
When a string is provided, it is assumed to be readable by an
astropy.time.Time instance.
**Returns**
An array containing the geocentric cartesian ITRS unit vectors
corresponding to the icrs_position_rad at obstime.
**Examples**
>>> itrs_from_icrs((array([358.7928571, 89.91033405])*u.deg).to(u.rad),
... obstime='2015-01-01 00:00:00')
array([[ 1.25385143e-12, 6.93857340e-12, 1.00000000e+00]])
>>> itrs_from_icrs((array([[358.7928571, 89.91033405],
... [90,-20],
... [30, 60]])*u.deg).to(u.rad),
... obstime='2015-01-01 00:00:00')
array([[ 1.25385143e-12, 6.93857340e-12, 1.00000000e+00],
[ 9.24930853e-01, -1.65788673e-01, -3.42077526e-01],
[ 1.70157684e-01, -4.68937638e-01, 8.66685557e-01]])
'''
ra, dec = array(icrs_position_rad, dtype='float64').T*u.rad
icrs = SkyCoord(ra, dec,frame='icrs',
obstime=obstime, equinox='J2000')
itrs = icrs.itrs
return array([itrs.x, itrs.y, itrs.z], dtype=float64).T
def pqr_from_icrs(icrs_rad, obstime, pqr_to_itrs_matrix):
r'''
Compute geocentric station-local PQR coordinates of a certain ICRS
direction. Geocentric here means that parallax between the centre
of the Earth and the station's phase centre is not taken into
account.
**Parameters**
icrs_rad : numpy.array
An ICRS position in radians. Either one vector of length 2, or
an array of shape (N, 2) containing N ICRS [RA, Dec] pairs.
obstime : string or astropy.time.Time
The date/time of observation.
pqr_to_itrs_matrix : 3x3 numpy.array
The rotation matrix that is used to convert a direction in the
station's PQR system into an ITRS direction. This matrix is
found in the /opt/lofar/etc/AntennaField.conf files at the
stations. These files are also found in the antenna-fields/
directory of this project.
**Returns**
A numpy.array instance with shape (3,) or (N, 3) containing the
ITRS directions.
**Example**
>>> import astropy
>>> core_pqr_itrs_mat = array([[ -1.19595000e-01, -7.91954000e-01, 5.98753000e-01],
... [ 9.92823000e-01, -9.54190000e-02, 7.20990000e-02],
... [ 3.30000000e-05, 6.03078000e-01, 7.97682000e-01]],
... dtype=float64)
>>> obstime='2015-06-19 13:50:00'
>>> target_3c196 = SkyCoord('08h13m36.0561s', '+48d13m02.636s', frame='icrs')
>>> target_icrs_rad = array([target_3c196.icrs.ra.rad, target_3c196.icrs.dec.rad])
>>> pqr_from_icrs(target_icrs_rad, obstime, core_pqr_itrs_mat)
array([ 0.02131259, -0.08235505, 0.99637513])
>>> pqr_from_icrs(array([target_icrs_rad, target_icrs_rad]),
... obstime, core_pqr_itrs_mat)
array([[ 0.02131259, -0.08235505, 0.99637513],
[ 0.02131259, -0.08235505, 0.99637513]])
>>> pqr_from_icrs(target_icrs_rad, astropy.time.Time(obstime)+7*u.minute +18*u.second,
... core_pqr_itrs_mat)
array([ 1.49721138e-05, -8.26257777e-02, 9.96580636e-01])
'''
return dot(pqr_to_itrs_matrix.T, itrs_from_icrs(icrs_rad, obstime).T).T.squeeze()
def icrs_from_pqr(pqr, obstime, pqr_to_itrs_matrix):
r'''
Convert directions from the station-local PQR system into
geocentric ICRS position.
**Parameters**
pqr : numpy array
The direction to convert it is assumed to be in the ITRS
system. Either a numpy.array of 3 floats, or an (N,3)
numpy.array contining a collection of PQR positions to
convert.
obstime : astropy.time.Time or string
When a string is provided, it is assumed to be readable by an
astropy.time.Time instance.
pqr_to_itrs_matrix : 3x3 numpy.array
The rotation matrix that is used to convert a direction in the
station's PQR system into an ITRS direction. This matrix is
found in the /opt/lofar/etc/AntennaField.conf files at the
stations. These files are also found in the antenna-fields/
directory of this project.
**Returns**
A SkyCoord containing the corresponding ICRS position(s).
**Examples**
>>> import astropy
>>> core_pqr_itrs_mat = array([[ -1.19595000e-01, -7.91954000e-01, 5.98753000e-01],
... [ 9.92823000e-01, -9.54190000e-02, 7.20990000e-02],
... [ 3.30000000e-05, 6.03078000e-01, 7.97682000e-01]],
... dtype=float64)
>>> obstime='2015-06-19 13:50:00'
>>> icrs_from_pqr(array([[0, -1, 0], [0, 0, 1], [0, 1, 0]]),
... obstime, core_pqr_itrs_mat)
<SkyCoord (ICRS): (ra, dec) in deg
[(121.70435893, -37.0409395), (121.55293128, 52.95456766),
(301.69271378, 37.04486433)]>
>>> icrs_from_pqr(array([[0, -1, 0], [0, 0, 1], [0, 1, 0]]),
... astropy.time.Time(obstime)+4*u.minute, core_pqr_itrs_mat)
<SkyCoord (ICRS): (ra, dec) in deg
[(122.70638666, -37.03964013), (122.55669495, 52.95575629),
(302.6948803, 37.04366179)]>
>>> source_3c196 = SkyCoord('08h13m36.0561s', '+48d13m02.636s', frame='icrs')
>>> pqr = pqr_from_icrs([source_3c196.icrs.ra.rad, source_3c196.icrs.dec.rad],
... obstime, core_pqr_itrs_mat)
>>> pqr
array([ 0.02131259, -0.08235505, 0.99637513])
>>> icrs_from_pqr(pqr, obstime, core_pqr_itrs_mat)
<SkyCoord (ICRS): (ra, dec) in deg
(123.40023895, 48.21740382)>
Let's check if that is indeed 3C 196:
>>> source_3c196
<SkyCoord (ICRS): (ra, dec) in deg
(123.40023375, 48.21739889)>
And it is not. The reason for this is the limited precision with
which the rotation matrices are listed in the AntennaField.conf
files. Sure enough:
>>> dot(core_pqr_itrs_mat, core_pqr_itrs_mat.T)
array([[ 9.99999257e-01, 2.84588000e-07, 5.09499000e-07],
[ 2.84588000e-07, 1.00000056e+00, -2.62005000e-07],
[ 5.09499000e-07, -2.62005000e-07, 9.99999648e-01]])
A 2e-7 differences correspond to about 1e-5 degrees, which isthe
level at which we see the difference.
'''
return icrs_from_itrs(dot(pqr_to_itrs_matrix, pqr.T).T, obstime=obstime)
from classes import *
from utils import *
import numpy
import argparse
from casacore.tables import table as MS_Table
def main():
cla_parser = specify_command_line_arguments()
arguments = parse_command_line_arguments(cla_parser)
bsfs, mss = read_holography_datasets(arguments.input_path)
target_station = bsfs[0].target_station_names[0]
reference_station = bsfs[0].reference_station_names[0]
write_to_numpy_array(bsfs[0].station_specification_map[target_station],
target_station,
reference_station,
mss.items()[0][1], 'dummy')
write_to_numpy_array(bsfs[0].station_specification_map[target_station],
target_station,
reference_station,
mss.items()[0][1], 'dummy2')
def parse_command_line_arguments(parser):
return parser.parse_args()
def specify_command_line_arguments():
parser = argparse.ArgumentParser(description='This program is meant for convert the Holography'
' observation\'s data into an holography dataset')
parser.add_argument('input_path', help='path to the holography observation data')
parser.add_argument('--holography_bsf', help='override default path for the holography beam'
' specification file')
parser.add_argument('--holography_ms', help='override default path for the holography'
' observation MS files')
return parser
def read_holography_datasets(holography_observation_path,
holography_bsf_path=None,
holography_ms_path=None):
if holography_ms_path is not None:
raise NotImplementedError()
if holography_bsf_path is not None:
raise NotImplementedError()
list_of_bsf_files = list_bsf_files_in_path(holography_observation_path)
list_of_ms = list_ms_files_in_path(holography_observation_path)
for bsf in list_of_bsf_files:
bsf.read_file()
return list_of_bsf_files, list_of_ms
def write_to_numpy_array(bsf_specifications, target, reference, ms_per_sub_band, location):
ra_dec_list, crosscorrelation_list,\
timestamp_list = extract_crosscorrelation_per_ra_dec_timestamp(bsf_specifications, target,
reference,
ms_per_sub_band)
ra_dec, timestamp, matrix = create_crosscorrelation_matrix(ra_dec_list,
crosscorrelation_list,
timestamp_list)
numpy.savez(location, ra_dec=numpy.array(list(ra_dec)), timestamp=timestamp, crosscorrelation=matrix)
def extract_crosscorrelation_per_ra_dec_timestamp(bsf_specifications,
target,
reference,
ms_per_sub_band):
ra_dec_list = []
crosscorrelation_list = []
timestamp_list = []
for bsf_specification in bsf_specifications:
beamlet = int(bsf_specification['beamlets'])
ms = ms_per_sub_band[beamlet]
timestamp, crosscorrelation = ms.read_cross_correlation_per_station_names(target,
reference)
virtual_pointing = bsf_specification['virtual_pointing']
ra = virtual_pointing['ra']
dec = virtual_pointing['dec']
ra_dec_list += [(ra, dec)]
crosscorrelation_list += list(crosscorrelation)
timestamp_list += list(timestamp)
return ra_dec_list, numpy.array(crosscorrelation_list), timestamp_list
def create_crosscorrelation_matrix(ra_dec_list, crosscorrelation_list, timestamp_list):
unique_timestamp = numpy.unique(timestamp_list)
unique_ra_dec = set(ra_dec_list)
polarizations = crosscorrelation_list[0].shape[1]
print('polarizations', polarizations)
number_of_pointings = len(ra_dec_list)
number_of_timestamps = len(unique_timestamp)
shape = (number_of_pointings, number_of_timestamps, polarizations)
print('shape', shape)
cross_correlation_matrix = numpy.zeros(shape, dtype=numpy.complex64)
for time_index, timestamp in enumerate(unique_timestamp):
selected_indexes = list(numpy.where(timestamp_list == timestamp)[0])
print(selected_indexes)
selected_crosscorrelations = crosscorrelation_list[selected_indexes]
for ra_dec_index, cross_correlation in enumerate(selected_crosscorrelations):
cross_correlation_matrix[ra_dec_index, time_index, :] = cross_correlation
return unique_ra_dec, unique_timestamp, cross_correlation_matrix
if __name__ == '__main__':
main()
\ No newline at end of file
from classes import *
import os
from glob import glob
def list_ms_files_in_path(path):
ms_dir_name_pattern = 'L(?P<sas_id>\d{6})'
ms_dirs_path_pattern = '^' + os.path.join(path, ms_dir_name_pattern, 'uv$')
matched_ms_paths_per_sas_id = {}
for root, dirnames, filenames in os.walk(path):
match = re.match(ms_dirs_path_pattern, root)
if match:
matched_ms_paths_per_sas_id[int(match.group('sas_id'))] =\
MeasurementSet.create_ms_dict_from_ms_name_list_and_path(dirnames, root)
return matched_ms_paths_per_sas_id
def list_bsf_files_in_path(path):
bsf_files_name_pattern = 'Holog-*.txt'
bsf_files_path_pattern = os.path.join(path, bsf_files_name_pattern)
matched_path_list = glob(bsf_files_path_pattern)
matched_file_path_list = list(filter(lambda path_i : os.path.isfile(path_i), matched_path_list))
matched_file_name_list = list(map(lambda path_i: os.path.basename(path_i),
matched_file_path_list))
holography_beam_specification = HolographySpecification.create_hs_list_from_name_list_and_path(
matched_file_name_list, path)
return holography_beam_specification
\ No newline at end of file
......@@ -21,5 +21,7 @@ include(LofarCTest)
lofar_add_test(t_holography_ms_class)
lofar_add_test(t_holography_beam_specification_class)
lofar_add_test(t_holography_dataset_class)
import logging
import unittest
from lofar.calibration.common.classes import HolographyDataset
logger = logging.getLogger('t_holography_dataset_class')
class TestHolographyDatasetClass(unittest.TestCase):
def test_one(self):
self.assertTrue(True)
\ No newline at end of file
#!/bin/bash
# Copyright (C) 2018 ASTRON (Netherlands Institute for Radio Astronomy)
# P.O. Box 2, 7990 AA Dwingeloo, The Netherlands
#
# This file is part of the LOFAR software suite.
# The LOFAR software suite is free software: you can redistribute it and/or
# modify it under the terms of the GNU General Public License as published
# by the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# The LOFAR software suite is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License along
# with the LOFAR software suite. If not, see <http://www.gnu.org/licenses/>.
# Run the unit test
source python-coverage.sh
python_coverage_test "*CAL*" t_holography_dataset_class.py
#!/bin/sh
# Copyright (C) 2018 ASTRON (Netherlands Institute for Radio Astronomy)
# P.O. Box 2, 7990 AA Dwingeloo, The Netherlands
#
# This file is part of the LOFAR software suite.
# The LOFAR software suite is free software: you can redistribute it and/or
# modify it under the terms of the GNU General Public License as published
# by the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# The LOFAR software suite is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License along
# with the LOFAR software suite. If not, see <http://www.gnu.org/licenses/>.
./runctest.sh t_holography_dataset_class
#!/usr/bin/env bash
# Copyright (C) 2018 ASTRON (Netherlands Institute for Radio Astronomy)
# P.O. Box 2, 7990 AA Dwingeloo, The Netherlands
#
# This file is part of the LOFAR software suite.
# The LOFAR software suite is free software: you can redistribute it and/or
# modify it under the terms of the GNU General Public License as published
# by the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# The LOFAR software suite is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License along
# with the LOFAR software suite. If not, see <http://www.gnu.org/licenses/>.
# $Id$
install(DIRECTORY
antenna-fields DESTINATION var/share/lofar_antenna_fields)
\ No newline at end of file
# Blitz-0.10 formatted
#
# AntennaPositions for CS001
# ITRF2005 target_date = 2015.5
# Created: 2015-01-20 21:24:05
#
NORMAL_VECTOR LBA
(0,2) [ 0.598753 0.072099 0.797682 ]
ROTATION_MATRIX LBA
(0,2) x (0,2) [
-0.1195950000 -0.7919540000 0.5987530000
0.9928230000 -0.0954190000 0.0720990000
0.0000330000 0.6030780000 0.7976820000
]
LBA
(0,2) [ 3826923.502750000 460915.488115000 5064643.517 ]
(0,95) x (0,1) x (0,2) [
0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
-2.019000 -0.243000 1.538000 -2.019000 -0.243000 1.538000
-1.338000 2.105000 0.814000 -1.338000 2.105000 0.814000
0.800000 2.363000 -0.814000 0.800000 2.363000 -0.814000
2.019000 0.243000 -1.538000 2.019000 0.243000 -1.538000
1.338000 -2.105000 -0.814000 1.338000 -2.105000 -0.814000
-0.800000 -2.363000 0.814000 -0.800000 -2.363000 0.814000
-0.210000 5.784000 -0.365000 -0.210000 5.784000 -0.365000
2.822000 4.398000 -2.516000 2.822000 4.398000 -2.516000
4.534000 0.954001 -3.489000 4.534000 0.954001 -3.489000
4.124000 -2.936999 -2.830000 4.124000 -2.936999 -2.830000
1.784000 -5.454000 -0.846000 1.784000 -5.454000 -0.846000
-1.390000 -5.418000 1.533000 -1.390000 -5.418000 1.533000
-3.914000 -2.848000 3.196000 -3.914000 -2.848000 3.196000
-4.606000 1.056000 3.362000 -4.606000 1.056000 3.362000
-3.143000 4.464000 1.956000 -3.143000 4.464000 1.956000
-4.555000 7.669000 2.726000 -4.555000 7.669000 2.726000
-2.165000 9.615000 0.756000 -2.165000 9.615000 0.756000
1.846000 9.900000 -2.281000 1.846000 9.900000 -2.281000
6.594000 5.379001 -5.436000 6.594000 5.379001 -5.436000
7.783000 2.445001 -6.063000 7.783000 2.445001 -6.063000
8.013000 -1.837999 -5.849000 8.013000 -1.837999 -5.849000
6.247000 -5.347999 -4.206000 6.247000 -5.347999 -4.206000
2.528000 -8.547000 -1.125000 2.528000 -8.547000 -1.125000
-1.137000 -9.695000 1.730000 -1.137000 -9.695000 1.730000
-4.464000 -8.074000 4.081000 -4.464000 -8.074000 4.081000
-6.224000 -4.693000 5.096000 -6.224000 -4.693000 5.096000
-7.236000 0.039000 5.428000 -7.236000 0.039000 5.428000
-6.583000 4.570000 4.528000 -6.583000 4.570000 4.528000
-6.093000 12.525000 3.442000 -6.093000 12.525000 3.442000
-3.475000 14.313000 1.315000 -3.475000 14.313000 1.315000
-0.972000 14.115000 -0.546000 -0.972000 14.115000 -0.546000
5.095000 11.170001 -4.834000 5.095000 11.170001 -4.834000
7.477000 11.350001 -6.638000 7.477000 11.350001 -6.638000
9.876000 4.350001 -7.806000 9.876000 4.350001 -7.806000
10.504000 -0.687999 -7.823000 10.504000 -0.687999 -7.823000
11.373000 -5.926999 -8.001000 11.373000 -5.926999 -8.001000
9.181000 -7.718999 -6.194000 9.181000 -7.718999 -6.194000
6.374000 -10.725999 -3.815000 6.374000 -10.725999 -3.815000
1.229000 -13.161000 0.267000 1.229000 -13.161000 0.267000
-0.957000 -13.579000 1.946000 -0.957000 -13.579000 1.946000
-4.932000 -11.373000 4.730000 -4.932000 -11.373000 4.730000
-8.215000 -10.491001 7.115000 -8.215000 -10.491001 7.115000
-11.072000 -2.214001 8.511000 -11.072000 -2.214001 8.511000
-11.116000 1.756999 8.185000 -11.116000 1.756999 8.185000
-10.397000 7.095999 7.163000 -10.397000 7.095999 7.163000
-7.780010 64.586999 0.002000 -7.780010 64.586999 0.002000
51.520000 6.207005 -39.233000 51.520000 6.207005 -39.233000
-11.566000 11.404999 7.650000 -11.566000 11.404999 7.650000
-11.932000 14.456999 7.650000 -11.932000 14.456999 7.650000
-0.701000 20.337000 -1.312000 -0.701000 20.337000 -1.312000
11.181000 13.262001 -9.591000 11.181000 13.262001 -9.591000
14.056000 9.704002 -11.428000 14.056000 9.704002 -11.428000
14.637000 -2.400998 -10.770000 14.637000 -2.400998 -10.770000
9.698000 -13.786999 -6.033000 9.698000 -13.786999 -6.033000
5.548000 -19.529999 -2.399000 5.548000 -19.529999 -2.399000
-5.804000 -16.092000 5.811000 -5.804000 -16.092000 5.811000
-11.813000 -12.457001 9.993000 -11.813000 -12.457001 9.993000
-14.654000 -3.007001 11.271000 -14.654000 -3.007001 11.271000
-19.883000 2.274998 14.719000 -19.883000 2.274998 14.719000
-16.889000 14.487999 11.367000 -16.889000 14.487999 11.367000
-9.526000 22.971999 5.074000 -9.526000 22.971999 5.074000
-4.898000 25.092000 1.408000 -4.898000 25.092000 1.408000
1.709000 26.419000 -3.671000 1.709000 26.419000 -3.671000
13.531000 16.011002 -11.604000 13.531000 16.011002 -11.604000
17.573000 9.587002 -14.057000 17.573000 9.587002 -14.057000
20.864000 -3.812998 -15.316000 20.864000 -3.812998 -15.316000
16.001000 -16.283998 -10.539000 16.001000 -16.283998 -10.539000
10.981000 -22.798999 -6.182000 10.981000 -22.798999 -6.182000
1.805000 -24.699000 0.878000 1.805000 -24.699000 0.878000
-11.054000 -22.397001 10.322000 -11.054000 -22.397001 10.322000
-20.215000 -8.270002 15.921000 -20.215000 -8.270002 15.921000
-27.744000 2.318998 20.616000 -27.744000 2.318998 20.616000
-22.547000 15.330998 15.538000 -22.547000 15.330998 15.538000
-15.013010 25.966999 8.922000 -15.013010 25.966999 8.922000
-10.063010 30.703999 4.779000 -10.063010 30.703999 4.779000
5.060000 30.482001 -6.553000 5.060000 30.482001 -6.553000
18.696000 23.593002 -16.166000 18.696000 23.593002 -16.166000
24.503000 13.316003 -19.596000 24.503000 13.316003 -19.596000
23.586000 4.116002 -18.076000 23.586000 4.116002 -18.076000
26.543000 -10.202997 -19.002000 26.543000 -10.202997 -19.002000
18.631000 -25.851998 -11.648000 18.631000 -25.851998 -11.648000
12.738000 -30.559999 -6.799000 12.738000 -30.559999 -6.799000
0.562000 -34.573000 2.703000 0.562000 -34.573000 2.703000
-5.545000 -31.083000 6.972000 -5.545000 -31.083000 6.972000
-20.102000 -20.428002 16.936000 -20.102000 -20.428002 16.936000
-21.281000 -13.823002 17.224000 -21.281000 -13.823002 17.224000
-26.161000 -15.007002 20.993000 -26.161000 -15.007002 20.993000
-28.304010 24.091998 19.068000 -28.304010 24.091998 19.068000
-11.638010 37.911999 5.309000 -11.638010 37.911999 5.309000
4.049990 41.232001 -6.767000 4.049990 41.232001 -6.767000
15.825000 33.938002 -14.946000 15.825000 33.938002 -14.946000
32.762000 -6.838997 -23.973000 32.762000 -6.838997 -23.973000
22.668000 -24.574998 -14.794000 22.668000 -24.574998 -14.794000
15.365000 -37.625998 -8.133000 15.365000 -37.625998 -8.133000
-18.877000 -29.400002 16.827000 -18.877000 -29.400002 16.827000
]
HBA
(0,2) [ 3826937.766750000 460938.573117000 5064630.724 ]
(0,47) x (0,1) x (0,2) [
-52.064010 42.569995 35.232000 -52.064010 42.569995 35.232000
-50.969010 47.440995 33.970000 -50.969010 47.440995 33.970000
-49.184010 36.068995 33.658000 -49.184010 36.068995 33.658000
-48.088010 40.939995 32.395000 -48.088010 40.939995 32.395000
-46.992010 45.809995 31.132000 -46.992010 45.809995 31.132000
-45.895010 50.680995 29.869000 -45.895010 50.680995 29.869000
-46.304010 29.565995 32.084000 -46.304010 29.565995 32.084000
-45.207010 34.437995 30.820000 -45.207010 34.437995 30.820000
-44.112010 39.308995 29.558000 -44.112010 39.308995 29.558000
-43.015010 44.179995 28.294000 -43.015010 44.179995 28.294000
-41.919010 49.049996 27.031000 -41.919010 49.049996 27.031000
-40.823010 53.920996 25.768000 -40.823010 53.920996 25.768000
-42.327010 27.935996 29.246000 -42.327010 27.935996 29.246000
-41.232010 32.806996 27.983000 -41.232010 32.806996 27.983000
-40.135010 37.677996 26.720000 -40.135010 37.677996 26.720000
-39.039010 42.548996 25.457000 -39.039010 42.548996 25.457000
-37.943010 47.419996 24.194000 -37.943010 47.419996 24.194000
-36.847010 52.290996 22.931000 -36.847010 52.290996 22.931000
-37.255010 31.175996 25.146000 -37.255010 31.175996 25.146000
-36.158010 36.047996 23.883000 -36.158010 36.047996 23.883000
-35.063010 40.917996 22.620000 -35.063010 40.917996 22.620000
-33.966010 45.788996 21.357000 -33.966010 45.788996 21.357000
-32.182010 34.416996 21.045000 -32.182010 34.416996 21.045000
-31.086010 39.287997 19.782000 -31.086010 39.287997 19.782000
31.084000 -39.287997 -19.782000 31.084000 -39.287997 -19.782000
32.180000 -34.416997 -21.044000 32.180000 -34.416997 -21.044000
33.965000 -45.788997 -21.356000 33.965000 -45.788997 -21.356000
35.061000 -40.917997 -22.619000 35.061000 -40.917997 -22.619000
36.157000 -36.047997 -23.882000 36.157000 -36.047997 -23.882000
37.253000 -31.176997 -25.145000 37.253000 -31.176997 -25.145000
36.845000 -52.290997 -22.930000 36.845000 -52.290997 -22.930000
37.941000 -47.419997 -24.194000 37.941000 -47.419997 -24.194000
39.037000 -42.548997 -25.456000 39.037000 -42.548997 -25.456000
40.133000 -37.677997 -26.720000 40.133000 -37.677997 -26.720000
41.230000 -32.807996 -27.983000 41.230000 -32.807996 -27.983000
42.326000 -27.936996 -29.246000 42.326000 -27.936996 -29.246000
40.821000 -53.920997 -25.768000 40.821000 -53.920997 -25.768000
41.917000 -49.049996 -27.031000 41.917000 -49.049996 -27.031000
43.014000 -44.179996 -28.294000 43.014000 -44.179996 -28.294000
44.110000 -39.308996 -29.557000 44.110000 -39.308996 -29.557000
45.206000 -34.437996 -30.820000 45.206000 -34.437996 -30.820000
46.302000 -29.566996 -32.083000 46.302000 -29.566996 -32.083000
45.894000 -50.680996 -29.868000 45.894000 -50.680996 -29.868000
46.990000 -45.809996 -31.131000 46.990000 -45.809996 -31.131000
48.086000 -40.939996 -32.394000 48.086000 -40.939996 -32.394000
49.182000 -36.068996 -33.657000 49.182000 -36.068996 -33.657000
50.967000 -47.440996 -33.969000 50.967000 -47.440996 -33.969000
52.063000 -42.569995 -35.232000 52.063000 -42.569995 -35.232000
]
NORMAL_VECTOR HBA0
(0,2) [ 0.598753 0.072099 0.797682 ]
ROTATION_MATRIX HBA0
(0,2) x (0,2) [
-0.1195950000 -0.7919540000 0.5987530000
0.9928230000 -0.0954190000 0.0720990000
0.0000330000 0.6030780000 0.7976820000
]
HBA0
(0,2) [ 3826896.191740000 460979.502113000 5064658.231 ]
NORMAL_VECTOR HBA1
(0,2) [ 0.598753 0.072099 0.797682 ]
ROTATION_MATRIX HBA1
(0,2) x (0,2) [
-0.1195950000 -0.7919540000 0.5987530000
0.9928230000 -0.0954190000 0.0720990000
0.0000330000 0.6030780000 0.7976820000
]
HBA1
(0,2) [ 3826979.340750000 460897.644121000 5064603.217 ]
#
# HBADeltas for CS001
# Created: 2015-01-20 21:24:05
#
HBADeltas
(0,31) x (0,2) [
-1.847 -1.180 1.493
-1.581 0.003 1.186
-1.315 1.185 0.880
-1.049 2.367 0.573
-0.882 -1.575 0.804
-0.616 -0.393 0.498
-0.350 0.789 0.191
-0.083 1.971 -0.116
0.083 -1.971 0.116
0.350 -0.789 -0.191
0.616 0.393 -0.498
0.882 1.575 -0.804
1.049 -2.367 -0.573
1.315 -1.185 -0.880
1.581 -0.003 -1.186
1.847 1.180 -1.493
-1.847 -1.180 1.493
-1.581 0.003 1.186
-1.315 1.185 0.880
-1.049 2.367 0.573
-0.882 -1.575 0.804
-0.616 -0.393 0.498
-0.350 0.789 0.191
-0.083 1.971 -0.116
0.083 -1.971 0.116
0.350 -0.789 -0.191
0.616 0.393 -0.498
0.882 1.575 -0.804
1.049 -2.367 -0.573
1.315 -1.185 -0.880
1.581 -0.003 -1.186
1.847 1.180 -1.493
]
# Blitz-0.10 formatted
#
# AntennaPositions for CS002
# ITRF2005 target_date = 2015.5
# Created: 2015-01-20 21:24:05
#
NORMAL_VECTOR LBA
(0,2) [ 0.598753 0.072099 0.797682 ]
ROTATION_MATRIX LBA
(0,2) x (0,2) [
-0.1195950000 -0.7919540000 0.5987530000
0.9928230000 -0.0954190000 0.0720990000
0.0000330000 0.6030780000 0.7976820000
]
LBA
(0,2) [ 3826577.022720000 461022.995082000 5064892.814 ]
(0,95) x (0,1) x (0,2) [
0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
-2.019000 -0.243000 1.538000 -2.019000 -0.243000 1.538000
-1.338000 2.105000 0.814000 -1.338000 2.105000 0.814000
0.800000 2.363000 -0.814000 0.800000 2.363000 -0.814000
2.019000 0.243000 -1.538000 2.019000 0.243000 -1.538000
1.338000 -2.105000 -0.814000 1.338000 -2.105000 -0.814000
-0.800000 -2.363000 0.814000 -0.800000 -2.363000 0.814000
-3.073000 -4.125000 2.679000 -3.073000 -4.125000 2.679000
-4.593000 -0.553000 3.497990 -4.593000 -0.553000 3.497990
-3.965000 3.277000 2.680000 -3.965000 3.277000 2.680000
-1.481000 5.575000 0.607000 -1.481000 5.575000 0.607000
1.696000 5.264000 -1.749000 1.696000 5.264000 -1.749000
4.079000 2.490001 -3.287000 4.079000 2.490001 -3.287000
4.553000 -1.449999 -3.287000 4.553000 -1.449999 -3.287000
2.897000 -4.710000 -1.749000 2.897000 -4.710000 -1.749000
-0.114000 -5.767000 0.607000 -0.114000 -5.767000 0.607000
-0.837000 -9.174000 1.457000 -0.837000 -9.174000 1.457000
-3.733000 -8.431000 3.564000 -3.733000 -8.431000 3.564000
-6.928000 -5.151000 5.665990 -6.928000 -5.151000 5.665990
-7.997000 2.044999 5.817990 -7.997000 2.044999 5.817990
-7.257000 5.076000 4.987990 -7.257000 5.076000 4.987990
-5.039000 8.195000 3.041990 -5.039000 8.195000 3.041990
-1.741000 9.059000 0.488000 -1.741000 9.059000 0.488000
2.862000 8.017000 -2.873000 2.862000 8.017000 -2.873000
6.276000 5.624001 -5.219000 6.276000 5.624001 -5.219000
7.890000 1.637001 -6.070000 7.890000 1.637001 -6.070000
7.333000 -2.192999 -5.306000 7.333000 -2.192999 -5.306000
5.455000 -6.294999 -3.526000 5.455000 -6.294999 -3.526000
2.430000 -8.819000 -1.027000 2.430000 -8.819000 -1.027000
-2.385000 -13.817000 3.039000 -2.385000 -13.817000 3.039000
-5.365000 -12.768000 5.181000 -5.365000 -12.768000 5.181000
-7.149000 -10.465000 6.311990 -7.149000 -10.465000 6.311990
-10.096000 -3.201001 7.867990 -10.096000 -3.201001 7.867990
-12.000000 -1.260001 9.121990 -12.000000 -1.260001 9.121990
-9.905000 5.589999 6.929990 -9.905000 5.589999 6.929990
-7.567000 9.569999 4.814990 -7.567000 9.569999 4.814990
-5.297000 13.893000 2.719990 -5.297000 13.893000 2.719990
-2.638000 13.216000 0.786000 -2.638000 13.216000 0.786000
1.167000 12.834000 -2.036000 1.167000 12.834000 -2.036000
6.422000 10.036001 -5.728000 6.422000 10.036001 -5.728000
8.310000 8.429001 -7.000000 8.310000 8.429001 -7.000000
10.086000 3.481001 -7.885000 10.086000 3.481001 -7.885000
12.079000 0.035001 -9.070000 12.079000 0.035001 -9.070000
9.618000 -8.082999 -6.489000 9.618000 -8.082999 -6.489000
7.433000 -10.827999 -4.601000 7.433000 -10.827999 -4.601000
3.905000 -13.844999 -1.680000 3.905000 -13.844999 -1.680000
26.346010 -45.625997 -15.651000 26.346010 -45.625997 -15.651000
35.480000 30.151004 -29.357000 35.480000 30.151004 -29.357000
2.383000 -17.795000 -0.180000 2.383000 -17.795000 -0.180000
0.956000 -20.194000 1.108000 0.956000 -20.194000 1.108000
-10.831000 -14.471001 9.437990 -10.831000 -14.471001 9.437990
-15.871000 0.644999 11.854990 -15.871000 0.644999 11.854990
-16.061000 5.561999 11.552990 -16.061000 5.561999 11.552990
-9.738000 14.317999 6.015990 -9.738000 14.317999 6.015990
0.362000 17.801000 -1.880000 0.362000 17.801000 -1.880000
6.711000 18.121001 -6.675000 6.711000 18.121001 -6.675000
13.383000 5.942001 -10.583000 13.383000 5.942001 -10.583000
15.901000 -1.739998 -11.778000 15.901000 -1.739998 -11.778000
12.772000 -10.645999 -8.624000 12.772000 -10.645999 -8.624000
13.779000 -18.775998 -8.646000 13.779000 -18.775998 -8.646000
4.690000 -24.507999 -1.305000 4.690000 -24.507999 -1.305000
-5.623000 -23.914000 6.382000 -5.623000 -23.914000 6.382000
-10.310000 -21.350001 9.667990 -10.310000 -21.350001 9.667990
-16.053000 -16.531001 13.542990 -16.053000 -16.531001 13.542990
-19.186000 0.806998 14.328990 -19.186000 0.806998 14.328990
-18.657000 8.687998 13.218990 -18.657000 8.687998 13.218990
-13.662000 20.673999 8.386990 -13.662000 20.673999 8.386990
-3.015000 24.964000 0.006990 -3.015000 24.964000 0.006990
4.425000 25.056001 -5.586000 4.425000 25.056001 -5.586000
12.433000 18.402001 -10.995000 12.433000 18.402001 -10.995000
20.879000 5.693002 -16.187000 20.879000 5.693002 -16.187000
19.922000 -11.873998 -13.880000 19.922000 -11.873998 -13.880000
19.705000 -25.616998 -12.476000 19.705000 -25.616998 -12.476000
8.502000 -29.983999 -3.671000 8.502000 -29.983999 -3.671000
-3.142000 -30.709000 5.134000 -3.142000 -30.709000 5.134000
-9.535000 -29.650001 9.836990 -9.535000 -29.650001 9.836990
-20.859000 -16.397002 17.138990 -20.859000 -16.397002 17.138990
-27.332000 0.111998 20.504990 -27.332000 0.111998 20.504990
-25.986000 12.149998 18.406990 -25.986000 12.149998 18.406990
-20.153000 17.626998 13.533990 -20.153000 17.626998 13.533990
-14.392000 29.950999 8.095990 -14.392000 29.950999 8.095990
0.340000 33.765000 -3.307000 0.340000 33.765000 -3.307000
7.430000 31.869001 -8.458000 7.430000 31.869001 -8.458000
18.888000 24.057002 -16.352000 18.888000 24.057002 -16.352000
21.562000 16.387002 -17.666000 21.562000 16.387002 -17.666000
26.629000 -3.488997 -19.673000 26.629000 -3.488997 -19.673000
23.831000 -9.012998 -17.073000 23.831000 -9.012998 -17.073000
28.186000 -12.432997 -20.033000 28.186000 -12.432997 -20.033000
7.964000 -40.944999 -2.277000 7.964000 -40.944999 -2.277000
-12.370000 -35.929001 12.532990 -12.370000 -35.929001 12.532990
-26.100000 -24.601002 21.813990 -26.100000 -24.601002 21.813990
-30.937000 -9.427003 24.073990 -30.937000 -9.427003 24.073990
-20.978000 33.044998 12.758990 -20.978000 33.044998 12.758990
-3.429000 36.392000 -0.715010 -3.429000 36.392000 -0.715010
9.389000 38.963001 -10.569000 9.389000 38.963001 -10.569000
30.713000 3.690003 -23.387000 30.713000 3.690003 -23.387000
]
HBA
(0,2) [ 3826583.234730000 460955.803083000 5064894.225 ]
(0,47) x (0,1) x (0,2) [
8.775000 5.195001 -7.056000 8.775000 5.195001 -7.056000
11.609000 8.730001 -9.504000 11.609000 8.730001 -9.504000
8.937000 -2.066999 -6.522000 8.937000 -2.066999 -6.522000
11.771000 1.468001 -8.969000 11.771000 1.468001 -8.969000
14.606000 5.004001 -11.416000 14.606000 5.004001 -11.416000
17.440000 8.538001 -13.863000 17.440000 8.538001 -13.863000
9.098000 -9.328999 -5.986000 9.098000 -9.328999 -5.986000
11.932000 -5.793999 -8.434000 11.932000 -5.793999 -8.434000
14.768000 -2.257999 -10.881000 14.768000 -2.257999 -10.881000
17.602000 1.276002 -13.329000 17.602000 1.276002 -13.329000
20.437000 4.812002 -15.776000 20.437000 4.812002 -15.776000
23.271000 8.347002 -18.223000 23.271000 8.347002 -18.223000
12.095000 -13.054999 -7.899000 12.095000 -13.054999 -7.899000
14.929000 -9.519999 -10.346000 14.929000 -9.519999 -10.346000
17.764000 -5.983998 -12.793000 17.764000 -5.983998 -12.793000
20.598000 -2.449998 -15.241000 20.598000 -2.449998 -15.241000
23.434000 1.086002 -17.688000 23.434000 1.086002 -17.688000
26.268000 4.621002 -20.136000 26.268000 4.621002 -20.136000
17.926000 -13.245998 -12.259000 17.926000 -13.245998 -12.259000
20.760000 -9.711998 -14.706000 20.760000 -9.711998 -14.706000
23.595000 -6.175998 -17.153000 23.595000 -6.175998 -17.153000
26.429000 -2.640998 -19.600000 26.429000 -2.640998 -19.600000
23.757000 -13.437998 -16.618000 23.757000 -13.437998 -16.618000
26.591000 -9.902998 -19.066000 26.591000 -9.902998 -19.066000
-27.572000 -1.431003 20.825990 -27.572000 -1.431003 20.825990
-28.188000 3.681997 20.825990 -28.188000 3.681997 20.825990
-22.878000 -6.053002 17.719990 -22.878000 -6.053002 17.719990
-23.494000 -0.940002 17.719990 -23.494000 -0.940002 17.719990
-24.110000 4.172998 17.719990 -24.110000 4.172998 17.719990
-24.726010 9.285997 17.719990 -24.726010 9.285997 17.719990
-18.183000 -10.674002 14.614000 -18.183000 -10.674002 14.614000
-18.799000 -5.561002 14.613990 -18.799000 -5.561002 14.613990
-19.415000 -0.448002 14.613990 -19.415000 -0.448002 14.613990
-20.031000 4.664998 14.613990 -20.031000 4.664998 14.613990
-20.647000 9.777998 14.613990 -20.647000 9.777998 14.613990
-21.263010 14.890998 14.613990 -21.263010 14.890998 14.613990
-14.105000 -10.183002 11.508000 -14.105000 -10.183002 11.508000
-14.721000 -5.070002 11.508000 -14.721000 -5.070002 11.508000
-15.337000 0.042998 11.508000 -15.337000 0.042998 11.508000
-15.953000 5.155998 11.508000 -15.953000 5.155998 11.508000
-16.569000 10.268998 11.507990 -16.569000 10.268998 11.507990
-17.185010 15.381998 11.507990 -17.185010 15.381998 11.507990
-10.642000 -4.578001 8.402000 -10.642000 -4.578001 8.402000
-11.258000 0.534999 8.402000 -11.258000 0.534999 8.402000
-11.874000 5.647999 8.402000 -11.874000 5.647999 8.402000
-12.490000 10.760999 8.402000 -12.490000 10.760999 8.402000
-7.180000 1.025999 5.296000 -7.180000 1.025999 5.296000
-7.796000 6.138999 5.296000 -7.796000 6.138999 5.296000
]
NORMAL_VECTOR HBA0
(0,2) [ 0.598753 0.072099 0.797682 ]
ROTATION_MATRIX HBA0
(0,2) x (0,2) [
-0.1195950000 -0.7919540000 0.5987530000
0.9928230000 -0.0954190000 0.0720990000
0.0000330000 0.6030780000 0.7976820000
]
HBA0
(0,2) [ 3826600.917730000 460953.449085000 5064881.164 ]
NORMAL_VECTOR HBA1
(0,2) [ 0.598753 0.072099 0.797682 ]
ROTATION_MATRIX HBA1
(0,2) x (0,2) [
-0.1195950000 -0.7919540000 0.5987530000
0.9928230000 -0.0954190000 0.0720990000
0.0000330000 0.6030780000 0.7976820000
]
HBA1
(0,2) [ 3826565.550730000 460958.157081000 5064907.286 ]
#
# HBADeltas for CS002
# Created: 2015-01-20 21:24:05
#
HBADeltas
(0,31) x (0,2) [
-2.123 0.070 1.587
-1.435 0.928 0.993
-0.747 1.786 0.399
-0.059 2.644 -0.195
-1.396 -0.835 1.123
-0.708 0.023 0.529
-0.020 0.881 -0.065
0.668 1.739 -0.659
-0.668 -1.739 0.659
0.020 -0.881 0.065
0.708 -0.023 -0.529
1.396 0.835 -1.123
0.059 -2.644 0.195
0.747 -1.786 -0.399
1.435 -0.928 -0.993
2.123 -0.070 -1.587
-1.261 -2.040 1.131
-1.410 -0.799 1.131
-1.560 0.442 1.131
-1.709 1.683 1.131
-0.271 -1.921 0.377
-0.420 -0.680 0.377
-0.570 0.561 0.377
-0.719 1.802 0.377
0.719 -1.802 -0.377
0.570 -0.561 -0.377
0.420 0.680 -0.377
0.271 1.921 -0.377
1.709 -1.683 -1.131
1.560 -0.442 -1.131
1.410 0.799 -1.131
1.261 2.040 -1.131
]
# Blitz-0.10 formatted
#
# AntennaPositions for CS003
# ITRF2005 target_date = 2015.5
# Created: 2015-01-20 21:24:06
#
NORMAL_VECTOR LBA
(0,2) [ 0.598753 0.072099 0.797682 ]
ROTATION_MATRIX LBA
(0,2) x (0,2) [
-0.1195950000 -0.7919540000 0.5987530000
0.9928230000 -0.0954190000 0.0720990000
0.0000330000 0.6030780000 0.7976820000
]
LBA
(0,2) [ 3826516.704730000 460930.113077000 5064946.485 ]
(0,95) x (0,1) x (0,2) [
0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
-2.020000 -0.244001 1.538000 -2.020000 -0.244001 1.538000
-1.339000 2.104999 0.814000 -1.339000 2.104999 0.814000
0.800000 2.362000 -0.814000 0.800000 2.362000 -0.814000
2.019000 0.243000 -1.538000 2.019000 0.243000 -1.538000
1.338000 -2.105000 -0.814000 1.338000 -2.105000 -0.814000
-0.800000 -2.363001 0.814000 -0.800000 -2.363001 0.814000
4.553000 -1.450000 -3.287000 4.553000 -1.450000 -3.287000
2.897000 -4.710000 -1.749000 2.897000 -4.710000 -1.749000
-0.115000 -5.767000 0.607000 -0.115000 -5.767000 0.607000
-3.073000 -4.125001 2.680000 -3.073000 -4.125001 2.680000
-4.594000 -0.554001 3.498000 -4.594000 -0.554001 3.498000
-3.965000 3.276999 2.680000 -3.965000 3.276999 2.680000
-1.481000 5.574999 0.608000 -1.481000 5.574999 0.608000
1.696000 5.263000 -1.749000 1.696000 5.263000 -1.749000
4.079000 2.490000 -3.287000 4.079000 2.490000 -3.287000
6.847000 3.253000 -5.434000 6.847000 3.253000 -5.434000
7.913000 -0.211000 -5.920000 7.913000 -0.211000 -5.920000
7.375000 -5.149000 -5.071000 7.375000 -5.149000 -5.071000
2.948000 -9.564000 -1.349000 2.948000 -9.564000 -1.349000
0.433000 -10.138000 0.591000 0.433000 -10.138000 0.591000
-2.952000 -9.156001 3.043000 -2.952000 -9.156001 3.043000
-5.357000 -5.986001 4.562000 -5.357000 -5.986001 4.562000
-7.148000 -0.542001 5.415000 -7.148000 -0.542001 5.415000
-7.350000 4.229999 5.134000 -7.350000 4.229999 5.134000
-5.455000 7.778999 3.392000 -5.455000 7.778999 3.392000
-2.483000 8.912999 1.058000 -2.483000 8.912999 1.058000
1.400000 8.743000 -1.841000 1.400000 8.743000 -1.841000
4.813000 6.620000 -4.211000 4.813000 6.620000 -4.211000
10.927000 3.685001 -8.535000 10.927000 3.685001 -8.535000
11.827000 -0.009999 -8.877000 11.827000 -0.009999 -8.877000
11.198000 -2.978999 -8.136000 11.198000 -2.978999 -8.136000
7.750000 -9.454000 -4.963000 7.750000 -9.454000 -4.963000
7.439000 -12.390000 -4.464000 7.439000 -12.390000 -4.464000
1.523000 -13.231000 0.053000 1.523000 -13.231000 0.053000
-2.527000 -12.509001 3.027000 -2.527000 -12.509001 3.027000
-6.779000 -12.016001 6.175000 -6.779000 -12.016001 6.175000
-7.762000 -8.838001 6.625000 -7.762000 -8.838001 6.625000
-9.577000 -4.555001 7.600000 -9.577000 -4.555001 7.600000
-10.502000 2.387999 7.667000 -10.502000 2.387999 7.667000
-10.415000 5.155999 7.351000 -10.415000 5.155999 7.351000
-7.941010 9.314999 5.118000 -7.941010 9.314999 5.118000
-6.630010 13.028999 3.799000 -6.630010 13.028999 3.799000
0.369990 14.050000 -1.547000 0.369990 14.050000 -1.547000
3.476000 12.934000 -3.778000 3.476000 12.934000 -3.778000
7.506000 10.491000 -6.583000 7.506000 10.491000 -6.583000
-38.650000 -23.644004 31.149000 -38.650000 -23.644004 31.149000
-21.152010 49.524997 11.401000 -21.152010 49.524997 11.401000
11.090000 10.637001 -9.286000 11.090000 10.637001 -9.286000
13.541000 10.183001 -11.085000 13.541000 10.183001 -11.085000
16.002000 -5.139999 -11.547000 16.002000 -5.139999 -11.547000
8.230000 -17.434000 -4.602000 8.230000 -17.434000 -4.602000
4.910000 -19.866000 -1.890000 4.910000 -19.866000 -1.890000
-4.647000 -17.006001 5.025000 -4.647000 -17.006001 5.025000
-12.595000 -7.676002 10.148000 -12.595000 -7.676002 10.148000
-16.291000 -0.964002 12.315000 -16.291000 -0.964002 12.315000
-11.457010 11.761998 7.537000 -11.457010 11.761998 7.537000
-7.484010 17.962999 3.994000 -7.484010 17.962999 3.994000
0.428990 18.618000 -2.005000 0.428990 18.618000 -2.005000
5.539990 23.390000 -6.273000 5.539990 23.390000 -6.273000
14.504000 16.170001 -12.349000 14.504000 16.170001 -12.349000
19.730000 4.763001 -15.240000 19.730000 4.763001 -15.240000
20.507000 -1.460999 -15.261000 20.507000 -1.460999 -15.261000
20.291000 -9.846999 -14.341000 20.291000 -9.846999 -14.341000
9.930000 -21.088000 -5.548000 9.930000 -21.088000 -5.548000
4.152000 -24.087000 -0.940000 4.152000 -24.087000 -0.940000
-6.928000 -24.124001 7.380000 -6.928000 -24.124001 7.380000
-15.737000 -14.568002 13.129000 -15.737000 -14.568002 13.129000
-19.871000 -6.575002 15.509000 -19.871000 -6.575002 15.509000
-19.616010 5.088998 14.264000 -19.616010 5.088998 14.264000
-15.383010 19.969998 9.742000 -15.383010 19.969998 9.742000
-2.625010 26.896999 -0.461000 -2.625010 26.896999 -0.461000
7.063990 32.890000 -8.275000 7.063990 32.890000 -8.275000
16.233990 22.768001 -14.243000 16.233990 22.768001 -14.243000
23.106000 10.521002 -18.295000 23.106000 10.521002 -18.295000
25.864000 3.137002 -19.698000 25.864000 3.137002 -19.698000
22.826000 -15.097998 -15.769000 22.826000 -15.097998 -15.769000
14.869000 -29.569999 -8.488000 14.869000 -29.569999 -8.488000
5.748000 -33.571000 -1.281000 5.748000 -33.571000 -1.281000
-1.256000 -29.753001 3.632000 -1.256000 -29.753001 3.632000
-12.990000 -29.115002 12.382000 -12.990000 -29.115002 12.382000
-23.702000 -14.934003 19.141000 -23.702000 -14.934003 19.141000
-26.259000 -6.416003 20.290000 -26.259000 -6.416003 20.290000
-27.085010 9.498997 19.472000 -27.085010 9.498997 19.472000
-23.205010 15.861997 15.984000 -23.205010 15.861997 15.984000
-12.133010 30.339998 6.365000 -12.133010 30.339998 6.365000
-6.755010 29.821999 2.375000 -6.755010 29.821999 2.375000
-6.755010 36.074999 1.810000 -6.755010 36.074999 1.810000
24.160990 27.156002 -20.590000 24.160990 27.156002 -20.590000
31.787000 2.920003 -24.124000 31.787000 2.920003 -24.124000
31.406000 -17.041997 -22.034000 31.406000 -17.041997 -22.034000
23.484000 -29.140998 -14.993000 23.484000 -29.140998 -14.993000
-11.543000 -37.630002 12.065000 -11.543000 -37.630002 12.065000
-23.470000 -20.195003 19.442000 -23.470000 -20.195003 19.442000
-32.270000 -7.515004 24.902000 -32.270000 -7.515004 24.902000
-19.366010 31.498998 11.690000 -19.366010 31.498998 11.690000
]
HBA
(0,2) [ 3826494.536720000 461017.745074000 5064955.204 ]
(0,47) x (0,1) x (0,2) [
-33.585000 -14.570003 26.526000 -33.585000 -14.570003 26.526000
-32.080000 -9.895003 24.973000 -32.080000 -9.895003 24.973000
-31.251000 -21.375003 25.389000 -31.251000 -21.375003 25.389000
-29.745000 -16.701002 23.836000 -29.745000 -16.701002 23.836000
-28.240000 -12.026002 22.284000 -28.240000 -12.026002 22.284000
-26.734000 -7.353002 20.731000 -26.734000 -7.353002 20.731000
-28.917000 -28.179002 24.252000 -28.917000 -28.179002 24.252000
-27.411000 -23.506002 22.699000 -27.411000 -23.506002 22.699000
-25.905000 -18.832002 21.147000 -25.905000 -18.832002 21.147000
-24.400000 -14.157002 19.594000 -24.400000 -14.157002 19.594000
-22.894000 -9.484002 18.041000 -22.894000 -9.484002 18.041000
-21.388000 -4.810002 16.488000 -21.388000 -4.810002 16.488000
-25.076000 -30.310002 21.562000 -25.076000 -30.310002 21.562000
-23.570000 -25.636002 20.009000 -23.570000 -25.636002 20.009000
-22.064000 -20.963002 18.456000 -22.064000 -20.963002 18.456000
-20.559000 -16.288002 16.903000 -20.559000 -16.288002 16.903000
-19.053000 -11.614001 15.351000 -19.053000 -11.614001 15.351000
-17.547000 -6.941001 13.798000 -17.547000 -6.941001 13.798000
-19.730000 -27.767002 17.319000 -19.730000 -27.767002 17.319000
-18.224000 -23.094001 15.766000 -18.224000 -23.094001 15.766000
-16.719000 -18.419001 14.214000 -16.719000 -18.419001 14.214000
-15.213000 -13.745001 12.661000 -15.213000 -13.745001 12.661000
-14.384000 -25.225001 13.077000 -14.384000 -25.225001 13.077000
-12.879000 -20.550001 11.524000 -12.879000 -20.550001 11.524000
14.324000 25.109002 -13.021000 14.324000 25.109002 -13.021000
17.157990 28.644002 -15.469000 17.157990 28.644002 -15.469000
14.486000 17.847002 -12.487000 14.486000 17.847002 -12.487000
17.320000 21.382002 -14.934000 17.320000 21.382002 -14.934000
20.155000 24.918002 -17.381000 20.155000 24.918002 -17.381000
22.989000 28.452003 -19.828000 22.989000 28.452003 -19.828000
14.647000 10.585002 -11.951000 14.647000 10.585002 -11.951000
17.481000 14.120002 -14.399000 17.481000 14.120002 -14.399000
20.317000 17.656002 -16.846000 20.317000 17.656002 -16.846000
23.151000 21.190003 -19.294000 23.151000 21.190003 -19.294000
25.986000 24.726003 -21.741000 25.986000 24.726003 -21.741000
28.820000 28.261003 -24.188000 28.820000 28.261003 -24.188000
17.644000 6.859002 -13.864000 17.644000 6.859002 -13.864000
20.478000 10.394002 -16.311000 20.478000 10.394002 -16.311000
23.313000 13.930003 -18.758000 23.313000 13.930003 -18.758000
26.147000 17.464003 -21.206000 26.147000 17.464003 -21.206000
28.983000 21.000003 -23.653000 28.983000 21.000003 -23.653000
31.817000 24.535003 -26.101000 31.817000 24.535003 -26.101000
23.475000 6.668003 -18.224000 23.475000 6.668003 -18.224000
26.309000 10.202003 -20.671000 26.309000 10.202003 -20.671000
29.144000 13.738003 -23.118000 29.144000 13.738003 -23.118000
31.978000 17.273003 -25.565000 31.978000 17.273003 -25.565000
29.306000 6.476003 -22.583000 29.306000 6.476003 -22.583000
32.140000 10.011003 -25.031000 32.140000 10.011003 -25.031000
]
NORMAL_VECTOR HBA0
(0,2) [ 0.598753 0.072099 0.797682 ]
ROTATION_MATRIX HBA0
(0,2) x (0,2) [
-0.1195950000 -0.7919540000 0.5987530000
0.9928230000 -0.0954190000 0.0720990000
0.0000330000 0.6030780000 0.7976820000
]
HBA0
(0,2) [ 3826471.304720000 461000.185072000 5064974.229 ]
NORMAL_VECTOR HBA1
(0,2) [ 0.598753 0.072099 0.797682 ]
ROTATION_MATRIX HBA1
(0,2) x (0,2) [
-0.1195950000 -0.7919540000 0.5987530000
0.9928230000 -0.0954190000 0.0720990000
0.0000330000 0.6030780000 0.7976820000
]
HBA1
(0,2) [ 3826517.768720000 461035.305077000 5064936.178 ]
#
# HBADeltas for CS003
# Created: 2015-01-20 21:24:06
#
HBADeltas
(0,31) x (0,2) [
-1.946 -0.926 1.545
-1.581 0.209 1.168
-1.215 1.343 0.791
-0.850 2.477 0.414
-1.014 -1.443 0.892
-0.649 -0.309 0.515
-0.283 0.826 0.138
0.082 1.960 -0.239
-0.082 -1.960 0.239
0.283 -0.826 -0.138
0.649 0.309 -0.515
1.014 1.443 -0.892
0.850 -2.477 -0.414
1.215 -1.343 -0.791
1.581 -0.209 -1.168
1.946 0.926 -1.545
-2.123 0.070 1.587
-1.435 0.928 0.993
-0.747 1.786 0.399
-0.059 2.644 -0.195
-1.396 -0.835 1.123
-0.708 0.023 0.529
-0.020 0.881 -0.065
0.668 1.739 -0.659
-0.668 -1.739 0.659
0.020 -0.881 0.065
0.708 -0.023 -0.529
1.396 0.835 -1.123
0.059 -2.644 0.195
0.747 -1.786 -0.399
1.435 -0.928 -0.993
2.123 -0.070 -1.587
]
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment