From e403b3f9b248818625f74e1e2bd82254cb4c22f9 Mon Sep 17 00:00:00 2001
From: jay <jlaura@usgs.gov>
Date: Fri, 24 Aug 2018 07:53:39 -0700
Subject: [PATCH] Removes dirs
---
plio/camera/__init__.py | 20 -
plio/camera/csm.py | 37 --
plio/camera/tests/__init__.py | 0
plio/camera/tests/test_csm.py | 45 --
plio/spatial/__init__.py | 0
plio/spatial/footprint.py | 59 ---
plio/spatial/tests/test_transformations.py | 74 ----
plio/spatial/transformations.py | 469 ---------------------
8 files changed, 704 deletions(-)
delete mode 100644 plio/camera/__init__.py
delete mode 100644 plio/camera/csm.py
delete mode 100644 plio/camera/tests/__init__.py
delete mode 100644 plio/camera/tests/test_csm.py
delete mode 100644 plio/spatial/__init__.py
delete mode 100644 plio/spatial/footprint.py
delete mode 100644 plio/spatial/tests/test_transformations.py
delete mode 100644 plio/spatial/transformations.py
diff --git a/plio/camera/__init__.py b/plio/camera/__init__.py
deleted file mode 100644
index 3f4336c..0000000
--- a/plio/camera/__init__.py
+++ /dev/null
@@ -1,20 +0,0 @@
-import importlib
-import warnings
-
-cam = importlib.find_loader('usgscam')
-cycsm_isd = importlib.find_loader('cycsm.isd')
-
-if cam:
- cam = cam.load_module()
-
-if cycsm_isd:
- cycsm_isd = cycsm_isd.load_module()
-
-def conditional_cameras(func):
- def cam_check(*args, **kwargs):
- if cam:
- return func(*args, **kwargs)
- else:
- warning.warn('Trying to call a camera method, but usgscam is not installed.')
- return None
- return cam_check
\ No newline at end of file
diff --git a/plio/camera/csm.py b/plio/camera/csm.py
deleted file mode 100644
index 6e20e79..0000000
--- a/plio/camera/csm.py
+++ /dev/null
@@ -1,37 +0,0 @@
-import datetime
-import json
-
-import requests
-
-from plio.utils.utils import find_in_dict
-from plio.io.io_json import NumpyEncoder
-from plio.camera import conditional_cameras, cam, cycsm_isd
-
-def data_from_cube(header):
- data = {}
- data['START_TIME'] = find_in_dict(header, 'StartTime')
- data['SPACECRAFT_NAME'] = find_in_dict(header, 'SpacecraftName')
- data['INSTRUMENT_NAME'] = find_in_dict(header, 'InstrumentId')
- data['SAMPLING_FACTOR'] = find_in_dict(header, 'SpatialSumming')
- data['SAMPLE_FIRST_PIXEL'] = find_in_dict(header, 'SampleFirstPixel')
- data['IMAGE'] = {}
- data['IMAGE']['LINES'] = find_in_dict(header, 'Lines')
- data['IMAGE']['SAMPLES'] = find_in_dict(header, 'Samples')
- data['TARGET_NAME'] = find_in_dict(header, 'TargetName')
- data['LINE_EXPOSURE_DURATION'] = find_in_dict(header, 'LineExposureDuration')
- data['SPACECRAFT_CLOCK_START_COUNT'] = find_in_dict(header, 'SpacecraftClockCount')
- return data
-
-@conditional_cameras
-def create_camera(obj, url='http://smalls:8002/api/1.0/missions/mars_reconnaissance_orbiter/csm_isd'):
-
- data = json.dumps(data_from_cube(obj.metadata), cls=NumpyEncoder)
- r = requests.post(url, data=data)
-
- # Get the ISD back and instantiate a local ISD for the image
- isd = r.json()['data']['isd']
- i = cycsm_isd.Isd.loads(isd)
-
- # Create the plugin and camera as usual
- plugin = cam.genericls.Plugin()
- return plugin.from_isd(i, plugin.modelname(0))
\ No newline at end of file
diff --git a/plio/camera/tests/__init__.py b/plio/camera/tests/__init__.py
deleted file mode 100644
index e69de29..0000000
diff --git a/plio/camera/tests/test_csm.py b/plio/camera/tests/test_csm.py
deleted file mode 100644
index 7f2cd90..0000000
--- a/plio/camera/tests/test_csm.py
+++ /dev/null
@@ -1,45 +0,0 @@
-import datetime
-import json
-from unittest import mock
-
-import pytest
-import pvl
-
-from plio.camera import csm, cam, conditional_cameras
-from plio.examples import get_path
-
-def mock_requests_post(*args, **kwargs):
- class MockResponse:
- def __init__(self, json_data, status_code):
- self.json_data = json_data
- self.status_code = status_code
-
- def json(self):
- return self.json_data
-
- if 'mars_reconnaissance_orbiter' in args[0]:
- with open(get_path('ctx.response'), 'r') as f:
- resp = json.load(f)
- return MockResponse(resp, 200)
-
- return MockResponse(None, 404)
-
-
-@pytest.fixture
-def header():
- return pvl.load(get_path('ctx.pvl'))
-
-@pytest.fixture
-def req_obj():
- return
-
-@pytest.mark.skipif(cam is None, reason="Cameras not installed")
-def test_data_from_cube(header):
- data = csm.data_from_cube(header)
- assert data['START_TIME'] == datetime.datetime(2008, 9, 17, 5, 8, 10, 820000)
-
-@pytest.mark.skipif(cam is None, reason="Cameras not installed")
-@mock.patch('requests.post', side_effect=mock_requests_post)
-def test_create_camera(header):
- created_camera = csm.create_camera(header)
- assert isinstance(create_camera, cam.genericls.SensorModel)
diff --git a/plio/spatial/__init__.py b/plio/spatial/__init__.py
deleted file mode 100644
index e69de29..0000000
diff --git a/plio/spatial/footprint.py b/plio/spatial/footprint.py
deleted file mode 100644
index 43ac649..0000000
--- a/plio/spatial/footprint.py
+++ /dev/null
@@ -1,59 +0,0 @@
-import numpy as np
-import pyproj
-import ogr
-
-def generate_gcps(camera, nnodes=5, semi_major=3396190, semi_minor=3376200):
- ecef = pyproj.Proj(proj='geocent', a=semi_major, b=semi_minor)
- lla = pyproj.Proj(proj='latlon', a=semi_major, b=semi_minor)
-
- isize = camera.imagesize[::-1]
- x = np.linspace(0,isize[1], 10)
- y = np.linspace(0,isize[0], 10)
- boundary = [(i,0.) for i in x] + [(isize[1], i) for i in y[1:]] +\
- [(i, isize[0]) for i in x[::-1][1:]] + [(0.,i) for i in y[::-1][1:]]
- gnds = np.empty((len(boundary), 3))
- for i, b in enumerate(boundary):
- gnds[i] = camera.imageToGround(*b, 0)
- lons, lats, alts = pyproj.transform(ecef, lla, gnds[:,0], gnds[:,1], gnds[:,2])
- lla = np.vstack((lons, lats, alts)).T
-
- tr = zip(boundary, lla)
-
- gcps = []
- for i, t in enumerate(tr):
- l = '<GCP Id="{}" Info="{}" Pixel="{}" Line="{}" X="{}" Y="{}" Z="{}" />'.format(i, i, t[0][1], t[0][0], t[1][0], t[1][1], t[1][2])
- gcps.append(l)
-
- return gcps
-
-def generate_latlon_footprint(camera, nnodes=5, semi_major=3396190, semi_minor=3376200):
- ecef = pyproj.Proj(proj='geocent', a=semi_major, b=semi_minor)
- lla = pyproj.Proj(proj='latlon', a=semi_major, b=semi_minor)
-
- isize = camera.imagesize[::-1]
- x = np.linspace(0,isize[1], 10)
- y = np.linspace(0,isize[0], 10)
- boundary = [(i,0.) for i in x] + [(isize[1], i) for i in y[1:]] +\
- [(i, isize[0]) for i in x[::-1][1:]] + [(0.,i) for i in y[::-1][1:]]
- ring = ogr.Geometry(ogr.wkbLinearRing)
- for i in boundary:
- gnd = camera.imageToGround(*i, 0)
- lons, lats, alts = pyproj.transform(ecef, lla, gnd[0], gnd[1], gnd[2])
- ring.AddPoint(lons, lats)
- poly = ogr.Geometry(ogr.wkbPolygon)
- poly.AddGeometry(ring)
- return poly
-
-def generate_bodyfixed_footprint(camera, nnodes=5):
- isize = camera.imagesize[::-1]
- x = np.linspace(0,isize[1], 10)
- y = np.linspace(0,isize[0], 10)
- boundary = [(i,0.) for i in x] + [(isize[1], i) for i in y[1:]] +\
- [(i, isize[0]) for i in x[::-1][1:]] + [(0.,i) for i in y[::-1][1:]]
- ring = ogr.Geometry(ogr.wkbLinearRing)
- for i in boundary:
- gnd = camera.imageToGround(*i, 0)
- ring.AddPoint(gnd[0], gnd[1], gnd[2])
- poly = ogr.Geometry(ogr.wkbPolygon)
- poly.AddGeometry(ring)
- return poly
diff --git a/plio/spatial/tests/test_transformations.py b/plio/spatial/tests/test_transformations.py
deleted file mode 100644
index 5c37fa9..0000000
--- a/plio/spatial/tests/test_transformations.py
+++ /dev/null
@@ -1,74 +0,0 @@
-
-import pytest
-from plio.examples import get_path
-from plio.spatial import transformations as trans
-
-@pytest.mark.parametrize("stat, expected",
- [({'stat':0}, True),
- ({'stat':1}, False),
- ({'stat':True}, False), # Is this the desired result?
- ({'stat':False}, True)])
-def test_stat_toggle(stat, expected):
- assert trans.stat_toggle(stat) == expected
-
-@pytest.mark.parametrize("stat, expected",
- [({'stat':True}, 0),
- ({'stat':False}, 1),
- ({'stat':False}, True), # Is this the desired result?
- ({'stat':True}, False)])
-def test_ignore_toggle(stat, expected):
- assert trans.ignore_toggle(stat) == expected
-
-def test_get_axis():
- fname = get_path('CTX_Athabasca_Middle.prj')
- erad, prad = trans.get_axis(fname)
- assert erad == 3.39619000000000e+006
- assert prad == 3.3762000000000866e+006
-
-@pytest.mark.parametrize("og, major, minor, oc",
- [(0, 3396190, 3376200, 0),
- (90, 3396190, 3376200, 90),
- (-90, 3396190, 3376200, -90),
- (45, 3396190, 3376200, 44.6617680)])
-def test_og2oc(og, major, minor, oc):
- assert trans.og2oc(og, major, minor) == pytest.approx(oc)
-
-@pytest.mark.parametrize("og, major, minor, oc",
- [(0, 3396190, 3376200, 0),
- (90, 3396190, 3376200, 90),
- (-90, 3396190, 3376200, -90),
- (45.338231, 3396190, 3376200, 45)])
-def test_oc2og(og, major, minor, oc):
- assert trans.oc2og(oc, major, minor) == pytest.approx(og)
-
-
-@pytest.mark.parametrize("known, expected",
- [({'known':0},'Free'),
- ({'known':1},'Constrained'),
- ({'known':2},'Constrained'),
- ({'known':3},'Constrained'),
- ({'known':False},'Free'), # Is this the desired result?
- ({'known':True},'Constrained')])
-def test_known(known, expected):
- assert trans.known(known) == expected
-
-
-@pytest.mark.parametrize("known, expected",
- [({'known':0},0),
- ({'known':2},0),
- ({'known':1},3),
- ({'known':3},3),
- ({'known':4},3),
- ({'known':False},0), # Is this the desired result?
- ({'known':True},3)])
-def test_known(known, expected):
- assert trans.reverse_known(known) == expected
-
-@pytest.mark.parametrize("num, expected",
- [(0, 0),
- (-180, 180),
- (370, 10),
- (-90, 270),
- (90, 90)])
-def test_to_360(num, expected):
- assert trans.to_360(num) == expected
\ No newline at end of file
diff --git a/plio/spatial/transformations.py b/plio/spatial/transformations.py
deleted file mode 100644
index 9dc4a67..0000000
--- a/plio/spatial/transformations.py
+++ /dev/null
@@ -1,469 +0,0 @@
-import os
-import pvl
-import math
-import pyproj
-
-import numpy as np
-
-import plio.io.isis_serial_number as sn
-from plio.utils.utils import find_in_dict
-
-def line_sample_size(record, path):
- """
- Converts columns l. and s. to sample size, line size, and generates an
- image index
-
- Parameters
- ----------
- record : object
- Pandas series object
-
- path : str
- Path to the associated sup files for a socet project
-
- Returns
- -------
- : list
- A list of sample_size, line_size, and img_index
- """
- with open(os.path.join(path, record['ipf_file'] + '.sup')) as f:
- for i, line in enumerate(f):
- if i == 2:
- img_index = line.split('\\')
- img_index = img_index[-1].strip()
- img_index = img_index.split('.')[0]
-
- if i == 3:
- line_size = line.split(' ')
- line_size = line_size[-1].strip()
- assert int(line_size) > 0, "Line number {} from {} is a negative number: Invalid Data".format(line_size, record['ipf_file'])
-
- if i == 4:
- sample_size = line.split(' ')
- sample_size = sample_size[-1].strip()
- assert int(sample_size) > 0, "Sample number {} from {} is a negative number: Invalid Data".format(sample_size, record['ipf_file'])
- break
-
-
- line_size = int(line_size)/2.0 + record['l.'] + 1
- sample_size = int(sample_size)/2.0 + record['s.'] + 1
- return sample_size, line_size, img_index
-
-def known(record):
- """
- Converts the known field from a socet dataframe into the
- isis point_type column
-
- Parameters
- ----------
- record : object
- Pandas series object
-
- Returns
- -------
- : str
- String representation of a known field
- """
-
- lookup = {0: 'Free',
- 1: 'Constrained',
- 2: 'Constrained',
- 3: 'Constrained'}
- return lookup[record['known']]
-
-def reverse_known(record):
- """
- Converts the known field from an isis dataframe into the
- socet known column
-
- Parameters
- ----------
- record : object
- Pandas series object
-
- Returns
- -------
- : str
- String representation of a known field
- """
- lookup = {0:0,
- 2:0,
- 1:3,
- 3:3,
- 4:3}
- record_type = record['known']
- return lookup[record_type]
-
-def to_360(num):
- """
- Transforms a given number into 0 - 360 space
-
- Parameters
- ----------
- num : int
- A given integer
-
- Returns
- -------
- : int
- num moduloed by 360
- """
- return num % 360
-
-def oc2og(dlat, dMajorRadius, dMinorRadius):
- """
- Ocentric to ographic latitudes
-
- Parameters
- ----------
- dlat : float
- Latitude to convert
-
- dMajorRadius : float
- Radius from the center of the body to the equater
-
- dMinorRadius : float
- Radius from the pole to the center of mass
-
- Returns
- -------
- dlat : float
- Converted latitude into ographic space
- """
- try:
- dlat = math.radians(dlat)
- dlat = math.atan(((dMajorRadius / dMinorRadius)**2) * (math.tan(dlat)))
- dlat = math.degrees(dlat)
- except:
- print ("Error in oc2og conversion")
- return dlat
-
-def og2oc(dlat, dMajorRadius, dMinorRadius):
- """
- Ographic to ocentric latitudes
-
- Parameters
- ----------
- dlat : float
- Latitude to convert
-
- dMajorRadius : float
- Radius from the center of the body to the equater
-
- dMinorRadius : float
- Radius from the pole to the center of mass
-
- Returns
- -------
- dlat : float
- Converted latitude into ocentric space
- """
- try:
- dlat = math.radians(dlat)
- dlat = math.atan((math.tan(dlat) / ((dMajorRadius / dMinorRadius)**2)))
- dlat = math.degrees(dlat)
- except:
- print ("Error in og2oc conversion")
- return dlat
-
-def get_axis(file):
- """
- Gets eRadius and pRadius from a .prj file
-
- Parameters
- ----------
- file : str
- file with path to a given socet project file
-
- Returns
- -------
- : list
- A list of the eRadius and pRadius of the project file
- """
- with open(file) as f:
- from collections import defaultdict
-
- files = defaultdict(list)
-
- for line in f:
-
- ext = line.strip().split(' ')
- files[ext[0]].append(ext[-1])
-
- eRadius = float(files['A_EARTH'][0])
- pRadius = eRadius * math.sqrt(1 - (float(files['E_EARTH'][0]) ** 2))
-
- return eRadius, pRadius
-
-def reproject(record, semi_major, semi_minor, source_proj, dest_proj, **kwargs):
- """
- Thin wrapper around PyProj's Transform() function to transform 1 or more three-dimensional
- point from one coordinate system to another. If converting between Cartesian
- body-centered body-fixed (BCBF) coordinates and Longitude/Latitude/Altitude coordinates,
- the values input for semi-major and semi-minor axes determine whether latitudes are
- planetographic or planetocentric and determine the shape of the datum for altitudes.
- If semi_major == semi_minor, then latitudes are interpreted/created as planetocentric
- and altitudes are interpreted/created as referenced to a spherical datum.
- If semi_major != semi_minor, then latitudes are interpreted/created as planetographic
- and altitudes are interpreted/created as referenced to an ellipsoidal datum.
-
- Parameters
- ----------
- record : object
- Pandas series object
-
- semi_major : float
- Radius from the center of the body to the equater
-
- semi_minor : float
- Radius from the pole to the center of mass
-
- source_proj : str
- Pyproj string that defines a projection space ie. 'geocent'
-
- dest_proj : str
- Pyproj string that defines a project space ie. 'latlon'
-
- Returns
- -------
- : list
- Transformed coordinates as y, x, z
-
- """
- source_pyproj = pyproj.Proj(proj = source_proj, a = semi_major, b = semi_minor)
- dest_pyproj = pyproj.Proj(proj = dest_proj, a = semi_major, b = semi_minor)
-
- y, x, z = pyproj.transform(source_pyproj, dest_pyproj, record[0], record[1], record[2], **kwargs)
- return y, x, z
-
-def stat_toggle(record):
- if record['stat'] == 0:
- return True
- else:
- return False
-
-def apply_isis_transformations(df, eRadius, pRadius, serial_dict, cub_dict):
- """
- Takes a atf dictionary and a socet dataframe and applies the necessary
- transformations to convert that dataframe into a isis compatible
- dataframe
-
- Parameters
- ----------
- df : object
- Pandas dataframe object
-
- eRadius : float
- Equitorial radius
-
- pRadius : float
- Polar radius
-
- serial_dict : dict
- Dictionary mapping serials as keys to images as the values
-
- extension : str
- String extension of all cubes being used
-
- cub_path : str
- Path to all cubes being used
-
- """
- # Convert from geocentered coords (x, y, z), to lat lon coords (latitude, longitude, alltitude)
- ecef = np.array([[df['long_X_East']], [df['lat_Y_North']], [df['ht']]])
- lla = reproject(ecef, semi_major = eRadius, semi_minor = pRadius,
- source_proj = 'latlon', dest_proj = 'geocent')
-
- df['long_X_East'], df['lat_Y_North'], df['ht'] = lla[0][0], lla[1][0], lla[2][0]
-
- # Update the stat fields and add the val field as it is just a clone of stat
- df['stat'] = df.apply(ignore_toggle, axis = 1)
- df['val'] = df['stat']
-
- # Update the known field, add the ipf_file field for saving, and
- # update the line, sample using data from the cubes
- df['known'] = df.apply(reverse_known, axis = 1)
- df['ipf_file'] = df['serialnumber'].apply(lambda serial_number: serial_dict[serial_number])
- df['l.'], df['s.'] = zip(*df.apply(fix_sample_line, serial_dict = serial_dict,
- cub_dict = cub_dict, axis = 1))
-
- # Add dummy for generic value setting
- x_dummy = lambda x: np.full(len(df), x)
-
- df['sig0'] = x_dummy(1)
- df['sig1'] = x_dummy(1)
- df['sig2'] = x_dummy(1)
-
- df['res0'] = x_dummy(0)
- df['res1'] = x_dummy(0)
- df['res2'] = x_dummy(0)
-
- df['fid_x'] = x_dummy(0)
- df['fid_y'] = x_dummy(0)
-
- df['no_obs'] = x_dummy(1)
- df['fid_val'] = x_dummy(0)
-
-def apply_socet_transformations(atf_dict, df):
- """
- Takes a atf dictionary and a socet dataframe and applies the necessary
- transformations to convert that dataframe into a isis compatible
- dataframe
-
- Parameters
- ----------
- atf_dict : dict
- Dictionary containing information from an atf file
-
- df : object
- Pandas dataframe object
-
- """
- prj_file = os.path.join(atf_dict['PATH'], atf_dict['PROJECT'])
-
- eRadius, pRadius = get_axis(prj_file)
-
- lla = np.array([[df['long_X_East']], [df['lat_Y_North']], [df['ht']]])
-
- ecef = reproject(lla, semi_major = eRadius, semi_minor = pRadius,
- source_proj = 'geocent', dest_proj = 'latlon')
-
- df['s.'], df['l.'], df['image_index'] = (zip(*df.apply(line_sample_size, path = atf_dict['PATH'], axis=1)))
- df['known'] = df.apply(known, axis=1)
- df['long_X_East'] = ecef[0][0]
- df['lat_Y_North'] = ecef[1][0]
- df['ht'] = ecef[2][0]
- df['aprioriCovar'] = df.apply(compute_cov_matrix, semimajor_axis = eRadius, axis=1)
- df['stat'] = df.apply(stat_toggle, axis=1)
-
-# TODO: Does isis cnet need a convariance matrix for sigmas? Even with a static matrix of 1,1,1,1
-def compute_sigma_covariance_matrix(lat, lon, rad, latsigma, lonsigma, radsigma, semimajor_axis):
-
- """
- Given geospatial coordinates, desired accuracy sigmas, and an equitorial radius, compute a 2x3
- sigma covariange matrix.
- Parameters
- ----------
- lat : float
- A point's latitude in degrees
-
- lon : float
- A point's longitude in degrees
-
- rad : float
- The radius (z-value) of the point in meters
-
- latsigma : float
- The desired latitude accuracy in meters (Default 10.0)
-
- lonsigma : float
- The desired longitude accuracy in meters (Default 10.0)
-
- radsigma : float
- The desired radius accuracy in meters (Defualt: 15.0)
-
- semimajor_axis : float
- The semi-major or equitorial radius in meters (Default: 1737400.0 - Moon)
- Returns
- -------
- rectcov : ndarray
- (2,3) covariance matrix
- """
- lat = math.radians(lat)
- lon = math.radians(lon)
-
- # SetSphericalSigmasDistance
- scaled_lat_sigma = latsigma / semimajor_axis
-
- # This is specific to each lon.
- scaled_lon_sigma = lonsigma * math.cos(lat) / semimajor_axis
-
- # SetSphericalSigmas
- cov = np.eye(3,3)
- cov[0,0] = math.radians(scaled_lat_sigma) ** 2
- cov[1,1] = math.radians(scaled_lon_sigma) ** 2
- cov[2,2] = radsigma ** 2
-
- # Approximate the Jacobian
- j = np.zeros((3,3))
- cosphi = math.cos(lat)
- sinphi = math.sin(lat)
- cos_lmbda = math.cos(lon)
- sin_lmbda = math.sin(lon)
- rcosphi = rad * cosphi
- rsinphi = rad * sinphi
- j[0,0] = -rsinphi * cos_lmbda
- j[0,1] = -rcosphi * sin_lmbda
- j[0,2] = cosphi * cos_lmbda
- j[1,0] = -rsinphi * sin_lmbda
- j[1,1] = rcosphi * cos_lmbda
- j[1,2] = cosphi * sin_lmbda
- j[2,0] = rcosphi
- j[2,1] = 0.
- j[2,2] = sinphi
- mat = j.dot(cov)
- mat = mat.dot(j.T)
- rectcov = np.zeros((2,3))
- rectcov[0,0] = mat[0,0]
- rectcov[0,1] = mat[0,1]
- rectcov[0,2] = mat[0,2]
- rectcov[1,0] = mat[1,1]
- rectcov[1,1] = mat[1,2]
- rectcov[1,2] = mat[2,2]
-
- return rectcov
-
-def compute_cov_matrix(record, semimajor_axis):
- cov_matrix = compute_sigma_covariance_matrix(record['lat_Y_North'], record['long_X_East'], record['ht'], record['sig0'], record['sig1'], record['sig2'], semimajor_axis)
- return cov_matrix.ravel().tolist()
-
-
-
-def fix_sample_line(record, serial_dict, cub_dict):
- """
- Extracts the sample, line data from a cube and computes deviation from the
- center of the image
-
- Parameters
- ----------
- record : dict
- Dict containing the key serialnumber, l., and s.
-
- serial_dict : dict
- Maps serial numbers to images
-
- cub_dict : dict
- Maps basic cub names to their assocated absoluate path cubs
-
- Returns
- -------
- new_line : int
- new line deviation from the center
-
- new_sample : int
- new sample deviation from the center
-
- """
- # Cube location to load
- cube = pvl.load(cub_dict[serial_dict[record['serialnumber']]])
- line_size = find_in_dict(cube, 'Lines')
- sample_size = find_in_dict(cube, 'Samples')
-
- new_line = record['l.'] - (int(line_size / 2.0)) - 1
- new_sample = record['s.'] - (int(sample_size / 2.0)) - 1
-
- return new_line, new_sample
-
-def ignore_toggle(record):
- """
- Maps the stat column in a record to 0 or 1 based on True or False
-
- Parameters
- ----------
- record : dict
- Dict containing the key stat
- """
- if record['stat'] == True:
- return 0
- else:
- return 1
--
GitLab