Frames Module / フレームモジュール¶
The wandas.frames module provides immutable Frame families for time-domain,
frequency-domain, time-frequency, cepstral, octave-band, and roughness results.
wandas.framesは時間領域、周波数領域、時間周波数、ケプストラム、オクターブ帯域、
ラフネス結果を表すimmutableなFrameファミリーを提供します。
For callable input/output contracts and a complete extension example, see Add your own processing. 独自関数の入出力契約と拡張の完全な例は、独自処理のHow-toを参照してください。
wandas.frames.channel.ChannelFrame
¶
Bases: BaseFrame[NDArrayReal], ChannelProcessingMixin, ChannelTransformMixin
Channel-based data frame for handling audio signals and time series data.
This frame represents channel-based data such as audio signals and time series data, with each channel containing data samples in the time domain.
Source code in wandas/frames/channel.py
388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 839 840 841 842 843 844 845 846 847 848 849 850 851 852 853 854 855 856 857 858 859 860 861 862 863 864 865 866 867 868 869 870 871 872 873 874 875 876 877 878 879 880 881 882 883 884 885 886 887 888 889 890 891 892 893 894 895 896 897 898 899 900 901 902 903 904 905 906 907 908 909 910 911 912 913 914 915 916 917 918 919 920 921 922 923 924 925 926 927 928 929 930 931 932 933 934 935 936 937 938 939 940 941 942 943 944 945 946 947 948 949 950 951 952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967 968 969 970 971 972 973 974 975 976 977 978 979 980 981 982 983 984 985 986 987 988 989 990 991 992 993 994 995 996 997 998 999 1000 1001 1002 1003 1004 1005 1006 1007 1008 1009 1010 1011 1012 1013 1014 1015 1016 1017 1018 1019 1020 1021 1022 1023 1024 1025 1026 1027 1028 1029 1030 1031 1032 1033 1034 1035 1036 1037 1038 1039 1040 1041 1042 1043 1044 1045 1046 1047 1048 1049 1050 1051 1052 1053 1054 1055 1056 1057 1058 1059 1060 1061 1062 1063 1064 1065 1066 1067 1068 1069 1070 1071 1072 1073 1074 1075 1076 1077 1078 1079 1080 1081 1082 1083 1084 1085 1086 1087 1088 1089 1090 1091 1092 1093 1094 1095 1096 1097 1098 1099 1100 1101 1102 1103 1104 1105 1106 1107 1108 1109 1110 1111 1112 1113 1114 1115 1116 1117 1118 1119 1120 1121 1122 1123 1124 1125 1126 1127 1128 1129 1130 1131 1132 1133 1134 1135 1136 1137 1138 1139 1140 1141 1142 1143 1144 1145 1146 1147 1148 1149 1150 1151 1152 1153 1154 1155 1156 1157 1158 1159 1160 1161 1162 1163 1164 1165 1166 1167 1168 1169 1170 1171 1172 1173 1174 1175 1176 1177 1178 1179 1180 1181 1182 1183 1184 1185 1186 1187 1188 1189 1190 1191 1192 1193 1194 1195 1196 1197 1198 1199 1200 1201 1202 1203 1204 1205 1206 1207 1208 1209 1210 1211 1212 1213 1214 1215 1216 1217 1218 1219 1220 1221 1222 1223 1224 1225 1226 1227 1228 1229 1230 1231 1232 1233 1234 1235 1236 1237 1238 1239 1240 1241 1242 1243 1244 1245 1246 1247 1248 1249 1250 1251 1252 1253 1254 1255 1256 1257 1258 1259 1260 1261 1262 1263 1264 1265 1266 1267 1268 1269 1270 1271 1272 1273 1274 1275 1276 1277 1278 1279 1280 1281 1282 1283 1284 1285 1286 1287 1288 1289 1290 1291 1292 1293 1294 1295 1296 1297 1298 1299 1300 1301 1302 1303 1304 1305 1306 1307 1308 1309 1310 1311 1312 1313 1314 1315 1316 1317 1318 1319 1320 1321 1322 1323 1324 1325 1326 1327 1328 1329 1330 1331 1332 1333 1334 1335 1336 1337 1338 1339 1340 1341 1342 1343 1344 1345 1346 1347 1348 1349 1350 1351 1352 1353 1354 1355 1356 1357 1358 1359 1360 1361 1362 1363 1364 1365 1366 1367 1368 1369 1370 1371 1372 1373 1374 1375 1376 1377 1378 1379 1380 1381 1382 1383 1384 1385 1386 1387 1388 1389 1390 1391 1392 1393 1394 1395 1396 1397 1398 1399 1400 1401 1402 1403 1404 1405 1406 1407 1408 1409 1410 1411 1412 1413 1414 1415 1416 1417 1418 1419 1420 1421 1422 1423 1424 1425 1426 1427 1428 1429 1430 1431 1432 1433 1434 1435 1436 1437 1438 1439 1440 1441 1442 1443 1444 1445 1446 1447 1448 1449 1450 1451 1452 1453 1454 1455 1456 1457 1458 1459 1460 1461 1462 1463 1464 1465 1466 1467 1468 1469 1470 1471 1472 1473 1474 1475 1476 1477 1478 1479 1480 1481 1482 1483 1484 1485 1486 1487 1488 1489 1490 1491 1492 1493 1494 1495 1496 1497 1498 1499 1500 1501 1502 1503 1504 1505 1506 1507 1508 1509 1510 1511 1512 1513 1514 1515 1516 1517 1518 1519 1520 1521 1522 1523 1524 1525 1526 1527 1528 1529 1530 1531 1532 1533 1534 1535 1536 1537 1538 1539 1540 1541 1542 1543 1544 1545 1546 1547 1548 1549 1550 1551 1552 1553 1554 1555 1556 1557 1558 1559 1560 1561 1562 1563 1564 1565 1566 1567 1568 1569 1570 1571 1572 1573 1574 1575 1576 1577 1578 1579 1580 1581 1582 1583 1584 1585 1586 1587 1588 1589 1590 1591 1592 1593 1594 1595 1596 1597 1598 1599 1600 1601 1602 1603 1604 1605 1606 1607 1608 1609 1610 1611 1612 1613 1614 1615 1616 1617 1618 1619 1620 1621 1622 1623 1624 1625 1626 1627 1628 1629 1630 1631 1632 1633 1634 1635 1636 1637 1638 1639 1640 1641 1642 1643 1644 1645 1646 1647 1648 1649 1650 1651 1652 1653 1654 1655 1656 1657 1658 1659 1660 1661 1662 1663 1664 1665 1666 1667 1668 1669 1670 1671 1672 1673 1674 1675 1676 1677 1678 1679 1680 1681 1682 1683 1684 1685 1686 1687 1688 1689 1690 1691 1692 1693 1694 1695 1696 1697 1698 1699 1700 1701 1702 1703 1704 1705 1706 1707 1708 1709 1710 1711 1712 1713 1714 1715 1716 1717 1718 1719 1720 1721 1722 1723 1724 1725 1726 1727 1728 1729 1730 1731 1732 1733 1734 1735 1736 1737 1738 1739 1740 1741 1742 1743 1744 1745 1746 1747 1748 1749 1750 1751 1752 1753 1754 1755 1756 1757 1758 1759 1760 1761 1762 1763 1764 1765 1766 1767 1768 1769 1770 1771 1772 1773 1774 1775 1776 1777 1778 1779 1780 1781 1782 1783 1784 1785 1786 1787 1788 1789 1790 1791 1792 1793 1794 1795 1796 1797 1798 1799 1800 1801 1802 1803 1804 1805 1806 1807 1808 1809 1810 1811 1812 1813 1814 1815 1816 1817 1818 1819 1820 1821 1822 1823 1824 1825 1826 1827 1828 1829 1830 1831 1832 1833 1834 1835 1836 1837 1838 1839 1840 1841 1842 1843 1844 1845 1846 1847 1848 1849 1850 1851 1852 1853 1854 1855 1856 1857 1858 1859 1860 1861 1862 1863 1864 1865 1866 1867 1868 1869 1870 1871 1872 1873 1874 1875 1876 1877 1878 1879 1880 1881 1882 1883 1884 1885 1886 1887 1888 1889 1890 1891 1892 1893 1894 1895 1896 1897 1898 1899 1900 1901 1902 1903 1904 1905 1906 1907 1908 1909 1910 1911 1912 1913 1914 1915 1916 1917 1918 1919 1920 1921 1922 1923 1924 1925 1926 1927 1928 1929 1930 1931 1932 1933 1934 1935 1936 1937 1938 1939 1940 1941 1942 1943 1944 1945 1946 1947 1948 1949 1950 1951 1952 1953 1954 1955 1956 1957 1958 1959 1960 1961 1962 1963 1964 1965 1966 1967 1968 1969 1970 1971 1972 1973 1974 1975 1976 1977 1978 1979 1980 1981 1982 1983 1984 1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 2019 2020 2021 2022 2023 2024 2025 2026 2027 2028 2029 2030 2031 2032 2033 2034 2035 2036 2037 2038 2039 2040 2041 2042 2043 2044 2045 2046 2047 2048 2049 2050 2051 2052 2053 2054 2055 2056 2057 2058 2059 2060 2061 2062 2063 2064 2065 2066 2067 2068 2069 2070 2071 2072 2073 2074 2075 2076 2077 2078 2079 2080 2081 2082 2083 2084 2085 2086 2087 2088 2089 2090 2091 2092 2093 2094 2095 2096 2097 2098 2099 2100 2101 2102 2103 2104 2105 2106 2107 2108 2109 2110 2111 2112 2113 2114 2115 | |
Attributes¶
time
property
¶
Get time array for the signal.
The time array represents the start time of each sample, calculated as sample_index / sampling_rate. This provides a uniform, evenly-spaced time axis that is consistent across all frame types in wandas.
For frames resulting from windowed analysis operations (e.g., FFT, loudness, roughness), each time point corresponds to the start of the analysis window, not the center. This differs from some libraries (e.g., MoSQITo) which use window center times, but does not affect the calculated values themselves.
Returns:
| Type | Description |
|---|---|
NDArrayReal
|
Array of time points in seconds, starting from 0.0. |
Examples:
>>> import wandas as wd
>>> signal = wd.read("audio.wav")
>>> time = signal.time
>>> print(f"Duration: {time[-1]:.3f}s")
>>> print(f"Time step: {time[1] - time[0]:.6f}s")
source_time
property
¶
Get sample times relative to the original source timeline.
n_samples
property
¶
Returns the number of samples.
duration
property
¶
Returns the duration in seconds.
rms
property
¶
Calculate one linear RMS amplitude for each channel.
This is a scalar reduction: it computes one value per channel and
triggers immediate computation of the underlying Dask graph. The
result is a plain NumPy array and does not produce a new frame,
so no runtime lineage or operation history view entry is created.
Per-channel calibration factors are applied before the reduction, so
each result uses that channel's physical unit. A calibrated Pa
channel therefore returns RMS pressure in Pa. This property never
performs logarithmic conversion and must not be labeled dB or dB SPL.
The RMS is defined as::
rms[i] = sqrt(mean(x[i] ** 2))
where x[i] is the sample array for channel i.
Returns:
| Type | Description |
|---|---|
NDArrayReal
|
NDArrayReal of shape |
Examples:
>>> import wandas as wd
>>> cf = wd.read("audio.wav")
>>> rms_values = cf.rms
>>> print(f"RMS values: {rms_values}")
>>> # Select channels with RMS > threshold
>>> active_channels = cf[cf.rms > 0.5]
crest_factor
property
¶
Calculate the crest factor (peak-to-RMS ratio) for each channel.
This is a scalar reduction: it computes one value per channel and triggers immediate computation of the underlying Dask graph. The result is a plain NumPy array and does not produce a new frame, so no runtime lineage or operation history view entry is created.
The crest factor is defined as::
crest_factor[i] = max(|x[i]|) / sqrt(mean(x[i] ** 2))
where x[i] is the sample array for channel i.
For a pure sine wave the theoretical continuous-time crest factor is sqrt(2) ≈ 1.414; in discrete-time this implementation typically yields a value close to this, and exactly equal only when the sampled waveform contains its true peaks. Channels with zero RMS (all-zero signals) return 1.0 (defined by convention; no division by zero is performed).
Returns:
| Type | Description |
|---|---|
NDArrayReal
|
NDArrayReal of shape |
Examples:
>>> import wandas as wd
>>> cf = wd.read("audio.wav")
>>> cf_values = cf.crest_factor
>>> print(f"Crest factors: {cf_values}")
>>> # Select channels with crest factor above threshold
>>> impulsive_channels = cf[cf.crest_factor > 3.0]
Functions¶
__init__(data, sampling_rate, label=None, metadata=None, channel_metadata=None, channel_ids=None, previous=None, source_time_offset=0.0, lineage=None, operation_history_prefix=())
¶
Initialize a ChannelFrame.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
data
|
Array
|
Dask array containing channel data. Shape should be (n_channels, n_samples). |
required |
sampling_rate
|
float
|
The sampling rate of the data in Hz. Must be a positive value. |
required |
label
|
str | None
|
A label for the frame. |
None
|
metadata
|
dict[str, Any] | None
|
Optional metadata dictionary. |
None
|
lineage
|
Any | None
|
Runtime operation lineage for this frame. This is the provenance source for executable replay and derived history views. |
None
|
channel_metadata
|
Sequence[ChannelMetadata | dict[str, Any]] | None
|
Metadata for each channel. |
None
|
previous
|
BaseFrame[Any] | None
|
Immediate receiver Frame for process-local data comparison. For multi-input operations, follows only the left/base receiver. Not persisted in WDF. |
None
|
operation_history_prefix
|
Sequence[Mapping[str, Any]]
|
Display history restored at a persistence boundary. |
()
|
Raises:
| Type | Description |
|---|---|
ValueError
|
If data has more than 2 dimensions, or if sampling_rate is not positive. |
Source code in wandas/frames/channel.py
397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 | |
with_calibration(values)
¶
Return a frame configured with replacement per-channel calibrations.
A sequence or one-dimensional NumPy array fully replaces factors in current
channel order. A mapping partially updates labels and/or call-time indices.
Numeric values replace only the factor; :class:ChannelCalibration replaces
factor, unit, and ref. Stored samples stay raw and multiplication remains lazy.
Source code in wandas/frames/channel.py
680 681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 | |
derive_calibration(*, target_rms=None, target_level=None, unit)
¶
Derive absolute per-channel calibration from this reference event.
Exactly one known physical scalar is broadcast to every channel. The
frame is not changed and no operation is added to its history.
target_rms is a linear value in unit. target_level is an
amplitude level using 20 * log10(target_rms / ref) and the default
reference for unit; for unit="Pa" that reference is 2e-5 Pa.
Source code in wandas/frames/channel.py
710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 | |
trigger(threshold, duration, pre=0.0, channel=0, slope='rising', holdoff=None)
¶
Cut out every event where one channel crosses a trigger level.
For transient measurements such as impact hammer tests: each time the
trigger channel crosses threshold, every channel is cut out from
pre seconds before the crossing, for duration seconds. Only the
trigger channel is computed to find the events; the segments stay lazy.
The level is compared with the calibrated channel data, so
threshold is in that channel's physical unit (for example N for a
calibrated force channel). A rising trigger fires at the first sample
that reaches the threshold after a sample below it, and a falling
trigger at the first sample that reaches it after a sample above it.
The crossing sample is usually after the start of the pulse, so give
pre to keep the rise in the segment; plot() on the result shows
whether it is kept. To trigger on a dB level, detect on a level Frame
such as cf.rms_trend(dB=True) and pass its
trigger(...).trigger_times to :meth:segment of this Frame.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
threshold
|
float
|
Trigger level in the trigger channel's physical unit. |
required |
duration
|
float
|
Segment length in seconds, including |
required |
pre
|
float
|
Pre-trigger time in seconds kept before each event. Default is 0.0. |
0.0
|
channel
|
int | str
|
Index or label of the trigger channel. Default is 0. |
0
|
slope
|
TriggerSlope
|
|
'rising'
|
holdoff
|
float | None
|
Minimum time in seconds between two triggers. Crossings
within |
None
|
Returns:
| Type | Description |
|---|---|
ChannelFrameSegments
|
ChannelFrameSegments with one segment per event, as returned by
:meth: |
Raises:
| Type | Description |
|---|---|
IndexError
|
If |
KeyError
|
If |
ValueError
|
If |
Examples:
>>> import wandas as wd
>>> hits = (
... wd.read("hammer_test.wav")
... .high_pass_filter(5)
... .trigger(10.0, duration=1.0, pre=0.01, channel="force")
... )
>>> hits.plot()
>>> frf = hits.drop([2]).transfer_function()
Source code in wandas/frames/channel.py
853 854 855 856 857 858 859 860 861 862 863 864 865 866 867 868 869 870 871 872 873 874 875 876 877 878 879 880 881 882 883 884 885 886 887 888 889 890 891 892 893 894 895 896 897 898 899 900 901 902 903 904 905 906 907 908 909 910 911 912 913 914 915 916 917 918 919 920 921 922 | |
segment(times, duration, pre=0.0)
¶
Cut fixed-length segments of all channels around the given times.
Use this when the event times come from elsewhere, such as a log or a
detection on another Frame; :meth:trigger finds them from a level.
Each segment starts pre seconds before its time and lasts
duration seconds. Every segment is a lazy time slice of this Frame
with the same contract as :meth:trim: channel labels, calibration,
metadata and channel IDs are preserved, and source_time_offset
advances to the first sample of the segment, which records where the
event lies on the original recording. No data is computed.
Each segment's event ID is the position of its time in times.
Times whose segment would start before the first sample or end after
the last one are skipped with a UserWarning, so every segment has
exactly round(duration * sampling_rate) samples and the IDs of the
kept segments may have gaps.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
times
|
float | Sequence[float] | NDArrayReal
|
Event times in seconds from the start of this Frame. |
required |
duration
|
float
|
Segment length in seconds, including |
required |
pre
|
float
|
Pre-trigger time in seconds kept before each event. Default is 0.0. |
0.0
|
Returns:
| Type | Description |
|---|---|
ChannelFrameSegments
|
ChannelFrameSegments with one segment per kept event, in the order
of |
Raises:
| Type | Description |
|---|---|
ValueError
|
If |
Examples:
>>> import wandas as wd
>>> cf = wd.read("hammer_test.wav")
>>> hits = cf.segment([0.5, 1.7, 2.9], duration=1.0, pre=0.01)
>>> frf = hits.transfer_function(input="force", output="accel")
Source code in wandas/frames/channel.py
947 948 949 950 951 952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967 968 969 970 971 972 973 974 975 976 977 978 979 980 981 982 983 984 985 986 987 988 989 990 991 992 993 994 995 996 997 998 | |
info()
¶
Display comprehensive information about the ChannelFrame.
This method prints a summary of the frame's properties including: - Number of channels - Sampling rate - Duration - Number of samples - Channel labels
This is a convenience method to view all key properties at once, similar to pandas DataFrame.info().
Examples:
>>> import wandas as wd
>>> cf = wd.read("audio.wav")
>>> cf.info()
Channels: 2
Sampling rate: 44100 Hz
Duration: 1.0 s
Samples: 44100
Channel labels: ['ch0', 'ch1']
Source code in wandas/frames/channel.py
1036 1037 1038 1039 1040 1041 1042 1043 1044 1045 1046 1047 1048 1049 1050 1051 1052 1053 1054 1055 1056 1057 1058 1059 1060 1061 1062 1063 1064 1065 | |
mix(other, *, align='strict', snr_db=None)
¶
Mix another signal lazily by array index.
The base frame owns output length, channels, metadata, labels, and
source_time_offset. pad accepts only a shorter other signal;
truncate accepts only a longer one.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
other
|
ChannelFrame | NDArrayReal | Array
|
A ChannelFrame or one-/two-dimensional NumPy or Dask array. A single channel broadcasts across the base channels; otherwise channel counts must match. |
required |
align
|
Literal['strict', 'pad', 'truncate']
|
Length policy. |
'strict'
|
snr_db
|
float | None
|
Optional signal-to-noise ratio in decibels used to scale |
None
|
Returns:
| Type | Description |
|---|---|
ChannelFrame
|
A new lazy ChannelFrame with the base frame's structure and metadata. |
Raises:
| Type | Description |
|---|---|
TypeError
|
If |
ValueError
|
If sampling rates, dimensions, channel counts, lengths, or
|
Notes
Source-time offsets describe provenance and do not shift array positions. Signals from different source-time regions can therefore be mixed directly.
Source code in wandas/frames/channel.py
1075 1076 1077 1078 1079 1080 1081 1082 1083 1084 1085 1086 1087 1088 1089 1090 1091 1092 1093 1094 1095 1096 1097 1098 1099 1100 1101 1102 1103 1104 1105 1106 1107 1108 1109 1110 1111 1112 1113 1114 1115 1116 1117 1118 1119 1120 1121 1122 1123 1124 1125 1126 1127 1128 1129 1130 1131 1132 1133 1134 1135 1136 1137 1138 1139 1140 1141 1142 1143 1144 1145 1146 1147 1148 1149 1150 1151 1152 1153 1154 1155 1156 1157 1158 1159 1160 1161 1162 1163 1164 1165 1166 1167 | |
plot(plot_type='waveform', ax=None, title=None, overlay=False, xlabel=None, ylabel=None, alpha=1.0, xlim=None, ylim=None, **kwargs)
¶
Plot the frame data.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
plot_type
|
str
|
Type of plot. Default is "waveform". |
'waveform'
|
ax
|
Axes | None
|
Optional matplotlib axes for plotting. |
None
|
title
|
str | None
|
Title for the plot. If None, uses the frame label. |
None
|
overlay
|
bool
|
Whether to overlay all channels on a single plot (True) or create separate subplots for each channel (False). |
False
|
xlabel
|
str | None
|
Label for the x-axis. If None, uses default based on plot type. |
None
|
ylabel
|
str | None
|
Label for the y-axis. If None, uses default based on plot type. |
None
|
alpha
|
float
|
Transparency level for the plot lines (0.0 to 1.0). |
1.0
|
xlim
|
tuple[float, float] | None
|
Limits for the x-axis as (min, max) tuple. |
None
|
ylim
|
tuple[float, float] | None
|
Limits for the y-axis as (min, max) tuple. |
None
|
**kwargs
|
Any
|
Additional matplotlib Line2D parameters
(e.g., color, linewidth, linestyle).
These are passed to the underlying matplotlib plot functions.
The waveform plot also accepts |
{}
|
Returns:
| Type | Description |
|---|---|
Axes | Iterator[Axes]
|
Single Axes object or iterator of Axes objects. |
Examples:
>>> import wandas as wd
>>> cf = wd.read("audio.wav")
>>> # Basic plot
>>> cf.plot()
>>> # Overlay all channels
>>> cf.plot(overlay=True, alpha=0.7)
>>> # Custom styling
>>> cf.plot(title="My Signal", ylabel="Voltage [V]", color="red")
>>> # Draw every sample instead of the envelope fitted to the axes
>>> cf.plot(max_points=None)
Source code in wandas/frames/channel.py
1169 1170 1171 1172 1173 1174 1175 1176 1177 1178 1179 1180 1181 1182 1183 1184 1185 1186 1187 1188 1189 1190 1191 1192 1193 1194 1195 1196 1197 1198 1199 1200 1201 1202 1203 1204 1205 1206 1207 1208 1209 1210 1211 1212 1213 1214 1215 1216 1217 1218 1219 1220 1221 1222 1223 1224 1225 1226 1227 1228 1229 1230 1231 1232 1233 1234 1235 1236 1237 1238 1239 1240 1241 1242 1243 1244 1245 1246 1247 1248 1249 1250 1251 1252 1253 1254 | |
rms_plot(ax=None, title=None, overlay=True, Aw=False, **kwargs)
¶
Plot a windowed RMS amplitude level in decibels.
rms_plot() calls :meth:rms_trend with dB=True. It is not a
plot of the linear :attr:rms scalar. Values use
20 * log10(window_rms / channel_ref); Aw=True applies the
implemented digital A-weighting filter before the RMS calculation.
The implementation has not been validated as a conforming sound-level
meter.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
ax
|
Axes | None
|
Optional matplotlib axes for plotting. |
None
|
title
|
str | None
|
Title for the plot. |
None
|
overlay
|
bool
|
Whether to overlay the plot on the existing axis. |
True
|
Aw
|
bool
|
Apply the implemented A-frequency-weighting filter. |
False
|
**kwargs
|
Any
|
Additional arguments passed to the plot() method. Accepts the same arguments as plot() including xlabel, ylabel, alpha, xlim, ylim, and matplotlib Line2D parameters. |
{}
|
Returns:
| Type | Description |
|---|---|
Axes | Iterator[Axes]
|
Single Axes object or iterator of Axes objects. |
Examples:
>>> cf = wd.read("audio.wav")
>>> # Basic RMS plot
>>> cf.rms_plot()
>>> # With A-weighting
>>> cf.rms_plot(Aw=True)
>>> # Custom styling
>>> cf.rms_plot(ylabel="RMS level [dB re channel reference]", alpha=0.8, color="blue")
Source code in wandas/frames/channel.py
1256 1257 1258 1259 1260 1261 1262 1263 1264 1265 1266 1267 1268 1269 1270 1271 1272 1273 1274 1275 1276 1277 1278 1279 1280 1281 1282 1283 1284 1285 1286 1287 1288 1289 1290 1291 1292 1293 1294 1295 1296 1297 1298 1299 1300 | |
describe(normalize=True, is_close=True, *, fmin=0, fmax=None, cmap='jet', vmin=None, vmax=None, xlim=None, ylim=None, Aw=False, waveform=None, spectral=None, image_save=None, **kwargs)
¶
Display visual and audio representation of the frame.
This method creates a comprehensive visualization with three plots: 1. Time-domain waveform (top) 2. Spectrogram (bottom-left) 3. Frequency spectrum via Welch method (bottom-right)
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
normalize
|
bool
|
Whether to normalize the audio data for playback. Default: True |
True
|
is_close
|
bool
|
Whether to close the figure after displaying.
Default: True. In a notebook (Jupyter, VS Code, Colab or
marimo) with IPython, each figure is shown with its audio
player. In a PyScript page whose running code has an output
element, each figure and its audio player are added to that
element with |
True
|
fmin
|
float
|
Minimum frequency to display in the spectrogram (Hz). Default: 0 |
0
|
fmax
|
float | None
|
Maximum frequency to display in the spectrogram (Hz). Default: Nyquist frequency (sampling_rate / 2) |
None
|
cmap
|
str
|
Colormap for the spectrogram. Default: 'jet' |
'jet'
|
vmin
|
float | None
|
Minimum value for spectrogram color scale (dB). Auto-calculated if None. |
None
|
vmax
|
float | None
|
Maximum value for spectrogram color scale (dB). Auto-calculated if None. |
None
|
xlim
|
tuple[float, float] | None
|
Time axis limits (seconds) for all time-based plots. Format: (start_time, end_time) |
None
|
ylim
|
tuple[float, float] | None
|
Frequency axis limits (Hz) for frequency-based plots. Format: (min_freq, max_freq) |
None
|
Aw
|
bool
|
Apply A-weighting to the frequency analysis. Default: False |
False
|
waveform
|
dict[str, Any] | None
|
Additional configuration dict for waveform subplot. Can include 'xlabel', 'ylabel', 'xlim', 'ylim'. |
None
|
spectral
|
dict[str, Any] | None
|
Additional configuration dict for spectral subplot. Can include 'xlabel', 'ylabel', 'xlim', 'ylim'. |
None
|
image_save
|
str | Path | None
|
Path to save the figure as an image file. If provided, the figure will be saved before closing. File format is determined from the extension (e.g., '.png', '.jpg', '.pdf'). For multi-channel frames, the channel index is appended to the filename stem (e.g., 'output_0.png', 'output_1.png'). Default: None. |
None
|
**kwargs
|
Any
|
Deprecated parameters for backward compatibility only. - axis_config: Old configuration format (use waveform/spectral instead) - cbar_config: Old colorbar configuration (use vmin/vmax instead) |
{}
|
Returns:
| Type | Description |
|---|---|
list[Figure] | None
|
None (default). When |
Examples:
>>> cf = wd.read("audio.wav")
>>> # Basic usage
>>> cf.describe()
>>>
>>> # Custom frequency range
>>> cf.describe(fmin=100, fmax=5000)
>>>
>>> # Custom color scale
>>> cf.describe(vmin=-80, vmax=-20, cmap="viridis")
>>>
>>> # A-weighted analysis
>>> cf.describe(Aw=True)
>>>
>>> # Custom time range
>>> cf.describe(xlim=(0, 5)) # Show first 5 seconds
>>>
>>> # Custom waveform subplot settings
>>> cf.describe(waveform={"ylabel": "Custom Label"})
>>>
>>> # Save the figure to a file
>>> cf.describe(image_save="output.png")
>>>
>>> # Get Figure objects for further manipulation (is_close=False)
>>> figures = cf.describe(is_close=False)
>>> fig = figures[0]
>>> fig.savefig("custom_output.png") # Custom save with modifications
Source code in wandas/frames/channel.py
1302 1303 1304 1305 1306 1307 1308 1309 1310 1311 1312 1313 1314 1315 1316 1317 1318 1319 1320 1321 1322 1323 1324 1325 1326 1327 1328 1329 1330 1331 1332 1333 1334 1335 1336 1337 1338 1339 1340 1341 1342 1343 1344 1345 1346 1347 1348 1349 1350 1351 1352 1353 1354 1355 1356 1357 1358 1359 1360 1361 1362 1363 1364 1365 1366 1367 1368 1369 1370 1371 1372 1373 1374 1375 1376 1377 1378 1379 1380 1381 1382 1383 1384 1385 1386 1387 1388 1389 1390 1391 1392 1393 1394 1395 1396 1397 1398 1399 1400 1401 1402 1403 1404 1405 1406 1407 1408 1409 1410 1411 1412 1413 1414 1415 1416 1417 1418 1419 1420 1421 | |
from_numpy(data, sampling_rate, label=None, metadata=None, ch_labels=None, ch_units=None)
classmethod
¶
Create a ChannelFrame from a NumPy array.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
data
|
NDArrayReal
|
NumPy array containing channel data. |
required |
sampling_rate
|
float
|
The sampling rate in Hz. |
required |
label
|
str | None
|
A label for the frame. |
None
|
metadata
|
dict[str, Any] | None
|
Optional metadata dictionary. |
None
|
ch_labels
|
list[str] | None
|
Labels for each channel. |
None
|
ch_units
|
list[str] | str | None
|
Units for each channel. |
None
|
Returns:
| Type | Description |
|---|---|
ChannelFrame
|
A new ChannelFrame containing the NumPy data. |
Source code in wandas/frames/channel.py
1423 1424 1425 1426 1427 1428 1429 1430 1431 1432 1433 1434 1435 1436 1437 1438 1439 1440 1441 1442 1443 1444 1445 1446 1447 1448 1449 1450 1451 1452 1453 1454 1455 1456 1457 1458 1459 1460 1461 1462 1463 1464 1465 1466 1467 1468 | |
from_ndarray(array, sampling_rate, labels=None, unit=None, frame_label=None, metadata=None)
classmethod
¶
Create a ChannelFrame from a NumPy array.
This method is deprecated. Use from_numpy instead.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
array
|
NDArrayReal
|
Signal data. Each row corresponds to a channel. |
required |
sampling_rate
|
float
|
Sampling rate (Hz). |
required |
labels
|
list[str] | None
|
Labels for each channel. |
None
|
unit
|
list[str] | str | None
|
Unit of the signal. |
None
|
frame_label
|
str | None
|
Label for the frame. |
None
|
metadata
|
dict[str, Any] | None
|
Optional metadata dictionary. |
None
|
Returns:
| Type | Description |
|---|---|
ChannelFrame
|
A new ChannelFrame containing the data. |
Source code in wandas/frames/channel.py
1470 1471 1472 1473 1474 1475 1476 1477 1478 1479 1480 1481 1482 1483 1484 1485 1486 1487 1488 1489 1490 1491 1492 1493 1494 1495 1496 1497 1498 1499 1500 1501 1502 1503 1504 1505 1506 1507 | |
from_file(path, channel=None, start=None, end=None, ch_labels=None, time_column=0, delimiter=',', header=0, file_type=None, source_name=None, timeout=10.0, dtype=None)
classmethod
¶
Create a ChannelFrame from an audio file or URL.
Note
The chunk_size parameter has been removed. ChannelFrame uses
channel-wise chunking by default (chunks=(1, -1)). Use .rechunk(...)
on the returned frame for custom sample-axis chunking.
Audio sample decoding is deferred through the returned Dask array. CSV metadata inspection synchronously parses the complete table to determine shape and sampling rate; the sample table is parsed again when the Dask data is computed.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
path
|
str | Path | bytes | bytearray | memoryview | BinaryIO
|
Path to the audio file, in-memory bytes/stream, or an HTTP/HTTPS
URL. When a URL is given it is streamed into a temporary file
before processing, subject to the maximum download size
enforced by |
required |
channel
|
int | list[int] | None
|
Channel(s) to load. None loads all channels. |
None
|
start
|
float | None
|
Start time in seconds. |
None
|
end
|
float | None
|
End time in seconds. |
None
|
ch_labels
|
list[str] | None
|
Labels for each channel. |
None
|
time_column
|
int | str
|
For CSV files, index or name of the time column. Default is 0 (first column). |
0
|
delimiter
|
str
|
For CSV files, delimiter character. Default is ",". |
','
|
header
|
int | None
|
For CSV files, row number to use as header. Default is 0 (first row). Set to None if no header. |
0
|
file_type
|
str | None
|
File extension for in-memory data or URLs without a recognisable extension (e.g. ".wav", ".csv"). |
None
|
source_name
|
str | None
|
Optional source name for in-memory data. Used in metadata. |
None
|
timeout
|
float
|
Timeout in seconds for HTTP/HTTPS URL downloads. Default is 10.0 seconds. Has no effect for local files or in-memory data. |
10.0
|
dtype
|
DTypeLike | None
|
|
None
|
Returns:
| Type | Description |
|---|---|
ChannelFrame
|
A new ChannelFrame containing the loaded data. SoundFile-backed
audio channels use the explicit linear unit |
Raises:
| Type | Description |
|---|---|
ValueError
|
If channel specification is invalid or file cannot be read. Error message includes absolute path, current directory, and troubleshooting suggestions. |
Examples:
>>> import wandas as wd
>>> # Load WAV file as full-scale float64 audio
>>> cf = wd.read("audio.wav")
>>> # Decode as float32 to halve the memory
>>> cf = wd.read("audio.wav", dtype="float32")
>>> # Load specific channels
>>> cf = wd.read("audio.wav", channel=[0, 2])
>>> # Load CSV file
>>> cf = wd.read("data.csv", time_column=0, delimiter=",", header=0)
>>> # Load from a URL
>>> cf = wd.read("https://example.com/audio.wav")
Source code in wandas/frames/channel.py
1509 1510 1511 1512 1513 1514 1515 1516 1517 1518 1519 1520 1521 1522 1523 1524 1525 1526 1527 1528 1529 1530 1531 1532 1533 1534 1535 1536 1537 1538 1539 1540 1541 1542 1543 1544 1545 1546 1547 1548 1549 1550 1551 1552 1553 1554 1555 1556 1557 1558 1559 1560 1561 1562 1563 1564 1565 1566 1567 1568 1569 1570 1571 1572 1573 1574 1575 1576 1577 1578 1579 1580 1581 1582 1583 1584 1585 1586 1587 1588 1589 1590 1591 1592 1593 1594 1595 1596 1597 1598 1599 1600 1601 1602 1603 1604 1605 1606 1607 1608 1609 1610 1611 1612 1613 1614 1615 1616 1617 1618 1619 1620 1621 1622 1623 1624 1625 1626 1627 1628 1629 1630 1631 1632 1633 1634 1635 1636 1637 1638 1639 1640 1641 1642 1643 1644 1645 1646 1647 1648 1649 1650 1651 1652 1653 1654 1655 1656 1657 1658 1659 1660 1661 1662 1663 1664 1665 1666 1667 1668 1669 1670 1671 1672 1673 1674 1675 1676 1677 1678 1679 1680 1681 1682 1683 1684 1685 1686 1687 1688 1689 1690 1691 1692 1693 1694 1695 1696 1697 1698 1699 1700 1701 1702 1703 1704 1705 1706 1707 1708 1709 1710 1711 1712 1713 1714 1715 1716 1717 1718 1719 1720 1721 1722 1723 1724 1725 1726 1727 1728 1729 1730 1731 1732 1733 1734 1735 1736 1737 1738 1739 1740 1741 1742 1743 1744 1745 1746 1747 1748 1749 1750 1751 1752 1753 1754 1755 1756 1757 1758 1759 1760 1761 1762 1763 1764 1765 1766 1767 1768 1769 | |
read_wav(filename, labels=None)
classmethod
¶
Utility method to read a WAV file.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
filename
|
str | Path | bytes | bytearray | memoryview | BinaryIO
|
Path to the WAV file or in-memory bytes/stream. |
required |
labels
|
list[str] | None
|
Labels to set for each channel. |
None
|
Returns:
| Type | Description |
|---|---|
ChannelFrame
|
A new ChannelFrame containing the data (lazy loading). |
Source code in wandas/frames/channel.py
1771 1772 1773 1774 1775 1776 1777 1778 1779 1780 1781 1782 1783 1784 1785 1786 1787 1788 1789 1790 1791 1792 1793 1794 1795 1796 1797 1798 | |
read_csv(filename, time_column=0, labels=None, delimiter=',', header=0)
classmethod
¶
Utility method to read a CSV file.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
filename
|
str
|
Path to the CSV file. |
required |
time_column
|
int | str
|
Index or name of the time column. |
0
|
labels
|
list[str] | None
|
Labels to set for each channel. |
None
|
delimiter
|
str
|
Delimiter character. |
','
|
header
|
int | None
|
Row number to use as header. |
0
|
Returns:
| Type | Description |
|---|---|
ChannelFrame
|
A new ChannelFrame containing Dask-backed sample data. CSV metadata inspection occurs synchronously before the Frame is returned. |
Examples:
>>> # Read CSV with default settings
>>> cf = wd.read("data.csv")
>>> # Read CSV with custom delimiter
>>> cf = wd.read("data.csv", delimiter=";")
>>> # Read CSV without header
>>> cf = wd.read("data.csv", header=None)
Source code in wandas/frames/channel.py
1800 1801 1802 1803 1804 1805 1806 1807 1808 1809 1810 1811 1812 1813 1814 1815 1816 1817 1818 1819 1820 1821 1822 1823 1824 1825 1826 1827 1828 1829 1830 1831 1832 1833 1834 1835 1836 1837 1838 1839 | |
to_wav(path, format=None)
¶
Save the audio data to a WAV file.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
path
|
str | PathLike[str] | BinaryIO
|
Destination path or writable binary stream, such as BytesIO. |
required |
format
|
str | None
|
Required for streams. Paths infer the format from their extension. |
None
|
Raises:
| Type | Description |
|---|---|
ValueError
|
If a stream is supplied without an explicit format. |
Source code in wandas/frames/channel.py
1841 1842 1843 1844 1845 1846 1847 1848 1849 1850 1851 1852 1853 | |
load(path)
classmethod
¶
Load a ChannelFrame from a WDF (Wandas Data File) file.
This loads data saved with the save() method, preserving all channel data, metadata, labels, and units.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
path
|
str | Path
|
Path to the WDF file |
required |
Returns:
| Type | Description |
|---|---|
ChannelFrame
|
A new ChannelFrame with all data and metadata loaded |
Raises:
| Type | Description |
|---|---|
FileNotFoundError
|
If the file doesn't exist |
Examples:
>>> import wandas as wd
>>> cf = wd.load("audio_analysis.wdf")
Source code in wandas/frames/channel.py
1855 1856 1857 1858 1859 1860 1861 1862 1863 1864 1865 1866 1867 1868 1869 1870 1871 1872 1873 1874 1875 1876 1877 1878 1879 1880 1881 1882 1883 1884 1885 | |
add_channel(data, label=None, align='strict', suffix_on_dup=None, source_time_offset=None)
¶
Add a new channel to the frame.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
data
|
ndarray[Any, Any] | Array
|
NumPy or Dask data for exactly one channel. The accepted shapes
are |
required |
label
|
str | None
|
Label for the new channel. If None, generates a default label. |
None
|
align
|
Literal['strict', 'pad', 'truncate']
|
How to handle length mismatches: - "strict": Raise error if lengths don't match - "pad": Pad shorter data with zeros - "truncate": Truncate longer data to match |
'strict'
|
suffix_on_dup
|
str | None
|
Suffix to add to duplicate labels. If None, raises error. |
None
|
source_time_offset
|
float | Sequence[float] | NDArrayReal | None
|
Offset in seconds for the new channel. Accepts a finite real scalar or a one-item 1-D sequence/NumPy array. If None, the new channel uses 0.0. Accepted forms are normalized to a built-in float before execution and Recipe capture. |
None
|
Returns:
| Type | Description |
|---|---|
ChannelFrame
|
A new ChannelFrame. |
Raises:
| Type | Description |
|---|---|
ValueError
|
If data length doesn't match and align="strict", or if label is duplicate and suffix_on_dup is None. |
TypeError
|
If data is not a NumPy or Dask array. Pass another
ChannelFrame to :meth: |
Examples:
>>> cf = wd.read("audio.wav")
>>> # Add a numpy array as a new channel
>>> new_data = np.sin(2 * np.pi * 440 * cf.time)
>>> cf_new = cf.add_channel(new_data, label="sine_440Hz")
>>> # Concatenate another ChannelFrame's channels
>>> cf2 = wd.read("audio2.wav")
>>> cf_combined = cf.concat_frame(cf2)
Source code in wandas/frames/channel.py
1887 1888 1889 1890 1891 1892 1893 1894 1895 1896 1897 1898 1899 1900 1901 1902 1903 1904 1905 1906 1907 1908 1909 1910 1911 1912 1913 1914 1915 1916 1917 1918 1919 1920 1921 1922 1923 1924 1925 1926 1927 1928 1929 1930 1931 1932 1933 1934 1935 1936 1937 1938 1939 1940 1941 1942 1943 1944 1945 1946 1947 1948 1949 1950 1951 1952 1953 1954 1955 1956 1957 1958 1959 1960 1961 1962 1963 1964 1965 1966 1967 1968 1969 1970 1971 1972 1973 1974 1975 1976 1977 1978 1979 1980 1981 1982 1983 1984 1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 | |
concat_frame(other, label_prefix=None, align='strict', suffix_on_dup=None)
¶
Concatenate all channels from another frame along the channel axis.
The result preserves other channel metadata, calibration, and
source-time offsets. Neither input frame is changed.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
other
|
ChannelFrame
|
ChannelFrame whose channels are appended in their current order. |
required |
label_prefix
|
str | None
|
Optional prefix for appended labels, producing
|
None
|
align
|
Literal['strict', 'pad', 'truncate']
|
|
'strict'
|
suffix_on_dup
|
str | None
|
Suffix for duplicate labels. If None, duplicates raise. |
None
|
Returns:
| Type | Description |
|---|---|
ChannelFrame
|
A new lazy ChannelFrame containing both channel collections. |
Raises:
| Type | Description |
|---|---|
TypeError
|
If other is not a ChannelFrame. |
ValueError
|
If sampling rates, lengths, or labels are incompatible. |
Source code in wandas/frames/channel.py
1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 2019 2020 2021 2022 2023 2024 2025 2026 2027 2028 2029 2030 2031 2032 2033 2034 2035 2036 2037 2038 2039 2040 2041 2042 2043 2044 2045 2046 2047 2048 2049 2050 2051 2052 2053 2054 2055 2056 2057 2058 2059 2060 2061 2062 2063 2064 2065 2066 2067 2068 2069 2070 2071 2072 2073 2074 | |
remove_channel(key)
¶
Return a new frame without one channel.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
key
|
int | str
|
Zero-based channel index or exact channel label to remove. |
required |
Returns:
| Type | Description |
|---|---|
ChannelFrame
|
A lazy ChannelFrame preserving the remaining channels' metadata, stable channel identifiers, source-time offsets, and semantic lineage. |
Raises:
| Type | Description |
|---|---|
IndexError
|
If an integer index is outside the channel range. |
KeyError
|
If a string label does not exist. |
Source code in wandas/frames/channel.py
2076 2077 2078 2079 2080 2081 2082 2083 2084 2085 2086 2087 2088 2089 2090 2091 2092 2093 2094 2095 2096 2097 2098 2099 2100 2101 2102 2103 2104 2105 2106 2107 2108 2109 2110 | |
wandas.frames.spectral.SpectralFrame
¶
Bases: SpectralPropertiesMixin, BaseFrame[NDArrayComplex]
Class for handling frequency-domain signal data.
This class represents spectral data, providing methods for spectral analysis, manipulation, and visualization. It handles complex-valued frequency domain data obtained through operations like FFT.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
data
|
Array
|
DaArray. The spectral data. Must be a dask array with shape: - (channels, frequency_bins) for multi-channel data - (frequency_bins,) for single-channel data, which will be reshaped to (1, frequency_bins) |
required |
sampling_rate
|
float
|
float. The sampling rate of the original time-domain signal in Hz. |
required |
n_fft
|
int
|
int. Required. The FFT size used to generate this spectral data. Must be a
positive integer and must be the same FFT size that was used to create
the complete one-sided spectrum (for example, 512 or 1024). The data must
contain exactly |
required |
window
|
str
|
str, default="hann". The window function used in the FFT. |
'hann'
|
win_length
|
int | None
|
int, optional. Length of the analysis window applied before
zero-padding to |
None
|
label
|
str | None
|
str, optional. A label for the frame. |
None
|
metadata
|
dict[str, Any] | None
|
dict, optional. Additional metadata for the frame. |
None
|
lineage
|
Any | None
|
LineageNode, optional. Constructor override for the runtime lineage. When omitted, a source node is
created. |
None
|
channel_metadata
|
Sequence[ChannelMetadata | dict[str, Any]] | None
|
list[ChannelMetadata], optional. Metadata for each channel in the frame. |
None
|
previous
|
BaseFrame[Any] | None
|
BaseFrame, optional. Immediate receiver Frame for process-local data comparison. For multi-input operations, follows only the left/base receiver. Not persisted in WDF. |
None
|
Attributes:
| Name | Type | Description |
|---|---|---|
magnitude |
NDArrayReal
|
NDArrayReal. Absolute value of the stored spectral quantity. FFT and Welch results are amplitudes in the input channel unit. |
phase |
NDArrayReal
|
NDArrayReal. The phase spectrum in radians. |
unwrapped_phase |
NDArrayReal
|
NDArrayReal. The unwrapped phase spectrum in radians. |
power |
NDArrayReal
|
NDArrayReal. Squared magnitude. This compatibility property is not necessarily physical power or power spectral density. |
dB |
NDArrayReal
|
NDArrayReal. Magnitude level, |
dBA |
NDArrayReal
|
NDArrayReal. A-weighted magnitude level. For FFT and Welch results this is an A-weighted amplitude level. |
freqs |
NDArrayReal
|
NDArrayReal. The frequency axis values in Hz. |
Examples:
Create a SpectralFrame from FFT:
>>> signal = ChannelFrame.from_numpy(data, sampling_rate=44100)
>>> spectrum = signal.fft(n_fft=2048)
Plot the amplitude level spectrum:
>>> spectrum.plot()
Perform binary operations:
>>> scaled = spectrum * 2.0
>>> summed = spectrum1 + spectrum2 # Must have matching sampling rates
Convert back to time domain:
>>> time_signal = spectrum.ifft()
Notes
- All operations are performed lazily using dask arrays for efficient memory usage.
- Binary operations (+, -, *, /) can be performed between SpectralFrames or with
scalar values. - The class maintains runtime lineage and metadata through all operations.
Source code in wandas/frames/spectral.py
33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 | |
Attributes¶
n_fft
property
¶
Return the read-only FFT size defining the frequency axis.
window
property
¶
Return the read-only analysis window; rerun FFT to change it.
win_length
property
¶
Return the read-only length of the analysis window before zero-padding.
unwrapped_phase
property
¶
Get the unwrapped phase spectrum.
The unwrapped phase removes discontinuities of 2π radians, providing continuous phase values across frequency bins.
Returns:
| Name | Type | Description |
|---|---|---|
NDArrayReal |
NDArrayReal
|
The unwrapped phase angles of the complex spectrum in radians. |
freqs
property
¶
Get the frequency axis values in Hz.
Values are derived on access from sampling_rate and n_fft using the
canonical one-sided real-FFT grid. They are not duplicated in Frame state.
Returns:
| Name | Type | Description |
|---|---|---|
NDArrayReal |
NDArrayReal
|
Array of frequency values corresponding to each frequency bin. |
Functions¶
__init__(data, sampling_rate, n_fft, window='hann', label=None, metadata=None, channel_metadata=None, channel_ids=None, previous=None, source_time_offset=0.0, lineage=None, operation_history_prefix=(), win_length=None)
¶
Initialize a complete canonical one-sided spectrum.
See the class docstring for parameter descriptions. The frequency axis is
derived from sampling_rate and n_fft rather than stored as mutable
coordinate state.
Source code in wandas/frames/spectral.py
102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 | |
plot(plot_type='frequency', ax=None, title=None, overlay=False, xlabel=None, ylabel=None, alpha=1.0, xlim=None, ylim=None, Aw=False, **kwargs)
¶
Plot the spectral data using various visualization strategies.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
plot_type
|
str
|
str, default="frequency". Type of plot to create. Options include: - "frequency": Standard frequency plot - "matrix": Matrix plot for comparing channels - Other types as defined by available plot strategies |
'frequency'
|
ax
|
Axes | None
|
matplotlib.axes.Axes, optional. Axes to plot on. If None, creates new axes. |
None
|
title
|
str | None
|
str, optional. Title for the plot. If None, uses the frame label. |
None
|
overlay
|
bool
|
bool, default=False. Whether to overlay all channels on a single plot (True) or create separate subplots for each channel (False). |
False
|
xlabel
|
str | None
|
str, optional. Label for the x-axis. If None, uses default "Frequency [Hz]". |
None
|
ylabel
|
str | None
|
str, optional. Label for the y-axis. If None, uses default based on data type. |
None
|
alpha
|
float
|
float, default=1.0. Transparency level for the plot lines (0.0 to 1.0). |
1.0
|
xlim
|
tuple[float, float] | None
|
tuple[float, float], optional. Limits for the x-axis as (min, max) tuple. |
None
|
ylim
|
tuple[float, float] | None
|
tuple[float, float], optional. Limits for the y-axis as (min, max) tuple. |
None
|
Aw
|
bool
|
bool, default=False. Whether to apply A-weighting to the data. |
False
|
**kwargs
|
Any
|
dict. Additional matplotlib Line2D parameters (e.g., color, linewidth, linestyle). |
{}
|
Returns:
| Type | Description |
|---|---|
Axes | Iterator[Axes]
|
Union[Axes, Iterator[Axes]]: The matplotlib axes containing the plot, or an iterator of axes for multi-plot outputs. |
Examples:
>>> spectrum = cf.fft()
>>> # Basic frequency plot
>>> spectrum.plot()
>>> # Overlay with A-weighting
>>> spectrum.plot(overlay=True, Aw=True)
>>> # Custom styling
>>> spectrum.plot(title="Frequency Spectrum", color="red", linewidth=2)
Source code in wandas/frames/spectral.py
238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 | |
ifft()
¶
Invert Wandas FFT normalization to a windowed time-domain signal.
For a spectrum returned by ChannelFrame.fft() with matching stored
n_fft, window and win_length, this returns the analysed
samples multiplied by the FFT window and zero-padded to n_fft. With
window="boxcar",
that prepared input is reconstructed exactly. A tapered window such as
the default Hann window is not divided out because its zero-valued
samples cannot be recovered. Graph construction remains lazy.
Returns:
| Name | Type | Description |
|---|---|---|
ChannelFrame |
ChannelFrame
|
A new ChannelFrame containing the windowed time-domain signal in the original channel unit. |
Source code in wandas/frames/spectral.py
318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 | |
noct_synthesis(fmin, fmax, n=3, G=10, fr=1000)
¶
Synthesize N-octave band spectrum.
This method combines frequency components into N-octave bands according to
standard acoustical band definitions. This is commonly used in noise and
vibration analysis. The authoritative FFT size and analysis window are
read from this frame's n_fft, window and win_length; they are
not caller-supplied parameters. The peak-amplitude bins are converted to
RMS and divided by the square root of the window's equivalent noise
bandwidth, so band levels match :meth:ChannelFrame.noct_spectrum.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
fmin
|
float
|
float. Lower frequency bound in Hz. |
required |
fmax
|
float
|
float. Upper frequency bound in Hz. |
required |
n
|
int
|
int, default=3. Number of bands per octave (e.g., 3 for third-octave bands). |
3
|
G
|
int
|
int, default=10. Exact center-frequency ratio convention. Use 10 for base
|
10
|
fr
|
int
|
int, default=1000. Reference frequency in Hz. |
1000
|
Returns:
| Name | Type | Description |
|---|---|---|
NOctFrame |
NOctFrame
|
A new NOctFrame containing the N-octave band spectrum. |
Raises:
| Type | Description |
|---|---|
ValueError
|
If the sampling rate is not 48000 Hz. |
TypeError
|
If |
Source code in wandas/frames/spectral.py
400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 | |
plot_matrix(plot_type='matrix', **kwargs)
¶
Plot channel relationships in matrix format.
This method creates a matrix plot showing relationships between channels, such as coherence, transfer functions, or cross-spectral density.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
plot_type
|
str
|
str, default="matrix". Type of matrix plot to create. |
'matrix'
|
**kwargs
|
Any
|
dict. Additional plot parameters: - vmin, vmax: Color scale limits - cmap: Colormap name - title: Plot title |
{}
|
Returns:
| Type | Description |
|---|---|
Axes | Iterator[Axes]
|
Union[Axes, Iterator[Axes]]: The matplotlib axes containing the plot. |
Source code in wandas/frames/spectral.py
538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 | |
info()
¶
Display comprehensive information about the SpectralFrame.
This method prints a summary of the frame's properties including: - Number of channels - Sampling rate - FFT size - Frequency range - Number of frequency bins - Frequency resolution (ΔF) - Channel labels
This is a convenience method to view all key properties at once, similar to pandas DataFrame.info().
Examples:
>>> spectrum = cf.fft()
>>> spectrum.info()
SpectralFrame Information:
Channels: 2
Sampling rate: 44100 Hz
FFT size: 2048
Frequency range: 0.0 - 22050.0 Hz
Frequency bins: 1025
Frequency resolution (ΔF): 21.5 Hz
Channel labels: ['ch0', 'ch1']
Operations Applied: 1
Source code in wandas/frames/spectral.py
573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 | |
wandas.frames.pairwise.CoherenceFrame
¶
Bases: PairwiseSpectralFrame
Typed storage for dimensionless magnitude-squared coherence.
data is real numeric, rank one or two, and uses flattened
(pair, frequency) storage internally. A single pair is exposed with
the normal single-channel public shape. n_fft, window, and
pair_state are required constructor state; frequency_indices may
retain a selected, ordered subset of the n_fft // 2 + 1 frequency
bins. pair_state supplies the output/input roles, source identity,
dimensionless domain, and row order; labels and lineage do not define
quantity meaning.
Magnitude-squared coherence produced by :meth:ChannelFrame.coherence is
mathematically in [0, 1], with NaN for undefined bins. This constructor
validates array structure and typed state, but does not scan, clip, or
otherwise validate array values from direct construction, WDF decoding, or
external Dask inputs. The constructor and all public operations remain lazy
for Dask input. Pair selection, frequency slicing, metadata changes, and
annotation copies preserve this concrete type and the corresponding typed
rows; arithmetic, amplitude-level APIs, inverse FFT, synthesis, and
A-weighting are not defined for this Frame.
Source code in wandas/frames/pairwise.py
765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 839 840 841 842 843 | |
Attributes¶
coherence
property
¶
Return raw magnitude-squared coherence values in the public shape.
Functions¶
__init__(data, sampling_rate, n_fft, window, pair_state, frequency_indices=None, source_channel_ids=None, label=None, metadata=None, channel_metadata=None, channel_ids=None, previous=None, source_time_offset=0.0, lineage=None, operation_history_prefix=())
¶
Source code in wandas/frames/pairwise.py
790 791 792 793 794 795 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 | |
wandas.frames.pairwise.CrossSpectralFrame
¶
Bases: PairwiseSpectralFrame
Complex cross-spectral values P_out_in with typed pair domains.
data is complex numeric, rank one or two, with flattened
(pair, frequency) storage internally. The required constructor state
is sampling_rate, n_fft, window, scaling, and immutable
pair_state; frequency_indices selects an ordered subset of the
n_fft // 2 + 1 bins. Each pair represents
conj(X_input) * X_output and its typed domain determines the channel
unit, reference, and level denominator. scaling is either
"spectrum" or "density"; density domains include /Hz.
magnitude, phase, and level_db are the quantity-specific
public projections, where level_db uses 10 * log10 of the
magnitude/reference ratio. Dask input remains lazy. Pair selection,
frequency slicing, metadata changes, and annotation copies preserve the
concrete type and row-matched pair state; arithmetic and A-weighting are
rejected rather than inferred from labels or history.
Source code in wandas/frames/pairwise.py
846 847 848 849 850 851 852 853 854 855 856 857 858 859 860 861 862 863 864 865 866 867 868 869 870 871 872 873 874 875 876 877 878 879 880 881 882 883 884 885 886 887 888 889 890 891 892 893 894 895 896 897 898 899 900 901 902 903 904 905 906 907 908 909 910 911 912 913 914 915 916 917 918 919 920 921 922 923 924 925 926 927 928 929 930 931 932 933 934 935 936 937 938 939 940 941 942 943 944 945 946 947 948 949 950 951 952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967 968 969 970 971 972 973 974 975 976 977 978 979 980 981 982 983 984 985 986 987 988 989 990 991 992 993 994 995 996 997 998 999 1000 1001 1002 1003 1004 1005 1006 1007 1008 1009 1010 1011 1012 1013 1014 1015 1016 1017 1018 1019 1020 1021 1022 1023 1024 1025 1026 1027 1028 1029 1030 1031 1032 1033 1034 1035 1036 1037 1038 1039 1040 1041 1042 1043 1044 1045 1046 1047 1048 1049 1050 1051 1052 1053 1054 1055 1056 1057 1058 1059 1060 1061 1062 1063 1064 1065 1066 1067 1068 1069 1070 1071 1072 1073 1074 1075 1076 1077 1078 1079 1080 1081 1082 1083 1084 1085 1086 1087 1088 1089 1090 1091 1092 1093 1094 1095 1096 1097 1098 1099 1100 1101 1102 1103 1104 1105 1106 1107 1108 1109 1110 1111 1112 1113 | |
Attributes¶
scaling
property
¶
Return the immutable CSD scaling contract.
magnitude
property
¶
Return abs(P_out_in) in the public shape.
phase
property
¶
Return angle(P_out_in) in radians in the public shape.
level_db
property
¶
Return CSD level 10 * log10(abs(P_out_in) / pair_reference).
Functions¶
__init__(data, sampling_rate, n_fft, window, pair_state, *, scaling, frequency_indices=None, source_channel_ids=None, label=None, metadata=None, channel_metadata=None, channel_ids=None, previous=None, source_time_offset=0.0, lineage=None, operation_history_prefix=())
¶
Source code in wandas/frames/pairwise.py
868 869 870 871 872 873 874 875 876 877 878 879 880 881 882 883 884 885 886 887 888 889 890 891 892 893 894 895 896 897 898 899 900 901 902 903 904 905 906 907 | |
transfer_function()
¶
Derive the canonical transfer-function matrix from this CSD without recomputing spectra.
Each row is H_out_in = P_out_in / P_in_in using the stored cross and
auto spectra, so the result equals :meth:ChannelFrame.transfer_function
with the same Welch parameters while reusing the averaging already done
for this Frame. scaling, n_fft, window, the represented
frequency bins, source identity, and per-pair source_time_offset are
preserved. The result uses canonical output-major/input-minor row order
and stays Dask-lazy.
Returns:
| Type | Description |
|---|---|
TransferFunctionFrame
|
TransferFunctionFrame with every output/input pair and |
TransferFunctionFrame
|
|
TransferFunctionFrame
|
meth: |
Raises:
| Type | Description |
|---|---|
KeyError
|
If this Frame no longer holds every output/input pair, for
example after :meth: |
Example
h = frame.csd(n_fft=4096).cache().transfer_function().select_pair(1, 0)
Source code in wandas/frames/pairwise.py
940 941 942 943 944 945 946 947 948 949 950 951 952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967 968 969 970 971 972 | |
coherence()
¶
Derive magnitude-squared coherence from this CSD without recomputing spectra.
Each row is abs(P_out_in)**2 / (P_out_out * P_in_in); bins with a
zero auto-spectrum product are NaN. Coherence does not depend on the
CSD scaling, so a mean-averaged CSD gives the same values as
:meth:ChannelFrame.coherence with matching Welch parameters. The
segment averaging (average) and detrending of this CSD are reused as
they were computed. A median-averaged CSD does not give SciPy's
coherence definition, and its values can exceed 1; deriving coherence
from a CSD whose lineage records average="median" emits a
UserWarning. n_fft,
window, represented bins, source identity, and per-pair
source_time_offset are preserved, rows use canonical
output-major/input-minor order, and the result stays Dask-lazy.
Returns:
| Type | Description |
|---|---|
CoherenceFrame
|
CoherenceFrame with every output/input pair. |
Raises:
| Type | Description |
|---|---|
KeyError
|
If this Frame no longer holds every output/input pair, for
example after :meth: |
Example
g = frame.csd(n_fft=4096).cache().coherence().select_pair(1, 0)
Source code in wandas/frames/pairwise.py
974 975 976 977 978 979 980 981 982 983 984 985 986 987 988 989 990 991 992 993 994 995 996 997 998 999 1000 1001 1002 1003 1004 1005 1006 1007 1008 1009 1010 1011 1012 1013 1014 | |
wandas.frames.pairwise.TransferFunctionFrame
¶
Bases: PairwiseSpectralFrame
Complex transfer values with truthful denominator and reference state.
data is complex numeric, rank one or two, with flattened
(pair, frequency) storage internally. The required constructor state
is sampling_rate, n_fft, window, scaling,
denominator_role, and immutable pair_state; frequency_indices
selects an ordered subset of the n_fft // 2 + 1 bins. The canonical
denominator_role="input" stores
H_out_in = P_out_in / P_in_in. denominator_role="output" exists only
to read WDF files that store the older output-denominator values and is
never relabeled as canonical.
gain, phase, gain_db, and transfer_level_db are the
quantity-specific public projections. gain_db is available only when
every selected pair is dimensionless; transfer_level_db uses each
typed output/input reference ratio. Zero-denominator non-finite values
remain observable. Dask input remains lazy, and selection, slicing,
annotation, and metadata operations preserve the concrete type and typed
rows. Arithmetic and A-weighting are rejected rather than inferred from
labels or operation history.
Source code in wandas/frames/pairwise.py
1116 1117 1118 1119 1120 1121 1122 1123 1124 1125 1126 1127 1128 1129 1130 1131 1132 1133 1134 1135 1136 1137 1138 1139 1140 1141 1142 1143 1144 1145 1146 1147 1148 1149 1150 1151 1152 1153 1154 1155 1156 1157 1158 1159 1160 1161 1162 1163 1164 1165 1166 1167 1168 1169 1170 1171 1172 1173 1174 1175 1176 1177 1178 1179 1180 1181 1182 1183 1184 1185 1186 1187 1188 1189 1190 1191 1192 1193 1194 1195 1196 1197 1198 1199 1200 1201 1202 1203 1204 1205 1206 1207 1208 1209 1210 1211 1212 1213 1214 1215 1216 1217 1218 1219 1220 1221 1222 1223 1224 1225 1226 1227 1228 1229 1230 1231 1232 1233 1234 1235 1236 1237 1238 1239 1240 1241 1242 1243 1244 1245 1246 1247 1248 1249 1250 1251 1252 1253 1254 1255 1256 1257 1258 1259 1260 1261 1262 1263 1264 1265 1266 1267 1268 1269 1270 1271 1272 1273 | |
Attributes¶
scaling
property
¶
Return the immutable transfer scaling contract.
denominator_role
property
¶
Return the immutable transfer denominator definition.
definition
property
¶
Return the persisted numerical-definition identifier.
gain
property
¶
Return linear transfer gain abs(H_out_in).
phase
property
¶
Return transfer phase angle(H_out_in) in radians.
gain_db
property
¶
Return 20 * log10(abs(H)) for dimensionless selected pairs only.
transfer_level_db
property
¶
Return transfer level relative to each typed output/input reference ratio.
Functions¶
__init__(data, sampling_rate, n_fft, window, pair_state, *, scaling, denominator_role='input', frequency_indices=None, source_channel_ids=None, label=None, metadata=None, channel_metadata=None, channel_ids=None, previous=None, source_time_offset=0.0, lineage=None, operation_history_prefix=())
¶
Source code in wandas/frames/pairwise.py
1141 1142 1143 1144 1145 1146 1147 1148 1149 1150 1151 1152 1153 1154 1155 1156 1157 1158 1159 1160 1161 1162 1163 1164 1165 1166 1167 1168 1169 1170 1171 1172 1173 1174 1175 1176 1177 1178 1179 1180 1181 1182 1183 1184 | |
wandas.frames.spectrogram.SpectrogramFrame
¶
Bases: SpectralPropertiesMixin, BaseFrame[NDArrayComplex]
Class for handling time-frequency domain data (spectrograms).
This class represents spectrogram data obtained through Short-Time Fourier Transform (STFT) or similar time-frequency analysis methods. It provides methods for visualization, manipulation, and conversion back to time domain.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
data
|
Array
|
DaArray. The spectrogram data. Must be a dask array with shape: - (channels, frequency_bins, time_frames) for multi-channel data - (frequency_bins, time_frames) for single-channel data, which will be reshaped to (1, frequency_bins, time_frames) |
required |
sampling_rate
|
float
|
float. The sampling rate of the original time-domain signal in Hz. |
required |
n_fft
|
int
|
int. The FFT size used to generate this spectrogram. The frequency dimension must
contain exactly |
required |
hop_length
|
int
|
int. Number of samples between successive frames. |
required |
allow_sparse
|
bool
|
bool, default=False. Permit hops larger than the window for display sampling. |
False
|
win_length
|
int | None
|
int, optional. The window length in samples. If None, defaults to n_fft. |
None
|
window
|
str
|
str, default="hann". The window function to use (e.g., "hann", "hamming", "blackman"). |
'hann'
|
label
|
str | None
|
str, optional. A label for the frame. |
None
|
metadata
|
dict[str, Any] | None
|
dict, optional. Additional metadata for the frame. |
None
|
lineage
|
Any | None
|
LineageNode, optional. Constructor override for the runtime lineage. When omitted, a source node is
created. |
None
|
channel_metadata
|
Sequence[ChannelMetadata | dict[str, Any]] | None
|
list[ChannelMetadata], optional. Metadata for each channel in the frame. |
None
|
previous
|
BaseFrame[Any] | None
|
BaseFrame, optional. Immediate receiver Frame for process-local data comparison. For multi-input operations, follows only the left/base receiver. Not persisted in WDF. |
None
|
frame_time_origin
|
float | None
|
float, optional. First physical frame center relative to input start, in seconds. None explicitly means unknown; STFT supplies its padding-dependent origin. Existing local time axes are unchanged. |
None
|
Attributes:
| Name | Type | Description |
|---|---|---|
magnitude |
NDArrayReal
|
NDArrayReal. The magnitude spectrogram. |
phase |
NDArrayReal
|
NDArrayReal. The phase spectrogram in radians. |
power |
NDArrayReal
|
NDArrayReal. The power spectrogram. |
dB |
NDArrayReal
|
NDArrayReal. The spectrogram in decibels relative to channel reference values. |
dBA |
NDArrayReal
|
NDArrayReal. The A-weighted spectrogram in decibels. |
n_frames |
int
|
int. Number of time frames. |
n_freq_bins |
int
|
int. Number of frequency bins. |
freqs |
NDArrayReal
|
NDArrayReal. The frequency axis values in Hz. |
times |
NDArrayReal
|
NDArrayReal. The time axis values in seconds. |
Examples:
Create a spectrogram from a time-domain signal:
>>> signal = ChannelFrame.from_wav("audio.wav")
>>> spectrogram = signal.stft(n_fft=2048, hop_length=512)
Extract a specific time frame:
>>> frame_at_1s = spectrogram.get_frame_at(int(1.0 * sampling_rate / hop_length))
Convert back to time domain:
>>> reconstructed = spectrogram.to_channel_frame()
Plot the spectrogram:
>>> spectrogram.plot()
Source code in wandas/frames/spectrogram.py
29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 | |
Attributes¶
n_fft
property
¶
Return the read-only FFT size defining the frequency axis.
hop_length
property
¶
Return the read-only hop length defining the time spacing.
win_length
property
¶
Return the read-only analysis window length.
window
property
¶
Return the read-only analysis window; rerun STFT to change it.
n_frames
property
¶
Get the number of time frames.
Returns:
| Name | Type | Description |
|---|---|---|
int |
int
|
The number of time frames in the spectrogram. |
n_freq_bins
property
¶
Get the number of frequency bins.
Returns:
| Name | Type | Description |
|---|---|---|
int |
int
|
The number of frequency bins (n_fft // 2 + 1). |
freqs
property
¶
Get the frequency axis values in Hz.
Values are derived on access from sampling_rate and n_fft using the
canonical one-sided real-FFT grid.
Returns:
| Name | Type | Description |
|---|---|---|
NDArrayReal |
NDArrayReal
|
Array of frequency values corresponding to each frequency bin. |
times
property
¶
Get the time axis values in seconds.
This is a zero-based local axis derived from hop_length and
sampling_rate. Absolute placement belongs to source_time_offset.
Returns:
| Name | Type | Description |
|---|---|---|
NDArrayReal |
NDArrayReal
|
Array of time values corresponding to each time frame. |
frame_time_origin
property
¶
First physical frame center relative to input start; None means unknown.
This read-only constructor state includes the STFT padding displacement. It is independent of channel-specific source offsets and is retained through same-domain transforms and contiguous time slicing.
frame_center_times
property
¶
Return actual frame centers on every channel's source timeline, in seconds.
Returns:
| Type | Description |
|---|---|
NDArrayReal
|
A new array shaped |
NDArrayReal
|
padding centers. Only axis metadata is computed; PCM and STFT values |
NDArrayReal
|
remain lazy. Legacy |
NDArrayReal
|
Binary arithmetic retains the left frame's time placement, as it |
NDArrayReal
|
does for |
Raises:
| Type | Description |
|---|---|
ValueError
|
The physical origin is unknown, for example in an old
WDF or a manually constructed frame without |
Examples:
>>> spectrum = recording.stft(n_fft=2048, hop_length=512)
>>> centers = spectrum.frame_center_times[0]
source_times
property
¶
Get frame times relative to the original source timeline.
Functions¶
__init__(data, sampling_rate, n_fft, hop_length, win_length=None, window='hann', label=None, metadata=None, lineage=None, channel_metadata=None, channel_ids=None, previous=None, source_time_offset=0.0, operation_history_prefix=(), *, frame_time_origin=None, allow_sparse=False)
¶
Initialize a complete canonical one-sided spectrogram.
See the class docstring for parameter descriptions. Frequency and local time axes are derived from the analysis parameters rather than stored as mutable coordinate state.
Source code in wandas/frames/spectrogram.py
92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 | |
plot(plot_type='spectrogram', ax=None, title=None, cmap='jet', vmin=None, vmax=None, fmin=0, fmax=None, xlim=None, ylim=None, Aw=False, overlay=False, **kwargs)
¶
Plot the spectrogram using various visualization strategies.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
plot_type
|
str
|
str, default="spectrogram". Type of plot to create. |
'spectrogram'
|
ax
|
Axes | None
|
matplotlib.axes.Axes, optional. Axes to plot on. If None, creates new axes. |
None
|
title
|
str | None
|
str, optional. Title for the plot. If None, uses the frame label. |
None
|
cmap
|
str
|
str, default="jet". Colormap name for the spectrogram visualization. |
'jet'
|
vmin
|
float | None
|
float, optional. Minimum value for colormap scaling (dB). Auto-calculated if None. |
None
|
vmax
|
float | None
|
float, optional. Maximum value for colormap scaling (dB). Auto-calculated if None. |
None
|
fmin
|
float
|
float, default=0. Minimum frequency to display (Hz). |
0
|
fmax
|
float | None
|
float, optional. Maximum frequency to display (Hz). If None, uses Nyquist frequency. |
None
|
xlim
|
tuple[float, float] | None
|
tuple[float, float], optional. Time axis limits as (start_time, end_time) in seconds. |
None
|
ylim
|
tuple[float, float] | None
|
tuple[float, float], optional. Frequency axis limits as (min_freq, max_freq) in Hz. |
None
|
Aw
|
bool
|
bool, default=False. Whether to apply A-weighting to the spectrogram. |
False
|
overlay
|
bool
|
bool, default=False. Whether to overlay channels on a single axes. |
False
|
**kwargs
|
Any
|
dict. Additional keyword arguments passed to Matplotlib plotting methods. |
{}
|
Returns:
| Type | Description |
|---|---|
Axes | Iterator[Axes]
|
Union[Axes, Iterator[Axes]]: The matplotlib axes containing the plot, or an iterator of axes for multi-plot outputs. |
Examples:
>>> stft = cf.stft()
>>> # Basic spectrogram
>>> stft.plot()
>>> # Custom color scale and frequency range
>>> stft.plot(vmin=-80, vmax=-20, fmin=100, fmax=5000)
>>> # A-weighted spectrogram
>>> stft.plot(Aw=True, cmap="viridis")
Source code in wandas/frames/spectrogram.py
284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 | |
plot_Aw(plot_type='spectrogram', ax=None, **kwargs)
¶
Plot the A-weighted spectrogram.
A convenience method that calls plot() with Aw=True, applying A-weighting to the spectrogram before plotting.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
plot_type
|
str
|
str, default="spectrogram". Type of plot to create. |
'spectrogram'
|
ax
|
Axes | None
|
matplotlib.axes.Axes, optional. Axes to plot on. If None, creates new axes. |
None
|
**kwargs
|
Any
|
dict. Additional keyword arguments passed to plot(). Accepts all parameters from plot() except Aw (which is set to True). |
{}
|
Returns:
| Type | Description |
|---|---|
Axes | Iterator[Axes]
|
Union[Axes, Iterator[Axes]]: The matplotlib axes containing the plot. |
Examples:
>>> stft = cf.stft()
>>> # A-weighted spectrogram with custom settings
>>> stft.plot_Aw(vmin=-60, vmax=-10, cmap="magma")
Source code in wandas/frames/spectrogram.py
361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 | |
cepstrum(floor=1e-12)
¶
Calculate a real cepstrum independently at every time frame.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
floor
|
float
|
float, default=1e-12. Positive finite floor applied to normalized STFT magnitude before taking the logarithm. |
1e-12
|
Returns:
| Name | Type | Description |
|---|---|---|
CepstrogramFrame |
CepstrogramFrame
|
New lazy coefficients shaped |
Raises:
| Type | Description |
|---|---|
TypeError
|
If |
ValueError
|
If |
Notes
The source SpectrogramFrame already contains normalized one-sided
STFT amplitudes. This method computes
irfft(log(max(abs(stft), floor))) along its frequency axis without
recomputing the time-domain STFT. It only builds a Dask graph.
Examples:
>>> cepstrogram = frame.stft(n_fft=2048).cepstrum()
>>> envelope = cepstrogram.lifter(0.002).to_spectral_envelope()
Source code in wandas/frames/spectrogram.py
389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 | |
abs()
¶
Compute the absolute value (magnitude) of the complex spectrogram.
This method calculates the magnitude of each complex value in the spectrogram, converting the complex-valued data to real-valued magnitude data. The result remains a SpectrogramFrame but carries a real numeric dtype.
Returns:
| Name | Type | Description |
|---|---|---|
SpectrogramFrame |
SpectrogramFrame
|
A new SpectrogramFrame containing real-valued magnitudes. |
Examples:
>>> signal = ChannelFrame.from_wav("audio.wav")
>>> spectrogram = signal.stft(n_fft=2048, hop_length=512)
>>> magnitude_spectrogram = spectrogram.abs()
>>> # The magnitude can be accessed via the magnitude property or data
>>> print(magnitude_spectrogram.magnitude.shape)
Source code in wandas/frames/spectrogram.py
446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 | |
get_frame_at(time_idx)
¶
Extract spectral data at a specific time frame.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
time_idx
|
int
|
int. Index of the time frame to extract. |
required |
Returns:
| Name | Type | Description |
|---|---|---|
SpectralFrame |
SpectralFrame
|
A new SpectralFrame containing the spectral data at the specified time. |
Raises:
| Type | Description |
|---|---|
IndexError
|
If time_idx is out of range. |
Source code in wandas/frames/spectrogram.py
482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 | |
to_channel_frame()
¶
Convert the spectrogram back to time domain using inverse STFT.
This method performs an inverse Short-Time Fourier Transform (ISTFT) to reconstruct the time-domain signal from the spectrogram.
Returns:
| Name | Type | Description |
|---|---|---|
ChannelFrame |
ChannelFrame
|
A new ChannelFrame containing the reconstructed time-domain signal. |
See Also
istft : Alias for this method with more intuitive naming.
Source code in wandas/frames/spectrogram.py
524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 | |
istft()
¶
Convert the spectrogram back to time domain using inverse STFT.
This is an alias for to_channel_frame() with a more intuitive name.
It performs an inverse Short-Time Fourier Transform (ISTFT) to
reconstruct the time-domain signal from the spectrogram.
Returns:
| Name | Type | Description |
|---|---|---|
ChannelFrame |
ChannelFrame
|
A new ChannelFrame containing the reconstructed time-domain signal. |
See Also
to_channel_frame : The underlying implementation.
Examples:
>>> signal = ChannelFrame.from_wav("audio.wav")
>>> spectrogram = signal.stft(n_fft=2048, hop_length=512)
>>> reconstructed = spectrogram.istft()
Source code in wandas/frames/spectrogram.py
572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 | |
csd(scaling='spectrum')
¶
Calculate the cross-spectral matrix averaged over every time frame.
The result is the same :class:CrossSpectralFrame that
:meth:ChannelFrame.csd returns, with every output/input pair in
canonical order, so transfer_function(), coherence() and
select_pair() work on it unchanged. Reusing an STFT this way avoids a
second Welch pass over the signal.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
scaling
|
Literal['spectrum', 'density']
|
|
'spectrum'
|
Returns:
| Type | Description |
|---|---|
CrossSpectralFrame
|
CrossSpectralFrame whose rows store |
CrossSpectralFrame
|
converted to SciPy CSD scaling for this frame's window. |
Raises:
| Type | Description |
|---|---|
ValueError
|
If |
Notes
Frames are mean-averaged. Detrending happens in the STFT, so the
result equals ChannelFrame.csd(detrend=...) over the same
segments when the STFT used the same detrend. The default
stft() does not detrend, so it matches csd(detrend=False),
not the csd() default of "constant". The STFT pads the
signal, so its first and last frames include zeros; slice the time
axis (spec[:, :, a:b]) to choose which frames are averaged. Each
pair's source_time_offset is the start of the
first averaged window, as in :meth:ChannelFrame.csd; it stays the
spectrogram's offset when frame_time_origin is unknown.
Construction stays lazy.
Examples:
>>> spec = frame.stft(n_fft=4096, hop_length=1024).cache()
>>> h = spec[:, :, 4:-4].csd().transfer_function().select_pair(1, 0)
Source code in wandas/frames/spectrogram.py
593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 | |
to_dataframe()
¶
DataFrame conversion is not supported for SpectrogramFrame.
SpectrogramFrame contains 3D data (channels, frequency_bins, time_frames) which cannot be directly converted to a 2D DataFrame. Consider using get_frame_at() to extract a specific time frame as a SpectralFrame, then convert that to a DataFrame.
Raises:
| Type | Description |
|---|---|
NotImplementedError
|
Always raised as DataFrame conversion is not supported. |
Source code in wandas/frames/spectrogram.py
681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 | |
info()
¶
Display comprehensive information about the SpectrogramFrame.
This method prints a summary of the frame's properties including: - Number of channels - Sampling rate - FFT size - Hop length - Window length - Window function - Frequency range - Number of frequency bins - Frequency resolution (ΔF) - Number of time frames - Time resolution (ΔT) - Total duration - Channel labels - Number of operations applied
This is a convenience method to view all key properties at once, similar to pandas DataFrame.info().
Examples:
>>> signal = ChannelFrame.from_wav("audio.wav")
>>> spectrogram = signal.stft(n_fft=2048, hop_length=512)
>>> spectrogram.info()
SpectrogramFrame Information:
Channels: 2
Sampling rate: 44100 Hz
FFT size: 2048
Hop length: 512 samples
Window length: 2048 samples
Window: hann
Frequency range: 0.0 - 22050.0 Hz
Frequency bins: 1025
Frequency resolution (ΔF): 21.5 Hz
Time frames: 100
Time resolution (ΔT): 11.6 ms
Total duration: 1.16 s
Channel labels: ['ch0', 'ch1']
Operations Applied: 1
Source code in wandas/frames/spectrogram.py
698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755 756 757 758 759 | |
from_numpy(data, sampling_rate, n_fft, hop_length, win_length=None, window='hann', label=None, metadata=None, lineage=None, channel_metadata=None, channel_ids=None, previous=None, *, frame_time_origin=None)
classmethod
¶
Create a SpectrogramFrame from a NumPy array.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
data
|
NDArrayComplex
|
NumPy array containing spectrogram data. Shape should be (n_channels, n_freq_bins, n_time_frames) or (n_freq_bins, n_time_frames) for single channel. |
required |
sampling_rate
|
float
|
The sampling rate in Hz. |
required |
n_fft
|
int
|
The FFT size used to generate this spectrogram. |
required |
hop_length
|
int
|
Number of samples between successive frames. |
required |
win_length
|
int | None
|
The window length in samples. If None, defaults to n_fft. |
None
|
window
|
str
|
The window function used (e.g., "hann", "hamming"). |
'hann'
|
label
|
str | None
|
A label for the frame. |
None
|
metadata
|
dict[str, Any] | None
|
Optional metadata dictionary. |
None
|
frame_time_origin
|
float | None
|
Known first physical frame center in seconds, including STFT padding displacement. None leaves physical centers unknown. |
None
|
lineage
|
Any | None
|
Runtime operation lineage for this frame. |
None
|
channel_metadata
|
Sequence[ChannelMetadata | dict[str, Any]] | None
|
Metadata for each channel. |
None
|
previous
|
BaseFrame[Any] | None
|
Immediate receiver Frame for process-local data comparison. For multi-input operations, follows only the left/base receiver. Not persisted in WDF. |
None
|
Returns:
| Type | Description |
|---|---|
SpectrogramFrame
|
A new SpectrogramFrame containing the NumPy data. |
Source code in wandas/frames/spectrogram.py
761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 | |
wandas.frames.cepstral.CepstralFrame
¶
Bases: BaseFrame[NDArrayReal]
Immutable, lazy real-cepstrum data on a quefrency axis.
Data is rank two with dimensions (channel, quefrency) and a real dtype.
n_fft is the circular period of the complete cepstrum; sliced frames keep
that period so methods can reject incomplete axes. Quefrency is measured in
seconds at spacing 1 / sampling_rate. The sampling rate is immutable
because changing it would reinterpret the stored axis.
lifter() preserves this frame family. to_spectral_envelope() returns a
:class:~wandas.frames.spectral.SpectralFrame. Both operations remain Dask
backed and preserve channel identity, metadata, source-time offsets, and
semantic lineage. Arithmetic between two CepstralFrames requires the same
n_fft, window, win_length and quefrency coordinates.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
data
|
Array
|
dask.array.Array. Real coefficients shaped |
required |
sampling_rate
|
float
|
float. Sampling rate in Hz that defines the quefrency-bin spacing. |
required |
n_fft
|
int
|
int. Positive FFT size of the complete cepstrum. Sliced data may contain fewer bins but never more than this value. |
required |
window
|
str
|
str, default="hann". Window used by the originating cepstrum analysis. |
'hann'
|
label
|
str | None
|
str, optional. Human-readable frame label. |
None
|
metadata
|
dict[str, Any] | None
|
dict, optional. User and recording metadata, copied on construction. |
None
|
channel_metadata
|
Sequence[ChannelMetadata | dict[str, Any]] | None
|
sequence, optional. Metadata aligned with the channel axis. |
None
|
channel_ids
|
list[str] | None
|
list[str], optional. Stable identifiers aligned with the channel axis. |
None
|
previous
|
BaseFrame[Any] | None
|
BaseFrame, optional. Immediate receiver Frame for process-local data comparison. For multi-input operations, follows only the left/base receiver. Not persisted in WDF. |
None
|
source_time_offset
|
float | Sequence[float] | NDArrayReal
|
float or sequence, default=0.0. Per-channel source timeline offsets, preserved across domain changes. |
0.0
|
lineage
|
LineageNode | None
|
LineageNode, optional. Authoritative runtime semantic lineage. |
None
|
operation_history_prefix
|
Sequence[Mapping[str, Any]]
|
sequence, default=(). Persisted display history for a new source frame. |
()
|
win_length
|
int | None
|
int, optional. Length of the originating analysis window, |
None
|
Raises:
| Type | Description |
|---|---|
TypeError
|
If coefficients are complex, |
ValueError
|
If rank, FFT size, or coefficient count violates the domain contract. |
Examples:
>>> cepstrum = frame.cepstrum(n_fft=2048)
>>> envelope = cepstrum.lifter(0.002).to_spectral_envelope()
Source code in wandas/frames/cepstral.py
27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 | |
Attributes¶
n_fft
property
¶
Return the immutable circular period of the complete cepstrum.
window
property
¶
Return the immutable originating analysis-window name.
win_length
property
¶
Return the length of the originating analysis window in samples.
sampling_rate
property
¶
Return the immutable rate that defines the quefrency spacing.
quefrencies
property
¶
Return a defensive copy of the represented quefrency bins in seconds.
Functions¶
__init__(data, sampling_rate, n_fft, window='hann', label=None, metadata=None, channel_metadata=None, channel_ids=None, previous=None, source_time_offset=0.0, lineage=None, operation_history_prefix=(), win_length=None)
¶
Source code in wandas/frames/cepstral.py
77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 | |
lifter(cutoff, mode='low')
¶
Keep low- or high-quefrency coefficients.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
cutoff
|
float
|
float. Positive quefrency boundary in seconds. It must reach at least one represented bin and remain below half the complete cepstrum. |
required |
mode
|
Literal['low', 'high']
|
{"low", "high"}, default="low". |
'low'
|
Returns:
| Name | Type | Description |
|---|---|---|
CepstralFrame |
CepstralFrame
|
A new lazy frame with the same axes and metadata. |
Raises:
| Type | Description |
|---|---|
ValueError
|
If this frame has been sliced on the quefrency axis or the lifter parameters cannot be represented. |
Notes
The method builds a Dask graph and does not compute coefficients.
Examples:
>>> smooth = frame.cepstrum().lifter(0.002, mode="low")
Source code in wandas/frames/cepstral.py
241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 | |
to_spectral_envelope()
¶
Convert a complete real cepstrum to a smooth spectral envelope.
Returns:
| Name | Type | Description |
|---|---|---|
SpectralFrame |
SpectralFrame
|
New lazy complex-valued frequency data with zero phase, the original
|
Raises:
| Type | Description |
|---|---|
ValueError
|
If this frame has been sliced on the quefrency axis. Asymmetric concrete coefficients raise when the lazy result is computed. |
Notes
This method builds a Dask graph. It does not compute the envelope.
Examples:
>>> envelope = frame.cepstrum().lifter(0.002).to_spectral_envelope()
Source code in wandas/frames/cepstral.py
274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 | |
plot(plot_type='quefrency', ax=None, *, title=None, xlabel='Quefrency [s]', ylabel='Real cepstrum', **kwargs)
¶
Plot real coefficients against quefrency.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
plot_type
|
str
|
str, default="quefrency". Only |
'quefrency'
|
ax
|
Axes | None
|
matplotlib.axes.Axes, optional. Existing axes. A new figure and axes are created when omitted. |
None
|
title
|
str | None
|
str, optional. Plot title; defaults to the frame label. |
None
|
xlabel
|
str
|
str. Horizontal axis label. |
'Quefrency [s]'
|
ylabel
|
str
|
str. Vertical axis label. |
'Real cepstrum'
|
**kwargs
|
Any
|
Any. Keyword arguments passed to |
{}
|
Returns:
| Type | Description |
|---|---|
Axes | Iterator[Axes]
|
matplotlib.axes.Axes: Axes containing one line per channel. |
Raises:
| Type | Description |
|---|---|
ValueError
|
If another plot type is requested. |
Notes
Plotting is an explicit compute boundary and materializes coefficients.
Examples:
>>> cepstrum.plot()
Source code in wandas/frames/cepstral.py
333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 | |
wandas.frames.cepstrogram.CepstrogramFrame
¶
Bases: BaseFrame[NDArrayReal]
Immutable, lazy real cepstrum evolving over STFT time frames.
Data is rank three with dimensions (channel, quefrency, time).
n_fft defines the complete circular quefrency axis, while
hop_length defines the spacing of the retained STFT time frames.
lifter() preserves this frame family and
to_spectral_envelope() returns a
:class:~wandas.frames.spectrogram.SpectrogramFrame.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
data
|
Array
|
dask.array.Array. Real coefficients shaped |
required |
sampling_rate
|
float
|
float. Sampling rate in Hz defining both axis spacings. |
required |
n_fft
|
int
|
int. Positive FFT size of the complete cepstrum. |
required |
hop_length
|
int
|
int. Positive sample distance between adjacent time frames. |
required |
allow_sparse
|
bool
|
bool. Permit gaps between windows; default False. Other values raise TypeError. Sparse conversion preserves time axes. |
False
|
win_length
|
int | None
|
int, optional. Analysis-window length inherited from the source spectrogram. Defaults
to |
None
|
window
|
str
|
str, default="hann". Analysis-window name inherited from the source spectrogram. |
'hann'
|
label
|
str | None
|
str, optional. Human-readable frame label. |
None
|
metadata
|
dict[str, Any] | None
|
dict, optional. User and recording metadata, copied on construction. |
None
|
channel_metadata
|
Sequence[ChannelMetadata | dict[str, Any]] | None
|
sequence, optional. Metadata aligned with the channel axis. |
None
|
channel_ids
|
list[str] | None
|
list[str], optional. Stable identifiers aligned with the channel axis. |
None
|
previous
|
BaseFrame[Any] | None
|
BaseFrame, optional. Immediate receiver Frame for process-local data comparison. For multi-input operations, follows only the left/base receiver. Not persisted in WDF. |
None
|
source_time_offset
|
float | Sequence[float] | NDArrayReal
|
float or sequence, default=0.0. Per-channel source timeline offsets. |
0.0
|
frame_time_origin
|
float | None
|
float, optional. First physical frame center relative to input start, in seconds. None explicitly means unknown; STFT supplies its padding-dependent origin. Existing local time axes are unchanged. |
None
|
lineage
|
LineageNode | None
|
LineageNode, optional. Authoritative runtime semantic lineage. |
None
|
operation_history_prefix
|
Sequence[Mapping[str, Any]]
|
sequence, default=(). Persisted display history for a new source frame. |
()
|
Raises:
| Type | Description |
|---|---|
TypeError
|
If coefficients are complex or domain parameters have invalid types. |
ValueError
|
If rank, FFT size, coefficient count, or time-analysis parameters are invalid. |
Examples:
>>> cepstrogram = frame.stft(n_fft=2048).cepstrum()
>>> envelope = cepstrogram.lifter(0.002).to_spectral_envelope()
Source code in wandas/frames/cepstrogram.py
29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 | |
Attributes¶
n_fft
property
¶
Return the immutable circular period of the cepstrum.
hop_length
property
¶
Return the immutable sample spacing between time frames.
win_length
property
¶
Return the immutable originating analysis-window length.
window
property
¶
Return the immutable originating analysis-window name.
sampling_rate
property
¶
Return the immutable rate defining quefrency and time spacing.
n_quefrency_bins
property
¶
Return the represented quefrency-bin count.
n_frames
property
¶
Return the time-frame count.
quefrencies
property
¶
Return a defensive copy of represented quefrencies in seconds.
times
property
¶
Return a defensive copy of frame times in seconds.
frame_time_origin
property
¶
First physical frame center relative to input start; None means unknown.
This read-only constructor state includes the STFT padding displacement. It is independent of channel-specific source offsets and is retained through same-domain transforms and contiguous time slicing.
frame_center_times
property
¶
Return actual frame centers on every channel's source timeline, in seconds.
Returns:
| Type | Description |
|---|---|
NDArrayReal
|
A new array shaped |
NDArrayReal
|
padding centers. Only axis metadata is computed; PCM and STFT values |
NDArrayReal
|
remain lazy. Legacy |
NDArrayReal
|
Binary arithmetic retains the left frame's time placement, as it |
NDArrayReal
|
does for |
Raises:
| Type | Description |
|---|---|
ValueError
|
The physical origin is unknown, for example in an old
WDF or a manually constructed frame without |
Examples:
>>> spectrum = recording.stft(n_fft=2048, hop_length=512)
>>> centers = spectrum.frame_center_times[0]
source_times
property
¶
Return frame times on each channel's original source timeline.
Functions¶
__init__(data, sampling_rate, n_fft, hop_length, win_length=None, window='hann', label=None, metadata=None, channel_metadata=None, channel_ids=None, previous=None, source_time_offset=0.0, lineage=None, operation_history_prefix=(), *, frame_time_origin=None, allow_sparse=False)
¶
Source code in wandas/frames/cepstrogram.py
76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 | |
lifter(cutoff, mode='low')
¶
Keep low- or high-quefrency coefficients at every time frame.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
cutoff
|
float
|
float. Positive quefrency boundary in seconds. It must reach at least one bin and remain below half of the complete cepstrum. |
required |
mode
|
Literal['low', 'high']
|
{"low", "high"}, default="low". |
'low'
|
Returns:
| Name | Type | Description |
|---|---|---|
CepstrogramFrame |
CepstrogramFrame
|
New lazy coefficients with unchanged time and channel axes. |
Raises:
| Type | Description |
|---|---|
ValueError
|
If the quefrency axis was sliced or the cutoff is not representable. |
Notes
This method only builds a Dask graph.
Source code in wandas/frames/cepstrogram.py
343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 | |
to_spectral_envelope()
¶
Reconstruct a smooth magnitude spectrogram with zero phase.
Returns:
| Name | Type | Description |
|---|---|---|
SpectrogramFrame |
SpectrogramFrame
|
New lazy frequency-time data preserving the original STFT analysis state, channels, metadata, and source-time offsets. |
Raises:
| Type | Description |
|---|---|
ValueError
|
If the quefrency axis was sliced. Asymmetric concrete coefficients raise when the lazy result is computed. |
Notes
This method only builds a Dask graph.
Source code in wandas/frames/cepstrogram.py
377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 | |
to_dataframe()
¶
Reject conversion because the frame has three semantic dimensions.
Raises:
| Type | Description |
|---|---|
NotImplementedError
|
Always raised. Materialize selected data when a tabular representation is required. |
Source code in wandas/frames/cepstrogram.py
447 448 449 450 451 452 453 454 | |
plot(plot_type='cepstrogram', ax=None, *, title=None, xlabel='Time [s]', ylabel='Quefrency [s]', cmap='RdBu_r', qmin=0.0, qmax=None, vmin=None, vmax=None, **kwargs)
¶
Plot real coefficients over time and quefrency.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
plot_type
|
str
|
str, default="cepstrogram". Only |
'cepstrogram'
|
ax
|
Axes | None
|
matplotlib.axes.Axes, optional. Existing axes for a single-channel frame. Multi-channel frames create one axes per channel when omitted. |
None
|
title
|
str | None
|
str, optional. Plot title prefix; defaults to the frame label. |
None
|
xlabel
|
str
|
str. Horizontal axis label. |
'Time [s]'
|
ylabel
|
str
|
str. Vertical axis label. |
'Quefrency [s]'
|
cmap
|
str
|
str, default="RdBu_r". Matplotlib colormap for signed real coefficients. |
'RdBu_r'
|
qmin
|
float
|
float, optional. Lower display bound on the quefrency axis in seconds. |
0.0
|
qmax
|
float | None
|
float, optional. Upper display bound on the quefrency axis in seconds. |
None
|
vmin
|
float | None
|
float, optional. Lower shared color limit. When omitted, a symmetric robust range is estimated from the displayed coefficients except the dominant zero-quefrency row. |
None
|
vmax
|
float | None
|
float, optional. Upper shared color limit. |
None
|
**kwargs
|
Any
|
Any. Additional keyword arguments passed to |
{}
|
Returns:
| Type | Description |
|---|---|
Axes | Iterator[Axes]
|
matplotlib.axes.Axes or Iterator[matplotlib.axes.Axes]: One axes for mono data or an iterator for multiple channels. |
Notes
Plotting is an explicit compute boundary.
Source code in wandas/frames/cepstrogram.py
456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 | |
wandas.frames.noct.NOctFrame
¶
Bases: BaseFrame[NDArrayReal]
Class for handling N-octave band analysis data.
This class represents frequency data analyzed in fractional octave bands, typically used in acoustic and vibration analysis. It handles real-valued data representing RMS amplitude in each frequency band, following standard acoustical band definitions. Values retain the input channel's physical unit.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
data
|
Array
|
DaArray. The N-octave band data. Must be a dask array with shape: - (channels, frequency_bins) for multi-channel data - (frequency_bins,) for single-channel data, which will be reshaped to (1, frequency_bins) |
required |
sampling_rate
|
float
|
float. The sampling rate of the original time-domain signal in Hz. |
required |
fmin
|
float
|
float, default=0. Lower frequency bound in Hz. |
0
|
fmax
|
float
|
float, default=0. Upper frequency bound in Hz. |
0
|
n
|
int
|
int, default=3. Number of bands per octave (e.g., 3 for third-octave bands). |
3
|
G
|
int
|
int, default=10. Exact center-frequency ratio convention: 10 selects base
|
10
|
fr
|
int
|
int, default=1000. Reference frequency in Hz, typically 1000 Hz for acoustic analysis. |
1000
|
label
|
str | None
|
str, optional. A label for the frame. |
None
|
metadata
|
dict[str, Any] | None
|
dict, optional. Additional metadata for the frame. |
None
|
lineage
|
Any | None
|
LineageNode, optional. Constructor override for the runtime lineage. When omitted, a source node is
created. |
None
|
channel_metadata
|
Sequence[ChannelMetadata | dict[str, Any]] | None
|
list[ChannelMetadata], optional. Metadata for each channel in the frame. |
None
|
previous
|
BaseFrame[Any] | None
|
BaseFrame, optional. Immediate receiver Frame for process-local data comparison. For multi-input operations, follows only the left/base receiver. Not persisted in WDF. |
None
|
Attributes:
| Name | Type | Description |
|---|---|---|
freqs |
NDArrayReal
|
NDArrayReal. The center frequencies of each band in Hz, calculated according to the standard fractional octave band definitions. |
dB |
NDArrayReal
|
NDArrayReal. Band RMS amplitude level,
|
dBA |
NDArrayReal
|
NDArrayReal. The A-weighted spectrum in decibels, applying frequency weighting for better correlation with perceived loudness. |
fmin |
float
|
float. Lower frequency bound in Hz. |
fmax |
float
|
float. Upper frequency bound in Hz. |
n |
int
|
int. Number of bands per octave. |
G |
int
|
int. Exact center-frequency ratio convention. |
fr |
int
|
int. Reference frequency in Hz. |
Examples:
Create an N-octave band spectrum from a time-domain signal:
>>> signal = ChannelFrame.from_wav("audio.wav")
>>> spectrum = signal.noct_spectrum(fmin=20, fmax=20000, n=3)
Plot the N-octave band spectrum:
>>> spectrum.plot()
Plot with A-weighting applied:
>>> spectrum.plot(Aw=True)
Notes
- Binary operations (addition, multiplication, etc.) are not currently
supported for N-octave band data. - The actual frequency bands are determined by the parameters n, G, and fr according to IEC 61260-1:2014 standard for fractional octave band filters. - The class follows acoustic standards for band definitions and analysis, making it suitable for noise measurements and sound level analysis. - A-weighting is available for better correlation with human hearing perception, following IEC 61672-1:2013.
Source code in wandas/frames/noct.py
33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 | |
Attributes¶
n = n
instance-attribute
¶
G = G
instance-attribute
¶
fr = fr
instance-attribute
¶
fmin = fmin
instance-attribute
¶
fmax = fmax
instance-attribute
¶
dB
property
¶
Get band RMS amplitude levels relative to each channel reference.
The conversion is
20 * log10(max(band_rms / channel_ref, 1e-12)). The reference has
the same physical unit as the band RMS value and is specified by each
channel's calibration metadata.
Returns:
| Name | Type | Description |
|---|---|---|
NDArrayReal |
NDArrayReal
|
The spectrum in decibels. Shape matches the input data shape: (channels, frequency_bins). |
dBA
property
¶
Get the A-weighted spectrum in decibels.
A-weighting applies a frequency-dependent weighting filter that approximates the human ear's response to different frequencies. This is particularly useful for analyzing noise and acoustic measurements as it provides a better correlation with perceived loudness.
The weighting is applied according to IEC 61672-1:2013 standard.
Returns:
| Name | Type | Description |
|---|---|---|
NDArrayReal |
NDArrayReal
|
The A-weighted spectrum in decibels. Shape matches the input data shape: (channels, frequency_bins). |
freqs
property
¶
Get the center frequencies of each band in Hz.
These frequencies are calculated based on the N-octave band parameters (n, G, fr) and the frequency bounds (fmin, fmax) according to IEC 61260-1:2014 standard for fractional octave band filters.
Bands in MoSQITo's nominal table (25 Hz to 20 kHz for 1/3 octave,
31.5 Hz to 16 kHz for octave) are labelled with their nominal
frequency, such as 125 Hz. Bands outside it are labelled with their
exact center frequency, such as 19.95 Hz. A frame from
noct_synthesis omits the top bands MoSQITo drops above the Nyquist
frequency.
Returns:
| Name | Type | Description |
|---|---|---|
NDArrayReal |
NDArrayReal
|
Array of center frequencies for each frequency band. |
Raises:
| Type | Description |
|---|---|
ValueError
|
If the center frequencies cannot be calculated or the result is not a numpy array. |
Functions¶
__init__(data, sampling_rate, fmin=0, fmax=0, n=3, G=10, fr=1000, label=None, metadata=None, channel_metadata=None, channel_ids=None, previous=None, source_time_offset=0.0, lineage=None, operation_history_prefix=())
¶
Initialize a NOctFrame instance.
Sets up N-octave band analysis parameters and prepares the frame for storing band-filtered data. Data shape is validated to ensure compatibility with N-octave band analysis.
See class docstring for parameter descriptions.
Source code in wandas/frames/noct.py
106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 | |
plot(plot_type='noct', ax=None, title=None, overlay=False, xlabel=None, ylabel=None, alpha=1.0, xlim=None, ylim=None, Aw=False, **kwargs)
¶
Plot the N-octave band data using various visualization strategies.
Supports standard plotting configurations for acoustic analysis, including decibel scales and A-weighting.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
plot_type
|
str
|
str, default="noct". Type of plot to create. The default "noct" type creates a step plot suitable for displaying N-octave band data. |
'noct'
|
ax
|
Axes | None
|
matplotlib.axes.Axes, optional. Axes to plot on. If None, creates new axes. |
None
|
title
|
str | None
|
str, optional. Title for the plot. If None, uses a default title with band specification. |
None
|
overlay
|
bool
|
bool, default=False. Whether to overlay all channels on a single plot (True) or create separate subplots for each channel (False). |
False
|
xlabel
|
str | None
|
str, optional. Label for the x-axis. If None, uses default "Center frequency [Hz]". |
None
|
ylabel
|
str | None
|
str, optional. Label for the y-axis. If None, uses default based on data type. |
None
|
alpha
|
float
|
float, default=1.0. Transparency level for the plot lines (0.0 to 1.0). |
1.0
|
xlim
|
tuple[float, float] | None
|
tuple[float, float], optional. Limits for the x-axis as (min, max) tuple. |
None
|
ylim
|
tuple[float, float] | None
|
tuple[float, float], optional. Limits for the y-axis as (min, max) tuple. |
None
|
Aw
|
bool
|
bool, default=False. Whether to apply A-weighting to the data. |
False
|
**kwargs
|
Any
|
dict. Additional matplotlib Line2D parameters (e.g., color, linewidth, linestyle). |
{}
|
Returns:
| Type | Description |
|---|---|
Axes | Iterator[Axes]
|
Union[Axes, Iterator[Axes]]: The matplotlib axes containing the plot, or an iterator of axes for multi-plot outputs. |
Examples:
>>> noct = spectrum.noct(n=3)
>>> # Basic 1/3-octave plot
>>> noct.plot()
>>> # Overlay with A-weighting
>>> noct.plot(overlay=True, Aw=True)
>>> # Custom styling
>>> noct.plot(title="1/3-Octave Spectrum", color="blue", linewidth=2)
Source code in wandas/frames/noct.py
244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 | |
wandas.frames.roughness.RoughnessFrame
¶
Bases: BaseFrame[NDArrayReal]
Frame for detailed roughness analysis with Bark-band information.
This frame contains specific roughness (R_spec) data organized by Bark frequency bands over time, calculated using the Daniel & Weber (1997) method.
The relationship between total roughness and specific roughness follows: R = 0.25 * sum(R_spec, axis=bark_bands)
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
data
|
Array
|
da.Array. Specific roughness data with shape: - (n_bark_bands, n_time) for mono signals, stored as (1, n_bark_bands, n_time) - (n_channels, n_bark_bands, n_time) for multi-channel signals where n_bark_bands is always 47. |
required |
sampling_rate
|
float
|
float. Sampling rate of the roughness time series in Hz. For overlap=0.5, this is approximately 10 Hz (100ms hop). For overlap=0.0, this is approximately 5 Hz (200ms hop). |
required |
bark_axis
|
NDArrayReal
|
NDArrayReal. Bark frequency axis with 47 values from 0.5 to 23.5 Bark. |
required |
overlap
|
float
|
float. Overlap coefficient used in the calculation (0.0 to 1.0). |
required |
label
|
str | None
|
str, optional. Frame label. Defaults to "roughness_spec". |
None
|
metadata
|
dict[str, Any] | None
|
dict, optional. Additional metadata. |
None
|
lineage
|
Any | None
|
LineageNode, optional. Constructor override for the runtime lineage. When omitted, a source node is
created. |
None
|
channel_metadata
|
Sequence[ChannelMetadata | dict[str, Any]] | None
|
list[ChannelMetadata], optional. Metadata for each channel. |
None
|
previous
|
BaseFrame[Any] | None
|
BaseFrame, optional. Immediate receiver Frame for process-local data comparison. For multi-input operations, follows only the left/base receiver. Not persisted in WDF. |
None
|
Attributes:
| Name | Type | Description |
|---|---|---|
bark_axis |
NDArrayReal
|
NDArrayReal. Frequency axis in Bark scale. |
n_bark_bands |
int
|
int. Number of Bark bands (always 47). |
n_time_points |
int
|
int. Number of time points. |
time |
NDArrayReal
|
NDArrayReal. Time axis based on sampling rate. |
overlap |
float
|
float. Overlap coefficient used (0.0 to 1.0). |
Examples:
Create a roughness frame from a signal:
>>> import wandas as wd
>>> signal = wd.read("motor.wav")
>>> roughness_spec = signal.roughness_dw_spec(overlap=0.5)
>>>
>>> # Plot Bark-Time heatmap
>>> roughness_spec.plot()
>>>
>>> # Find dominant Bark band
>>> dominant_idx = roughness_spec.data.mean(axis=1).argmax()
>>> dominant_bark = roughness_spec.bark_axis[dominant_idx]
>>> print(f"Dominant frequency: {dominant_bark:.1f} Bark")
>>>
>>> # Extract specific Bark band
>>> bark_10_idx = np.argmin(np.abs(roughness_spec.bark_axis - 10.0))
>>> roughness_at_10bark = roughness_spec.data[bark_10_idx, :]
The Daniel & Weber (1997) roughness model calculates specific roughness for 47 critical bands (Bark scale) over time, then integrates them to produce the total roughness:
.. math:: R = 0.25 \sum_{i=1}^{47} R'_i
where R'_i is the specific roughness in the i-th Bark band.
References
.. [1] Daniel, P., & Weber, R. (1997). "Psychoacoustical roughness: Implementation of an optimized model". Acta Acustica united with Acustica, 83(1), 113-123.
Source code in wandas/frames/roughness.py
23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 | |
Attributes¶
bark_axis = normalized_bark_axis.copy()
instance-attribute
¶
overlap = overlap
instance-attribute
¶
n_bark_bands
property
¶
Number of Bark bands.
Returns:
| Name | Type | Description |
|---|---|---|
int |
int
|
Always 47 for the Daniel & Weber model. |
n_time_points
property
¶
Number of time points in the roughness time series.
Returns:
| Name | Type | Description |
|---|---|---|
int |
int
|
Number of time frames in the analysis. |
time
property
¶
Time axis based on sampling rate.
Returns:
| Name | Type | Description |
|---|---|---|
NDArrayReal |
NDArrayReal
|
Time values in seconds for each frame. |
source_time
property
¶
Return roughness analysis time points on the source timeline.
Functions¶
__init__(data, sampling_rate, bark_axis, overlap, label=None, metadata=None, channel_metadata=None, channel_ids=None, previous=None, source_time_offset=0.0, lineage=None, operation_history_prefix=())
¶
Initialize a roughness tensor and its exact analysis state.
See the class docstring for parameter descriptions. The constructor requires
one finite Bark coordinate for each of the 47 model bands and an overlap in
the closed interval [0.0, 1.0].
Source code in wandas/frames/roughness.py
98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 | |
to_dataframe()
¶
DataFrame conversion is not supported for RoughnessFrame.
RoughnessFrame contains 3D data (channels, bark_bands, time_frames) which cannot be directly converted to a 2D DataFrame.
Raises:
| Type | Description |
|---|---|
NotImplementedError
|
Always raised as DataFrame conversion is not supported. |
Source code in wandas/frames/roughness.py
217 218 219 220 221 222 223 224 225 226 | |
plot(plot_type='heatmap', ax=None, title=None, cmap='viridis', vmin=None, vmax=None, xlabel='Time [s]', ylabel='Frequency [Bark]', colorbar_label='Specific Roughness [Asper/Bark]', **kwargs)
¶
Plot Bark-Time-Roughness heatmap.
For multi-channel signals, the mean across channels is plotted.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
plot_type
|
Literal['heatmap']
|
{"heatmap"}, default="heatmap". Plot strategy. Only the Bark-time heatmap is supported. |
'heatmap'
|
ax
|
Axes | None
|
Axes, optional. Matplotlib axes to plot on. If None, a new figure is created. |
None
|
title
|
str | None
|
str, optional. Plot title. If None, a default title is used. |
None
|
cmap
|
str
|
str, default="viridis". Colormap name for the heatmap. |
'viridis'
|
vmin
|
float | None
|
float, optional. Lower color scale limit. If None, automatic scaling is used. |
None
|
vmax
|
float | None
|
float, optional. Upper color scale limit. If None, automatic scaling is used. |
None
|
xlabel
|
str
|
str, default="Time [s]". Label for the x-axis. |
'Time [s]'
|
ylabel
|
str
|
str, default="Frequency [Bark]". Label for the y-axis. |
'Frequency [Bark]'
|
colorbar_label
|
str
|
str, default="Specific Roughness [Asper/Bark]". Label for the colorbar. |
'Specific Roughness [Asper/Bark]'
|
**kwargs
|
Any
|
Any. Additional keyword arguments passed to pcolormesh. |
{}
|
Returns:
| Name | Type | Description |
|---|---|---|
Axes |
Axes
|
The matplotlib axes object containing the plot. |
Raises:
| Type | Description |
|---|---|
ValueError
|
If |
Notes
Plotting is an explicit compute boundary.
Examples:
>>> import wandas as wd
>>> signal = wd.read("motor.wav")
>>> roughness_spec = signal.roughness_dw_spec(overlap=0.5)
>>> roughness_spec.plot(cmap="hot", title="Motor Roughness Analysis")
Source code in wandas/frames/roughness.py
235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 | |
wandas.frames.segments.FrameSegments
¶
Bases: Generic[F]
Equal-shaped Frames cut from recordings, one per event.
Segments are created by :meth:ChannelFrame.trigger or
:meth:ChannelFrame.segment, or by :meth:apply on existing segments.
Every segment keeps a stable event ID, its position among the detected
events, so IDs stay the same after selecting or dropping segments.
Indexing with an integer returns that segment's Frame. Indexing with a
slice, a list of positions or a boolean mask returns a new collection of the
same type. The collection is immutable and keeps every Frame lazy; only
:meth:to_numpy and plotting compute data.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
items
|
Sequence[F]
|
Frames of the same type, shape and sampling rate. |
required |
event_ids
|
Sequence[int] | NDArray[integer[Any]]
|
Unique integer ID of each segment. |
required |
trigger_times
|
Sequence[float] | NDArray[floating[Any]]
|
Event time in seconds of each segment, from the start of the Frame it was cut from. |
required |
Raises:
| Type | Description |
|---|---|
TypeError
|
If the items are not Frames of one type. |
ValueError
|
If shapes, sampling rates or the number of IDs and times differ, or IDs repeat. |
Source code in wandas/frames/segments.py
81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 | |
Attributes¶
event_ids
property
¶
Stable event ID of each segment, in segment order (read-only).
trigger_times
property
¶
Event time in seconds of each segment, from the start of the Frame it was cut from (read-only).
For segments cut by :meth:ChannelFrame.trigger or
:meth:ChannelFrame.segment, hits[i].source_time_offset + hits.pre
gives the event time on the original recording per channel, to the
nearest sample.
Functions¶
__init__(items, event_ids, trigger_times)
¶
Source code in wandas/frames/segments.py
106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 | |
__len__()
¶
Source code in wandas/frames/segments.py
160 161 | |
__iter__()
¶
Source code in wandas/frames/segments.py
163 164 | |
__getitem__(key)
¶
__getitem__(key: int) -> F
__getitem__(key: slice | Sequence[int] | npt.NDArray[np.integer[Any]] | npt.NDArray[np.bool_]) -> SegmentsT
Return one segment's Frame, or a subset of the segments.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
key
|
Any
|
Integer position for one Frame; a slice, a sequence of integer
positions or a boolean mask of length |
required |
Raises:
| Type | Description |
|---|---|
IndexError
|
If a position is out of range. |
TypeError
|
If the key is of another type. |
ValueError
|
If a mask has the wrong length or positions repeat. |
Source code in wandas/frames/segments.py
174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 | |
drop(event_ids)
¶
Return the segments without the given event IDs.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
event_ids
|
int | Sequence[int]
|
Event ID or IDs to exclude, as listed in
:attr: |
required |
Raises:
| Type | Description |
|---|---|
KeyError
|
If an ID is not among the segments. |
ValueError
|
If an ID repeats. |
Source code in wandas/frames/segments.py
213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 | |
apply(func)
¶
Apply a function to every segment and keep the results as segments.
Use this to process each segment separately, for example
hits.apply(lambda hit: hit.high_pass_filter(5)). Filtering each
segment differs from filtering the recording before cutting it, because
each segment starts its own filter state; choose the order that matches
the analysis. Event IDs and trigger times are kept, and each result
records its own lineage. For ChannelFrame results, the trigger channel
follows its channel ID through reordering, and pre follows a moved
segment start (read from source_time_offset).
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
func
|
Callable[[F], BaseFrame[Any]]
|
Function that takes one segment's Frame and returns a Frame. Every result must have the same type and shape. |
required |
Returns:
| Type | Description |
|---|---|
FrameSegments[Any]
|
class: |
Source code in wandas/frames/segments.py
237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 | |
to_numpy()
¶
Compute every segment and stack them as (segment, channel, ...).
The channel axis is kept even for single-channel segments, so a
ChannelFrame collection gives (segment, channel, sample). Values are
calibrated, as in each Frame's data. This computes the data.
Raises:
| Type | Description |
|---|---|
ValueError
|
If there are no segments. |
Source code in wandas/frames/segments.py
275 276 277 278 279 280 281 282 283 284 285 286 287 288 | |
__repr__()
¶
Source code in wandas/frames/segments.py
300 301 302 303 | |
wandas.frames.segments.ChannelFrameSegments
¶
Bases: FrameSegments[ChannelFrame]
Time segments of all channels, one per event, such as the hits of an impact test.
Every segment is a ChannelFrame of the same channels and length. Index it
to look at one hit, select or :meth:drop hits to exclude bad ones, and
average the remaining hits with :meth:transfer_function,
:meth:coherence and :meth:average.
Detection by :meth:ChannelFrame.trigger computes the trigger channel;
segments, selections and averages stay lazy until data is accessed.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
items
|
Sequence[ChannelFrame]
|
ChannelFrames with the same channels, units, calibration, length and sampling rate. |
required |
event_ids
|
Sequence[int] | NDArray[integer[Any]]
|
Unique integer ID of each segment. |
required |
trigger_times
|
Sequence[float] | NDArray[floating[Any]]
|
Event time in seconds of each segment, from the start of the Frame it was cut from. |
required |
pre
|
float
|
Time in seconds from the start of each segment to its event. |
0.0
|
trigger_channel
|
int | None
|
Index of the channel the events were detected on, or
|
None
|
Raises:
| Type | Description |
|---|---|
ValueError
|
If the segments differ in channels, units or calibration,
or |
Source code in wandas/frames/segments.py
306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 | |
Attributes¶
pre
property
¶
Time in seconds from the start of each segment to its event.
trigger_channel
property
¶
Label of the channel the events were detected on, if known.
sampling_rate
property
¶
Sampling rate of the segments in Hz.
labels
property
¶
Channel labels of the segments.
n_samples
property
¶
Samples per segment.
duration
property
¶
Segment length in seconds.
Functions¶
__init__(items, event_ids, trigger_times, pre=0.0, trigger_channel=None)
¶
Source code in wandas/frames/segments.py
332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 | |
transfer_function(input=None, output=None, n_fft=None, window='boxcar', detrend=None)
¶
Average the segments into one frequency response function (H1).
The estimate is H1 = sum(G_xy) / sum(G_xx) over the segments,
where x is the input channel and y an output channel: every
segment is one spectral average with equal weight, so this is not the
mean of per-segment transfer functions. Each segment is transformed
whole, without overlap. The defaults follow impact testing (rectangular
window, no detrending), unlike :meth:ChannelFrame.transfer_function,
which averages overlapping Hann windows within one recording. Bins where
the input has no energy are NaN.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
input
|
int | str | None
|
Excitation channel (index or label). |
None
|
output
|
int | str | Sequence[int | str] | None
|
Response channel or channels. |
None
|
n_fft
|
int | None
|
FFT length. |
None
|
window
|
str
|
SciPy window name applied to each segment. Default
|
'boxcar'
|
detrend
|
Literal['constant', 'linear'] | None
|
|
None
|
Returns:
| Type | Description |
|---|---|
TransferFunctionFrame
|
TransferFunctionFrame with one output/input pair per output channel. |
Raises:
| Type | Description |
|---|---|
ValueError
|
If there are no segments, the channels are invalid or
|
Source code in wandas/frames/segments.py
391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 | |
coherence(input=None, output=None, n_fft=None, window='boxcar', detrend=None)
¶
Coherence between the input and each output over the segments.
Computed from the same averaged cross and auto spectra as
:meth:transfer_function, with the same arguments. Values near 1 mean
the output follows the input linearly in every segment. At least two
segments are needed, because with one segment the coherence is always
1.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
input
|
int | str | None
|
Excitation channel (index or label). |
None
|
output
|
int | str | Sequence[int | str] | None
|
Response channel or channels. |
None
|
n_fft
|
int | None
|
FFT length. |
None
|
window
|
str
|
SciPy window name applied to each segment. Default
|
'boxcar'
|
detrend
|
Literal['constant', 'linear'] | None
|
|
None
|
Returns:
| Type | Description |
|---|---|
CoherenceFrame
|
CoherenceFrame with one output/input pair per output channel. |
Raises:
| Type | Description |
|---|---|
ValueError
|
If there are fewer than two segments, or as in
:meth: |
Source code in wandas/frames/segments.py
439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 | |
average()
¶
Average the segments sample by sample, aligned on their events.
This is the time-domain mean of the waveforms (synchronous averaging),
not an average spectrum. The result keeps the channels of the segments
and is lazy. Its source_time_offset is that of the first segment,
since an average has no single position on the recording.
Raises:
| Type | Description |
|---|---|
ValueError
|
If there are no segments. |
Source code in wandas/frames/segments.py
489 490 491 492 493 494 495 496 497 498 499 500 501 502 | |
plot(channel=None, ax=None)
¶
Overlay every segment, one axes per channel, labelled by event ID.
The time axis is in seconds from the event, so the trigger is at 0 and the pre-trigger part is negative. This computes the data.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
channel
|
int | str | None
|
Channel to plot (index or label). |
None
|
ax
|
Axes | None
|
Axes to draw on; only allowed with one channel. |
None
|
Returns:
| Type | Description |
|---|---|
Axes | Iterator[Axes]
|
The Axes for one channel, otherwise an iterator over the Axes. |
Raises:
| Type | Description |
|---|---|
ValueError
|
If there are no segments, or |
Source code in wandas/frames/segments.py
504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 | |
__repr__()
¶
Source code in wandas/frames/segments.py
642 643 644 645 646 647 648 649 | |
wandas.frames.mixins.channel_processing_mixin.ChannelProcessingMixin
¶
Mixin that provides methods related to signal processing.
This mixin provides processing methods applied to audio signals and other time-series data, such as signal processing filters and transformation operations.
Source code in wandas/frames/mixins/channel_processing_mixin.py
61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 839 840 841 842 843 844 845 846 847 848 849 850 851 852 853 854 855 856 857 858 859 860 861 862 863 864 865 866 867 868 869 870 871 872 873 874 875 876 877 878 879 880 881 882 883 884 885 886 887 888 889 890 891 892 893 894 895 896 897 898 899 900 901 902 903 904 905 906 907 908 909 910 911 912 913 914 915 916 917 918 919 920 921 922 923 924 925 926 927 928 929 930 931 932 933 934 935 936 937 938 939 940 941 942 943 944 945 946 947 948 949 950 951 952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967 968 969 970 971 972 973 974 975 976 977 978 979 980 981 982 983 984 985 986 987 988 989 990 991 992 993 994 995 996 997 998 999 1000 1001 1002 1003 1004 1005 1006 1007 1008 1009 1010 1011 1012 1013 1014 1015 1016 1017 1018 1019 1020 1021 1022 1023 1024 1025 1026 1027 1028 1029 1030 1031 1032 1033 1034 1035 1036 1037 1038 1039 1040 1041 1042 1043 1044 1045 1046 1047 1048 1049 1050 1051 1052 1053 1054 1055 1056 1057 1058 1059 1060 1061 1062 1063 1064 1065 1066 1067 1068 1069 1070 1071 1072 1073 1074 1075 1076 1077 1078 1079 1080 1081 1082 1083 1084 1085 1086 1087 1088 1089 1090 1091 1092 1093 1094 1095 1096 1097 1098 1099 1100 1101 1102 1103 1104 1105 1106 1107 1108 1109 1110 1111 1112 1113 1114 1115 1116 1117 1118 1119 1120 1121 1122 1123 1124 1125 1126 1127 1128 1129 1130 1131 1132 1133 1134 1135 1136 1137 1138 1139 1140 1141 1142 1143 1144 1145 1146 1147 1148 1149 1150 1151 1152 1153 1154 1155 1156 1157 1158 1159 1160 1161 1162 1163 1164 1165 1166 1167 1168 1169 1170 1171 1172 1173 1174 1175 1176 1177 1178 1179 1180 1181 1182 1183 1184 1185 1186 1187 1188 1189 1190 1191 1192 1193 1194 1195 1196 1197 1198 1199 1200 1201 1202 1203 1204 1205 1206 1207 1208 1209 1210 1211 1212 1213 1214 1215 1216 1217 1218 1219 1220 1221 1222 1223 1224 1225 1226 1227 1228 1229 1230 1231 1232 1233 1234 1235 1236 1237 1238 1239 1240 1241 1242 1243 1244 1245 1246 1247 1248 1249 1250 1251 1252 1253 1254 1255 1256 1257 1258 1259 1260 1261 1262 1263 1264 1265 1266 1267 1268 1269 1270 1271 1272 1273 1274 1275 1276 1277 1278 1279 1280 1281 1282 1283 1284 1285 1286 1287 1288 1289 1290 1291 1292 1293 1294 1295 1296 1297 1298 1299 1300 1301 1302 1303 1304 1305 1306 1307 1308 1309 1310 1311 1312 1313 1314 1315 1316 1317 1318 1319 1320 1321 1322 1323 1324 1325 1326 1327 1328 1329 1330 1331 1332 1333 1334 1335 1336 1337 1338 1339 1340 1341 1342 1343 1344 1345 1346 1347 1348 1349 1350 1351 1352 1353 1354 1355 1356 1357 1358 1359 1360 1361 1362 1363 1364 1365 1366 1367 1368 1369 1370 | |
Functions¶
apply(func, output_shape_func=None, output_frame_class=None, output_frame_kwargs=None, *, dask_pure=True, **kwargs)
¶
apply(func: Callable[..., Any], output_shape_func: Callable[[tuple[int, ...]], tuple[int, ...]] | None = ..., output_frame_class: None = ..., output_frame_kwargs: dict[str, Any] | None = ..., *, dask_pure: bool = ..., **kwargs: Any) -> T_Processing
apply(func: Callable[..., Any], output_shape_func: Callable[[tuple[int, ...]], tuple[int, ...]] | None = ..., output_frame_class: type[T_OutputFrame] = ..., output_frame_kwargs: dict[str, Any] | None = ..., *, dask_pure: bool = ..., **kwargs: Any) -> T_OutputFrame
Apply a custom function to the signal.
Builds a lazy Dask graph without calling func. At materialization,
func receives the whole channel-first NumPy array, shaped
(channels, samples) for a ChannelFrame, including mono input.
Return a new array with the input dtype; this API has no output-dtype
callback. Supply output_shape_func for sample-axis shape changes
that preserve channel count and order. Channel-count changes also need
new channel IDs and metadata; use a standard channel operation or an
extension that constructs the output Frame explicitly.
The input Frame remains unchanged. Array slicing alone does not update
sampling rate or source-time offsets; use a semantic operation such as
trim() when appropriate. Runtime lineage is recorded, but arbitrary
callables cannot be extracted into a Recipe. See the
custom processing how-to for the
external Operation and Recipe extension path.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
func
|
Callable[..., Any]
|
Function to apply. |
required |
output_shape_func
|
Callable[[tuple[int, ...]], tuple[int, ...]] | None
|
Optional function to calculate output shape. |
None
|
output_frame_class
|
type[T_OutputFrame] | None
|
Optional frame class for the output. When
provided, the result is wrapped in this class instead of the
caller's type, enabling domain transitions (e.g.
|
None
|
output_frame_kwargs
|
dict[str, Any] | None
|
Extra constructor keyword arguments required
by output_frame_class (e.g. |
None
|
dask_pure
|
bool
|
Dask execution-control flag for delayed custom
operations. Set to |
True
|
**kwargs
|
Any
|
Additional arguments for the function. Formal |
{}
|
Returns:
| Type | Description |
|---|---|
Any
|
New frame with the custom function applied. |
Raises:
| Type | Description |
|---|---|
ValueError
|
If function kwargs use reserved |
Examples:
>>> import numpy as np
>>> import wandas as wd
>>> frame = wd.from_numpy(np.ones((2, 8)), sampling_rate=8000)
>>> scaled = frame.apply(lambda data, factor: data * factor, factor=2.0)
>>> np.allclose(scaled.to_numpy(), 2.0)
True
Source code in wandas/frames/mixins/channel_processing_mixin.py
171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 | |
high_pass_filter(cutoff, order=4)
¶
Apply a high-pass filter to the signal.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
cutoff
|
float
|
Filter cutoff frequency (Hz) |
required |
order
|
int
|
Filter order. Default is 4. |
4
|
Returns:
| Type | Description |
|---|---|
T_Processing
|
New ChannelFrame after filter application |
Raises:
| Type | Description |
|---|---|
ValueError
|
If the cutoff is outside |
Source code in wandas/frames/mixins/channel_processing_mixin.py
268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 | |
low_pass_filter(cutoff, order=4)
¶
Apply a low-pass filter to the signal.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
cutoff
|
float
|
Filter cutoff frequency (Hz) |
required |
order
|
int
|
Filter order. Default is 4. |
4
|
Returns:
| Type | Description |
|---|---|
T_Processing
|
New ChannelFrame after filter application |
Raises:
| Type | Description |
|---|---|
ValueError
|
If the cutoff is outside |
Source code in wandas/frames/mixins/channel_processing_mixin.py
288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 | |
band_pass_filter(low_cutoff, high_cutoff, order=4)
¶
Apply a band-pass filter to the signal.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
low_cutoff
|
float
|
Lower cutoff frequency (Hz) |
required |
high_cutoff
|
float
|
Higher cutoff frequency (Hz) |
required |
order
|
int
|
Order of the Butterworth prototype. Default is 4. The
band-pass filter itself has order |
4
|
Returns:
| Type | Description |
|---|---|
T_Processing
|
New ChannelFrame after filter application |
Raises:
| Type | Description |
|---|---|
ValueError
|
If either cutoff is outside |
Source code in wandas/frames/mixins/channel_processing_mixin.py
308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 | |
normalize(norm=float('inf'), axis=-1, threshold=None, fill=None)
¶
Normalize signal levels using NumPy-based normalization.
This method normalizes the signal amplitude according to the specified norm.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
norm
|
float | None
|
Norm type. Default is np.inf (maximum absolute value normalization). Supported values: - np.inf: Maximum absolute value normalization - -np.inf: Minimum absolute value normalization - 0: Peak normalization - float: Lp norm - None: No normalization |
float('inf')
|
axis
|
int | None
|
Axis along which to normalize. Default is -1 (time axis). - -1: Normalize along time axis (each channel independently) - None: Global normalization across all axes - int: Normalize along specified axis |
-1
|
threshold
|
float | None
|
Threshold below which values are considered zero. If None, no threshold is applied. |
None
|
fill
|
bool | None
|
Value to fill when the norm is zero. If None, the zero vector remains zero. |
None
|
Returns:
| Type | Description |
|---|---|
T_Processing
|
New ChannelFrame containing the normalized signal |
Examples:
>>> import wandas as wd
>>> signal = wd.read("audio.wav")
>>> # Normalize to maximum absolute value of 1.0 (per channel)
>>> normalized = signal.normalize()
>>> # Global normalization across all channels
>>> normalized_global = signal.normalize(axis=None)
>>> # L2 normalization
>>> normalized_l2 = signal.normalize(norm=2)
Source code in wandas/frames/mixins/channel_processing_mixin.py
369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 | |
remove_dc()
¶
Remove DC component (DC offset) from the signal.
This method removes the DC (direct current) component by subtracting the mean value from each channel. This is equivalent to centering the signal around zero.
Returns:
| Type | Description |
|---|---|
T_Processing
|
New ChannelFrame with DC component removed |
Examples:
>>> import wandas as wd
>>> import numpy as np
>>> # Create signal with DC offset
>>> signal = wd.read("audio.wav")
>>> signal_with_dc = signal + 2.0 # Add DC offset
>>> # Remove DC offset
>>> signal_clean = signal_with_dc.remove_dc()
>>> # Verify DC removal
>>> assert np.allclose(signal_clean.data.mean(axis=1), 0, atol=1e-10)
Notes
- This operation is performed per channel
- Equivalent to applying a high-pass filter with very low cutoff
- Useful for removing sensor drift or measurement offset
Source code in wandas/frames/mixins/channel_processing_mixin.py
415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 | |
a_weighting()
¶
Apply A-weighting filter to the signal.
A-weighting adjusts the frequency response to approximate human auditory perception using the implemented digital curve. This returns a weighted linear waveform in the input unit; it does not calculate RMS or convert to dB. No sound-level-meter conformance is implied.
Returns:
| Type | Description |
|---|---|
T_Processing
|
New ChannelFrame containing the A-weighted signal |
Source code in wandas/frames/mixins/channel_processing_mixin.py
446 447 448 449 450 451 452 453 454 455 456 457 458 459 | |
abs()
¶
Compute the absolute value of the signal.
Returns:
| Type | Description |
|---|---|
T_Processing
|
New ChannelFrame containing the absolute values |
Source code in wandas/frames/mixins/channel_processing_mixin.py
461 462 463 464 465 466 467 468 469 | |
power(exponent=2.0)
¶
Compute the power of the signal.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
exponent
|
float
|
Exponent to raise the signal to. Default is 2.0. |
2.0
|
Returns:
| Type | Description |
|---|---|
T_Processing
|
New ChannelFrame containing the powered signal |
Source code in wandas/frames/mixins/channel_processing_mixin.py
471 472 473 474 475 476 477 478 479 480 481 482 | |
sum()
¶
Sum all channels.
The result keeps the unit and reference shared by all channels and the first channel's source-time offset.
Returns:
| Type | Description |
|---|---|
T_Processing
|
A new ChannelFrame with summed signal. |
Source code in wandas/frames/mixins/channel_processing_mixin.py
531 532 533 534 535 536 537 538 539 540 541 | |
mean()
¶
Average all channels.
The result keeps the unit and reference shared by all channels and the first channel's source-time offset.
Returns:
| Type | Description |
|---|---|
T_Processing
|
A new ChannelFrame with averaged signal. |
Source code in wandas/frames/mixins/channel_processing_mixin.py
543 544 545 546 547 548 549 550 551 552 553 | |
trim(start=0, end=None)
¶
Trim the signal to the specified time range.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
start
|
float
|
Start time (seconds) |
0
|
end
|
float | None
|
End time (seconds) |
None
|
Returns:
| Type | Description |
|---|---|
T_Processing
|
New ChannelFrame containing the trimmed signal. The operation is a lazy structural time slice: channel labels, calibration, metadata, and channel IDs are preserved, while source-time offsets advance to the first selected sample. |
Raises:
| Type | Description |
|---|---|
ValueError
|
If either time is negative or end is earlier than start |
Source code in wandas/frames/mixins/channel_processing_mixin.py
591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 | |
fix_length(length=None, duration=None)
¶
Adjust the signal to the specified length.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
duration
|
float | None
|
Signal length in seconds |
None
|
length
|
int | None
|
Signal length in samples |
None
|
Returns:
| Type | Description |
|---|---|
T_Processing
|
New ChannelFrame containing the adjusted signal |
Source code in wandas/frames/mixins/channel_processing_mixin.py
621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 | |
rms_trend(frame_length=2048, hop_length=512, dB=False, Aw=False)
¶
Compute a linear RMS trend or an RMS amplitude level.
This method calculates root mean square over centered, zero-padded
sliding windows. Calibration is applied per channel. With dB=False
the output is linear and retains each channel's physical unit. With
dB=True the output is 20 * log10(max(window_rms / channel_ref,
1e-12)), bounded below by -240 dB. It is dB SPL only when the signal
is pressure in Pa and the reference is 2e-5 Pa.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
frame_length
|
int
|
Size of the sliding window in samples. Default is 2048. |
2048
|
hop_length
|
int
|
Hop length between windows in samples. Default is 512. |
512
|
dB
|
bool
|
Return amplitude level relative to each channel reference. Default is False. |
False
|
Aw
|
bool
|
Apply the implemented A-frequency-weighting filter before RMS. Default is False. |
False
|
Returns:
| Type | Description |
|---|---|
T_Processing
|
New lazy ChannelFrame with shape |
Raises:
| Type | Description |
|---|---|
ValueError
|
If the window parameters or channel references are invalid for the input Frame. |
Source code in wandas/frames/mixins/channel_processing_mixin.py
640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 | |
sound_level(freq_weighting='Z', time_weighting='Fast', dB=False)
¶
Compute a frequency- and time-weighted RMS or reference-relative level.
The selected frequency weighting is applied first. Squared samples are
then smoothed by a first-order exponential filter using 125 ms (Fast)
or 1 s (Slow). With dB=False the square root is returned in the
calibrated input unit. With dB=True the result is
10 * log10(max(smoothed_power / channel_ref**2, 1e-20)), bounded
below by -200 dB. A Pa channel whose reference is 2e-5 Pa yields
dB SPL; an uncalibrated channel yields relative dB re 1 input unit.
This method validates the implemented filters and time constants, not the complete tolerance, detector, calibration, or directional-response requirements of an IEC/JIS sound-level meter.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
freq_weighting
|
Literal['A', 'C', 'Z'] | None
|
Implemented frequency-weighting curve: |
'Z'
|
time_weighting
|
Literal['Fast', 'Slow']
|
Exponential time constant: |
'Fast'
|
dB
|
bool
|
Return level relative to the channel reference when |
False
|
Returns:
| Type | Description |
|---|---|
T_Processing
|
New lazy ChannelFrame with shape |
Raises:
| Type | Description |
|---|---|
ValueError
|
If the frequency/time weighting or channel references are invalid for the input Frame. |
Source code in wandas/frames/mixins/channel_processing_mixin.py
695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 | |
channel_difference(other_channel=0)
¶
Compute index-wise differences between channels.
channel_difference subtracts the selected reference channel from
each channel at the current array indices. It does not compare
per-channel source_time_offset values and does not perform
source-time alignment.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
other_channel
|
int | str
|
Index or label of the reference channel. Default is 0. |
0
|
Returns:
| Type | Description |
|---|---|
T_Processing
|
New ChannelFrame containing the channel difference with the input source-time offsets preserved. |
Source code in wandas/frames/mixins/channel_processing_mixin.py
752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 | |
resampling(target_sr, **kwargs)
¶
Resample audio data.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
target_sr
|
float
|
Target sampling rate (Hz) |
required |
**kwargs
|
Any
|
Additional resampling parameters |
{}
|
Returns:
| Type | Description |
|---|---|
T_Processing
|
Resampled ChannelFrame |
Source code in wandas/frames/mixins/channel_processing_mixin.py
775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 797 | |
hpss_harmonic(kernel_size=31, power=2, margin=1, n_fft=2048, hop_length=None, win_length=None, window='hann', center=True, pad_mode='constant')
¶
Extract harmonic components using HPSS (Harmonic-Percussive Source Separation).
This method separates the harmonic (tonal) components from the signal.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
kernel_size
|
HpssKernelSize
|
Median filter size for HPSS. |
31
|
power
|
HpssFloatLike
|
Exponent for the Weiner filter used in HPSS. |
2
|
margin
|
HpssMargin
|
Margin size for the separation. |
1
|
n_fft
|
HpssIntLike
|
Size of FFT window. |
2048
|
hop_length
|
HpssIntLike | None
|
Hop length for STFT. |
None
|
win_length
|
HpssIntLike | None
|
Window length for STFT. |
None
|
window
|
Any
|
Window type for STFT. |
'hann'
|
center
|
bool
|
If True, center the frames. |
True
|
pad_mode
|
str
|
Padding mode for STFT. |
'constant'
|
Returns:
| Type | Description |
|---|---|
T_Processing
|
A new ChannelFrame containing the harmonic components. |
Source code in wandas/frames/mixins/channel_processing_mixin.py
827 828 829 830 831 832 833 834 835 836 837 838 839 840 841 842 843 844 845 846 847 848 849 850 851 852 853 854 855 856 857 858 859 860 861 862 863 864 865 866 867 868 869 870 871 872 873 874 | |
hpss_percussive(kernel_size=31, power=2, margin=1, n_fft=2048, hop_length=None, win_length=None, window='hann', center=True, pad_mode='constant')
¶
Extract percussive components using HPSS (Harmonic-Percussive Source Separation).
This method separates the percussive (tonal) components from the signal.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
kernel_size
|
HpssKernelSize
|
Median filter size for HPSS. |
31
|
power
|
HpssFloatLike
|
Exponent for the Weiner filter used in HPSS. |
2
|
margin
|
HpssMargin
|
Margin size for the separation. |
1
|
Returns:
| Type | Description |
|---|---|
T_Processing
|
A new ChannelFrame containing the harmonic components. |
Source code in wandas/frames/mixins/channel_processing_mixin.py
876 877 878 879 880 881 882 883 884 885 886 887 888 889 890 891 892 893 894 895 896 897 898 899 900 901 902 903 904 905 906 907 908 909 910 911 912 913 914 915 916 917 | |
loudness_zwtv(field_type='free')
¶
Calculate time-varying loudness using Zwicker method (ISO 532-1:2017).
This method computes the loudness of non-stationary signals according to the Zwicker method, as specified in ISO 532-1:2017. The loudness is calculated in sones, where a doubling of sones corresponds to a doubling of perceived loudness.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
field_type
|
Literal['free', 'diffuse']
|
Type of sound field. Options: - 'free': Free field (sound from a specific direction) - 'diffuse': Diffuse field (sound from all directions) Default is 'free'. |
'free'
|
Returns:
| Type | Description |
|---|---|
T_Processing
|
New ChannelFrame containing time-varying loudness values in sones. Each channel is processed independently. The output sampling rate is adjusted based on the loudness calculation time resolution (typically ~500 Hz for 2ms steps). |
Raises:
| Type | Description |
|---|---|
ValueError
|
If field_type is not 'free' or 'diffuse' |
Examples:
Calculate loudness for a signal:
>>> import wandas as wd
>>> signal = wd.read("audio.wav")
>>> loudness = signal.loudness_zwtv(field_type="free")
>>> loudness.plot(title="Time-varying Loudness")
Compare free field and diffuse field:
>>> loudness_free = signal.loudness_zwtv(field_type="free")
>>> loudness_diffuse = signal.loudness_zwtv(field_type="diffuse")
Notes
- The output contains time-varying loudness values in sones
- Typical loudness: 1 sone ≈ 40 phon (loudness level)
- The time resolution is approximately 2ms (determined by the algorithm)
- For multi-channel signals, loudness is calculated per channel
- The output sampling rate is updated to reflect the time resolution
Time axis convention: The time axis in the returned frame represents the start time of each 2ms analysis step. This differs slightly from the MoSQITo library, which uses the center time of each step. For example:
- wandas time: [0.000s, 0.002s, 0.004s, ...] (step start)
- MoSQITo time: [0.001s, 0.003s, 0.005s, ...] (step center)
The difference is very small (~1ms) and does not affect the loudness values themselves. This design choice ensures consistency with wandas's time axis convention across all frame types.
References
ISO 532-1:2017, "Acoustics — Methods for calculating loudness — Part 1: Zwicker method"
Source code in wandas/frames/mixins/channel_processing_mixin.py
919 920 921 922 923 924 925 926 927 928 929 930 931 932 933 934 935 936 937 938 939 940 941 942 943 944 945 946 947 948 949 950 951 952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967 968 969 970 971 972 973 974 975 976 977 978 979 980 981 | |
loudness_zwst(field_type='free')
¶
Calculate steady-state loudness using Zwicker method (ISO 532-1:2017).
This method computes the loudness of stationary (steady) signals according to the Zwicker method, as specified in ISO 532-1:2017. The loudness is calculated in sones, where a doubling of sones corresponds to a doubling of perceived loudness.
This method is suitable for analyzing steady sounds such as fan noise, constant machinery sounds, or other stationary signals.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
field_type
|
Literal['free', 'diffuse']
|
Type of sound field. Options: - 'free': Free field (sound from a specific direction) - 'diffuse': Diffuse field (sound from all directions) Default is 'free'. |
'free'
|
Returns:
| Type | Description |
|---|---|
NDArrayReal
|
Loudness values in sones, one per channel. Shape: (n_channels,) |
Raises:
| Type | Description |
|---|---|
ValueError
|
If field_type is not 'free' or 'diffuse' |
Examples:
Calculate steady-state loudness for a fan noise:
>>> import wandas as wd
>>> signal = wd.read("fan_noise.wav")
>>> loudness = signal.loudness_zwst(field_type="free")
>>> print(f"Channel 0 loudness: {loudness[0]:.2f} sones")
>>> print(f"Mean loudness: {loudness.mean():.2f} sones")
Compare free field and diffuse field:
>>> loudness_free = signal.loudness_zwst(field_type="free")
>>> loudness_diffuse = signal.loudness_zwst(field_type="diffuse")
>>> print(f"Free field: {loudness_free[0]:.2f} sones")
>>> print(f"Diffuse field: {loudness_diffuse[0]:.2f} sones")
Notes
- Returns a 1D array with one loudness value per channel
- Typical loudness: 1 sone ≈ 40 phon (loudness level)
- For multi-channel signals, loudness is calculated independently per channel
- This method is designed for stationary signals (constant sounds)
- For time-varying signals, use loudness_zwtv() instead
- Similar to the rms property, returns NDArrayReal for consistency
References
ISO 532-1:2017, "Acoustics — Methods for calculating loudness — Part 1: Zwicker method"
Source code in wandas/frames/mixins/channel_processing_mixin.py
983 984 985 986 987 988 989 990 991 992 993 994 995 996 997 998 999 1000 1001 1002 1003 1004 1005 1006 1007 1008 1009 1010 1011 1012 1013 1014 1015 1016 1017 1018 1019 1020 1021 1022 1023 1024 1025 1026 1027 1028 1029 1030 1031 1032 1033 1034 1035 1036 1037 | |
roughness_dw(overlap=0.5)
¶
Calculate time-varying roughness using Daniel and Weber method.
Roughness is a psychoacoustic metric that quantifies the perceived harshness or roughness of a sound, measured in asper. This method implements the Daniel & Weber (1997) standard calculation.
The calculation follows the standard formula: R = 0.25 * sum(R'_i) for i=1 to 47 Bark bands
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
overlap
|
float
|
Overlapping coefficient for 200ms analysis windows (0.0 to 1.0). - overlap=0.5: 100ms hop → ~10 Hz output sampling rate - overlap=0.0: 200ms hop → ~5 Hz output sampling rate Default is 0.5. |
0.5
|
Returns:
| Type | Description |
|---|---|
T_Processing
|
New ChannelFrame containing time-varying roughness values in asper. The output sampling rate depends on the overlap parameter. |
Raises:
| Type | Description |
|---|---|
ValueError
|
If overlap is not in the range [0.0, 1.0] |
Examples:
Calculate roughness for a motor noise:
>>> import wandas as wd
>>> signal = wd.read("motor_noise.wav")
>>> roughness = signal.roughness_dw(overlap=0.5)
>>> roughness.plot(ylabel="Roughness [asper]")
Analyze roughness statistics:
>>> mean_roughness = roughness.data.mean()
>>> max_roughness = roughness.data.max()
>>> print(f"Mean: {mean_roughness:.2f} asper")
>>> print(f"Max: {max_roughness:.2f} asper")
Compare before and after modification:
>>> before = wd.read("motor_before.wav").roughness_dw()
>>> after = wd.read("motor_after.wav").roughness_dw()
>>> improvement = before.data.mean() - after.data.mean()
>>> print(f"Roughness reduction: {improvement:.2f} asper")
Notes
- Returns a ChannelFrame with time-varying roughness values
- Typical roughness values: 0-2 asper for most sounds
- Higher values indicate rougher, harsher sounds
- For multi-channel signals, roughness is calculated independently per channel
- This is the standard-compliant total roughness (R)
- For detailed Bark-band analysis, use roughness_dw_spec() instead
Time axis convention: The time axis in the returned frame represents the start time of each 200ms analysis window. This differs from the MoSQITo library, which uses the center time of each window. For example:
- wandas time: [0.0s, 0.1s, 0.2s, ...] (window start)
- MoSQITo time: [0.1s, 0.2s, 0.3s, ...] (window center)
The difference is constant (half the window duration = 100ms) and does not affect the roughness values themselves. This design choice ensures consistency with wandas's time axis convention across all frame types.
References
Daniel, P., & Weber, R. (1997). "Psychoacoustical roughness: Implementation of an optimized model." Acustica, 83, 113-123.
Source code in wandas/frames/mixins/channel_processing_mixin.py
1039 1040 1041 1042 1043 1044 1045 1046 1047 1048 1049 1050 1051 1052 1053 1054 1055 1056 1057 1058 1059 1060 1061 1062 1063 1064 1065 1066 1067 1068 1069 1070 1071 1072 1073 1074 1075 1076 1077 1078 1079 1080 1081 1082 1083 1084 1085 1086 1087 1088 1089 1090 1091 1092 1093 1094 1095 1096 1097 1098 1099 1100 1101 1102 1103 1104 1105 1106 1107 1108 1109 1110 | |
roughness_dw_spec(overlap=0.5)
¶
Calculate specific roughness with Bark-band frequency information.
This method returns detailed roughness analysis data organized by Bark frequency bands over time, allowing for frequency-specific roughness analysis. It uses the Daniel & Weber (1997) method.
The relationship between total roughness and specific roughness: R = 0.25 * sum(R'_i) for i=1 to 47 Bark bands
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
overlap
|
float
|
Overlapping coefficient for 200ms analysis windows (0.0 to 1.0). - overlap=0.5: 100ms hop → ~10 Hz output sampling rate - overlap=0.0: 200ms hop → ~5 Hz output sampling rate Default is 0.5. |
0.5
|
Returns:
| Type | Description |
|---|---|
RoughnessFrame
|
RoughnessFrame containing: - data: Specific roughness by Bark band, shape (47, n_time) for mono or (n_channels, 47, n_time) for multi-channel - bark_axis: Frequency axis in Bark scale (47 values, 0.5-23.5) - time: Time axis for each analysis frame - overlap: Overlap coefficient used - plot(): Method for Bark-Time heatmap visualization |
Raises:
| Type | Description |
|---|---|
ValueError
|
If overlap is not in the range [0.0, 1.0] |
Examples:
Analyze frequency-specific roughness:
>>> import wandas as wd
>>> import numpy as np
>>> signal = wd.read("motor.wav")
>>> roughness_spec = signal.roughness_dw_spec(overlap=0.5)
>>>
>>> # Plot Bark-Time heatmap
>>> roughness_spec.plot(cmap="viridis", title="Roughness Analysis")
>>>
>>> # Find dominant Bark band
>>> dominant_idx = roughness_spec.data.mean(axis=1).argmax()
>>> dominant_bark = roughness_spec.bark_axis[dominant_idx]
>>> print(f"Most contributing band: {dominant_bark:.1f} Bark")
>>>
>>> # Extract specific Bark band time series
>>> bark_10_idx = np.argmin(np.abs(roughness_spec.bark_axis - 10.0))
>>> roughness_at_10bark = roughness_spec.data[bark_10_idx, :]
>>>
>>> # Verify standard formula
>>> total_roughness = 0.25 * roughness_spec.data.sum(axis=-2)
>>> # This should match signal.roughness_dw(overlap=0.5).data
Notes
- Returns a RoughnessFrame (not ChannelFrame)
- Contains 47 Bark bands from 0.5 to 23.5 Bark
- Each Bark band corresponds to a critical band of hearing
- Useful for identifying which frequencies contribute most to roughness
- The specific roughness can be integrated to obtain total roughness
- For simple time-series analysis, use roughness_dw() instead
Time axis convention: The time axis represents the start time of each 200ms analysis window, consistent with roughness_dw() and other wandas methods.
References
Daniel, P., & Weber, R. (1997). "Psychoacoustical roughness: Implementation of an optimized model." Acustica, 83, 113-123.
Source code in wandas/frames/mixins/channel_processing_mixin.py
1112 1113 1114 1115 1116 1117 1118 1119 1120 1121 1122 1123 1124 1125 1126 1127 1128 1129 1130 1131 1132 1133 1134 1135 1136 1137 1138 1139 1140 1141 1142 1143 1144 1145 1146 1147 1148 1149 1150 1151 1152 1153 1154 1155 1156 1157 1158 1159 1160 1161 1162 1163 1164 1165 1166 1167 1168 1169 1170 1171 1172 1173 1174 1175 1176 1177 1178 1179 1180 1181 1182 1183 1184 1185 1186 1187 1188 1189 1190 1191 1192 1193 1194 1195 1196 1197 1198 1199 1200 1201 1202 1203 1204 1205 1206 1207 1208 1209 1210 1211 1212 1213 1214 1215 1216 1217 1218 1219 1220 1221 1222 1223 | |
fade(fade_ms=50)
¶
Apply symmetric fade-in and fade-out to the signal using Tukey window.
This method applies a symmetric fade-in and fade-out envelope to the signal using a Tukey (tapered cosine) window. The fade duration is the same for both the beginning and end of the signal.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
fade_ms
|
float
|
Fade duration in milliseconds for each end of the signal. The total fade duration is 2 * fade_ms. Default is 50 ms. Must be positive and less than half the signal duration. |
50
|
Returns:
| Type | Description |
|---|---|
T_Processing
|
New ChannelFrame containing the faded signal |
Raises:
| Type | Description |
|---|---|
ValueError
|
If fade_ms is negative or too long for the signal |
Examples:
>>> import wandas as wd
>>> signal = wd.read("audio.wav")
>>> # Apply 10ms fade-in and fade-out
>>> faded = signal.fade(fade_ms=10.0)
>>> # Apply very short fade (almost no effect)
>>> faded_short = signal.fade(fade_ms=0.1)
Notes
- Uses SciPy's Tukey window for smooth fade transitions
- Fade is applied symmetrically to both ends of the signal
- The Tukey window alpha parameter is computed automatically based on the fade duration and signal length
- For multi-channel signals, the same fade envelope is applied to all channels
- Lazy evaluation is preserved - computation occurs only when needed
Source code in wandas/frames/mixins/channel_processing_mixin.py
1225 1226 1227 1228 1229 1230 1231 1232 1233 1234 1235 1236 1237 1238 1239 1240 1241 1242 1243 1244 1245 1246 1247 1248 1249 1250 1251 1252 1253 1254 1255 1256 1257 1258 1259 1260 1261 1262 1263 | |
sharpness_din(weighting='din', field_type='free')
¶
Calculate sharpness using DIN 45692 method.
This method computes the time-varying sharpness of the signal according to DIN 45692 standard, which quantifies the perceived sharpness of sounds.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
weighting
|
Literal['din', 'aures', 'bismarck', 'fastl']
|
Default "din". Weighting type for sharpness calculation. Options: - 'din': DIN 45692 method - 'aures': Aures method - 'bismarck': Bismarck method - 'fastl': Fastl method |
'din'
|
field_type
|
Literal['free', 'diffuse']
|
Default "free". Type of sound field. Options: - 'free': Free field (sound from a specific direction) - 'diffuse': Diffuse field (sound from all directions) |
'free'
|
Returns:
| Name | Type | Description |
|---|---|---|
T_Processing |
T_Processing
|
New ChannelFrame containing sharpness time series in acum. The output sampling rate is approximately 500 Hz (2ms time steps). |
Raises:
| Type | Description |
|---|---|
ValueError
|
If the signal sampling rate is not supported by the algorithm. |
Examples:
>>> import wandas as wd
>>> signal = wd.read("sharp_sound.wav")
>>> sharpness = signal.sharpness_din(weighting="din", field_type="free")
>>> print(f"Mean sharpness: {sharpness.data.mean():.2f} acum")
Notes
- Sharpness is measured in acum (acum = 1 when the sound has the
same sharpness as a 2 kHz narrow-band noise at 60 dB SPL) - The calculation uses MoSQITo's implementation of DIN 45692 - Output sampling rate is fixed at 500 Hz regardless of input rate - For multi-channel signals, sharpness is calculated per channel
References
.. [1] DIN 45692:2009, "Measurement technique for the simulation of the auditory sensation of sharpness"
Source code in wandas/frames/mixins/channel_processing_mixin.py
1265 1266 1267 1268 1269 1270 1271 1272 1273 1274 1275 1276 1277 1278 1279 1280 1281 1282 1283 1284 1285 1286 1287 1288 1289 1290 1291 1292 1293 1294 1295 1296 1297 1298 1299 1300 1301 1302 1303 1304 1305 1306 1307 1308 1309 1310 1311 1312 1313 1314 1315 1316 1317 1318 1319 1320 1321 | |
sharpness_din_st(weighting='din', field_type='free')
¶
Calculate steady-state sharpness using DIN 45692 method.
This method computes the steady-state sharpness of the signal according to DIN 45692 standard, which quantifies the perceived sharpness of stationary sounds.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
weighting
|
Literal['din', 'aures', 'bismarck', 'fastl']
|
Default "din". Weighting type for sharpness calculation. Options: - 'din': DIN 45692 method - 'aures': Aures method - 'bismarck': Bismarck method - 'fastl': Fastl method |
'din'
|
field_type
|
Literal['free', 'diffuse']
|
Default "free". Type of sound field. Options: - 'free': Free field (sound from a specific direction) - 'diffuse': Diffuse field (sound from all directions) |
'free'
|
Returns:
| Name | Type | Description |
|---|---|---|
NDArrayReal |
NDArrayReal
|
Sharpness values in acum, one per channel. Shape: (n_channels,) |
Raises:
| Type | Description |
|---|---|
ValueError
|
If the signal sampling rate is not supported by the algorithm. |
Examples:
>>> import wandas as wd
>>> signal = wd.read("constant_tone.wav")
>>> sharpness = signal.sharpness_din_st(weighting="din", field_type="free")
>>> print(f"Steady-state sharpness: {sharpness[0]:.2f} acum")
Notes
- Sharpness is measured in acum (acum = 1 when the sound has the
same sharpness as a 2 kHz narrow-band noise at 60 dB SPL) - The calculation uses MoSQITo's implementation of DIN 45692 - Output is a single value per channel, suitable for stationary signals - For multi-channel signals, sharpness is calculated per channel
References
.. [1] DIN 45692:2009, "Measurement technique for the simulation of the auditory sensation of sharpness"
Source code in wandas/frames/mixins/channel_processing_mixin.py
1323 1324 1325 1326 1327 1328 1329 1330 1331 1332 1333 1334 1335 1336 1337 1338 1339 1340 1341 1342 1343 1344 1345 1346 1347 1348 1349 1350 1351 1352 1353 1354 1355 1356 1357 1358 1359 1360 1361 1362 1363 1364 1365 1366 1367 1368 1369 1370 | |
wandas.frames.mixins.channel_transform_mixin.ChannelTransformMixin
¶
Mixin providing methods related to frequency transformations.
This mixin provides operations related to frequency analysis and transformations such as FFT, STFT, and Welch method.
Source code in wandas/frames/mixins/channel_transform_mixin.py
186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 797 798 799 | |
Functions¶
cepstrum(n_fft=None, window='hann', floor=1e-12)
¶
Calculate the normalized real cepstrum of each channel.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
n_fft
|
int | None
|
int, optional. FFT size. |
None
|
window
|
str
|
str, default="hann". SciPy window name applied to the analysed samples. |
'hann'
|
floor
|
float
|
float, default=1e-12. Positive finite floor applied to normalized magnitude before |
1e-12
|
Returns:
| Name | Type | Description |
|---|---|---|
CepstralFrame |
CepstralFrame
|
New lazy real coefficients with dimensions
|
Raises:
| Type | Description |
|---|---|
TypeError
|
If the input is complex or a parameter has the wrong type. |
ValueError
|
If |
Notes
The method only builds a Dask graph. Accessing data, calling
compute(), or plotting materializes the coefficients.
Examples:
>>> cepstrum = frame.cepstrum(n_fft=2048, window="hann")
>>> envelope = cepstrum.lifter(0.002).to_spectral_envelope()
Source code in wandas/frames/mixins/channel_transform_mixin.py
198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 | |
fft(n_fft=None, window='hann')
¶
Calculate a one-sided peak-amplitude FFT spectrum.
The signal is truncated to n_fft, the window is applied to the
remaining min(n_samples, n_fft) samples, and the result is
zero-padded to n_fft and normalized by that window's coherent gain,
so a short signal keeps its amplitude. The window length is stored as
SpectralFrame.win_length. Values retain each channel's
physical unit. Positive-frequency bins other than Nyquist are doubled,
so an on-bin sinusoid's magnitude equals its peak amplitude. Graph
construction remains lazy.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
n_fft
|
int | None
|
Number of FFT points. By default, use the current sample count exactly. |
None
|
window
|
str
|
Window type. Default is "hann". |
'hann'
|
Returns:
| Type | Description |
|---|---|
SpectralFrame
|
A lazy SpectralFrame containing complex peak-amplitude values. |
Source code in wandas/frames/mixins/channel_transform_mixin.py
292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 | |
welch(n_fft=2048, hop_length=None, win_length=None, window='hann', average='mean')
¶
Calculate a Welch-averaged one-sided peak-amplitude spectrum.
Segment power spectra are averaged and converted to peak amplitude.
Values retain each channel's physical unit and are not power spectral
density or expressed per hertz. SpectralFrame.dB therefore uses the
amplitude rule 20 * log10(amplitude / channel_ref). Graph
construction remains lazy.
With average="mean", a long channel is passed to SciPy in runs of
whole segments and the run means are weighted by their segment counts,
so peak memory stays near the input size. average="median" needs all
segments at once and uses several times the input size.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
n_fft
|
int
|
Number of FFT points. Default is 2048. |
2048
|
hop_length
|
int | None
|
Number of samples between frames.
Default is |
None
|
win_length
|
int | None
|
Window length. Default is n_fft. |
None
|
window
|
str
|
Window type. Default is "hann". |
'hann'
|
average
|
Literal['mean', 'median']
|
Method for averaging segments. Default is "mean". |
'mean'
|
Returns:
| Type | Description |
|---|---|
SpectralFrame
|
A lazy SpectralFrame containing real peak-amplitude values. |
Source code in wandas/frames/mixins/channel_transform_mixin.py
366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 | |
noct_spectrum(fmin=25, fmax=20000, n=3, G=10, fr=1000)
¶
Calculate N-octave band spectrum.
Each output value is the RMS amplitude in one fractional-octave band
and retains the input channel's physical unit. NOctFrame.dB applies
20 * log10(band_rms / channel_ref).
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
fmin
|
float
|
Minimum center frequency (Hz). Default is 25 Hz. |
25
|
fmax
|
float
|
Maximum center frequency (Hz). Default is 20000 Hz. |
20000
|
n
|
int
|
Band division (1: octave, 3: 1/3 octave). Default is 3. |
3
|
G
|
int
|
Exact center-frequency ratio convention. Use 10 for base
|
10
|
fr
|
int
|
Reference frequency (Hz). Default is 1000 Hz. |
1000
|
Returns:
| Type | Description |
|---|---|
NOctFrame
|
A lazy NOctFrame containing per-band RMS amplitudes. |
Source code in wandas/frames/mixins/channel_transform_mixin.py
460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 | |
stft(n_fft=2048, hop_length=None, win_length=None, window='hann', *, allow_sparse=False, detrend=None)
¶
Calculate a one-sided peak-amplitude Short-Time Fourier Transform.
Each time frame is normalized by its window's coherent gain. Complex values retain the input physical unit; an on-bin sinusoid's magnitude is its peak amplitude.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
n_fft
|
int
|
Number of FFT points. Default is 2048. |
2048
|
hop_length
|
int | None
|
Number of samples between frames.
Default is |
None
|
win_length
|
int | None
|
Window length. Default is n_fft. |
None
|
window
|
str
|
Window type. Default is "hann". |
'hann'
|
allow_sparse
|
bool
|
Permit hop_length > win_length for sparse display sampling. Default is False; sparse results cannot reconstruct the full signal. |
False
|
detrend
|
Literal['constant', 'linear'] | None
|
Trend removed from each frame's samples before windowing:
None (default), |
None
|
Returns:
| Type | Description |
|---|---|
SpectrogramFrame
|
SpectrogramFrame containing STFT results |
Source code in wandas/frames/mixins/channel_transform_mixin.py
524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 | |
coherence(n_fft=2048, hop_length=None, win_length=None, window='hann', detrend='constant')
¶
Calculate typed magnitude-squared coherence for every channel pair.
The result is a :class:CoherenceFrame with flattened (pair,
frequency) storage and output-major/input-minor pair order. Its real
raw values are dimensionless and lie in [0, 1]; NaN is retained
for undefined zero-energy bins. This mathematical result contract is
enforced by the numerical operation, not by scanning values in the
:class:CoherenceFrame constructor. Pair roles, source identity,
domains, and row order are carried by immutable typed state, not labels
or operation history. See the spectral numerical contracts for the
canonical mathematical definition.
Sampling rate and user metadata are preserved. Each pair's
source_time_offset is derived from its input-role source offset, and
input calibration is consumed before the pairwise operation. Constructing
the result remains Dask-lazy; accessing data or plotting is the
materialization boundary. Invalid spectral parameters or structural
Frame state raise an actionable TypeError or ValueError instead
of being silently coerced.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
n_fft
|
int
|
Number of FFT points. Default is 2048. |
2048
|
hop_length
|
int | None
|
Number of samples between frames.
Default is |
None
|
win_length
|
int | None
|
Window length. Default is n_fft. |
None
|
window
|
str
|
Window type. Default is "hann". |
'hann'
|
detrend
|
str
|
Detrend method. Options: "constant", "linear", None. |
'constant'
|
Returns:
| Type | Description |
|---|---|
CoherenceFrame
|
CoherenceFrame whose public single-pair shape is |
CoherenceFrame
|
whose multi-pair shape is |
CoherenceFrame
|
for the quantity-specific raw values. |
Example
coherence = frame.coherence(n_fft=1024, window="hann")
Source code in wandas/frames/mixins/channel_transform_mixin.py
611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 | |
csd(n_fft=2048, hop_length=None, win_length=None, window='hann', detrend='constant', scaling='spectrum', average='mean')
¶
Calculate a typed cross-spectral density matrix.
The result is a :class:CrossSpectralFrame with flattened
(pair, frequency) storage and output-major/input-minor pair order.
Each raw complex row stores P_out_in = conj(X_input) * X_output;
pair domains provide the unit and reference, with /Hz included for
scaling="density". Pair roles and domains are immutable typed state;
labels and operation history are display/provenance views only. See the
spectral numerical contracts for the canonical definition and scaling.
Sampling rate and user metadata are preserved. Pair
source_time_offset uses the input-role source offset, and input
calibration is consumed before constructing output metadata. The result
stays Dask-lazy until data, a property, or a plot is materialized.
Invalid spectral parameters or domain/shape violations raise an actionable
TypeError or ValueError. Use the quantity-specific magnitude,
phase, and level_db properties; pairwise A-weighting is rejected.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
n_fft
|
int
|
Number of FFT points. Default is 2048. |
2048
|
hop_length
|
int | None
|
Number of samples between frames.
Default is |
None
|
win_length
|
int | None
|
Window length. Default is n_fft. |
None
|
window
|
str
|
Window type. Default is "hann". |
'hann'
|
detrend
|
str
|
Detrend method. Options: "constant", "linear", None. |
'constant'
|
scaling
|
Literal['spectrum', 'density']
|
Scaling method. Options: "spectrum", "density". |
'spectrum'
|
average
|
Literal['mean', 'median']
|
Method for averaging segments. Default is "mean". |
'mean'
|
Returns:
| Type | Description |
|---|---|
CrossSpectralFrame
|
CrossSpectralFrame whose public single-pair shape is |
CrossSpectralFrame
|
and whose multi-pair shape is |
Example
spectrum = frame.csd(n_fft=1024, scaling="density")
Source code in wandas/frames/mixins/channel_transform_mixin.py
672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 | |
transfer_function(n_fft=2048, hop_length=None, win_length=None, window='hann', detrend='constant', scaling='spectrum', average='mean')
¶
Calculate the canonical typed output/input transfer-function matrix.
The v2 result is a :class:TransferFunctionFrame with flattened
(pair, frequency) storage and output-major/input-minor pair order. It
stores H_out_in = P_out_in / P_in_in and carries the denominator
definition, pair roles, unit/reference domain, and row order as immutable
typed state. Labels and operation history do not define its meaning.
See the spectral numerical contracts for the canonical formulas.
Sampling rate and user metadata are preserved. Pair
source_time_offset uses the input-role source offset, and input
calibration is consumed before output metadata is derived. Construction
remains Dask-lazy; accessing data, a property, or a plot materializes the
requested values. Invalid spectral parameters or shape/domain violations
raise an actionable TypeError or ValueError. gain_db is
available only after selecting dimensionless pairs; transfer_level_db
uses each pair's explicit reference ratio. Pairwise A-weighting is
rejected.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
n_fft
|
int
|
Number of FFT points. Default is 2048. |
2048
|
hop_length
|
int | None
|
Number of samples between frames.
Default is |
None
|
win_length
|
int | None
|
Window length. Default is n_fft. |
None
|
window
|
str
|
Window type. Default is "hann". |
'hann'
|
detrend
|
str
|
Detrend method. Options: "constant", "linear", None. |
'constant'
|
scaling
|
Literal['spectrum', 'density']
|
Scaling method. Options: "spectrum", "density". |
'spectrum'
|
average
|
Literal['mean', 'median']
|
Method for averaging segments. Default is "mean". |
'mean'
|
Returns:
| Type | Description |
|---|---|
TransferFunctionFrame
|
TransferFunctionFrame whose public single-pair shape is |
TransferFunctionFrame
|
and whose multi-pair shape is |
Example
transfer = frame.transfer_function(n_fft=1024, scaling="spectrum")
Source code in wandas/frames/mixins/channel_transform_mixin.py
736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 797 798 799 | |
wandas.frames.mixins.spectral_properties_mixin.SpectralPropertiesMixin
¶
Shared magnitude, phase, squared-magnitude, and level properties.
Host classes must provide data (computed array),
_data (Dask array), _channel_metadata, and freqs.
The operation that created the host defines the stored quantity and unit.
NumPy properties use the same channel-axis convention as data: a
single-channel SpectralFrame returns (frequency,) and a
single-channel SpectrogramFrame returns (frequency, time); multiple
channels retain a leading channel axis. Plotting restores that axis only at
its boundary when it needs channel-first input. dBA uses the documented
internal _data.ndim (2 for spectra, 3 for spectrograms) to locate the
public frequency axis.
Source code in wandas/frames/mixins/spectral_properties_mixin.py
52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 | |
Attributes¶
magnitude
property
¶
Absolute magnitude of the stored spectral quantity.
phase
property
¶
Phase angles in radians.
power
property
¶
Squared magnitude, a compatibility property that is not a PSD.
dB
property
¶
Magnitude level: 20 * log10(magnitude / channel_ref).
For the canonical FFT, STFT, and Welch amplitude quantities, this is an amplitude level.
dBA
property
¶
A-weighted magnitude level relative to each channel reference.
For the canonical FFT, STFT, and Welch amplitude quantities, this is an A-weighted amplitude level.