30:T190d,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 } }); //////////////// //Testing, this sets up a feature layer using a combination of geometries //and symbols to use with CIMSymbols. This is a potential option to //explore for the arrows and error ellipses, becuase it allows custom //symbols to be plotted in the same layer, and the ellipses may be plotted at //the end of the arrow. const polygon = { type: "polygon", // autocasts as new Polygon() rings: [[-64.78, 32.3], [-66.07, 18.45], [-80.21, 25.78], [-64.78, 32.3]] }; const lineSymbol = { type: "simple-fill", // autocasts as SimpleLineSymbol() color: [226, 119, 40, 0], outline: { color: [226, 119, 40, 1], width: 1 } }; const polylineGraphic = [ { geometry: polygon, // Add the geometry created in step 4 attributes: { OBJECTID: 1 } } ]; const testLayer = new FeatureLayer({ source: polylineGraphic, objectIdField: "OBJECTID", //renderer: { // type: "simple", // field: "OBJECTID", // value: 1, // symbol: lineSymbol //} renderer: { type: "simple", symbol: { type: "cim", // CIM Line Symbol data: { type: "CIMSymbolReference", symbol: { type: "CIMLineSymbol", symbolLayers: [ { // white dashed layer at center of the line type: "CIMSolidStroke", effects: [ { type: "CIMGeometricEffectDashes", dashTemplate: [2, 2, 2, 2], // width of dashes and spacing between the dashes lineDashEnding: "NoConstraint", controlPointEnding: "NoConstraint" } ], enable: true, // must be set to true in order for the symbol layer to be visible capStyle: "Butt", joinStyle: "Round", width: 1, color: [255, 255, 255, 255] }, { // lighter green line layer that surrounds the dashes type: "CIMSolidStroke", enable: true, capStyle: "Butt", joinStyle: "Round", width: 3, color: [56, 168, 0, 255] }, { // darker green outline around the line symbol type: "CIMSolidStroke", enable: true, capStyle: "Butt", joinStyle: "Round", width: 6, color: [0, 115, 76, 255] } ] } } } } }); map.add(testLayer); /////////////////// // Defining layers // Horizontal velocity vectors const arrowLayer = new CSVLayer({ url: fileurl, latitudeField: ["Latitude"], longitudeField: ["Longitude"], popupTemplate: arrowTemplate, renderer: arrowRender, visible: options.includes(1) }); // Station Labels, with invisible placemark // (when placed in arrow layer, the labels were too far from the station) const stationLayer = new CSVLayer({ url: fileurl, latitudeField: ["Latitude"], longitudeField: ["Longitude"], popupTemplate: arrowTemplate, renderer: { type: "simple", symbol: { type: "simple-marker", color: [0, 0, 0, 0], outline: { color: [0, 0, 0, 0] } } }, labelingInfo: stationLabel, labelsVisible: options.includes(3) }); // This needs major work... but has a placeholder svg path const ellipseLayer = new CSVLayer({ url: fileurl, latitudeField: ["Latitude"], longitudeField: ["Longitude"], renderer: ellipseRender, visible: options.includes(2) }); // 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 }); // 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) }); // Vertical Velocity layer const upVelLayer = new CSVLayer({ url: fileurl, latitudeField: ["Latitude"], longitudeField: ["Longitude"], renderer: upVelRender, popupTemplate: arrowTemplate, visible: options.includes(8) }); // Vertical Velocity Error Layer const upErrorLayer = new CSVLayer({ url: fileurl, latitudeField: ["Latitude"], longitudeField: ["Longitude"], renderer: upErrorRender, visible: options.includes(9) }); /////////////////////////// // Adding the layers, legend, etc to the map map.add(pbLayer); map.add(faultLayer); map.add(volcLayer); map.add(earthquakeLayer); map.add(upVelLayer); map.add(upErrorLayer); map.add(ellipseLayer); map.add(stationLayer); map.add(arrowLayer); // This legend needs help const legend = new Legend({ // Need to Fix view: view, layerInfos: [ { layer: arrowLayer, title: "Speed" } ] }); // Add widget to the bottom right corner of the view // This adds the legend //view.ui.add(legend, "bottom-right"); }31:Tae3,arrowTemplate = Object({ title: "{Station} Info", // content: arrowTempText, content: [ { type: "fields", fieldInfos: [ { fieldName: "Latitude", visible: true, format: { places: 2, digitSeparator: true } }, { fieldName: "Longitude", visible: true, format: { places: 2, digitSeparator: true } }, { fieldName: "Speed", label: "Horizontal Speed (mm/yr)", visible: true, format: { places: 2, digitSeparator: true } }, { fieldName: "Azimuth", visible: true, format: { places: 2, digitSeparator: true } }, { fieldName: "NVel", label: "North Velocity (mm/yr)", visible: true, format: { places: 2, digitSeparator: true } }, { fieldName: "Evel", label: "East Velocity (mm/yr)", visible: true, format: { places: 2, digitSeparator: true } }, { fieldName: "UVel", label: "Up Velocity (mm/yr)", visible: true, format: { places: 2, digitSeparator: true } }, { fieldName: "sigN", label: "North Velocity Standard Deviation (mm/yr)", visible: true, format: { places: 2, digitSeparator: true } }, { fieldName: "sigE", label: "East Velocity Standard Deviation (mm/yr)", visible: true, format: { places: 2, digitSeparator: true } }, { fieldName: "sigU", label: "Up Velocity Standard Deviation (mm/yr)", visible: true, format: { places: 2, digitSeparator: true } }, { fieldName: "rEN", label: "East/North Velocity Correlation Coefficient", visible: true, format: { places: 2, digitSeparator: true } } ] }, { type: "text", text: " Click for NAM14 ASCII Time Series

" } ] })32:T8d9,md`### File Setup 1. (With awk or python) I'm sure there is a better way to do this, but I took the header info out of the velocity files, replaced tabs with commas, deleted repeat rows for the station (I kept the latest velocity info only), and uploaded it here. [Brooks, if you know how to do this in Javascript, and link it to a url, we can do it all in one step. I haven't figured it out yet] 2. I used the downloadFile cell below to tranform the column entries into a format that is useful for plotting on the map above. 3. Using the download csv button on the downloadFile cell, I downloaded the file, then re-uploaded it as the plottingFile (without letting it read the csv format into an array, because the arcgis mapping tool for csvs expects a url to the file) The complete list of files: The files ready for plotting are inmy github at https://github.com/cehanagan/VelocityViewer/tree/main/VelocityFiles cwu_final_nam14_editted.csv (unavco source, NAM14 reference [stable north america]) cwu_final_igs14_editted.csv (unavco source, IGS14 reference [global reference frame]) GSRM_GPS_AF_editted.csv (GSRM source, IGS08 system, Africa) GSRM_GPS_AN_editted.csv (GSRM source, IGS08 system, Antartica) GSRM_GPS_AR_editted.csv (GSRM source, IGS08 system, Arabia) GSRM_GPS_AU_editted.csv (GSRM source, IGS08 system, Australia) GSRM_GPS_CA_editted.csv (GSRM source, IGS08 system, Carribean) GSRM_GPS_CO_editted.csv (GSRM source, IGS08 system, Cocos) GSRM_GPS_EU_editted.csv (GSRM source, IGS08 system, Europe) GSRM_GPS_IN_editted.csv (GSRM source, IGS08 system, India) GSRM_GPS_JF_editted.csv (GSRM source, IGS08 system, Juan de Fuca) GSRM_GPS_NA_editted.csv (GSRM source, IGS08 system, North America) GSRM_GPS_NZ_editted.csv (GSRM source, IGS08 system, New Zealand) GSRM_GPS_PA_editted.csv (GSRM source, IGS08 system, Pacific) GSRM_GPS_PS_editted.csv (GSRM source, IGS08 system, Philippine Sea) GSRM_GPS_SA_editted.csv (GSRM source, IGS08 system, South America) GSRM_GPS_SC_editted.csv (GSRM source, IGS08 system, Scotia) GSRM_4degGRID_ITRF08_editted.csv (GSRM source, ITRF08, Global) GSRM_4degGRID_JF_editted.csv (GSRM source, Juan de Fuca) GSRM_4degGRID_NA_editted.csv (GSRM source, North America) GSRM_4degGRID_PA_editted.csv (GSRM source, Pacific) `33: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]; }34:[{"id":"936dc54a8f43f51e70fdc3856e4ffa4f77b675e8f4184a9f6242520f7f659aa0c1dacfd414c44b463cc0e1879dc37fe0be77d5982a1f052e1af093ba923a1b29","href":"https://static.observableusercontent.com/files/936dc54a8f43f51e70fdc3856e4ffa4f77b675e8f4184a9f6242520f7f659aa0c1dacfd414c44b463cc0e1879dc37fe0be77d5982a1f052e1af093ba923a1b29","name":"transformedFile (1).csv","type":"text/csv","lastModified":1627503038172.952,"size":1294926},{"id":"7b0a8b107631b0a167db5aa1cfb893261a1e8e19823394e60ce5eabe763af729ab4e05b79b70e3401056dbec971792da414464de664edf83d71e89d485d67fb4","href":"https://static.observableusercontent.com/files/7b0a8b107631b0a167db5aa1cfb893261a1e8e19823394e60ce5eabe763af729ab4e05b79b70e3401056dbec971792da414464de664edf83d71e89d485d67fb4","name":"cwu.final_nam14_cut.vel.csv","type":"text/csv","lastModified":1627659724226.881,"size":867981},{"id":"ce4e3a150fbc877ab2beac9b2768ce94894cf192e177b650a220feb0a22b68a8993a5b98447b2334d5253b44b1a4f5f9e481d2d1d7e9988ecd6da5c514b084f9","href":"https://static.observableusercontent.com/files/ce4e3a150fbc877ab2beac9b2768ce94894cf192e177b650a220feb0a22b68a8993a5b98447b2334d5253b44b1a4f5f9e481d2d1d7e9988ecd6da5c514b084f9","name":"downloadFile (1).csv","type":"text/csv","lastModified":1627659803686.986,"size":524803},{"id":"2179da05c088ef1f3c5b3c97fdc71cd293240529779226dbedec94b82ecdcce74d7751c8049933bcf56fe1a510bdc00d716c00190bf13dcaf7d7a2d4bcbba135","href":"https://static.observableusercontent.com/files/2179da05c088ef1f3c5b3c97fdc71cd293240529779226dbedec94b82ecdcce74d7751c8049933bcf56fe1a510bdc00d716c00190bf13dcaf7d7a2d4bcbba135","name":"GVP_Volcano_List_Holocene.csv","type":"text/csv","lastModified":1627570352968.358,"size":268439},{"id":"56cf6c14609697ae81713a791cf719c0f7ebd079c9cfbc835cc356b8c413cd18da27c4da5d2b41c44be4d75e8b956f8449dd0608e537948a7cba742331ef0a68","href":"https://static.observableusercontent.com/files/56cf6c14609697ae81713a791cf719c0f7ebd079c9cfbc835cc356b8c413cd18da27c4da5d2b41c44be4d75e8b956f8449dd0608e537948a7cba742331ef0a68","name":"GSRMv21_vel_1deg_Pacific_4deg_grid.gsrm.csv","type":"text/csv","lastModified":1627669773202.593,"size":185051},{"id":"7eabf31e262ada16d46e658366f2dcc99378a96d9220d19bd043da023619cafbc8e188493f42b19096975601e3be185f4c6547b5324ca6bd8fa16de7c754061f","href":"https://static.observableusercontent.com/files/7eabf31e262ada16d46e658366f2dcc99378a96d9220d19bd043da023619cafbc8e188493f42b19096975601e3be185f4c6547b5324ca6bd8fa16de7c754061f","name":"cwu_final_nam14_editted.csv","type":"text/csv","lastModified":1627670708242.388,"size":471228},{"id":"6d22fb6d966d0dbe0a9578d4ddd483d94936230796e033cb0690c681dec5c412622244ed66511acece98949c5177fd9d0e63d1c060df0fb7cb2995fb5295bb65","href":"https://static.observableusercontent.com/files/6d22fb6d966d0dbe0a9578d4ddd483d94936230796e033cb0690c681dec5c412622244ed66511acece98949c5177fd9d0e63d1c060df0fb7cb2995fb5295bb65","name":"cwu.final_igs14_cut.vel.csv","type":"text/csv","lastModified":1626199176075.12,"size":2375904}] 2e:{"title":"Velocity Viewer Support","cells":[{"id":0,"value":"md`# Velocity Viewer Support`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1593,"value":"viewof options = Inputs.checkbox(\n new Map([\n [\"Velocity Vectors (mm/yr)\", 1],\n [\"Error Ellipses (mm/yr)\", 2],\n [\"Station Labels\", 3],\n [\"Tectonic Plates\", 4],\n [\"Faults\", 5],\n [\"Earthquakes\", 6],\n [\"Volcanos\", 7],\n [\"Up Velocity (mm/yr)\", 8],\n [\"Up Velocity Error (mm/yr)\", 9]\n ]),\n {\n label: md`**Mapping Options**`,\n value: [1, 3]\n }\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":580,"value":"$30","pinned":false,"hidden":false,"mode":"ojs"},{"id":1953,"value":"md`### Everything for Velocity Plotting`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1980,"value":"arrowRender = Object({\n type: \"simple\",\n //field: \"Arrow\",\n //value: \"yes\",\n symbol: arrowSymbol,\n visualVariables: [\n { type: \"rotation\", field: \"Azimuth\", rotationType: \"geographic\" },\n {\n type: \"size\",\n field: \"Speed\",\n minDataValue: 0,\n maxDataValue: maxspeed,\n maxSize: {\n type: \"size\",\n valueExpression: \"$view.scale\",\n stops: [{ size: 200, value: 288890 }, { size: 60, value: 147914381 }]\n },\n minSize: {\n type: \"size\",\n valueExpression: \"$view.scale\",\n stops: [{ size: 3, value: 288890 }, { size: 1, value: 147914381 }]\n }\n }\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1977,"value":"stationLabel = Object({\n symbol: {\n type: \"text\", // autocasts as new TextSymbol()\n color: \"magenta\",\n haloColor: \"black\",\n haloSize: 1,\n font: {\n size: 9,\n weight: \"bold\"\n }\n },\n labelPlacement: \"above-center\",\n labelExpressionInfo: {\n expression: \"$feature.Station\"\n }\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1962,"value":"arrowSymbol = Object({\n type: \"simple-marker\", // autocasts as new SimpleMarkerSymbol()\n // use an SVG path to create an arrow shape\n path: \"M-0.5 0L-0.5 -7L-2 -7L0 -9L2 -7L 0.5 -7L0.5 0ZM0 0L0 9\",\n color: \"yellow\",\n outline: {\n color: \"black\",\n width: 0.5\n }\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1951,"value":"$31","pinned":false,"hidden":false,"mode":"ojs"},{"id":2143,"value":"upVelSymbol = Object({\n type: \"simple-marker\", // autocasts as new SimpleMarkerSymbol()\n style: \"triangle\",\n color: [0, 0, 0, 0],\n outline: {\n color: \"black\",\n width: 0.5\n }\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":2159,"value":"downVelSymbol = Object({\n type: \"simple-marker\", // autocasts as new SimpleMarkerSymbol()\n style: \"triangle\",\n angle: 180,\n color: [0, 0, 0, 0],\n outline: {\n color: \"black\",\n width: 0.5\n }\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":2100,"value":"upVelRender = Object({\n type: \"unique-value\",\n valueExpression:\n \"When($feature.UVel < 0.0, 'down', $feature.UVel >= 0.0, 'up','down')\",\n defaultSymbol: upVelSymbol,\n defaultLabel: \"Up Velocity (mm/yr)\",\n uniqueValueInfos: [\n { value: 'down', symbol: downVelSymbol },\n { value: 'up', symbol: upVelSymbol }\n ],\n visualVariables: [\n {\n type: \"size\",\n valueExpression: \"$view.scale\",\n stops: [{ size: 15, value: 2800 }, { size: 0.5, value: 100000000 }]\n },\n {\n type: \"color\",\n field: \"UVel\",\n legendOptions: {\n title: \"Fill Me In!\"\n },\n stops: [\n {\n value: quantile5up,\n color: \"blue\",\n label: \"S\" + String(quantile95up)\n }\n ]\n }\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":2302,"value":"upErrorRender = Object({\n type: \"simple\",\n symbol: {\n type: \"simple-marker\",\n color: [0, 0, 0, 0],\n outline: {\n color: \"orange\",\n width: 0.5\n }\n },\n visualVariables: [\n {\n type: \"size\",\n field: \"sigU\",\n minDataValue: 0,\n maxDataValue: maxuperror,\n maxSize: {\n type: \"size\",\n valueExpression: \"$view.scale\",\n stops: [{ size: 15, value: 2800 }, { size: 0.5, value: 100000000 }]\n },\n minSize: {\n type: \"size\",\n valueExpression: \"$view.scale\",\n stops: [{ size: 5, value: 2800 }, { size: 0.05, value: 100000000 }]\n }\n }\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1988,"value":"md`### Everything for Error Ellipse Plotting`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1986,"value":"ellipseSymbol = Object({\n type: \"simple-marker\", // autocasts as new SimpleMarkerSymbol()\n // use an SVG path to create an arrow shape\n path: \"M 0 0 a 1 1 0 1 0 0.001 0 z\",\n color: \"yellow\",\n outline: {\n color: \"black\",\n width: 0.5\n }\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1992,"value":"ellipseRender = Object({\n type: \"simple\",\n field: \"Ellipse\",\n symbol: ellipseSymbol\n})","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 =\n \"https://earthquake.usgs.gov/earthquakes/feed/v1.0/summary/all_month.geojson\"\n\n \n ","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":1717,"value":"$32","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":"$33","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":"$W34","stdlib":"1","titling":"auto"}