# Geo

**URL:** https://discourse.julialang.org/c/domain/geo/18.md

[Latest](https://discourse.julialang.org/latest.md) · [Categories](https://discourse.julialang.org/categories.md) · [Tags](https://discourse.julialang.org/tags.md)

---

## [About the Geo category](https://discourse.julialang.org/t/about-the-geo-category/241)

<div class="topic-metadata">

**Author:** [@StefanKarpinski](https://discourse.julialang.org/u/StefanKarpinski)\
**Replies:** 0

</div>

This subcategory is for geo-related discussions. This includes using Julia in geosciences, or dealing with geospatial data. Besides using Julia for the above, the development of an ecosystem of packages that allow typic…

---

## [Triangulate.jl - Generating a triangle grid using Triangulate](https://discourse.julialang.org/t/triangulate-jl-generating-a-triangle-grid-using-triangulate/139023)

<div class="topic-metadata">

**Author:** [@omgplshelp](https://discourse.julialang.org/u/omgplshelp)\
**Replies:** 3\
**Last updated:** [August 25, 2026, 3:26pm UTC](https://discourse.julialang.org/t/triangulate-jl-generating-a-triangle-grid-using-triangulate/139023 "2026-08-25T15:26:59Z")

</div>

Hi, I’m trying to generate an unstructured triangle grid using Triangulate.jl with coordinate data specifying layer boundaries and measuring points within those boundaries. I’ve tried just about every bit of example cod…

---

## [GMT.jl : Coloring lakes, keeping oceans bathymetry](https://discourse.julialang.org/t/gmt-jl-coloring-lakes-keeping-oceans-bathymetry/135163)

<div class="topic-metadata">

**Author:** [@Balinus](https://discourse.julialang.org/u/Balinus)\
**Replies:** 6\
**Last updated:** [January 22, 2026, 2:06pm UTC](https://discourse.julialang.org/t/gmt-jl-coloring-lakes-keeping-oceans-bathymetry/135163 "2026-01-22T14:06:22Z")

</div>

Hello! I’m trying to plot some points over a map. I would like to keep the default topography/bathymetry visuals, but the lakes are “green”, being over 0m at sea level (my guess). Here’s my current code. I tried to use t…

---

## [Plot density map of points in GeoMakie](https://discourse.julialang.org/t/plot-density-map-of-points-in-geomakie/135107)

<div class="topic-metadata">

**Author:** [@mihalybaci](https://discourse.julialang.org/u/mihalybaci)\
**Replies:** 4\
**Last updated:** [January 18, 2026, 2:01am UTC](https://discourse.julialang.org/t/plot-density-map-of-points-in-geomakie/135107 "2026-01-18T02:01:09Z")

</div>

I am trying to figure out how to plot a heat map / density map of geolocations (lon, lat) and overlay it on a geographical map. I think I am nearly there using GeoMakie.jl, but I just can’t get the final step. I am able …

---

## [Possible to geocode from file using GMT?](https://discourse.julialang.org/t/possible-to-geocode-from-file-using-gmt/135105)

<div class="topic-metadata">

**Author:** [@mihalybaci](https://discourse.julialang.org/u/mihalybaci)\
**Replies:** 5\
**Last updated:** [January 18, 2026, 1:51am UTC](https://discourse.julialang.org/t/possible-to-geocode-from-file-using-gmt/135105 "2026-01-18T01:51:08Z")

</div>

I am looking to gecode some address information (convert addresses to lon/lat pairs). I found the geodcode function in GMT.jl, and it looks simple enough. But, is it possible to avoid sending the addresses to a web servi…

---

## [GMT.jl: How to specify several river styles e.g. in coast?](https://discourse.julialang.org/t/gmt-jl-how-to-specify-several-river-styles-e-g-in-coast/133834)

<div class="topic-metadata">

**Author:** [@trg818](https://discourse.julialang.org/u/trg818)\
**Replies:** 3\
**Last updated:** [November 13, 2025, 1:20pm UTC](https://discourse.julialang.org/t/gmt-jl-how-to-specify-several-river-styles-e-g-in-coast/133834 "2025-11-13T13:20:10Z")

</div>

Hi, basically, my question is how to translate the regular GMT coast options -I0/thicker,lightblue -I1/thick,lightblue -I2/thin,lightblue to the keyword argument style. I tried rivers=(type=(0,1,2),style=("thicker","…

---

## [How to extract location from Strided GeoTIFF raster?](https://discourse.julialang.org/t/how-to-extract-location-from-strided-geotiff-raster/124973)

<div class="topic-metadata">

**Author:** [@cshen](https://discourse.julialang.org/u/cshen)\
**Replies:** 14\
**Last updated:** [October 30, 2025, 4:32pm UTC](https://discourse.julialang.org/t/how-to-extract-location-from-strided-geotiff-raster/124973 "2025-10-30T16:32:44Z")

</div>

Hello, I’ve been trying to read some geotiff datasets. Reading it with GeoTIFF.jl results in a GeoTIFF.GeoTIFFImageIterator{TiffImages.StridedTaggedImage{...}}}. This case is not really documented anywhere neither in G…

---

## [Area interpolation between Shapefile (country borders) and WKT (multi)polygon elements in a CSV file](https://discourse.julialang.org/t/area-interpolation-between-shapefile-country-borders-and-wkt-multi-polygon-elements-in-a-csv-file/133458)

<div class="topic-metadata">

**Author:** [@sylvaticus](https://discourse.julialang.org/u/sylvaticus)\
**Replies:** 15\
**Last updated:** [October 30, 2025, 1:01am UTC](https://discourse.julialang.org/t/area-interpolation-between-shapefile-country-borders-and-wkt-multi-polygon-elements-in-a-csv-file/133458 "2025-10-30T01:01:00Z")

</div>

I have a shapefile of administrative country borders and a CSV file where one field is “geometry” and embed polygon and multipolygon data as "Well Known Text (WKT). Both are in EPSG:4326. How can I retrieve for each ele…

---

## [Status of geodesy packages in Julia (DGGS experiments)](https://discourse.julialang.org/t/status-of-geodesy-packages-in-julia-dggs-experiments/133379)

<div class="topic-metadata">

**Author:** [@allixender](https://discourse.julialang.org/u/allixender)\
**Replies:** 12\
**Last updated:** [October 25, 2025, 10:27am UTC](https://discourse.julialang.org/t/status-of-geodesy-packages-in-julia-dggs-experiments/133379 "2025-10-25T10:27:13Z")

</div>

Hi all, I would like to instigate a short review discussion on geodesy-related packages in Julia. My background is GIS, geospatial, and recently Discrete Global Grid Systems (DGGS) in particular. If you have heard of Ub…

---

## [Heatmap! in GeoMakie](https://discourse.julialang.org/t/heatmap-in-geomakie/133249)

<div class="topic-metadata">

**Author:** [@mreichMPI-BGC](https://discourse.julialang.org/u/mreichMPI-BGC)\
**Replies:** 3\
**Last updated:** [October 18, 2025, 1:32pm UTC](https://discourse.julialang.org/t/heatmap-in-geomakie/133249 "2025-10-18T13:32:02Z")

</div>

The projection in heatmap! does not seem to work (using GeoMAkie 0.7.15), the example using GeoMakie, GLMakie lons = -180:180 lats = -90:90 field = \[exp(cosd(l)) + 3(y/90) for l in lons, y in lats\] fig = Figure() ax = G…

---

## [PooledArray.PooledVector causing StackOverflowError in GeoStats when using viz/plot, expected?](https://discourse.julialang.org/t/pooledarray-pooledvector-causing-stackoverflowerror-in-geostats-when-using-viz-plot-expected/132962)

<div class="topic-metadata">

**Author:** [@diadora77](https://discourse.julialang.org/u/diadora77)\
**Replies:** 8\
**Last updated:** [October 10, 2025, 12:11am UTC](https://discourse.julialang.org/t/pooledarray-pooledvector-causing-stackoverflowerror-in-geostats-when-using-viz-plot-expected/132962 "2025-10-10T00:11:44Z")

</div>

I ran into a situation that had me going in circle in that GeoStats viz/plot does not play well with PooledArray.PooledVector of anything(?). Situation was reading in a csv file into a dataframe then using GeoStats viz …

---

## [I want to build a hydrogeology model program](https://discourse.julialang.org/t/i-want-to-build-a-hydrogeology-model-program/132129)

<div class="topic-metadata">

**Author:** [@GeodeticR](https://discourse.julialang.org/u/GeodeticR)\
**Replies:** 14\
**Last updated:** [September 12, 2025, 9:51am UTC](https://discourse.julialang.org/t/i-want-to-build-a-hydrogeology-model-program/132129 "2025-09-12T09:51:08Z")

</div>

Hi everyone. I’m not really sure how to preface this. so I guess I’ll just provide my background, and my reason for wanting to try and build something useful in Julia for my usage. I recently finished my master’s degre…

---

## [Let's Choroplethize the US](https://discourse.julialang.org/t/lets-choroplethize-the-us/131254)

<div class="topic-metadata">

**Author:** [@joa-quim](https://discourse.julialang.org/u/joa-quim)\
**Replies:** 18\
**Last updated:** [August 31, 2025, 10:34pm UTC](https://discourse.julialang.org/t/lets-choroplethize-the-us/131254 "2025-08-31T22:34:53Z")

</div>

Hi, Making choropleth maps of US seems to a popular request/example here. Advanced solution often lack, IMO, the necessary map quality. In these last days before vacations I spent some time preparing these examples that…

---

## [How to extract the x,y coordinated otf the center of pixel 200, 300 of a Rasters.Raster object?](https://discourse.julialang.org/t/how-to-extract-the-x-y-coordinated-otf-the-center-of-pixel-200-300-of-a-rasters-raster-object/131936)

<div class="topic-metadata">

**Author:** [@sylvaticus](https://discourse.julialang.org/u/sylvaticus)\
**Replies:** 5\
**Last updated:** [August 30, 2025, 10:03am UTC](https://discourse.julialang.org/t/how-to-extract-the-x-y-coordinated-otf-the-center-of-pixel-200-300-of-a-rasters-raster-object/131936 "2025-08-30T10:03:42Z")

</div>

Provided I have a Rasters.Raster object foo, how can I extract the coordinated at the centre of a given pixel, in the raster coordinates system ?

---

## [\[NEW BOOK\] Geospatial Data Science with Julia](https://discourse.julialang.org/t/new-book-geospatial-data-science-with-julia/103364)

<div class="topic-metadata">

**Author:** [@juliohm](https://discourse.julialang.org/u/juliohm)\
**Replies:** 13\
**Last updated:** [August 8, 2025, 12:15pm UTC](https://discourse.julialang.org/t/new-book-geospatial-data-science-with-julia/103364 "2025-08-08T12:15:13Z")

</div>

The book will be published soon, and will contain practical applications of the framework we built over the years. If you would like to see an application in a specific area (e.g., agriculture, groundwater, environmen…

---

## [Plotting in GeoStats](https://discourse.julialang.org/t/plotting-in-geostats/130770)

<div class="topic-metadata">

**Author:** [@beppi](https://discourse.julialang.org/u/beppi)\
**Replies:** 27\
**Last updated:** [July 29, 2025, 8:39pm UTC](https://discourse.julialang.org/t/plotting-in-geostats/130770 "2025-07-29T20:39:52Z")

</div>

Hey, I’ve seen the page here: Visualization · GeoStats.jl However as a newbie I might be unable to understand how to do this. I have a GeoTable of points and I want to plot the points coloured by a column in the GeoTab…

---

## [Add new column to GeoTable from Vector](https://discourse.julialang.org/t/add-new-column-to-geotable-from-vector/130735)

<div class="topic-metadata">

**Author:** [@beppi](https://discourse.julialang.org/u/beppi)\
**Replies:** 3\
**Last updated:** [July 16, 2025, 2:53am UTC](https://discourse.julialang.org/t/add-new-column-to-geotable-from-vector/130735 "2025-07-16T02:53:16Z")

</div>

Hey guys, new to Julia I’ve searched so much and can’t seem to find an answer to this simple question.. How do I take a Vector, and add it as a new column to a GeoTable? I see there is a way to create a new column (@t…

---

## [How to get value and matrix index at a given (x;y) point?](https://discourse.julialang.org/t/how-to-get-value-and-matrix-index-at-a-given-x-y-point/129692)

<div class="topic-metadata">

**Author:** [@sylvaticus](https://discourse.julialang.org/u/sylvaticus)\
**Replies:** 9\
**Last updated:** [June 27, 2025, 5:04am UTC](https://discourse.julialang.org/t/how-to-get-value-and-matrix-index-at-a-given-x-y-point/129692 "2025-06-27T05:04:10Z")

</div>

I have a set of rasters with the same grid and a large number of points and need to create a CZV with these points + information from the rasters (so my questions is similar to this one, but I am using the Rasters.jl pac…

---

## [GeoMakie GeoAxis Subticks?](https://discourse.julialang.org/t/geomakie-geoaxis-subticks/130024)

<div class="topic-metadata">

**Author:** [@consumere](https://discourse.julialang.org/u/consumere)\
**Replies:** 2\
**Last updated:** [June 19, 2025, 1:16pm UTC](https://discourse.julialang.org/t/geomakie-geoaxis-subticks/130024 "2025-06-19T13:16:34Z")

</div>

If i run this code, i get wrong axis labels. using GeoMakie set\_theme!(theme\_latexfonts()) fig = Figure(fontsize=12) lonlims = (10,11) latlims = (49,51) ax = GeoAxis(fig\[1, 1\], limits=(lonlims, latlims)) fig Q…

---

## [Checking that a position is inside a AREA works in EPSG:25382 projection, but not in CRS84](https://discourse.julialang.org/t/checking-that-a-position-is-inside-a-area-works-in-epsg-25382-projection-but-not-in-crs84/129471)

<div class="topic-metadata">

**Author:** [@sfuerst](https://discourse.julialang.org/u/sfuerst)\
**Replies:** 1\
**Last updated:** [May 30, 2025, 11:52am UTC](https://discourse.julialang.org/t/checking-that-a-position-is-inside-a-area-works-in-epsg-25382-projection-but-not-in-crs84/129471 "2025-05-30T11:52:00Z")

</div>

I have first to say, that I’m a newbie when it comes to geo-related stuff, so please excuse me if I ask or formulate something stupid here, but I’ve been biting my teeth out on this problem for hours. First a little con…

---

## [Plotting a georeferenced field on a 3d globe](https://discourse.julialang.org/t/plotting-a-georeferenced-field-on-a-3d-globe/128895)

<div class="topic-metadata">

**Author:** [@aasdelat](https://discourse.julialang.org/u/aasdelat)\
**Replies:** 5\
**Last updated:** [May 17, 2025, 4:25pm UTC](https://discourse.julialang.org/t/plotting-a-georeferenced-field-on-a-3d-globe/128895 "2025-05-17T16:25:53Z")

</div>

I have georeferenced data in longitude/latitude, and I want to plot them on a three-dimensional globe with the WGS84 datum, using the GLMakie and GeoMakie libraries. You can download the data from here. The grid is irr…

---

## [Satellogic data in pure Julia](https://discourse.julialang.org/t/satellogic-data-in-pure-julia/127909)

<div class="topic-metadata">

**Author:** [@Zgabalj](https://discourse.julialang.org/u/Zgabalj)\
**Replies:** 3\
**Last updated:** [April 10, 2025, 11:41am UTC](https://discourse.julialang.org/t/satellogic-data-in-pure-julia/127909 "2025-04-10T11:41:33Z")

</div>

Mark Litwintschik has a blog (cannot add link) showing how to obtain satellite images from Satellogic using python. How would one start to implement this in Julia? The author uses GDAL and duckdb; both GDAL.jl and Duck…

---

## [\[ANN\] GeoIDs.jl: a package to create a postgres/postgis database of TIGER county boundaries](https://discourse.julialang.org/t/ann-geoids-jl-a-package-to-create-a-postgres-postgis-database-of-tiger-county-boundaries/127778)

<div class="topic-metadata">

**Author:** [@technocrat](https://discourse.julialang.org/u/technocrat)\
**Replies:** 3\
**Last updated:** [April 7, 2025, 10:40am UTC](https://discourse.julialang.org/t/ann-geoids-jl-a-package-to-create-a-postgres-postgis-database-of-tiger-county-boundaries/127778 "2025-04-07T10:40:40Z")

</div>

GeoIDs.jl provides a comprehensive solution for working with Census GEOID identifiers, particularly at the county level. The package enables you to: Create and Manage Geographic Area Sets: Define sets of counties for …

---

## [Best library for calculating slope in an elevation raster?](https://discourse.julialang.org/t/best-library-for-calculating-slope-in-an-elevation-raster/114154)

<div class="topic-metadata">

**Author:** [@rholland4julia](https://discourse.julialang.org/u/rholland4julia)\
**Replies:** 13\
**Last updated:** [April 6, 2025, 4:14pm UTC](https://discourse.julialang.org/t/best-library-for-calculating-slope-in-an-elevation-raster/114154 "2025-04-06T16:14:46Z")

</div>

I am trying to calculate the slope from a large number of DEM rasters. It seems like Geomorphometry.jl has an implementation of the Horn (1981) method, but are there other library options?

---

## [GeoTiff files retrieving metadata by using Rasters.jl as in rasterio.py](https://discourse.julialang.org/t/geotiff-files-retrieving-metadata-by-using-rasters-jl-as-in-rasterio-py/127587)

<div class="topic-metadata">

**Author:** [@JosephPollacco](https://discourse.julialang.org/u/JosephPollacco)\
**Replies:** 12\
**Last updated:** [April 2, 2025, 8:45am UTC](https://discourse.julialang.org/t/geotiff-files-retrieving-metadata-by-using-rasters-jl-as-in-rasterio-py/127587 "2025-04-02T08:45:44Z")

</div>

Hi, I am looking how to retrive metadata from GeoTiff files by using Julia language. I can obtain the relevant information by using the rasterio.jl python package, but I am struggling to perform the same commands by u…

---

## [Retrieving CRS from vector and raster data](https://discourse.julialang.org/t/retrieving-crs-from-vector-and-raster-data/107903)

<div class="topic-metadata">

**Author:** [@Bojas](https://discourse.julialang.org/u/Bojas)\
**Replies:** 7\
**Last updated:** [April 1, 2025, 2:15pm UTC](https://discourse.julialang.org/t/retrieving-crs-from-vector-and-raster-data/107903 "2025-04-01T14:15:00Z")

</div>

Hello, I’m new to Julia, as for now I was using mainly Python. What I want to do is to crop part of raster based on polygon. Both polygon and raster are in different spatial references and I have to convert either one of…

---

## [FE Discretisation for embedded cables](https://discourse.julialang.org/t/fe-discretisation-for-embedded-cables/127025)

<div class="topic-metadata">

**Author:** [@Br0BOLE](https://discourse.julialang.org/u/Br0BOLE)\
**Replies:** 13\
**Last updated:** [March 18, 2025, 10:12pm UTC](https://discourse.julialang.org/t/fe-discretisation-for-embedded-cables/127025 "2025-03-18T22:12:38Z")

</div>

Hi there, I’m trying to automate the discritization of FE meshes for High voltage cables embedded in the ground. This is to conduct thermal analysis and parameterise the geometry of the cable (cables with tri cores). I…

---

## [Geostatistics in Julia](https://discourse.julialang.org/t/geostatistics-in-julia/126862)

<div class="topic-metadata">

**Author:** [@juliohm](https://discourse.julialang.org/u/juliohm)\
**Replies:** 0\
**Last updated:** [March 12, 2025, 10:59am UTC](https://discourse.julialang.org/t/geostatistics-in-julia/126862 "2025-03-12T10:59:34Z")

</div>

In case you are not following the latest updates in our discourse thread, here is a comparison table that you can share with colleagues when they ask about the status of geo + stats software in Julia:

---

## [How to load a raster of a NC file with multiple monthly bands](https://discourse.julialang.org/t/how-to-load-a-raster-of-a-nc-file-with-multiple-monthly-bands/126487)

<div class="topic-metadata">

**Author:** [@sylvaticus](https://discourse.julialang.org/u/sylvaticus)\
**Replies:** 4\
**Last updated:** [March 3, 2025, 9:18am UTC](https://discourse.julialang.org/t/how-to-load-a-raster-of-a-nc-file-with-multiple-monthly-bands/126487 "2025-03-03T09:18:30Z")

</div>

I am using the chelsa\_cmip6 script to download climatic data, and the output is a NC file with the data embedded in monthly series, like in this (small) example of precipitations. How do I load it into a vector of raste…

---

## [Can't get the expected resolution using Rasters.resample](https://discourse.julialang.org/t/cant-get-the-expected-resolution-using-rasters-resample/126118)

<div class="topic-metadata">

**Author:** [@MJulia](https://discourse.julialang.org/u/MJulia)\
**Replies:** 2\
**Last updated:** [February 21, 2025, 2:25am UTC](https://discourse.julialang.org/t/cant-get-the-expected-resolution-using-rasters-resample/126118 "2025-02-21T02:25:04Z")

</div>

I tried the following code snippet to resample high-resolution data to a lower resolution. The first method works, but the second one doesn’t. It seems I missed some parameters for the “to” keyword in the resample functi…

[Next page](https://discourse.julialang.org/c/domain/geo/18.md?page=1)