water-column-sonar-annotation 26.1.8__py3-none-any.whl

tests/astronomical/test_astronomical_manager.py

@@ -0,0 +1,148 @@
+import numpy as np
+
+from water_column_sonar_annotation.astronomical import AstronomicalManager
+
+
+#######################################################
+def setup_module():
+    print("setup")
+
+
+def teardown_module():
+    print("teardown")
+
+
+def test_get_solar_azimuth():
+    astronomical_manager = AstronomicalManager()
+    # https://www.suncalc.org/#/39.9812,-105.2495,13/2026.01.26/11:52/1/3
+    azimuth_noon = astronomical_manager.get_solar_azimuth(
+        iso_time="2026-01-26T19:00:00Z",  # noon
+        latitude=39.9674884,  # Boulder
+        longitude=-105.2532602,
+    )
+    assert np.isclose(azimuth_noon, 31.38)  # 15)
+
+    azimuth_sunset = astronomical_manager.get_solar_azimuth(
+        iso_time="2026-01-26T00:00:00Z",  # sunset
+        latitude=39.9674884,
+        longitude=-105.2532602,
+    )
+    assert np.isclose(azimuth_sunset, 1.25)  # 27)
+
+
+def test_get_solar_azimuth_boulder_2pm():
+    # 2026-01-29 @2pm is UTC: "2026-01-29T21:02:23Z"
+    astronomical_manager = AstronomicalManager()
+    # https://www.suncalc.org/#/39.9812,-105.2495,13/2026.01.26/11:52/1/3
+    azimuth_noon = astronomical_manager.get_solar_azimuth(
+        iso_time="2026-01-29T21:02:23Z",  # 2pm
+        latitude=39.9674884,  # Boulder
+        longitude=-105.2532602,
+    )
+    assert np.isclose(azimuth_noon, 27.01)
+
+
+### PHASE OF DAY ###
+# { 'dawn': 1, 'day': 2, 'dusk': 3, 'night': 4 }
+def test_phase_of_day_at_noon():
+    astronomical_manager = AstronomicalManager()
+    phase = astronomical_manager.phase_of_day(
+        iso_time="2026-01-27T19:00:00Z",  # noon
+        latitude=39.9674884,  # Boulder
+        longitude=-105.2532602,
+    )
+    assert phase == 2
+
+
+def test_phase_of_day_at_midnight():
+    astronomical_manager = AstronomicalManager()
+    phase = astronomical_manager.phase_of_day(
+        iso_time="2026-01-28T07:00:00Z",
+        latitude=39.9674884,  # Boulder
+        longitude=-105.2532602,
+    )
+    assert phase == 4
+
+
+def test_phase_of_day_before_sunset():
+    astronomical_manager = AstronomicalManager()
+    phase = astronomical_manager.phase_of_day(
+        # sunset is at 5:13pm on jan 27th, per https://psl.noaa.gov/boulder/boulder.sunset.html
+        iso_time="2026-01-28T00:09:00Z",  # sunset @ 5:13pm, nautical sunset @6:16pm
+        latitude=39.9674884,  # Boulder
+        longitude=-105.2532602,
+    )
+    assert phase == 2  # day
+
+
+def test_phase_of_day_after_sunset():
+    astronomical_manager = AstronomicalManager()
+    phase = astronomical_manager.phase_of_day(
+        # sunset is at 5:13pm on jan 27th, per https://psl.noaa.gov/boulder/boulder.sunset.html
+        iso_time="2026-01-28T00:10:00Z",  # sunset @ 5:13pm, nautical sunset @6:16pm
+        latitude=39.9674884,  # Boulder
+        longitude=-105.2532602,
+    )
+    assert phase == 3  # dusk
+
+
+def test_phase_of_day_before_nautical_sunset():
+    astronomical_manager = AstronomicalManager()
+    # an hour'ish later is nautical sunset
+    phase_before_nautical_sunset = astronomical_manager.phase_of_day(
+        iso_time="2026-01-28T01:15:00Z",  # sunset @5:13pm, nautical sunset @6:16pm
+        latitude=39.9674884,  # Boulder
+        longitude=-105.2532602,
+    )
+    assert phase_before_nautical_sunset == 3  # dusk
+
+
+def test_phase_of_day_after_nautical_sunset():
+    astronomical_manager = AstronomicalManager()
+    phase_after_nautical_sunset = astronomical_manager.phase_of_day(
+        iso_time="2026-01-28T01:17:00Z",  # sunset @5:13pm, nautical sunset @6:16pm
+        latitude=39.9674884,  # Boulder
+        longitude=-105.2532602,
+    )
+    assert phase_after_nautical_sunset == 4  # night
+
+
+def test_phase_of_day_before_sunrise():
+    astronomical_manager = AstronomicalManager()
+    phase_at_sunrise = astronomical_manager.phase_of_day(
+        iso_time="2026-01-27T14:13:00Z",  # sunrise @7:13am, nautical sunrise @6:12am
+        latitude=39.9674884,  # Boulder
+        longitude=-105.2532602,
+    )
+    assert phase_at_sunrise == 1  # dusk
+
+
+def test_phase_of_day_after_sunrise():
+    astronomical_manager = AstronomicalManager()
+    phase_at_sunrise = astronomical_manager.phase_of_day(
+        iso_time="2026-01-27T14:20:00Z",  # sunrise @7:13am, nautical sunrise @6:12am
+        latitude=39.9674884,  # Boulder
+        longitude=-105.2532602,
+    )
+    assert phase_at_sunrise == 2  # day
+
+
+def test_phase_of_day_before_nautical_sunrise():
+    astronomical_manager = AstronomicalManager()
+    # about an hour before is nautical sunrise
+    phase_before_nautical_sunrise = astronomical_manager.phase_of_day(
+        iso_time="2026-01-27T13:12:00Z",  # sunrise @7:13am, nautical sunrise @6:12am
+        latitude=39.9674884,  # Boulder
+        longitude=-105.2532602,
+    )
+    assert phase_before_nautical_sunrise == 4  # night
+
+
+def test_phase_of_day_after_nautical_sunrise():
+    astronomical_manager = AstronomicalManager()
+    phase = astronomical_manager.phase_of_day(
+        iso_time="2026-01-27T14:13:00Z",  # sunrise @7:13am, nautical sunrise @6:12am
+        latitude=39.9674884,  # Boulder
+        longitude=-105.2532602,
+    )
+    assert phase == 1

tests/conftest.py

@@ -0,0 +1,60 @@
+from pathlib import Path
+
+import pooch
+import pytest
+
+HERE = Path(__file__).parent.absolute()
+TEST_DATA_FOLDER = HERE / "test_resources"
+
+HB1906_DATA = pooch.create(
+    path=pooch.os_cache("water-column-sonar-annotation"),
+    base_url="https://github.com/CI-CMG/water-column-sonar-annotation/releases/download/v26.1.0/",
+    retry_if_failed=1,
+    registry={
+        "HB201906_BOTTOMS.zip": "sha256:20609581493ea3326c1084b6868e02aafbb6c0eae871d946f30b8b5f0e7ba059",
+        "HB201906_EVR.zip": "sha256:ceed912a25301be8f1b8f91e134d0ca4cff717f52b6623a58677832fd60c2990",
+        #
+        # "ne_50m_coastline.shp": "sha256:797d675af9613f80b51ab6049fa32e589974d7a97c6497ca56772965f179ed26",
+        # "ne_50m_coastline.shx": "sha256:0ff1792f2d16b58246d074215edd9d12fa280880ecaad61a91b9382fee854065",
+        #
+        "ne_10m_coastline.shp": "sha256:459a4a97c09db19aadf5244026612de9d43748be27f83a360242b99f7fabb3c1",
+        "ne_10m_coastline.shx": "sha256:f873afee7f56779ce52253f740ec251c2f12244aea911dc40f0a85d75de8d5f2",
+    },
+)
+
+
+def fetch_raw_files():
+    HB1906_DATA.fetch(fname="HB201906_BOTTOMS.zip", progressbar=True)
+    HB1906_DATA.fetch(fname="HB201906_EVR.zip", progressbar=True)
+
+    # HB1906_DATA.fetch(fname="ne_50m_coastline.shp", progressbar=True)
+    # HB1906_DATA.fetch(fname="ne_50m_coastline.shx", progressbar=True)
+
+    HB1906_DATA.fetch(fname="ne_10m_coastline.shp", progressbar=True)
+    HB1906_DATA.fetch(fname="ne_10m_coastline.shx", progressbar=True)
+
+    file_name = HB1906_DATA.fetch(fname="HB201906_EVR.zip", progressbar=True)
+
+    """
+    water-column-sonar-annotation user$ ls /Users/user/Library/Caches/water-column-sonar-annotation
+    HB201906_BOTTOMS.zip	HB201906_EVR.zip	ne_10m_coastline.shp	ne_10m_coastline.shx
+    """
+    return Path(file_name).parent
+
+
+@pytest.fixture(scope="session")
+def test_path():
+    return {
+        "DATA_TEST_PATH": fetch_raw_files(),
+    }
+
+
+# """
+# Folder locations in mac and windows:
+#
+# Windows
+# C:\Users\<user>\AppData\Local\echopype\Cache\2024.12.23.10.10
+#
+# MacOS
+# /Users//Library/Caches/echopype/2024.12.23.10.10
+# """

tests/cruise/test_cruise_manager.py

@@ -0,0 +1,80 @@
+import numpy as np
+import pytest
+
+from water_column_sonar_annotation.cruise import CruiseManager
+
+
+#######################################################
+def setup_module():
+    print("setup")
+
+
+def teardown_module():
+    print("teardown")
+
+
+@pytest.fixture
+def process_cruise_path(test_path):
+    return test_path["DATA_TEST_PATH"]
+
+
+#######################################################
+# def test_get_cruise_bottom_nan(process_cruise_path, tmp_path):
+#     cruise_manager = CruiseManager()
+#     cruise = cruise_manager.get_cruise()
+#     # count of non-nan values np.count_nonzero(~np.isnan(cruise.bottom.values)) / cruise.Sv.shape[1]
+#     assert len(cruise.Sv.shape) == 3
+
+
+def test_get_cruise(process_cruise_path, tmp_path):
+    cruise_manager = CruiseManager()
+    cruise = cruise_manager.get_cruise()
+    assert len(cruise.Sv.shape) == 3
+
+
+def test_get_coordinates():
+    """This only gets the depth over the interval, need to calculate the 'altitude'"""
+    cruise_manager = CruiseManager()
+    lat, lon = cruise_manager.get_coordinates(
+        start_time="2019-10-16T16:20:00",
+        end_time="2019-10-16T16:30:00",
+    )
+    assert np.isclose(lat, 41.48177)
+    assert np.isclose(lon, -68.50478)
+
+
+def test_get_depth():
+    """This only gets the depth over the interval, need to calculate the 'altitude'"""
+    cruise_manager = CruiseManager()
+    depth_value = cruise_manager.get_depth(
+        start_time="2019-10-16T16:20:00",
+        end_time="2019-10-16T16:50:00",
+    )
+    assert np.isclose(depth_value, 96.36)  # 96.356674
+
+
+def test_get_altitude():
+    """This gets the distance from EVR to the bottom"""
+    cruise_manager = CruiseManager()
+    altitude_value = cruise_manager.get_altitude(
+        start_time="2019-10-16T16:20:00",
+        end_time="2019-10-16T16:50:00",
+        bbox_max=80.0,
+    )
+    # bottom is at 96.356674
+    # setting the bbox at 80.
+    assert np.isclose(altitude_value, 16.36)  # 96.356674
+
+
+def test_get_altitude_nan_bottom():
+    """This gets the distance from EVR to the bottom when there are nan values"""
+    cruise_manager = CruiseManager()
+    altitude_value = cruise_manager.get_altitude(
+        start_time="2019-09-26T10:02:39.03",
+        end_time="2019-09-26T10:02:40.00",
+        bbox_max=10.0,
+    )
+    assert np.isclose(altitude_value, 0.0)
+
+
+# get_gps

tests/geospatial/test_geospatial_manager.py

@@ -0,0 +1,86 @@
+import numpy as np
+import pytest
+
+from water_column_sonar_annotation.geospatial import GeospatialManager
+
+
+#######################################################
+def setup_module():
+    print("setup")
+
+
+def teardown_module():
+    print("teardown")
+
+
+@pytest.fixture
+def process_check_distance_from_coastline(test_path):
+    return test_path["DATA_TEST_PATH"]
+
+
+#######################################################
+def test_check_distance_from_coastline(process_check_distance_from_coastline, tmp_path):
+    geospatial_manager = GeospatialManager()
+    # Point in middle of atlantic https://wktmap.com/?ab28cbae
+    distance = geospatial_manager.check_distance_from_coastline(
+        latitude=51.508742,
+        longitude=-30.410156,
+        shapefile_path=process_check_distance_from_coastline,
+    )
+    # assert np.isclose(distance, 1_236_212.37356)  # 1,200 km
+    assert np.isclose(distance, 1_233_910.720702243)
+
+
+def test_check_distance_from_coastline_woods_hole(
+    process_check_distance_from_coastline, tmp_path
+):
+    geospatial_manager = GeospatialManager()
+    # Point in middle of woods hole vineyard sound: https://wktmap.com/?9b405aa9
+    distance = geospatial_manager.check_distance_from_coastline(
+        latitude=41.494692,
+        longitude=-70.647926,
+        shapefile_path=process_check_distance_from_coastline,
+    )
+    # The sound is 5 km across
+    # assert np.isclose(distance, 4_457.0347)  # 4.5 km --> should be 2.5 km?
+    assert np.isclose(distance, 3_093)  # 3_093.3015 km is close enough
+
+
+def test_get_local_time():
+    geospatial_manager = GeospatialManager()
+    local_time = geospatial_manager.get_local_time(
+        iso_time="2026-01-26T20:35:00Z",
+        latitude=51.508742,
+        longitude=-30.410156,
+    )
+    assert local_time == "2026-01-26T18:35:00-02:00"
+
+
+def test_get_local_hour_of_day():
+    geospatial_manager = GeospatialManager()
+    local_hour_of_day = geospatial_manager.get_local_hour_of_day(
+        iso_time="2026-01-26T20:35:00Z",
+        latitude=51.508742,
+        longitude=-30.410156,
+    )
+    assert local_hour_of_day == 18
+
+
+def test_get_month_of_year_january():
+    geospatial_manager = GeospatialManager()
+    month_of_year = geospatial_manager.get_month_of_year(
+        iso_time="2026-01-26T20:35:00Z",
+        latitude=51.508742,
+        longitude=-30.410156,
+    )
+    assert month_of_year == 1
+
+
+def test_get_month_of_year_july():
+    geospatial_manager = GeospatialManager()
+    month_of_year = geospatial_manager.get_month_of_year(
+        iso_time="2026-07-13T12:00:00Z",
+        latitude=51.508742,
+        longitude=-30.410156,
+    )
+    assert month_of_year == 7

tests/record/test_echoview_record_manager.py

@@ -0,0 +1,160 @@
+import pytest
+
+from water_column_sonar_annotation.record import EchoviewRecordManager
+
+
+#######################################################
+def setup_module():
+    print("setup")
+
+
+def teardown_module():
+    print("teardown")
+
+
+@pytest.fixture
+def process_evr_file_path(test_path):
+    return test_path["DATA_TEST_PATH"]
+
+
+#######################################################
+# @pytest.mark.skip(reason="todo implement this")
+# def test_process_evr_file(process_evr_file_path, tmp_path):
+#     echoview_record_manager = EchoviewRecordManager()
+#     echoview_record_manager.process_point_data()
+#     assert True
+
+record_header = """EVRG 7 11.0.244.39215
+11"""
+
+possible_herring_example = """13 4 7 0 3 -1 1 20190925 2247242130  24.1290795746 20190925 2247362460  35.1668500183
+0
+0
+possible_herring
+20190925 2247242130 24.1290795746 20190925 2247242130 35.1668500183 20190925 2247362460 35.1668500183 20190925 2247362460 24.1290795746 1
+100"""
+
+fish_school_example = """13 30 8 0 7 -1 1 20190925 1749451605  20.2800268439 20190925 1749501645  26.3028953087
+0
+10
+School detected with:
+Minimum data threshold:  -66.00
+Maximum data threshold: (none)
+Distance mode: GPS distance
+Minimum total school height (meters):   4.00
+Minimum candidate length (meters):   1.00
+Minimum candidate height (meters):   2.00
+Maximum vertical linking distance (meters):   2.00
+Maximum horizontal linking distance (meters):  20.00
+Minimum total school length (meters):   4.00
+fish_school
+20190925 1749451605 22.0286015595 20190925 1749451605 25.1371788317 20190925 1749461355 25.1371788317 20190925 1749461605 25.1420359837 20190925 1749461605 25.7200370702 20190925 1749471370 25.7200370702 20190925 1749471620 25.7248942222 20190925 1749471620 26.1086092292 20190925 1749481385 26.1086092292 20190925 1749481635 26.1134663812 20190925 1749481635 26.3028953087 20190925 1749491640 26.3028953087 20190925 1749491640 25.1371788317 20190925 1749501645 25.1371788317 20190925 1749501645 24.5543205932 20190925 1749491640 24.5543205932 20190925 1749491640 24.3600345136 20190925 1749501645 24.3600345136 20190925 1749501645 21.6400294005 20190925 1749491640 21.6400294005 20190925 1749491640 21.0571711620 20190925 1749481635 21.0571711620 20190925 1749481635 20.2800268439 20190925 1749471620 20.2800268439 20190925 1749471620 20.8580279305 20190925 1749471370 20.8628850825 20190925 1749461605 20.8628850825 20190925 1749461605 21.8294583280 20190925 1749461355 21.8343154800 20190925 1749451605 21.8343154800 1
+Region 8"""
+
+unclassified_regions_example = """13 12 1 0 2 -1 1 20190925 2053458953  9.2818 20190925 2054119318  11.5333
+0
+0
+Unclassified regions
+20190925 2053458953 9.6034489515 20190925 2053521545 11.1197829964 20190925 2054046730 11.5333286451 20190925 2054064248 11.5333286451 20190925 2054079263 11.4414296120 20190925 2054116810 10.8440858974 20190925 2054119318 10.4764897652 20190925 2054116810 9.6953479845 20190925 2054111800 9.4196508854 20190925 2054091775 9.2818023359 20190925 2054076760 9.2818023359 20190925 2053483995 9.5574994350 0
+100"""
+
+krill_schools_example = """13 23 49 0 7 -1 1 20191020 0536541560  13.0000000000 20191020 0536591655  17.0000000000
+0
+10
+School detected with:
+Minimum data threshold:  -80.00
+Maximum data threshold: (none)
+Distance mode: GPS distance
+Minimum total school height (meters):   4.00
+Minimum candidate length (meters):   1.00
+Minimum candidate height (meters):   2.00
+Maximum vertical linking distance (meters):   2.00
+Maximum horizontal linking distance (meters):  20.00
+Minimum total school length (meters):   4.00
+krill_schools
+20191020 0536551605 15.0250000000 20191020 0536551605 17.0000000000 20191020 0536561605 17.0000000000 20191020 0536561605 16.0000000000 20191020 0536571635 16.0000000000 20191020 0536571635 15.0000000000 20191020 0536581405 15.0000000000 20191020 0536581655 15.0250000000 20191020 0536581655 16.0000000000 20191020 0536591655 16.0000000000 20191020 0536591655 15.0000000000 20191020 0536581655 15.0000000000 20191020 0536581655 14.0000000000 20191020 0536571635 14.0000000000 20191020 0536571635 13.0000000000 20191020 0536561605 13.0000000000 20191020 0536561605 14.9750000000 20191020 0536561355 15.0000000000 20191020 0536551605 15.0000000000 20191020 0536551605 14.0000000000 20191020 0536541560 14.0000000000 20191020 0536541560 15.0000000000 20191020 0536551355 15.0000000000 1
+Region 49"""
+
+ah_school_example1 = """13 5 23 0 7 -1 1 20191106 1314583780  25.3008882713 20191106 1314593790  30.2941528987
+0
+10
+School detected with:
+Minimum data threshold:  -66.00
+Maximum data threshold: (none)
+Distance mode: GPS distance
+Minimum total school height (meters):   4.00
+Minimum candidate length (meters):   1.00
+Minimum candidate height (meters):   2.00
+Maximum vertical linking distance (meters):   2.00
+Maximum horizontal linking distance (meters):  20.00
+Minimum total school length (meters):   4.00
+AH_School
+20191106 1314583780 25.4929369108 20191106 1314583780 30.2941528987 20191106 1314593790 30.2941528987 20191106 1314593790 25.3008882713 20191106 1314583780 25.3008882713 1
+Region 23"""
+
+
+ah_school_example2 = """13 16 28 0 7 -1 1 20191106 1317305715  31.8305420148 20191106 1317335745  35.8635634446
+0
+10
+School detected with:
+Minimum data threshold:  -66.00
+Maximum data threshold: (none)
+Distance mode: GPS distance
+Minimum total school height (meters):   4.00
+Minimum candidate length (meters):   1.00
+Minimum candidate height (meters):   2.00
+Maximum vertical linking distance (meters):   2.00
+Maximum horizontal linking distance (meters):  20.00
+Minimum total school length (meters):   4.00
+AH_School
+20191106 1317315725 34.3319755445 20191106 1317315725 35.8635634446 20191106 1317325735 35.8635634446 20191106 1317325735 35.4794661656 20191106 1317335745 35.4794661656 20191106 1317335745 34.5192229680 20191106 1317325735 34.5192229680 20191106 1317325735 34.3271743285 20191106 1317335745 34.3271743285 20191106 1317335745 31.8305420148 20191106 1317315725 31.8305420148 20191106 1317315725 34.1303244730 20191106 1317315475 34.1351256890 20191106 1317305715 34.1351256890 20191106 1317305715 34.3271743285 20191106 1317315475 34.3271743285 1
+Region 28"""
+
+atlantic_herring_example = """13 70 2 0 2 -1 1 20191014 0011274130  26.8060139347 20191014 0017599670  76.3659699303
+0
+0
+atlantic_herring
+20191014 0011274130 32.7282728273 20191014 0011339205 48.0014668133 20191014 0011474370 54.2354235424 20191014 0011499420 54.2354235424 20191014 0011544520 54.2354235424 20191014 0012260060 57.6640997433 20191014 0012310115 57.6640997433 20191014 0012350165 57.6640997433 20191014 0012390210 57.6640997433 20191014 0013370930 70.7554088742 20191014 0013420980 71.6905023836 20191014 0013506080 71.6905023836 20191014 0013541115 71.3788045471 20191014 0013586170 71.3788045471 20191014 0014011185 70.7554088742 20191014 0015067420 75.7425742574 20191014 0015107470 76.0542720939 20191014 0015187560 76.3659699303 20191014 0015282670 76.3659699303 20191014 0015307700 76.0542720939 20191014 0015337730 75.7425742574 20191014 0016278420 75.1191785845 20191014 0016383575 75.1191785845 20191014 0016473695 75.1191785845 20191014 0016528780 74.4957829116 20191014 0016568850 73.2489915658 20191014 0017129040 68.2618261826 20191014 0017154065 66.7033370004 20191014 0017184115 64.8331499817 20191014 0017199150 62.9629629630 20191014 0017339330 56.4173083975 20191014 0017384400 55.7939127246 20191014 0017439460 55.1705170517 20191014 0017479510 53.6120278695 20191014 0017599670 49.5599559956 20191014 0017599670 48.0014668133 20191014 0017584660 44.8844884488 20191014 0017559630 44.5727906124 20191014 0017309300 44.5727906124 20191014 0017284270 44.2610927759 20191014 0017259240 44.2610927759 20191014 0017214180 43.3259992666 20191014 0017179110 43.0143014301 20191014 0016323505 43.0143014301 20191014 0016288450 43.0143014301 20191014 0016223340 42.0792079208 20191014 0016188310 42.0792079208 20191014 0016143230 41.1441144114 20191014 0015467870 39.2739273927 20191014 0015402790 38.0271360469 20191014 0015327710 36.7803447011 20191014 0014351845 34.9101576824 20191014 0014306720 34.5984598460 20191014 0014261660 33.9750641731 20191014 0014221575 33.3516685002 20191014 0014191530 33.3516685002 20191014 0014161465 33.3516685002 20191014 0014116410 33.3516685002 20191014 0013150670 30.2346901357 20191014 0013060570 28.6762009534 20191014 0013020530 27.4294096076 20191014 0012555450 26.8060139347 20191014 0012485320 26.8060139347 20191014 0012455280 26.8060139347 20191014 0012134870 29.6112944628 20191014 0012074750 29.9229922992 20191014 0011534490 31.4814814815 20191014 0011429315 32.4165749908 20191014 0011339205 32.7282728273 20191014 0011314180 32.7282728273 1
+100"""
+
+
+def test_process_evr_record_possible_herring():
+    echoview_record_manager = EchoviewRecordManager()
+    echoview_record_manager.process_evr_record(
+        evr_record=possible_herring_example, filename="test.evr"
+    )
+
+
+def test_process_evr_record_fish_school():
+    echoview_record_manager = EchoviewRecordManager()
+    echoview_record_manager.process_evr_record(
+        evr_record=fish_school_example, filename="test.evr"
+    )
+
+
+def test_process_evr_record_unclassified_regions():
+    echoview_record_manager = EchoviewRecordManager()
+    echoview_record_manager.process_evr_record(
+        evr_record=unclassified_regions_example, filename="test.evr"
+    )
+
+
+def test_process_evr_record_krill_schools():
+    echoview_record_manager = EchoviewRecordManager()
+    echoview_record_manager.process_evr_record(
+        evr_record=krill_schools_example, filename="test.evr"
+    )
+
+
+def test_process_evr_record_ah_school_1():
+    echoview_record_manager = EchoviewRecordManager()
+    echoview_record_manager.process_evr_record(
+        evr_record=ah_school_example1, filename="test.evr"
+    )
+
+
+def test_process_evr_record_ah_school_2():
+    echoview_record_manager = EchoviewRecordManager()
+    echoview_record_manager.process_evr_record(
+        evr_record=ah_school_example2, filename="test.evr"
+    )

water_column_sonar_annotation/__init__.py

@@ -0,0 +1,5 @@
+from __future__ import absolute_import
+
+from . import astronomical, geospatial, record, shape
+
+__all__ = ["astronomical", "record", "geospatial", "shape"]

water_column_sonar_annotation/astronomical/__init__.py

@@ -0,0 +1,5 @@
+from .astronomical_manager import AstronomicalManager
+
+__all__ = [
+    AstronomicalManager,
+]

water_column_sonar_annotation/astronomical/astronomical_manager.py

@@ -0,0 +1,82 @@
+import numpy as np
+import pandas as pd
+import pvlib
+
+from water_column_sonar_annotation.geospatial import GeospatialManager
+
+
+class AstronomicalManager:
+    #######################################################
+    def __init__(
+        self,
+    ):
+        self.DECIMAL_PRECISION = 6
+        # https://github.com/CI-CMG/water-column-sonar-annotation/issues/6
+        self.SUNRISE_DEGREES = 0.0
+        self.CIVIL_TWILIGHT_DEGREES = -6.0
+        self.NAUTICAL_TWILIGHT_DEGREES = -12.0  # Requested metric to calculate
+        self.ASTRONOMICAL_TWILIGHT_DEGREES = -18.0
+
+    @staticmethod
+    def get_solar_azimuth(
+        iso_time: str = "2026-01-26T00:06:00Z",
+        latitude: float = 39.9674884,  # boulder gps coordinates
+        longitude: float = -105.2532602,
+    ) -> float:
+        """
+        Good reference for calculating: https://www.suncalc.org/#/39.9812,-105.2495,13/2026.01.26/11:52/1/3
+        utc time now: '2026-01-25T18:42:00Z' # 11:43am
+            7:14 am↑ (sunrise)
+                (Timestamp('2026-01-25 18:42:00+0000', tz='UTC'), Timestamp('2026-01-25 14:15:07.145030400+0000', tz='UTC'))
+            5:13 pm↑ (sunset)
+                (Timestamp('2026-01-25 18:42:00+0000', tz='UTC'), Timestamp('2026-01-26 00:10:51.244243200+0000', tz='UTC'))
+            solar altitude should be: 31.26°, azimuth should be: 174.01°
+        """
+        solar_position = pvlib.solarposition.get_solarposition(
+            time=pd.DatetimeIndex([iso_time]),
+            latitude=latitude,
+            longitude=longitude,
+        )
+        # 'elevation' is analogous to 'altitude' in suncalc
+        elevation = solar_position.elevation.iloc[0]
+        ### The altitude aka elevation is the angle between horizon and the center of the sun including refraction ###
+        return np.round(elevation, 2)
+
+    def phase_of_day(
+        self,
+        iso_time: str,
+        latitude: float,
+        longitude: float,
+    ) -> int:
+        """
+        Returns whether the time/gps references a Nautical Daylight time
+        Going to need to verify the az is correctly computed
+        { 'night': 4, 'dawn': 1, 'day': 2, 'dusk': 3 }
+        """
+        # categories = {"night": 4, "dawn": 1, "day": 2, "dusk": 3}
+        solar_azimuth = self.get_solar_azimuth(iso_time, latitude, longitude)
+        geospatial_manager = GeospatialManager()
+        local_hour = geospatial_manager.get_local_hour_of_day(
+            iso_time=iso_time,
+            latitude=latitude,
+            longitude=longitude,
+        )
+        if solar_azimuth < self.NAUTICAL_TWILIGHT_DEGREES:
+            return 4  # night
+        if solar_azimuth >= 0.0:
+            return 2  # day
+        if local_hour < 12:
+            return 1  # dawn
+        return 3  # dusk
+
+    # def get_moon_phase(self):
+    #     # TODO: add method for getting the moon phase
+    #     pass
+
+    # TODO: calculate moonrise and moonset
+
+
+# if __name__ == "__main__":
+#     astronomical_manager = AstronomicalManager()
+#     azimuth = astronomical_manager.get_solar_azimuth()
+#     print(azimuth)

water_column_sonar_annotation/cruise/__init__.py

@@ -0,0 +1,5 @@
+from .cruise_manager import CruiseManager
+
+__all__ = [
+    CruiseManager,
+]