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ファミリーを提供します。
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
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 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 | |
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
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 | |
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
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 | |
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
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 | |
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
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 | |
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
|
str
|
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
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 | |
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. |
{}
|
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")
Source code in wandas/frames/channel.py
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 | |
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
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 | |
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 |
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
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 | |
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
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 | |
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
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 | |
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)
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
|
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")
>>> # 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
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 | |
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
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 | |
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
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 | |
to_wav(path, format=None)
¶
Save the audio data to a WAV file.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
path
|
str | Path
|
Path to save the file. |
required |
format
|
str | None
|
File format. If None, determined from file extension. |
None
|
Source code in wandas/frames/channel.py
1597 1598 1599 1600 1601 1602 1603 1604 1605 1606 | |
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
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 | |
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
|
str
|
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
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 | |
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
|
str
|
|
'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
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 | |
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
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 | |
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'
|
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 | |
Attributes¶
n_fft = n_fft
instance-attribute
¶
window = window
instance-attribute
¶
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=())
¶
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
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 | |
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
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 | |
ifft()
¶
Invert Wandas FFT normalization to a windowed time-domain signal.
For a spectrum returned by ChannelFrame.fft() with matching stored
n_fft and window, this returns the truncated-or-zero-padded
analysis input multiplied by the FFT window. 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
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 | |
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 is read from this
SpectralFrame.n_fft; it is not a new caller-supplied parameter.
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
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 | |
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
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 | |
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
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 | |
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
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 | |
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
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 | |
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
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 | |
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
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 | |
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. The truthful legacy v1 value uses
denominator_role="output" and is never relabeled as canonical v2.
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
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 | |
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
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 | |
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 |
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
|
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
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 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 | |
Attributes¶
n_fft = n_fft
instance-attribute
¶
hop_length = hop_length
instance-attribute
¶
win_length = resolved_win_length
instance-attribute
¶
window = window
instance-attribute
¶
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. |
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=())
¶
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
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 | |
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
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 | |
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
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 | |
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
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 | |
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
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 | |
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
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 | |
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
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 | |
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
506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 | |
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
548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 | |
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
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 | |
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)
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
|
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
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 | |
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.
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. |
()
|
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
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 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 | |
Attributes¶
n_fft
property
¶
Return the immutable circular period of the complete cepstrum.
window
property
¶
Return the immutable originating analysis-window name.
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=())
¶
Source code in wandas/frames/cepstral.py
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 | |
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
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 | |
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
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 | |
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
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 | |
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 |
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
|
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
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 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 | |
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.
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=())
¶
Source code in wandas/frames/cepstrogram.py
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 | |
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
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 | |
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
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 | |
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
392 393 394 395 396 397 398 399 | |
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
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 | |
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
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 | |
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.
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
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 | |
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
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 | |
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 - (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
22 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 319 320 321 322 323 324 325 326 327 328 | |
Attributes¶
bark_axis = normalized_bark_axis.copy()
instance-attribute
¶
overlap = overlap
instance-attribute
¶
data
property
¶
Returns the computed data without squeezing.
For RoughnessFrame, even mono signals have 2D shape (47, n_time) so we don't squeeze the channel dimension.
Returns:
| Name | Type | Description |
|---|---|---|
NDArrayReal |
NDArrayReal
|
Computed data array. |
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
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 | |
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
228 229 230 231 232 233 234 235 236 237 | |
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
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 | |
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 | |
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.
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. |
{}
|
Returns:
| Type | Description |
|---|---|
Any
|
New frame with the custom function applied. |
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 | |
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 |
Source code in wandas/frames/mixins/channel_processing_mixin.py
236 237 238 239 240 241 242 243 244 245 246 247 248 249 | |
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 |
Source code in wandas/frames/mixins/channel_processing_mixin.py
251 252 253 254 255 256 257 258 259 260 261 262 263 264 | |
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
|
Filter order. Default is 4. |
4
|
Returns:
| Type | Description |
|---|---|
T_Processing
|
New ChannelFrame after filter application |
Source code in wandas/frames/mixins/channel_processing_mixin.py
266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 | |
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
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 | |
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
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 | |
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
366 367 368 369 370 371 372 373 374 375 376 377 378 379 | |
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
381 382 383 384 385 386 387 388 389 | |
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
391 392 393 394 395 396 397 398 399 400 401 402 | |
sum()
¶
Sum all channels.
Returns:
| Type | Description |
|---|---|
T_Processing
|
A new ChannelFrame with summed signal. |
Source code in wandas/frames/mixins/channel_processing_mixin.py
442 443 444 445 446 447 448 449 | |
mean()
¶
Average all channels.
Returns:
| Type | Description |
|---|---|
T_Processing
|
A new ChannelFrame with averaged signal. |
Source code in wandas/frames/mixins/channel_processing_mixin.py
451 452 453 454 455 456 457 458 | |
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
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 | |
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
516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 | |
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
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 | |
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
|
str | None
|
Implemented frequency-weighting curve: |
'Z'
|
time_weighting
|
str
|
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
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 | |
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
699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 | |
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
722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 | |
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
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 | |
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
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 | |
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
|
str
|
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
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 | |
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
|
str
|
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
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 | |
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
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 | |
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
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 | |
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
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 | |
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
|
str
|
str, default="din". Weighting type for sharpness calculation. Options: - 'din': DIN 45692 method - 'aures': Aures method - 'bismarck': Bismarck method - 'fastl': Fastl method |
'din'
|
field_type
|
str
|
str, 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
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 | |
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
|
str
|
str, default="din". Weighting type for sharpness calculation. Options: - 'din': DIN 45692 method - 'aures': Aures method - 'bismarck': Bismarck method - 'fastl': Fastl method |
'din'
|
field_type
|
str
|
str, 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
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 | |
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
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 | |
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 before the FFT. |
'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
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 | |
fft(n_fft=None, window='hann')
¶
Calculate a one-sided peak-amplitude FFT spectrum.
The signal is truncated or zero-padded to n_fft, windowed, and
normalized by the window's coherent gain. 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
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 | |
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.
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
|
str
|
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
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 | |
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
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 | |
stft(n_fft=2048, hop_length=None, win_length=None, window='hann')
¶
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'
|
Returns:
| Type | Description |
|---|---|
SpectrogramFrame
|
SpectrogramFrame containing STFT results |
Source code in wandas/frames/mixins/channel_transform_mixin.py
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 | |
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 n_fft//4. |
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
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 | |
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 n_fft//4. |
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
|
str
|
Scaling method. Options: "spectrum", "density". |
'spectrum'
|
average
|
str
|
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
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 | |
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; the
released v1 denominator contract is replayed separately by the v1 Recipe
handler. 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 n_fft//4. |
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
|
str
|
Scaling method. Options: "spectrum", "density". |
'spectrum'
|
average
|
str
|
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
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 | |
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.