private Vector3f interpolate_vertex(final Vector3f vertex0, final Vector3f vertex1) {
float value0 = this.substitute_plane_equation(vertex0);
float value1 = this.substitute_plane_equation(vertex1);
float ratio = kvs.core.util.Math.abs(value0 / (value1 - value0));
return (vertex0.mul(1.0f - ratio).add(vertex1.mul(ratio)));
}
public final boolean isIntersected(final Vector3f v0, final Vector3f v1, final Vector3f v2) {
final Vector3f v01 = v1.sub(v0);
final Vector3f v02 = v2.sub(v0);
final Vector3f pvec = this.direction().cross(v02);
final float det = v01.dot(pvec);
if (det > 1e-6f) {
final Vector3f tvec = this.from().sub(v0);
final float u = tvec.dot(pvec);
if (u < 0.0f || u > det) {
return (false);
}
final Vector3f qvec = tvec.cross(v01);
final float v = this.direction().dot(qvec);
if (v < 0.0f || u + v > det) {
return (false);
}
final float t = v02.dot(qvec) / det;
this.setT(t);
return (true); // hit!!
}
return (false);
}
/**
* このオブジェクトの文字列表現を返します。
*
* @return このオブジェクトの文字列表現
*/
@Override
public String toString() {
String crlf = System.getProperty("line.separator");
StringBuilder sb = new StringBuilder();
sb.append("from : ");
sb.append(m_from.toString());
sb.append(crlf);
sb.append("direction : ");
sb.append(m_direction.toString());
sb.append(crlf);
sb.append("t : ");
sb.append(m_t);
return sb.toString();
}
public void setOrigin(final int win_x, final int win_y) {
final float[] in = {
win_x * m_delta.get(0) + m_constant.get(0), win_y * m_delta.get(1) + m_constant.get(1)
};
final float[] center = {
m_inverse.get(0, 0) * in[0] + m_inverse.get(0, 1) * in[1] + m_inverse.get(0, 3),
m_inverse.get(1, 0) * in[0] + m_inverse.get(1, 1) * in[1] + m_inverse.get(1, 3),
m_inverse.get(2, 0) * in[0] + m_inverse.get(2, 1) * in[1] + m_inverse.get(2, 3),
m_inverse.get(3, 0) * in[0] + m_inverse.get(3, 1) * in[1] + m_inverse.get(3, 3)
};
float[] front = {
center[0] - m_inverse.get(0, 2),
center[1] - m_inverse.get(1, 2),
center[2] - m_inverse.get(2, 2),
center[3] - m_inverse.get(3, 2)
};
final float front_inv = 1.0f / front[3];
front[0] *= front_inv;
front[1] *= front_inv;
front[2] *= front_inv;
float[] back = {
center[0] + m_inverse.get(0, 2),
center[1] + m_inverse.get(1, 2),
center[2] + m_inverse.get(2, 2),
center[3] + m_inverse.get(3, 2)
};
final float back_inv = 1.0f / back[3];
back[0] *= back_inv;
back[1] *= back_inv;
back[2] *= back_inv;
m_t = 0.0f;
m_from = new Vector3f(front[0], front[1], front[2]);
Vector3f direction = new Vector3f(back[0] - front[0], back[1] - front[1], back[2] - front[2]);
m_direction = direction.normalize();
}
public float depth() {
Vector3f point = this.point();
final float view2 =
point.getX() * m_combined.get(2, 0)
+ point.getY() * m_combined.get(2, 1)
+ point.getZ() * m_combined.get(2, 2)
+ m_combined.get(2, 3);
final float view3 =
point.getX() * m_combined.get(3, 0)
+ point.getY() * m_combined.get(3, 1)
+ point.getZ() * m_combined.get(3, 2)
+ m_combined.get(3, 3);
return ((1.0f + view2 / view3) * 0.5f);
}
private double interpolate_value(
final StructuredVolumeObject volume, final Vector3f vertex0, final Vector3f vertex1)
throws KVSException {
Buffer buf = volume.values();
final int line_size = volume.nnodesPerLine();
final int slice_size = volume.nnodesPerSlice();
final float value0 = this.substitute_plane_equation(vertex0);
final float value1 = this.substitute_plane_equation(vertex1);
final float ratio = kvs.core.util.Math.abs(value0 / (value1 - value0));
final int index0 =
(int) (vertex0.getX() + vertex0.getY() * line_size + vertex0.getZ() * slice_size);
final int index1 =
(int) (vertex1.getX() + vertex1.getY() * line_size + vertex1.getZ() * slice_size);
if (buf instanceof ByteBuffer) {
ByteBuffer values = (ByteBuffer) buf;
return (values.get(index0) + ratio * (values.get(index1) - values.get(index0)));
} else if (buf instanceof ShortBuffer) {
ShortBuffer values = (ShortBuffer) buf;
return (values.get(index0) + ratio * (values.get(index1) - values.get(index0)));
} else if (buf instanceof IntBuffer) {
IntBuffer values = (IntBuffer) buf;
return (values.get(index0) + ratio * (values.get(index1) - values.get(index0)));
} else if (buf instanceof LongBuffer) {
LongBuffer values = (LongBuffer) buf;
return (values.get(index0) + ratio * (values.get(index1) - values.get(index0)));
} else if (buf instanceof FloatBuffer) {
FloatBuffer values = (FloatBuffer) buf;
return (values.get(index0) + ratio * (values.get(index1) - values.get(index0)));
} else if (buf instanceof DoubleBuffer) {
DoubleBuffer values = (DoubleBuffer) buf;
return (values.get(index0) + ratio * (values.get(index1) - values.get(index0)));
} else if (buf instanceof CharBuffer) {
CharBuffer values = (CharBuffer) buf;
return (values.get(index0) + ratio * (values.get(index1) - values.get(index0)));
} else {
throw new KVSException("Unsupported data type");
}
}
public void wheel(final float value) {
Vector3f scale = new Vector3f(1.0f);
switch (getScalingType()) {
case ScalingXYZ:
scale = scale.mul(value);
break;
case ScalingX:
scale = new Vector3f(scale.getX() * value, scale.getY(), scale.getZ());
break;
case ScalingY:
scale = new Vector3f(scale.getX(), scale.getY() * value, scale.getZ());
break;
case ScalingZ:
scale = new Vector3f(scale.getX(), scale.getY(), scale.getZ() * value);
break;
case ScalingXY:
scale = new Vector3f(scale.getX() * value, scale.getY() * value, scale.getZ());
break;
case ScalingYZ:
scale = new Vector3f(scale.getX(), scale.getY() * value, scale.getZ() * value);
break;
case ScalingZX:
scale = new Vector3f(scale.getX() * value, scale.getY(), scale.getZ() * value);
break;
default:
break;
}
m_scale = scale;
scale(m_old, m_new, getScalingType());
m_old = m_new;
}
public Vector3f point() {
return (m_from.add(m_direction.mul(m_t)));
}
private void extract_plane(final StructuredVolumeObject volume) throws KVSException {
// Calculated the coordinate data array and the normal vector array.
ArrayList<Float> coords = new ArrayList<Float>();
ArrayList<Float> normals = new ArrayList<Float>();
ArrayList<Integer> colors = new ArrayList<Integer>();
// Calculate min/max values of the node data.
if (!volume.hasMinMaxValues()) {
volume.updateMinMaxValues();
}
// Calculate a normalize_factor.
double min_value = volume.minValue();
double max_value = volume.maxValue();
double normalize_factor = 255.0 / (max_value - min_value);
Vector3i ncells = volume.resolution().sub(new Vector3i(1));
int line_size = volume.nnodesPerLine();
ColorMap color_map = transferFunction().colorMap();
// Extract surfaces.
int index = 0;
for (int z = 0; z < ncells.getZ(); ++z) {
for (int y = 0; y < ncells.getY(); ++y) {
for (int x = 0; x < ncells.getX(); ++x) {
// Calculate the index of the reference table.
int table_index = this.calculate_table_index(x, y, z);
if (table_index == 0) {
continue;
}
if (table_index == 255) {
continue;
}
// Calculate the triangle polygons.
for (int i = 0; MarchingCubesTable.TriangleID[table_index][i] != -1; i += 3) {
// Refer the edge IDs from the TriangleTable by using the table_index.
int e0 = MarchingCubesTable.TriangleID[table_index][i];
int e1 = MarchingCubesTable.TriangleID[table_index][i + 2];
int e2 = MarchingCubesTable.TriangleID[table_index][i + 1];
// Determine vertices for each edge.
Vector3f v0 =
new Vector3f(
x + MarchingCubesTable.VertexID[e0][0][0],
y + MarchingCubesTable.VertexID[e0][0][1],
z + MarchingCubesTable.VertexID[e0][0][2]);
Vector3f v1 =
new Vector3f(
x + MarchingCubesTable.VertexID[e0][1][0],
y + MarchingCubesTable.VertexID[e0][1][1],
z + MarchingCubesTable.VertexID[e0][1][2]);
Vector3f v2 =
new Vector3f(
x + MarchingCubesTable.VertexID[e1][0][0],
y + MarchingCubesTable.VertexID[e1][0][1],
z + MarchingCubesTable.VertexID[e1][0][2]);
Vector3f v3 =
new Vector3f(
x + MarchingCubesTable.VertexID[e1][1][0],
y + MarchingCubesTable.VertexID[e1][1][1],
z + MarchingCubesTable.VertexID[e1][1][2]);
Vector3f v4 =
new Vector3f(
x + MarchingCubesTable.VertexID[e2][0][0],
y + MarchingCubesTable.VertexID[e2][0][1],
z + MarchingCubesTable.VertexID[e2][0][2]);
Vector3f v5 =
new Vector3f(
x + MarchingCubesTable.VertexID[e2][1][0],
y + MarchingCubesTable.VertexID[e2][1][1],
z + MarchingCubesTable.VertexID[e2][1][2]);
// Calculate coordinates of the vertices which are composed
// of the triangle polygon.
Vector3f vertex0 = this.interpolate_vertex(v0, v1);
coords.add(vertex0.getX());
coords.add(vertex0.getY());
coords.add(vertex0.getZ());
final Vector3f vertex1 = this.interpolate_vertex(v2, v3);
coords.add(vertex1.getX());
coords.add(vertex1.getY());
coords.add(vertex1.getZ());
final Vector3f vertex2 = this.interpolate_vertex(v4, v5);
coords.add(vertex2.getX());
coords.add(vertex2.getY());
coords.add(vertex2.getZ());
final double value0 = this.interpolate_value(volume, v0, v1);
final double value1 = this.interpolate_value(volume, v2, v3);
final double value2 = this.interpolate_value(volume, v4, v5);
final int color0 = (int) (normalize_factor * (value0 - min_value));
colors.add(color_map.getAt(color0).getRed());
colors.add(color_map.getAt(color0).getGreen());
colors.add(color_map.getAt(color0).getBlue());
final int color1 = (int) (normalize_factor * (value1 - min_value));
colors.add(color_map.getAt(color1).getRed());
colors.add(color_map.getAt(color1).getGreen());
colors.add(color_map.getAt(color1).getBlue());
final int color2 = (int) (normalize_factor * (value2 - min_value));
colors.add(color_map.getAt(color2).getRed());
colors.add(color_map.getAt(color2).getGreen());
colors.add(color_map.getAt(color2).getBlue());
// Calculate a normal vector for the triangle polygon.
final Vector3f normal = (vertex2.sub(vertex0)).cross(vertex1.sub(vertex0));
normals.add(normal.getX());
normals.add(normal.getY());
normals.add(normal.getZ());
} // end of loop-triangle
} // end of loop-x
++index;
} // end of loop-y
index += line_size;
} // end of loop-z
m_object.setCoords(Utility.ListToFloatArray(coords));
m_object.setColors(Utility.ListToIntArray(colors));
m_object.setNormals(Utility.ListToFloatArray(normals));
m_object.setOpacity((byte) 255);
m_object.setPolygonType(PolygonObject.PolygonType.Triangle);
m_object.setColorType(PolygonObject.ColorType.VertexColor);
m_object.setNormalType(PolygonObject.NormalType.PolygonNormal);
}
public SlicePlane(Vector3f point, Vector3f norm, TransferFunction transfer_function) {
super(transfer_function);
m_coefficients = new Vector4f(norm, -point.dot(norm));
}
public SlicePlane(Vector3f point, Vector3f norm) {
super();
m_coefficients = new Vector4f(norm, -point.dot(norm));
}
private float substitute_plane_equation(final Vector3f vertex) {
return (m_coefficients.x() * vertex.getX()
+ m_coefficients.y() * vertex.getY()
+ m_coefficients.z() * vertex.getZ()
+ m_coefficients.w());
}
private void extract_plane(final UnstructuredVolumeObject volume) throws KVSException {
// Calculated the coordinate data array and the normal vector array.
ArrayList<Float> coords = new ArrayList<Float>();
ArrayList<Float> normals = new ArrayList<Float>();
ArrayList<Integer> colors = new ArrayList<Integer>();
// Calculate min/max values of the node data.
if (!volume.hasMinMaxValues()) {
volume.updateMinMaxValues();
}
// Calculate a normalize factor.
double min_value = volume.minValue();
double max_value = volume.maxValue();
double normalize_factor = 255.0 / (max_value - min_value);
// Refer the parameters of the unstructured volume object.
final float[] volume_coords = volume.coords();
final int[] volume_connections = volume.connections();
final int ncells = volume.ncells();
final ColorMap color_map = transferFunction().colorMap();
// Extract surfaces.
int index = 0;
int[] local_index = new int[4];
for (int cell = 0; cell < ncells; ++cell, index += 4) {
// Calculate the indices of the target cell.
local_index[0] = volume_connections[index];
local_index[1] = volume_connections[index + 1];
local_index[2] = volume_connections[index + 2];
local_index[3] = volume_connections[index + 3];
// Calculate the index of the reference table.
final int table_index = this.calculate_table_index(local_index);
if (table_index == 0) {
continue;
}
if (table_index == 15) {
continue;
}
// Calculate the triangle polygons.
for (int i = 0; MarchingTetrahedraTable.TriangleID[table_index][i] != -1; i += 3) {
// Refer the edge IDs from the TriangleTable using the table_index.
final int e0 = MarchingTetrahedraTable.TriangleID[table_index][i];
final int e1 = MarchingTetrahedraTable.TriangleID[table_index][i + 1];
final int e2 = MarchingTetrahedraTable.TriangleID[table_index][i + 2];
// Refer indices of the coordinate array from the VertexTable using the edgeIDs.
final int c0 = local_index[MarchingTetrahedraTable.VertexID[e0][0]];
final int c1 = local_index[MarchingTetrahedraTable.VertexID[e0][1]];
final int c2 = local_index[MarchingTetrahedraTable.VertexID[e1][0]];
final int c3 = local_index[MarchingTetrahedraTable.VertexID[e1][1]];
final int c4 = local_index[MarchingTetrahedraTable.VertexID[e2][0]];
final int c5 = local_index[MarchingTetrahedraTable.VertexID[e2][1]];
// Determine vertices for each edge.
final Vector3f v0 = new Vector3f(volume_coords, 3 * c0);
final Vector3f v1 = new Vector3f(volume_coords, 3 * c1);
final Vector3f v2 = new Vector3f(volume_coords, 3 * c2);
final Vector3f v3 = new Vector3f(volume_coords, 3 * c3);
final Vector3f v4 = new Vector3f(volume_coords, 3 * c4);
final Vector3f v5 = new Vector3f(volume_coords, 3 * c5);
// Calculate coordinates of the vertices which are composed
// of the triangle polygon.
final Vector3f vertex0 = this.interpolate_vertex(v0, v1);
coords.add(vertex0.getX());
coords.add(vertex0.getY());
coords.add(vertex0.getZ());
final Vector3f vertex1 = this.interpolate_vertex(v2, v3);
coords.add(vertex1.getX());
coords.add(vertex1.getY());
coords.add(vertex1.getZ());
final Vector3f vertex2 = this.interpolate_vertex(v4, v5);
coords.add(vertex2.getX());
coords.add(vertex2.getY());
coords.add(vertex2.getZ());
final double value0 = this.interpolate_value(volume, c0, c1);
final double value1 = this.interpolate_value(volume, c2, c3);
final double value2 = this.interpolate_value(volume, c4, c5);
final int color0 = (int) (normalize_factor * (value0 - min_value));
colors.add(color_map.getAt(color0).getRed());
colors.add(color_map.getAt(color0).getGreen()); // blue?
colors.add(color_map.getAt(color0).getBlue());
final int color1 = (int) (normalize_factor * (value1 - min_value));
colors.add(color_map.getAt(color1).getRed());
colors.add(color_map.getAt(color1).getGreen());
colors.add(color_map.getAt(color1).getBlue());
final int color2 = (int) (normalize_factor * (value2 - min_value));
colors.add(color_map.getAt(color2).getRed());
colors.add(color_map.getAt(color2).getGreen());
colors.add(color_map.getAt(color2).getBlue());
// Calculate a normal vector for the triangle polygon.
final Vector3f normal = (vertex2.sub(vertex0)).cross(vertex1.sub(vertex0));
normals.add(normal.getX());
normals.add(normal.getY());
normals.add(normal.getZ());
} // end of loop-triangle
} // end of loop-cell
m_object.setCoords(Utility.ListToFloatArray(coords));
m_object.setColors(Utility.ListToIntArray(colors));
m_object.setNormals(Utility.ListToFloatArray(normals));
m_object.setOpacity((byte) 255);
m_object.setPolygonType(PolygonObject.PolygonType.Triangle);
m_object.setColorType(PolygonObject.ColorType.PolygonColor);
m_object.setNormalType(PolygonObject.NormalType.PolygonNormal);
}