void Formant_scatterPlot (Formant me, Graphics g, double tmin, double tmax,
int iformant1, double fmin1, double fmax1, int iformant2, double fmin2, double fmax2,
double size_mm, const char32 *mark, int garnish)
{
if (iformant1 < 1 || iformant2 < 1) return;
if (tmax <= tmin) { tmin = my xmin; tmax = my xmax; }
long itmin, itmax;
if (! Sampled_getWindowSamples (me, tmin, tmax, & itmin, & itmax)) return;
if (fmax1 == fmin1)
Formant_getExtrema (me, iformant1, tmin, tmax, & fmin1, & fmax1);
if (fmax1 == fmin1) return;
if (fmax2 == fmin2)
Formant_getExtrema (me, iformant2, tmin, tmax, & fmin2, & fmax2);
if (fmax2 == fmin2) return;
Graphics_setInner (g);
Graphics_setWindow (g, fmin1, fmax1, fmin2, fmax2);
for (long iframe = itmin; iframe <= itmax; iframe ++) {
Formant_Frame frame = & my d_frames [iframe];
if (iformant1 > frame -> nFormants || iformant2 > frame -> nFormants) continue;
double x = frame -> formant [iformant1]. frequency;
double y = frame -> formant [iformant2]. frequency;
if (x == 0.0 || y == 0.0) continue;
Graphics_mark (g, x, y, size_mm, mark);
}
Graphics_unsetInner (g);
if (garnish) {
Graphics_drawInnerBox (g);
Graphics_textBottom (g, 1, Melder_cat (U"%%F_", iformant1, U" (Hz)"));
Graphics_textLeft (g, 1, Melder_cat (U"%%F_", iformant2, U" (Hz)"));
Graphics_marksBottom (g, 2, 1, 1, 0);
Graphics_marksLeft (g, 2, 1, 1, 0);
}
}
void Polygon_drawMarks (Polygon me, Graphics g, double xmin, double xmax,
double ymin, double ymax, double size_mm, const wchar_t *mark) {
Graphics_setInner (g);
setWindow (me, g, xmin, xmax, ymin, ymax);
for (long i = 1; i <= my numberOfPoints; i++) {
Graphics_mark (g, my x[i], my y[i], size_mm, mark);
}
Graphics_unsetInner (g);
}
void Matrix_scatterPlot (Matrix me, Graphics g, long icx, long icy, double xmin, double xmax, double ymin, double ymax, double size_mm, const char32 *mark, int garnish) {
long ix = labs (icx), iy = labs (icy);
if (ix < 1 || ix > my nx || iy < 1 || iy > my nx) {
return;
}
if (xmax <= xmin) {
(void) Matrix_getWindowExtrema (me, ix, ix, 1, my ny, & xmin, & xmax);
if (xmax <= xmin) {
xmin -= 0.5; xmax += 0.5;
}
}
if (ymax <= ymin) {
(void) Matrix_getWindowExtrema (me, iy, iy, 1, my ny, & ymin, & ymax);
if (ymax <= ymin) {
ymin -= 0.5; ymax += 0.5;
}
}
Graphics_setInner (g);
if (icx < 0) {
double t = xmin;
xmin = xmax;
xmax = t;
}
if (icy < 0) {
double t = ymin;
ymin = ymax;
ymax = t;
}
Graphics_setWindow (g, xmin, xmax, ymin, ymax);
for (long i = 1; i <= my ny; i++) {
if (my z[i][ix] >= xmin && my z[i][ix] <= xmax && my z[i][iy] >= ymin && my z[i][iy] <= ymax) {
Graphics_mark (g, my z[i][ix], my z[i][iy], size_mm, mark);
}
}
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);
}
}
}