public SpectrumPlot(double[] periods) {
Application app = Application.getInstance();
if (app.isUseFrequencies()) {
xValues = new double[periods.length];
for (int i = 0; i < periods.length; ++i) {
xValues[xValues.length - i - 1] = 1 / periods[i];
}
} else {
xValues = periods;
}
this.reference1 = readReference(Application.PropertySpecRef1);
this.reference2 = readReference(Application.PropertySpecRef2);
this.logScale = Application.getInstance().getProperty(Application.PropertySpecLogScale, false);
// this.setBackground(BackgroundColor);
}
@Override
protected void paintComponent(Graphics g) {
super.paintComponent(g);
Graphics2D g2 = (Graphics2D) g;
g2.setRenderingHint(RenderingHints.KEY_ANTIALIASING, RenderingHints.VALUE_ANTIALIAS_ON);
int width = getWidth() - Border.left - Border.right;
int height = getHeight() - Border.top - Border.bottom;
g2.translate(Border.left, Border.top);
// create x and y axes
g2.setColor(Color.black);
g2.drawLine(0, height, width, height);
g2.drawLine(0, 0, 0, height);
g2.setFont(g2.getFont().deriveFont(10.0f));
FontMetrics fm = g2.getFontMetrics();
if (xValues.length == 0 || target == null || xValues[xValues.length - 1] <= 0) {
return;
}
double factor;
String yText;
Application app = Application.getInstance();
Shaking.Type t = app.getSpectrumParameter();
if (t == Shaking.Type.PSA) {
factor = Application.EarthAcceleration1;
yText = t.toString().toUpperCase() + ", g";
} else if (t == Shaking.Type.DRS) {
factor = 100;
yText = t.toString().toUpperCase() + ", cm";
} else {
LOG.warn("unsupported spectrum parameter: " + t.toString());
return;
}
synchronized (target) {
// determine min/max of x and y axis
double xMin;
double xMax;
double dx = 0;
if (logScale) {
xMin = Math.log10(xValues[0]);
xMax = Math.log10(xValues[xValues.length - 1]);
} else {
xMin = 0.0;
xMax = xValues[xValues.length - 1];
}
if (xMax > xMin) {
dx = (double) width / (xMax - xMin);
}
double yMin = 0;
double yMax = 0;
double dy = 0;
boolean first = true;
for (double v : reference1) {
if (first) {
first = false;
/*yMin =*/ yMax = v;
}
/*yMin = Math.min(yMin, v);*/
yMax = Math.max(yMax, v);
}
for (double v : reference2) {
if (first) {
first = false;
/*yMin =*/ yMax = v;
}
/*yMin = Math.min(yMin, v);*/
yMax = Math.max(yMax, v);
}
for (Shaking s : target.spectralValues) {
if (first) {
first = false;
/*yMin =*/ yMax = s.expectedSI * factor;
} else {
/*yMin = Math.min(yMin, s.expectedSI * factor);*/
yMax = Math.max(yMax, s.expectedSI * factor);
}
/*
yMin = Math.min(yMin, s.expectedSI * factor);
yMin = Math.min(yMin, s.percentile84 * factor);
yMin = Math.min(yMin, s.percentile16 * factor);
*/
yMax = Math.max(yMax, s.expectedSI * factor);
yMax = Math.max(yMax, s.percentile84 * factor);
yMax = Math.max(yMax, s.percentile16 * factor);
}
if (yMax > yMin) {
dy = (double) height / (yMax - yMin);
}
List<Point> expectedPoints = new ArrayList();
List<Point> percentile84Points = new ArrayList();
List<Point> percentile16Points = new ArrayList();
List<Point> ref1Points = new ArrayList();
List<Point> ref2Points = new ArrayList();
int x, y, iV;
boolean isFreq = app.isUseFrequencies();
for (int i = 0; i < xValues.length; i++) {
if (logScale) {
x = (int) ((Math.log10(xValues[i]) - xMin) / (xMax - xMin) * width);
} else {
x = (int) ((xValues[i] - xMin) / (xMax - xMin) * width);
}
iV = isFreq ? xValues.length - i - 1 : i;
if (iV < target.spectralValues.size()) {
Shaking s = target.spectralValues.get(iV);
y = (int) ((yMax - s.expectedSI * factor) * dy);
expectedPoints.add(new Point(x, y));
y = (int) ((yMax - s.percentile84 * factor) * dy);
percentile84Points.add(new Point(x, y));
y = (int) ((yMax - s.percentile16 * factor) * dy);
percentile16Points.add(new Point(x, y));
}
if (iV < reference1.size()) {
y = (int) ((yMax - reference1.get(iV)) * dy);
ref1Points.add(new Point(x, y));
}
if (iV < reference2.size()) {
y = (int) ((yMax - reference2.get(iV)) * dy);
ref2Points.add(new Point(x, y));
}
}
// Y-axis
String text;
int halfAscent = fm.getAscent() / 2 - 1;
double q = Math.log10(2 * (yMax - yMin) * 2 * fm.getHeight() / height);
double rx = q - Math.floor(q);
int d = rx < 0.3 ? 1 : rx > 0.7 ? 5 : 2;
double tickStep = d * Math.pow(10, Math.floor(q - rx));
for (double v = 0; v < yMax; v += tickStep) {
y = height - (int) (v * dy) - 1;
g2.drawLine(0, y, -TickLength, y);
text = Double.toString(((int) (v * 1000.0 + 0.5)) / 1000.0);
int w = (int) fm.getStringBounds(text, null).getWidth();
g2.drawString(text, -TickLength - 3 - w, y + halfAscent);
}
// X-axis
q = Math.log10(2 * (xMax - xMin) * 50 / width);
rx = q - Math.floor(q);
d = rx < 0.3 ? 1 : rx > 0.7 ? 5 : 2;
tickStep = d * Math.pow(10, Math.floor(q - rx));
for (double v = xMin; v < xMax; v += tickStep) {
x = (int) ((v - xMin) * dx);
g2.drawLine(x, height, x, height + TickLength);
text = String.format("%.2f", logScale ? Math.pow(10, v) : v);
// String text = Double.toString(logScale ? Math.pow(10, v) : v);//((int) (v * 1000.0 +
// 0.5)) / 1000.0);
int w = (int) fm.getStringBounds(text, null).getWidth();
g2.drawString(text, x - w / 2, height + TickLength + 3 + fm.getAscent());
}
// axis label
g2.setFont(g2.getFont().deriveFont(14.0f));
AffineTransform orig = g2.getTransform();
int w = (int) ((height + fm.getStringBounds(yText, null).getWidth()) / 2);
g2.translate(-Border.left + fm.getAscent() + 7, w);
g2.rotate(-Math.PI / 2);
g2.drawString(yText, 0, 0);
g2.setTransform(orig);
text = "vibration " + (isFreq ? "frequency, Hz" : "period, s");
g2.drawString(
text,
(int) ((width - fm.getStringBounds(text, null).getWidth()) / 2),
height + TickLength + 3 + 2 * fm.getHeight());
g2.setClip(0, 0, width, height);
g2.setStroke(DefaultStroke);
g2.setColor(Ref1Color);
drawGraph(g2, ref1Points);
g2.setColor(Ref2Color);
drawGraph(g2, ref2Points);
g2.setStroke(PercentileStroke);
g2.setColor(ShakingColor);
drawGraph(g2, percentile16Points);
drawGraph(g2, percentile84Points);
g2.setStroke(DefaultStroke);
drawGraph(g2, expectedPoints);
}
}
