\n",
+ "
"
+ ],
+ "text/plain": [
+ "<xarray.Dataset> Size: 309MB\n", + "Dimensions: (channel: 1, ping_time: 2667,\n", + " range_sample: 2894, filenames: 1)\n", + "Coordinates:\n", + " * channel (channel) <U25 100B 'WBT 741862-15 ES38-7_ES'\n", + " * ping_time (ping_time) datetime64[ns] 21kB 2023-06-26...\n", + " * range_sample (range_sample) int32 12kB 0 1 2 ... 2892 2893\n", + " * filenames (filenames) int32 4B 0\n", + "Data variables: (12/28)\n", + " Sv (channel, ping_time, range_sample) float64 62MB ...\n", + " echo_range (channel, ping_time, range_sample) float64 62MB ...\n", + " tau_effective (channel) float64 8B 0.000832\n", + " frequency_nominal (channel) float64 8B 3.8e+04\n", + " sound_speed float64 8B 1.44e+03\n", + " sound_absorption (channel) float64 8B 0.003039\n", + " ... ...\n", + " receiver_sampling_frequency (channel) float64 8B 1.5e+06\n", + " source_filenames (filenames) <U80 320B '..\\\\..\\\\..\\\\run_ech...\n", + " water_level float64 8B 0.0\n", + " depth (channel, ping_time, range_sample) float64 62MB ...\n", + " angle_alongship (channel, ping_time, range_sample) float64 62MB ...\n", + " angle_athwartship (channel, ping_time, range_sample) float64 62MB ...\n", + "Attributes:\n", + " processing_software_name: echopype\n", + " processing_software_version: 0.11.1.dev26+gebe1d44.d20251123\n", + " processing_time: 2026-02-24T08:37:06+00:00\n", + " processing_function: calibrate.compute_Sv
<xarray.DataArray 'bottom_depth' (ping_time: 2667)> Size: 21kB\n", + "array([31.8731256, 31.8731256, 31.8731256, ..., 37.0913232, 37.0913232,\n", + " 37.0913232])\n", + "Coordinates:\n", + " * ping_time (ping_time) datetime64[ns] 21kB 2023-06-26T23:58:35.018728 ......\n", + "Attributes:\n", + " detector: blackwell\n", + " threshold_Sv: -40.0\n", + " threshold_angle_major: 0.4\n", + " threshold_angle_minor: 0.4\n", + " offset_m: 0.3\n", + " channel: WBT 741862-15 ES38-7_ES" + ], + "text/plain": [ + "
\n",
+ "
"
+ ],
+ "text/plain": [
+ "<xarray.Dataset> Size: 154MB\n", + "Dimensions: (channel: 1, ping_time: 2667,\n", + " range_sample: 1447, filenames: 1)\n", + "Coordinates:\n", + " * channel (channel) <U25 100B 'WBT 741862-15 ES38-7_ES'\n", + " * ping_time (ping_time) datetime64[ns] 21kB 2023-06-26...\n", + " * range_sample (range_sample) int32 6kB 0 1 2 ... 1445 1446\n", + " * filenames (filenames) int32 4B 0\n", + "Data variables: (12/28)\n", + " Sv (channel, ping_time, range_sample) float64 31MB ...\n", + " echo_range (channel, ping_time, range_sample) float64 31MB ...\n", + " tau_effective (channel) float64 8B 0.000832\n", + " frequency_nominal (channel) float64 8B 3.8e+04\n", + " sound_speed float64 8B 1.44e+03\n", + " sound_absorption (channel) float64 8B 0.003039\n", + " ... ...\n", + " receiver_sampling_frequency (channel) float64 8B 1.5e+06\n", + " source_filenames (filenames) <U80 320B '..\\\\..\\\\..\\\\run_ech...\n", + " water_level float64 8B 0.0\n", + " depth (channel, ping_time, range_sample) float64 31MB ...\n", + " angle_alongship (channel, ping_time, range_sample) float64 31MB ...\n", + " angle_athwartship (channel, ping_time, range_sample) float64 31MB ...\n", + "Attributes:\n", + " processing_software_name: echopype\n", + " processing_software_version: 0.11.1.dev26+gebe1d44.d20251123\n", + " processing_time: 2026-02-24T08:37:06+00:00\n", + " processing_function: calibrate.compute_Sv
\n",
+ "
"
+ ],
+ "text/plain": [
+ "<xarray.Dataset> Size: 1GB\n", + "Dimensions: (channel: 6, ping_time: 2667,\n", + " range_sample: 2894, filenames: 1)\n", + "Coordinates:\n", + " * channel (channel) <U25 600B 'WBT 741110-15 ES18_ES...\n", + " * ping_time (ping_time) datetime64[ns] 21kB 2023-06-26...\n", + " * range_sample (range_sample) int32 12kB 0 1 2 ... 2892 2893\n", + " * filenames (filenames) int32 4B 0\n", + "Data variables: (12/26)\n", + " TS (channel, ping_time, range_sample) float64 370MB ...\n", + " echo_range (channel, ping_time, range_sample) float64 370MB ...\n", + " frequency_nominal (channel) float64 48B 1.8e+04 ... 2e+05\n", + " sound_speed float64 8B 1.44e+03\n", + " sound_absorption (channel) float64 48B 0.0008724 ... 0.02312\n", + " temperature float64 8B 5.4\n", + " ... ...\n", + " impedance_transceiver (channel) int32 24B 5400 5400 ... 5400 5400\n", + " receiver_sampling_frequency (channel) float64 48B 1.5e+06 ... 1.5e+06\n", + " source_filenames (filenames) <U80 320B '..\\\\..\\\\..\\\\run_ech...\n", + " water_level float64 8B 0.0\n", + " angle_alongship (channel, ping_time, range_sample) float64 370MB ...\n", + " angle_athwartship (channel, ping_time, range_sample) float64 370MB ...\n", + "Attributes:\n", + " processing_software_name: echopype\n", + " processing_software_version: 0.11.1.dev26+gebe1d44.d20251123\n", + " processing_time: 2026-02-24T08:41:32+00:00\n", + " processing_function: calibrate.compute_TS
<xarray.Dataset> Size: 247MB\n", + "Dimensions: (channel: 1, ping_time: 2667,\n", + " range_sample: 2894, filenames: 1)\n", + "Coordinates:\n", + " * channel (channel) <U25 100B 'WBT 741862-15 ES38-7_ES'\n", + " * ping_time (ping_time) datetime64[ns] 21kB 2023-06-26...\n", + " * range_sample (range_sample) int32 12kB 0 1 2 ... 2892 2893\n", + " * filenames (filenames) int32 4B 0\n", + " depth (range_sample) float64 23kB 0.0 ... nan\n", + "Data variables: (12/26)\n", + " TS (channel, ping_time, range_sample) float64 62MB ...\n", + " echo_range (channel, ping_time, range_sample) float64 62MB ...\n", + " frequency_nominal (channel) float64 8B 3.8e+04\n", + " sound_speed float64 8B 1.44e+03\n", + " sound_absorption (channel) float64 8B 0.003039\n", + " temperature float64 8B 5.4\n", + " ... ...\n", + " impedance_transceiver (channel) int32 4B 5400\n", + " receiver_sampling_frequency (channel) float64 8B 1.5e+06\n", + " source_filenames (filenames) <U80 320B '..\\\\..\\\\..\\\\run_ech...\n", + " water_level float64 8B 0.0\n", + " angle_alongship (channel, ping_time, range_sample) float64 62MB ...\n", + " angle_athwartship (channel, ping_time, range_sample) float64 62MB ...\n", + "Attributes:\n", + " processing_software_name: echopype\n", + " processing_software_version: 0.11.1.dev26+gebe1d44.d20251123\n", + " processing_time: 2026-02-24T08:41:32+00:00\n", + " processing_function: calibrate.compute_TS" + ], + "text/plain": [ + "
\n",
+ "
"
+ ],
+ "text/plain": [
+ "<xarray.Dataset> Size: 247MB\n", + "Dimensions: (channel: 1, ping_time: 2667,\n", + " range_sample: 2894, filenames: 1)\n", + "Coordinates:\n", + " * channel (channel) object 8B 'WBT 741862-15 ES38-7_ES'\n", + " * ping_time (ping_time) datetime64[ns] 21kB 2023-06-26...\n", + " * range_sample (range_sample) int32 12kB 0 1 2 ... 2892 2893\n", + " * filenames (filenames) int32 4B 0\n", + " depth (range_sample) float64 23kB 0.0 ... nan\n", + "Data variables: (12/26)\n", + " TS (channel, ping_time, range_sample) float64 62MB ...\n", + " echo_range (channel, ping_time, range_sample) float64 62MB ...\n", + " frequency_nominal (channel) float64 8B 3.8e+04\n", + " sound_speed float64 8B 1.44e+03\n", + " sound_absorption (channel) float64 8B 0.003039\n", + " temperature float64 8B 5.4\n", + " ... ...\n", + " impedance_transceiver (channel) int32 4B 5400\n", + " receiver_sampling_frequency (channel) float64 8B 1.5e+06\n", + " source_filenames (filenames) <U80 320B '..\\\\..\\\\..\\\\run_ech...\n", + " water_level float64 8B 0.0\n", + " angle_alongship (channel, ping_time, range_sample) float64 62MB ...\n", + " angle_athwartship (channel, ping_time, range_sample) float64 62MB ...\n", + "Attributes:\n", + " processing_software_name: echopype\n", + " processing_software_version: 0.11.1.dev26+gebe1d44.d20251123\n", + " processing_time: 2026-02-24T08:41:32+00:00\n", + " processing_function: calibrate.compute_TS\n", + " mask_software_name: echopype\n", + " mask_software_version: 0.11.1.dev26+gebe1d44.d20251123\n", + " mask_time: 2026-02-24T08:42:00+00:00\n", + " mask_function: mask.apply_mask
\n",
+ "
"
+ ],
+ "text/plain": [
+ "<xarray.Dataset> Size: 247MB\n", + "Dimensions: (channel: 1, ping_time: 2667,\n", + " range_sample: 2894, filenames: 1)\n", + "Coordinates:\n", + " * channel (channel) object 8B 'WBT 741862-15 ES38-7_ES'\n", + " * ping_time (ping_time) datetime64[ns] 21kB 2023-06-26...\n", + " * range_sample (range_sample) int32 12kB 0 1 2 ... 2892 2893\n", + " * filenames (filenames) int32 4B 0\n", + " depth (range_sample) float64 23kB 0.0 ... nan\n", + "Data variables: (12/29)\n", + " TS (channel, ping_time, range_sample) float64 62MB ...\n", + " echo_range (channel, ping_time, range_sample) float64 62MB ...\n", + " frequency_nominal (channel) float64 8B 3.8e+04\n", + " sound_speed float64 8B 1.44e+03\n", + " sound_absorption (channel) float64 8B 0.003039\n", + " temperature float64 8B 5.4\n", + " ... ...\n", + " water_level float64 8B 0.0\n", + " angle_alongship (channel, ping_time, range_sample) float64 62MB ...\n", + " angle_athwartship (channel, ping_time, range_sample) float64 62MB ...\n", + " tau_effective (channel) float64 8B 0.000832\n", + " sample_interval (channel, ping_time) float64 21kB 4.8e-05 ...\n", + " beam_type (channel) int32 4B 65\n", + "Attributes:\n", + " processing_software_name: echopype\n", + " processing_software_version: 0.11.1.dev26+gebe1d44.d20251123\n", + " processing_time: 2026-02-24T08:41:32+00:00\n", + " processing_function: calibrate.compute_TS\n", + " mask_software_name: echopype\n", + " mask_software_version: 0.11.1.dev26+gebe1d44.d20251123\n", + " mask_time: 2026-02-24T08:42:00+00:00\n", + " mask_function: mask.apply_mask
\n",
+ " \n",
+ " 
\n",
+ "
\n",
+ " "
+ ],
+ "text/plain": [
+ "Canvas(toolbar=Toolbar(toolitems=[('Home', 'Reset original view', 'home', 'home'), ('Back', 'Back to previous …"
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "import numpy as np\n",
+ "import matplotlib.pyplot as plt\n",
+ "import matplotlib as mpl\n",
+ "from matplotlib.widgets import SpanSelector\n",
+ "import echopype.colormap\n",
+ "\n",
+ "vmin_ts, vmax_ts = -80, -20\n",
+ "cmap_ts = mpl.colormaps[\"ep.ek500\"]\n",
+ "ang_cmap = \"RdBu_r\"\n",
+ "ymax = 35\n",
+ "\n",
+ "ds = ds_m2.sortby(\"ping_time\")\n",
+ "\n",
+ "depth1d = ds[\"depth\"].values.astype(float)\n",
+ "valid = np.isfinite(depth1d)\n",
+ "idx = np.where(valid)[0]\n",
+ "depth1d = depth1d[valid]\n",
+ "\n",
+ "def to_depth(da):\n",
+ " return (\n",
+ " da.isel(range_sample=idx)\n",
+ " .assign_coords(depth=(\"range_sample\", depth1d))\n",
+ " .swap_dims({\"range_sample\": \"depth\"})\n",
+ " .sortby(\"depth\")\n",
+ " )\n",
+ "\n",
+ "TS = to_depth(ds[\"TS\"])\n",
+ "Al = to_depth(ds[\"angle_alongship\"])\n",
+ "At = to_depth(ds[\"angle_athwartship\"])\n",
+ "\n",
+ "bottom = ds[\"bottom\"] if \"bottom\" in ds else None\n",
+ "\n",
+ "# detections\n",
+ "ev_t = st_ev[\"ping_time\"].values.astype(\"datetime64[ns]\")\n",
+ "ev_z = st_ev[\"target_range\"].values.astype(float) # <\n",
+ "