31:T69c,// Credit to RE:earth // https://github.com/reearth/reearth/blob/main/web/src/beta/core/mantle/data/shapefile.ts function parseDbf(dbf) { const dataView = new DataView(new Uint8Array(dbf).buffer); let idx = 4; const numberOfRecords = dataView.getInt32(idx, true); idx += 28; let end = false; const fields = []; while (!end) { const field = {}; const nameArray = []; for (let i = 0; i < 10; i++) { const letter = dataView.getUint8(idx); if (letter !== 0) { nameArray.push(String.fromCharCode(letter)); } idx += 1; } field.name = nameArray.join(""); idx += 1; field.type = String.fromCharCode(dataView.getUint8(idx)); idx += 5; field.fieldLength = dataView.getUint8(idx); idx += 16; fields.push(field); if (dataView.getUint8(idx) === 0x0d) { break; } } idx += 1; const propertiesArray = []; for (let i = 0; i < numberOfRecords; i++) { const properties = {}; if (!end) { try { idx += 1; for (let j = 0; j < fields.length; j++) { let str = ""; const charString = []; for (let h = 0; h < fields[j].fieldLength; h++) { charString.push(String.fromCharCode(dataView.getUint8(idx))); idx += 1; } str = charString.join(""); str = str.trim(); const number = parseFloat(str); if (isNaN(number)) { properties[fields[j].name] = str; } else { properties[fields[j].name] = number; } } } catch (err) { end = true; } } propertiesArray.push(properties); } return propertiesArray; }32:Ta9c,// Credit to RE:earth // https://github.com/reearth/reearth/blob/main/web/src/beta/core/mantle/data/shapefile.ts function parseShapefile(shp) { const dataView = new DataView(shp); let idx = 0; const wordLength = dataView.getInt32((idx += 6 * 4), false); const byteLength = wordLength * 2; idx += 4; //version idx += 4; //shapeType idx += 4; //minX, minY idx += 8 * 8; //min(Y, Z, M),max(X, Y, Z, M) const features = []; while (idx < byteLength) { const feature = {}; const length = dataView.getInt32((idx += 4), false); const type = dataView.getInt32((idx += 4), true); let idxFeature = idx + 4; let numberOfParts, nbpoints, numberOfPoints, nbpartsPoint = []; switch (type) { case 1: case 11: case 21: feature.type = "Point"; feature.coordinates = [ dataView.getFloat64(idxFeature, true), dataView.getFloat64(idxFeature + 8, true) ]; break; case 3: case 13: case 23: case 5: case 15: case 25: if (type === 3 || type === 13 || type === 23) { feature.type = "MultiLineString"; } else if (type === 5 || type === 15 || type === 25) { feature.type = "Polygon"; } numberOfParts = dataView.getInt32(idxFeature + 32, true); nbpoints = dataView.getInt32(idxFeature + 36, true); idxFeature += 40; nbpartsPoint = new Array(numberOfParts).fill(0).map(() => { const result = dataView.getInt32(idxFeature, true); idxFeature += 4; return result; }); feature.coordinates = new Array(numberOfParts).fill(0).map((_, i) => { const idstart = nbpartsPoint[i]; const idend = (i < numberOfParts - 1 ? nbpartsPoint[i + 1] : nbpoints) - 1; const part = []; for (let j = idstart; j <= idend; j++) { part.push([ dataView.getFloat64(idxFeature, true), dataView.getFloat64(idxFeature + 8, true) ]); idxFeature += 16; } return part; }); break; case 8: case 18: case 28: feature.type = "MultiPoint"; numberOfPoints = dataView.getInt32(idxFeature + 32, true); idxFeature += 36; feature.coordinates = new Array(numberOfPoints).fill(0).map(() => { const result = [ dataView.getFloat64(idxFeature, true), dataView.getFloat64(idxFeature + 8, true) ]; idxFeature += 16; return result; }); break; } idx += length * 2; features.push(feature); } return features; }33:[{"id":"93096a641bed29b9a3b2044d57e30b8857eeb888dcb3d6cc03bb0b419fb0e1048164478827614c9d9f48df18563a6db6e677065a166089326d92e3a56209cec5","href":"https://static.observableusercontent.com/files/93096a641bed29b9a3b2044d57e30b8857eeb888dcb3d6cc03bb0b419fb0e1048164478827614c9d9f48df18563a6db6e677065a166089326d92e3a56209cec5","name":"UGB_2021.zip","type":"application/zip","lastModified":1686103458190.097,"size":1228799}] 2f:{"title":"shp2geojson","cells":[{"id":0,"value":"# shp2geojson\n\nConvert Shapefile to GeoJSON. \n\n*NOTE: Requires .shp, .dbf, and .proj files.*\n\n### Usage\n~~~js\nimport {zip, zipreader} from \"@fil/jszip\"\n\n// Fetch your ZIP file\nshp_zip = await fetch(\"\").then((z) => zipreader(z, { type: \"arraybuffer\" }))\n\n// Convert to GeoJSON\nimport { geojson } with { my_shp_zip as shp_zip } from \"@rjsheperd/shp2geojson\"\n\n// TODO: Add example using FileAttachment\n// shp_zip = await FileAttachment(\".zip\").zip()\n~~~","pinned":false,"hidden":false,"mode":"md"},{"id":153,"value":"### Fetch Example","pinned":false,"hidden":false,"mode":"md"},{"id":12,"value":"shp_zip = await fetchp(\n \"https://ftp.gis.oregon.gov/framework/adminbound/UGB_2021.zip\"\n).then((z) => zipreader(z, { type: \"arraybuffer\" }))","pinned":false,"hidden":false,"mode":"ojs"},{"id":126,"value":"shp = await shp_zip.get(withExt(shp_zip, \"shp\"))","pinned":false,"hidden":false,"mode":"ojs"},{"id":144,"value":"dbf = await shp_zip.get(withExt(shp_zip, \"dbf\"))","pinned":false,"hidden":false,"mode":"ojs"},{"id":43,"value":"geojson = shp2geojson(shp, dbf, projection, proj4.WGS84)","pinned":false,"hidden":false,"mode":"ojs"},{"id":49,"value":"projection = {\n return new TextDecoder().decode(\n new Uint8Array(await shp_zip.get(withExt(shp_zip, \"prj\")))\n );\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":61,"value":"// Adapted from RE:earth\n// https://github.com/reearth/reearth/blob/main/web/src/beta/core/mantle/data/shapefile.ts\nfunction shp2geojson(shp, dbf, from_proj, to_proj) {\n const features = parseShapefile(shp);\n const propertiesArray = parseDbf(dbf);\n\n const geojson = {\n type: \"FeatureCollection\",\n features: []\n };\n for (let i = 0; i < Math.min(features.length, propertiesArray.length); i++) {\n geojson.features.push({\n type: \"Feature\",\n geometry: features[i],\n properties: propertiesArray[i]\n });\n }\n return reproject.reproject(geojson, from_proj, to_proj || proj4.WGS84);\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":29,"value":"download(new Blob([JSON.stringify(geojson)]), \"ugb.geojson\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":111,"value":"### Functions","pinned":false,"hidden":false,"mode":"md"},{"id":133,"value":"withExt = (zip, ext) => Array.from(zip.keys()).filter((f) => f.endsWith(ext))[0]","pinned":false,"hidden":false,"mode":"ojs"},{"id":58,"value":"$31","pinned":false,"hidden":false,"mode":"ojs"},{"id":64,"value":"$32","pinned":false,"hidden":false,"mode":"ojs"},{"id":15,"value":"### Imports","pinned":false,"hidden":false,"mode":"md"},{"id":3,"value":"import { zip, zipreader } from \"@fil/jszip\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":52,"value":"reproject = import(\"https://cdn.skypack.dev/reproject@1.2.7?min\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":18,"value":"import { fetchp } from \"@tomlarkworthy/fetchp\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":33,"value":"proj4 = require(\"proj4@2.9.0/dist/proj4-src.js\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":68,"value":"import { download } from \"@rjsheperd/gdal\"","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W33","stdlib":"1","titling":"auto"}