/**
* CircleTool
* Copyright 2007 by Michael Christen
* First released 22.05.2007 at http://yacy.net
*
* This library is free software; you can redistribute it and/or
* modify it under the terms of the GNU Lesser General Public
* License as published by the Free Software Foundation; either
* version 2.1 of the License, or (at your option) any later version.
*
* This library is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* Lesser General Public License for more details.
*
* You should have received a copy of the GNU Lesser General Public License
* along with this program in the file lgpl21.txt
* If not, see <http://www.gnu.org/licenses/>.
*/
package net.yacy.visualization;
import java.util.ArrayList;
import java.util.HashSet;
import java.util.List;
import java.util.Set;
public class CircleTool {
private static List<int[]> circles = new ArrayList<>();
public static void clearcache() {
circles.clear();
}
private static int[] getCircleCoords(final short radius) {
if (radius - 1 < circles.size()) return circles.get(radius - 1);
// read some lines from known circles
Set<Integer> crds = new HashSet<>();
Integer co;
for (short i = (short) Math.max(0, circles.size() - 1); i < circles.size(); i++) {
int[] circle = circles.get(i);
for (int c: circle) crds.add(c);
}
// compute more lines in new circles
short x, y;
List<Integer> crc;
short r1;
int rc = radius < 200 ? (radius < 100 ? 100 : radius + 9) : radius;
for (short r = (short) circles.size(); r < rc; r++) {
r1 = (short) (r + 1);
crc = new ArrayList<>();
for (short a = 0; a < 2 * r1; a++) {
double h = RasterPlotter.PI4 * a / r1;
double cosh = Math.cos(h);
double sinh = Math.sin(h);
x = (short) (r1 * cosh);
y = (short) (r1 * sinh);
co = x << 16 | y;
if (!(crds.contains(co))) {
crc.add(co);
crds.add(co);
}
x = (short) ((r + 0.5) * cosh);
y = (short) ((r + 0.5) * sinh);
co = x << 16 | y;
if (!(crds.contains(co))) {
crc.add(co);
crds.add(co);
}
}
// put coordinates into array
//System.out.print("Radius " + r + " => " + crc.size() + " points: ");
int[] newCircle = new int[crc.size() - 1];
int coords;
for (short i = 0; i < crc.size() - 1; i++) {
coords = crc.get(i);
newCircle[i] = coords;
}
circles.add(newCircle);
}
crc = null;
crds = null;
// finally return wanted slice
return circles.get(radius - 1);
}
public static void circle(final RasterPlotter matrix, final int xc, final int yc, final int radius, final int intensity) {
if (radius == 0) {
//matrix.plot(xc, yc, 100);
} else {
final int[] c = getCircleCoords((short) radius);
short x, y;
short limit = (short) c.length;
int co;
for (short i = 0; i < limit; i++) {
co = c[i];
x = (short) (0xffff & (co >> 16));
y = (short) (0xffff & co);
matrix.plot(xc + x , yc - y - 1, intensity); // quadrant 1
matrix.plot(xc - x + 1, yc - y - 1, intensity); // quadrant 2
matrix.plot(xc + x , yc + y , intensity); // quadrant 4
matrix.plot(xc - x + 1, yc + y , intensity); // quadrant 3
}
}
}
public static void circle(final RasterPlotter matrix, final int xc, final int yc, final int radius, int fromArc, int toArc) {
// draws only a part of a circle
// arc is given in degree
while (fromArc > 360) fromArc -=360;
while (fromArc < 0 ) fromArc +=360;
while ( toArc > 360) toArc -=360;
while ( toArc < 0 ) toArc +=360;
if (radius == 0) {
//matrix.plot(xc, yc, 100);
} else {
int[] c = getCircleCoords((short) radius);
if (c == null) c = getCircleCoords((short) radius);
final short q = (short) c.length;
final short q2 = (short) (q * 2);
final short q3 = (short) (q * 3);
final short q4 = (short) (q * 4);
final short[] c4x = new short[q4];
final short[] c4y = new short[q4];
short a0, a1, a2, a3, b0, b1;
int co;
for (short i = 0; i < q; i++) {
b0 = i;
b1 = (short) (q - 1 - i);
co = c[b0];
a0 = (short) (0xffff & (co >> 16));
a1 = (short) (0xffff & co);
co = c[b1];
a2 = (short) (0xffff & (co >> 16));
a3 = (short) (0xffff & co);
c4x[i ] = a0 ; // quadrant 1
c4y[i ] = (short) (-a1 - 1); // quadrant 1
c4x[i + q ] = (short) ( 1 - a2); // quadrant 2
c4y[i + q ] = (short) (-a3 - 1); // quadrant 2
c4x[i + q2] = (short) ( 1 - a0); // quadrant 3
c4y[i + q2] = a1 ; // quadrant 3
c4x[i + q3] = a2 ; // quadrant 4
c4y[i + q3] = a3 ; // quadrant 4
}
if (fromArc == toArc) {
int i = q4 * fromArc / 360;
matrix.plot(xc + c4x[i], yc + c4y[i], 100);
} else if (fromArc > toArc) {
// draw two parts
int from = q4 * fromArc / 360;
int to = q4 * toArc / 360;
for (int i = from; i < q4; i++) {
matrix.plot(xc + c4x[i], yc + c4y[i], 100);
}
for (int i = 0; i < to; i++) {
matrix.plot(xc + c4x[i], yc + c4y[i], 100);
}
} else {
// can be drawn in one part
int from = q4 * fromArc / 360;
int to = q4 * toArc / 360;
for (int i = from; i < to; i++) {
matrix.plot(xc + c4x[i], yc + c4y[i], 100);
}
}
}
}
}