Many resources are needed to download a project. Please understand that we have to compensate our server costs. Thank you in advance.
Project price only 1 $
You can buy this project and download/modify it how often you want.
org.math.plot.utils.Histogram Maven / Gradle / Ivy
/*
* Created on 21 juil. 2005 by richet
*/
package org.math.plot.utils;
import static org.math.plot.utils.Array.*;
public class Histogram {
// histograms functions
public static double[][] histogram_classes(double[] values, double[] bounds) {
return mergeColumns(centers(bounds), histogram(values, bounds));
}
public static double[][] histogram_classes(double[] values, double min, double max, int n) {
double[] bounds = bounds(values, min, max, n);
return mergeColumns(centers(bounds), histogram(values, bounds));
}
public static double[][] histogram_classes(double[] values, int n) {
double[] bounds = bounds(values, n);
return mergeColumns(centers(bounds), histogram(values, bounds));
}
public static double[] histogram(double[] values, double[] bounds) {
double[] h = new double[bounds.length - 1];
for (int i = 0; i < values.length; i++) {
for (int j = 0; j < h.length; j++) {
if (((bounds[j + 1] - values[i]) * (bounds[j] - values[i]) < 0) || ((bounds[j] == values[i])))
h[j]++;
}
}
return h;
}
public static double[] histogram(double[] values, double min, double max, int n) {
double[] bounds = bounds(values, min, max, n);
return histogram(values, bounds);
}
public static double[] histogram(double[] values, int n) {
return histogram(values, n);
}
private static double[] bounds(double[] values, int n) {
double min = min(values);
double max = max(values);
return bounds(values, min, max, n);
}
private static double[] bounds(double[] values, double min, double max, int n) {
double[] bounds = new double[n + 1];
for (int i = 0; i < bounds.length; i++) {
bounds[i] = min + (max - min) * i / (double) n;
}
return bounds;
}
private static double[] centers(double[] bounds) {
double[] center = new double[bounds.length - 1];
for (int i = 0; i < center.length; i++)
center[i] = (bounds[i] + bounds[i + 1]) / 2;
return center;
}
// histograms 2D functions
public static double[][] histogram_classes_2D(double[][] values, double[] boundsX, double[] boundsY) {
return insertColumn(centers_2D(boundsX, boundsY), histogram_2D(values, boundsX, boundsY), 2);
}
public static double[][] histogram_classes_2D(double[][] values, double minX, double maxX, int nX, double minY, double maxY, int nY) {
double[] valuesX = getColumnCopy(values, 0);
double[] valuesY = getColumnCopy(values, 1);
double[] boundsX = bounds(valuesX, minX, maxX, nX);
double[] boundsY = bounds(valuesY, minY, maxY, nY);
return insertColumn(centers_2D(boundsX, boundsY), histogram_2D(values, boundsX, boundsY), 2);
}
public static double[][] histogram_classes_2D(double[][] values, int nX, int nY) {
double[] valuesX = getColumnCopy(values, 0);
double[] valuesY = getColumnCopy(values, 1);
double[] boundsX = bounds(valuesX, nX);
double[] boundsY = bounds(valuesY, nY);
return insertColumn(centers_2D(boundsX, boundsY), histogram_2D(values, boundsX, boundsY), 2);
}
public static double[] histogram_2D(double[][] values, double[] boundsX, double[] boundsY) {
double[] h = new double[(boundsX.length - 1) * (boundsY.length - 1)];
for (int n = 0; n < values.length; n++) {
for (int i = 0; i < boundsX.length - 1; i++) {
for (int j = 0; j < boundsY.length - 1; j++) {
if ((((boundsX[i + 1] - values[n][0]) * (boundsX[i] - values[n][0]) < 0) || ((boundsX[i] == values[n][0])))
&& (((boundsY[j + 1] - values[n][1]) * (boundsY[j] - values[n][1]) < 0) || ((boundsY[j] == values[n][1]))))
h[index2(i, j, boundsX.length - 1)]++;
}
}
}
return h;
}
public static double[] histogram_2D(double[][] values, double minX, double maxX, int nX, double minY, double maxY, int nY) {
double[] valuesX = getColumnCopy(values, 0);
double[] valuesY = getColumnCopy(values, 1);
double[] boundsX = bounds(valuesX, minX, maxX, nX);
double[] boundsY = bounds(valuesY, minY, maxY, nY);
return histogram_2D(values, boundsX, boundsY);
}
public static double[] histogram_2D(double[][] values, int nX, int nY) {
double[] valuesX = getColumnCopy(values, 0);
double[] valuesY = getColumnCopy(values, 1);
double[] boundsX = bounds(valuesX, nX);
double[] boundsY = bounds(valuesY, nY);
return histogram_2D(values, boundsX, boundsY);
}
private static double[][] centers_2D(double[] boundsX, double[] boundsY) {
int nb_centers = (boundsX.length - 1) * (boundsY.length - 1);
double[][] center = new double[nb_centers][2];
for (int i = 0; i < boundsX.length - 1; i++) {
for (int j = 0; j < boundsY.length - 1; j++) {
int k = index2(i, j, boundsX.length - 1);
center[k][0] = (boundsX[i] + boundsX[i + 1]) / 2;
center[k][1] = (boundsY[j] + boundsY[j + 1]) / 2;
}
}
return center;
}
private static int index2(int i, int j, int imax) {
return i + imax * j;
}
// histograms 3D functions
public static double[][] histogram_classes_3D(double[][] values, int nX, int nY, int nZ) {
double[] valuesX = getColumnCopy(values, 0);
double[] valuesY = getColumnCopy(values, 1);
double[] valuesZ = getColumnCopy(values, 2);
double[] boundsX = bounds(valuesX, nX);
double[] boundsY = bounds(valuesY, nY);
double[] boundsZ = bounds(valuesZ, nZ);
return insertColumn(centers_3D(boundsX, boundsY, boundsZ), histogram_3D(values, boundsX, boundsY, boundsZ), 3);
}
public static double[] histogram_3D(double[][] values, double[] boundsX, double[] boundsY, double[] boundsZ) {
double[] h = new double[(boundsX.length - 1) * (boundsY.length - 1)* (boundsZ.length - 1)];
for (int n = 0; n < values.length; n++) {
for (int i = 0; i < boundsX.length - 1; i++) {
for (int j = 0; j < boundsY.length - 1; j++) {
for (int k = 0; k < boundsZ.length - 1; k++) {
if ((((boundsX[i + 1] - values[n][0]) * (boundsX[i] - values[n][0]) < 0) || ((boundsX[i] == values[n][0])))
&& (((boundsY[j + 1] - values[n][1]) * (boundsY[j] - values[n][1]) < 0) || ((boundsY[j] == values[n][1])))&& (((boundsZ[k + 1] - values[n][2]) * (boundsZ[k] - values[n][2]) < 0) || ((boundsZ[k] == values[n][2]))))
h[index3(i, j,k, boundsX.length - 1, boundsY.length - 1)]++;
}}
}
}
return h;
}
private static double[][] centers_3D(double[] boundsX, double[] boundsY, double[] boundsZ) {
int nb_centers = (boundsX.length - 1) * (boundsY.length - 1)* (boundsZ.length - 1);
double[][] center = new double[nb_centers][3];
for (int i = 0; i < boundsX.length - 1; i++) {
for (int j = 0; j < boundsY.length - 1; j++) {
for (int k = 0; k < boundsZ.length - 1; k++) {
int l = index3(i, j,k, boundsX.length - 1, boundsY.length - 1);
center[l][0] = (boundsX[i] + boundsX[i + 1]) / 2;
center[l][1] = (boundsY[j] + boundsY[j + 1]) / 2;
center[l][2] = (boundsZ[k] + boundsZ[k + 1]) / 2;
}
}
}
return center;
}
private static int index3(int i, int j, int k,int imax, int jmax) {
return i + imax * j + imax*jmax*k;
}
}
