itowns 2.43.2-next.0 → 2.43.2-next.10

package/lib/Parser/LASLoader.js

@@ -15,6 +15,10 @@ import { Las } from 'copc';
  * xOffset, zOffset]`) added to the scaled X, Y, Z point record values.
  */
 
+function defaultColorEncoding(header) {
+  return header.majorVersion === 1 && header.minorVersion <= 2 ? 8 : 16;
+}
+
 /**
  * @classdesc
  * Loader for LAS and LAZ (LASZip) point clouds. It uses the copc.js library and
@@ -117,6 +121,46 @@ class LASLoader {
     this._wasmPromise = null;
   }
 
+  /**
+   * Parses a LAS or LAZ (LASZip) chunk. Note that this function is
+   * **CPU-bound** and shall be parallelised in a dedicated worker.
+   * @param {Uint8Array} data - File chunk data.
+   * @param {Object} options - Parsing options.
+   * @param {Header} options.header - Partial LAS header.
+   * @param {number} options.pointCount - Number of points encoded in this
+   * data chunk.
+   * @param {Las.ExtraBytes[]} [options.eb] - Extra bytes LAS VLRs
+   * headers.
+   * @param {8 | 16} [options.colorDepth] - Color depth encoding (in bits).
+   * Either 8 or 16 bits. Defaults to 8 bits for LAS 1.2 and 16 bits for later
+   * versions (as mandatory by the specification).
+   */
+  async parseChunk(data, options) {
+    const {
+      header,
+      eb,
+      pointCount
+    } = options;
+    const {
+      pointDataRecordFormat,
+      pointDataRecordLength
+    } = header;
+    const colorDepth = options.colorDepth ?? defaultColorEncoding(header);
+    const bytes = new Uint8Array(data);
+    const pointData = await Las.PointData.decompressChunk(bytes, {
+      pointCount,
+      pointDataRecordFormat,
+      pointDataRecordLength
+    }, this._initDecoder());
+    const view = Las.View.create(pointData, header, eb);
+    const attributes = this._parseView(view, {
+      colorDepth
+    });
+    return {
+      attributes
+    };
+  }
+
   /**
    * Parses a LAS or LAZ (LASZip) file. Note that this function is
    * **CPU-bound** and shall be parallelised in a dedicated worker.
@@ -131,7 +175,7 @@ class LASLoader {
     const bytes = new Uint8Array(data);
     const pointData = await Las.PointData.decompressFile(bytes, this._initDecoder());
     const header = Las.Header.parse(bytes);
-    const colorDepth = options.colorDepth ?? (header.majorVersion === 1 && header.minorVersion <= 2 ? 8 : 16);
+    const colorDepth = options.colorDepth ?? defaultColorEncoding(header);
     const getter = async (begin, end) => bytes.slice(begin, end);
     const vlrs = await Las.Vlr.walk(getter, header);
     const ebVlr = Las.Vlr.find(vlrs, 'LASF_Spec', 4);

package/lib/Parser/LASParser.js

@@ -1,6 +1,29 @@
 import * as THREE from 'three';
 import LASLoader from "./LASLoader.js";
 const lasLoader = new LASLoader();
+function buildBufferGeometry(attributes) {
+  const geometry = new THREE.BufferGeometry();
+  const positionBuffer = new THREE.BufferAttribute(attributes.position, 3);
+  geometry.setAttribute('position', positionBuffer);
+  const intensityBuffer = new THREE.BufferAttribute(attributes.intensity, 1);
+  geometry.setAttribute('intensity', intensityBuffer);
+  const returnNumber = new THREE.BufferAttribute(attributes.returnNumber, 1);
+  geometry.setAttribute('returnNumber', returnNumber);
+  const numberOfReturns = new THREE.BufferAttribute(attributes.numberOfReturns, 1);
+  geometry.setAttribute('numberOfReturns', numberOfReturns);
+  const classBuffer = new THREE.BufferAttribute(attributes.classification, 1);
+  geometry.setAttribute('classification', classBuffer);
+  const pointSourceID = new THREE.BufferAttribute(attributes.pointSourceID, 1);
+  geometry.setAttribute('pointSourceID', pointSourceID);
+  if (attributes.color) {
+    const colorBuffer = new THREE.BufferAttribute(attributes.color, 4, true);
+    geometry.setAttribute('color', colorBuffer);
+  }
+  const scanAngle = new THREE.BufferAttribute(attributes.scanAngle, 1);
+  geometry.setAttribute('scanAngle', scanAngle);
+  geometry.userData.origin = new THREE.Vector3().fromArray(attributes.origin);
+  return geometry;
+}
 
 /** The LASParser module provides a [parse]{@link
  * module:LASParser.parse} method that takes a LAS or LAZ (LASZip) file in, and
@@ -21,6 +44,36 @@ export default {
     }
     lasLoader.lazPerf = path;
   },
+  /**
+   * Parses a chunk of a LAS or LAZ (LASZip) and returns the corresponding
+   * `THREE.BufferGeometry`.
+   *
+   * @param {ArrayBuffer} data - The file content to parse.
+   * @param {Object} options
+   * @param {Object} options.in - Options to give to the parser.
+   * @param {number} options.in.pointCount - Number of points encoded in this
+   * data chunk.
+   * @param {Object} options.in.header - Partial LAS file header.
+   * @param {number} options.in.header.pointDataRecordFormat - Type of Point
+   * Data Record contained in the LAS file.
+   * @param {number} options.in.header.pointDataRecordLength - Size (in bytes)
+   * of the Point Data Record.
+   * @param {Object} [options.eb] - Extra bytes LAS VLRs headers.
+   * @param { 8 | 16 } [options.in.colorDepth] - Color depth (in bits).
+   * Defaults to 8 bits for LAS 1.2 and 16 bits for later versions
+   * (as mandatory by the specification)
+   *
+   * @return {Promise<THREE.BufferGeometry>} A promise resolving with a
+   * `THREE.BufferGeometry`.
+   */
+  parseChunk(data) {
+    let options = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : {};
+    return lasLoader.parseChunk(data, options.in).then(parsedData => {
+      const geometry = buildBufferGeometry(parsedData.attributes);
+      geometry.computeBoundingBox();
+      return geometry;
+    });
+  },
   /**
    * Parses a LAS file or a LAZ (LASZip) file and return the corresponding
    * `THREE.BufferGeometry`.
@@ -36,37 +89,16 @@ export default {
    * header of the file is contained in `userData`.
    */
   parse(data) {
-    var _options$out, _options$in;
     let options = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : {};
-    if ((_options$out = options.out) !== null && _options$out !== void 0 && _options$out.skip) {
+    if (options.out?.skip) {
       console.warn("Warning: options 'skip' not supported anymore");
     }
     return lasLoader.parseFile(data, {
-      colorDepth: (_options$in = options.in) === null || _options$in === void 0 ? void 0 : _options$in.colorDepth
+      colorDepth: options.in?.colorDepth
     }).then(parsedData => {
-      const geometry = new THREE.BufferGeometry();
-      const attributes = parsedData.attributes;
-      geometry.userData = parsedData.header;
-      const positionBuffer = new THREE.BufferAttribute(attributes.position, 3);
-      geometry.setAttribute('position', positionBuffer);
-      const intensityBuffer = new THREE.BufferAttribute(attributes.intensity, 1);
-      geometry.setAttribute('intensity', intensityBuffer);
-      const returnNumber = new THREE.BufferAttribute(attributes.returnNumber, 1);
-      geometry.setAttribute('returnNumber', returnNumber);
-      const numberOfReturns = new THREE.BufferAttribute(attributes.numberOfReturns, 1);
-      geometry.setAttribute('numberOfReturns', numberOfReturns);
-      const classBuffer = new THREE.BufferAttribute(attributes.classification, 1);
-      geometry.setAttribute('classification', classBuffer);
-      const pointSourceID = new THREE.BufferAttribute(attributes.pointSourceID, 1);
-      geometry.setAttribute('pointSourceID', pointSourceID);
-      if (attributes.color) {
-        const colorBuffer = new THREE.BufferAttribute(attributes.color, 4, true);
-        geometry.setAttribute('color', colorBuffer);
-      }
-      const scanAngle = new THREE.BufferAttribute(attributes.scanAngle, 1);
-      geometry.setAttribute('scanAngle', scanAngle);
+      const geometry = buildBufferGeometry(parsedData.attributes);
+      geometry.userData.header = parsedData.header;
       geometry.computeBoundingBox();
-      geometry.userData.origin = new THREE.Vector3().fromArray(attributes.origin);
       return geometry;
     });
   }

package/lib/Parser/deprecated/LegacyGLTFLoader.js

@@ -321,7 +321,6 @@ threeExamples.LegacyGLTFLoader = function () {
     1029: THREE.FrontSide // Culling back
     //1032: THREE.NoSide   // Culling front and back, what to do?
   };
-
   var WEBGL_DEPTH_FUNCS = {
     512: THREE.NeverDepth,
     513: THREE.LessDepth,
@@ -355,7 +354,6 @@ threeExamples.LegacyGLTFLoader = function () {
     //32771: CONSTANT_ALPHA,
     //32772: ONE_MINUS_CONSTANT_COLOR
   };
-
   var WEBGL_TYPE_SIZES = {
     'SCALAR': 1,
     'VEC2': 2,
@@ -1116,6 +1114,7 @@ threeExamples.LegacyGLTFLoader = function () {
 
         // According to COLLADA spec...
         // aspectRatio = xfov / yfov
+
         var _camera = new THREE.PerspectiveCamera(THREE.MathUtils.radToDeg(yfov * aspectRatio), aspectRatio, camera.perspective.znear || 1, camera.perspective.zfar || 2e6);
         if (camera.name !== undefined) _camera.name = camera.name;
         if (camera.extras) _camera.userData = camera.extras;

package/lib/Process/FeatureProcessing.js

@@ -17,9 +17,8 @@ export default {
     if (node.layerUpdateState[layer.id] === undefined) {
       node.layerUpdateState[layer.id] = new LayerUpdateState();
     } else if (!node.layerUpdateState[layer.id].canTryUpdate()) {
-      var _node$link$layer$id;
       // toggle visibility features
-      (_node$link$layer$id = node.link[layer.id]) === null || _node$link$layer$id === void 0 ? void 0 : _node$link$layer$id.forEach(f => {
+      node.link[layer.id]?.forEach(f => {
         f.layer.object3d.add(f);
         f.meshes.position.z = geoidLayerIsVisible(layer.parent) ? node.geoidHeight : 0;
         f.meshes.updateMatrixWorld();

package/lib/Process/LayeredMaterialNodeProcessing.js

@@ -1,8 +1,6 @@
 import { chooseNextLevelToFetch } from "../Layer/LayerUpdateStrategy.js";
 import LayerUpdateState from "../Layer/LayerUpdateState.js";
 import handlingError from "./handlerNodeError.js";
-export const SIZE_TEXTURE_TILE = 256;
-export const SIZE_DIAGONAL_TEXTURE = (2 * (SIZE_TEXTURE_TILE * SIZE_TEXTURE_TILE)) ** 0.5;
 function materialCommandQueuePriorityFunction(material) {
   // We know that 'node' is visible because commands can only be
   // issued for visible nodes.
@@ -62,14 +60,12 @@ export function updateLayeredMaterialNodeImagery(context, layer, node, parent) {
         return;
       } // ok, we're going to inherit our parent's texture
     }
-
     if (!nodeLayer) {
-      var _parent$material;
       // Create new raster node
       nodeLayer = layer.setupRasterNode(node);
 
       // Init the node by parent
-      const parentLayer = (_parent$material = parent.material) === null || _parent$material === void 0 ? void 0 : _parent$material.getLayer(layer.id);
+      const parentLayer = parent.material?.getLayer(layer.id);
       nodeLayer.initFromParent(parentLayer, extentsDestination);
     }
 
@@ -158,9 +154,8 @@ export function updateLayeredMaterialNodeElevation(context, layer, node, parent)
     nodeLayer = layer.setupRasterNode(node);
   }
   if (node.layerUpdateState[layer.id] === undefined) {
-    var _parent$material2;
     node.layerUpdateState[layer.id] = new LayerUpdateState();
-    const parentLayer = (_parent$material2 = parent.material) === null || _parent$material2 === void 0 ? void 0 : _parent$material2.getLayer(layer.id);
+    const parentLayer = parent.material?.getLayer(layer.id);
     nodeLayer.initFromParent(parentLayer, extentsDestination);
     if (nodeLayer.level >= layer.source.zoom.min) {
       context.view.notifyChange(node, false);
@@ -210,8 +205,7 @@ export function removeLayeredMaterialNodeLayer(layerId) {
    * @param {TileMesh} node - The node to udpate.
    */
   return function (node) {
-    var _node$material;
-    if ((_node$material = node.material) !== null && _node$material !== void 0 && _node$material.removeLayer) {
+    if (node.material?.removeLayer) {
       if (node.material.elevationLayerIds.indexOf(layerId) > -1) {
         node.setBBoxZ({
           min: 0,

package/lib/Process/ObjectRemovalHelper.js

@@ -28,7 +28,6 @@ export default {
         // see https://github.com/iTowns/itowns/issues/869
         // obj.geometry = null;
       }
-
       if (obj.material) {
         if (Array.isArray(obj.material)) {
           for (const material of obj.material) {
@@ -82,7 +81,7 @@ export default {
     // of the objects which have their own removal logic
     let toRemove = obj.children.filter(c => c.layer && c.layer.id === layer.id);
     const linked = obj.link && obj.link[layer.id];
-    if (linked !== null && linked !== void 0 && linked.children.length) {
+    if (linked?.children.length) {
       toRemove = toRemove.concat(linked.children);
       delete obj.link[layer.id];
     }

package/lib/Provider/3dTilesProvider.js

@@ -114,6 +114,7 @@ function executeCommand(command) {
         } else if (magic == 'b3dm') {
           func = supportedFormats.b3dm;
         } else if (magic == 'pnts') {
+          layer.hasPnts = true;
           func = supportedFormats.pnts;
         } else if (magic == 'glTF') {
           func = supportedFormats.gltf;

package/lib/Renderer/ColorLayersOrdering.js

@@ -2,8 +2,7 @@ import { ImageryLayers } from "../Layer/Layer.js";
 function updateLayersOrdering(geometryLayer, imageryLayers) {
   const sequence = ImageryLayers.getColorLayersIdOrderedBySequence(imageryLayers);
   const cO = function (object) {
-    var _object$material;
-    if ((_object$material = object.material) !== null && _object$material !== void 0 && _object$material.setSequence) {
+    if (object.material?.setSequence) {
       object.material.setSequence(sequence);
     }
   };

package/lib/Renderer/Label2DRenderer.js

@@ -157,10 +157,7 @@ class Label2DRenderer {
       }
     });
     labelLayers.forEach(labelLayer => {
-      labelLayer.toHide.children.forEach(labelsNode => {
-        var _labelsNode$domElemen;
-        return (_labelsNode$domElemen = labelsNode.domElements) === null || _labelsNode$domElemen === void 0 ? void 0 : _labelsNode$domElemen.labels.hide();
-      });
+      labelLayer.toHide.children.forEach(labelsNode => labelsNode.domElements?.labels.hide());
       labelLayer.toHide.clear();
     });
   }

package/lib/Renderer/PointsMaterial.js

@@ -1,6 +1,6 @@
 import * as THREE from 'three';
 /* babel-plugin-inline-import './Shader/PointsVS.glsl' */
-const PointsVS = "#include <itowns/precision_qualifier>\n#if defined(USE_TEXTURES_PROJECTIVE)\n#include <itowns/projective_texturing_pars_vertex>\n#endif\n#include <common>\n#include <logdepthbuf_pars_vertex>\n\n#define NB_CLASS 8.\n\nuniform float size;\nuniform float scale;\n\nuniform bool picking;\nuniform int mode;\nuniform float opacity;\nuniform vec4 overlayColor;\n\nuniform vec2 elevationRange;\nuniform vec2 intensityRange;\nuniform vec2 angleRange;\n\nuniform bool applyOpacityClassication;\n\nuniform sampler2D classificationTexture;\nuniform sampler2D discreteTexture;\nuniform sampler2D gradientTexture;\nuniform int sizeMode;\nuniform float minAttenuatedSize;\nuniform float maxAttenuatedSize;\n\nattribute vec3 color;\nattribute vec2 range;\nattribute vec4 unique_id;\nattribute float intensity;\nattribute float classification;\nattribute float pointSourceID;\n\nattribute float returnNumber;\nattribute float numberOfReturns;\nattribute float scanAngle;\n\n#if defined(NORMAL_OCT16)\nattribute vec2 oct16Normal;\n#elif defined(NORMAL_SPHEREMAPPED)\nattribute vec2 sphereMappedNormal;\n#endif\n\nvarying vec4 vColor;\n\n// see https://web.archive.org/web/20150303053317/http://lgdv.cs.fau.de/get/1602\n// and implementation in PotreeConverter (BINPointReader.cpp) and potree (BinaryDecoderWorker.js)\n#if defined(NORMAL_OCT16)\nvec3 decodeOct16Normal(vec2 encodedNormal) {\n    vec2 nNorm = 2. * (encodedNormal / 255.) - 1.;\n    vec3 n;\n    n.z = 1. - abs(nNorm.x) - abs(nNorm.y);\n    if (n.z >= 0.) {\n        n.x = nNorm.x;\n        n.y = nNorm.y;\n    } else {\n        n.x = sign(nNorm.x) - sign(nNorm.x) * sign(nNorm.y) * nNorm.y;\n        n.y = sign(nNorm.y) - sign(nNorm.y) * sign(nNorm.x) * nNorm.x;\n    }\n    return normalize(n);\n}\n#elif defined(NORMAL_SPHEREMAPPED)\n// see http://aras-p.info/texts/CompactNormalStorage.html method #4\n// or see potree's implementation in BINPointReader.cpp\nvec3 decodeSphereMappedNormal(vec2 encodedNormal) {\n    vec2 fenc = 2. * encodedNormal / 255. - 1.;\n    float f = dot(fenc,fenc);\n    float g = 2. * sqrt(1. - f);\n    vec3 n;\n    n.xy = fenc * g;\n    n.z = 1. - 2. * f;\n    return n;\n}\n#endif\n\nvoid main() {\n\n#if defined(NORMAL_OCT16)\n    vec3  normal = decodeOct16Normal(oct16Normal);\n#elif defined(NORMAL_SPHEREMAPPED)\n    vec3 normal = decodeSphereMappedNormal(sphereMappedNormal);\n#elif defined(NORMAL)\n    // nothing to do\n#else\n    // default to color\n    vec3 normal = color;\n#endif\n\n    if (picking) {\n        vColor = unique_id;\n    } else {\n        vColor.a = opacity;\n        if (applyOpacityClassication || mode == PNTS_MODE_CLASSIFICATION) {\n            vec2 uv = vec2(classification/255., 0.5);\n            vColor = texture2D(classificationTexture, uv);\n            vColor.a *= opacity;\n        }\n\n        if (mode == PNTS_MODE_NORMAL) {\n            vColor.rgb = abs(normal);\n        } else if (mode == PNTS_MODE_COLOR) {\n            // default to color mode\n            vColor.rgb = mix(color, overlayColor.rgb, overlayColor.a);\n        } else if (mode == PNTS_MODE_RETURN_NUMBER) {\n            vec2 uv = vec2(returnNumber/255., 0.5);\n            vColor = texture2D(discreteTexture, uv);\n        } else if (mode == PNTS_MODE_RETURN_TYPE) {\n            float returnType;\n            if (returnNumber > numberOfReturns) {\n                returnType = 4.;\n            } else if (returnNumber == 1.) {\n                if (numberOfReturns == 1.) {\n                    // single\n                    returnType = 0.;\n                } else {\n                    // first\n                    returnType = 1.;\n                }\n            } else {\n                if (returnNumber == numberOfReturns) {\n                    // last\n                    returnType = 3.;\n                } else {\n                    // intermediate\n                    returnType = 2.;\n                }\n            }\n            vec2 uv = vec2(returnType/255., 0.5);\n            vColor = texture2D(discreteTexture, uv);\n        } else if (mode == PNTS_MODE_RETURN_COUNT) {\n            vec2 uv = vec2(numberOfReturns/255., 0.5);\n            vColor = texture2D(discreteTexture, uv);\n        } else if (mode == PNTS_MODE_POINT_SOURCE_ID) {\n            vec2 uv = vec2(mod(pointSourceID, NB_CLASS)/255., 0.5);\n            vColor = texture2D(discreteTexture, uv);\n        } else if (mode == PNTS_MODE_SCAN_ANGLE) {\n            float i = (scanAngle - angleRange.x) / (angleRange.y - angleRange.x);\n            vec2 uv = vec2(i, (1. - i));\n            vColor = texture2D(gradientTexture, uv);\n        } else if (mode == PNTS_MODE_INTENSITY) {\n            float i = (intensity - intensityRange.x) / (intensityRange.y - intensityRange.x);\n            vec2 uv = vec2(i, (1. - i));\n            vColor = texture2D(gradientTexture, uv);\n        } else if (mode == PNTS_MODE_ELEVATION) {\n            float i = (position.z - elevationRange.x) / (elevationRange.y - elevationRange.x);\n            vec2 uv = vec2(i, (1. - i));\n            vColor = texture2D(gradientTexture, uv);\n        }\n    }\n\n    #include <begin_vertex>\n    #include <project_vertex>\n\n    gl_PointSize = size;\n\n    if (sizeMode == PNTS_SIZE_MODE_ATTENUATED) {\n        bool isPerspective = isPerspectiveMatrix(projectionMatrix);\n\n        if (isPerspective) {\n            gl_PointSize *= scale / -mvPosition.z;\n            gl_PointSize = clamp(gl_PointSize, minAttenuatedSize, maxAttenuatedSize);\n        }\n    }\n\n#if defined(USE_TEXTURES_PROJECTIVE)\n    #include <itowns/projective_texturing_vertex>\n#endif\n    #include <logdepthbuf_vertex>\n}\n";
+const PointsVS = "#include <itowns/precision_qualifier>\n#if defined(USE_TEXTURES_PROJECTIVE)\n#include <itowns/projective_texturing_pars_vertex>\n#endif\n#include <common>\n#include <logdepthbuf_pars_vertex>\n\n#define NB_CLASS 8.\n\nuniform float size;\nuniform float scale;\n\nuniform bool picking;\nuniform int mode;\nuniform float opacity;\nuniform vec4 overlayColor;\n\nuniform vec2 elevationRange;\nuniform vec2 intensityRange;\nuniform vec2 angleRange;\n\nuniform sampler2D classificationTexture;\nuniform sampler2D discreteTexture;\nuniform sampler2D gradientTexture;\nuniform int sizeMode;\nuniform float minAttenuatedSize;\nuniform float maxAttenuatedSize;\n\nattribute vec3 color;\nattribute vec2 range;\nattribute vec4 unique_id;\nattribute float intensity;\nattribute float classification;\nattribute float pointSourceID;\n\nattribute float returnNumber;\nattribute float numberOfReturns;\nattribute float scanAngle;\n\n#if defined(NORMAL_OCT16)\nattribute vec2 oct16Normal;\n#elif defined(NORMAL_SPHEREMAPPED)\nattribute vec2 sphereMappedNormal;\n#endif\n\nvarying vec4 vColor;\n\n// see https://web.archive.org/web/20150303053317/http://lgdv.cs.fau.de/get/1602\n// and implementation in PotreeConverter (BINPointReader.cpp) and potree (BinaryDecoderWorker.js)\n#if defined(NORMAL_OCT16)\nvec3 decodeOct16Normal(vec2 encodedNormal) {\n    vec2 nNorm = 2. * (encodedNormal / 255.) - 1.;\n    vec3 n;\n    n.z = 1. - abs(nNorm.x) - abs(nNorm.y);\n    if (n.z >= 0.) {\n        n.x = nNorm.x;\n        n.y = nNorm.y;\n    } else {\n        n.x = sign(nNorm.x) - sign(nNorm.x) * sign(nNorm.y) * nNorm.y;\n        n.y = sign(nNorm.y) - sign(nNorm.y) * sign(nNorm.x) * nNorm.x;\n    }\n    return normalize(n);\n}\n#elif defined(NORMAL_SPHEREMAPPED)\n// see http://aras-p.info/texts/CompactNormalStorage.html method #4\n// or see potree's implementation in BINPointReader.cpp\nvec3 decodeSphereMappedNormal(vec2 encodedNormal) {\n    vec2 fenc = 2. * encodedNormal / 255. - 1.;\n    float f = dot(fenc,fenc);\n    float g = 2. * sqrt(1. - f);\n    vec3 n;\n    n.xy = fenc * g;\n    n.z = 1. - 2. * f;\n    return n;\n}\n#endif\n\nvoid main() {\n\n#if defined(NORMAL_OCT16)\n    vec3  normal = decodeOct16Normal(oct16Normal);\n#elif defined(NORMAL_SPHEREMAPPED)\n    vec3 normal = decodeSphereMappedNormal(sphereMappedNormal);\n#elif defined(NORMAL)\n    // nothing to do\n#else\n    // default to color\n    vec3 normal = color;\n#endif\n\n    if (picking) {\n        vColor = unique_id;\n    } else {\n        vColor.a = 1.0;\n        if (mode == PNTS_MODE_CLASSIFICATION) {\n            vec2 uv = vec2(classification/255., 0.5);\n            vColor = texture2D(classificationTexture, uv);\n        } else if (mode == PNTS_MODE_NORMAL) {\n            vColor.rgb = abs(normal);\n        } else if (mode == PNTS_MODE_COLOR) {\n            // default to color mode\n            vColor.rgb = mix(color, overlayColor.rgb, overlayColor.a);\n        } else if (mode == PNTS_MODE_RETURN_NUMBER) {\n            vec2 uv = vec2(returnNumber/255., 0.5);\n            vColor = texture2D(discreteTexture, uv);\n        } else if (mode == PNTS_MODE_RETURN_TYPE) {\n            float returnType;\n            if (returnNumber > numberOfReturns) {\n                returnType = 4.;\n            } else if (returnNumber == 1.) {\n                if (numberOfReturns == 1.) {\n                    // single\n                    returnType = 0.;\n                } else {\n                    // first\n                    returnType = 1.;\n                }\n            } else {\n                if (returnNumber == numberOfReturns) {\n                    // last\n                    returnType = 3.;\n                } else {\n                    // intermediate\n                    returnType = 2.;\n                }\n            }\n            vec2 uv = vec2(returnType/255., 0.5);\n            vColor = texture2D(discreteTexture, uv);\n        } else if (mode == PNTS_MODE_RETURN_COUNT) {\n            vec2 uv = vec2(numberOfReturns/255., 0.5);\n            vColor = texture2D(discreteTexture, uv);\n        } else if (mode == PNTS_MODE_POINT_SOURCE_ID) {\n            vec2 uv = vec2(mod(pointSourceID, NB_CLASS)/255., 0.5);\n            vColor = texture2D(discreteTexture, uv);\n        } else if (mode == PNTS_MODE_SCAN_ANGLE) {\n            float i = (scanAngle - angleRange.x) / (angleRange.y - angleRange.x);\n            vec2 uv = vec2(i, (1. - i));\n            vColor = texture2D(gradientTexture, uv);\n        } else if (mode == PNTS_MODE_INTENSITY) {\n            float i = (intensity - intensityRange.x) / (intensityRange.y - intensityRange.x);\n            vec2 uv = vec2(i, (1. - i));\n            vColor = texture2D(gradientTexture, uv);\n        } else if (mode == PNTS_MODE_ELEVATION) {\n            float i = (position.z - elevationRange.x) / (elevationRange.y - elevationRange.x);\n            vec2 uv = vec2(i, (1. - i));\n            vColor = texture2D(gradientTexture, uv);\n        }\n\n        vColor.a *= opacity;\n    }\n\n    #include <begin_vertex>\n    #include <project_vertex>\n\n    gl_PointSize = size;\n\n    if (sizeMode == PNTS_SIZE_MODE_ATTENUATED) {\n        bool isPerspective = isPerspectiveMatrix(projectionMatrix);\n\n        if (isPerspective) {\n            gl_PointSize *= scale / -mvPosition.z;\n            gl_PointSize = clamp(gl_PointSize, minAttenuatedSize, maxAttenuatedSize);\n        }\n    }\n\n#if defined(USE_TEXTURES_PROJECTIVE)\n    #include <itowns/projective_texturing_vertex>\n#endif\n    #include <logdepthbuf_vertex>\n}\n";
 /* babel-plugin-inline-import './Shader/PointsFS.glsl' */
 const PointsFS = "#include <itowns/precision_qualifier>\n#include <logdepthbuf_pars_fragment>\n#if defined(USE_TEXTURES_PROJECTIVE)\n#include <itowns/projective_texturing_pars_fragment>\n#endif\n\nvarying vec4 vColor;\nuniform bool picking;\nuniform int shape;\n\nvoid main() {\n    #include <logdepthbuf_fragment>\n    //square shape does not require any change.\n    if (shape == PNTS_SHAPE_CIRCLE) {\n        //circular rendering in glsl\n        if ((length(gl_PointCoord - 0.5) > 0.5) || (vColor.a == 0.0)) {\n            discard;\n        }\n    }\n\n#if defined(USE_TEXTURES_PROJECTIVE)\n    vec4 color = vColor;\n    if (!picking) {\n        #pragma unroll_loop\n        for (int i = 0; i < ORIENTED_IMAGES_COUNT; i++) {\n            color = projectiveTextureColor(projectiveTextureCoords[ ORIENTED_IMAGES_COUNT - 1 - i ], projectiveTextureDistortion[ ORIENTED_IMAGES_COUNT - 1 - i ], projectiveTexture[ ORIENTED_IMAGES_COUNT - 1 - i ], mask[ORIENTED_IMAGES_COUNT - 1 - i], color);\n        }\n        gl_FragColor = vec4(color.rgb, color.a * opacity);\n    } else {\n        gl_FragColor = color;\n    }\n#else\n    gl_FragColor = vColor;\n#endif\n}\n";
 import ShaderUtils from "./Shader/ShaderUtils.js";
@@ -127,7 +127,7 @@ export const ClassificationScheme = {
       visible: true,
       name: 'default',
       color: new THREE.Color(0.3, 0.6, 0.6),
-      opacity: 0.5
+      opacity: 1.0
     }
   }
 };
@@ -185,7 +185,7 @@ const DiscreteScheme = {
       visible: true,
       name: 'default',
       color: white,
-      opacity: 0.5
+      opacity: 1.0
     }
   }
 };
@@ -220,6 +220,7 @@ function generateGradientTexture(gradient) {
   return texture;
 }
 function recomputeTexture(scheme, texture, nbClass) {
+  let needTransparency;
   const data = texture.image.data;
   const width = texture.image.width;
   if (!nbClass) {
@@ -250,8 +251,10 @@ function recomputeTexture(scheme, texture, nbClass) {
     data[j + 1] = parseInt(255 * color.g, 10);
     data[j + 2] = parseInt(255 * color.b, 10);
     data[j + 3] = visible ? parseInt(255 * opacity, 10) : 0;
+    needTransparency = needTransparency || opacity < 1;
   }
   texture.needsUpdate = true;
+  return needTransparency;
 }
 class PointsMaterial extends THREE.ShaderMaterial {
   /**
@@ -264,7 +267,6 @@ class PointsMaterial extends THREE.ShaderMaterial {
    * @param      {THREE.Vector2}  [options.intensityRange=new THREE.Vector2(1, 65536)]  intensity range.
    * @param      {THREE.Vector2}  [options.elevationRange=new THREE.Vector2(0, 1000)]  elevation range.
    * @param      {THREE.Vector2}  [options.angleRange=new THREE.Vector2(-90, 90)]  scan angle range.
-   * @param      {boolean}  [options.applyOpacityClassication=false]  apply opacity classification on all display mode.
    * @param      {Scheme}  [options.classification]  LUT for point classification colorization.
    * @param      {Scheme}  [options.discreteScheme]  LUT for other discret point values colorization.
    * @param      {string}  [options.gradient]  Descrition of the gradient to use for continuous point values.
@@ -293,7 +295,6 @@ class PointsMaterial extends THREE.ShaderMaterial {
     const oiMaterial = options.orientedImageMaterial;
     const classificationScheme = options.classification || ClassificationScheme.DEFAULT;
     const discreteScheme = options.discreteScheme || DiscreteScheme.DEFAULT;
-    const applyOpacityClassication = options.applyOpacityClassication == undefined ? false : options.applyOpacityClassication;
     const size = options.size || 0;
     const mode = options.mode || PNTS_MODE.COLOR;
     const shape = options.shape || PNTS_SHAPE.CIRCLE;
@@ -313,7 +314,6 @@ class PointsMaterial extends THREE.ShaderMaterial {
     delete options.orientedImageMaterial;
     delete options.classification;
     delete options.discreteScheme;
-    delete options.applyOpacityClassication;
     delete options.size;
     delete options.mode;
     delete options.shape;
@@ -322,6 +322,7 @@ class PointsMaterial extends THREE.ShaderMaterial {
     delete options.maxAttenuatedSize;
     delete options.gradient;
     super(options);
+    this.userData.needTransparency = {};
     this.gradients = gradients;
     this.gradientTexture = new THREE.CanvasTexture();
     this.vertexShader = PointsVS;
@@ -339,7 +340,6 @@ class PointsMaterial extends THREE.ShaderMaterial {
     CommonMaterial.setUniformProperty(this, 'intensityRange', intensityRange);
     CommonMaterial.setUniformProperty(this, 'elevationRange', elevationRange);
     CommonMaterial.setUniformProperty(this, 'angleRange', angleRange);
-    CommonMaterial.setUniformProperty(this, 'applyOpacityClassication', applyOpacityClassication);
     CommonMaterial.setUniformProperty(this, 'sizeMode', sizeMode);
     CommonMaterial.setUniformProperty(this, 'scale', scale);
     CommonMaterial.setUniformProperty(this, 'minAttenuatedSize', minAttenuatedSize);
@@ -390,14 +390,19 @@ class PointsMaterial extends THREE.ShaderMaterial {
     }
   }
   recomputeClassification() {
-    recomputeTexture(this.classificationScheme, this.classificationTexture, 32);
+    const needTransparency = recomputeTexture(this.classificationScheme, this.classificationTexture, 256);
+    this.userData.needTransparency[PNTS_MODE.CLASSIFICATION] = needTransparency;
     this.dispatchEvent({
       type: 'material_property_changed',
       target: this.uniforms
     });
   }
   recomputeDiscreteTexture() {
-    recomputeTexture(this.discreteScheme, this.discreteTexture);
+    const needTransparency = recomputeTexture(this.discreteScheme, this.discreteTexture);
+    this.userData.needTransparency[PNTS_MODE.RETURN_NUMBER] = needTransparency;
+    this.userData.needTransparency[PNTS_MODE.RETURN_TYPE] = needTransparency;
+    this.userData.needTransparency[PNTS_MODE.RETURN_COUNT] = needTransparency;
+    this.userData.needTransparency[PNTS_MODE.POINT_SOURCE_ID] = needTransparency;
     this.dispatchEvent({
       type: 'material_property_changed',
       target: this.uniforms

package/lib/Renderer/RenderMode.js

@@ -5,7 +5,6 @@ const MODES = {
   // depth buffer
   ID: 2 // id object
 };
-
 function push(object3d, mode) {
   const _mode = object3d.mode ?? MODES.FINAL;
   if (_mode == mode) {

package/lib/Source/CopcSource.js

@@ -0,0 +1,118 @@
+import { Binary, Info, Las } from 'copc';
+import Extent from "../Core/Geographic/Extent.js";
+import Fetcher from "../Provider/Fetcher.js";
+import LASParser from "../Parser/LASParser.js";
+import Source from "./Source.js";
+import * as THREE from 'three';
+
+/**
+ * @param {function(number, number):Promise<Uint8Array>} fetcher
+ */
+async function getHeaders(fetcher) {
+  const header = Las.Header.parse(await fetcher(0, Las.Constants.minHeaderLength));
+  const vlrs = await Las.Vlr.walk(fetcher, header);
+
+  // info VLR: required by COPC
+  const infoVlr = Las.Vlr.find(vlrs, 'copc', 1);
+  if (!infoVlr) {
+    return Promise.reject('COPC info VLR is required');
+  }
+  const info = Info.parse(await Las.Vlr.fetch(fetcher, infoVlr));
+
+  // OGC Coordinate System WKT: required by LAS1.4
+  const wktVlr = Las.Vlr.find(vlrs, 'LASF_Projection', 2112);
+  if (!wktVlr) {
+    return Promise.reject('LAS1.4 WKT VLR is required');
+  }
+  const wkt = Binary.toCString(await Las.Vlr.fetch(fetcher, wktVlr));
+
+  // Extra bytes: optional by LAS1.4
+  const ebVlr = Las.Vlr.find(vlrs, 'LASF_Spec', 4);
+  const eb = ebVlr ? Las.ExtraBytes.parse(await Las.Vlr.fetch(fetcher, ebVlr)) : [];
+  return {
+    header,
+    info,
+    wkt,
+    eb
+  };
+}
+
+/**
+ * @classdesc
+ * A source for [Cloud Optimised Point Cloud](https://copc.io/) (COPC) data.
+ * Such data consists of a [LAZ 1.4](https://www.ogc.org/standard/las/) file
+ * that stores compressed points data organized in a clustered octree.
+ *
+ * A freshly created source fetches and parses portions of the file
+ * corresponding to the LAS 1.4 header, all the Variable Length Record (VLR)
+ * headers as well the following VLRs:
+ * - COPC [`info`](https://copc.io/#info-vlr) record (mandatory)
+ * - LAS 1.4 `OGC Coordinate System WKT` record (mandatory, see [Las 1.4
+ *   spec](https://portal.ogc.org/files/?artifact_id=74523))
+ * - LAS 1.4 `Extra Bytes` record (optional, see [Las 1.4
+ *   spec](https://portal.ogc.org/files/?artifact_id=74523))
+ *
+ * @extends {Source}
+ *
+ * @property {boolean} isCopcSource - Read-only flag to check that a given
+ * object is of type CopcSource.
+ * @property {Object} header - LAS header of the source.
+ * @property {Object[]} eb - List of headers of each Variable Length Records
+ * (VLRs).
+ * @property {Object} info - COPC `info` VLR.
+ * @property {number[]} info.cube - Bounding box of the octree as a 6-elements.
+ * tuple `[minX, minY, minZ, maxX, maxY, maxZ]`. Computed from `center_x`,
+ * `center_y`, `center_z` and `halfSize` properties.
+ * @property {Object} info.rootHierarchyPage - Hierarchy page of the root node.
+ * @property {number} info.rootHierarchyPage.pageOffset - Absolute Offset to the
+ * root node data chunk.
+ * @property {number} info.rootHierarchyPage.pageOffset - Size (in bytes) of the
+ * root node data chunk.
+ * @property {number[]} gpsTimeRange - A 2-element tuple denoting the minimum
+ * and maximum values of attribute `gpsTime`.
+ */
+class CopcSource extends Source {
+  /**
+   * @param {Object} config - Source configuration
+   * @param {string} config.url - URL of the COPC resource.
+   * @param {8 | 16} [config.colorDepth=16] - Encoding of the `color`
+   * attribute. Either `8` or `16` bits.
+   * @param {string} [config._lazPerfBaseUrl] - (experimental) Overrides base
+   * url of the `las-zip.wasm` file of the `laz-perf` library.
+   * @param {string} [config.crs='EPSG:4326'] - Native CRS of the COPC
+   * ressource. Note that this is not for now inferred from the COPC header.
+   * @param {RequestInit} [config.networkOptions] - Fetch options (passed
+   * directly to `fetch()`), see [the syntax for more information]{@link
+   * https://developer.mozilla.org/en-US/docs/Web/API/WindowOrWorkerGlobalScope/fetch#Syntax}.
+   * @param {Object} [config.attribution] - Attribution of the data.
+   *
+   * @constructor
+   */
+  constructor(config) {
+    super(config);
+    this.isCopcSource = true;
+    this.parser = LASParser.parseChunk;
+    this.fetcher = Fetcher.arrayBuffer;
+    this.colorDepth = config.colorDepth ?? 16;
+    const get = ( /** @type {number} */begin, /** @type {number} */end) => this.fetcher(this.url, {
+      ...this.networkOptions,
+      headers: {
+        ...this.networkOptions.headers,
+        range: `bytes=${begin}-${end - 1}`
+      }
+    }).then(buffer => new Uint8Array(buffer));
+    this.whenReady = getHeaders(get).then(metadata => {
+      this.header = metadata.header;
+      this.info = metadata.info;
+      this.eb = metadata.eb;
+      // TODO: use wkt definition in `metadata.wkt` to infer/define crs
+      this.crs = config.crs || 'EPSG:4326';
+      const bbox = new THREE.Box3();
+      bbox.min.fromArray(this.info.cube, 0);
+      bbox.max.fromArray(this.info.cube, 3);
+      this.extent = Extent.fromBox3(this.crs, bbox);
+      return this;
+    });
+  }
+}
+export default CopcSource;

package/lib/Source/Source.js

@@ -56,6 +56,8 @@ let uid = 0;
  * To extend a Source, it is necessary to implement two functions:
  * `urlFromExtent` and `extentInsideLimit`.
  *
+ * @extends InformationsData
+ *
  * @property {boolean} isSource - Used to checkout whether this source is a
  * Source. Default is true. You should not change this, as it is used internally
  * for optimisation.
@@ -97,7 +99,6 @@ class Source extends InformationsData {
    * Source. Only the `url` property is mandatory.
    *
    * @constructor
-   * @extends InformationsData
    */
   constructor(source) {
     super(source);
@@ -118,6 +119,7 @@ class Source extends InformationsData {
       crossOrigin: 'anonymous'
     };
     this.attribution = source.attribution;
+    /** @type {Promise<any>} */
     this.whenReady = Promise.resolve();
     this._featuresCaches = {};
     if (source.extent && !source.extent.isExtent) {

package/lib/ThreeExtended/loaders/DDSLoader.js

@@ -17,23 +17,29 @@ class DDSLoader extends CompressedTextureLoader {
 
     // All values and structures referenced from:
     // http://msdn.microsoft.com/en-us/library/bb943991.aspx/
+
     // const DDSD_CAPS = 0x1;
     // const DDSD_HEIGHT = 0x2;
     // const DDSD_WIDTH = 0x4;
     // const DDSD_PITCH = 0x8;
     // const DDSD_PIXELFORMAT = 0x1000;
+
     // const DDSD_LINEARSIZE = 0x80000;
     // const DDSD_DEPTH = 0x800000;
+
     // const DDSCAPS_COMPLEX = 0x8;
     // const DDSCAPS_MIPMAP = 0x400000;
     // const DDSCAPS_TEXTURE = 0x1000;
+
     // const DDSCAPS2_VOLUME = 0x200000;
+
     // const DDPF_ALPHAPIXELS = 0x1;
     // const DDPF_ALPHA = 0x2;
     // const DDPF_FOURCC = 0x4;
     // const DDPF_RGB = 0x40;
     // const DDPF_YUV = 0x200;
     // const DDPF_LUMINANCE = 0x20000;
+
     function fourCCToInt32(value) {
       return value.charCodeAt(0) + (value.charCodeAt(1) << 8) + (value.charCodeAt(2) << 16) + (value.charCodeAt(3) << 24);
     }
@@ -66,7 +72,6 @@ class DDSLoader extends CompressedTextureLoader {
           dst++; //a
         }
       }
-
       return byteArray;
     }
     const FOURCC_DXT1 = fourCCToInt32('DXT1');
@@ -78,11 +83,16 @@ class DDSLoader extends CompressedTextureLoader {
     const extendedHeaderLengthInt = 5; // The extended header length in 32 bit ints
 
     // Offsets into the header array
+
     // const off_pfFlags = 20;
+
     // const off_caps = 27;
+
     // const off_caps3 = 29;
     // const off_caps4 = 30;
+
     // If fourCC = DX10, the extended header starts after 32
+
     // Parse header
 
     const header = new Int32Array(buffer, 0, headerLengthInt);

package/lib/ThreeExtended/loaders/DRACOLoader.js

@@ -261,7 +261,6 @@ function DRACOWorker() {
           };
           DracoDecoderModule(decoderConfig); // eslint-disable-line no-undef
         });
-
         break;
       case 'decode':
         const buffer = message.buffer;

package/lib/ThreeExtended/loaders/GLTFLoader.js

@@ -2627,6 +2627,7 @@ class GLTFParser {
 
       // Removes dangling associations, associations that reference a node that
       // didn't make it into the scene.
+
       parser.associations = (node => {
         const reducedAssociations = new Map();
         for (const [key, value] of parser.associations) {