Created
August 10, 2023 22:01
-
-
Save scottyhq/687660d4b115868e7886c52fba47ccf0 to your computer and use it in GitHub Desktop.
icesat2 hackweek 2023 geopandas
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
{ | |
"cells": [ | |
{ | |
"cell_type": "code", | |
"execution_count": 1, | |
"id": "0a2951e5-1e9f-4338-867f-79f93de9c0d3", | |
"metadata": { | |
"tags": [] | |
}, | |
"outputs": [], | |
"source": [ | |
"# NOTE: REQUIRES PROJ>=9.2 and GEOPANDAS>=0.13\n", | |
"\n", | |
"#!pip install geopandas==0.13" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 2, | |
"id": "aff6d69e-07f9-4b6c-87d5-5ca7a13c08d0", | |
"metadata": { | |
"tags": [] | |
}, | |
"outputs": [ | |
{ | |
"name": "stdout", | |
"output_type": "stream", | |
"text": [ | |
"ON\n" | |
] | |
} | |
], | |
"source": [ | |
"import os\n", | |
"os.environ['USE_PYGEOS'] = '0'\n", | |
"os.environ['PROJ_DEBUG'] = '2' # For showing logs\n", | |
"#os.environ['PROJ_NETWORK'] = 'ON' # Ensure this is 'ON' to get shift grids over the internet\n", | |
"print(os.environ['PROJ_NETWORK']) \n", | |
"\n", | |
"import fiona\n", | |
"import geopandas as gpd\n", | |
"import shapely" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 3, | |
"id": "a2078edc-442b-4ee2-be8a-fd6f47ed0743", | |
"metadata": { | |
"tags": [] | |
}, | |
"outputs": [ | |
{ | |
"name": "stdout", | |
"output_type": "stream", | |
"text": [ | |
"\n", | |
"SYSTEM INFO\n", | |
"-----------\n", | |
"python : 3.10.12 | packaged by conda-forge | (main, Jun 23 2023, 22:40:32) [GCC 12.3.0]\n", | |
"executable : /srv/conda/envs/notebook/bin/python\n", | |
"machine : Linux-5.10.167-147.601.amzn2.x86_64-x86_64-with-glibc2.27\n", | |
"\n", | |
"GEOS, GDAL, PROJ INFO\n", | |
"---------------------\n", | |
"GEOS : 3.11.2\n", | |
"GEOS lib : None\n", | |
"GDAL : 3.7.0\n", | |
"GDAL data dir: /srv/conda/envs/notebook/share/gdal\n", | |
"PROJ : 9.2.1\n", | |
"PROJ data dir: /srv/conda/envs/notebook/share/proj\n", | |
"\n", | |
"PYTHON DEPENDENCIES\n", | |
"-------------------\n", | |
"geopandas : 0.13.0\n", | |
"numpy : 1.23.5\n", | |
"pandas : 1.5.1\n", | |
"pyproj : 3.6.0\n", | |
"shapely : 2.0.1\n", | |
"fiona : 1.9.4\n", | |
"geoalchemy2: None\n", | |
"geopy : None\n", | |
"matplotlib : 3.6.2\n", | |
"mapclassify: 2.5.0\n", | |
"pygeos : 0.14\n", | |
"pyogrio : None\n", | |
"psycopg2 : 2.9.6 (dt dec pq3 ext lo64)\n", | |
"pyarrow : 12.0.1\n", | |
"rtree : 1.0.1\n" | |
] | |
} | |
], | |
"source": [ | |
"gpd.show_versions()" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 4, | |
"id": "f4c5d69c-7a9b-48e3-a79e-13addc10ff35", | |
"metadata": { | |
"tags": [] | |
}, | |
"outputs": [], | |
"source": [ | |
"# from slidrule single granule, default parameters\n", | |
"#gf = icesat2.atl06p({}, resources=['ATL03_20181019224323_03250112_005_01.h5'])\n", | |
"gf = gpd.read_file('sliderule-output.gpkg',\n", | |
" rows=100) " | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 5, | |
"id": "20851bf4-27e5-4b53-a82d-610844621713", | |
"metadata": { | |
"tags": [] | |
}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"0 POINT (-100.63246 -78.98922)\n", | |
"1 POINT (-100.63265 -78.98904)\n", | |
"Name: geometry, dtype: geometry" | |
] | |
}, | |
"execution_count": 5, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"gf.geometry.head(2)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 6, | |
"id": "18957036-0c63-4c6d-b7ba-92d33d66c8ac", | |
"metadata": { | |
"tags": [] | |
}, | |
"outputs": [ | |
{ | |
"name": "stderr", | |
"output_type": "stream", | |
"text": [ | |
"/tmp/ipykernel_2695/4052444587.py:3: FutureWarning: In a future version, `df.iloc[:, i] = newvals` will attempt to set the values inplace instead of always setting a new array. To retain the old behavior, use either `df[df.columns[i]] = newvals` or, if columns are non-unique, `df.isetitem(i, newvals)`\n", | |
" points.loc[:, 'time'] = points.time.dt.strftime('%Y-%m-%d')\n" | |
] | |
}, | |
{ | |
"data": { | |
"text/html": [ | |
"<div style=\"width:100%;\"><div style=\"position:relative;width:100%;height:0;padding-bottom:60%;\"><span style=\"color:#565656\">Make this Notebook Trusted to load map: File -> Trust Notebook</span><iframe srcdoc=\"<!DOCTYPE html>\n", | |
"<html>\n", | |
"<head>\n", | |
" \n", | |
" <meta http-equiv="content-type" content="text/html; charset=UTF-8" />\n", | |
" \n", | |
" <script>\n", | |
" L_NO_TOUCH = false;\n", | |
" L_DISABLE_3D = false;\n", | |
" </script>\n", | |
" \n", | |
" <style>html, body {width: 100%;height: 100%;margin: 0;padding: 0;}</style>\n", | |
" <style>#map {position:absolute;top:0;bottom:0;right:0;left:0;}</style>\n", | |
" <script src="https://cdn.jsdelivr.net/npm/[email protected]/dist/leaflet.js"></script>\n", | |
" <script src="https://code.jquery.com/jquery-1.12.4.min.js"></script>\n", | |
" <script src="https://cdn.jsdelivr.net/npm/[email protected]/dist/js/bootstrap.bundle.min.js"></script>\n", | |
" <script src="https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.js"></script>\n", | |
" <link rel="stylesheet" href="https://cdn.jsdelivr.net/npm/[email protected]/dist/leaflet.css"/>\n", | |
" <link rel="stylesheet" href="https://cdn.jsdelivr.net/npm/[email protected]/dist/css/bootstrap.min.css"/>\n", | |
" <link rel="stylesheet" href="https://netdna.bootstrapcdn.com/bootstrap/3.0.0/css/bootstrap.min.css"/>\n", | |
" <link rel="stylesheet" href="https://cdn.jsdelivr.net/npm/@fortawesome/[email protected]/css/all.min.css"/>\n", | |
" <link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.css"/>\n", | |
" <link rel="stylesheet" href="https://cdn.jsdelivr.net/gh/python-visualization/folium/folium/templates/leaflet.awesome.rotate.min.css"/>\n", | |
" \n", | |
" <meta name="viewport" content="width=device-width,\n", | |
" initial-scale=1.0, maximum-scale=1.0, user-scalable=no" />\n", | |
" <style>\n", | |
" #map_8366bb0414240e77d3544f24a5d8ceb1 {\n", | |
" position: relative;\n", | |
" width: 100.0%;\n", | |
" height: 100.0%;\n", | |
" left: 0.0%;\n", | |
" top: 0.0%;\n", | |
" }\n", | |
" .leaflet-container { font-size: 1rem; }\n", | |
" </style>\n", | |
" \n", | |
" \n", | |
" <style>\n", | |
" .foliumtooltip {\n", | |
" \n", | |
" }\n", | |
" .foliumtooltip table{\n", | |
" margin: auto;\n", | |
" }\n", | |
" .foliumtooltip tr{\n", | |
" text-align: left;\n", | |
" }\n", | |
" .foliumtooltip th{\n", | |
" padding: 2px; padding-right: 8px;\n", | |
" }\n", | |
" </style>\n", | |
" \n", | |
"</head>\n", | |
"<body>\n", | |
" \n", | |
" \n", | |
" <div class="folium-map" id="map_8366bb0414240e77d3544f24a5d8ceb1" ></div>\n", | |
" \n", | |
"</body>\n", | |
"<script>\n", | |
" \n", | |
" \n", | |
" var map_8366bb0414240e77d3544f24a5d8ceb1 = L.map(\n", | |
" "map_8366bb0414240e77d3544f24a5d8ceb1",\n", | |
" {\n", | |
" center: [-78.99172577412558, -100.7896232307395],\n", | |
" crs: L.CRS.EPSG3857,\n", | |
" zoom: 10,\n", | |
" zoomControl: true,\n", | |
" preferCanvas: false,\n", | |
" }\n", | |
" );\n", | |
" L.control.scale().addTo(map_8366bb0414240e77d3544f24a5d8ceb1);\n", | |
"\n", | |
" \n", | |
"\n", | |
" \n", | |
" \n", | |
" var tile_layer_a06b64598dd95de7331cf01f01c72ffe = L.tileLayer(\n", | |
" "https://{s}.tile.openstreetmap.org/{z}/{x}/{y}.png",\n", | |
" {"attribution": "Data by \\u0026copy; \\u003ca target=\\"_blank\\" href=\\"http://openstreetmap.org\\"\\u003eOpenStreetMap\\u003c/a\\u003e, under \\u003ca target=\\"_blank\\" href=\\"http://www.openstreetmap.org/copyright\\"\\u003eODbL\\u003c/a\\u003e.", "detectRetina": false, "maxNativeZoom": 18, "maxZoom": 18, "minZoom": 0, "noWrap": false, "opacity": 1, "subdomains": "abc", "tms": false}\n", | |
" ).addTo(map_8366bb0414240e77d3544f24a5d8ceb1);\n", | |
" \n", | |
" \n", | |
" map_8366bb0414240e77d3544f24a5d8ceb1.fitBounds(\n", | |
" [[-79.00104942913823, -100.94678306984892], [-78.98240211911293, -100.63246339163007]],\n", | |
" {}\n", | |
" );\n", | |
" \n", | |
" \n", | |
" function geo_json_617164272ebf824d35c18b195035293f_styler(feature) {\n", | |
" switch(feature.id) {\n", | |
" default:\n", | |
" return {"fillOpacity": 0.5, "weight": 2};\n", | |
" }\n", | |
" }\n", | |
" function geo_json_617164272ebf824d35c18b195035293f_highlighter(feature) {\n", | |
" switch(feature.id) {\n", | |
" default:\n", | |
" return {"fillOpacity": 0.75};\n", | |
" }\n", | |
" }\n", | |
" function geo_json_617164272ebf824d35c18b195035293f_pointToLayer(feature, latlng) {\n", | |
" var opts = {"bubblingMouseEvents": true, "color": "#3388ff", "dashArray": null, "dashOffset": null, "fill": true, "fillColor": "#3388ff", "fillOpacity": 0.2, "fillRule": "evenodd", "lineCap": "round", "lineJoin": "round", "opacity": 1.0, "radius": 2, "stroke": true, "weight": 3};\n", | |
" \n", | |
" let style = geo_json_617164272ebf824d35c18b195035293f_styler(feature)\n", | |
" Object.assign(opts, style)\n", | |
" \n", | |
" return new L.CircleMarker(latlng, opts)\n", | |
" }\n", | |
"\n", | |
" function geo_json_617164272ebf824d35c18b195035293f_onEachFeature(feature, layer) {\n", | |
" layer.on({\n", | |
" mouseout: function(e) {\n", | |
" if(typeof e.target.setStyle === "function"){\n", | |
" geo_json_617164272ebf824d35c18b195035293f.resetStyle(e.target);\n", | |
" }\n", | |
" },\n", | |
" mouseover: function(e) {\n", | |
" if(typeof e.target.setStyle === "function"){\n", | |
" const highlightStyle = geo_json_617164272ebf824d35c18b195035293f_highlighter(e.target.feature)\n", | |
" e.target.setStyle(highlightStyle);\n", | |
" }\n", | |
" },\n", | |
" });\n", | |
" };\n", | |
" var geo_json_617164272ebf824d35c18b195035293f = L.geoJson(null, {\n", | |
" onEachFeature: geo_json_617164272ebf824d35c18b195035293f_onEachFeature,\n", | |
" \n", | |
" style: geo_json_617164272ebf824d35c18b195035293f_styler,\n", | |
" pointToLayer: geo_json_617164272ebf824d35c18b195035293f_pointToLayer\n", | |
" });\n", | |
"\n", | |
" function geo_json_617164272ebf824d35c18b195035293f_add (data) {\n", | |
" geo_json_617164272ebf824d35c18b195035293f\n", | |
" .addData(data)\n", | |
" .addTo(map_8366bb0414240e77d3544f24a5d8ceb1);\n", | |
" }\n", | |
" geo_json_617164272ebf824d35c18b195035293f_add({"bbox": [-100.94678306984892, -79.00104942913823, -100.63246339163007, -78.98240211911293], "features": [{"bbox": [-100.63246339163007, -78.98921540514897, -100.63246339163007, -78.98921540514897], "geometry": {"coordinates": [-100.63246339163007, -78.98921540514897], "type": "Point"}, "id": "0", "properties": {"cycle": 1, "dh_fit_dx": 0.010870090490449251, "dh_fit_dy": 0.0, "distance": 31318165.597002406, "gt": 50, "h_mean": 1875.2172759600721, "h_sigma": 0.06681861493225781, "index": 0, "n_fit_photons": 85, "pflags": 4, "rgt": 325, "rms_misfit": 0.6123660041359527, "segment_id": 1565343, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 5.361178154641397}, "type": "Feature"}, {"bbox": [-100.63264780169405, -78.98904067650072, -100.63264780169405, -78.98904067650072], "geometry": {"coordinates": [-100.63264780169405, -78.98904067650072], "type": "Point"}, "id": "1", "properties": {"cycle": 1, "dh_fit_dx": 0.0011063828145405745, "dh_fit_dy": 0.0, "distance": 31318185.49600201, "gt": 50, "h_mean": 1875.3668984219225, "h_sigma": 0.03348806441098648, "index": 1, "n_fit_photons": 64, "pflags": 0, "rgt": 325, "rms_misfit": 0.2618516925056621, "segment_id": 1565344, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 4.32712149857403}, "type": "Feature"}, {"bbox": [-100.63283317956532, -78.98886598730077, -100.63283317956532, -78.98886598730077], "geometry": {"coordinates": [-100.63283317956532, -78.98886598730077], "type": "Point"}, "id": "2", "properties": {"cycle": 1, "dh_fit_dx": 0.020058689100974306, "dh_fit_dy": 0.0, "distance": 31318205.395001765, "gt": 50, "h_mean": 1875.2522535069245, "h_sigma": 0.06875128723901988, "index": 2, "n_fit_photons": 47, "pflags": 0, "rgt": 325, "rms_misfit": 0.29810067949800356, "segment_id": 1565345, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0625393981159528}, "type": "Feature"}, {"bbox": [-100.63302050675973, -78.9886913739691, -100.63302050675973, -78.9886913739691], "geometry": {"coordinates": [-100.63302050675973, -78.9886913739691], "type": "Point"}, "id": "3", "properties": {"cycle": 1, "dh_fit_dx": 0.002709077806056982, "dh_fit_dy": 0.0, "distance": 31318225.29400167, "gt": 50, "h_mean": 1875.5041817571055, "h_sigma": 0.020000326170482034, "index": 3, "n_fit_photons": 86, "pflags": 0, "rgt": 325, "rms_misfit": 0.18224103535297453, "segment_id": 1565346, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.63320809674873, -78.988516768308, -100.63320809674873, -78.988516768308], "geometry": {"coordinates": [-100.63320809674873, -78.988516768308], "type": "Point"}, "id": "4", "properties": {"cycle": 1, "dh_fit_dx": -0.0008054574414121906, "dh_fit_dy": 0.0, "distance": 31318245.19300172, "gt": 50, "h_mean": 1875.5264100114541, "h_sigma": 0.02722500758751929, "index": 4, "n_fit_photons": 100, "pflags": 0, "rgt": 325, "rms_misfit": 0.27177678343134465, "segment_id": 1565347, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.63339585229686, -78.9883421712285, -100.63339585229686, -78.9883421712285], "geometry": {"coordinates": [-100.63339585229686, -78.9883421712285], "type": "Point"}, "id": "5", "properties": {"cycle": 1, "dh_fit_dx": 0.0001297665120547375, "dh_fit_dy": 0.0, "distance": 31318265.092001922, "gt": 50, "h_mean": 1875.564067193102, "h_sigma": 0.03173536141256648, "index": 5, "n_fit_photons": 94, "pflags": 0, "rgt": 325, "rms_misfit": 0.30602311835506507, "segment_id": 1565348, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.63358346280299, -78.98816756798207, -100.63358346280299, -78.98816756798207], "geometry": {"coordinates": [-100.63358346280299, -78.98816756798207], "type": "Point"}, "id": "6", "properties": {"cycle": 1, "dh_fit_dx": -0.0017951471206665418, "dh_fit_dy": 0.0, "distance": 31318284.991002273, "gt": 50, "h_mean": 1875.4769320515989, "h_sigma": 0.05632849352510723, "index": 6, "n_fit_photons": 86, "pflags": 4, "rgt": 325, "rms_misfit": 0.5220323142618515, "segment_id": 1565349, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.5537283289121704}, "type": "Feature"}, {"bbox": [-100.633770818648, -78.98799295503935, -100.633770818648, -78.98799295503935], "geometry": {"coordinates": [-100.633770818648, -78.98799295503935], "type": "Point"}, "id": "7", "properties": {"cycle": 1, "dh_fit_dx": 0.010707852730142896, "dh_fit_dy": 0.0, "distance": 31318304.890002772, "gt": 50, "h_mean": 1875.525429338919, "h_sigma": 0.07587017121151803, "index": 7, "n_fit_photons": 88, "pflags": 0, "rgt": 325, "rms_misfit": 0.7115993788608449, "segment_id": 1565350, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 6.490640771006121}, "type": "Feature"}, {"bbox": [-100.7871943205637, -78.99517154328777, -100.7871943205637, -78.99517154328777], "geometry": {"coordinates": [-100.7871943205637, -78.99517154328777], "type": "Point"}, "id": "8", "properties": {"cycle": 1, "dh_fit_dx": -0.02702510792175261, "dh_fit_dy": 0.0, "distance": 31318165.597002406, "gt": 30, "h_mean": 1877.151904385501, "h_sigma": 0.0, "index": 8, "n_fit_photons": 157, "pflags": 2, "rgt": 325, "rms_misfit": 0.0, "segment_id": 1565343, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.63395771337754, -78.98781832540189, -100.63395771337754, -78.98781832540189], "geometry": {"coordinates": [-100.63395771337754, -78.98781832540189], "type": "Point"}, "id": "9", "properties": {"cycle": 1, "dh_fit_dx": -0.0012187921232758919, "dh_fit_dy": 0.0, "distance": 31318324.78900342, "gt": 50, "h_mean": 1875.5895346741133, "h_sigma": 0.04438870671584975, "index": 9, "n_fit_photons": 80, "pflags": 4, "rgt": 325, "rms_misfit": 0.3968501834819791, "segment_id": 1565351, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.78737480400298, -78.99499667409697, -100.78737480400298, -78.99499667409697], "geometry": {"coordinates": [-100.78737480400298, -78.99499667409697], "type": "Point"}, "id": "10", "properties": {"cycle": 1, "dh_fit_dx": -0.02366348155912469, "dh_fit_dy": 0.0, "distance": 31318185.49600201, "gt": 30, "h_mean": 1874.180242030695, "h_sigma": 0.08808078980038779, "index": 10, "n_fit_photons": 87, "pflags": 4, "rgt": 325, "rms_misfit": 0.8177833483011702, "segment_id": 1565344, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 4.676488127248978}, "type": "Feature"}, {"bbox": [-100.6341437708202, -78.98764366207739, -100.6341437708202, -78.98764366207739], "geometry": {"coordinates": [-100.6341437708202, -78.98764366207739], "type": "Point"}, "id": "11", "properties": {"cycle": 1, "dh_fit_dx": 0.006857301337553129, "dh_fit_dy": 0.0, "distance": 31318344.688004218, "gt": 50, "h_mean": 1875.6631881334922, "h_sigma": 0.02630952391692351, "index": 11, "n_fit_photons": 79, "pflags": 0, "rgt": 325, "rms_misfit": 0.2328872576278779, "segment_id": 1565352, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.7875573009986, -78.99482188222775, -100.7875573009986, -78.99482188222775], "geometry": {"coordinates": [-100.7875573009986, -78.99482188222775], "type": "Point"}, "id": "12", "properties": {"cycle": 1, "dh_fit_dx": -0.00029498537467771187, "dh_fit_dy": 0.0, "distance": 31318205.395001765, "gt": 30, "h_mean": 1874.3258263176474, "h_sigma": 0.031169179538033054, "index": 12, "n_fit_photons": 71, "pflags": 0, "rgt": 325, "rms_misfit": 0.26197110800383255, "segment_id": 1565345, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.63432880334332, -78.98746896027504, -100.63432880334332, -78.98746896027504], "geometry": {"coordinates": [-100.63432880334332, -78.98746896027504], "type": "Point"}, "id": "13", "properties": {"cycle": 1, "dh_fit_dx": -0.004941196744140535, "dh_fit_dy": 0.0, "distance": 31318364.587005164, "gt": 50, "h_mean": 1875.825610649132, "h_sigma": 0.0319001820458606, "index": 13, "n_fit_photons": 62, "pflags": 0, "rgt": 325, "rms_misfit": 0.1953497960002299, "segment_id": 1565353, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.78774118399565, -78.9946471431624, -100.78774118399565, -78.9946471431624], "geometry": {"coordinates": [-100.78774118399565, -78.9946471431624], "type": "Point"}, "id": "14", "properties": {"cycle": 1, "dh_fit_dx": 0.0019812957768892366, "dh_fit_dy": 0.0, "distance": 31318225.29400167, "gt": 30, "h_mean": 1874.4415734970635, "h_sigma": 0.06275460628887168, "index": 14, "n_fit_photons": 72, "pflags": 4, "rgt": 325, "rms_misfit": 0.5310937040188418, "segment_id": 1565346, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.3791155465055382}, "type": "Feature"}, {"bbox": [-100.63451283135224, -78.98729421671737, -100.63451283135224, -78.98729421671737], "geometry": {"coordinates": [-100.63451283135224, -78.98729421671737], "type": "Point"}, "id": "15", "properties": {"cycle": 1, "dh_fit_dx": 0.0036175428121803044, "dh_fit_dy": 0.0, "distance": 31318384.48600626, "gt": 50, "h_mean": 1875.8015832716728, "h_sigma": 0.02403797518064921, "index": 15, "n_fit_photons": 63, "pflags": 0, "rgt": 325, "rms_misfit": 0.15723462204799432, "segment_id": 1565354, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.78792567012488, -78.99447242718041, -100.78792567012488, -78.99447242718041], "geometry": {"coordinates": [-100.78792567012488, -78.99447242718041], "type": "Point"}, "id": "16", "properties": {"cycle": 1, "dh_fit_dx": -0.0011488995301435474, "dh_fit_dy": 0.0, "distance": 31318245.19300172, "gt": 30, "h_mean": 1874.3486912236879, "h_sigma": 0.03815236458853898, "index": 16, "n_fit_photons": 66, "pflags": 0, "rgt": 325, "rms_misfit": 0.3099335342653288, "segment_id": 1565347, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.63469626166224, -78.98711945286234, -100.63469626166224, -78.98711945286234], "geometry": {"coordinates": [-100.63469626166224, -78.98711945286234], "type": "Point"}, "id": "17", "properties": {"cycle": 1, "dh_fit_dx": 0.005940338825529322, "dh_fit_dy": 0.0, "distance": 31318404.385007504, "gt": 50, "h_mean": 1875.8041803064864, "h_sigma": 0.030391611172854816, "index": 17, "n_fit_photons": 94, "pflags": 4, "rgt": 325, "rms_misfit": 0.28646023324938075, "segment_id": 1565355, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.78811042213263, -78.99429772161223, -100.78811042213263, -78.99429772161223], "geometry": {"coordinates": [-100.78811042213263, -78.99429772161223], "type": "Point"}, "id": "18", "properties": {"cycle": 1, "dh_fit_dx": -0.002710615407636506, "dh_fit_dy": 0.0, "distance": 31318265.092001922, "gt": 30, "h_mean": 1874.2438278113248, "h_sigma": 0.03430718027891223, "index": 18, "n_fit_photons": 73, "pflags": 0, "rgt": 325, "rms_misfit": 0.2872312499166884, "segment_id": 1565348, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.78829526013388, -78.99412301953608, -100.78829526013388, -78.99412301953608], "geometry": {"coordinates": [-100.78829526013388, -78.99412301953608], "type": "Point"}, "id": "19", "properties": {"cycle": 1, "dh_fit_dx": -0.0021863053041103146, "dh_fit_dy": 0.0, "distance": 31318284.991002273, "gt": 30, "h_mean": 1874.1823642164595, "h_sigma": 0.04735446489618689, "index": 19, "n_fit_photons": 77, "pflags": 0, "rgt": 325, "rms_misfit": 0.4117800537204691, "segment_id": 1565349, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 4.172236553519962}, "type": "Feature"}, {"bbox": [-100.63487941412191, -78.98694467951341, -100.63487941412191, -78.98694467951341], "geometry": {"coordinates": [-100.63487941412191, -78.98694467951341], "type": "Point"}, "id": "20", "properties": {"cycle": 1, "dh_fit_dx": -0.009999775615746536, "dh_fit_dy": 0.0, "distance": 31318424.284008898, "gt": 50, "h_mean": 1875.9775858176865, "h_sigma": 0.02788116574369504, "index": 20, "n_fit_photons": 65, "pflags": 0, "rgt": 325, "rms_misfit": 0.20407602247223217, "segment_id": 1565356, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.94209365462089, -79.00104942913823, -100.94209365462089, -79.00104942913823], "geometry": {"coordinates": [-100.94209365462089, -79.00104942913823], "type": "Point"}, "id": "21", "properties": {"cycle": 1, "dh_fit_dx": 0.06489782643660054, "dh_fit_dy": 0.0, "distance": 31318165.597002406, "gt": 10, "h_mean": 1870.5539657756367, "h_sigma": 0.7138667849245016, "index": 21, "n_fit_photons": 62, "pflags": 0, "rgt": 325, "rms_misfit": 5.572152190986255, "segment_id": 1565343, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 21.06562988708786}, "type": "Feature"}, {"bbox": [-100.78847995598332, -78.99394831161932, -100.78847995598332, -78.99394831161932], "geometry": {"coordinates": [-100.78847995598332, -78.99394831161932], "type": "Point"}, "id": "22", "properties": {"cycle": 1, "dh_fit_dx": -0.004118221980817997, "dh_fit_dy": 0.0, "distance": 31318304.890002772, "gt": 30, "h_mean": 1874.114927612472, "h_sigma": 0.05836417621953996, "index": 22, "n_fit_photons": 70, "pflags": 4, "rgt": 325, "rms_misfit": 0.4877136137483986, "segment_id": 1565350, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.970378585345956}, "type": "Feature"}, {"bbox": [-100.63506209289945, -78.98676988470486, -100.63506209289945, -78.98676988470486], "geometry": {"coordinates": [-100.63506209289945, -78.98676988470486], "type": "Point"}, "id": "23", "properties": {"cycle": 1, "dh_fit_dx": -0.0007573229390307668, "dh_fit_dy": 0.0, "distance": 31318444.18301044, "gt": 50, "h_mean": 1875.840876596566, "h_sigma": 0.023306853767456254, "index": 23, "n_fit_photons": 87, "pflags": 0, "rgt": 325, "rms_misfit": 0.2173435727046234, "segment_id": 1565357, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.94227496732512, -79.000874601023, -100.94227496732512, -79.000874601023], "geometry": {"coordinates": [-100.94227496732512, -79.000874601023], "type": "Point"}, "id": "24", "properties": {"cycle": 1, "dh_fit_dx": -0.11892988263679345, "dh_fit_dy": 0.0, "distance": 31318185.49600201, "gt": 10, "h_mean": 1869.8129360696862, "h_sigma": 0.689543215525079, "index": 24, "n_fit_photons": 88, "pflags": 0, "rgt": 325, "rms_misfit": 6.458700011007781, "segment_id": 1565344, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 29.381894307150006}, "type": "Feature"}, {"bbox": [-100.78866419118486, -78.99377358613661, -100.78866419118486, -78.99377358613661], "geometry": {"coordinates": [-100.78866419118486, -78.99377358613661], "type": "Point"}, "id": "25", "properties": {"cycle": 1, "dh_fit_dx": -0.04930973488031595, "dh_fit_dy": 0.0, "distance": 31318324.78900342, "gt": 30, "h_mean": 1873.2657837677148, "h_sigma": 0.11217788030285256, "index": 25, "n_fit_photons": 15, "pflags": 0, "rgt": 325, "rms_misfit": 0.2847612706238475, "segment_id": 1565351, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.63524482634762, -78.98659509357977, -100.63524482634762, -78.98659509357977], "geometry": {"coordinates": [-100.63524482634762, -78.98659509357977], "type": "Point"}, "id": "26", "properties": {"cycle": 1, "dh_fit_dx": 0.0010281862941789882, "dh_fit_dy": 0.0, "distance": 31318464.08201213, "gt": 50, "h_mean": 1875.7550723709128, "h_sigma": 0.04064644724894023, "index": 26, "n_fit_photons": 90, "pflags": 4, "rgt": 325, "rms_misfit": 0.38557829694368995, "segment_id": 1565358, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.94245534595221, -79.0006997376683, -100.94245534595221, -79.0006997376683], "geometry": {"coordinates": [-100.94245534595221, -79.0006997376683], "type": "Point"}, "id": "27", "properties": {"cycle": 1, "dh_fit_dx": 0.0025370215332371784, "dh_fit_dy": 0.0, "distance": 31318205.395001765, "gt": 10, "h_mean": 1871.0954725817817, "h_sigma": 0.07599974503070754, "index": 27, "n_fit_photons": 36, "pflags": 0, "rgt": 325, "rms_misfit": 0.4496833047054977, "segment_id": 1565345, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.78884821672014, -78.99359885393359, -100.78884821672014, -78.99359885393359], "geometry": {"coordinates": [-100.78884821672014, -78.99359885393359], "type": "Point"}, "id": "28", "properties": {"cycle": 1, "dh_fit_dx": 0.003282121969478211, "dh_fit_dy": 0.0, "distance": 31318344.688004218, "gt": 30, "h_mean": 1873.9942219535285, "h_sigma": 0.03127353722992504, "index": 28, "n_fit_photons": 83, "pflags": 0, "rgt": 325, "rms_misfit": 0.28370083001805146, "segment_id": 1565352, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.63542714994745, -78.98642028635697, -100.63542714994745, -78.98642028635697], "geometry": {"coordinates": [-100.63542714994745, -78.98642028635697], "type": "Point"}, "id": "29", "properties": {"cycle": 1, "dh_fit_dx": -0.0551520852500177, "dh_fit_dy": 0.0, "distance": 31318483.981013972, "gt": 50, "h_mean": 1874.991421513761, "h_sigma": 0.08836128148685711, "index": 29, "n_fit_photons": 37, "pflags": 0, "rgt": 325, "rms_misfit": 0.4049911277905729, "segment_id": 1565359, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.94263397054134, -79.00052480814938, -100.94263397054134, -79.00052480814938], "geometry": {"coordinates": [-100.94263397054134, -79.00052480814938], "type": "Point"}, "id": "30", "properties": {"cycle": 1, "dh_fit_dx": 0.024342964232492083, "dh_fit_dy": 0.0, "distance": 31318225.29400167, "gt": 10, "h_mean": 1871.019163018047, "h_sigma": 0.282960359901167, "index": 30, "n_fit_photons": 58, "pflags": 4, "rgt": 325, "rms_misfit": 2.1534579390213624, "segment_id": 1565346, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 9.4921875}, "type": "Feature"}, {"bbox": [-100.78903100596598, -78.9934240741356, -100.78903100596598, -78.9934240741356], "geometry": {"coordinates": [-100.78903100596598, -78.9934240741356], "type": "Point"}, "id": "31", "properties": {"cycle": 1, "dh_fit_dx": -0.0028271820045435597, "dh_fit_dy": 0.0, "distance": 31318364.587005164, "gt": 30, "h_mean": 1873.9894885017964, "h_sigma": 0.026503030046553043, "index": 31, "n_fit_photons": 76, "pflags": 0, "rgt": 325, "rms_misfit": 0.22904287861397776, "segment_id": 1565353, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.63561032502378, -78.98624551236507, -100.63561032502378, -78.98624551236507], "geometry": {"coordinates": [-100.63561032502378, -78.98624551236507], "type": "Point"}, "id": "32", "properties": {"cycle": 1, "dh_fit_dx": 0.01707093305775969, "dh_fit_dy": 0.0, "distance": 31318503.88001596, "gt": 50, "h_mean": 1875.7432101824697, "h_sigma": 0.04146776388137765, "index": 32, "n_fit_photons": 73, "pflags": 0, "rgt": 325, "rms_misfit": 0.32054719948321936, "segment_id": 1565360, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.94281144870938, -79.00034983568705, -100.94281144870938, -79.00034983568705], "geometry": {"coordinates": [-100.94281144870938, -79.00034983568705], "type": "Point"}, "id": "33", "properties": {"cycle": 1, "dh_fit_dx": -0.11661534059311761, "dh_fit_dy": 0.0, "distance": 31318245.19300172, "gt": 10, "h_mean": 1870.1074370381295, "h_sigma": 0.18845157267757237, "index": 33, "n_fit_photons": 33, "pflags": 0, "rgt": 325, "rms_misfit": 1.045818352563494, "segment_id": 1565347, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 7.846335398149742}, "type": "Feature"}, {"bbox": [-100.78921321465148, -78.99324927103645, -100.78921321465148, -78.99324927103645], "geometry": {"coordinates": [-100.78921321465148, -78.99324927103645], "type": "Point"}, "id": "34", "properties": {"cycle": 1, "dh_fit_dx": -0.002476850074770127, "dh_fit_dy": 0.0, "distance": 31318384.48600626, "gt": 30, "h_mean": 1873.8334831709058, "h_sigma": 0.08776367928389069, "index": 34, "n_fit_photons": 80, "pflags": 4, "rgt": 325, "rms_misfit": 0.7745330983152842, "segment_id": 1565354, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 5.0097907178544325}, "type": "Feature"}, {"bbox": [-100.63579369670691, -78.98607074733596, -100.63579369670691, -78.98607074733596], "geometry": {"coordinates": [-100.63579369670691, -78.98607074733596], "type": "Point"}, "id": "35", "properties": {"cycle": 1, "dh_fit_dx": 0.007137363227544392, "dh_fit_dy": 0.0, "distance": 31318523.7790181, "gt": 50, "h_mean": 1875.8171851265106, "h_sigma": 0.041625631345336256, "index": 35, "n_fit_photons": 80, "pflags": 0, "rgt": 325, "rms_misfit": 0.35463913341795156, "segment_id": 1565361, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.942989171836, -79.00017487257404, -100.942989171836, -79.00017487257404], "geometry": {"coordinates": [-100.942989171836, -79.00017487257404], "type": "Point"}, "id": "36", "properties": {"cycle": 1, "dh_fit_dx": -0.015471007677448653, "dh_fit_dy": 0.0, "distance": 31318265.092001922, "gt": 10, "h_mean": 1867.8170125617714, "h_sigma": 0.0, "index": 36, "n_fit_photons": 195, "pflags": 2, "rgt": 325, "rms_misfit": 0.0, "segment_id": 1565348, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.78939530019755, -78.99307446470668, -100.78939530019755, -78.99307446470668], "geometry": {"coordinates": [-100.78939530019755, -78.99307446470668], "type": "Point"}, "id": "37", "properties": {"cycle": 1, "dh_fit_dx": -0.008409875496191078, "dh_fit_dy": 0.0, "distance": 31318404.385007504, "gt": 30, "h_mean": 1873.7289580748693, "h_sigma": 0.07266829836859198, "index": 37, "n_fit_photons": 79, "pflags": 4, "rgt": 325, "rms_misfit": 0.6390425742258685, "segment_id": 1565355, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.958958090565245}, "type": "Feature"}, {"bbox": [-100.6359779931904, -78.98589601707094, -100.6359779931904, -78.98589601707094], "geometry": {"coordinates": [-100.6359779931904, -78.98589601707094], "type": "Point"}, "id": "38", "properties": {"cycle": 1, "dh_fit_dx": -0.0028571455143197983, "dh_fit_dy": 0.0, "distance": 31318543.67802039, "gt": 50, "h_mean": 1875.8550527610769, "h_sigma": 0.03787919274180721, "index": 38, "n_fit_photons": 94, "pflags": 0, "rgt": 325, "rms_misfit": 0.36700612083125095, "segment_id": 1565362, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.94316709305949, -78.99999991686799, -100.94316709305949, -78.99999991686799], "geometry": {"coordinates": [-100.94316709305949, -78.99999991686799], "type": "Point"}, "id": "39", "properties": {"cycle": 1, "dh_fit_dx": -0.10341691971591796, "dh_fit_dy": 0.0, "distance": 31318284.991002273, "gt": 10, "h_mean": 1869.9850352792416, "h_sigma": 0.20747155923406332, "index": 39, "n_fit_photons": 54, "pflags": 4, "rgt": 325, "rms_misfit": 1.1438327165046744, "segment_id": 1565349, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 6.079430986365825}, "type": "Feature"}, {"bbox": [-100.78957747929813, -78.99289966131913, -100.78957747929813, -78.99289966131913], "geometry": {"coordinates": [-100.78957747929813, -78.99289966131913], "type": "Point"}, "id": "40", "properties": {"cycle": 1, "dh_fit_dx": 0.005598364475377973, "dh_fit_dy": 0.0, "distance": 31318424.284008898, "gt": 30, "h_mean": 1873.8673152163728, "h_sigma": 0.04186702906895921, "index": 40, "n_fit_photons": 72, "pflags": 0, "rgt": 325, "rms_misfit": 0.3552199926220937, "segment_id": 1565356, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.1428385121921583}, "type": "Feature"}, {"bbox": [-100.63616323799279, -78.98572132378759, -100.63616323799279, -78.98572132378759], "geometry": {"coordinates": [-100.63616323799279, -78.98572132378759], "type": "Point"}, "id": "41", "properties": {"cycle": 1, "dh_fit_dx": 0.004567016861487032, "dh_fit_dy": 0.0, "distance": 31318563.57702283, "gt": 50, "h_mean": 1875.839722147026, "h_sigma": 0.03581468206076469, "index": 41, "n_fit_photons": 95, "pflags": 0, "rgt": 325, "rms_misfit": 0.3466546954478622, "segment_id": 1565363, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.94334659502681, -78.99982502128023, -100.94334659502681, -78.99982502128023], "geometry": {"coordinates": [-100.94334659502681, -78.99982502128023], "type": "Point"}, "id": "42", "properties": {"cycle": 1, "dh_fit_dx": -0.05954291857901705, "dh_fit_dy": 0.0, "distance": 31318304.890002772, "gt": 10, "h_mean": 1866.0808921696978, "h_sigma": 0.0, "index": 42, "n_fit_photons": 199, "pflags": 2, "rgt": 325, "rms_misfit": 0.0, "segment_id": 1565350, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 5.014720537295542}, "type": "Feature"}, {"bbox": [-100.63634961125874, -78.98554667405912, -100.63634961125874, -78.98554667405912], "geometry": {"coordinates": [-100.63634961125874, -78.98554667405912], "type": "Point"}, "id": "43", "properties": {"cycle": 1, "dh_fit_dx": 0.0002857643931784028, "dh_fit_dy": 0.0, "distance": 31318583.476025417, "gt": 50, "h_mean": 1875.908673998327, "h_sigma": 0.02263678353648522, "index": 43, "n_fit_photons": 101, "pflags": 0, "rgt": 325, "rms_misfit": 0.22710134068094762, "segment_id": 1565364, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.78975983206166, -78.99272486488178, -100.78975983206166, -78.99272486488178], "geometry": {"coordinates": [-100.78975983206166, -78.99272486488178], "type": "Point"}, "id": "44", "properties": {"cycle": 1, "dh_fit_dx": 0.010844705527191287, "dh_fit_dy": 0.0, "distance": 31318444.18301044, "gt": 30, "h_mean": 1873.9199881475724, "h_sigma": 0.05844579099143813, "index": 44, "n_fit_photons": 87, "pflags": 0, "rgt": 325, "rms_misfit": 0.5447669063974334, "segment_id": 1565357, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 6.122397446540845}, "type": "Feature"}, {"bbox": [-100.94352688338655, -78.99965015513469, -100.94352688338655, -78.99965015513469], "geometry": {"coordinates": [-100.94352688338655, -78.99965015513469], "type": "Point"}, "id": "45", "properties": {"cycle": 1, "dh_fit_dx": 0.03216445383142637, "dh_fit_dy": 0.0, "distance": 31318324.78900342, "gt": 10, "h_mean": 1866.8024807501063, "h_sigma": 0.0, "index": 45, "n_fit_photons": 179, "pflags": 2, "rgt": 325, "rms_misfit": 0.0, "segment_id": 1565351, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.63653687768243, -78.98537205927262, -100.63653687768243, -78.98537205927262], "geometry": {"coordinates": [-100.63653687768243, -78.98537205927262], "type": "Point"}, "id": "46", "properties": {"cycle": 1, "dh_fit_dx": -0.002127145761420209, "dh_fit_dy": 0.0, "distance": 31318603.375028156, "gt": 50, "h_mean": 1875.8399196060107, "h_sigma": 0.032368385024775986, "index": 46, "n_fit_photons": 96, "pflags": 0, "rgt": 325, "rms_misfit": 0.31496666329100487, "segment_id": 1565365, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.78994244515671, -78.99255007903287, -100.78994244515671, -78.99255007903287], "geometry": {"coordinates": [-100.78994244515671, -78.99255007903287], "type": "Point"}, "id": "47", "properties": {"cycle": 1, "dh_fit_dx": -0.0012754978616258486, "dh_fit_dy": 0.0, "distance": 31318464.08201213, "gt": 30, "h_mean": 1873.9612452287138, "h_sigma": 0.054727856630136594, "index": 47, "n_fit_photons": 89, "pflags": 4, "rgt": 325, "rms_misfit": 0.5161662274383139, "segment_id": 1565358, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.685583392251897}, "type": "Feature"}, {"bbox": [-100.9437081253888, -78.99947532501601, -100.9437081253888, -78.99947532501601], "geometry": {"coordinates": [-100.9437081253888, -78.99947532501601], "type": "Point"}, "id": "48", "properties": {"cycle": 1, "dh_fit_dx": -0.0031012578411041595, "dh_fit_dy": 0.0, "distance": 31318344.688004218, "gt": 10, "h_mean": 1870.9365366137642, "h_sigma": 0.036032485447342624, "index": 48, "n_fit_photons": 53, "pflags": 0, "rgt": 325, "rms_misfit": 0.259545849404715, "segment_id": 1565352, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.63672560494659, -78.9851975009372, -100.63672560494659, -78.9851975009372], "geometry": {"coordinates": [-100.63672560494659, -78.9851975009372], "type": "Point"}, "id": "49", "properties": {"cycle": 1, "dh_fit_dx": 0.0037611397175001835, "dh_fit_dy": 0.0, "distance": 31318623.274031043, "gt": 50, "h_mean": 1875.8064164133948, "h_sigma": 0.10505394360064503, "index": 49, "n_fit_photons": 90, "pflags": 0, "rgt": 325, "rms_misfit": 0.9944381532063316, "segment_id": 1565366, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 9.45345621944209}, "type": "Feature"}, {"bbox": [-100.79012591915667, -78.99237532754324, -100.79012591915667, -78.99237532754324], "geometry": {"coordinates": [-100.79012591915667, -78.99237532754324], "type": "Point"}, "id": "50", "properties": {"cycle": 1, "dh_fit_dx": 0.11957927706894707, "dh_fit_dy": 0.0, "distance": 31318483.981013972, "gt": 30, "h_mean": 1875.9022707133909, "h_sigma": 0.3015615093505047, "index": 50, "n_fit_photons": 27, "pflags": 0, "rgt": 325, "rms_misfit": 0.8562231842107929, "segment_id": 1565359, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 4.4685388933939}, "type": "Feature"}, {"bbox": [-100.94389248916063, -78.99930061219007, -100.94389248916063, -78.99930061219007], "geometry": {"coordinates": [-100.94389248916063, -78.99930061219007], "type": "Point"}, "id": "51", "properties": {"cycle": 1, "dh_fit_dx": 0.09649452411249727, "dh_fit_dy": 0.0, "distance": 31318364.587005164, "gt": 10, "h_mean": 1872.205041285209, "h_sigma": 0.13441537345268678, "index": 51, "n_fit_photons": 35, "pflags": 0, "rgt": 325, "rms_misfit": 0.5988306136246216, "segment_id": 1565353, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 3.527910759612125}, "type": "Feature"}, {"bbox": [-100.63691577715035, -78.98502299820642, -100.63691577715035, -78.98502299820642], "geometry": {"coordinates": [-100.63691577715035, -78.98502299820642], "type": "Point"}, "id": "52", "properties": {"cycle": 1, "dh_fit_dx": -0.04328354425021256, "dh_fit_dy": 0.0, "distance": 31318643.17303408, "gt": 50, "h_mean": 1875.1817977142725, "h_sigma": 0.0732778870932537, "index": 52, "n_fit_photons": 23, "pflags": 0, "rgt": 325, "rms_misfit": 0.21252272938650826, "segment_id": 1565367, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.79030948299747, -78.99220057657068, -100.79030948299747, -78.99220057657068], "geometry": {"coordinates": [-100.79030948299747, -78.99220057657068], "type": "Point"}, "id": "53", "properties": {"cycle": 1, "dh_fit_dx": -0.012445771592666629, "dh_fit_dy": 0.0, "distance": 31318503.88001596, "gt": 30, "h_mean": 1873.8485191352036, "h_sigma": 0.04037026761982368, "index": 53, "n_fit_photons": 60, "pflags": 0, "rgt": 325, "rms_misfit": 0.28730599756080943, "segment_id": 1565360, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.53685528734785}, "type": "Feature"}, {"bbox": [-100.94407926790902, -78.9991259912435, -100.94407926790902, -78.9991259912435], "geometry": {"coordinates": [-100.94407926790902, -78.9991259912435], "type": "Point"}, "id": "54", "properties": {"cycle": 1, "dh_fit_dx": 0.019585295760252596, "dh_fit_dy": 0.0, "distance": 31318384.48600626, "gt": 10, "h_mean": 1871.1506193713833, "h_sigma": 0.07032039759191158, "index": 54, "n_fit_photons": 32, "pflags": 0, "rgt": 325, "rms_misfit": 0.3202554725877171, "segment_id": 1565354, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.63710739707766, -78.98484855408772, -100.63710739707766, -78.98484855408772], "geometry": {"coordinates": [-100.63710739707766, -78.98484855408772], "type": "Point"}, "id": "55", "properties": {"cycle": 1, "dh_fit_dx": 0.0018878282534452453, "dh_fit_dy": 0.0, "distance": 31318663.072037265, "gt": 50, "h_mean": 1875.9114521180265, "h_sigma": 0.036757150557396916, "index": 55, "n_fit_photons": 84, "pflags": 4, "rgt": 325, "rms_misfit": 0.3364453508856395, "segment_id": 1565368, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.79049450773643, -78.99202588237486, -100.79049450773643, -78.99202588237486], "geometry": {"coordinates": [-100.79049450773643, -78.99202588237486], "type": "Point"}, "id": "56", "properties": {"cycle": 1, "dh_fit_dx": 0.03196625074798625, "dh_fit_dy": 0.0, "distance": 31318523.7790181, "gt": 30, "h_mean": 1876.103772148735, "h_sigma": 0.0, "index": 56, "n_fit_photons": 142, "pflags": 2, "rgt": 325, "rms_misfit": 0.0, "segment_id": 1565361, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.9442663925819, -78.9989513832374, -100.9442663925819, -78.9989513832374], "geometry": {"coordinates": [-100.9442663925819, -78.9989513832374], "type": "Point"}, "id": "57", "properties": {"cycle": 1, "dh_fit_dx": -0.0024252795841737712, "dh_fit_dy": 0.0, "distance": 31318404.385007504, "gt": 10, "h_mean": 1870.7437341815887, "h_sigma": 0.07222830736841135, "index": 57, "n_fit_photons": 54, "pflags": 4, "rgt": 325, "rms_misfit": 0.5203456675612699, "segment_id": 1565355, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.63729826448628, -78.98467407918245, -100.63729826448628, -78.98467407918245], "geometry": {"coordinates": [-100.63729826448628, -78.98467407918245], "type": "Point"}, "id": "58", "properties": {"cycle": 1, "dh_fit_dx": 0.016588558633418123, "dh_fit_dy": 0.0, "distance": 31318682.9710406, "gt": 50, "h_mean": 1875.9780448009851, "h_sigma": 0.08304162251533909, "index": 58, "n_fit_photons": 95, "pflags": 4, "rgt": 325, "rms_misfit": 0.8093252747544806, "segment_id": 1565369, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 5.011445263428664}, "type": "Feature"}, {"bbox": [-100.79067910110065, -78.99185117152926, -100.79067910110065, -78.99185117152926], "geometry": {"coordinates": [-100.79067910110065, -78.99185117152926], "type": "Point"}, "id": "59", "properties": {"cycle": 1, "dh_fit_dx": -0.11697492242663, "dh_fit_dy": 0.0, "distance": 31318543.67802039, "gt": 30, "h_mean": 1875.908555937835, "h_sigma": 0.2043880760333796, "index": 59, "n_fit_photons": 12, "pflags": 0, "rgt": 325, "rms_misfit": 0.34312810977377733, "segment_id": 1565362, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.13912337625456}, "type": "Feature"}, {"bbox": [-100.94445335703321, -78.99877676968856, -100.94445335703321, -78.99877676968856], "geometry": {"coordinates": [-100.94445335703321, -78.99877676968856], "type": "Point"}, "id": "60", "properties": {"cycle": 1, "dh_fit_dx": 0.031748444022672725, "dh_fit_dy": 0.0, "distance": 31318424.284008898, "gt": 10, "h_mean": 1870.388295245465, "h_sigma": 0.10074885426823686, "index": 60, "n_fit_photons": 37, "pflags": 0, "rgt": 325, "rms_misfit": 0.4163886693243068, "segment_id": 1565356, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.7908633096352, -78.99167644755238, -100.7908633096352, -78.99167644755238], "geometry": {"coordinates": [-100.7908633096352, -78.99167644755238], "type": "Point"}, "id": "61", "properties": {"cycle": 1, "dh_fit_dx": 0.02552290117845439, "dh_fit_dy": 0.0, "distance": 31318563.57702283, "gt": 30, "h_mean": 1873.7801869139018, "h_sigma": 0.09957487510744628, "index": 61, "n_fit_photons": 91, "pflags": 4, "rgt": 325, "rms_misfit": 0.9483944242678065, "segment_id": 1565363, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 6.121030148052542}, "type": "Feature"}, {"bbox": [-100.63748824926877, -78.98449956985039, -100.63748824926877, -78.98449956985039], "geometry": {"coordinates": [-100.63748824926877, -78.98449956985039], "type": "Point"}, "id": "62", "properties": {"cycle": 1, "dh_fit_dx": 7.635820462387591e-05, "dh_fit_dy": 0.0, "distance": 31318702.870044082, "gt": 50, "h_mean": 1876.0678504624484, "h_sigma": 0.05982414892123962, "index": 62, "n_fit_photons": 87, "pflags": 4, "rgt": 325, "rms_misfit": 0.5502270161270527, "segment_id": 1565370, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 4.49166926459711}, "type": "Feature"}, {"bbox": [-100.94463929671117, -78.998602116338, -100.94463929671117, -78.998602116338], "geometry": {"coordinates": [-100.94463929671117, -78.998602116338], "type": "Point"}, "id": "63", "properties": {"cycle": 1, "dh_fit_dx": 0.008153498127616851, "dh_fit_dy": 0.0, "distance": 31318444.18301044, "gt": 10, "h_mean": 1870.574906615831, "h_sigma": 0.07194001605440659, "index": 63, "n_fit_photons": 66, "pflags": 4, "rgt": 325, "rms_misfit": 0.5776529134400445, "segment_id": 1565357, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 3.570545993276083}, "type": "Feature"}, {"bbox": [-100.79104800391258, -78.99150174132265, -100.79104800391258, -78.99150174132265], "geometry": {"coordinates": [-100.79104800391258, -78.99150174132265], "type": "Point"}, "id": "64", "properties": {"cycle": 1, "dh_fit_dx": -0.006050729308555293, "dh_fit_dy": 0.0, "distance": 31318583.476025417, "gt": 30, "h_mean": 1873.7108019631228, "h_sigma": 0.0349283665605208, "index": 64, "n_fit_photons": 72, "pflags": 0, "rgt": 325, "rms_misfit": 0.29629030547735324, "segment_id": 1565364, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.63767680398878, -78.98432500502778, -100.63767680398878, -78.98432500502778], "geometry": {"coordinates": [-100.63767680398878, -78.98432500502778], "type": "Point"}, "id": "65", "properties": {"cycle": 1, "dh_fit_dx": -0.00166783191201688, "dh_fit_dy": 0.0, "distance": 31318722.769047715, "gt": 50, "h_mean": 1875.9335259979853, "h_sigma": 0.03266338888364193, "index": 65, "n_fit_photons": 79, "pflags": 0, "rgt": 325, "rms_misfit": 0.2903174053005659, "segment_id": 1565371, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.94482404682856, -78.99842741896596, -100.94482404682856, -78.99842741896596], "geometry": {"coordinates": [-100.94482404682856, -78.99842741896596], "type": "Point"}, "id": "66", "properties": {"cycle": 1, "dh_fit_dx": 0.009114998485159587, "dh_fit_dy": 0.0, "distance": 31318464.08201213, "gt": 10, "h_mean": 1870.6862339231263, "h_sigma": 0.05748924168245941, "index": 66, "n_fit_photons": 64, "pflags": 4, "rgt": 325, "rms_misfit": 0.4569204930962826, "segment_id": 1565358, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 3.7983619785808536}, "type": "Feature"}, {"bbox": [-100.79123246311664, -78.99132702610947, -100.79123246311664, -78.99132702610947], "geometry": {"coordinates": [-100.79123246311664, -78.99132702610947], "type": "Point"}, "id": "67", "properties": {"cycle": 1, "dh_fit_dx": -0.005752119961143618, "dh_fit_dy": 0.0, "distance": 31318603.375028156, "gt": 30, "h_mean": 1873.532934949737, "h_sigma": 0.05345940748454325, "index": 67, "n_fit_photons": 82, "pflags": 0, "rgt": 325, "rms_misfit": 0.48349691971466663, "segment_id": 1565365, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 6.535347388911264}, "type": "Feature"}, {"bbox": [-100.63786398107916, -78.98415038928214, -100.63786398107916, -78.98415038928214], "geometry": {"coordinates": [-100.63786398107916, -78.98415038928214], "type": "Point"}, "id": "68", "properties": {"cycle": 1, "dh_fit_dx": 0.003829363619506022, "dh_fit_dy": 0.0, "distance": 31318742.668051496, "gt": 50, "h_mean": 1875.8078850679387, "h_sigma": 0.1104363994118082, "index": 68, "n_fit_photons": 92, "pflags": 4, "rgt": 325, "rms_misfit": 1.0569492354946115, "segment_id": 1565372, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 8.071697579018593}, "type": "Feature"}, {"bbox": [-100.94500811141845, -78.99825269638282, -100.94500811141845, -78.99825269638282], "geometry": {"coordinates": [-100.94500811141845, -78.99825269638282], "type": "Point"}, "id": "69", "properties": {"cycle": 1, "dh_fit_dx": 0.09054971795456165, "dh_fit_dy": 0.0, "distance": 31318483.981013972, "gt": 10, "h_mean": 1869.2936646070868, "h_sigma": 0.15531369348852814, "index": 69, "n_fit_photons": 37, "pflags": 0, "rgt": 325, "rms_misfit": 0.5329911109120381, "segment_id": 1565359, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 3.4469129750127534}, "type": "Feature"}, {"bbox": [-100.7914166656317, -78.99115230153978, -100.7914166656317, -78.99115230153978], "geometry": {"coordinates": [-100.7914166656317, -78.99115230153978], "type": "Point"}, "id": "70", "properties": {"cycle": 1, "dh_fit_dx": -0.001978709916820555, "dh_fit_dy": 0.0, "distance": 31318623.274031043, "gt": 30, "h_mean": 1873.5241716078626, "h_sigma": 0.05038792372031999, "index": 70, "n_fit_photons": 83, "pflags": 0, "rgt": 325, "rms_misfit": 0.4562162591352633, "segment_id": 1565366, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.906667584474688}, "type": "Feature"}, {"bbox": [-100.63804783471791, -78.98397564397, -100.63804783471791, -78.98397564397], "geometry": {"coordinates": [-100.63804783471791, -78.98397564397], "type": "Point"}, "id": "71", "properties": {"cycle": 1, "dh_fit_dx": 0.04310239626925391, "dh_fit_dy": 0.0, "distance": 31318762.567055427, "gt": 50, "h_mean": 1873.7419069430414, "h_sigma": 0.0, "index": 71, "n_fit_photons": 146, "pflags": 2, "rgt": 325, "rms_misfit": 0.0, "segment_id": 1565373, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.060244224387975}, "type": "Feature"}, {"bbox": [-100.94518783811398, -78.99807780956287, -100.94518783811398, -78.99807780956287], "geometry": {"coordinates": [-100.94518783811398, -78.99807780956287], "type": "Point"}, "id": "72", "properties": {"cycle": 1, "dh_fit_dx": 0.006890036083376788, "dh_fit_dy": 0.0, "distance": 31318503.88001596, "gt": 10, "h_mean": 1870.5725796590023, "h_sigma": 0.05278928460648562, "index": 72, "n_fit_photons": 74, "pflags": 4, "rgt": 325, "rms_misfit": 0.44504654105662395, "segment_id": 1565360, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 3.4021275648330245}, "type": "Feature"}, {"bbox": [-100.7916005274896, -78.99097756409802, -100.7916005274896, -78.99097756409802], "geometry": {"coordinates": [-100.7916005274896, -78.99097756409802], "type": "Point"}, "id": "73", "properties": {"cycle": 1, "dh_fit_dx": -0.0021585747856443208, "dh_fit_dy": 0.0, "distance": 31318643.17303408, "gt": 30, "h_mean": 1873.458947557445, "h_sigma": 0.029768318340703398, "index": 73, "n_fit_photons": 75, "pflags": 0, "rgt": 325, "rms_misfit": 0.257223163709322, "segment_id": 1565367, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.63823090024293, -78.98380086621383, -100.63823090024293, -78.98380086621383], "geometry": {"coordinates": [-100.63823090024293, -78.98380086621383], "type": "Point"}, "id": "74", "properties": {"cycle": 1, "dh_fit_dx": 0.004387620659518638, "dh_fit_dy": 0.0, "distance": 31318782.466059506, "gt": 50, "h_mean": 1875.8988284853817, "h_sigma": 0.031557304997857846, "index": 74, "n_fit_photons": 75, "pflags": 0, "rgt": 325, "rms_misfit": 0.2700658488003202, "segment_id": 1565374, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.94536479254263, -78.99790281886217, -100.94536479254263, -78.99790281886217], "geometry": {"coordinates": [-100.94536479254263, -78.99790281886217], "type": "Point"}, "id": "75", "properties": {"cycle": 1, "dh_fit_dx": -0.03147509394814602, "dh_fit_dy": 0.0, "distance": 31318523.7790181, "gt": 10, "h_mean": 1873.1679706989605, "h_sigma": 0.0, "index": 75, "n_fit_photons": 193, "pflags": 2, "rgt": 325, "rms_misfit": 0.0, "segment_id": 1565361, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.79178388691884, -78.99080280770679, -100.79178388691884, -78.99080280770679], "geometry": {"coordinates": [-100.79178388691884, -78.99080280770679], "type": "Point"}, "id": "76", "properties": {"cycle": 1, "dh_fit_dx": -0.0004629755265847285, "dh_fit_dy": 0.0, "distance": 31318663.072037265, "gt": 30, "h_mean": 1873.389419169239, "h_sigma": 0.03576189274704694, "index": 76, "n_fit_photons": 81, "pflags": 0, "rgt": 325, "rms_misfit": 0.3177509743692562, "segment_id": 1565368, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.63841092802602, -78.98362597299744, -100.63841092802602, -78.98362597299744], "geometry": {"coordinates": [-100.63841092802602, -78.98362597299744], "type": "Point"}, "id": "77", "properties": {"cycle": 1, "dh_fit_dx": 0.0029848899929703727, "dh_fit_dy": 0.0, "distance": 31318802.365063734, "gt": 50, "h_mean": 1875.9642488556465, "h_sigma": 0.01828130010359642, "index": 77, "n_fit_photons": 83, "pflags": 0, "rgt": 325, "rms_misfit": 0.1664700069656454, "segment_id": 1565375, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.94554108789045, -78.9977278026363, -100.94554108789045, -78.9977278026363], "geometry": {"coordinates": [-100.94554108789045, -78.9977278026363], "type": "Point"}, "id": "78", "properties": {"cycle": 1, "dh_fit_dx": 0.052572798125603934, "dh_fit_dy": 0.0, "distance": 31318543.67802039, "gt": 10, "h_mean": 1871.1686367072175, "h_sigma": 0.0, "index": 78, "n_fit_photons": 29, "pflags": 1, "rgt": 325, "rms_misfit": 0.0, "segment_id": 1565362, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 3.1969602269424717}, "type": "Feature"}, {"bbox": [-100.79196648568356, -78.99062802215549, -100.79196648568356, -78.99062802215549], "geometry": {"coordinates": [-100.79196648568356, -78.99062802215549], "type": "Point"}, "id": "79", "properties": {"cycle": 1, "dh_fit_dx": 0.00233154516861674, "dh_fit_dy": 0.0, "distance": 31318682.9710406, "gt": 30, "h_mean": 1873.3646428218324, "h_sigma": 0.03710459759697346, "index": 79, "n_fit_photons": 67, "pflags": 0, "rgt": 325, "rms_misfit": 0.2943828509823545, "segment_id": 1565369, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.63858869998279, -78.98345099078989, -100.63858869998279, -78.98345099078989], "geometry": {"coordinates": [-100.63858869998279, -78.98345099078989], "type": "Point"}, "id": "80", "properties": {"cycle": 1, "dh_fit_dx": 0.0012462814109295947, "dh_fit_dy": 0.0, "distance": 31318822.26406811, "gt": 50, "h_mean": 1873.4207017470897, "h_sigma": 0.0, "index": 80, "n_fit_photons": 172, "pflags": 2, "rgt": 325, "rms_misfit": 0.0, "segment_id": 1565376, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.94571634598397, -78.99755274909245, -100.94571634598397, -78.99755274909245], "geometry": {"coordinates": [-100.94571634598397, -78.99755274909245], "type": "Point"}, "id": "81", "properties": {"cycle": 1, "dh_fit_dx": 0.14177031117903285, "dh_fit_dy": 0.0, "distance": 31318563.57702283, "gt": 10, "h_mean": 1872.8896521179129, "h_sigma": 0.2749271847349587, "index": 81, "n_fit_photons": 20, "pflags": 0, "rgt": 325, "rms_misfit": 1.0010259431403379, "segment_id": 1565363, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 17.07078843051403}, "type": "Feature"}, {"bbox": [-100.79214869065842, -78.99045322121317, -100.79214869065842, -78.99045322121317], "geometry": {"coordinates": [-100.79214869065842, -78.99045322121317], "type": "Point"}, "id": "82", "properties": {"cycle": 1, "dh_fit_dx": -0.0014982825136973812, "dh_fit_dy": 0.0, "distance": 31318702.870044082, "gt": 30, "h_mean": 1875.0125260962689, "h_sigma": 0.0, "index": 82, "n_fit_photons": 172, "pflags": 2, "rgt": 325, "rms_misfit": 0.0, "segment_id": 1565370, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.63876928019641, -78.98327611754638, -100.63876928019641, -78.98327611754638], "geometry": {"coordinates": [-100.63876928019641, -78.98327611754638], "type": "Point"}, "id": "83", "properties": {"cycle": 1, "dh_fit_dx": -0.0013734875445921351, "dh_fit_dy": 0.0, "distance": 31318842.16307264, "gt": 50, "h_mean": 1875.8800519167144, "h_sigma": 0.03839800524852741, "index": 83, "n_fit_photons": 98, "pflags": 4, "rgt": 325, "rms_misfit": 0.37720731370445776, "segment_id": 1565377, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.6864178991526506}, "type": "Feature"}, {"bbox": [-100.94589172139008, -78.9973776988299, -100.94589172139008, -78.9973776988299], "geometry": {"coordinates": [-100.94589172139008, -78.9973776988299], "type": "Point"}, "id": "84", "properties": {"cycle": 1, "dh_fit_dx": 0.13691795021679432, "dh_fit_dy": 0.0, "distance": 31318583.476025417, "gt": 10, "h_mean": 1872.4755011366858, "h_sigma": 0.28232263534572305, "index": 84, "n_fit_photons": 32, "pflags": 0, "rgt": 325, "rms_misfit": 1.290109364542204, "segment_id": 1565364, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 8.267331515360313}, "type": "Feature"}, {"bbox": [-100.79232918689905, -78.99027835515245, -100.79232918689905, -78.99027835515245], "geometry": {"coordinates": [-100.79232918689905, -78.99027835515245], "type": "Point"}, "id": "85", "properties": {"cycle": 1, "dh_fit_dx": -0.0011798297675438356, "dh_fit_dy": 0.0, "distance": 31318722.769047715, "gt": 30, "h_mean": 1873.2694177540454, "h_sigma": 0.06356799233426397, "index": 85, "n_fit_photons": 81, "pflags": 4, "rgt": 325, "rms_misfit": 0.5662471512113643, "segment_id": 1565371, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.4958448711889103}, "type": "Feature"}, {"bbox": [-100.63894887454863, -78.98310120655609, -100.63894887454863, -78.98310120655609], "geometry": {"coordinates": [-100.63894887454863, -78.98310120655609], "type": "Point"}, "id": "86", "properties": {"cycle": 1, "dh_fit_dx": 0.0028696703607674207, "dh_fit_dy": 0.0, "distance": 31318862.062077314, "gt": 50, "h_mean": 1875.8836573933709, "h_sigma": 0.027184544768953958, "index": 86, "n_fit_photons": 91, "pflags": 0, "rgt": 325, "rms_misfit": 0.2591678786168654, "segment_id": 1565378, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.94606749886299, -78.99720266442357, -100.94606749886299, -78.99720266442357], "geometry": {"coordinates": [-100.94606749886299, -78.99720266442357], "type": "Point"}, "id": "87", "properties": {"cycle": 1, "dh_fit_dx": -0.020887477702237887, "dh_fit_dy": 0.0, "distance": 31318603.375028156, "gt": 10, "h_mean": 1870.8459017480823, "h_sigma": 0.2995837182946662, "index": 87, "n_fit_photons": 56, "pflags": 4, "rgt": 325, "rms_misfit": 2.241711310845137, "segment_id": 1565365, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 9.4921875}, "type": "Feature"}, {"bbox": [-100.79250944084454, -78.9901034802122, -100.79250944084454, -78.9901034802122], "geometry": {"coordinates": [-100.79250944084454, -78.9901034802122], "type": "Point"}, "id": "88", "properties": {"cycle": 1, "dh_fit_dx": -0.005225478313106019, "dh_fit_dy": 0.0, "distance": 31318742.668051496, "gt": 30, "h_mean": 1873.2298704413781, "h_sigma": 0.024102234224253578, "index": 88, "n_fit_photons": 93, "pflags": 0, "rgt": 325, "rms_misfit": 0.2317817143661702, "segment_id": 1565372, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.63912916015498, -78.9829263232985, -100.63912916015498, -78.9829263232985], "geometry": {"coordinates": [-100.63912916015498, -78.9829263232985], "type": "Point"}, "id": "89", "properties": {"cycle": 1, "dh_fit_dx": -0.0035630456074394854, "dh_fit_dy": 0.0, "distance": 31318881.961082138, "gt": 50, "h_mean": 1875.9371177535304, "h_sigma": 0.04886434931134959, "index": 89, "n_fit_photons": 88, "pflags": 4, "rgt": 325, "rms_misfit": 0.45759760663692384, "segment_id": 1565379, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0927920841038086}, "type": "Feature"}, {"bbox": [-100.94624356992408, -78.9970276410372, -100.94624356992408, -78.9970276410372], "geometry": {"coordinates": [-100.94624356992408, -78.9970276410372], "type": "Point"}, "id": "90", "properties": {"cycle": 1, "dh_fit_dx": -0.16593783554023117, "dh_fit_dy": 0.0, "distance": 31318623.274031043, "gt": 10, "h_mean": 1871.0760209187074, "h_sigma": 0.40685201684858996, "index": 90, "n_fit_photons": 65, "pflags": 0, "rgt": 325, "rms_misfit": 3.273179431838088, "segment_id": 1565366, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 15.916794462464775}, "type": "Feature"}, {"bbox": [-100.79268881015877, -78.98992857134729, -100.79268881015877, -78.98992857134729], "geometry": {"coordinates": [-100.79268881015877, -78.98992857134729], "type": "Point"}, "id": "91", "properties": {"cycle": 1, "dh_fit_dx": -0.006488222276581484, "dh_fit_dy": 0.0, "distance": 31318762.567055427, "gt": 30, "h_mean": 1873.1260988093325, "h_sigma": 0.03348158220034802, "index": 91, "n_fit_photons": 93, "pflags": 0, "rgt": 325, "rms_misfit": 0.32216706066672113, "segment_id": 1565373, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 4.446523733842005}, "type": "Feature"}, {"bbox": [-100.63931101876791, -78.98275150100493, -100.63931101876791, -78.98275150100493], "geometry": {"coordinates": [-100.63931101876791, -78.98275150100493], "type": "Point"}, "id": "92", "properties": {"cycle": 1, "dh_fit_dx": -0.002420611823758191, "dh_fit_dy": 0.0, "distance": 31318901.86008711, "gt": 50, "h_mean": 1875.8483576044264, "h_sigma": 0.06421419999860602, "index": 92, "n_fit_photons": 87, "pflags": 4, "rgt": 325, "rms_misfit": 0.5988319038113967, "segment_id": 1565380, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.8997472350947904}, "type": "Feature"}, {"bbox": [-100.9464215815247, -78.99685269067514, -100.9464215815247, -78.99685269067514], "geometry": {"coordinates": [-100.9464215815247, -78.99685269067514], "type": "Point"}, "id": "93", "properties": {"cycle": 1, "dh_fit_dx": -0.002317418735671036, "dh_fit_dy": 0.0, "distance": 31318643.17303408, "gt": 10, "h_mean": 1870.2783495656931, "h_sigma": 0.04112922469482952, "index": 93, "n_fit_photons": 40, "pflags": 0, "rgt": 325, "rms_misfit": 0.2596849972865863, "segment_id": 1565367, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.79286687132232, -78.98975361263965, -100.79286687132232, -78.98975361263965], "geometry": {"coordinates": [-100.79286687132232, -78.98975361263965], "type": "Point"}, "id": "94", "properties": {"cycle": 1, "dh_fit_dx": -0.0002734077418122724, "dh_fit_dy": 0.0, "distance": 31318782.466059506, "gt": 30, "h_mean": 1873.0809181797488, "h_sigma": 0.038546619593472574, "index": 94, "n_fit_photons": 93, "pflags": 4, "rgt": 325, "rms_misfit": 0.3716872948213734, "segment_id": 1565374, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.63949519685809, -78.98257676860815, -100.63949519685809, -78.98257676860815], "geometry": {"coordinates": [-100.63949519685809, -78.98257676860815], "type": "Point"}, "id": "95", "properties": {"cycle": 1, "dh_fit_dx": -0.006443081948132504, "dh_fit_dy": 0.0, "distance": 31318921.759092234, "gt": 50, "h_mean": 1875.8117113884512, "h_sigma": 0.03171404330804392, "index": 95, "n_fit_photons": 88, "pflags": 0, "rgt": 325, "rms_misfit": 0.29729394948791193, "segment_id": 1565381, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.94660208212215, -78.99667783389975, -100.94660208212215, -78.99667783389975], "geometry": {"coordinates": [-100.94660208212215, -78.99667783389975], "type": "Point"}, "id": "96", "properties": {"cycle": 1, "dh_fit_dx": 0.062059165892559436, "dh_fit_dy": 0.0, "distance": 31318663.072037265, "gt": 10, "h_mean": 1871.0069817810402, "h_sigma": 0.13589371159410515, "index": 96, "n_fit_photons": 25, "pflags": 0, "rgt": 325, "rms_misfit": 0.5010211217118069, "segment_id": 1565368, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 3.4493517012761963}, "type": "Feature"}, {"bbox": [-100.7930448654226, -78.98957865183182, -100.7930448654226, -78.98957865183182], "geometry": {"coordinates": [-100.7930448654226, -78.98957865183182], "type": "Point"}, "id": "97", "properties": {"cycle": 1, "dh_fit_dx": 0.006728857799350986, "dh_fit_dy": 0.0, "distance": 31318802.365063734, "gt": 30, "h_mean": 1875.209336129865, "h_sigma": 0.0, "index": 97, "n_fit_photons": 200, "pflags": 2, "rgt": 325, "rms_misfit": 0.0, "segment_id": 1565375, "spot": 4, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.6396815562078, -78.98240211911293, -100.6396815562078, -78.98240211911293], "geometry": {"coordinates": [-100.6396815562078, -78.98240211911293], "type": "Point"}, "id": "98", "properties": {"cycle": 1, "dh_fit_dx": -0.0017539921477270407, "dh_fit_dy": 0.0, "distance": 31318941.658097506, "gt": 50, "h_mean": 1875.7605125634502, "h_sigma": 0.030065908325334036, "index": 98, "n_fit_photons": 91, "pflags": 0, "rgt": 325, "rms_misfit": 0.28650359319263563, "segment_id": 1565382, "spot": 2, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}, {"bbox": [-100.94678306984892, -78.99650299624977, -100.94678306984892, -78.99650299624977], "geometry": {"coordinates": [-100.94678306984892, -78.99650299624977], "type": "Point"}, "id": "99", "properties": {"cycle": 1, "dh_fit_dx": -0.03415061432999793, "dh_fit_dy": 0.0, "distance": 31318682.9710406, "gt": 10, "h_mean": 1870.229100068489, "h_sigma": 0.09669085009589629, "index": 99, "n_fit_photons": 28, "pflags": 4, "rgt": 325, "rms_misfit": 0.5036232303623501, "segment_id": 1565369, "spot": 6, "time": "2018-10-19", "w_surface_window_final": 3.0}, "type": "Feature"}], "type": "FeatureCollection"});\n", | |
"\n", | |
" \n", | |
" \n", | |
" geo_json_617164272ebf824d35c18b195035293f.bindTooltip(\n", | |
" function(layer){\n", | |
" let div = L.DomUtil.create('div');\n", | |
" \n", | |
" let handleObject = feature=>typeof(feature)=='object' ? JSON.stringify(feature) : feature;\n", | |
" let fields = ["index", "time", "rgt", "rms_misfit", "dh_fit_dy", "pflags", "h_sigma", "h_mean", "dh_fit_dx", "cycle", "segment_id", "gt", "distance", "w_surface_window_final", "spot", "n_fit_photons"];\n", | |
" let aliases = ["index", "time", "rgt", "rms_misfit", "dh_fit_dy", "pflags", "h_sigma", "h_mean", "dh_fit_dx", "cycle", "segment_id", "gt", "distance", "w_surface_window_final", "spot", "n_fit_photons"];\n", | |
" let table = '<table>' +\n", | |
" String(\n", | |
" fields.map(\n", | |
" (v,i)=>\n", | |
" `<tr>\n", | |
" <th>${aliases[i]}</th>\n", | |
" \n", | |
" <td>${handleObject(layer.feature.properties[v])}</td>\n", | |
" </tr>`).join(''))\n", | |
" +'</table>';\n", | |
" div.innerHTML=table;\n", | |
" \n", | |
" return div\n", | |
" }\n", | |
" ,{"className": "foliumtooltip", "sticky": true});\n", | |
" \n", | |
"</script>\n", | |
"</html>\" style=\"position:absolute;width:100%;height:100%;left:0;top:0;border:none !important;\" allowfullscreen webkitallowfullscreen mozallowfullscreen></iframe></div></div>" | |
], | |
"text/plain": [ | |
"<folium.folium.Map at 0x7f535b702b90>" | |
] | |
}, | |
"execution_count": 6, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"# Always check wehere you data is!\n", | |
"points = gf.reset_index()\n", | |
"points.loc[:, 'time'] = points.time.dt.strftime('%Y-%m-%d')\n", | |
"points.explore()" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 7, | |
"id": "46d8f194-ae4e-44d8-b0c7-6f820ed6d5d3", | |
"metadata": { | |
"tags": [] | |
}, | |
"outputs": [ | |
{ | |
"name": "stdout", | |
"output_type": "stream", | |
"text": [ | |
"-100.94678306984892 -79.00104942913823 -100.63246339163007 -78.98240211911293\n" | |
] | |
} | |
], | |
"source": [ | |
"# Get bounding box of all of our points\n", | |
"w,s,e,n = gf.unary_union.bounds #W, S, E, N\n", | |
"print(w,s,e,n)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 8, | |
"id": "e02040c2-aa7d-4962-affc-fc6430e76aae", | |
"metadata": { | |
"tags": [] | |
}, | |
"outputs": [ | |
{ | |
"name": "stdout", | |
"output_type": "stream", | |
"text": [ | |
"pj_open_lib(proj.ini): call fopen(/srv/conda/envs/notebook/share/proj/proj.ini) - succeeded\n", | |
"pj_open_lib(proj.db): call fopen(/srv/conda/envs/notebook/share/proj/proj.db) - succeeded\n", | |
"pj_open_lib(us_noaa_nadcon5_nad27_nad83_1986_conus.tif): call fopen(/srv/conda/envs/notebook/share/proj/us_noaa_nadcon5_nad27_nad83_1986_conus.tif) - failed\n", | |
"pj_open_lib(proj.db): call fopen(/srv/conda/envs/notebook/share/proj/proj.db) - succeeded\n", | |
"pj_open_lib(us_noaa_nadcon5_nad83_1986_nad83_harn_conus.tif): call fopen(/srv/conda/envs/notebook/share/proj/us_noaa_nadcon5_nad83_1986_nad83_harn_conus.tif) - failed\n", | |
"pj_open_lib(us_noaa_nadcon5_nad83_harn_nad83_fbn_conus.tif): call fopen(/srv/conda/envs/notebook/share/proj/us_noaa_nadcon5_nad83_harn_nad83_fbn_conus.tif) - failed\n", | |
"pj_open_lib(us_noaa_nadcon5_nad83_fbn_nad83_2007_conus.tif): call fopen(/srv/conda/envs/notebook/share/proj/us_noaa_nadcon5_nad83_fbn_nad83_2007_conus.tif) - failed\n", | |
"pj_open_lib(us_noaa_nadcon5_nad83_2007_nad83_2011_conus.tif): call fopen(/srv/conda/envs/notebook/share/proj/us_noaa_nadcon5_nad83_2007_nad83_2011_conus.tif) - failed\n", | |
"pj_open_lib(us_noaa_alhpgn.tif): call fopen(/srv/conda/envs/notebook/share/proj/us_noaa_alhpgn.tif) - failed\n", | |
"pj_open_lib(alhpgn.gsb): call fopen(/srv/conda/envs/notebook/share/proj/alhpgn.gsb) - failed\n", | |
"pj_open_lib(us_noaa_conus.tif): call fopen(/srv/conda/envs/notebook/share/proj/us_noaa_conus.tif) - failed\n", | |
"pj_open_lib(conus): call fopen(/srv/conda/envs/notebook/share/proj/conus) - failed\n", | |
"pj_open_lib(us_noaa_FL.tif): call fopen(/srv/conda/envs/notebook/share/proj/us_noaa_FL.tif) - failed\n", | |
"pj_open_lib(FL): call fopen(/srv/conda/envs/notebook/share/proj/FL) - failed\n", | |
"pj_open_lib(us_noaa_lahpgn.tif): call fopen(/srv/conda/envs/notebook/share/proj/us_noaa_lahpgn.tif) - failed\n", | |
"pj_open_lib(lahpgn.gsb): call fopen(/srv/conda/envs/notebook/share/proj/lahpgn.gsb) - failed\n", | |
"pj_open_lib(us_noaa_mshpgn.tif): call fopen(/srv/conda/envs/notebook/share/proj/us_noaa_mshpgn.tif) - failed\n", | |
"pj_open_lib(mshpgn.gsb): call fopen(/srv/conda/envs/notebook/share/proj/mshpgn.gsb) - failed\n", | |
"pj_open_lib(us_noaa_ethpgn.tif): call fopen(/srv/conda/envs/notebook/share/proj/us_noaa_ethpgn.tif) - failed\n", | |
"pj_open_lib(ethpgn.gsb): call fopen(/srv/conda/envs/notebook/share/proj/ethpgn.gsb) - failed\n", | |
"pj_open_lib(au_icsm_GDA94_GDA2020_conformal_and_distortion.tif): call fopen(/srv/conda/envs/notebook/share/proj/au_icsm_GDA94_GDA2020_conformal_and_distortion.tif) - failed\n", | |
"pj_open_lib(GDA94_GDA2020_conformal_and_distortion.gsb): call fopen(/srv/conda/envs/notebook/share/proj/GDA94_GDA2020_conformal_and_distortion.gsb) - failed\n", | |
"pj_open_lib(eur_nkg_nkgrf17vel.tif): call fopen(/srv/conda/envs/notebook/share/proj/eur_nkg_nkgrf17vel.tif) - failed\n", | |
"pj_open_lib(no_kv_NKGETRF14_EPSG7922_2000.tif): call fopen(/srv/conda/envs/notebook/share/proj/no_kv_NKGETRF14_EPSG7922_2000.tif) - failed\n", | |
"pj_open_lib(Und_min1x1_egm2008_isw=82_WGS84_TideFree.gz): call fopen(/srv/conda/envs/notebook/share/proj/Und_min1x1_egm2008_isw=82_WGS84_TideFree.gz) - failed\n", | |
"pj_open_lib(us_nga_egm08_25.tif): call fopen(/srv/conda/envs/notebook/share/proj/us_nga_egm08_25.tif) - failed\n", | |
"pj_open_lib(egm08_25.gtx): call fopen(/srv/conda/envs/notebook/share/proj/egm08_25.gtx) - failed\n", | |
"pj_open_lib(us_nga_egm08_25.tif): call fopen(/srv/conda/envs/notebook/share/proj/us_nga_egm08_25.tif) - failed\n", | |
"pj_open_lib(egm08_25.gtx): call fopen(/srv/conda/envs/notebook/share/proj/egm08_25.gtx) - failed\n", | |
"Using https://cdn.proj.org/us_nga_egm08_25.tif\n", | |
"Candidate operations found: 1\n", | |
"-------------------------------------\n", | |
"Operation No. 1:\n", | |
"\n", | |
"pj_open_lib(us_nga_egm08_25.tif): call fopen(/srv/conda/envs/notebook/share/proj/us_nga_egm08_25.tif) - failed\n", | |
"pj_open_lib(egm08_25.gtx): call fopen(/srv/conda/envs/notebook/share/proj/egm08_25.gtx) - failed\n", | |
"Using https://cdn.proj.org/us_nga_egm08_25.tif\n", | |
"unknown id, Inverse of Conversion from ITRF2014 (geocentric) to ITRF2014 (geog3D) + Inverse of WGS 84 (G2139) to ITRF2014 (1) + Inverse of Conversion from WGS 84 (G2139) (geog2D) to WGS 84 (G2139) (geocentric) + Inverse of WGS 84 to WGS 84 (G2139) + Inverse of Null geographic offset from WGS 84 (geog3D) to WGS 84 (geog2D) + WGS 84 to EGM2008 height (1), 3.01 m, World\n", | |
"\n", | |
"PROJ string:\n", | |
"+proj=pipeline\n", | |
" +step +proj=axisswap +order=2,1\n", | |
" +step +proj=unitconvert +xy_in=deg +xy_out=rad\n", | |
" +step +inv +proj=vgridshift +grids=us_nga_egm08_25.tif +multiplier=1\n", | |
" +step +proj=unitconvert +xy_in=rad +xy_out=deg\n", | |
" +step +proj=axisswap +order=2,1\n" | |
] | |
} | |
], | |
"source": [ | |
"#!projinfo -s EPSG:7912 -t EPSG:9518 -o PROJ --hide-ballpark --spatial-test intersects\n", | |
"!projinfo -s EPSG:7912 -t EPSG:9518 -o PROJ --bbox {w},{s},{e},{n} --hide-ballpark --spatial-test intersects" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 9, | |
"id": "4006bcc0-14dd-4bbd-aa8f-9e37df3ff1fc", | |
"metadata": { | |
"tags": [] | |
}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/html": [ | |
"<div>\n", | |
"<style scoped>\n", | |
" .dataframe tbody tr th:only-of-type {\n", | |
" vertical-align: middle;\n", | |
" }\n", | |
"\n", | |
" .dataframe tbody tr th {\n", | |
" vertical-align: top;\n", | |
" }\n", | |
"\n", | |
" .dataframe thead th {\n", | |
" text-align: right;\n", | |
" }\n", | |
"</style>\n", | |
"<table border=\"1\" class=\"dataframe\">\n", | |
" <thead>\n", | |
" <tr style=\"text-align: right;\">\n", | |
" <th></th>\n", | |
" <th>geometry</th>\n", | |
" </tr>\n", | |
" </thead>\n", | |
" <tbody>\n", | |
" <tr>\n", | |
" <th>0</th>\n", | |
" <td>POINT Z (-100.63246 -78.98922 1875.21728)</td>\n", | |
" </tr>\n", | |
" <tr>\n", | |
" <th>1</th>\n", | |
" <td>POINT Z (-100.63265 -78.98904 1875.36690)</td>\n", | |
" </tr>\n", | |
" </tbody>\n", | |
"</table>\n", | |
"</div>" | |
], | |
"text/plain": [ | |
" geometry\n", | |
"0 POINT Z (-100.63246 -78.98922 1875.21728)\n", | |
"1 POINT Z (-100.63265 -78.98904 1875.36690)" | |
] | |
}, | |
"execution_count": 9, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"# NOTE: only vshift so only expect vertical coordinate change\n", | |
"points3D = gpd.points_from_xy(gf.geometry.x, gf.geometry.y, gf.h_mean)\n", | |
"gf3D = gpd.GeoDataFrame(geometry=points3D, crs='EPSG:7912')\n", | |
"gf3D.head(2)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 10, | |
"id": "38a740de-3c30-4434-a865-22374c7a7f19", | |
"metadata": { | |
"tags": [] | |
}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"<Geographic 3D CRS: EPSG:7912>\n", | |
"Name: ITRF2014\n", | |
"Axis Info [ellipsoidal]:\n", | |
"- Lat[north]: Geodetic latitude (degree)\n", | |
"- Lon[east]: Geodetic longitude (degree)\n", | |
"- h[up]: Ellipsoidal height (metre)\n", | |
"Area of Use:\n", | |
"- name: World.\n", | |
"- bounds: (-180.0, -90.0, 180.0, 90.0)\n", | |
"Datum: International Terrestrial Reference Frame 2014\n", | |
"- Ellipsoid: GRS 1980\n", | |
"- Prime Meridian: Greenwich" | |
] | |
}, | |
"execution_count": 10, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"gf3D.crs" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 11, | |
"id": "af6a3674-4552-4242-8a62-5753a0ed1be4", | |
"metadata": { | |
"tags": [] | |
}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"array([[-100.63246339, -78.98921541, 1875.21727596]])" | |
] | |
}, | |
"execution_count": 11, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"shapely.get_coordinates(gf3D.iloc[0], include_z=True)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 12, | |
"id": "84ea693d-64b1-463c-889d-2a79295dbaf5", | |
"metadata": { | |
"tags": [] | |
}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"array([[-100.63246339, -78.98921541, 1904.67380708]])" | |
] | |
}, | |
"execution_count": 12, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"gfGeoid = gf3D.to_crs(epsg=9518)\n", | |
"shapely.get_coordinates(gfGeoid.iloc[0], include_z=True)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 13, | |
"id": "612f8f95-b2ad-4ef1-8b08-df0b0f37eb6c", | |
"metadata": { | |
"tags": [] | |
}, | |
"outputs": [], | |
"source": [ | |
"#gfGeoid.head()" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 14, | |
"id": "30c8b1dc-da22-4556-8f7b-fec65ff10d3c", | |
"metadata": { | |
"tags": [] | |
}, | |
"outputs": [ | |
{ | |
"name": "stderr", | |
"output_type": "stream", | |
"text": [ | |
"DEBUG:fiona._env:GDAL_DATA found in environment.\n", | |
"DEBUG:fiona._env:PROJ_DATA found in environment.\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(us_noaa_nadcon5_nad83_harn_nad83_fbn_conus.tif): call fopen(/home/jovyan/.local/share/proj/us_noaa_nadcon5_nad83_harn_nad83_fbn_conus.tif) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(us_noaa_nadcon5_nad83_fbn_nad83_2007_conus.tif): call fopen(/home/jovyan/.local/share/proj/us_noaa_nadcon5_nad83_fbn_nad83_2007_conus.tif) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(us_noaa_nadcon5_nad83_2007_nad83_2011_conus.tif): call fopen(/home/jovyan/.local/share/proj/us_noaa_nadcon5_nad83_2007_nad83_2011_conus.tif) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(us_noaa_g2018u0.tif): call fopen(/home/jovyan/.local/share/proj/us_noaa_g2018u0.tif) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(g2018u0.gtx): call fopen(/home/jovyan/.local/share/proj/g2018u0.gtx) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(us_noaa_geoid03_conus.tif): call fopen(/home/jovyan/.local/share/proj/us_noaa_geoid03_conus.tif) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(geoid03_conus.gtx): call fopen(/home/jovyan/.local/share/proj/geoid03_conus.gtx) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(us_noaa_geoid09_conus.tif): call fopen(/home/jovyan/.local/share/proj/us_noaa_geoid09_conus.tif) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(geoid09_conus.gtx): call fopen(/home/jovyan/.local/share/proj/geoid09_conus.gtx) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(us_noaa_vertconw.tif): call fopen(/home/jovyan/.local/share/proj/us_noaa_vertconw.tif) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(vertconw.gtx): call fopen(/home/jovyan/.local/share/proj/vertconw.gtx) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(us_noaa_g1999u01.tif): call fopen(/home/jovyan/.local/share/proj/us_noaa_g1999u01.tif) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(g1999u01.gtx): call fopen(/home/jovyan/.local/share/proj/g1999u01.gtx) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(us_noaa_g2018u0.tif): call fopen(/home/jovyan/.local/share/proj/us_noaa_g2018u0.tif) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(g2018u0.gtx): call fopen(/home/jovyan/.local/share/proj/g2018u0.gtx) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: Using https://cdn.proj.org/us_noaa_g2018u0.tif\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(us_noaa_nadcon5_nad83_1993_nad83_1997_prvi.tif): call fopen(/home/jovyan/.local/share/proj/us_noaa_nadcon5_nad83_1993_nad83_1997_prvi.tif) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(us_noaa_nadcon5_nad83_1997_nad83_2002_prvi.tif): call fopen(/home/jovyan/.local/share/proj/us_noaa_nadcon5_nad83_1997_nad83_2002_prvi.tif) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(us_noaa_nadcon5_nad83_2007_nad83_2011_conus.tif): call fopen(/home/jovyan/.local/share/proj/us_noaa_nadcon5_nad83_2007_nad83_2011_conus.tif) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: Using https://cdn.proj.org/us_noaa_nadcon5_nad83_2007_nad83_2011_conus.tif\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(us_noaa_nadcon5_nad83_fbn_nad83_2007_conus.tif): call fopen(/home/jovyan/.local/share/proj/us_noaa_nadcon5_nad83_fbn_nad83_2007_conus.tif) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: Using https://cdn.proj.org/us_noaa_nadcon5_nad83_fbn_nad83_2007_conus.tif\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(us_noaa_nadcon5_nad83_harn_nad83_fbn_conus.tif): call fopen(/home/jovyan/.local/share/proj/us_noaa_nadcon5_nad83_harn_nad83_fbn_conus.tif) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: Using https://cdn.proj.org/us_noaa_nadcon5_nad83_harn_nad83_fbn_conus.tif\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(us_noaa_geoid03_conus.tif): call fopen(/home/jovyan/.local/share/proj/us_noaa_geoid03_conus.tif) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(geoid03_conus.gtx): call fopen(/home/jovyan/.local/share/proj/geoid03_conus.gtx) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: Using https://cdn.proj.org/us_noaa_geoid03_conus.tif\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(us_noaa_geoid09_conus.tif): call fopen(/home/jovyan/.local/share/proj/us_noaa_geoid09_conus.tif) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(geoid09_conus.gtx): call fopen(/home/jovyan/.local/share/proj/geoid09_conus.gtx) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: Using https://cdn.proj.org/us_noaa_geoid09_conus.tif\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(us_noaa_g1999u01.tif): call fopen(/home/jovyan/.local/share/proj/us_noaa_g1999u01.tif) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(g1999u01.gtx): call fopen(/home/jovyan/.local/share/proj/g1999u01.gtx) - failed\n", | |
"DEBUG:pyproj:PROJ_DEBUG: Using https://cdn.proj.org/us_noaa_g1999u01.tif\n", | |
"DEBUG:pyproj:PROJ_DEBUG: Using coordinate operation Inverse of Transformation from NAVD88 height to ITRF2014 (ballpark vertical transformation, without ellipsoid height to vertical height correction) + Inverse of Ballpark geographic offset from NAD83(HARN) to ITRF2014 + SPCS83 Washington South zone (US Survey feet)\n" | |
] | |
} | |
], | |
"source": [ | |
"# SHow debug log\n", | |
"\n", | |
"# CAREFUL, geopandas will happily do bogus transforms for you\n", | |
"\n", | |
"import logging\n", | |
"logging.basicConfig(level=logging.DEBUG)\n", | |
"\n", | |
"with fiona.Env(CPL_DEBUG=True):\n", | |
" test = gf3D.to_crs(epsg=\"2927+5703\")" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 15, | |
"id": "daf65bfc-2dff-4a88-a0dd-d2ef04c1fa76", | |
"metadata": { | |
"tags": [] | |
}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/html": [ | |
"<div>\n", | |
"<style scoped>\n", | |
" .dataframe tbody tr th:only-of-type {\n", | |
" vertical-align: middle;\n", | |
" }\n", | |
"\n", | |
" .dataframe tbody tr th {\n", | |
" vertical-align: top;\n", | |
" }\n", | |
"\n", | |
" .dataframe thead th {\n", | |
" text-align: right;\n", | |
" }\n", | |
"</style>\n", | |
"<table border=\"1\" class=\"dataframe\">\n", | |
" <thead>\n", | |
" <tr style=\"text-align: right;\">\n", | |
" <th></th>\n", | |
" <th>geometry</th>\n", | |
" </tr>\n", | |
" </thead>\n", | |
" <tbody>\n", | |
" <tr>\n", | |
" <th>0</th>\n", | |
" <td>POINT Z (53991548.388 -183138438.398 1875.217)</td>\n", | |
" </tr>\n", | |
" <tr>\n", | |
" <th>1</th>\n", | |
" <td>POINT Z (53990465.763 -183136201.989 1875.367)</td>\n", | |
" </tr>\n", | |
" <tr>\n", | |
" <th>2</th>\n", | |
" <td>POINT Z (53989380.806 -183133966.815 1875.252)</td>\n", | |
" </tr>\n", | |
" <tr>\n", | |
" <th>3</th>\n", | |
" <td>POINT Z (53988291.113 -183131734.019 1875.504)</td>\n", | |
" </tr>\n", | |
" <tr>\n", | |
" <th>4</th>\n", | |
" <td>POINT Z (53987200.796 -183129501.561 1875.526)</td>\n", | |
" </tr>\n", | |
" </tbody>\n", | |
"</table>\n", | |
"</div>" | |
], | |
"text/plain": [ | |
" geometry\n", | |
"0 POINT Z (53991548.388 -183138438.398 1875.217)\n", | |
"1 POINT Z (53990465.763 -183136201.989 1875.367)\n", | |
"2 POINT Z (53989380.806 -183133966.815 1875.252)\n", | |
"3 POINT Z (53988291.113 -183131734.019 1875.504)\n", | |
"4 POINT Z (53987200.796 -183129501.561 1875.526)" | |
] | |
}, | |
"execution_count": 15, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"test.head()" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": null, | |
"id": "2308788e-9479-450d-b678-dfe000a9a70c", | |
"metadata": {}, | |
"outputs": [], | |
"source": [] | |
} | |
], | |
"metadata": { | |
"kernelspec": { | |
"display_name": "Python 3 (ipykernel)", | |
"language": "python", | |
"name": "python3" | |
}, | |
"language_info": { | |
"codemirror_mode": { | |
"name": "ipython", | |
"version": 3 | |
}, | |
"file_extension": ".py", | |
"mimetype": "text/x-python", | |
"name": "python", | |
"nbconvert_exporter": "python", | |
"pygments_lexer": "ipython3", | |
"version": "3.10.12" | |
} | |
}, | |
"nbformat": 4, | |
"nbformat_minor": 5 | |
} |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment