Skip to content

Handle streamtube coordinates in string format and handle different data orders in isosurface and volume #4431

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
merged 11 commits into from
Dec 30, 2019
Merged

Handle streamtube coordinates in string format and handle different data orders in isosurface and volume #4431

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
11 commits
Select commit Hold shift + click to select a range
07ce3d2
fix issue 4427 - streamtube should handle string values
archmoj Dec 18, 2019
942f92f
streamtube - no need for casting typed arrays to numbers
archmoj Dec 21, 2019
d015b57
create a reusable process grid function
archmoj Dec 24, 2019
5a95dab
Apply streamtube processGrid in isosurface and volume
archmoj Dec 24, 2019
1f0ee0f
improve test image for isosurface
archmoj Dec 25, 2019
990cdda
match isosurface value in the mock to be equal to sqrt(u2+v2+w2)
archmoj Dec 25, 2019
b71677a
pass len to distictVals and filter to help improve the performance
archmoj Dec 27, 2019
6d06229
add an image test for volume with different fill orders
archmoj Dec 27, 2019
df9b4c4
cast x, y, z, u, v, w, values to numbers at the start of calc and sto…
archmoj Dec 27, 2019
41504aa
remove extra cast operations
archmoj Dec 29, 2019
c82da07
rewrite without forEach
archmoj Dec 30, 2019
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
27 changes: 22 additions & 5 deletions src/traces/isosurface/calc.js
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -9,22 +9,39 @@
'use strict';

var colorscaleCalc = require('../../components/colorscale/calc');
var processGrid = require('../streamtube/calc').processGrid;
var filter = require('../streamtube/calc').filter;

module.exports = function calc(gd, trace) {
trace._len = Math.min(trace.x.length, trace.y.length, trace.z.length, trace.value.length);
trace._len = Math.min(
trace.x.length,
trace.y.length,
trace.z.length,
trace.value.length
);

trace._x = filter(trace.x, trace._len);
trace._y = filter(trace.y, trace._len);
trace._z = filter(trace.z, trace._len);
trace._value = filter(trace.value, trace._len);

var grid = processGrid(trace);
trace._gridFill = grid.fill;
trace._Xs = grid.Xs;
trace._Ys = grid.Ys;
trace._Zs = grid.Zs;
trace._len = grid.len;

var min = Infinity;
var max = -Infinity;
var len = trace.value.length;
for(var i = 0; i < len; i++) {
var v = trace.value[i];
for(var i = 0; i < trace._len; i++) {
var v = trace._value[i];
min = Math.min(min, v);
max = Math.max(max, v);
}

trace._minValues = min;
trace._maxValues = max;

trace._vMin = (trace.isomin === undefined || trace.isomin === null) ? min : trace.isomin;
trace._vMax = (trace.isomax === undefined || trace.isomin === null) ? max : trace.isomax;

Expand Down
96 changes: 53 additions & 43 deletions src/traces/isosurface/convert.js
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -9,18 +9,11 @@
'use strict';

var createMesh = require('gl-mesh3d');

var Lib = require('../../lib');

var parseColorScale = require('../../lib/gl_format_color').parseColorScale;
var str2RgbaArray = require('../../lib/str2rgbarray');
var extractOpts = require('../../components/colorscale').extractOpts;
var zip3 = require('../../plots/gl3d/zip3');

function distinctVals(col) {
return Lib.distinctVals(col).vals;
}

var findNearestOnAxis = function(w, arr) {
for(var q = arr.length - 1; q > 0; q--) {
var min = Math.min(arr[q], arr[q - 1]);
Expand Down Expand Up @@ -53,9 +46,9 @@ proto.handlePick = function(selection) {
if(selection.object === this.mesh) {
var rawId = selection.data.index;

var x = this.data._x[rawId];
var y = this.data._y[rawId];
var z = this.data._z[rawId];
var x = this.data._meshX[rawId];
var y = this.data._meshY[rawId];
var z = this.data._meshZ[rawId];

var height = this.data._Ys.length;
var depth = this.data._Zs.length;
Expand All @@ -67,10 +60,10 @@ proto.handlePick = function(selection) {
var selectIndex = selection.index = k + depth * j + depth * height * i;

selection.traceCoordinate = [
this.data._x[selectIndex],
this.data._y[selectIndex],
this.data._z[selectIndex],
this.data.value[selectIndex]
this.data._meshX[selectIndex],
this.data._meshY[selectIndex],
this.data._meshZ[selectIndex],
this.data._value[selectIndex]
];

var text = this.data.hovertext || this.data.text;
Expand Down Expand Up @@ -98,11 +91,11 @@ proto.update = function(data) {
}

var positions = zip3(
toDataCoords(layout.xaxis, data._x, scene.dataScale[0], data.xcalendar),
toDataCoords(layout.yaxis, data._y, scene.dataScale[1], data.ycalendar),
toDataCoords(layout.zaxis, data._z, scene.dataScale[2], data.zcalendar));
toDataCoords(layout.xaxis, data._meshX, scene.dataScale[0], data.xcalendar),
toDataCoords(layout.yaxis, data._meshY, scene.dataScale[1], data.ycalendar),
toDataCoords(layout.zaxis, data._meshZ, scene.dataScale[2], data.zcalendar));

var cells = zip3(data._i, data._j, data._k);
var cells = zip3(data._meshI, data._meshJ, data._meshK);

var config = {
positions: positions,
Expand All @@ -123,7 +116,7 @@ proto.update = function(data) {
};

var cOpts = extractOpts(data);
config.vertexIntensity = data._intensity;
config.vertexIntensity = data._meshIntensity;
config.vertexIntensityBounds = [cOpts.min, cOpts.max];
config.colormap = parseColorScale(data);

Expand All @@ -136,10 +129,12 @@ proto.dispose = function() {
this.mesh.dispose();
};

var GRID_TYPES = ['xyz', 'xzy', 'yxz', 'yzx', 'zxy', 'zyx'];

function generateIsoMeshes(data) {
data._i = [];
data._j = [];
data._k = [];
data._meshI = [];
data._meshJ = [];
data._meshK = [];

var showSurface = data.surface.show;
var showSpaceframe = data.spaceframe.show;
Expand All @@ -154,17 +149,32 @@ function generateIsoMeshes(data) {
var numVertices;
var beginVertextLength;

var Xs = distinctVals(data.x.slice(0, data._len));
var Ys = distinctVals(data.y.slice(0, data._len));
var Zs = distinctVals(data.z.slice(0, data._len));
var Xs = data._Xs;
var Ys = data._Ys;
var Zs = data._Zs;

var width = Xs.length;
var height = Ys.length;
var depth = Zs.length;

function getIndex(i, j, k) {
return k + depth * j + depth * height * i;
}
var filled = GRID_TYPES.indexOf(data._gridFill.replace(/-/g, '').replace(/\+/g, ''));

var getIndex = function(i, j, k) {
switch(filled) {
case 5: // 'zyx'
return k + depth * j + depth * height * i;
case 4: // 'zxy'
return k + depth * i + depth * width * j;
case 3: // 'yzx'
return j + height * k + height * depth * i;
case 2: // 'yxz'
return j + height * i + height * width * k;
case 1: // 'xzy'
return i + width * k + width * depth * j;
default: // case 0: // 'xyz'
return i + width * j + width * height * k;
}
};

var minValues = data._minValues;
var maxValues = data._maxValues;
Expand Down Expand Up @@ -218,9 +228,9 @@ function generateIsoMeshes(data) {
}

function addFace(a, b, c) {
data._i.push(a);
data._j.push(b);
data._k.push(c);
data._meshI.push(a);
data._meshJ.push(b);
data._meshK.push(c);
numFaces++;

return numFaces - 1;
Expand Down Expand Up @@ -379,10 +389,10 @@ function generateIsoMeshes(data) {
var index = indecies[q];
xyzv.push(
[
data.x[index],
data.y[index],
data.z[index],
data.value[index]
data._x[index],
data._y[index],
data._z[index],
data._value[index]
]
);
}
Expand Down Expand Up @@ -870,10 +880,10 @@ function generateIsoMeshes(data) {
for(var k = 0; k < depth; k++) {
var index = getIndex(i, j, k);
addVertex(
data.x[index],
data.y[index],
data.z[index],
data.value[index]
data._x[index],
data._y[index],
data._z[index],
data._value[index]
);
}
}
Expand Down Expand Up @@ -1011,10 +1021,10 @@ function generateIsoMeshes(data) {
emptyVertices();
}

data._x = allXs;
data._y = allYs;
data._z = allZs;
data._intensity = allVs;
data._meshX = allXs;
data._meshY = allYs;
data._meshZ = allZs;
data._meshIntensity = allVs;

data._Xs = Xs;
data._Ys = Ys;
Expand Down
Loading