31:Tddf,render = { var dx = d; var dy = d; var vertexArray = new Float64Array((360 * 180 * 4 * 3) / (dx * dy)); var polysArray = new Uint32Array((360 * 180 * 5) / (dx * dy)); var valuesArray1 = new Float64Array((360 * 180) / (dx * dy)); var valuesArray2 = new Float64Array((360 * 180) / (dx * dy)); var id = 0, n = 0, v = 0, p = 0; for (var y = -90 + dy / 2; y <= 90 - dy / 2; y = y + dy) { for (var x = -180 + dx / 2; x <= 180 - dx / 2; x = x + dx) { var cellBounds = [ [x - dx / 2, y - dy / 2], [x - dx / 2, y + dy / 2], [x + dx / 2, y + dy / 2], [x + dx / 2, y - dy / 2] ]; for (var i = 0; i < 4; i++) { var xyz = vertex(cellBounds[i], 1.0); vertexArray[v] = xyz[0]; vertexArray[v + 1] = xyz[1]; vertexArray[v + 2] = xyz[2]; v = v + 3; } polysArray[p] = 4; polysArray[p + 1] = id; polysArray[p + 2] = id + 1; polysArray[p + 3] = id + 2; polysArray[p + 4] = id + 3; valuesArray1[n] = Math.random() * 10; valuesArray2[n] = Math.random() * 10; id = id + 4; p = p + 5; n = n + 1; } } var polydata = vtk.Common.DataModel.vtkPolyData.newInstance(); polydata.getPoints().setData(vertexArray, 3); polydata.getPolys().setData(polysArray); console.log(polydata.getCellData().getArrays()); polydata.getCellData().setScalars( vtk.Common.Core.vtkDataArray.newInstance({ name: 'var1', values: valuesArray1 }) ); console.log('Nb of points: ', polydata.getPoints().getData().length); console.log( 'Nb of cells: ', polydata .getCellData() .getScalars() .getData().length ); var graticule = createLineSegments(graticules(deltaGraticule)); var graticulepolygonmapper = vtk.Rendering.Core.vtkMapper.newInstance(); graticulepolygonmapper.setInputData(graticule); var graticulepolygonactor = vtk.Rendering.Core.vtkActor.newInstance(); graticulepolygonactor.getProperty().setColor(1, 0, 0); graticulepolygonactor.getProperty().setAmbient(1); graticulepolygonactor.getProperty().setDiffuse(0); graticulepolygonactor.setMapper(graticulepolygonmapper); var continentspolygonmapper = vtk.Rendering.Core.vtkMapper.newInstance(); continentspolygonmapper.setInputData(continents); var continentspolygonactor = vtk.Rendering.Core.vtkActor.newInstance(); continentspolygonactor.getProperty().setColor(0, 0, 0); continentspolygonactor.getProperty().setAmbient(1); continentspolygonactor.getProperty().setDiffuse(0); continentspolygonactor.setMapper(continentspolygonmapper); var lookup = vtk.Rendering.Core.vtkColorTransferFunction.newInstance(); lookup.addRGBPoint(0, 1.0, 0.0, 1.0); //red lookup.addRGBPoint(1, 1.0, 1.0, 1.0); //white lookup.addRGBPoint(2, 0.0, 1.0, 0.0); //green lookup.build(); var mapper = vtk.Rendering.Core.vtkMapper.newInstance(); mapper.setInputData(polydata); mapper.setLookupTable(lookup); mapper.setScalarRange(0, 10); var actor = vtk.Rendering.Core.vtkActor.newInstance(); actor.setMapper(mapper); // Needed when object changes invalidation.then(() => { actor.delete(); graticulepolygonactor.delete(); mapper.delete(); polydata.delete(); graticule.delete(); }); var renderer = fullScreenRenderer.getRenderer(); renderer.addActor(actor); renderer.addActor(graticulepolygonactor); renderer.addActor(continentspolygonactor); renderer.resetCamera(); return fullScreenRenderer.getRenderWindow().render; }32:T531,function graticules(delta) { var epsilon = 1e-6, x1 = 180, x0 = -x1, y1 = 80, y0 = -y1, dx = delta, dy = delta, X1 = 180, X0 = -X1, Y1 = 90, Y0 = -Y1, DX = 90, DY = 360, x = graticuleX(y0, y1, 2.5), y = graticuleY(x0, x1, 2.5), X = graticuleX(Y0, Y1, 2.5), Y = graticuleY(X0, X1, 2.5); function graticuleX(y0, y1, dy) { var y = d3.range(y0, y1 - epsilon, dy).concat(y1); return function(x) { return y.map(function(y) { return [x, y]; }); }; } function graticuleY(x0, x1, dx) { var x = d3.range(x0, x1 - epsilon, dx).concat(x1); return function(y) { return x.map(function(x) { return [x, y]; }); }; } return { type: "MultiLineString", coordinates: d3 .range(Math.ceil(X0 / DX) * DX, X1, DX) .map(X) .concat(d3.range(Math.ceil(Y0 / DY) * DY, Y1, DY).map(Y)) .concat( d3 .range(Math.ceil(x0 / dx) * dx, x1, dx) .filter(function(x) { return Math.abs(x % DX) > epsilon; }) .map(x) ) .concat( d3 .range(Math.ceil(y0 / dy) * dy, y1 + epsilon, dy) .filter(function(y) { return Math.abs(y % DY) > epsilon; }) .map(y) ) }; }33:[{"id":"b2896732353d3c4292a82cff2335cc93f29bdf765080ad00892ecaad960da2c715460151db9a17e86812df136f01a868db4bbecac2b36036be448d4c67e4538e","href":"https://static.observableusercontent.com/files/b2896732353d3c4292a82cff2335cc93f29bdf765080ad00892ecaad960da2c715460151db9a17e86812df136f01a868db4bbecac2b36036be448d4c67e4538e","name":"50m@1.json","type":"application/json","lastModified":1599239411245.621,"size":747178}] 2f:{"title":"VTK Mapper is back !","cells":[{"id":0,"value":"md`# VTK Mapper is back !\nInitial work long time ago from https://github.com/PBrockmann/VTK_Mapper/blob/master/PRISM_Report19.pdf \nWould be nice to have it from an Observable notebook !`","pinned":false,"hidden":false,"mode":"ojs"},{"id":12,"value":"viewof fullScreenRenderer = {\n const view = html`

`;\n yield view;\n\n view.value = vtk.Rendering.Misc.vtkFullScreenRenderWindow.newInstance({\n rootContainer: view,\n containerStyle: { width: `${width}px`, height: \"600px\" }\n });\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":160,"value":"viewof deltaGraticule = html``","pinned":false,"hidden":false,"mode":"ojs"},{"id":287,"value":"html`DeltaGraticule: ${deltaGraticule}`","pinned":false,"hidden":false,"mode":"ojs"},{"id":260,"value":"viewof d = html``","pinned":false,"hidden":false,"mode":"ojs"},{"id":262,"value":"html`DeltaPolydata: ${d}`","pinned":false,"hidden":false,"mode":"ojs"},{"id":239,"value":"html`

`","pinned":false,"hidden":false,"mode":"ojs"},{"id":323,"value":"render()","pinned":true,"hidden":false,"mode":"ojs"},{"id":2,"value":"vtk = require(\"vtk.js\").then(() => window.vtk)","pinned":false,"hidden":false,"mode":"ojs"},{"id":27,"value":"function vertex(point, scaleFactor) {\n var radius = 100;\n var radius1 = radius * scaleFactor;\n var lambda = (point[0] * Math.PI) / 180,\n phi = (point[1] * Math.PI) / 180,\n cosPhi = Math.cos(phi);\n return [\n radius1 * cosPhi * Math.cos(lambda),\n radius1 * cosPhi * Math.sin(lambda),\n radius1 * Math.sin(phi)\n ];\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":16,"value":"$31","pinned":false,"hidden":false,"mode":"ojs"},{"id":131,"value":"continents = createLineSegments(topojson.mesh(world, world.objects.land))","pinned":false,"hidden":false,"mode":"ojs"},{"id":75,"value":"world = FileAttachment(\"50m@1.json\").json()","pinned":false,"hidden":false,"mode":"ojs"},{"id":40,"value":"$32","pinned":false,"hidden":false,"mode":"ojs"},{"id":46,"value":"function createLineSegments(multilinestring) {\n var polydata = vtk.Common.DataModel.vtkPolyData.newInstance();\n var vertexArray = [];\n var lineArray = [];\n var id = 0;\n for (var p = 0; p < multilinestring.coordinates.length; p++) {\n var line = multilinestring.coordinates[p];\n lineArray.push(line.length);\n for (var i = 0; i < line.length; i++) {\n var xyz = vertex(line[i], 1.005);\n vertexArray.push(xyz[0], xyz[1], xyz[2]);\n lineArray.push(id++);\n }\n polydata.getPoints().setData(Float32Array.from(vertexArray), 3);\n polydata.getLines().setData(Uint16Array.from(lineArray));\n }\n return polydata;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":42,"value":"d3 = require(\"d3\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":44,"value":"d3.version\n","pinned":false,"hidden":false,"mode":"ojs"},{"id":68,"value":"topojson = require(\"topojson\")","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W33","stdlib":"1","titling":"auto"}