2d:null 2e:["$","title",null,{"children":"VTK Mapper is back ! | Patrick Brockmann | Observable"}] 33: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; }34: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) ) }; }35:[{"id":"b2896732353d3c4292a82cff2335cc93f29bdf765080ad00892ecaad960da2c715460151db9a17e86812df136f01a868db4bbecac2b36036be448d4c67e4538e","href":"https://static.observableusercontent.com/files/b2896732353d3c4292a82cff2335cc93f29bdf765080ad00892ecaad960da2c715460151db9a17e86812df136f01a868db4bbecac2b36036be448d4c67e4538e","name":"50m@1.json","type":"application/json","lastModified":1599239411245.621,"size":747178}] 31:{"id":"42482f44a1477d56","slug":"vtk-mapper-is-back","title":"VTK Mapper is back !","thumbnail":"4d59f3845d63ba17bc5727f7668fd2442f95535c0121ce6737574e52437bbc38","access_level":"public","create_time":"$D2020-09-04T13:06:29.763Z","update_time":"$D2020-09-10T19:32:07.545Z","view_time":null,"first_public_time":"$D2020-09-04T17:26:07.356Z","has_unpublished_edits":null,"fork":false,"stars":1,"starred":false,"owner":{"type":"team","id":"63e96f00663997c7","login":"pbrockmann","name":"Patrick Brockmann","avatar_url":"https://avatars-next.observableusercontent.com/25f3e1483d611e61cc2dc87b5091004972f4ad8dbf4792cd225809597ddb2463","tier":"starter_2024"},"version":343,"publish_version":343,"role":null,"authored":false,"authors":[{"type":"individual","id":"68e18af682792012","login":"pbrockmann","name":"Patrick Brockmann","avatar_url":"https://avatars-next.observableusercontent.com/25f3e1483d611e61cc2dc87b5091004972f4ad8dbf4792cd225809597ddb2463","tier":null}],"members":[{"type":"individual","id":"68e18af682792012","login":"pbrockmann","name":"Patrick Brockmann","avatar_url":"https://avatars-next.observableusercontent.com/25f3e1483d611e61cc2dc87b5091004972f4ad8dbf4792cd225809597ddb2463","tier":null,"role":"editor"}],"comments":[{"id":"12c160ac3ec4f4fb","reply_to_id":null,"content":"Get a \"render = TypeError: Cannot read property 'getBounds' of undefined\" while updating delta of the grid of the polydata.\nDo not understand why.","create_time":"$D2020-09-09T20:30:11.745Z","update_time":null,"resolved":true,"user":{"type":"individual","id":"68e18af682792012","login":"pbrockmann","name":"Patrick Brockmann","avatar_url":"https://avatars-next.observableusercontent.com/25f3e1483d611e61cc2dc87b5091004972f4ad8dbf4792cd225809597ddb2463","tier":null}},{"id":"43aa6ec861cfa6e2","reply_to_id":"12c160ac3ec4f4fb","content":"I think it's because you deleted graticule which is an outside object, not one created internally. Comment it out:\n// graticule.delete();","create_time":"$D2020-09-09T20:40:43.993Z","update_time":null,"resolved":true,"user":{"type":"individual","id":"45a379fcfcb14253","login":"fil","name":"Fil","avatar_url":"https://avatars-next.observableusercontent.com/d5b1777907117a2bdf2a07eb98600d7a9beb21639fde42b9753e95c88c7731a6","tier":null}}],"custom_thumbnail":false,"trashed":false,"lapsed":null,"subscribed":null,"forks":1,"collection":null,"type":"notebook","body":{"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","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"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","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":160,"value":"viewof deltaGraticule = html``","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":287,"value":"html`DeltaGraticule: ${deltaGraticule}`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":260,"value":"viewof d = html``","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":262,"value":"html`DeltaPolydata: ${d}`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":239,"value":"html`

`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":323,"value":"render()","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2,"value":"vtk = require(\"vtk.js\").then(() => window.vtk)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"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","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":16,"value":"$33","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":131,"value":"continents = createLineSegments(topojson.mesh(world, world.objects.land))","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":75,"value":"world = FileAttachment(\"50m@1.json\").json()","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":40,"value":"$34","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"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","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":42,"value":"d3 = require(\"d3\")","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":44,"value":"d3.version\n","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":68,"value":"topojson = require(\"topojson\")","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"}],"files":"$W35","stdlib":"1","titling":"auto"}} 32:"$31:body"