/* rotate counterclockwise w.r.t. (xc,yc) */
void Polygon_rotate (Polygon me, double alpha, double xc, double yc) {
double f = alpha * NUMpi / 180, cosa = cos (f), sina = sin (f);
Polygon_translate (me, -xc, -yc);
for (long i = 1; i <= my numberOfPoints; i++) {
double x = my x[i];
my x[i] = cosa * my x[i] - sina * my y[i];
my y[i] = sina * x + cosa * my y[i];
}
Polygon_translate (me, xc, yc);
}
void Confusion_Matrix_draw (Confusion me, Matrix thee, Graphics g, long index, double lowerPercentage, double xmin, double xmax, double ymin, double ymax, int garnish) {
long ib = 1, ie = my numberOfRows;
if (index > 0 && index <= my numberOfColumns) {
ib = ie = index;
}
if (thy ny != my numberOfRows) {
Melder_throw (U"Wrong number of positions.");
}
if (xmax <= xmin) {
(void) Matrix_getWindowExtrema (thee, 1, 1, 1, thy ny, &xmin, &xmax);
}
if (xmax <= xmin) {
return;
}
if (ymax <= ymin) {
(void) Matrix_getWindowExtrema (thee, 2, 2, 1, thy ny, &ymin, &ymax);
}
if (ymax <= ymin) {
return;
}
double rmax = fabs (xmax - xmin) / 10.0;
double rmin = rmax / 10;
Graphics_setInner (g);
Graphics_setWindow (g, xmin - rmax, xmax + rmax, ymin - rmax, ymax + rmax);
Graphics_setTextAlignment (g, Graphics_CENTRE, Graphics_HALF);
for (long i = 1; i <= my numberOfRows; i++) {
Graphics_text (g, thy z[i][1], thy z[i][2], my rowLabels[i]);
}
for (long i = ib; i <= ie; i++) {
double xSum = 0.0;
for (long j = 1; j <= my numberOfColumns; j++) {
xSum += my data[i][j];
}
if (xSum <= 0.0) {
continue; /* no confusions */
}
double x1 = thy z[i][1];
double y1 = thy z[i][2];
double r = rmax * my data[i][i] / xSum;
Graphics_circle (g, x1, y1, r > rmin ? r : rmin);
for (long j = 1; j <= my numberOfColumns; j++) {
double x2 = thy z[j][1], y2 = thy z[j][2];
double perc = 100.0 * my data[i][j] / xSum;
double dx = x2 - x1, dy = y2 - y1;
double alpha = atan2 (dy, dx);
if (perc == 0.0 || perc < lowerPercentage || j == i) {
continue;
}
xmin = x1; xmax = x2;
if (x2 < x1) {
xmin = x2; xmax = x1;
}
ymin = y1; xmax = y2;
if (y2 < y1) {
ymin = y2; ymax = y1;
}
autoPolygon p = Polygon_createPointer();
double xs = sqrt (dx * dx + dy * dy) - 2.2 * r;
if (xs < 0.0) {
xs = 0.0;
}
double ys = perc * rmax / 100.0;
Polygon_scale (p.get(), xs, ys);
Polygon_translate (p.get(), x1, y1 - ys / 2);
Polygon_rotate (p.get(), alpha, x1, y1);
Polygon_translate (p.get(), 1.1 * r * cos (alpha), 1.1 * r * sin (alpha));
Polygon_drawInside (p.get(), g);
}
}
Graphics_unsetInner (g);
if (garnish) {
Graphics_drawInnerBox (g);
Graphics_marksBottom (g, 2, true, true, false);
if (ymin * ymax < 0.0) {
Graphics_markLeft (g, 0.0, true, true, true, nullptr);
}
Graphics_marksLeft (g, 2, true, true, false);
if (xmin * xmax < 0.0) {
Graphics_markBottom (g, 0.0, true, true, true, nullptr);
}
}
}
