void Matrix_drawSliceY (Matrix me, Graphics g, double x, double ymin, double ymax, double min, double max) {
if (x < my xmin || x > my xmax) {
return;
}
long ix = Matrix_xToNearestColumn (me, x);
if (ymax <= ymin) {
ymin = my ymin;
ymax = my ymax;
}
long iymin, iymax;
long ny = Matrix_getWindowSamplesY (me, ymin, ymax, &iymin, &iymax);
if (ny < 1) {
return;
}
if (max <= min) {
Matrix_getWindowExtrema (me, ix, ix, iymin, iymax, &min, &max);
}
if (max <= min) {
min -= 0.5; max += 0.5;
}
autoNUMvector<double> y (iymin, iymax);
Graphics_setWindow (g, ymin, ymax, min, max);
Graphics_setInner (g);
for (long i = iymin; i <= iymax; i++) {
y[i] = my z[i][ix];
}
Graphics_function (g, y.peek(), iymin, iymax, Matrix_rowToY (me, iymin), Matrix_rowToY (me, iymax));
Graphics_unsetInner (g);
}
void PowerCepstrogram_paint (PowerCepstrogram me, Graphics g, double tmin, double tmax, double qmin, double qmax, double dBmaximum, int autoscaling, double dynamicRangedB, double dynamicCompression, int garnish) {
if (tmax <= tmin) { tmin = my xmin; tmax = my xmax; }
if (qmax <= qmin) { qmin = my ymin; qmax = my ymax; }
long itmin, itmax, ifmin, ifmax;
if (! Matrix_getWindowSamplesX (me, tmin - 0.49999 * my dx, tmax + 0.49999 * my dx, & itmin, & itmax) ||
! Matrix_getWindowSamplesY (me, qmin - 0.49999 * my dy, qmax + 0.49999 * my dy, & ifmin, & ifmax)) {
return;
}
autoMatrix thee = Data_copy (me);
double min = 1e308, max = -min;
for (long i = 1; i <= my ny; i++) {
for (long j = 1; j <= my nx; j++) {
double val = TO10LOG (my z[i][j]);
min = val < min ? val : min;
max = val > max ? val : max;
thy z[i][j] = val;
}
}
double dBminimum = dBmaximum - dynamicRangedB;
if (autoscaling) {
dBminimum = min; dBmaximum = max;
}
for (long j = 1; j <= my nx; j++) {
double lmax = thy z[1][j];
for (long i = 2; i <= my ny; i++) {
if (thy z[i][j] > lmax) {
lmax = thy z[i][j];
}
}
double factor = dynamicCompression * (max - lmax);
for (long i = 1; i <= my ny; i++) {
thy z[i][j] += factor;
}
}
Graphics_setInner (g);
Graphics_setWindow (g, tmin, tmax, qmin, qmax);
Graphics_image (g, thy z,
itmin, itmax,
Matrix_columnToX (thee.get(), itmin - 0.5),
Matrix_columnToX (thee.get(), itmax + 0.5),
ifmin, ifmax,
Matrix_rowToY (thee.get(), ifmin - 0.5),
Matrix_rowToY (thee.get(), ifmax + 0.5),
dBminimum, dBmaximum);
Graphics_unsetInner (g);
if (garnish) {
Graphics_drawInnerBox (g);
Graphics_textBottom (g, true, U"Time (s)");
Graphics_marksBottom (g, 2, true, true, false);
Graphics_marksLeft (g, 2, true, true, false);
Graphics_textLeft (g, true, U"Quefrency (s)");
}
}
void Spectrogram_paintInside (Spectrogram me, Graphics g, double tmin, double tmax, double fmin, double fmax,
double maximum, int autoscaling, double dynamic, double preemphasis, double dynamicCompression)
{
if (tmax <= tmin) { tmin = my xmin; tmax = my xmax; }
if (fmax <= fmin) { fmin = my ymin; fmax = my ymax; }
long itmin, itmax, ifmin, ifmax;
if (! Matrix_getWindowSamplesX (me, tmin - 0.49999 * my dx, tmax + 0.49999 * my dx, & itmin, & itmax) ||
! Matrix_getWindowSamplesY (me, fmin - 0.49999 * my dy, fmax + 0.49999 * my dy, & ifmin, & ifmax))
return;
Graphics_setWindow (g, tmin, tmax, fmin, fmax);
autoNUMvector <double> preemphasisFactor (ifmin, ifmax);
autoNUMvector <double> dynamicFactor (itmin, itmax);
/* Pre-emphasis in place; also compute maximum after pre-emphasis. */
for (long ifreq = ifmin; ifreq <= ifmax; ifreq ++) {
preemphasisFactor [ifreq] = (preemphasis / NUMln2) * log (ifreq * my dy / 1000.0);
for (long itime = itmin; itime <= itmax; itime ++) {
double value = my z [ifreq] [itime]; /* Power. */
value = (10.0/NUMln10) * log ((value + 1e-30) / 4.0e-10) + preemphasisFactor [ifreq]; /* dB */
if (value > dynamicFactor [itime]) dynamicFactor [itime] = value; /* Local maximum. */
my z [ifreq] [itime] = value;
}
}
/* Compute global maximum. */
if (autoscaling) {
maximum = 0.0;
for (long itime = itmin; itime <= itmax; itime ++)
if (dynamicFactor [itime] > maximum) maximum = dynamicFactor [itime];
}
/* Dynamic compression in place. */
for (long itime = itmin; itime <= itmax; itime ++) {
dynamicFactor [itime] = dynamicCompression * (maximum - dynamicFactor [itime]);
for (long ifreq = ifmin; ifreq <= ifmax; ifreq ++)
my z [ifreq] [itime] += dynamicFactor [itime];
}
Graphics_image (g, my z,
itmin, itmax,
Matrix_columnToX (me, itmin - 0.5),
Matrix_columnToX (me, itmax + 0.5),
ifmin, ifmax,
Matrix_rowToY (me, ifmin - 0.5),
Matrix_rowToY (me, ifmax + 0.5),
maximum - dynamic, maximum);
for (long ifreq = ifmin; ifreq <= ifmax; ifreq ++)
for (long itime = itmin; itime <= itmax; itime ++) {
double value = 4.0e-10 * exp ((my z [ifreq] [itime] - dynamicFactor [itime]
- preemphasisFactor [ifreq]) * (NUMln10 / 10.0)) - 1e-30;
my z [ifreq] [itime] = value > 0.0 ? value : 0.0;
}
}
void Matrix_paintSurface (Matrix me, Graphics g, double xmin, double xmax, double ymin, double ymax,
double minimum, double maximum, double elevation, double azimuth)
{
if (xmax <= xmin) { xmin = my xmin; xmax = my xmax; }
if (ymax <= ymin) { ymin = my ymin; ymax = my ymax; }
long ixmin, ixmax, iymin, iymax;
(void) Matrix_getWindowSamplesX (me, xmin, xmax, & ixmin, & ixmax);
(void) Matrix_getWindowSamplesY (me, ymin, ymax, & iymin, & iymax);
if (maximum <= minimum)
(void) Matrix_getWindowExtrema (me, ixmin, ixmax, iymin, iymax, & minimum, & maximum);
if (maximum <= minimum) { minimum -= 1.0; maximum += 1.0; }
Graphics_setInner (g);
Graphics_setWindow (g, -1.0, 1.0, minimum, maximum);
Graphics_surface (g, my z,
ixmin, ixmax, Matrix_columnToX (me, ixmin), Matrix_columnToX (me, ixmax),
iymin, iymax, Matrix_rowToY (me, iymin), Matrix_rowToY (me, iymax),
minimum, maximum, elevation, azimuth);
Graphics_unsetInner (g);
}
void BarkFilter_drawSekeyHansonFilterFunctions (BarkFilter me, Graphics g,
int toFreqScale, int fromFilter, int toFilter, double zmin, double zmax,
int dbScale, double ymin, double ymax, int garnish) {
if (! checkLimits (me, FilterBank_BARK, toFreqScale, & fromFilter, & toFilter,
& zmin, & zmax, dbScale, & ymin, & ymax)) {
return;
}
long n = 1000;
autoNUMvector<double> a (1, n);
Graphics_setInner (g);
Graphics_setWindow (g, zmin, zmax, ymin, ymax);
for (long j = fromFilter; j <= toFilter; j++) {
double df = (zmax - zmin) / (n - 1);
double zMid = Matrix_rowToY (me, j);
long ibegin, iend;
for (long i = 1; i <= n; i++) {
double f = zmin + (i - 1) * df;
double z = scaleFrequency (f, toFreqScale, FilterBank_BARK);
if (z == NUMundefined) {
a[i] = NUMundefined;
} else {
z -= zMid + 0.215;
a[i] = 7 - 7.5 * z - 17.5 * sqrt (0.196 + z * z);
if (! dbScale) {
a[i] = pow (10, a[i]);
}
}
}
setDrawingLimits (a.peek(), n, ymin, ymax, &ibegin, &iend);
if (ibegin <= iend) {
double fmin = zmin + (ibegin - 1) * df;
double fmax = zmax - (n - iend) * df;
Graphics_function (g, a.peek(), ibegin, iend, fmin, fmax);
}
}
Graphics_unsetInner (g);
if (garnish) {
double distance = dbScale ? 10 : 1;
const char32 *ytext = dbScale ? U"Amplitude (dB)" : U"Amplitude";
Graphics_drawInnerBox (g);
Graphics_marksBottom (g, 2, 1, 1, 0);
Graphics_marksLeftEvery (g, 1, distance, 1, 1, 0);
Graphics_textLeft (g, 1, ytext);
Graphics_textBottom (g, 1, GetFreqScaleText (toFreqScale));
}
}
void PowerCepstrogram_paint (PowerCepstrogram me, Graphics g, double tmin, double tmax, double qmin, double qmax, double dBminimum, double dBmaximum, int garnish) {
if (tmax <= tmin) { tmin = my xmin; tmax = my xmax; }
if (qmax <= qmin) { qmin = my ymin; qmax = my ymax; }
long itmin, itmax, ifmin, ifmax;
if (! Matrix_getWindowSamplesX (me, tmin - 0.49999 * my dx, tmax + 0.49999 * my dx, & itmin, & itmax) ||
! Matrix_getWindowSamplesY (me, qmin - 0.49999 * my dy, qmax + 0.49999 * my dy, & ifmin, & ifmax)) {
return;
}
autoMatrix thee = (Matrix) Data_copy (me);
double min = 1e38, max = -min;
for (long i = 1; i <= my ny; i++) {
for (long j = 1; j <= my nx; j++) {
double val = TO10LOG (my z[i][j]);
min = val < min ? val : min;
max = val > max ? val : max;
thy z[i][j] = val;
}
}
if (dBmaximum <= dBminimum) {
dBminimum = min; dBmaximum = max;
}
Graphics_setInner (g);
Graphics_setWindow (g, tmin, tmax, qmin, qmax);
Graphics_image (g, thy z,
itmin, itmax,
Matrix_columnToX (thee.peek(), itmin - 0.5),
Matrix_columnToX (thee.peek(), itmax + 0.5),
ifmin, ifmax,
Matrix_rowToY (thee.peek(), ifmin - 0.5),
Matrix_rowToY (thee.peek(), ifmax + 0.5),
dBminimum, dBmaximum);
Graphics_unsetInner (g);
if (garnish) {
Graphics_drawInnerBox (g);
Graphics_textBottom (g, 1, L"Time (s)");
Graphics_marksBottom (g, 2, 1, 1, 0);
Graphics_marksLeft (g, 2, 1, 1, 0);
Graphics_textLeft (g, 1, L"Quefrency (s)");
}
}
void Matrix_drawOneContour (Matrix me, Graphics g, double xmin, double xmax, double ymin, double ymax,
double height)
{
bool xreversed = xmin > xmax, yreversed = ymin > ymax;
if (xmax == xmin) { xmin = my xmin; xmax = my xmax; }
if (ymax == ymin) { ymin = my ymin; ymax = my ymax; }
if (xreversed) { double temp = xmin; xmin = xmax; xmax = temp; }
if (yreversed) { double temp = ymin; ymin = ymax; ymax = temp; }
long ixmin, ixmax, iymin, iymax;
(void) Matrix_getWindowSamplesX (me, xmin, xmax, & ixmin, & ixmax);
(void) Matrix_getWindowSamplesY (me, ymin, ymax, & iymin, & iymax);
if (xmin == xmax || ymin == ymax) return;
Graphics_setInner (g);
Graphics_setWindow (g, xreversed ? xmax : xmin, xreversed ? xmin : xmax, yreversed ? ymax : ymin, yreversed ? ymin : ymax);
Graphics_contour (g, my z,
ixmin, ixmax, Matrix_columnToX (me, ixmin), Matrix_columnToX (me, ixmax),
iymin, iymax, Matrix_rowToY (me, iymin), Matrix_rowToY (me, iymax),
height);
Graphics_rectangle (g, xmin, xmax, ymin, ymax);
Graphics_unsetInner (g);
}
void Matrix_drawAsSquares (Matrix me, Graphics g, double xmin, double xmax, double ymin, double ymax, int garnish) {
Graphics_Colour colour = Graphics_inqColour (g);
long ixmin, ixmax, iymin, iymax;
if (xmax <= xmin) {
xmin = my xmin;
xmax = my xmax;
}
long nx = Matrix_getWindowSamplesX (me, xmin, xmax, &ixmin, &ixmax);
if (ymax <= ymin) {
ymin = my ymin;
ymax = my ymax;
}
long ny = Matrix_getWindowSamplesY (me, ymin, ymax, &iymin, &iymax);
double min, max = nx > ny ? nx : ny;
double dx = (xmax - xmin) / max, dy = (ymax - ymin) / max;
Graphics_setInner (g);
Graphics_setWindow (g, xmin, xmax, ymin, ymax);
Matrix_getWindowExtrema (me, ixmin, ixmax, iymin, iymax, & min, & max);
double wAbsMax = fabs (max) > fabs (min) ? fabs (max) : fabs (min);
for (long i = iymin; i <= iymax; i++) {
double y = Matrix_rowToY (me, i);
for (long j = ixmin; j <= ixmax; j++) {
double x = Matrix_columnToX (me, j);
double d = 0.95 * sqrt (fabs (my z[i][j]) / wAbsMax);
if (d > 0) {
double x1WC = x - d * dx / 2, x2WC = x + d * dx / 2;
double y1WC = y - d * dy / 2, y2WC = y + d * dy / 2;
if (my z[i][j] > 0) {
Graphics_setColour (g, Graphics_WHITE);
}
Graphics_fillRectangle (g, x1WC, x2WC, y1WC, y2WC);
Graphics_setColour (g, colour);
Graphics_rectangle (g, x1WC, x2WC , y1WC, y2WC);
}
}
}
Graphics_setGrey (g, 0.0);
Graphics_unsetInner (g);
if (garnish) {
Graphics_drawInnerBox (g);
Graphics_marksLeft (g, 2, true, true, false);
if (ymin * ymax < 0.0) {
Graphics_markLeft (g, 0.0, true, true, true, nullptr);
}
Graphics_marksBottom (g, 2, true, true, false);
if (xmin * xmax < 0.0) {
Graphics_markBottom (g, 0.0, true, true, true, nullptr);
}
}
}
void Matrix_paintContours (Matrix me, Graphics g, double xmin, double xmax, double ymin, double ymax,
double minimum, double maximum)
{
double border [1 + 30];
if (xmax <= xmin) { xmin = my xmin; xmax = my xmax; }
if (ymax <= ymin) { ymin = my ymin; ymax = my ymax; }
long ixmin, ixmax, iymin, iymax, iborder;
(void) Matrix_getWindowSamplesX (me, xmin, xmax, & ixmin, & ixmax);
(void) Matrix_getWindowSamplesY (me, ymin, ymax, & iymin, & iymax);
if (maximum <= minimum)
(void) Matrix_getWindowExtrema (me, ixmin, ixmax, iymin, iymax, & minimum, & maximum);
if (maximum <= minimum) { minimum -= 1.0; maximum += 1.0; }
for (iborder = 1; iborder <= 30; iborder ++)
border [iborder] = minimum + iborder * (maximum - minimum) / (30 + 1);
if (xmin >= xmax || ymin >= ymax) return;
Graphics_setInner (g);
Graphics_setWindow (g, xmin, xmax, ymin, ymax);
Graphics_grey (g, my z,
ixmin, ixmax, Matrix_columnToX (me, ixmin), Matrix_columnToX (me, ixmax),
iymin, iymax, Matrix_rowToY (me, iymin), Matrix_rowToY (me, iymax),
30, border);
Graphics_rectangle (g, xmin, xmax, ymin, ymax);
Graphics_unsetInner (g);
}
autoSound Spectrogram_to_Sound (Spectrogram me, double fsamp) {
try {
double dt = 1 / fsamp;
long n = (long) floor ((my xmax - my xmin) / dt);
if (n < 0) return autoSound ();
autoSound thee = Sound_create (1, my xmin, my xmax, n, dt, 0.5 * dt);
for (long i = 1; i <= n; i ++) {
double t = Sampled_indexToX (thee.peek(), i);
double rframe = Sampled_xToIndex (me, t), phase, value = 0.0;
long leftFrame, rightFrame;
if (rframe < 1 || rframe >= my nx) continue;
leftFrame = (long) floor (rframe), rightFrame = leftFrame + 1, phase = rframe - leftFrame;
for (long j = 1; j <= my ny; j ++) {
double f = Matrix_rowToY (me, j);
double power = my z [j] [leftFrame] * (1 - phase) + my z [j] [rightFrame] * phase;
value += sqrt (power) * sin (2 * NUMpi * f * t);
}
thy z [1] [i] = value;
}
return thee;
} catch (MelderError) {
Melder_throw (me, U": not converted to Sound.");
}
}
void MelSpectrogram_drawTriangularFilterFunctions (MelSpectrogram me, Graphics g, bool xIsHertz, int fromFilter, int toFilter, double zmin, double zmax, bool yscale_dB, double ymin, double ymax, int garnish) {
double xmin = zmin, xmax = zmax;
if (zmin >= zmax) {
zmin = my ymin; zmax = my ymax; // mel
xmin = xIsHertz ? my v_frequencyToHertz (zmin) : zmin;
xmax = xIsHertz ? my v_frequencyToHertz (zmax) : zmax;
}
if (xIsHertz) {
zmin = my v_hertzToFrequency (xmin); zmax = my v_hertzToFrequency (xmax);
}
if (ymin >= ymax) {
ymin = yscale_dB ? -60 : 0;
ymax = yscale_dB ? 0 : 1;
}
fromFilter = fromFilter <= 0 ? 1 : fromFilter;
toFilter = toFilter <= 0 || toFilter > my ny ? my ny : toFilter;
if (fromFilter > toFilter) {
fromFilter = 1; toFilter = my ny;
}
long n = xIsHertz ? 1000 : 500;
autoNUMvector<double> xz (1, n), xhz (1,n), y (1, n);
Graphics_setInner (g);
Graphics_setWindow (g, xmin, xmax, ymin, ymax);
double dz = (zmax - zmin) / (n - 1);
for (long iz = 1; iz <= n; iz++) {
double f = zmin + (iz - 1) * dz;
xz[iz] = f;
xhz[iz] = my v_frequencyToHertz (f); // just in case we need the linear scale
}
for (long ifilter = fromFilter; ifilter <= toFilter; ifilter++) {
double zc = Matrix_rowToY (me, ifilter), zl = zc - my dy, zh = zc + my dy;
double xo1, yo1, xo2, yo2;
if (yscale_dB) {
for (long iz = 1; iz <= n; iz++) {
double z = xz[iz];
double amp = NUMtriangularfilter_amplitude (zl, zc, zh, z);
y[iz] = yscale_dB ? (amp > 0 ? 20 * log10 (amp) : ymin - 10) : amp;
}
double x1 = xIsHertz ? xhz[1] : xz[1], y1 = y[1];
if (NUMdefined (y1)) {
for (long iz = 1; iz <= n; iz++) {
double x2 = xIsHertz ? xhz[iz] : xz[iz], y2 = y[iz];
if (NUMdefined (y2)) {
if (NUMclipLineWithinRectangle (x1, y1, x2, y2, xmin, ymin, xmax, ymax, &xo1, &yo1, &xo2, &yo2)) {
Graphics_line (g, xo1, yo1, xo2, yo2);
}
}
x1 = x2; y1 = y2;
}
}
} else {
double x1 = xIsHertz ? my v_frequencyToHertz (zl) : zl;
double x2 = xIsHertz ? my v_frequencyToHertz (zc) : zc;
if (NUMclipLineWithinRectangle (x1, 0, x2, 1, xmin, ymin, xmax, ymax, &xo1, &yo1, &xo2, &yo2)) {
Graphics_line (g, xo1, yo1, xo2, yo2);
}
double x3 = xIsHertz ? my v_frequencyToHertz (zh) : zh;
if (NUMclipLineWithinRectangle (x2, 1, x3, 0, xmin, ymin, xmax, ymax, &xo1, &yo1, &xo2, &yo2)) {
Graphics_line (g, xo1, yo1, xo2, yo2);
}
}
}
Graphics_unsetInner (g);
if (garnish) {
Graphics_drawInnerBox (g);
Graphics_marksBottom (g, 2, 1, 1, 0);
Graphics_marksLeftEvery (g, 1, yscale_dB ? 10 : 0.5, 1, 1, 0);
Graphics_textLeft (g, 1, yscale_dB ? U"Amplitude (dB)" : U"Amplitude");
Graphics_textBottom (g, 1, Melder_cat (U"Frequency (", ( xIsHertz ? U"Hz" : my v_getFrequencyUnit () ), U")"));
}
}
void BarkSpectrogram_drawSekeyHansonFilterFunctions (BarkSpectrogram me, Graphics g, bool xIsHertz, int fromFilter, int toFilter, double zmin, double zmax, bool yscale_dB, double ymin, double ymax, int garnish) {
double xmin = zmin, xmax = zmax;
if (zmin >= zmax) {
zmin = my ymin; zmax = my ymax;
xmin = xIsHertz ? my v_frequencyToHertz (zmin) : zmin;
xmax = xIsHertz ? my v_frequencyToHertz (zmax) : zmax;
}
if (xIsHertz) {
zmin = my v_hertzToFrequency (xmin); zmax = my v_hertzToFrequency (xmax);
}
if (ymin >= ymax) {
ymin = yscale_dB ? -60 : 0;
ymax = yscale_dB ? 0 : 1;
}
fromFilter = fromFilter <= 0 ? 1 : fromFilter;
toFilter = toFilter <= 0 || toFilter > my ny ? my ny : toFilter;
if (fromFilter > toFilter) {
fromFilter = 1; toFilter = my ny;
}
long n = xIsHertz ? 1000 : 500;
autoNUMvector<double> xz (1, n), xhz (1,n), y (1, n);
Graphics_setInner (g);
Graphics_setWindow (g, xmin, xmax, ymin, ymax);
double dz = (zmax - zmin) / (n - 1);
for (long iz = 1; iz <= n; iz++) {
double f = zmin + (iz - 1) * dz;
xz[iz] = f;
xhz[iz] = my v_frequencyToHertz (f); // just in case we need the linear scale
}
for (long ifilter = fromFilter; ifilter <= toFilter; ifilter++) {
double zMid = Matrix_rowToY (me, ifilter);
for (long iz = 1; iz <= n; iz++) {
double z = xz[iz] - (zMid - 0.215);
double amp = 7 - 7.5 * z - 17.5 * sqrt (0.196 + z * z);
y[iz] = yscale_dB ? amp : pow (10, amp / 10);
}
// the drawing
double x1 = xIsHertz ? xhz[1] : xz[1], y1 = y[1];
for (long iz = 2; iz <= n; iz++) {
double x2 = xIsHertz ? xhz[iz] : xz[iz], y2 = y[iz];
if (NUMdefined (x1) && NUMdefined (x2)) {
double xo1, yo1, xo2, yo2;
if (NUMclipLineWithinRectangle (x1, y1, x2, y2, xmin, ymin, xmax, ymax, &xo1, &yo1, &xo2, &yo2)) {
Graphics_line (g, xo1, yo1, xo2, yo2);
}
}
x1 = x2; y1 = y2;
}
}
Graphics_unsetInner (g);
if (garnish) {
double distance = yscale_dB ? 10 : 0.5;
Graphics_drawInnerBox (g);
Graphics_marksBottom (g, 2, 1, 1, 0);
Graphics_marksLeftEvery (g, 1, distance, 1, 1, 0);
Graphics_textLeft (g, 1, yscale_dB ? U"Amplitude (dB)" : U"Amplitude");
Graphics_textBottom (g, 1, Melder_cat (U"Frequency (", xIsHertz ? U"Hz" : my v_getFrequencyUnit (), U")"));
}
}
