Clean up comparison notebooks and fix vis.py

README.md

@@ -59,21 +59,3 @@ To build the docs:
 1. Browse to the `index.html` file in `docs/public` to open the docs.  If needed, run a utility like `http-server` in the `docs/public` directory to serve the docs on localhost.
 
 <hr>
-
-## Status Report - November 2019
-
-For full testing reports and example usage, please see the linked example Jupyter notebooks in the table below. 
-
-|       Instrument      |                      Jupyter Notebooks                     |  Production Ready |    Difference CSM -> ISIS (in pixels)   |    Difference ISIS -> CSM (in pixels)   |
-|:---------------------:|:-------------------------------------------------:|:-----------------:|:---------------------------------------:|:---------------------------------------:|
-|       MRO HiRISE      |     [link](examples/mro_hirise_isis_cmp.ipynb)    |     sub-pixel; in testing for production    | sample mean=-2.0e-05; line mean=2.5e-08 | sample mean=-3.0e-08; line mean=1.2e-04 |
-|        MRO CTX        |      [link](examples/mro_ctx_isis_cmp.ipynb)      | nearly sub-pixel; still in research |  gross error in line  |     sample mean=0.0002; line mean=-0.07     |
-|        MEX HRSC       |      [link](examples/mex_hrsc_isis_cmp.ipynb)     |         sub-pixel; in testing for production        |      sample mean=0.000038; line mean=-0.000072      |       sample mean=-0.000038 ; line mean=-7.512e-05       |
-|        LROC NAC       |      [link](examples/lrocnac_isis_cmp.ipynb)      | sub-pixel; in testing for production |      sample mean=-0.003; line mean=-0.0006      |                   sample mean=0.0005	line mean=0.003                   |
-| Kaguya Terrain Camera |     [link](examples/kaguya_tc_isis_cmp.ipynb)     | barely sub-pixel; in testing |     sample mean=0.0001; line mean=0.00003    |      sample mean=0.009; line mean=-1.242      |
-|   Messenger MDIS NAC  | [link](examples/messenger_mdisnac_isis_cmp.ipynb) |  sub-pixel; in testing for production |      sample mean=-0.01; line mean=-0.003     |       sample mean=0.01; line mean=0.003      |
-|  Cassini ISS NAC      |      [link](examples/cassini_isis_nac_cmp.ipynb)      |  sub-pixel; in testing for production  |    sample mean=-0.001; line mean=0.01    |     sample mean=0.001; line mean=-0.01     |
-|  Cassini ISS WAC      |      [link](examples/cassini_isis_wac_cmp.ipynb)      |  sub-pixel; in testing for production  |    sample mean=0.001; line mean=0.004    |     sample mean=-0.001; line mean=-0.004     |
-|  Dawn Framing Camera  |      [link](examples/dawn_fc_isis_cmp.ipynb)      |     sub-pixel; in testing for production    | sample mean=-0.02; line mean=0.003 | sample mean=0.02; line mean=-0.003 |
-
-The Difference column (CSM -> ISIS) represents the mean difference in pixels from running usgscsm's *image2ground* and then back to the camera using ISIS3's *campt* (ground2image). The Difference column (ISIS -> CSM) is simply the reverse starting with ISIS3 first with *campt* (image2ground) and then usgscsm's *ground2image*.

examples/data/N1702360370_1.LBL

-PDS_VERSION_ID = PDS3
-
-/* FILE CHARACTERISTICS */
-
-RECORD_TYPE = FIXED_LENGTH
-RECORD_BYTES = 1048
-FILE_RECORDS = 1028
-
-/* POINTERS TO DATA OBJECTS */
-
-^IMAGE_HEADER = ("N1702360370_1.IMG",1)
-^TELEMETRY_TABLE = ("N1702360370_1.IMG",4)
-^LINE_PREFIX_TABLE = ("N1702360370_1.IMG",5)
-^IMAGE = ("N1702360370_1.IMG",5)
-
-/* IDENTIFICATION DATA ELEMENTS */
-
-ANTIBLOOMING_STATE_FLAG = "OFF"
-BIAS_STRIP_MEAN = 8.850293
-CALIBRATION_LAMP_STATE_FLAG = "N/A"
-COMMAND_FILE_NAME = "trigger_3618_1.ioi"
-COMMAND_SEQUENCE_NUMBER = 3618
-DARK_STRIP_MEAN = 0.269221
-DATA_CONVERSION_TYPE = "TABLE"
-DATA_SET_ID = "CO-S-ISSNA/ISSWA-2-EDR-V1.0"
-DELAYED_READOUT_FLAG = "NO"
-DESCRIPTION = "Incomplete product finalized due to truncated lines."
-DETECTOR_TEMPERATURE = -89.243546 <DEGC>
-EARTH_RECEIVED_START_TIME = 2011-346T22:30:08.981
-EARTH_RECEIVED_STOP_TIME = 2011-346T22:30:49.765
-ELECTRONICS_BIAS = 112
-EXPECTED_MAXIMUM = (62.996498,69.454498)
-EXPECTED_PACKETS = 576
-EXPOSURE_DURATION = 4600.000000
-FILTER_NAME = ("CL1","UV3")
-FILTER_TEMPERATURE = 0.248629
-FLIGHT_SOFTWARE_VERSION_ID = "1.4"
-GAIN_MODE_ID = "29 ELECTRONS PER DN"
-IMAGE_MID_TIME = 2011-346T05:02:22.073
-IMAGE_NUMBER = "1702360370"
-IMAGE_OBSERVATION_TYPE = {"SCIENCE"}
-IMAGE_TIME = 2011-346T05:02:24.373
-INSTRUMENT_DATA_RATE = 182.783997
-INSTRUMENT_HOST_NAME = "CASSINI ORBITER"
-INSTRUMENT_ID = "ISSNA"
-INSTRUMENT_MODE_ID = "FULL"
-INSTRUMENT_NAME = "IMAGING SCIENCE SUBSYSTEM - NARROW ANGLE"
-INST_CMPRS_PARAM = ("N/A","N/A","N/A","N/A")
-INST_CMPRS_RATE = (5.333330,3.631307)
-INST_CMPRS_RATIO = 2.203063
-INST_CMPRS_TYPE = "LOSSLESS"
-LIGHT_FLOOD_STATE_FLAG = "ON"
-METHOD_DESC = "ISSPT2.7;Enceladus;ISS_158EN_ENCEL001_PRIME"
-MISSING_LINES = 31
-MISSING_PACKET_FLAG = "NO"
-MISSION_NAME = "CASSINI-HUYGENS"
-MISSION_PHASE_NAME = "EXTENDED-EXTENDED MISSION"
-OBSERVATION_ID = "ISS_158EN_ENCEL001_PRIME"
-OPTICS_TEMPERATURE = (0.627499,1.905708)
-ORDER_NUMBER = 16
-PARALLEL_CLOCK_VOLTAGE_INDEX = 9
-PREPARE_CYCLE_INDEX = 4
-PRODUCT_CREATION_TIME = 2011-346T15:34:07.000
-PRODUCT_ID = "1_N1702360370.120"
-PRODUCT_VERSION_TYPE = "FINAL"
-READOUT_CYCLE_INDEX = 6
-RECEIVED_PACKETS = 523
-SENSOR_HEAD_ELEC_TEMPERATURE = 1.633024
-SEQUENCE_ID = "S71"
-SEQUENCE_NUMBER = 16
-SEQUENCE_TITLE = "ISS_158EN_ENCEL001_PRIME"
-SHUTTER_MODE_ID = "NACONLY"
-SHUTTER_STATE_ID = "ENABLED"
-SOFTWARE_VERSION_ID = "ISS 11.00 05-24-2006"
-SPACECRAFT_CLOCK_CNT_PARTITION = 1
-SPACECRAFT_CLOCK_START_COUNT = "1702360365.220"
-SPACECRAFT_CLOCK_STOP_COUNT = "1702360370.120"
-START_TIME = 2011-346T05:02:19.773
-STOP_TIME = 2011-346T05:02:24.373
-TARGET_DESC = "Enceladus"
-TARGET_LIST = "N/A"
-TARGET_NAME = "ENCELADUS"
-TELEMETRY_FORMAT_ID = "S&ER3"
-VALID_MAXIMUM = (4095,4095)
-OBJECT = IMAGE_HEADER
-      INTERCHANGE_FORMAT = ASCII
-      HEADER_TYPE = VICAR2
-       BYTES = 3144
-      RECORDS = 1
-      ^DESCRIPTION = "../../label/vicar2.txt"
-END_OBJECT = IMAGE_HEADER
-OBJECT = TELEMETRY_TABLE
-      INTERCHANGE_FORMAT = BINARY
-      ROWS = 1
-      COLUMNS = 2
-      ROW_BYTES = 1048
-      ^STRUCTURE = "../../label/tlmtab.fmt"
-      OBJECT = COLUMN
-            NAME = NULL_PADDING
-            DATA_TYPE = MSB_UNSIGNED_INTEGER
-            START_BYTE = 61
-            BYTES = 987
-      END_OBJECT = COLUMN
-END_OBJECT = TELEMETRY_TABLE
-OBJECT = LINE_PREFIX_TABLE
-      INTERCHANGE_FORMAT = BINARY
-      ROWS = 1024
-      COLUMNS = 7
-      ROW_BYTES = 24
-      ROW_SUFFIX_BYTES = 1024
-      ^LINE_PREFIX_STRUCTURE = "../../label/prefix3.fmt"
-END_OBJECT = LINE_PREFIX_TABLE
-OBJECT = IMAGE
-      LINES = 1024
-      LINE_SAMPLES = 1024
-      SAMPLE_BITS = 8
-      SAMPLE_TYPE = SUN_INTEGER
-      LINE_PREFIX_BYTES = 24
-END_OBJECT = IMAGE
-END

knoten/vis.py

@@ -134,8 +134,6 @@ def point_info(cube_path, x, y, point_type, allow_outside=False):
 
 
 def plot_diff(data, title='diff plot', colx='x', coly='y', coldx='diffx', coldy='diffy', colmag='magnitude', width=500, height=500):
-    import matplotlib.cm as cm
-    from matplotlib.colors import Normalize
 
     fig = make_subplots(rows=2, cols=2, column_widths=[0.9, .1], row_width=[.9, .1],
                         shared_xaxes=True, shared_yaxes=True, horizontal_spacing = 0.01, vertical_spacing = 0.01)
@@ -175,7 +173,7 @@ def plot_diff(data, title='diff plot', colx='x', coly='y', coldx='diffx', coldy=
                         colorscale="viridis",
                 )), row=2, col=1)
 
-    xavg = data.groupby(colx).apply(np.mean)
+    xavg = data.groupby(colx).mean()
     fig.add_trace(go.Scatter(x=xavg.index, y=xavg[colmag],
         customdata=xavg,
         name=f'{colx} mean error',
@@ -187,7 +185,7 @@ def plot_diff(data, title='diff plot', colx='x', coly='y', coldx='diffx', coldy=
                 colorscale="viridis",
         )), row=1, col=1)
 
-    yavg = data.groupby(coly).apply(np.mean)
+    yavg = data.groupby(coly).mean()
     fig.add_trace(go.Scatter(x=yavg[colmag],y=yavg.index,
         customdata=yavg,
         name=f'{coly} mean error',

recipe/meta.yaml

@@ -35,6 +35,8 @@ requirements:
     - scipy
     - shapely
     - usgscsm>=2.0
+    - kalasiris
+    - plotly
 
 test:
   imports: