UNPKG

@kibeo/loaders.gl-geopackage

Version:

GeoPackage data loaders

395 lines (320 loc) 12.3 kB
"use strict"; var _interopRequireDefault = require("@babel/runtime/helpers/interopRequireDefault"); Object.defineProperty(exports, "__esModule", { value: true }); exports.default = parseGeoPackage; var _regenerator = _interopRequireDefault(require("@babel/runtime/regenerator")); var _slicedToArray2 = _interopRequireDefault(require("@babel/runtime/helpers/slicedToArray")); var _asyncToGenerator2 = _interopRequireDefault(require("@babel/runtime/helpers/asyncToGenerator")); var _sql = _interopRequireDefault(require("sql.js")); var _loaders = require("@kibeo/loaders.gl-wkt"); var _loaders2 = require("@kibeo/loaders.gl-schema"); var _loaders3 = require("@kibeo/loaders.gl-gis"); var _proj = require("@math.gl/proj4"); function _createForOfIteratorHelper(o, allowArrayLike) { var it = typeof Symbol !== "undefined" && o[Symbol.iterator] || o["@@iterator"]; if (!it) { if (Array.isArray(o) || (it = _unsupportedIterableToArray(o)) || allowArrayLike && o && typeof o.length === "number") { if (it) o = it; var i = 0; var F = function F() {}; return { s: F, n: function n() { if (i >= o.length) return { done: true }; return { done: false, value: o[i++] }; }, e: function e(_e) { throw _e; }, f: F }; } throw new TypeError("Invalid attempt to iterate non-iterable instance.\nIn order to be iterable, non-array objects must have a [Symbol.iterator]() method."); } var normalCompletion = true, didErr = false, err; return { s: function s() { it = it.call(o); }, n: function n() { var step = it.next(); normalCompletion = step.done; return step; }, e: function e(_e2) { didErr = true; err = _e2; }, f: function f() { try { if (!normalCompletion && it.return != null) it.return(); } finally { if (didErr) throw err; } } }; } function _unsupportedIterableToArray(o, minLen) { if (!o) return; if (typeof o === "string") return _arrayLikeToArray(o, minLen); var n = Object.prototype.toString.call(o).slice(8, -1); if (n === "Object" && o.constructor) n = o.constructor.name; if (n === "Map" || n === "Set") return Array.from(o); if (n === "Arguments" || /^(?:Ui|I)nt(?:8|16|32)(?:Clamped)?Array$/.test(n)) return _arrayLikeToArray(o, minLen); } function _arrayLikeToArray(arr, len) { if (len == null || len > arr.length) len = arr.length; for (var i = 0, arr2 = new Array(len); i < len; i++) { arr2[i] = arr[i]; } return arr2; } var ENVELOPE_BYTE_LENGTHS = { 0: 0, 1: 32, 2: 48, 3: 48, 4: 64, 5: 0, 6: 0, 7: 0 }; var SQL_TYPE_MAPPING = { BOOLEAN: _loaders2.Bool, TINYINT: _loaders2.Int8, SMALLINT: _loaders2.Int16, MEDIUMINT: _loaders2.Int32, INT: _loaders2.Int32, INTEGER: _loaders2.Int32, FLOAT: _loaders2.Float32, DOUBLE: _loaders2.Float64, REAL: _loaders2.Float64, TEXT: _loaders2.Utf8, BLOB: _loaders2.Binary, DATE: _loaders2.Utf8, DATETIME: _loaders2.Utf8, GEOMETRY: _loaders2.Binary }; function parseGeoPackage(_x, _x2) { return _parseGeoPackage.apply(this, arguments); } function _parseGeoPackage() { _parseGeoPackage = (0, _asyncToGenerator2.default)(_regenerator.default.mark(function _callee(arrayBuffer, options) { var _ref2, _ref2$sqlJsCDN, sqlJsCDN, _ref3, _ref3$reproject, reproject, _ref3$_targetCrs, _targetCrs, db, tables, projections, result, _iterator2, _step2, table, tableName; return _regenerator.default.wrap(function _callee$(_context) { while (1) { switch (_context.prev = _context.next) { case 0: _ref2 = (options === null || options === void 0 ? void 0 : options.geopackage) || {}, _ref2$sqlJsCDN = _ref2.sqlJsCDN, sqlJsCDN = _ref2$sqlJsCDN === void 0 ? 'https://sql.js.org/dist/' : _ref2$sqlJsCDN; _ref3 = (options === null || options === void 0 ? void 0 : options.gis) || {}, _ref3$reproject = _ref3.reproject, reproject = _ref3$reproject === void 0 ? false : _ref3$reproject, _ref3$_targetCrs = _ref3._targetCrs, _targetCrs = _ref3$_targetCrs === void 0 ? 'WGS84' : _ref3$_targetCrs; _context.next = 4; return loadDatabase(arrayBuffer, sqlJsCDN); case 4: db = _context.sent; tables = listVectorTables(db); projections = getProjections(db); result = {}; _iterator2 = _createForOfIteratorHelper(tables); try { for (_iterator2.s(); !(_step2 = _iterator2.n()).done;) { table = _step2.value; tableName = table.table_name; result[tableName] = getVectorTable(db, tableName, projections, { reproject: reproject, _targetCrs: _targetCrs }); } } catch (err) { _iterator2.e(err); } finally { _iterator2.f(); } return _context.abrupt("return", result); case 11: case "end": return _context.stop(); } } }, _callee); })); return _parseGeoPackage.apply(this, arguments); } function loadDatabase(_x3, _x4) { return _loadDatabase.apply(this, arguments); } function _loadDatabase() { _loadDatabase = (0, _asyncToGenerator2.default)(_regenerator.default.mark(function _callee2(arrayBuffer, sqlJsCDN) { var SQL; return _regenerator.default.wrap(function _callee2$(_context2) { while (1) { switch (_context2.prev = _context2.next) { case 0: if (!sqlJsCDN) { _context2.next = 6; break; } _context2.next = 3; return (0, _sql.default)({ locateFile: function locateFile(file) { return "".concat(sqlJsCDN).concat(file); } }); case 3: SQL = _context2.sent; _context2.next = 9; break; case 6: _context2.next = 8; return (0, _sql.default)(); case 8: SQL = _context2.sent; case 9: return _context2.abrupt("return", new SQL.Database(new Uint8Array(arrayBuffer))); case 10: case "end": return _context2.stop(); } } }, _callee2); })); return _loadDatabase.apply(this, arguments); } function listVectorTables(db) { var stmt = db.prepare("SELECT * FROM gpkg_contents WHERE data_type='features';"); var vectorTablesInfo = []; while (stmt.step()) { var vectorTableInfo = stmt.getAsObject(); vectorTablesInfo.push(vectorTableInfo); } return vectorTablesInfo; } function getVectorTable(db, tableName, projections, _ref) { var reproject = _ref.reproject, _targetCrs = _ref._targetCrs; var dataColumns = getDataColumns(db, tableName); var geomColumn = getGeometryColumn(db, tableName); var featureIdColumn = getFeatureIdName(db, tableName); var _db$exec$ = db.exec("SELECT * FROM `".concat(tableName, "`;"))[0], columns = _db$exec$.columns, values = _db$exec$.values; var projection; if (reproject) { var geomColumnProjStr = projections[geomColumn.srs_id]; projection = new _proj.Proj4Projection({ from: geomColumnProjStr, to: _targetCrs }); } var geojsonFeatures = []; var _iterator = _createForOfIteratorHelper(values), _step; try { for (_iterator.s(); !(_step = _iterator.n()).done;) { var row = _step.value; var geojsonFeature = constructGeoJsonFeature(columns, row, geomColumn, dataColumns, featureIdColumn); geojsonFeatures.push(geojsonFeature); } } catch (err) { _iterator.e(err); } finally { _iterator.f(); } var schema = getArrowSchema(db, tableName); if (projection) { return { geojsonFeatures: (0, _loaders3.transformGeoJsonCoords)(geojsonFeatures, projection.project), schema: schema }; } return { geojsonFeatures: geojsonFeatures, schema: schema }; } function getProjections(db) { var stmt = db.prepare('SELECT * FROM gpkg_spatial_ref_sys;'); var projectionMapping = {}; while (stmt.step()) { var srsInfo = stmt.getAsObject(); var srs_id = srsInfo.srs_id, definition = srsInfo.definition; projectionMapping[srs_id] = definition; } return projectionMapping; } function constructGeoJsonFeature(columns, row, geomColumn, dataColumns, featureIdColumn) { var idIdx = columns.indexOf(featureIdColumn); var id = row[idIdx]; var geomColumnIdx = columns.indexOf(geomColumn.column_name); var geometry = parseGeometry(row[geomColumnIdx].buffer); var properties = {}; if (dataColumns) { for (var _i = 0, _Object$entries = Object.entries(dataColumns); _i < _Object$entries.length; _i++) { var _Object$entries$_i = (0, _slicedToArray2.default)(_Object$entries[_i], 2), key = _Object$entries$_i[0], value = _Object$entries$_i[1]; var idx = columns.indexOf(key); properties[value] = row[idx]; } } else { for (var i = 0; i < columns.length; i++) { if (i === idIdx || i === geomColumnIdx) { continue; } var columnName = columns[i]; properties[columnName] = row[i]; } } return { id: id, type: 'Feature', geometry: geometry, properties: properties }; } function getGeopackageVersion(db) { var textDecoder = new TextDecoder(); var applicationIdQuery = db.exec('PRAGMA application_id;')[0]; var applicationId = applicationIdQuery.values[0][0]; var buffer = new ArrayBuffer(4); var view = new DataView(buffer); view.setInt32(0, Number(applicationId)); var versionString = textDecoder.decode(buffer); if (versionString === 'GP10') { return '1.0'; } if (versionString === 'GP11') { return '1.1'; } var userVersionQuery = db.exec('PRAGMA user_version;')[0]; var userVersionInt = userVersionQuery.values[0][0]; if (userVersionInt && userVersionInt < 10300) { return '1.2'; } return null; } function getFeatureIdName(db, tableName) { var stmt = db.prepare("PRAGMA table_info(`".concat(tableName, "`)")); while (stmt.step()) { var pragmaTableInfo = stmt.getAsObject(); var name = pragmaTableInfo.name, pk = pragmaTableInfo.pk; if (pk) { return name; } } return null; } function parseGeometry(arrayBuffer) { var view = new DataView(arrayBuffer); var _parseGeometryBitFlag = parseGeometryBitFlags(view.getUint8(3)), envelopeLength = _parseGeometryBitFlag.envelopeLength, emptyGeometry = _parseGeometryBitFlag.emptyGeometry; if (emptyGeometry) { return null; } var wkbOffset = 8 + envelopeLength; var binaryGeometry = _loaders.WKBLoader.parseSync(arrayBuffer.slice(wkbOffset)); return (0, _loaders3.binaryToGeometry)(binaryGeometry); } function parseGeometryBitFlags(byte) { var envelopeValue = (byte & 14) / 2; var envelopeLength = ENVELOPE_BYTE_LENGTHS[envelopeValue]; return { littleEndian: Boolean(byte & 1), envelopeLength: envelopeLength, emptyGeometry: Boolean(byte & 16), extendedGeometryType: Boolean(byte & 32) }; } function getGeometryColumn(db, tableName) { var stmt = db.prepare('SELECT * FROM gpkg_geometry_columns WHERE table_name=:tableName;'); stmt.bind({ ':tableName': tableName }); stmt.step(); var geometryColumn = stmt.getAsObject(); return geometryColumn; } function getDataColumns(db, tableName) { var stmt; try { stmt = db.prepare('SELECT * FROM gpkg_data_columns WHERE table_name=:tableName;'); } catch (error) { if (error.message.includes('no such table')) { return null; } throw error; } stmt.bind({ ':tableName': tableName }); var result = {}; while (stmt.step()) { var column = stmt.getAsObject(); var column_name = column.column_name, name = column.name; result[column_name] = name || null; } return result; } function getArrowSchema(db, tableName) { var stmt = db.prepare("PRAGMA table_info(`".concat(tableName, "`)")); var fields = []; while (stmt.step()) { var pragmaTableInfo = stmt.getAsObject(); var name = pragmaTableInfo.name, type = pragmaTableInfo.type, notnull = pragmaTableInfo.notnull; var field = new _loaders2.Field(name, new SQL_TYPE_MAPPING[type](), !notnull); fields.push(field); } return new _loaders2.Schema(fields); } //# sourceMappingURL=parse-geopackage.js.map