30:T891,map = { /// Setup the map and container let container = html`

`; yield container; /// This is the velocity file const fileurl = await plottingFile; /// Setup the map const map = new WebMap({ portalItem: { id: "86265e5a4bbb4187a59719cf134e0018" } }); /// Setup the view const view = new MapView({ container, map: map, center: [-100, 40], zoom: 4, navigation: { mouseWheelZoomEnabled: true } }); // This file is available at the url. KML formatting is not great in ArcGIS const pbLayer = new KMLLayer({ url: pbUrl, copyright: "USGS Plate Boundaries", visible: options.includes(4) }); // This is one of the only files native to the notebook const volcLayer = new CSVLayer({ url: volcFile, copyright: "Global Volcanism Program", latitudeField: ["Latitude"], longitudeField: ["Longitude"], renderer: volcRender, visible: options.includes(7), popupTemplate: volcTemplate }); const bsmLayer = new CSVLayer({ url: bsmFile, latitudeField: ["Latitude"], longitudeField: ["Longitude"], renderer: bsmRender, visible: options.includes(3), popupTemplate: bsmTemplate }); // Fault layer const faultLayer = new GeoJSONLayer({ url: faultUrl, copyright: "The GEM Global Active Faults Database", popupTemplate: { title: "{name}" }, visible: options.includes(5) }); // This is one option for plotting earthquakes that occured in the last month, // see the url and VV notebook for other potential earthquake files const earthquakeLayer = new GeoJSONLayer({ url: eqUrl, copyright: "USGS Earthquakes", popupTemplate: eqTemplate, renderer: eqRenderer, visible: options.includes(6) }); const lsLayer = new GeoJSONLayer({ url: lsUrl, copyright: "NASA Open Data Portal", popupTemplate: lsTemplate, renderer: lsRenderer, visible: options.includes(2) }); /////////////////////////// // Adding the layers, legend, etc to the map map.add(pbLayer); map.add(faultLayer); map.add(volcLayer); map.add(lsLayer); map.add(earthquakeLayer); map.add(bsmLayer); }31:T5b1,function ellipseInfo(sigE, sigN, rEN) { /// Returns ellipse semi-major and semi-minor axis /// and angle of major axis CCW from E (x) axis //ellipse scale for 95% confidence var s = 5.991; //build covariance matrix var a = sigE ** 2; var b = rEN * sigE * sigN; var c = rEN * sigE * sigN; var d = sigN ** 2; //eigenvectors V var tmp = Math.sqrt((a - d) * (a - d) + 4 * b * c); var denom = 2 * c; var denom = 0 ? denom : 2; // special because c can be 0, don't divide by 0! var V = [ [-(tmp - a + d) / denom, (tmp + a - d) / denom], [1, 1] ]; //sqrt of eigenvalues D*scale var eigv1 = (a + d - tmp) / 2; var eigv2 = (a + d + tmp) / 2; var sqrtD = [Math.sqrt(s * eigv1), Math.sqrt(s * eigv2)]; //normalized eigenvectors var norm1 = Math.hypot(V[0][0], 1); var norm2 = Math.hypot(V[0][1], 1); V[0][0] /= norm1; V[1][0] /= norm1; V[0][1] /= norm2; V[1][1] /= norm2; //index of larger eigenvalue var ndx = sqrtD[0] < sqrtD[1] ? 1 : 0; //scaled larger eigenvector var x1 = V[0][ndx] * sqrtD[ndx]; var y1 = V[1][ndx] * sqrtD[ndx]; //scaled smaller eigenvector var x2 = V[0][1 - ndx] * sqrtD[1 - ndx]; var y2 = V[1][1 - ndx] * sqrtD[1 - ndx]; // angle of larger eigenvector, // semi-major, semi-minor length var ang = (Math.atan2(y1, x1) * 180) / Math.PI; var ang = NaN ? 0 : ang; var major = Math.hypot(x1, y1); var minor = Math.hypot(x2, y2); return [major, minor, ang]; }32:[{"id":"56cf6c14609697ae81713a791cf719c0f7ebd079c9cfbc835cc356b8c413cd18da27c4da5d2b41c44be4d75e8b956f8449dd0608e537948a7cba742331ef0a68","href":"https://static.observableusercontent.com/files/56cf6c14609697ae81713a791cf719c0f7ebd079c9cfbc835cc356b8c413cd18da27c4da5d2b41c44be4d75e8b956f8449dd0608e537948a7cba742331ef0a68","name":"GSRMv21_vel_1deg_Pacific_4deg_grid.gsrm.csv","type":"text/csv","lastModified":1648562188982.012,"size":185051},{"id":"2179da05c088ef1f3c5b3c97fdc71cd293240529779226dbedec94b82ecdcce74d7751c8049933bcf56fe1a510bdc00d716c00190bf13dcaf7d7a2d4bcbba135","href":"https://static.observableusercontent.com/files/2179da05c088ef1f3c5b3c97fdc71cd293240529779226dbedec94b82ecdcce74d7751c8049933bcf56fe1a510bdc00d716c00190bf13dcaf7d7a2d4bcbba135","name":"GVP_Volcano_List_Holocene.csv","type":"text/csv","lastModified":1648562188982.012,"size":268439}] 2e:{"title":"ArcGIS Plotting","cells":[{"id":0,"value":"md`# ArcGIS Plotting`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1593,"value":"viewof options = Inputs.checkbox(\n new Map([\n [\"Landslides\", 2],\n [\"Borehole Strainmeters\", 3],\n [\"Tectonic Plates\", 4],\n [\"Faults\", 5],\n [\"Earthquakes\", 6],\n [\"Volcanos\", 7],\n ]),\n {\n label: md`**Mapping Options**`,\n value: [3,7]\n }\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":580,"value":"$30","pinned":false,"hidden":false,"mode":"ojs"},{"id":2001,"value":"md`### Everything for the Supporting Map Objects`","pinned":false,"hidden":false,"mode":"ojs"},{"id":2004,"value":"pbUrl = \"https://prd-wret.s3.us-west-2.amazonaws.com/assets/palladium/production/atoms/files/plate-boundaries.kmz\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":2011,"value":"volcFile = FileAttachment(\"GVP_Volcano_List_Holocene.csv\").url()\n\n","pinned":false,"hidden":false,"mode":"ojs"},{"id":2017,"value":"volcSymbol = Object({\n type: \"simple-marker\", // autocasts as new SimpleMarkerSymbol()\n style: \"triangle\",\n color: \"orange\",\n size: \"8px\", // pixels\n outline: {\n color: \"black\",\n width: 0.5 // points\n }\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":2019,"value":"volcTemplate = Object({\n title: \"{Volcano Name}\",\n content:\n \"Latitude: {Latitude}, Longitude: {Longitude}
Last known eruption: {Last Known Eruption}\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":2021,"value":"volcRender = Object({\n type: \"simple\",\n field: \"Volcano\",\n value: \"yes\",\n symbol: volcSymbol,\n visualVariables: [\n {\n type: \"size\",\n field: \"Volcano\",\n\n valueExpression: \"$view.scale\",\n stops: [{ size: 10, value: 2888900 }, { size: 0.5, value: 147914381 }]\n }\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":2033,"value":"faultUrl = \"https://raw.githubusercontent.com/GEMScienceTools/gem-global-active-faults/master/geojson/gem_active_faults.geojson\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":2038,"value":"eqUrl = \"https://earthquake.usgs.gov/earthquakes/feed/v1.0/summary/4.5_month.geojson\"\n// \"https://earthquake.usgs.gov/earthquakes/feed/v1.0/summary/all_month.geojson\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":2043,"value":"eqTemplate = Object({\n title: \"Earthquake Info\",\n content: \"Magnitude {mag} {type} hit {place} on {time}\",\n fieldInfos: [\n {\n fieldName: \"time\",\n format: {\n dateFormat: \"short-date-short-time\"\n }\n }\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":2039,"value":"eqRenderer = Object({\n type: \"simple\",\n field: \"mag\",\n symbol: {\n type: \"simple-marker\",\n color: \"blue\",\n outline: {\n color: \"white\"\n }\n },\n visualVariables: [\n {\n type: \"size\",\n field: \"mag\",\n stops: [\n {\n value: 2.5,\n size: \"4px\"\n },\n {\n value: 8,\n size: \"40px\"\n }\n ]\n }\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":2411,"value":"bsmFile = \"https://raw.githubusercontent.com/cehanagan/ObservableSupport/main/PBO_bsm_coords.csv\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":2413,"value":"bsmSymbol = Object({\n type: \"simple-marker\", // autocasts as new SimpleMarkerSymbol()\n style: \"circle\",\n color: \"magenta\",\n size: \"9px\", // pixels\n outline: {\n color: \"black\",\n width: 0.5 // points\n }\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":2418,"value":"bsmTemplate = Object({\n title:\n \"{ID} Info\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":2415,"value":"bsmRender = Object({\n type: \"simple\",\n field: \"BSM\",\n value: \"yes\",\n symbol: bsmSymbol,\n visualVariables: [\n {\n type: \"size\",\n field: \"BSM\",\n\n valueExpression: \"$view.scale\",\n stops: [{ size: 10, value: 2888900 }, { size: 0.5, value: 147914381 }]\n }\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":2456,"value":"lsUrl = \"https://raw.githubusercontent.com/cehanagan/ObservableSupport/main/GlobalLandslideCatalog.geojson\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":2464,"value":"lsTemplate = Object({\n title: \"Landslide Info\",\n content:\n \"{landslide1} landslide on
{date_}
{source_lin}\",\n fieldInfos: [\n {\n fieldName: \"date_\",\n format: {\n dateFormat: \"short-date-short-time\"\n }\n }\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":2461,"value":"lsRenderer = Object({\n type: \"simple\",\n field: \"ls\",\n symbol: {\n type: \"simple-marker\",\n color: \"green\",\n size: \"9px\",\n outline: {\n color: \"white\"\n }\n },\n visualVariables: [\n {\n type: \"size\",\n field: \"ls\",\n valueExpression: \"$view.scale\",\n stops: [\n { size: 8, value: 2888900 },\n { size: 0.4, value: 147914381 }\n ]\n }\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1717,"value":"md`### Velocity Files`","pinned":false,"hidden":false,"mode":"ojs"},{"id":307,"value":"file = FileAttachment(\"GSRMv21_vel_1deg_Pacific_4deg_grid.gsrm.csv\").csv() //FileAttachment(\"cwu.final_igs14_cut.vel.csv\").csv({ typed: true })","pinned":false,"hidden":false,"mode":"ojs"},{"id":1740,"value":"downloadFile = file.map((d) =>\n Object({\n Station: d[\"Sta\"], //Sta or *Dot#\n Longitude: d[\"Ref_Elong\"] > 180 ? d[\"Ref_Elong\"] - 360 : d[\"Ref_Elong\"],\n Latitude: d[\"Ref_Nlat\"],\n NVel: d[\"dN/dt\"], // mm/yr (everything)\n Evel: d[\"dE/dt\"],\n UVel: d[\"dU/dt\"],\n Azimuth: (Math.atan2(d[\"dE/dt\"], d[\"dN/dt\"]) * 180) / Math.PI,\n Speed: Math.sqrt(d[\"dE/dt\"] ** 2 + d[\"dN/dt\"] ** 2),\n sigN: d[\"SNd\"],\n sigE: d[\"SEd\"],\n sigU: d[\"SUd\"],\n rEN: d[\"Rne\"],\n EllipseSemiMaj: ellipseInfo(d[\"SEd\"], d[\"SNd\"], d[\"Rne\"])[0],\n EllipseSemiMin: ellipseInfo(d[\"SEd\"], d[\"SNd\"], d[\"Rne\"])[1],\n EllipseXAng: ellipseInfo(d[\"SEd\"], d[\"SNd\"], d[\"Rne\"])[2]\n })\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":1745,"value":"plottingFile = \"https://raw.githubusercontent.com/cehanagan/VelocityViewer/main/VelocityFiles/GSRM_GPS_PA_editted.csv\"\n\n//\"https://raw.githubusercontent.com/cehanagan/VelocityViewer/main/VelocityFiles/cwu_final_nam14_editted.csv\"\n\n//\"https://raw.githubusercontent.com/cehanagan/VelocityViewer/main/VelocityFiles/GSRM_GPS_AF_editted.csv\"\n\n//FileAttachment(\"cwu_final_nam14_editted.csv\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":1754,"value":"filecsv = d3.csv(plottingFile) // plottingFile.csv({ typed: true })","pinned":false,"hidden":false,"mode":"ojs"},{"id":1971,"value":"md`### Constants from File`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1757,"value":"maxspeed = Number(d3.max(filecsv.map((d) => d.Speed)))","pinned":false,"hidden":false,"mode":"ojs"},{"id":2102,"value":"maxup = d3.max(filecsv.map(d => d.UVel))","pinned":false,"hidden":false,"mode":"ojs"},{"id":2109,"value":"minup = d3.min(filecsv.map(d => d.UVel))","pinned":false,"hidden":false,"mode":"ojs"},{"id":2310,"value":"maxuperror = d3.max(filecsv.map(d => d.sigE))","pinned":false,"hidden":false,"mode":"ojs"},{"id":2248,"value":"quantile95up = d3.quantile(filecsv.map(d => d.UVel), 0.95)","pinned":false,"hidden":false,"mode":"ojs"},{"id":2256,"value":"quantile5up = d3.quantile(filecsv.map(d => d.UVel), 0.05)","pinned":false,"hidden":false,"mode":"ojs"},{"id":1714,"value":"md`### Functions`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1477,"value":"$31","pinned":false,"hidden":false,"mode":"ojs"},{"id":844,"value":"function createEllipse(a, b, x, y, rotation) {\n rotation = (rotation / 180) * Math.PI;\n var n = (n = Math.ceil(36 * Math.max(a / b, b / a))); // n sampling angles, more for more elongated ellipses\n var coords = [];\n\n for (var i = 0; i <= n; i++) {\n // get the current angle\n var θ = ((Math.PI * 2) / n) * i + rotation;\n\n // get the radius at that angle\n var r =\n (a * b) /\n Math.sqrt(\n a * a * Math.sin(θ) * Math.sin(θ) + b * b * Math.cos(θ) * Math.cos(θ)\n );\n\n // get the x,y coordinate that marks the ellipse at this angle\n const x1 = x + Math.cos(θ - rotation) * r;\n const y1 = y + Math.sin(θ - rotation) * r;\n\n coords.push([x1, y1]);\n }\n\n // return a geojson object:\n return { type: \"Polygon\", coordinates: [coords] };\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1700,"value":"md`### ArcGIS Class Imports and Document Setup`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1427,"value":"KMLLayer = (await import(\"https://js.arcgis.com/4.20/@arcgis/core/layers/KMLLayer.js\"))\n .default","pinned":false,"hidden":false,"mode":"ojs"},{"id":788,"value":"Graphic = (await import(\"https://js.arcgis.com/4.20/@arcgis/core/Graphic.js\"))\n .default","pinned":false,"hidden":false,"mode":"ojs"},{"id":983,"value":"CSVLayer = (\n await import(\"https://js.arcgis.com/4.20/@arcgis/core/layers/CSVLayer.js\")\n).default","pinned":false,"hidden":false,"mode":"ojs"},{"id":746,"value":"GeoJSONLayer = (await import(\"https://js.arcgis.com/4.20/@arcgis/core/layers/GeoJSONLayer.js\"))\n .default","pinned":false,"hidden":false,"mode":"ojs"},{"id":676,"value":"Legend = (await import(\"https://js.arcgis.com/4.20/@arcgis/core/widgets/Legend.js\"))\n .default","pinned":false,"hidden":false,"mode":"ojs"},{"id":565,"value":"FeatureLayer = (await import(\"https://js.arcgis.com/4.20/@arcgis/core/layers/FeatureLayer.js\"))\n .default","pinned":false,"hidden":false,"mode":"ojs"},{"id":22,"value":"WebMap = (await import(\"https://js.arcgis.com/4.20/@arcgis/core/WebMap.js\"))\n .default","pinned":false,"hidden":false,"mode":"ojs"},{"id":20,"value":"MapView = (await import(\"https://js.arcgis.com/4.20/@arcgis/core/views/MapView.js\"))\n .default","pinned":false,"hidden":false,"mode":"ojs"},{"id":14,"value":"document.head.appendChild(\n html`\n`\n)","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W32","stdlib":"1","titling":"auto"}