forked from iTowns/itowns
-
Notifications
You must be signed in to change notification settings - Fork 0
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
feat(LASParser): change lasparser package from loaders.gl to copc
- Loading branch information
1 parent
f89df8c
commit aa9d97e
Showing
5 changed files
with
257 additions
and
66 deletions.
There are no files selected for viewing
Binary file not shown.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,151 @@ | ||
| import { LazPerf } from 'laz-perf'; | ||
| import { Las } from 'copc'; | ||
|
|
||
| /** | ||
| * @typedef {Object} Header - Partial LAS header. | ||
| * @property {number} header.pointDataRecordFormat - Type of point data | ||
| * records contained by the buffer. | ||
| * @property {number} header.pointDataRecordLength - Size (in bytes) of the | ||
| * point data records. If the specified size is larger than implied by the | ||
| * point data record format (see above) the remaining bytes are user-specfic | ||
| * "extra bytes". Those are described by an Extra Bytes VLR. | ||
| * @property {number[]} header.scale - Scale factors (an array `[xScale, | ||
| * yScale, zScale]`) multiplied to the X, Y, Z point record values. | ||
| * @property {number[]} header.offset - Offsets (an array `[xOffset, | ||
| * xOffset, zOffset]`) added to the scaled X, Y, Z point record values. | ||
| */ | ||
|
|
||
| /** | ||
| * @classdesc | ||
| * Loader for LAS and LAZ (LASZip) point clouds. It uses the copc.js library and | ||
| * the laz-perf decoder under the hood. | ||
| * | ||
| * The laz-perf web assembly module is lazily fetched at runtime when a parsing | ||
| * request is initiated. Location of laz-perf wasm defaults to the unpkg | ||
| * repository. | ||
| */ | ||
| class LASLoader { | ||
| constructor() { | ||
| this._wasmPath = 'https://unpkg.com/[email protected]/lib/'; | ||
| this._wasmPromise = null; | ||
| } | ||
|
|
||
| _initDecoder() { | ||
| if (this._wasmPromise) { | ||
| return this._wasmPromise; | ||
| } | ||
|
|
||
| this._wasmPromise = LazPerf.create({ | ||
| locateFile: file => `${this._wasmPath}/${file}`, | ||
| }); | ||
|
|
||
| return this._wasmPromise; | ||
| } | ||
|
|
||
| _parseView(view, options) { | ||
| const colorDepth = options.colorDepth ?? 16; | ||
|
|
||
| const getPosition = ['X', 'Y', 'Z'].map(view.getter); | ||
| const getIntensity = view.getter('Intensity'); | ||
| const getReturnNumber = view.getter('ReturnNumber'); | ||
| const getNumberOfReturns = view.getter('NumberOfReturns'); | ||
| const getClassification = view.getter('Classification'); | ||
| const getPointSourceID = view.getter('PointSourceId'); | ||
| const getColor = view.dimensions.Red ? | ||
| ['Red', 'Green', 'Blue'].map(view.getter) : undefined; | ||
|
|
||
| const positions = new Float32Array(view.pointCount * 3); | ||
| const intensities = new Uint16Array(view.pointCount); | ||
| const returnNumbers = new Uint8Array(view.pointCount); | ||
| const numberOfReturns = new Uint8Array(view.pointCount); | ||
| const classifications = new Uint8Array(view.pointCount); | ||
| const pointSourceIDs = new Uint16Array(view.pointCount); | ||
| const colors = getColor ? new Uint8Array(view.pointCount * 4) : undefined; | ||
|
|
||
| for (let i = 0; i < view.pointCount; i++) { | ||
| // `getPosition` apply scale and offset transform to the X, Y, Z | ||
| // values. See https://github.com/connormanning/copc.js/blob/master/src/las/extractor.ts. | ||
| const [x, y, z] = getPosition.map(f => f(i)); | ||
| positions[i * 3] = x; | ||
| positions[i * 3 + 1] = y; | ||
| positions[i * 3 + 2] = z; | ||
|
|
||
| intensities[i] = getIntensity(i); | ||
| returnNumbers[i] = getReturnNumber(i); | ||
| numberOfReturns[i] = getNumberOfReturns(i); | ||
|
|
||
| if (getColor) { | ||
| // Note that we do not infer color depth as it is expensive | ||
| // (i.e. traverse the whole view to check if there exists a red, | ||
| // green or blue value > 255). | ||
| let [r, g, b] = getColor.map(f => f(i)); | ||
|
|
||
| if (colorDepth === 16) { | ||
| r /= 256; | ||
| g /= 256; | ||
| b /= 256; | ||
| } | ||
|
|
||
| colors[i * 4] = r; | ||
| colors[i * 4 + 1] = g; | ||
| colors[i * 4 + 2] = b; | ||
| colors[i * 4 + 3] = 255; | ||
| } | ||
|
|
||
| classifications[i] = getClassification(i); | ||
| pointSourceIDs[i] = getPointSourceID(i); | ||
| } | ||
|
|
||
| return { | ||
| position: positions, | ||
| intensity: intensities, | ||
| returnNumber: returnNumbers, | ||
| numberOfReturns, | ||
| classification: classifications, | ||
| pointSourceID: pointSourceIDs, | ||
| color: colors, | ||
| }; | ||
| } | ||
|
|
||
| /** | ||
| * Set LazPerf decoder path. | ||
| * @param {string} path - path to `laz-perf.wasm` folder. | ||
| */ | ||
| set lazPerf(path) { | ||
| this._wasmPath = path; | ||
| this._wasmPromise = null; | ||
| } | ||
|
|
||
| /** | ||
| * Parses a LAS or LAZ (LASZip) file. Note that this function is | ||
| * **CPU-bound** and shall be parallelised in a dedicated worker. | ||
| * @param {ArrayBuffer} data - Binary data to parse. | ||
| * @param {Object} [options] - Parsing options. | ||
| * @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 parseFile(data, options = {}) { | ||
| 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 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); | ||
| const eb = ebVlr && Las.ExtraBytes.parse(await Las.Vlr.fetch(getter, ebVlr)); | ||
|
|
||
| const view = Las.View.create(pointData, header, eb); | ||
| const attributes = this._parseView(view, { colorDepth }); | ||
| return { | ||
| header, | ||
| attributes, | ||
| }; | ||
| } | ||
| } | ||
|
|
||
| export default LASLoader; |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Oops, something went wrong.