31:Ted9,render = { var dx = 360 / isize; var dy = 180 / jsize; var kcurrent1 = Math.min(kcurrent, ksize); var lcurrent1 = Math.min(lcurrent, lsize); var vertexArray = new Float64Array(isize * jsize * 4 * 3); var polysArray = new Uint32Array(isize * jsize * 5); var id = 0, v = 0, p = 0; for (var j = 0; j < jsize; j++) { var y = axis2[j]; for (var i = 0; i < isize; i++) { var x = axis1[i]; 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 n = 0; n < cellBounds.length; n++) { var xyz = vertex(cellBounds[n], 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; id = id + 4; p = p + 5; } } var polydata = vtk.Common.DataModel.vtkPolyData.newInstance(); polydata.getPoints().setData(vertexArray, 3); polydata.getPolys().setData(polysArray); console.log(polydata.getCellData().getArrays()); var values = valuesAll .slice([lcurrent1 - 1, lcurrent1], [kcurrent1 - 1, kcurrent1], 0, 0) .flatten() .tolist(); polydata.getCellData().setScalars( vtk.Common.Core.vtkDataArray.newInstance({ name: varName, values: Float32Array.from(values) }) ); var graticule = createLineSegments(graticules(10)); var graticulepolygonmapper = vtk.Rendering.Core.vtkMapper.newInstance(); graticulepolygonmapper.setInputData(graticule); var graticulepolygonactor = vtk.Rendering.Core.vtkActor.newInstance(); graticulepolygonactor.getProperty().setColor(0.6, 0.6, 0.6); 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(); var preset = vtk.Rendering.Core.vtkColorTransferFunction.vtkColorMaps.getPresetByName( colorMap ); lookup.applyColorMap(preset); lookup.setNanColor(1.0, 1.0, 1.0, 1.0); lookup.build(); var mapper = vtk.Rendering.Core.vtkMapper.newInstance(); mapper.setInputData(polydata); mapper.setLookupTable(lookup); mapper.setScalarRange(levels[0], levels[1]); var actor = vtk.Rendering.Core.vtkActor.newInstance(); actor.setMapper(mapper); actor.getProperty().setEdgeVisibility(edge); // Needed when object changes invalidation.then(() => { actor.delete(); graticulepolygonactor.delete(); mapper.delete(); polydata.delete(); graticule.delete(); }); var camera = vtk.Rendering.Core.vtkCamera.newInstance(); camera.setClippingRange(282, 646); camera.setFocalPoint(0, 0, 0); camera.setPosition(327, -107, 307); camera.setViewUp(-0.64, 0.18, 0.75); var renderer = fullScreenRenderer.getRenderer(); renderer.addActor(actor); renderer.addActor(graticulepolygonactor); renderer.addActor(continentspolygonactor); renderer.setBackground(0.9, 0.9, 0.9); if (activeCamera) { renderer.setActiveCamera(camera); mutable activeCamera = false; } //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:T7c8,function readVariable(varName) { var varArrays = reader.variables; var dimIndex, varAttributes, varLongName, varMissingValue, varRecord, varCoordinates; varArrays.forEach(function(varArray) { if (varArray.name == varName) { dimIndex = varArray.dimensions; varAttributes = varArray.attributes; varLongName = varName; varMissingValue = 1E+34; varRecord = varArray.record; varCoordinates = 'None'; varAttributes.forEach(function(varAttribute) { if (varAttribute.name == 'long_name') { varLongName = varAttribute.value; } else if (varAttribute.name == 'missing_value') { varMissingValue = varAttribute.value; } else if (varAttribute.name == 'coordinates') { varCoordinates = varAttribute.value; } }); } }); var dim1 = reader.dimensions[dimIndex[dimIndex.length - 1]]; var isize = dim1.size; var dim2 = reader.dimensions[dimIndex[dimIndex.length - 2]]; var jsize = dim2.size; if (dimIndex.length == 2) { var ksize = 1; var lsize = 1; } else if (dimIndex.length == 3) { var dim3 = reader.dimensions[dimIndex[dimIndex.length - 3]]; if (dim3.size == 0) { dim3 = reader.recordDimension; ksize = 1; lsize = dim3.length; } else { ksize = dim3.size; lsize = 1; } } else if (dimIndex.length == 4) { dim3 = reader.dimensions[dimIndex[dimIndex.length - 3]]; ksize = dim3.size; var dim4 = reader.dimensions[dimIndex[dimIndex.length - 4]]; if (dim4.size == 0) { dim4 = reader.recordDimension; lsize = dim4.length; } else { lsize = dim4.size; } } var values = reader.getDataVariable(varName); var valuesAll = nj.array(values).reshape(lsize, ksize, jsize, isize); var axis1 = reader.getDataVariable(dim1.name); var axis2 = reader.getDataVariable(dim2.name); return { isize, jsize, ksize, lsize, varLongName, valuesAll, axis1, axis2 }; }34:[{"id":"b2896732353d3c4292a82cff2335cc93f29bdf765080ad00892ecaad960da2c715460151db9a17e86812df136f01a868db4bbecac2b36036be448d4c67e4538e","href":"https://static.observableusercontent.com/files/b2896732353d3c4292a82cff2335cc93f29bdf765080ad00892ecaad960da2c715460151db9a17e86812df136f01a868db4bbecac2b36036be448d4c67e4538e","name":"50m@1.json","type":"application/json","lastModified":1599766177030.199,"size":747178}] 2f:{"title":"VTK Mapper - Using vtks.js with netcdf input files","cells":[{"id":0,"value":"md`# VTK Mapper - Using vtks.js with netcdf input files\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":1054,"value":"viewof lcurrent = Player(Array.from({ length: lsize }, (_, i) => i + 1), {\n index: false,\n delay: 500,\n autoplay: false,\n format: s => s.toString().padStart(3, '0')\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":987,"value":"html`

Source

`","pinned":false,"hidden":false,"mode":"ojs"},{"id":351,"value":"viewof source = html``","pinned":false,"hidden":false,"mode":"ojs"},{"id":937,"value":"arrayVariables = reader.variables.filter(function(a) {\n if (\n a.dimensions.length >= 2 &&\n a.name.includes('bnds') == false &&\n a.name.includes('bounds') == false\n )\n return a;\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":750,"value":"html`

Variable

`","pinned":false,"hidden":false,"mode":"ojs"},{"id":739,"value":"viewof varName = html``","pinned":false,"hidden":false,"mode":"ojs"},{"id":869,"value":"html`

K index

from 1 to ${ksize}`","pinned":false,"hidden":false,"mode":"ojs"},{"id":872,"value":"viewof kcurrent = slider({\n min: 1,\n max: ksize,\n value: 1,\n step: 1\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":581,"value":"html`

ColorMap

`","pinned":false,"hidden":false,"mode":"ojs"},{"id":409,"value":"viewof colorMap = html``","pinned":false,"hidden":false,"mode":"ojs"},{"id":478,"value":"viewof levels = rangeSlider({\n min: -50,\n max: 50,\n step: 1,\n value: [-15, 15],\n title: '

Levels'\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1012,"value":"html`

Wireframe

`","pinned":false,"hidden":false,"mode":"ojs"},{"id":993,"value":"viewof edge = checkbox({\n options: [{ value: \"visible\", label: \"

Visible\" }],\n value: \"toggle\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":807,"value":"({\n min: valuesAll.min(),\n max: valuesAll.max()\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":323,"value":"render()","pinned":false,"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":561,"value":"$31","pinned":false,"hidden":false,"mode":"ojs"},{"id":1303,"value":"mutable activeCamera = true","pinned":true,"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":349,"value":"reader = new netcdf(\n await fetch(source).then(function(response) {\n return response.arrayBuffer();\n })\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":757,"value":"readVariable(varName)","pinned":false,"hidden":false,"mode":"ojs"},{"id":731,"value":"varLongName = readVariable(varName).varLongName","pinned":true,"hidden":false,"mode":"ojs"},{"id":603,"value":"valuesAll = readVariable(varName).valuesAll","pinned":true,"hidden":false,"mode":"ojs"},{"id":665,"value":"isize = readVariable(varName).isize","pinned":true,"hidden":false,"mode":"ojs"},{"id":668,"value":"jsize = readVariable(varName).jsize","pinned":true,"hidden":false,"mode":"ojs"},{"id":795,"value":"ksize = readVariable(varName).ksize","pinned":true,"hidden":false,"mode":"ojs"},{"id":797,"value":"lsize = readVariable(varName).lsize","pinned":true,"hidden":false,"mode":"ojs"},{"id":978,"value":"axis1 = readVariable(varName).axis1","pinned":true,"hidden":false,"mode":"ojs"},{"id":980,"value":"axis2 = readVariable(varName).axis2","pinned":true,"hidden":false,"mode":"ojs"},{"id":589,"value":"$33","pinned":false,"hidden":false,"mode":"ojs"},{"id":407,"value":"html``","pinned":false,"hidden":false,"mode":"ojs"},{"id":481,"value":"import { rangeSlider } from '@mootari/range-slider'","pinned":false,"hidden":false,"mode":"ojs"},{"id":655,"value":"nj = require('https://cdn.jsdelivr.net/gh/nicolaspanel/numjs@0.16/dist/numjs.min.js')","pinned":false,"hidden":false,"mode":"ojs"},{"id":329,"value":"netcdf = require('https://unpkg.com/netcdfjs@1.0.0/dist/netcdfjs.min.js')","pinned":false,"hidden":false,"mode":"ojs"},{"id":68,"value":"topojson = require(\"topojson\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":42,"value":"d3 = require(\"d3\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":842,"value":"import { slider, checkbox } from \"@jashkenas/inputs\"","pinned":true,"hidden":false,"mode":"ojs"},{"id":1019,"value":"import { animationControls } from \"@pbrockmann/animation-controls\"","pinned":true,"hidden":false,"mode":"ojs"},{"id":1055,"value":"import { Player } from \"@oscar6echo/player\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":1064,"value":"import { style as faStyle } from \"@airbornemint/fontawesome\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":1066,"value":"faStyle({ solid: true })","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W34","stdlib":"1","titling":"auto"}