Spilhaus in the browser

%23%20%2F%2F%2F%20script%0A%23%20requires-python%20%3D%20%22%3E%3D3.11%22%0A%23%20dependencies%20%3D%20%5B%0A%23%20%20%20%20%20%22marimo%3E%3D0.24%22%2C%0A%23%20%20%20%20%20%22numpy%22%2C%0A%23%20%20%20%20%20%22pyproj%22%2C%0A%23%20%20%20%20%20%22shapely%3E%3D2.0%22%2C%0A%23%20%20%20%20%20%22geopandas%3E%3D1.0%22%2C%0A%23%20%20%20%20%20%22matplotlib%22%2C%0A%23%20%5D%0A%23%20%2F%2F%2F%0A%22%22%22%0ASpilhaus%20in%20the%20browser%3A%20the%20vector-only%20half%20of%20spilhaus.py%2C%20exported%20with%0A%60marimo%20export%20html-wasm%60%20and%20run%20by%20Pyodide.%0A%0APyodide%20ships%20its%20own%20pyproj%20build%2C%20with%20its%20own%20PROJ.%20There%20is%20no%20GDAL%20in%0Athe%20browser%2C%20so%20no%20rasterio%2C%20pyogrio%20or%20%2Fvsicurl%20here%3A%20this%20page%20is%20exactly%0Athe%20%22which%20PROJ%20did%20I%20actually%20get%22%20question%2C%20asked%20of%20one%20more%20runtime.%0A%22%22%22%0A%0Aimport%20marimo%0A%0A__generated_with%20%3D%20%220.24.2%22%0Aapp%20%3D%20marimo.App(width%3D%22medium%22%2C%20app_title%3D%22Spilhaus%20in%20the%20browser%22)%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20import%20marimo%20as%20mo%0A%0A%20%20%20%20return%20(mo%2C)%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%20Spilhaus%20in%20the%20browser%0A%0A%20%20%20%20This%20page%20is%20Python%20running%20in%20your%20browser%20via%20**Pyodide**.%20It%20carries%20its%0A%20%20%20%20own%20%60pyproj%60%2C%20built%20against%20its%20own%20PROJ.%20Whether%20%60%2Bproj%3Dspilhaus%60%20works%20here%0A%20%20%20%20is%20decided%20by%20whichever%20PROJ%20the%20Pyodide%20maintainers%20built%20against%20-%20a%20fourth%0A%20%20%20%20or%20fifth%20answer%20to%20%22what%20PROJ%20do%20I%20have%22%2C%20alongside%20the%20ones%20in%20the%0A%20%20%20%20%5Bfull%20notebook%5D(..%2Fpython%2F)%20and%20%5Bthe%20R%20version%5D(..%2Fr%2F)%3B%20%5Bindex%5D(..%2F).%0A%0A%20%20%20%20There%20is%20no%20GDAL%20in%20the%20browser%2C%20so%20this%20page%20does%20the%20vector%20and%20point%20half%0A%20%20%20%20only%3A%20a%20hand-typed%20set%20of%20landmarks%2C%20and%20Natural%20Earth%20coastlines%20fetched%20as%0A%20%20%20%20GeoJSON%20with%20the%20browser's%20own%20%60fetch%60.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20import%20sys%0A%0A%20%20%20%20import%20numpy%20as%20np%0A%0A%20%20%20%20rows%20%3D%20%5B%5D%0A%20%20%20%20engine%20%3D%20None%0A%20%20%20%20try%3A%0A%20%20%20%20%20%20%20%20import%20pyproj%0A%0A%20%20%20%20%20%20%20%20row%20%3D%20dict(package%3D%22pyproj%22%2C%20version%3Dpyproj.__version__%2C%20proj%3Dpyproj.proj_version_str)%0A%20%20%20%20%20%20%20%20try%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20crs%20%3D%20pyproj.CRS.from_proj4(%22%2Bproj%3Dspilhaus%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20fwd%20%3D%20pyproj.Transformer.from_crs(4326%2C%20crs%2C%20always_xy%3DTrue)%0A%20%20%20%20%20%20%20%20%20%20%20%20engine%20%3D%20fwd%0A%20%20%20%20%20%20%20%20%20%20%20%20row%5B%22spilhaus%22%5D%20%3D%20%22ok%22%0A%20%20%20%20%20%20%20%20except%20Exception%20as%20e%3A%20%20%23%20noqa%3A%20BLE001%0A%20%20%20%20%20%20%20%20%20%20%20%20row%5B%22spilhaus%22%5D%20%3D%20f%22FAIL%3A%20%7Btype(e).__name__%7D%22%0A%20%20%20%20%20%20%20%20rows.append(row)%0A%20%20%20%20except%20ImportError%3A%0A%20%20%20%20%20%20%20%20rows.append(dict(package%3D%22pyproj%22%2C%20version%3DNone%2C%20proj%3DNone%2C%20spilhaus%3D%22not%20installed%22))%0A%0A%20%20%20%20try%3A%0A%20%20%20%20%20%20%20%20import%20shapely%20as%20_shp%0A%0A%20%20%20%20%20%20%20%20rows.append(dict(package%3D%22shapely%22%2C%20version%3D_shp.__version__%2C%20proj%3Df%22GEOS%20%7B_shp.geos_version_string%7D%22%2C%20spilhaus%3D%22n%2Fa%22))%0A%20%20%20%20except%20ImportError%3A%0A%20%20%20%20%20%20%20%20pass%0A%0A%20%20%20%20runtime%20%3D%20%22Pyodide%20(browser)%22%20if%20sys.platform%20%3D%3D%20%22emscripten%22%20else%20f%22CPython%20%7Bsys.version.split()%5B0%5D%7D%20(%7Bsys.platform%7D)%22%0A%0A%20%20%20%20mo.vstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(f%22Runtime%3A%20**%7Bruntime%7D**%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.ui.table(rows%2C%20selection%3DNone%2C%20label%3D%22Native%20library%20provenance%20in%20this%20runtime%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.callout(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22%60%2Bproj%3Dspilhaus%60%20is%20available%20here%3B%20the%20map%20below%20is%20drawn%20with%20it.%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20engine%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20else%20%22%60%2Bproj%3Dspilhaus%60%20is%20**not**%20available%20in%20this%20runtime's%20PROJ.%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22The%20table%20above%20is%20the%20whole%20result%3A%20nothing%20below%20can%20be%20drawn%20until%20Pyodide's%20pyproj%20is%20rebuilt%20against%20PROJ%20%3E%3D%209.6.%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20kind%3D%22success%22%20if%20engine%20else%20%22warn%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%20%20%20%20return%20engine%2C%20np%0A%0A%0A%40app.cell%0Adef%20_(engine%2C%20mo%2C%20np)%3A%0A%20%20%20%20mo.stop(engine%20is%20None%2C%20mo.md(%22_Skipping%20the%20map%3A%20no%20Spilhaus-capable%20PROJ%20in%20this%20runtime._%22))%0A%0A%20%20%20%20def%20forward(lon%2C%20lat)%3A%0A%20%20%20%20%20%20%20%20x%2C%20y%20%3D%20engine.transform(np.asarray(lon%2C%20float)%2C%20np.asarray(lat%2C%20float))%0A%20%20%20%20%20%20%20%20x%20%3D%20np.asarray(x%2C%20float)%0A%20%20%20%20%20%20%20%20y%20%3D%20np.asarray(y%2C%20float)%0A%20%20%20%20%20%20%20%20bad%20%3D%20~np.isfinite(x)%20%7C%20~np.isfinite(y)%20%7C%20(np.abs(x)%20%3E%201e15)%0A%20%20%20%20%20%20%20%20x%5Bbad%5D%20%3D%20np.nan%0A%20%20%20%20%20%20%20%20y%5Bbad%5D%20%3D%20np.nan%0A%20%20%20%20%20%20%20%20return%20x%2C%20y%0A%0A%20%20%20%20_lon%2C%20_lat%20%3D%20np.meshgrid(np.linspace(-180%2C%20180%2C%20721)%2C%20np.linspace(-90%2C%2090%2C%20361))%0A%20%20%20%20_x%2C%20_y%20%3D%20forward(_lon%2C%20_lat)%0A%20%20%20%20HALF%20%3D%20float(np.ceil(np.nanmax(np.abs(np.concatenate(%5B_x%2C%20_y%5D)))%20%2F%201000.0)%20*%201000.0)%0A%20%20%20%20return%20HALF%2C%20forward%0A%0A%0A%40app.cell%0Adef%20_(HALF%2C%20forward%2C%20mo%2C%20np)%3A%0A%20%20%20%20import%20json%0A%20%20%20%20import%20urllib.request%0A%0A%20%20%20%20NE_GEOJSON%20%3D%20(%0A%20%20%20%20%20%20%20%20%22https%3A%2F%2Fraw.githubusercontent.com%2Fnvkelso%2Fnatural-earth-vector%2Fmaster%2Fgeojson%2Fne_110m_coastline.geojson%22%0A%20%20%20%20)%0A%0A%20%20%20%20def%20_fetch_json(url)%3A%0A%20%20%20%20%20%20%20%20%23%20In%20Pyodide%20urllib%20is%20patched%20to%20use%20the%20browser's%20fetch%3B%20on%20CPython%20it%20is%20plain%20HTTP.%0A%20%20%20%20%20%20%20%20with%20urllib.request.urlopen(url%2C%20timeout%3D60)%20as%20r%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20json.loads(r.read().decode(%22utf-8%22))%0A%0A%20%20%20%20try%3A%0A%20%20%20%20%20%20%20%20_gj%20%3D%20_fetch_json(NE_GEOJSON)%0A%20%20%20%20%20%20%20%20_status%20%3D%20f%22fetched%20%7Blen(_gj%5B'features'%5D)%7D%20coastline%20features%22%0A%20%20%20%20except%20Exception%20as%20e%3A%20%20%23%20noqa%3A%20BLE001%0A%20%20%20%20%20%20%20%20_gj%20%3D%20%7B%22features%22%3A%20%5B%5D%7D%0A%20%20%20%20%20%20%20%20_status%20%3D%20f%22coastline%20fetch%20failed%20(%7Btype(e).__name__%7D)%3B%20drawing%20landmarks%20only%22%0A%0A%20%20%20%20import%20shapely%0A%20%20%20%20from%20shapely.geometry%20import%20shape%0A%0A%20%20%20%20coast_segments%20%3D%20%5B%5D%0A%20%20%20%20for%20_f%20in%20_gj%5B%22features%22%5D%3A%0A%20%20%20%20%20%20%20%20_g%20%3D%20shapely.segmentize(shape(_f%5B%22geometry%22%5D)%2C%200.5)%0A%20%20%20%20%20%20%20%20for%20_line%20in%20getattr(_g%2C%20%22geoms%22%2C%20%5B_g%5D)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_c%20%3D%20np.asarray(_line.coords)%0A%20%20%20%20%20%20%20%20%20%20%20%20_x%2C%20_y%20%3D%20forward(_c%5B%3A%2C%200%5D%2C%20_c%5B%3A%2C%201%5D)%0A%20%20%20%20%20%20%20%20%20%20%20%20_step%20%3D%20np.hypot(np.diff(_x)%2C%20np.diff(_y))%0A%20%20%20%20%20%20%20%20%20%20%20%20_ok%20%3D%20np.isfinite(_step)%20%26%20(_step%20%3C%20HALF%20%2F%208)%0A%20%20%20%20%20%20%20%20%20%20%20%20_start%20%3D%200%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20_i%20in%20list(np.flatnonzero(~_ok))%20%2B%20%5Blen(_step)%5D%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20_i%20-%20_start%20%3E%3D%201%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20coast_segments.append(np.c_%5B_x%2C%20_y%5D%5B_start%20%3A%20_i%20%2B%201%5D)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_start%20%3D%20_i%20%2B%201%0A%0A%20%20%20%20mo.md(f%22Coastlines%3A%20%7B_status%7D%3B%20%7Blen(coast_segments)%7D%20drawable%20runs.%22)%0A%20%20%20%20return%20(coast_segments%2C)%0A%0A%0A%40app.cell%0Adef%20_(HALF%2C%20coast_segments%2C%20forward%2C%20mo)%3A%0A%20%20%20%20import%20matplotlib%0A%0A%20%20%20%20matplotlib.use(%22Agg%22)%0A%20%20%20%20import%20matplotlib.pyplot%20as%20plt%0A%20%20%20%20from%20matplotlib%20import%20patheffects%0A%20%20%20%20from%20matplotlib.collections%20import%20LineCollection%0A%0A%20%20%20%20landmarks%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22Hobart%22%3A%20(147.33%2C%20-42.88)%2C%0A%20%20%20%20%20%20%20%20%22Tokyo%22%3A%20(139.69%2C%2035.69)%2C%0A%20%20%20%20%20%20%20%20%22Reykjavik%22%3A%20(-21.94%2C%2064.15)%2C%0A%20%20%20%20%20%20%20%20%22Cape%20of%20Good%20Hope%22%3A%20(18.47%2C%20-34.36)%2C%0A%20%20%20%20%20%20%20%20%22Bering%20Strait%22%3A%20(-169.0%2C%2065.8)%2C%0A%20%20%20%20%20%20%20%20%22Drake%20Passage%22%3A%20(-63.0%2C%20-59.0)%2C%0A%20%20%20%20%20%20%20%20%22Point%20Nemo%22%3A%20(-123.39%2C%20-48.88)%2C%0A%20%20%20%20%20%20%20%20%22Challenger%20Deep%22%3A%20(142.20%2C%2011.37)%2C%0A%20%20%20%20%20%20%20%20%22Puerto%20Rico%20Trench%22%3A%20(-66.5%2C%2019.7)%2C%0A%20%20%20%20%20%20%20%20%22Mid-Atlantic%20Ridge%20(Azores)%22%3A%20(-27.0%2C%2038.0)%2C%0A%20%20%20%20%20%20%20%20%22Macquarie%20Island%22%3A%20(158.94%2C%20-54.62)%2C%0A%20%20%20%20%20%20%20%20%22Mawson%20Station%22%3A%20(62.87%2C%20-67.60)%2C%0A%20%20%20%20%20%20%20%20%22Casey%20Station%22%3A%20(110.53%2C%20-66.28)%2C%0A%20%20%20%20%7D%0A%20%20%20%20_lx%2C%20_ly%20%3D%20forward(%5Bv%5B0%5D%20for%20v%20in%20landmarks.values()%5D%2C%20%5Bv%5B1%5D%20for%20v%20in%20landmarks.values()%5D)%0A%0A%20%20%20%20fig%2C%20ax%20%3D%20plt.subplots(figsize%3D(9%2C%209))%0A%20%20%20%20ax.set_facecolor(%22%23e9eef4%22)%0A%20%20%20%20ax.add_patch(plt.Rectangle((-HALF%2C%20-HALF)%2C%202%20*%20HALF%2C%202%20*%20HALF%2C%20fc%3D%22%23e9eef4%22%2C%20ec%3D%22%23333333%22%2C%20lw%3D0.8))%0A%20%20%20%20ax.add_collection(LineCollection(coast_segments%2C%20colors%3D%22%232d4a6d%22%2C%20linewidths%3D0.7))%0A%20%20%20%20ax.plot(_lx%2C%20_ly%2C%20%22o%22%2C%20ms%3D5%2C%20color%3D%22%23f2b134%22%2C%20mec%3D%22black%22%2C%20mew%3D0.6%2C%20zorder%3D5)%0A%20%20%20%20for%20_name%2C%20_x%2C%20_y%20in%20zip(landmarks%2C%20_lx%2C%20_ly)%3A%0A%20%20%20%20%20%20%20%20ax.annotate(%0A%20%20%20%20%20%20%20%20%20%20%20%20_name%2C%20(_x%2C%20_y)%2C%20xytext%3D(5%2C%204)%2C%20textcoords%3D%22offset%20points%22%2C%20fontsize%3D8.5%2C%20zorder%3D6%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20path_effects%3D%5Bpatheffects.withStroke(linewidth%3D2.5%2C%20foreground%3D%22white%22)%5D%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20ax.set_xlim(-HALF%2C%20HALF)%0A%20%20%20%20ax.set_ylim(-HALF%2C%20HALF)%0A%20%20%20%20ax.set_aspect(%22equal%22)%0A%20%20%20%20ax.set_axis_off()%0A%20%20%20%20ax.set_title(%22Spilhaus%2C%20drawn%20by%20pyproj%20running%20in%20your%20browser%22%2C%20pad%3D12)%0A%20%20%20%20fig.tight_layout()%0A%20%20%20%20mo.vstack(%5Bfig%5D)%0A%20%20%20%20return%0A%0A%0Aif%20__name__%20%3D%3D%20%22__main__%22%3A%0A%20%20%20%20app.run()%0A