Frames Module / フレームモジュール¶
The wandas.frames module provides various data frame classes for manipulating and representing audio data.
wandas.frames モジュールは、オーディオデータの操作と表現のための様々なデータフレームクラスを提供します。
ChannelFrame¶
ChannelFrame is the basic frame for handling time-domain waveform data. ChannelFrameは時間領域の波形データを扱うための基本的なフレームです。
Frame annotations are updated immutably with with_label(), with_metadata(),
or with_channel_extra(). Use with_source_time_offset()
for portable source-time intent and rename_channels() on any Frame family.
Direct mutation is unsupported; public state getters return detached snapshots.
String channel lookup always means a channel name. Use frame.channels.by_id()
for explicit stable-ID lookup.
ChannelCalibration.with_unit() preserves the factor and resets ref to the
new unit's default. Chain .with_unit(...).with_ref(...) for a custom reference.
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
258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 839 840 841 842 843 844 845 846 847 848 849 850 851 852 853 854 855 856 857 858 859 860 861 862 863 864 865 866 867 868 869 870 871 872 873 874 875 876 877 878 879 880 881 882 883 884 885 886 887 888 889 890 891 892 893 894 895 896 897 898 899 900 901 902 903 904 905 906 907 908 909 910 911 912 913 914 915 916 917 918 919 920 921 922 923 924 925 926 927 928 929 930 931 932 933 934 935 936 937 938 939 940 941 942 943 944 945 946 947 948 949 950 951 952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967 968 969 970 971 972 973 974 975 976 977 978 979 980 981 982 983 984 985 986 987 988 989 990 991 992 993 994 995 996 997 998 999 1000 1001 1002 1003 1004 1005 1006 1007 1008 1009 1010 1011 1012 1013 1014 1015 1016 1017 1018 1019 1020 1021 1022 1023 1024 1025 1026 1027 1028 1029 1030 1031 1032 1033 1034 1035 1036 1037 1038 1039 1040 1041 1042 1043 1044 1045 1046 1047 1048 1049 1050 1051 1052 1053 1054 1055 1056 1057 1058 1059 1060 1061 1062 1063 1064 1065 1066 1067 1068 1069 1070 1071 1072 1073 1074 1075 1076 1077 1078 1079 1080 1081 1082 1083 1084 1085 1086 1087 1088 1089 1090 1091 1092 1093 1094 1095 1096 1097 1098 1099 1100 1101 1102 1103 1104 1105 1106 1107 1108 1109 1110 1111 1112 1113 1114 1115 1116 1117 1118 1119 1120 1121 1122 1123 1124 1125 1126 1127 1128 1129 1130 1131 1132 1133 1134 1135 1136 1137 1138 1139 1140 1141 1142 1143 1144 1145 1146 1147 1148 1149 1150 1151 1152 1153 1154 1155 1156 1157 1158 1159 1160 1161 1162 1163 1164 1165 1166 1167 1168 1169 1170 1171 1172 1173 1174 1175 1176 1177 1178 1179 1180 1181 1182 1183 1184 1185 1186 1187 1188 1189 1190 1191 1192 1193 1194 1195 1196 1197 1198 1199 1200 1201 1202 1203 1204 1205 1206 1207 1208 1209 1210 1211 1212 1213 1214 1215 1216 1217 1218 1219 1220 1221 1222 1223 1224 1225 1226 1227 1228 1229 1230 1231 1232 1233 1234 1235 1236 1237 1238 1239 1240 1241 1242 1243 1244 1245 1246 1247 1248 1249 1250 1251 1252 1253 1254 1255 1256 1257 1258 1259 1260 1261 1262 1263 1264 1265 1266 1267 1268 1269 1270 1271 1272 1273 1274 1275 1276 1277 1278 1279 1280 1281 1282 1283 1284 1285 1286 1287 1288 1289 1290 1291 1292 1293 1294 1295 1296 1297 1298 1299 1300 1301 1302 1303 1304 1305 1306 1307 1308 1309 1310 1311 1312 1313 1314 1315 1316 1317 1318 1319 1320 1321 1322 1323 1324 1325 1326 1327 1328 1329 1330 1331 1332 1333 1334 1335 1336 1337 1338 1339 1340 1341 1342 1343 1344 1345 1346 1347 1348 1349 1350 1351 1352 1353 1354 1355 1356 1357 1358 1359 1360 1361 1362 1363 1364 1365 1366 1367 1368 1369 1370 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 | |
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 RMS (Root Mean Square) value 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 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 |
NDArrayReal
|
for each channel. |
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 |
NDArrayReal
|
for each channel. All-zero channels yield 1.0. |
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
|
Compatibility/debug pointer to the immediate prior frame; not the provenance source of truth. |
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
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 | |
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
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 | |
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.
Source code in wandas/frames/channel.py
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 | |
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
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 | |
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
|
Note
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
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 | |
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
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 | |
rms_plot(ax=None, title=None, overlay=True, Aw=False, **kwargs)
¶
Generate an RMS plot.
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 A-weighting. |
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 [V]", alpha=0.8, color="blue")
Source code in wandas/frames/channel.py
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 | |
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 |
list[Figure] | None
|
objects created for each channel. The list length equals the number of |
list[Figure] | None
|
channels in the frame. |
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
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 | |
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
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 | |
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
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 | |
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.
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 audio data. |
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
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 | |
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
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 | |
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 the data (lazy loading). |
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
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 | |
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
1434 1435 1436 1437 1438 1439 1440 1441 1442 1443 | |
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 |
Example
import wandas as wd cf = wd.load("audio_analysis.wdf")
Source code in wandas/frames/channel.py
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 | |
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 | ChannelFrame
|
Data to add as a new channel. Can be: - numpy array (1D or 2D) - dask array (1D or 2D) - ChannelFrame (channels will be added) |
required |
label
|
str | None
|
Label for the new channel. If None, generates a default label.
When data is a ChannelFrame, acts as a prefix: each channel in
the input frame is renamed to |
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 raw numpy or dask input. If None, raw input uses 0.0. When data is a ChannelFrame, offsets are taken from that frame and this argument must be None. |
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 type is not supported. |
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")
>>> # Add another ChannelFrame's channels
>>> cf2 = wd.read("audio2.wav")
>>> cf_combined = cf.add_channel(cf2)
Source code in wandas/frames/channel.py
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 | |
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 |
ChannelFrame
|
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
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 | |
Use the per-channel
calibration learning app to configure known conversion factors without
modifying the source frame. Calibrated physical values are available from
frame.data as a NumPy array; users do not need to manage the internal array
backend.
For a recorded reference event, use exactly one known scalar target:
microphone_calibration = microphone_reference.derive_calibration(
target_level=94.0,
unit="Pa",
)
acceleration_calibration = acceleration_reference.derive_calibration(
target_rms=1.0,
unit="m/s^2",
)
measurement = measurement.with_calibration(
{**microphone_calibration, **acceleration_calibration}
)
The scalar is broadcast to every channel in that reference event. References
with different targets, units, or recording times are separate events whose
label mappings can be combined. Derivation always uses the current
frame.data/frame.rms, including an existing factor, and does not inspect or
change operation history. It requires unique non-empty labels. The reference
and measurement recording chain—including amplifier gain—must represent the
same physical scale.
get_channel(..., validate_query_keys: bool = True) parameter¶
- validate_query_keys: When
True(default), dict-stylequeryarguments are validated against the known channel metadata fields and any existingextrakeys. Unknown keys raiseKeyErrorwith the message "Unknown channel metadata key". Set toFalseto skip this pre-validation and allow queries that reference keys not present on the model; in that case, normal matching proceeds and a no-match will raise the usualKeyErrorfor no results.
Source-time offsets and index-wise operations¶
source_time_offset records where each channel's local sample axis starts on
the original source timeline. Binary frame operators such as frame_a + frame_b
do not use this value for automatic alignment. They operate on the current array
indices after verifying that sampling rate, channel count, and shape match.
Different source_time_offset values are allowed. The result inherits the left
operand's source_time_offset, so frame_a + frame_b carries frame_a's
source timeline. channel_difference() follows the same index-wise principle
within one frame and preserves the input channel offsets.
When a workflow needs source-time alignment, trim or otherwise align frames explicitly before applying binary operators. A dedicated source-time alignment API may be added separately in the future.
SpectralFrame¶
SpectralFrame is a frame for handling frequency-domain data. SpectralFrameは周波数領域のデータを扱うためのフレームです。
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¶
data : 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)
sampling_rate : float
The sampling rate of the original time-domain signal in Hz.
n_fft : 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 n_fft // 2 + 1 frequency bins.
window : str, default="hann"
The window function used in the FFT.
label : str, optional
A label for the frame.
metadata : dict, optional
Additional metadata for the frame.
lineage : LineageNode, optional
Constructor override for the runtime lineage. When omitted, a source node is
created. operation_history is its public derived projection.
channel_metadata : list[ChannelMetadata], optional
Metadata for each channel in the frame.
previous : BaseFrame, optional
Compatibility/debug pointer to the immediate prior frame; not the
provenance source of truth.
Attributes¶
magnitude : NDArrayReal The magnitude spectrum of the data. phase : NDArrayReal The phase spectrum in radians. unwrapped_phase : NDArrayReal The unwrapped phase spectrum in radians. power : NDArrayReal The power spectrum (magnitude squared). dB : NDArrayReal The spectrum in decibels relative to channel reference values. dBA : NDArrayReal The A-weighted spectrum in decibels. freqs : 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 magnitude 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
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 | |
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¶
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¶
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
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 | |
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¶
plot_type : 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 ax : matplotlib.axes.Axes, optional Axes to plot on. If None, creates new axes. title : str, optional Title for the plot. If None, uses the frame label. overlay : bool, default=False Whether to overlay all channels on a single plot (True) or create separate subplots for each channel (False). xlabel : str, optional Label for the x-axis. If None, uses default "Frequency [Hz]". ylabel : str, optional Label for the y-axis. If None, uses default based on data type. alpha : float, default=1.0 Transparency level for the plot lines (0.0 to 1.0). xlim : tuple[float, float], optional Limits for the x-axis as (min, max) tuple. ylim : tuple[float, float], optional Limits for the y-axis as (min, max) tuple. Aw : bool, default=False Whether to apply A-weighting to the data. **kwargs : dict Additional matplotlib Line2D parameters (e.g., color, linewidth, linestyle).
Returns¶
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
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 | |
ifft()
¶
Compute the Inverse Fast Fourier Transform (IFFT) to return to time domain.
This method transforms the frequency-domain data back to the time domain using the inverse FFT operation. The window function used in the forward FFT is taken into account to ensure proper reconstruction.
Returns¶
ChannelFrame A new ChannelFrame containing the time-domain signal.
Source code in wandas/frames/spectral.py
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 | |
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.
Parameters¶
fmin : float Lower frequency bound in Hz. fmax : float Upper frequency bound in Hz. n : int, default=3 Number of bands per octave (e.g., 3 for third-octave bands). G : int, default=10 Reference band number. fr : int, default=1000 Reference frequency in Hz.
Returns¶
NOctFrame A new NOctFrame containing the N-octave band spectrum.
Raises¶
ValueError If the sampling rate is not 48000 Hz.
Source code in wandas/frames/spectral.py
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 | |
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¶
plot_type : str, default="matrix" Type of matrix plot to create. **kwargs : dict Additional plot parameters: - vmin, vmax: Color scale limits - cmap: Colormap name - title: Plot title
Returns¶
Union[Axes, Iterator[Axes]] The matplotlib axes containing the plot.
Source code in wandas/frames/spectral.py
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 | |
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
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 | |
CepstralFrame¶
CepstralFrame represents a normalized real cepstrum on a quefrency axis. Start with
the cepstral analysis guide for the typed
ChannelFrame -> CepstralFrame -> SpectralFrame workflow.
CepstralFrameは、ケフレンシー軸上の正規化された実ケプストラムを表します。型付きの
ChannelFrame -> CepstralFrame -> SpectralFrameワークフローは
ケプストラム解析ガイドを参照してください。
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¶
data : dask.array.Array
Real coefficients shaped (quefrency,) or
(channels, quefrency).
sampling_rate : float
Sampling rate in Hz that defines the quefrency-bin spacing.
n_fft : int
Positive FFT size of the complete cepstrum. Sliced data may contain fewer
bins but never more than this value.
window : str, default="hann"
Window used by the originating cepstrum analysis.
label : str, optional
Human-readable frame label.
metadata : dict, optional
User and recording metadata, copied on construction.
channel_metadata : sequence, optional
Metadata aligned with the channel axis.
channel_ids : list[str], optional
Stable identifiers aligned with the channel axis.
previous : BaseFrame, optional
Compatibility/debug pointer to the immediate prior frame.
source_time_offset : float or sequence, default=0.0
Per-channel source timeline offsets, preserved across domain changes.
lineage : LineageNode, optional
Authoritative runtime semantic lineage.
operation_history_prefix : sequence, default=()
Persisted display history for a new source frame.
Raises¶
TypeError
If coefficients are complex, n_fft is not integral, or window is
not a non-empty string.
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 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 | |
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
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 | |
lifter(cutoff, mode='low')
¶
Keep low- or high-quefrency coefficients.
Parameters¶
cutoff : float
Positive quefrency boundary in seconds. It must reach at least one
represented bin and remain below half the complete cepstrum.
mode : {"low", "high"}, default="low"
"low" keeps the smooth-envelope region; "high" keeps the
complementary fine structure.
Returns¶
CepstralFrame A new lazy frame with the same axes and metadata.
Raises¶
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
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 | |
to_spectral_envelope()
¶
Convert a complete real cepstrum to a smooth spectral envelope.
Returns¶
SpectralFrame
New lazy complex-valued frequency data with zero phase, the original
n_fft and window, and preserved metadata and source-time offsets.
Raises¶
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
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 | |
plot(plot_type='quefrency', ax=None, *, title=None, xlabel='Quefrency [s]', ylabel='Real cepstrum', **kwargs)
¶
Plot real coefficients against quefrency.
Parameters¶
plot_type : str, default="quefrency"
Only "quefrency" is supported.
ax : matplotlib.axes.Axes, optional
Existing axes. A new figure and axes are created when omitted.
title : str, optional
Plot title; defaults to the frame label.
xlabel, ylabel : str
Axis labels.
**kwargs : Any
Keyword arguments passed to Axes.plot.
Returns¶
matplotlib.axes.Axes Axes containing one line per channel.
Raises¶
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
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 | |
CepstrogramFrame¶
CepstrogramFrame represents a real cepstrum at every STFT time frame with dimensions
(channel, quefrency, time). It is created by SpectrogramFrame.cepstrum() and can
reconstruct a time-varying spectral envelope. See the
cepstral analysis guide.
CepstrogramFrameは各STFT時間フレームの実ケプストラムを
(channel, quefrency, time)で表します。SpectrogramFrame.cepstrum()から生成し、
時間変化するスペクトル包絡を再構成できます。
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¶
data : dask.array.Array
Real coefficients shaped (quefrency, time) or
(channels, quefrency, time).
sampling_rate : float
Sampling rate in Hz defining both axis spacings.
n_fft : int
Positive FFT size of the complete cepstrum.
hop_length : int
Positive sample distance between adjacent time frames.
win_length : int, optional
Analysis-window length inherited from the source spectrogram. Defaults
to n_fft.
window : str, default="hann"
Analysis-window name inherited from the source spectrogram.
label : str, optional
Human-readable frame label.
metadata : dict, optional
User and recording metadata, copied on construction.
channel_metadata : sequence, optional
Metadata aligned with the channel axis.
channel_ids : list[str], optional
Stable identifiers aligned with the channel axis.
previous : BaseFrame, optional
Compatibility/debug pointer to the immediate prior frame.
source_time_offset : float or sequence, default=0.0
Per-channel source timeline offsets.
lineage : LineageNode, optional
Authoritative runtime semantic lineage.
operation_history_prefix : sequence, default=()
Persisted display history for a new source frame.
Raises¶
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 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 | |
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
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 | |
lifter(cutoff, mode='low')
¶
Keep low- or high-quefrency coefficients at every time frame.
Parameters¶
cutoff : float
Positive quefrency boundary in seconds. It must reach at least one
bin and remain below half of the complete cepstrum.
mode : {"low", "high"}, default="low"
"low" keeps the smooth-envelope region; "high" keeps the
complementary fine structure.
Returns¶
CepstrogramFrame New lazy coefficients with unchanged time and channel axes.
Raises¶
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
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 | |
to_spectral_envelope()
¶
Reconstruct a smooth magnitude spectrogram with zero phase.
Returns¶
SpectrogramFrame New lazy frequency-time data preserving the original STFT analysis state, channels, metadata, and source-time offsets.
Raises¶
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
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 | |
to_dataframe()
¶
Reject conversion because the frame has three semantic dimensions.
Raises¶
NotImplementedError Always raised. Materialize selected data when a tabular representation is required.
Source code in wandas/frames/cepstrogram.py
422 423 424 425 426 427 428 429 430 431 | |
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¶
plot_type : str, default="cepstrogram"
Only "cepstrogram" is supported.
ax : matplotlib.axes.Axes, optional
Existing axes for a single-channel frame. Multi-channel frames
create one axes per channel when omitted.
title : str, optional
Plot title prefix; defaults to the frame label.
xlabel, ylabel : str
Axis labels.
cmap : str, default="RdBu_r"
Matplotlib colormap for signed real coefficients.
qmin, qmax : float, optional
Display range on the quefrency axis in seconds.
vmin, vmax : float, optional
Shared color limits. When omitted, a symmetric robust range is
estimated from the displayed coefficients except the dominant
zero-quefrency row.
**kwargs : Any
Additional keyword arguments passed to Axes.pcolormesh.
Returns¶
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
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 | |
SpectrogramFrame¶
SpectrogramFrame is a frame for handling time-frequency domain (spectrogram) data. SpectrogramFrameは時間-周波数領域(スペクトログラム)のデータを扱うフレームです。
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¶
data : 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)
sampling_rate : float
The sampling rate of the original time-domain signal in Hz.
n_fft : int
The FFT size used to generate this spectrogram. The frequency dimension must
contain exactly n_fft // 2 + 1 bins.
hop_length : int
Number of samples between successive frames.
win_length : int, optional
The window length in samples. If None, defaults to n_fft.
window : str, default="hann"
The window function to use (e.g., "hann", "hamming", "blackman").
label : str, optional
A label for the frame.
metadata : dict, optional
Additional metadata for the frame.
lineage : LineageNode, optional
Constructor override for the runtime lineage. When omitted, a source node is
created. operation_history is its public derived projection.
channel_metadata : list[ChannelMetadata], optional
Metadata for each channel in the frame.
previous : BaseFrame, optional
Compatibility/debug pointer to the immediate prior frame; not the
provenance source of truth.
Attributes¶
magnitude : NDArrayReal The magnitude spectrogram. phase : NDArrayReal The phase spectrogram in radians. power : NDArrayReal The power spectrogram. dB : NDArrayReal The spectrogram in decibels relative to channel reference values. dBA : NDArrayReal The A-weighted spectrogram in decibels. n_frames : int Number of time frames. n_freq_bins : int Number of frequency bins. freqs : NDArrayReal The frequency axis values in Hz. times : 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 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 | |
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
¶
n_freq_bins
property
¶
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¶
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¶
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
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 | |
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¶
plot_type : str, default="spectrogram" Type of plot to create. ax : matplotlib.axes.Axes, optional Axes to plot on. If None, creates new axes. title : str, optional Title for the plot. If None, uses the frame label. cmap : str, default="jet" Colormap name for the spectrogram visualization. vmin : float, optional Minimum value for colormap scaling (dB). Auto-calculated if None. vmax : float, optional Maximum value for colormap scaling (dB). Auto-calculated if None. fmin : float, default=0 Minimum frequency to display (Hz). fmax : float, optional Maximum frequency to display (Hz). If None, uses Nyquist frequency. xlim : tuple[float, float], optional Time axis limits as (start_time, end_time) in seconds. ylim : tuple[float, float], optional Frequency axis limits as (min_freq, max_freq) in Hz. Aw : bool, default=False Whether to apply A-weighting to the spectrogram. overlay : bool, default=False Whether to overlay channels on a single axes. **kwargs : dict Additional keyword arguments passed to Matplotlib plotting methods.
Returns¶
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
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 | |
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¶
plot_type : str, default="spectrogram" Type of plot to create. ax : matplotlib.axes.Axes, optional Axes to plot on. If None, creates new axes. **kwargs : dict Additional keyword arguments passed to plot(). Accepts all parameters from plot() except Aw (which is set to True).
Returns¶
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
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 | |
cepstrum(floor=1e-12)
¶
Calculate a real cepstrum independently at every time frame.
Parameters¶
floor : float, default=1e-12 Positive finite floor applied to normalized STFT magnitude before taking the logarithm.
Returns¶
CepstrogramFrame
New lazy coefficients shaped (channel, quefrency, time). The
source FFT size, hop length, window state, channels, metadata, and
source-time offsets are preserved.
Raises¶
TypeError
If floor is not a real number.
ValueError
If floor is non-positive or non-finite.
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
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 | |
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¶
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
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 | |
get_frame_at(time_idx)
¶
Extract spectral data at a specific time frame.
Parameters¶
time_idx : int Index of the time frame to extract.
Returns¶
SpectralFrame A new SpectralFrame containing the spectral data at the specified time.
Raises¶
IndexError If time_idx is out of range.
Source code in wandas/frames/spectrogram.py
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 | |
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¶
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
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 | |
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¶
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
579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 | |
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¶
NotImplementedError Always raised as DataFrame conversion is not supported.
Source code in wandas/frames/spectrogram.py
627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 | |
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
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 | |
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
|
Reference to the previous frame in the processing chain. |
None
|
Returns:
| Type | Description |
|---|---|
SpectrogramFrame
|
A new SpectrogramFrame containing the NumPy data. |
Source code in wandas/frames/spectrogram.py
710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 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 | |
NOctFrame¶
NOctFrame is a frame class for octave-band analysis. NOctFrameはオクターブバンド解析のためのフレームクラスです。
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 energy or power in each frequency band, following standard acoustical band definitions.
Parameters¶
data : 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)
sampling_rate : float
The sampling rate of the original time-domain signal in Hz.
fmin : float, default=0
Lower frequency bound in Hz.
fmax : float, default=0
Upper frequency bound in Hz.
n : int, default=3
Number of bands per octave (e.g., 3 for third-octave bands).
G : int, default=10
Reference band number according to IEC 61260-1:2014.
fr : int, default=1000
Reference frequency in Hz, typically 1000 Hz for acoustic analysis.
label : str, optional
A label for the frame.
metadata : dict, optional
Additional metadata for the frame.
lineage : LineageNode, optional
Constructor override for the runtime lineage. When omitted, a source node is
created. operation_history is its public derived projection.
channel_metadata : list[ChannelMetadata], optional
Metadata for each channel in the frame.
previous : BaseFrame, optional
Compatibility/debug pointer to the immediate prior frame; not the
provenance source of truth.
Attributes¶
freqs : NDArrayReal The center frequencies of each band in Hz, calculated according to the standard fractional octave band definitions. dB : NDArrayReal The spectrum in decibels relative to channel reference values. dBA : NDArrayReal The A-weighted spectrum in decibels, applying frequency weighting for better correlation with perceived loudness. fmin : float Lower frequency bound in Hz. fmax : float Upper frequency bound in Hz. n : int Number of bands per octave. G : int Reference band number. fr : 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 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 | |
Attributes¶
n = n
instance-attribute
¶
G = G
instance-attribute
¶
fr = fr
instance-attribute
¶
fmin = fmin
instance-attribute
¶
fmax = fmax
instance-attribute
¶
dB
property
¶
Get the spectrum in decibels relative to each channel's reference value.
The reference value for each channel is specified in its metadata. A minimum value of -120 dB is enforced to avoid numerical issues.
Returns¶
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¶
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¶
NDArrayReal Array of center frequencies for each frequency band.
Raises¶
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
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 | |
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¶
plot_type : str, default="noct" Type of plot to create. The default "noct" type creates a step plot suitable for displaying N-octave band data. ax : matplotlib.axes.Axes, optional Axes to plot on. If None, creates new axes. title : str, optional Title for the plot. If None, uses a default title with band specification. overlay : bool, default=False Whether to overlay all channels on a single plot (True) or create separate subplots for each channel (False). xlabel : str, optional Label for the x-axis. If None, uses default "Center frequency [Hz]". ylabel : str, optional Label for the y-axis. If None, uses default based on data type. alpha : float, default=1.0 Transparency level for the plot lines (0.0 to 1.0). xlim : tuple[float, float], optional Limits for the x-axis as (min, max) tuple. ylim : tuple[float, float], optional Limits for the y-axis as (min, max) tuple. Aw : bool, default=False Whether to apply A-weighting to the data. **kwargs : dict Additional matplotlib Line2D parameters (e.g., color, linewidth, linestyle).
Returns¶
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
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 | |
RoughnessFrame¶
RoughnessFrame is a frame class for psychoacoustic roughness analysis results. RoughnessFrameは心理音響ラフネス解析結果のためのフレームクラスです。
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¶
data : 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.
sampling_rate : 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).
bark_axis : NDArrayReal
Bark frequency axis with 47 values from 0.5 to 23.5 Bark.
overlap : float
Overlap coefficient used in the calculation (0.0 to 1.0).
label : str, optional
Frame label. Defaults to "roughness_spec".
metadata : dict, optional
Additional metadata.
lineage : LineageNode, optional
Constructor override for the runtime lineage. When omitted, a source node is
created. operation_history is its public derived projection.
channel_metadata : list[ChannelMetadata], optional
Metadata for each channel.
previous : BaseFrame, optional
Compatibility/debug pointer to the immediate prior frame; not the
provenance source of truth.
Attributes¶
bark_axis : NDArrayReal Frequency axis in Bark scale. n_bark_bands : int Number of Bark bands (always 47). n_time_points : int Number of time points. time : NDArrayReal Time axis based on sampling rate. overlap : 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 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 | |
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¶
NDArrayReal Computed data array.
n_time_points
property
¶
Number of time points in the roughness time series.
Returns¶
int Number of time frames in the analysis.
time
property
¶
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
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 | |
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¶
NotImplementedError Always raised as DataFrame conversion is not supported.
Source code in wandas/frames/roughness.py
255 256 257 258 259 260 261 262 263 264 265 266 | |
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¶
ax : Axes, optional Matplotlib axes to plot on. If None, a new figure is created. title : str, optional Plot title. If None, a default title is used. cmap : str, default="viridis" Colormap name for the heatmap. vmin, vmax : float, optional Color scale limits. If None, automatic scaling is used. xlabel : str, default="Time [s]" Label for the x-axis. ylabel : str, default="Frequency [Bark]" Label for the y-axis. colorbar_label : str, default="Specific Roughness [Asper/Bark]" Label for the colorbar. **kwargs : Any Additional keyword arguments passed to pcolormesh.
Returns¶
Axes The matplotlib axes object containing the plot.
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
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 | |
Mixins¶
Mixins for extending frame functionality. フレームの機能を拡張するためのミックスインです。
ChannelProcessingMixin¶
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
60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 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 | |
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
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 | |
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
205 206 207 208 209 210 211 212 213 214 215 216 217 218 | |
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
220 221 222 223 224 225 226 227 228 229 230 231 232 233 | |
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
235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 | |
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
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 | |
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
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 | |
a_weighting()
¶
Apply A-weighting filter to the signal.
A-weighting adjusts the frequency response to approximate human auditory perception, according to the IEC 61672-1:2013 standard.
Returns:
| Type | Description |
|---|---|
T_Processing
|
New ChannelFrame containing the A-weighted signal |
Source code in wandas/frames/mixins/channel_processing_mixin.py
335 336 337 338 339 340 341 342 343 344 345 346 | |
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
348 349 350 351 352 353 354 355 356 | |
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
358 359 360 361 362 363 364 365 366 367 368 369 | |
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
409 410 411 412 413 414 415 416 | |
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
418 419 420 421 422 423 424 425 | |
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 |
Raises:
| Type | Description |
|---|---|
ValueError
|
If end time is earlier than start time |
Source code in wandas/frames/mixins/channel_processing_mixin.py
427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 | |
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
452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 | |
rms_trend(frame_length=2048, hop_length=512, dB=False, Aw=False)
¶
Compute the RMS trend of the signal.
This method calculates the root mean square value over a sliding window.
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
|
Whether to return RMS values in decibels. Default is False. |
False
|
Aw
|
bool
|
Whether to apply A-weighting. Default is False. |
False
|
Returns:
| Type | Description |
|---|---|
T_Processing
|
New ChannelFrame containing the RMS trend |
Source code in wandas/frames/mixins/channel_processing_mixin.py
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 | |
sound_level(freq_weighting='Z', time_weighting='Fast', dB=False)
¶
Compute a time-weighted RMS trend or sound pressure level.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
freq_weighting
|
str | None
|
Frequency weighting curve. Supported values are
|
'Z'
|
time_weighting
|
str
|
Time weighting characteristic. Supported values are
|
'Fast'
|
dB
|
bool
|
When |
False
|
Returns:
| Type | Description |
|---|---|
T_Processing
|
New ChannelFrame containing the weighted time series. |
Source code in wandas/frames/mixins/channel_processing_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 | |
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 |
T_Processing
|
source-time offsets preserved. |
Source code in wandas/frames/mixins/channel_processing_mixin.py
540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 | |
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
563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 | |
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
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 | |
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
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 | |
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. |
T_Processing
|
Each channel is processed independently. |
T_Processing
|
The output sampling rate is adjusted based on the loudness |
T_Processing
|
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
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 | |
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
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 | |
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. |
T_Processing
|
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
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 | |
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
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 | |
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
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 | |
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¶
weighting : str, default="din" Weighting type for sharpness calculation. Options: - 'din': DIN 45692 method - 'aures': Aures method - 'bismarck': Bismarck method - 'fastl': Fastl method field_type : str, default="free" Type of sound field. Options: - 'free': Free field (sound from a specific direction) - 'diffuse': Diffuse field (sound from all directions)
Returns¶
T_Processing New ChannelFrame containing sharpness time series in acum. The output sampling rate is approximately 500 Hz (2ms time steps).
Raises¶
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
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 | |
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¶
weighting : str, default="din" Weighting type for sharpness calculation. Options: - 'din': DIN 45692 method - 'aures': Aures method - 'bismarck': Bismarck method - 'fastl': Fastl method field_type : str, default="free" Type of sound field. Options: - 'free': Free field (sound from a specific direction) - 'diffuse': Diffuse field (sound from all directions)
Returns¶
NDArrayReal Sharpness values in acum, one per channel. Shape: (n_channels,)
Raises¶
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
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 | |
ChannelTransformMixin¶
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
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 | |
Functions¶
cepstrum(n_fft=None, window='hann', floor=1e-12)
¶
Calculate the normalized real cepstrum of each channel.
Parameters¶
n_fft : int, optional
FFT size. None uses the current sample count. Smaller values
truncate and larger values zero-pad the analysis input.
window : str, default="hann"
SciPy window name applied before the FFT.
floor : float, default=1e-12
Positive finite floor applied to normalized magnitude before log.
Returns¶
CepstralFrame
New lazy real coefficients with dimensions
(channel, quefrency). Channel metadata, IDs, user metadata,
sampling rate, and source-time offsets are preserved.
Raises¶
TypeError
If the input is complex or a parameter has the wrong type.
ValueError
If n_fft or floor is invalid.
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
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 | |
fft(n_fft=None, window='hann')
¶
Calculate Fast Fourier Transform (FFT).
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
n_fft
|
int | None
|
Number of FFT points. Default is the next power of 2 of the data length. |
None
|
window
|
str
|
Window type. Default is "hann". |
'hann'
|
Returns:
| Type | Description |
|---|---|
SpectralFrame
|
SpectralFrame containing FFT results |
Source code in wandas/frames/mixins/channel_transform_mixin.py
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 | |
welch(n_fft=2048, hop_length=None, win_length=None, window='hann', average='mean')
¶
Calculate power spectral density using Welch's method.
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'
|
average
|
str
|
Method for averaging segments. Default is "mean". |
'mean'
|
Returns:
| Type | Description |
|---|---|
SpectralFrame
|
SpectralFrame containing power spectral density |
Source code in wandas/frames/mixins/channel_transform_mixin.py
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 | |
noct_spectrum(fmin=25, fmax=20000, n=3, G=10, fr=1000)
¶
Calculate N-octave band spectrum.
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
|
Reference gain (dB). Default is 10 dB. |
10
|
fr
|
int
|
Reference frequency (Hz). Default is 1000 Hz. |
1000
|
Returns:
| Type | Description |
|---|---|
NOctFrame
|
NOctFrame containing N-octave band spectrum |
Source code in wandas/frames/mixins/channel_transform_mixin.py
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 | |
stft(n_fft=2048, hop_length=None, win_length=None, window='hann')
¶
Calculate Short-Time Fourier Transform.
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'
|
Returns:
| Type | Description |
|---|---|
SpectrogramFrame
|
SpectrogramFrame containing STFT results |
Source code in wandas/frames/mixins/channel_transform_mixin.py
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 | |
coherence(n_fft=2048, hop_length=None, win_length=None, window='hann', detrend='constant')
¶
Calculate magnitude squared coherence.
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 |
|---|---|
SpectralFrame
|
SpectralFrame containing magnitude squared coherence |
Source code in wandas/frames/mixins/channel_transform_mixin.py
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 | |
csd(n_fft=2048, hop_length=None, win_length=None, window='hann', detrend='constant', scaling='spectrum', average='mean')
¶
Calculate cross-spectral density matrix.
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 |
|---|---|
SpectralFrame
|
SpectralFrame containing cross-spectral density matrix |
Source code in wandas/frames/mixins/channel_transform_mixin.py
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 | |
transfer_function(n_fft=2048, hop_length=None, win_length=None, window='hann', detrend='constant', scaling='spectrum', average='mean')
¶
Calculate transfer function matrix.
The transfer function represents the signal transfer characteristics between channels in the frequency domain and represents the input-output relationship of the system.
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 |
|---|---|
SpectralFrame
|
SpectralFrame containing transfer function matrix |
Source code in wandas/frames/mixins/channel_transform_mixin.py
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 | |