Last active
July 27, 2025 20:07
-
-
Save technocrat/96679a50eaceecc6ceda88353383216c to your computer and use it in GitHub Desktop.
Session log of a hello world
This file contains hidden or 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
| julia> import CairoMakie as Mke | |
| julia> using GeoIO | |
| [ Info: Precompiling GeoIO [f5a160d5-e41d-4189-8b61-d57781c419e3] (cache misses: wrong dep version loaded (4)) | |
| WARNING: Method definition apply_recipe(Base.AbstractDict{Symbol, Any}, Union{GeoInterface.MultiPointTrait, GeoInterface.PointTrait}, Any) in module GeoInterfaceRecipesBaseExt at /Users/technocrat/.julia/packages/RecipesBase/BRe07/src/RecipesBase.jl:296 overwritten in module GeoInterfaceRecipes on the same line (check for duplicate calls to `include`). | |
| ERROR: Method overwriting is not permitted during Module precompilation. Use `__precompile__(false)` to opt-out of precompilation. | |
| ┌ Info: Skipping precompilation due to precompilable error. Importing GeoIO [f5a160d5-e41d-4189-8b61-d57781c419e3]. | |
| └ exception = Error when precompiling module, potentially caused by a __precompile__(false) declaration in the module. | |
| [ Info: Precompiling Shapefile [8e980c4a-a4fe-5da2-b3a7-4b4b0353a2f4] | |
| WARNING: Method definition apply_recipe(Base.AbstractDict{Symbol, Any}, Union{GeoInterface.MultiPointTrait, GeoInterface.PointTrait}, Any) in module GeoInterfaceRecipesBaseExt at /Users/technocrat/.julia/packages/RecipesBase/BRe07/src/RecipesBase.jl:296 overwritten in module GeoInterfaceRecipes on the same line (check for duplicate calls to `include`). | |
| ERROR: Method overwriting is not permitted during Module precompilation. Use `__precompile__(false)` to opt-out of precompilation. | |
| ┌ Info: Skipping precompilation due to precompilable error. Importing Shapefile [8e980c4a-a4fe-5da2-b3a7-4b4b0353a2f4]. | |
| └ exception = Error when precompiling module, potentially caused by a __precompile__(false) declaration in the module. | |
| [ Info: Precompiling GeoInterfaceRecipes [0329782f-3d07-4b52-b9f6-d3137cf03c7a] (cache misses: wrong dep version loaded (2)) | |
| WARNING: Method definition apply_recipe(Base.AbstractDict{Symbol, Any}, Union{GeoInterface.MultiPointTrait, GeoInterface.PointTrait}, Any) in module GeoInterfaceRecipesBaseExt at /Users/technocrat/.julia/packages/RecipesBase/BRe07/src/RecipesBase.jl:296 overwritten in module GeoInterfaceRecipes on the same line (check for duplicate calls to `include`). | |
| ERROR: Method overwriting is not permitted during Module precompilation. Use `__precompile__(false)` to opt-out of precompilation. | |
| ┌ Info: Skipping precompilation due to precompilable error. Importing GeoInterfaceRecipes [0329782f-3d07-4b52-b9f6-d3137cf03c7a]. | |
| └ exception = Error when precompiling module, potentially caused by a __precompile__(false) declaration in the module. | |
| [ Info: Precompiling ShapefileMakieExt [ce22dd63-8704-5a0b-b9df-feab1218054b] (cache misses: wrong dep version loaded (4)) | |
| ┌ Warning: Module Shapefile with build ID ffffffff-ffff-ffff-0b81-ebd073943a75 is missing from the cache. | |
| │ This may mean Shapefile [8e980c4a-a4fe-5da2-b3a7-4b4b0353a2f4] does not support precompilation but is imported by a module that does. | |
| └ @ Base loading.jl:2541 | |
| ┌ Info: Skipping precompilation due to precompilable error. Importing ShapefileMakieExt [ce22dd63-8704-5a0b-b9df-feab1218054b]. | |
| └ exception = Error when precompiling module, potentially caused by a __precompile__(false) declaration in the module. | |
| [ Info: Precompiling ArchGDAL [c9ce4bd3-c3d5-55b8-8973-c0e20141b8c3] (cache misses: wrong dep version loaded (6)) | |
| ┌ Warning: Module GeoInterfaceRecipes with build ID ffffffff-ffff-ffff-47ee-98538561677b is missing from the cache. | |
| │ This may mean GeoInterfaceRecipes [0329782f-3d07-4b52-b9f6-d3137cf03c7a] does not support precompilation but is imported by a module that does. | |
| └ @ Base loading.jl:2541 | |
| ┌ Info: Skipping precompilation due to precompilable error. Importing ArchGDAL [c9ce4bd3-c3d5-55b8-8973-c0e20141b8c3]. | |
| └ exception = Error when precompiling module, potentially caused by a __precompile__(false) declaration in the module. | |
| [ Info: Precompiling ArchGDALMakieExt [ed506d0f-917e-5ac5-847b-cce2fb23f093] (cache misses: wrong dep version loaded (6)) | |
| ┌ Warning: Module ArchGDAL with build ID ffffffff-ffff-ffff-3467-d2cddbf9d1f5 is missing from the cache. | |
| │ This may mean ArchGDAL [c9ce4bd3-c3d5-55b8-8973-c0e20141b8c3] does not support precompilation but is imported by a module that does. | |
| └ @ Base loading.jl:2541 | |
| ┌ Info: Skipping precompilation due to precompilable error. Importing ArchGDALMakieExt [ed506d0f-917e-5ac5-847b-cce2fb23f093]. | |
| └ exception = Error when precompiling module, potentially caused by a __precompile__(false) declaration in the module. | |
| julia> gt = GeoIO.load("/Users/technocrat/projects/Maps/data/2024_shp/cb_2024_us_state_500k.shp") |> Filter(row -> row.NAME ∉ ["Alaska", "Hawaii", "Guam", "Puerto Rico", "American Samoa", "United States Virgin Islands", "Commonwealth of the Northern Mariana Islands"]) | |
| ERROR: UndefVarError: `Filter` not defined in `Main` | |
| Suggestion: check for spelling errors or missing imports. | |
| Hint: a global variable of this name may be made accessible by importing Transducers in the current active module Main | |
| Suggestion: check for spelling errors or missing imports. | |
| Hint: a global variable of this name may be made accessible by importing Transducers in the current active module Main | |
| Stacktrace: | |
| [1] top-level scope | |
| @ REPL[15]:1 | |
| julia> fig = Mke.Figure() | |
| # retry without the missing Filter, knowing that the bounding box is too large | |
| julia> gt = GeoIO.load("/Users/technocrat/projects/Maps/data/2024_shp/cb_2024_us_state_500k.shp") | |
| 56×10 GeoTable over 56 GeometrySet | |
| ┌─────────────┬─────────────┬─────────────┬─────────────┬─────────────┬──────────────┬─────────────┬──────────────┬─────── | |
| │ STATEFP │ STATENS │ GEOIDFQ │ GEOID │ STUSPS │ NAME │ LSAD │ ALAND │ AW ⋯ | |
| │ Categorical │ Categorical │ Categorical │ Categorical │ Categorical │ Categorical │ Categorical │ Categorical │ Cate ⋯ | |
| │ [NoUnits] │ [NoUnits] │ [NoUnits] │ [NoUnits] │ [NoUnits] │ [NoUnits] │ [NoUnits] │ [NoUnits] │ [No ⋯ | |
| ├─────────────┼─────────────┼─────────────┼─────────────┼─────────────┼──────────────┼─────────────┼──────────────┼─────── | |
| │ 35 │ 00897535 │ 0400000US35 │ 35 │ NM │ New Mexico │ 00 │ 314198519809 │ 726 ⋯ | |
| │ ⋮ │ ⋮ │ ⋮ │ ⋮ │ ⋮ │ ⋮ │ ⋮ │ ⋮ │ ⋱ | |
| └─────────────┴─────────────┴─────────────┴─────────────┴─────────────┴──────────────┴─────────────┴──────────────┴─────── | |
| 2 columns and 55 rows omitted | |
| julia> ax = Mke.Axis(fig[1, 1], title = "US 'Lower 48' States") | |
| Axis with 0 plots: | |
| julia> viz!(ax, gt.geometry, color=1:nrow(gt), segmentcolor="black", showsegments=true, segmentsize=0.3f0) | |
| # bring in viz! | |
| ERROR: UndefVarError: `nrow` not defined in `Main` | |
| Suggestion: check for spelling errors or missing imports. | |
| Hint: a global variable of this name may be made accessible by importing GeoTables in the current active module Main | |
| Hint: a global variable of this name may be made accessible by importing DataFrames in the current active module Main | |
| Suggestion: check for spelling errors or missing imports. | |
| Hint: a global variable of this name may be made accessible by importing GeoTables in the current active module Main | |
| Hint: a global variable of this name may be made accessible by importing DataFrames in the current active module Main | |
| Stacktrace: | |
| [1] top-level scope | |
| @ REPL[19]:1 | |
| julia> Pkg.add("GeoStats") | |
| Resolving package versions... | |
| Compat entries added for | |
| No Changes to `~/projects/testbed/Project.toml` | |
| No Changes to `~/projects/testbed/Manifest.toml` | |
| julia> using GeoStats | |
| julia> @time viz!(ax, gt.geometry, color=1:nrow(gt), segmentcolor="black", showsegments=true, segmentsize=0.3f0) | |
| # need GeoStats | |
| 103.097917 seconds (737.73 M allocations: 67.914 GiB, 7.64% gc time, 7.96% compilation time) | |
| Makie.Plot{viz, Tuple{GeometrySet{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}, MultiPolygon{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}, PolyArea{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}, Ring{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}, CircularArrays.CircularVector{Point{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}}, Vector{Point{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}}}}}, Vector{Ring{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}, CircularArrays.CircularVector{Point{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}}, Vector{Point{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}}}}}}}}}}} | |
| julia> Pkg.add("GeoStats") | |
| Resolving package versions... | |
| Compat entries added for | |
| No Changes to `~/projects/testbed/Project.toml` | |
| No Changes to `~/projects/testbed/Manifest.toml` | |
| julia> using GeoStats | |
| julia> gt = GeoIO.load("/Users/technocrat/projects/Maps/data/2024_shp/cb_2024_us_state_500k.shp") |> Filter(row -> row.NAME ∉ ["Alaska", "Hawaii", "Guam", "Puerto Rico", "American Samoa", "United States Virgin Islands", "Commonwealth of the Northern Mariana Islands"]) | |
| 49×10 GeoTable over 49 view(::GeometrySet, [1, 2, 3, 4, ..., 51, 52, 53, 56]) | |
| ┌─────────────┬─────────────┬─────────────┬─────────────┬─────────────┬──────────────┬─────────────┬──────────────┬─────── | |
| │ STATEFP │ STATENS │ GEOIDFQ │ GEOID │ STUSPS │ NAME │ LSAD │ ALAND │ AW ⋯ | |
| │ Categorical │ Categorical │ Categorical │ Categorical │ Categorical │ Categorical │ Categorical │ Categorical │ Cate ⋯ | |
| │ [NoUnits] │ [NoUnits] │ [NoUnits] │ [NoUnits] │ [NoUnits] │ [NoUnits] │ [NoUnits] │ [NoUnits] │ [No ⋯ | |
| ├─────────────┼─────────────┼─────────────┼─────────────┼─────────────┼──────────────┼─────────────┼──────────────┼─────── | |
| │ 35 │ 00897535 │ 0400000US35 │ 35 │ NM │ New Mexico │ 00 │ 314198519809 │ 726 ⋯ | |
| │ ⋮ │ ⋮ │ ⋮ │ ⋮ │ ⋮ │ ⋮ │ ⋮ │ ⋮ │ ⋱ | |
| └─────────────┴─────────────┴─────────────┴─────────────┴─────────────┴──────────────┴─────────────┴──────────────┴─────── | |
| 2 columns and 48 rows omitted | |
| julia> @time viz!(ax, gt.geometry, color=1:nrow(gt), segmentcolor="black", showsegments=true, segmentsize=0.3f0) | |
| 23.725956 seconds (543.23 M allocations: 47.472 GiB, 18.22% gc time, 2.21% compilation time) | |
| Makie.Plot{viz, Tuple{SubDomain{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}, GeometrySet{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}, MultiPolygon{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}, PolyArea{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}, Ring{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}, CircularArrays.CircularVector{Point{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}}, Vector{Point{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}}}}}, Vector{Ring{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}, CircularArrays.CircularVector{Point{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}}, Vector{Point{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}}}}}}}}}, Vector{Int64}}}} | |
| julia> ax.aspect = Mke.DataAspect() | |
| Makie.DataAspect() | |
| julia> Mke.display(fig) | |
| # find out why Filter did not remove the overseas territories and Alaska | |
| # ?GeoTables | |
| # online search | |
| CairoMakie.Screen{IMAGE} | |
| julia> unique(gt.NAME) | |
| # shows 49 states | |
| 49-element Vector{Union{Missing, String}}: | |
| "New Mexico" | |
| "South Dakota" | |
| "California" | |
| ⋮ | |
| "South Carolina" | |
| "Maine" | |
| "Rhode Island" | |
| julia> unique(gt.geometry) | |
| # try new session to see if the problem persists | |
| # same result | |
| # fix projection | |
| 49-element Vector{MultiPolygon{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}, PolyArea{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}, Ring{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}, CircularArrays.CircularVector{Point{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}}, Vector{Point{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}}}}}, Vector{Ring{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}, CircularArrays.CircularVector{Point{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}}, Vector{Point{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}}}}}}}}}: | |
| Multi(1×PolyArea) | |
| Multi(1×PolyArea) | |
| Multi(15×PolyArea) | |
| ⋮ | |
| Multi(4×PolyArea) | |
| Multi(212×PolyArea) | |
| Multi(14×PolyArea) | |
| julia> Pkg.add("GeoMakie") | |
| Resolving package versions... | |
| Compat entries added for | |
| No Changes to `~/projects/testbed/Project.toml` | |
| No Changes to `~/projects/testbed/Manifest.toml` | |
| julia> using GeoMakie | |
| WARNING: using GeoMakie.Point in module Main conflicts with an existing identifier. | |
| julia> fig = Mke.Figure() | |
| julia> ax = GeoAxis(fig[1, 1], title = "US 'Lower 48' States",dest = "+proj=aea +lat_0=37.5 +lon_0=-96 +lat_1=29.5 +lat_2=45.5 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs") | |
| GeoAxis() | |
| julia> @time viz!(ax, gt.geometry, color=1:nrow(gt), segmentcolor="black", showsegments=true, segmentsize=0.3f0) | |
| 22.784561 seconds (540.96 M allocations: 47.349 GiB, 16.33% gc time, 0.62% compilation time) | |
| Plot{viz, Tuple{SubDomain{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}, GeometrySet{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}, MultiPolygon{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}, PolyArea{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}, Ring{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}, CircularArrays.CircularVector{Point{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}}, Vector{Point{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}}}}}, Vector{Ring{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}, CircularArrays.CircularVector{Point{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}}, Vector{Point{🌐, GeodeticLatLon{NAD83, Quantity{Float64, NoDims, Unitful.FreeUnits{(°,), NoDims, nothing}}}}}}}}}}}, Vector{Int64}}}} | |
| julia> ax.aspect = GeoMakie.DataAspect() | |
| # GeoMakie.display(fig) | |
| # display hanging; kill REPL | |
| # restart REPL | |
| DataAspect() | |
| julia> using Pkg | |
| julia> Pkg.activate(@__DIR__) | |
| Activating project at `~/projects/testbed` | |
| julia> Pkg.status() | |
| Project testbed v0.1.0 | |
| Status `~/projects/testbed/Project.toml` | |
| [336ed68f] CSV v0.10.15 | |
| [13f3f980] CairoMakie v0.15.4 | |
| [e9467ef8] GLMakie v0.13.4 | |
| [5752ebe1] GMT v1.31.0 | |
| [db073c08] GeoMakie v0.7.15 | |
| [dcc97b0b] GeoStats v0.81.2 | |
| [e502b557] GeoTables v1.24.19 | |
| [e6f89c97] LoggingExtras v1.1.0 | |
| [56ddb016] Logging v1.11.0 | |
| julia> versioninfo() | |
| Julia Version 1.11.6 | |
| Commit 9615af0f269 (2025-07-09 12:58 UTC) | |
| Build Info: | |
| Official https://julialang.org/ release | |
| Platform Info: | |
| OS: macOS (arm64-apple-darwin24.0.0) | |
| CPU: 28 × Apple M3 Ultra | |
| WORD_SIZE: 64 | |
| LLVM: libLLVM-16.0.6 (ORCJIT, apple-m2) | |
| Threads: 1 default, 0 interactive, 1 GC (on 20 virtual cores) | |
| Environment: | |
| JULIA_EDITOR = lvim | |
| JULIA_DEPOT_PATH = /Users/technocrat/.julia | |
| JULIA_PKG_USE_CLI_GIT = true | |
| JULIA_PROMPT = > | |
| DYLD_LIBRARY_PATH = /opt/homebrew/opt/[email protected]/lib: | |
| DYLD_FALLBACK_LIBRARY_PATH = /opt/homebrew/lib: | |
Author
technocrat
commented
Jul 27, 2025
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment