# Geo

**URL:** https://discourse.julialang.org/c/domain/geo/18.md?page=1

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

**Page:** 2

---

## [Reprojecting maps](https://discourse.julialang.org/t/reprojecting-maps/97101)

<div class="topic-metadata">

**Author:** [@mreichMPI-BGC](https://discourse.julialang.org/u/mreichMPI-BGC)\
**Replies:** 14\
**Last updated:** [February 10, 2025, 8:47am UTC](https://discourse.julialang.org/t/reprojecting-maps/97101 "2025-02-10T08:47:57Z")

</div>

Dear community, does anyone have a recipe/example of how to reproject a map, e.g. from lat/lon to Robinson (for my use case just global is enough)? Please note I do not just want to plot (cf. GeoMakie), but get a Raster…

---

## [Add point coordinates to a GeoTable](https://discourse.julialang.org/t/add-point-coordinates-to-a-geotable/125126)

<div class="topic-metadata">

**Author:** [@phchavesmaia](https://discourse.julialang.org/u/phchavesmaia)\
**Replies:** 12\
**Last updated:** [February 7, 2025, 3:05pm UTC](https://discourse.julialang.org/t/add-point-coordinates-to-a-geotable/125126 "2025-02-07T15:05:51Z")

</div>

Hello everyone, I am trying to learn Julia and move some of my geostats projects to this language. I found that the GeoStats package is well regarded, hence I am using it. I have a grid and I want to separate the east f…

---

## [GeoStats.jl 2025 survey](https://discourse.julialang.org/t/geostats-jl-2025-survey/124362)

<div class="topic-metadata">

**Author:** [@juliohm](https://discourse.julialang.org/u/juliohm)\
**Replies:** 0\
**Last updated:** [January 2, 2025, 12:21pm UTC](https://discourse.julialang.org/t/geostats-jl-2025-survey/124362 "2025-01-02T12:21:04Z")

</div>

---

## [Intersection of Ball and GeometrySet in LatLon crs](https://discourse.julialang.org/t/intersection-of-ball-and-geometryset-in-latlon-crs/123992)

<div class="topic-metadata">

**Author:** [@stepanoslejsek](https://discourse.julialang.org/u/stepanoslejsek)\
**Replies:** 6\
**Last updated:** [December 20, 2024, 8:37am UTC](https://discourse.julialang.org/t/intersection-of-ball-and-geometryset-in-latlon-crs/123992 "2024-12-20T08:37:26Z")

</div>

Hi, I’d like to find an intersection of a Ball with GeometrySet in LatLon coordinate reference system. However, when I construct my Ball b like b = GeoStats.Ball(GeoStats.Point(LatLon(48.34, 18.45)), 10) and ropes r a…

---

## [Changing EPSG values from .shp file](https://discourse.julialang.org/t/changing-epsg-values-from-shp-file/122978)

<div class="topic-metadata">

**Author:** [@ivyc](https://discourse.julialang.org/u/ivyc)\
**Replies:** 7\
**Last updated:** [December 3, 2024, 3:42pm UTC](https://discourse.julialang.org/t/changing-epsg-values-from-shp-file/122978 "2024-12-03T15:42:00Z")

</div>

I am trying to convert my .shp file from an EPSG value of 26917 to a EPSG value of 4326. These seems to be an issue with my .shp file not having the EPSG built in. I have tried changing the EPSG using the code below but …

---

## [Spatial Join with GeometryOps is killing my Julia session](https://discourse.julialang.org/t/spatial-join-with-geometryops-is-killing-my-julia-session/122863)

<div class="topic-metadata">

**Author:** [@alejandromerchan](https://discourse.julialang.org/u/alejandromerchan)\
**Replies:** 15\
**Last updated:** [November 25, 2024, 4:13pm UTC](https://discourse.julialang.org/t/spatial-join-with-geometryops-is-killing-my-julia-session/122863 "2024-11-25T16:13:03Z")

</div>

I’m trying to perform a spatial join in Julia between two files that aren’t that large. One is the California PLSS sections shapefile, which is publicly available here and is 164627 x 17 after some cleaning. The other I …

---

## [GeoStats.jl workshop at FOSS4G 2024](https://discourse.julialang.org/t/geostats-jl-workshop-at-foss4g-2024/122311)

<div class="topic-metadata">

**Author:** [@juliohm](https://discourse.julialang.org/u/juliohm)\
**Replies:** 0\
**Last updated:** [November 6, 2024, 12:36am UTC](https://discourse.julialang.org/t/geostats-jl-workshop-at-foss4g-2024/122311 "2024-11-06T00:36:20Z")

</div>

Upcoming workshop at FOSS4G in Brazil this year:

---

## [Screening Interest in a unified Vector/Raster Package, akin to R Terra](https://discourse.julialang.org/t/screening-interest-in-a-unified-vector-raster-package-akin-to-r-terra/121613)

<div class="topic-metadata">

**Author:** [@souma](https://discourse.julialang.org/u/souma)\
**Replies:** 27\
**Last updated:** [October 26, 2024, 1:36pm UTC](https://discourse.julialang.org/t/screening-interest-in-a-unified-vector-raster-package-akin-to-r-terra/121613 "2024-10-26T13:36:33Z")

</div>

Hello all, I haven’t been active in Discourse but I’ve been working with Julia for scientific computing for quite a few years at this point. I want to know what the interest is in a JuliaGeo project/effort to harmonize…

---

## [DBSCAN clustering with Haversine metric](https://discourse.julialang.org/t/dbscan-clustering-with-haversine-metric/121685)

<div class="topic-metadata">

**Author:** [@miguelborrero](https://discourse.julialang.org/u/miguelborrero)\
**Replies:** 4\
**Last updated:** [October 24, 2024, 12:53pm UTC](https://discourse.julialang.org/t/dbscan-clustering-with-haversine-metric/121685 "2024-10-24T12:53:46Z")

</div>

Hi there, I have an array of points (latitude, longitude) in radians and I am trying to use the package Clustering.jl to cluster my points. In the documentation it seems like a metric can be given, the default being: Eu…

---

## [Recommended workflow and packages for computing a slope along a tramway track](https://discourse.julialang.org/t/recommended-workflow-and-packages-for-computing-a-slope-along-a-tramway-track/119975)

<div class="topic-metadata">

**Author:** [@zdenek\_hurak](https://discourse.julialang.org/u/zdenek_hurak)\
**Replies:** 39\
**Last updated:** [October 24, 2024, 9:00am UTC](https://discourse.julialang.org/t/recommended-workflow-and-packages-for-computing-a-slope-along-a-tramway-track/119975 "2024-10-24T09:00:37Z")

</div>

I would like to ask for a recommendation on a workflow for computing a slope along a tramway track using Julia’s geo ecosystem (organizations JuliaEarth, JuliaGeo, …). The tramway tracks are obtained from OpenStreetMaps…

---

## [Wavelet Transformation](https://discourse.julialang.org/t/wavelet-transformation/120737)

<div class="topic-metadata">

**Author:** [@Pradyumnya\_Dastidar](https://discourse.julialang.org/u/Pradyumnya_Dastidar)\
**Replies:** 4\
**Last updated:** [October 22, 2024, 1:47pm UTC](https://discourse.julialang.org/t/wavelet-transformation/120737 "2024-10-22T13:47:43Z")

</div>

Can someone please help me to write a julia code on Wavelet Transformation for a seismic signal.

---

## [Sampling random points in LINESTRING geometry](https://discourse.julialang.org/t/sampling-random-points-in-linestring-geometry/121429)

<div class="topic-metadata">

**Author:** [@miguelborrero](https://discourse.julialang.org/u/miguelborrero)\
**Replies:** 7\
**Last updated:** [October 18, 2024, 12:07am UTC](https://discourse.julialang.org/t/sampling-random-points-in-linestring-geometry/121429 "2024-10-18T00:07:33Z")

</div>

Hi there, I a new user of the geo-tools in general, not only within the Julia ecosystem and I am trying to sample randomly from a LINESTRING LibGEOS geometry but the only function that seemed potentially useful is point…

---

## [Issue with gridlines when cliping is used in GMT's grdimage](https://discourse.julialang.org/t/issue-with-gridlines-when-cliping-is-used-in-gmts-grdimage/120844)

<div class="topic-metadata">

**Author:** [@k945](https://discourse.julialang.org/u/k945)\
**Replies:** 0\
**Last updated:** [October 3, 2024, 8:50am UTC](https://discourse.julialang.org/t/issue-with-gridlines-when-cliping-is-used-in-gmts-grdimage/120844 "2024-10-03T08:50:20Z")

</div>

How do I mask out the ocean when using GMT? Doing this using grdimage does not appear to be working. Here is the image without clipping applied However when I set cliping in grdimage to :land I get the following …

---

## [How to map rectangular lat/lon region on non-rectangular projection with GeoMakie?](https://discourse.julialang.org/t/how-to-map-rectangular-lat-lon-region-on-non-rectangular-projection-with-geomakie/103140)

<div class="topic-metadata">

**Author:** [@haffner](https://discourse.julialang.org/u/haffner)\
**Replies:** 5\
**Last updated:** [September 30, 2024, 4:06pm UTC](https://discourse.julialang.org/t/how-to-map-rectangular-lat-lon-region-on-non-rectangular-projection-with-geomakie/103140 "2024-09-30T16:06:04Z")

</div>

I’m trying to map a rectangular lat/lon polygon on a non-rectangular map projection using GeoMakie, When I use four corner points to define the region, GeoMakie draws straight lines between them. The example below show…

---

## [Efficiently check if points are contained in polygons](https://discourse.julialang.org/t/efficiently-check-if-points-are-contained-in-polygons/74415)

<div class="topic-metadata">

**Author:** [@ahmoreira](https://discourse.julialang.org/u/ahmoreira)\
**Replies:** 33\
**Last updated:** [September 27, 2024, 9:15pm UTC](https://discourse.julialang.org/t/efficiently-check-if-points-are-contained-in-polygons/74415 "2024-09-27T21:15:34Z")

</div>

Hi all, I am new to Julia and I am trying to check if points (plant occurrence points) fall within polygons from a shapefile. The code I wrote works, but I would like to make it more efficient. Currently, the code takes …

---

## [Solving Geoscience questions using Julia:](https://discourse.julialang.org/t/solving-geoscience-questions-using-julia/119262)

<div class="topic-metadata">

**Author:** [@GeodeticR](https://discourse.julialang.org/u/GeodeticR)\
**Replies:** 2\
**Last updated:** [September 11, 2024, 9:08am UTC](https://discourse.julialang.org/t/solving-geoscience-questions-using-julia/119262 "2024-09-11T09:08:33Z")

</div>

Okay, so, I’m trying to become proficient at using Julia for various applications and not just cleaning and crunching datasets. One of my graduate classes assigned a few problems that we can do by hand, but we are also …

---

## ["Raster" syntax in GeoStats.jl](https://discourse.julialang.org/t/raster-syntax-in-geostats-jl/119214)

<div class="topic-metadata">

**Author:** [@juliohm](https://discourse.julialang.org/u/juliohm)\
**Replies:** 0\
**Last updated:** [September 9, 2024, 11:07am UTC](https://discourse.julialang.org/t/raster-syntax-in-geostats-jl/119214 "2024-09-09T11:07:09Z")

</div>

Given that most people in the geo community are not aware of the “raster” syntax in GeoStats.jl, I collect a few examples here to facilitate the transition. This post may become part of the official project documentation…

---

## [Geocoding in Julia? (getting latitude and longitude from addresses)](https://discourse.julialang.org/t/geocoding-in-julia-getting-latitude-and-longitude-from-addresses/118814)

<div class="topic-metadata">

**Author:** [@indymnv](https://discourse.julialang.org/u/indymnv)\
**Replies:** 5\
**Last updated:** [September 2, 2024, 1:10pm UTC](https://discourse.julialang.org/t/geocoding-in-julia-getting-latitude-and-longitude-from-addresses/118814 "2024-09-02T13:10:34Z")

</div>

I have a dataframe with several addresses and I want to extract the latitude and longitude potentially using OpenStreetMap, is there a package that help with that?

---

## [NCDatasets: performance issues deriving fast Minimum Maximum value as too slow](https://discourse.julialang.org/t/ncdatasets-performance-issues-deriving-fast-minimum-maximum-value-as-too-slow/118468)

<div class="topic-metadata">

**Author:** [@JosephPollacco](https://discourse.julialang.org/u/JosephPollacco)\
**Replies:** 12\
**Last updated:** [August 23, 2024, 4:14am UTC](https://discourse.julialang.org/t/ncdatasets-performance-issues-deriving-fast-minimum-maximum-value-as-too-slow/118468 "2024-08-23T04:14:54Z")

</div>

:turtle: :turtle: :turtle: :turtle: Finding minimum and maximum of a NetCDF files due to missing values is painfully slow. It takes me 2 minutes to perform this simple task: "NameOutput="SoilWater" Output\_NCDatasets =…

---

## [Shapefile.jl writing .shp without coordinates](https://discourse.julialang.org/t/shapefile-jl-writing-shp-without-coordinates/117434)

<div class="topic-metadata">

**Author:** [@alejandromerchan](https://discourse.julialang.org/u/alejandromerchan)\
**Replies:** 19\
**Last updated:** [August 4, 2024, 4:08pm UTC](https://discourse.julialang.org/t/shapefile-jl-writing-shp-without-coordinates/117434 "2024-08-04T16:08:36Z")

</div>

I loaded a Shapefile in Julia, transformed it to a datafame and then did some data cleaning, getting rid of some unneeded rows of data. I saved the shapefile again, with a new name, using “Shapefile.write” and tried to o…

---

## [Troubles with GeoStats/GeoIO coordinates with units](https://discourse.julialang.org/t/troubles-with-geostats-geoio-coordinates-with-units/117766)

<div class="topic-metadata">

**Author:** [@TimG](https://discourse.julialang.org/u/TimG)\
**Replies:** 5\
**Last updated:** [August 1, 2024, 10:39am UTC](https://discourse.julialang.org/t/troubles-with-geostats-geoio-coordinates-with-units/117766 "2024-08-01T10:39:09Z")

</div>

Thank you @juliohm! File reads fine now. New issue - too small to warrant new thread, perhaps. bb = boundingbox(gt.geometry) minx, miny = coordinates(bb.min)\[1:2\] coordinates function no longer works. ERROR: LoadErro…

---

## [ESRI code for British National Grid not known by GeoIO](https://discourse.julialang.org/t/esri-code-for-british-national-grid-not-known-by-geoio/117641)

<div class="topic-metadata">

**Author:** [@TimG](https://discourse.julialang.org/u/TimG)\
**Replies:** 22\
**Last updated:** [July 31, 2024, 2:19pm UTC](https://discourse.julialang.org/t/esri-code-for-british-national-grid-not-known-by-geoio/117641 "2024-07-31T14:19:56Z")

</div>

ERROR: LoadError: ArgumentError: EPSG/ESRI code for the ESRI ID "British\_National\_Grid" not found in dictionary. Please check https://github.com/JuliaEarth/CoordRefSystems.jl/blob/main/src/strings.jl If you know the EPSG…

---

## [Rasters: can I select certain time steps of a lazy rasterseries before or while reading it into memory?](https://discourse.julialang.org/t/rasters-can-i-select-certain-time-steps-of-a-lazy-rasterseries-before-or-while-reading-it-into-memory/116328)

<div class="topic-metadata">

**Author:** [@aasdelat](https://discourse.julialang.org/u/aasdelat)\
**Replies:** 6\
**Last updated:** [July 20, 2024, 8:38pm UTC](https://discourse.julialang.org/t/rasters-can-i-select-certain-time-steps-of-a-lazy-rasterseries-before-or-while-reading-it-into-memory/116328 "2024-07-20T20:38:28Z")

</div>

I have a bunch of grib files, each covering a different time period, but the same spatial extent. I have opened them as a RasterSeries (library Rasters.jl) with lazy=true, so that the data are not loaded into memory. da…

---

## [DimensionalData: How to use a dimension whose name is stored in a variable](https://discourse.julialang.org/t/dimensionaldata-how-to-use-a-dimension-whose-name-is-stored-in-a-variable/116098)

<div class="topic-metadata">

**Author:** [@aasdelat](https://discourse.julialang.org/u/aasdelat)\
**Replies:** 6\
**Last updated:** [June 27, 2024, 9:47pm UTC](https://discourse.julialang.org/t/dimensionaldata-how-to-use-a-dimension-whose-name-is-stored-in-a-variable/116098 "2024-06-27T21:47:51Z")

</div>

I have to use data with one time dimension and two spatial dimensions (longitude and latitude). I use the Rasters.jl library to manage them, that creates DimensionalData.jl arrays. Sometimes, I have a variable, say slp,…

---

## [ROMS.jl different error launching .in with romsS vs romsG](https://discourse.julialang.org/t/roms-jl-different-error-launching-in-with-romss-vs-romsg/115966)

<div class="topic-metadata">

**Author:** [@Unai\_Ganzedo](https://discourse.julialang.org/u/Unai_Ganzedo)\
**Replies:** 1\
**Last updated:** [June 22, 2024, 9:49pm UTC](https://discourse.julialang.org/t/roms-jl-different-error-launching-in-with-romss-vs-romsg/115966 "2024-06-22T21:49:44Z")

</div>

Hi; I have launched the lingurean example with romeS and with romsG and the questions are: Because it shows me different errors The solution to any of them to finish the example. Thank you Unai romsS: GET\_STATE\_NF…

---

## [How to transform CRS of entire shapefile (not a single point)?](https://discourse.julialang.org/t/how-to-transform-crs-of-entire-shapefile-not-a-single-point/84489)

<div class="topic-metadata">

**Author:** [@floswald](https://discourse.julialang.org/u/floswald)\
**Replies:** 6\
**Last updated:** [June 18, 2024, 6:45pm UTC](https://discourse.julialang.org/t/how-to-transform-crs-of-entire-shapefile-not-a-single-point/84489 "2024-06-18T18:45:04Z")

</div>

I could not find a way to transform an entire shapefile into another crs. There are too many options for a first time user. ArchGDAL.jl? Proj4? To me it seems reasonable that the Shapefile package should have a transform …

---

## [Polygonize a Raster](https://discourse.julialang.org/t/polygonize-a-raster/114189)

<div class="topic-metadata">

**Author:** [@consumere](https://discourse.julialang.org/u/consumere)\
**Replies:** 11\
**Last updated:** [June 18, 2024, 6:42pm UTC](https://discourse.julialang.org/t/polygonize-a-raster/114189 "2024-06-18T18:42:53Z")

</div>

Hello there, i know there is ArchGDAL.polygonize|\>methods \[1\] polygonize(f::Function, args...; kwargs...) @ \[2\] polygonize(geom::ArchGDAL.AbstractGeometry) and one could use PyCall or RCall to polygonize a ras…

---

## [ROMS.jl error with invalid dimension parameter](https://discourse.julialang.org/t/roms-jl-error-with-invalid-dimension-parameter/115468)

<div class="topic-metadata">

**Author:** [@Unai\_Ganzedo](https://discourse.julialang.org/u/Unai_Ganzedo)\
**Replies:** 1\
**Last updated:** [June 11, 2024, 11:49am UTC](https://discourse.julialang.org/t/roms-jl-error-with-invalid-dimension-parameter/115468 "2024-06-11T11:49:38Z")

</div>

Hi; I have a question regarding roms.in. I have followed the instructions and I have all the files prepared but when I launch the roms, I get this error..I have tried to search for it but I can’t find it, could you send …

---

## [How to rasterize a geopackage (Corine Land Cover)?](https://discourse.julialang.org/t/how-to-rasterize-a-geopackage-corine-land-cover/115316)

<div class="topic-metadata">

**Author:** [@sylvaticus](https://discourse.julialang.org/u/sylvaticus)\
**Replies:** 24\
**Last updated:** [June 8, 2024, 7:24pm UTC](https://discourse.julialang.org/t/how-to-rasterize-a-geopackage-corine-land-cover/115316 "2024-06-08T19:24:42Z")

</div>

Hello, being only marginally on the domain, I am a bit lost with the various GIS packages… I need to rasterize a geopackage (gpkg) vector file (Corine Land Cover). Visualizing the gpkg file in QGis I can see there are …

---

## [How to use Rasterize in GeoStats.jl](https://discourse.julialang.org/t/how-to-use-rasterize-in-geostats-jl/115338)

<div class="topic-metadata">

**Author:** [@sylvaticus](https://discourse.julialang.org/u/sylvaticus)\
**Replies:** 7\
**Last updated:** [June 7, 2024, 11:47pm UTC](https://discourse.julialang.org/t/how-to-use-rasterize-in-geostats-jl/115338 "2024-06-07T23:47:02Z")

</div>

(note: I am splitting this thread from this one in order to avoid misunderstanding) I can confirm that the latest version of GeoIO can now successfully import a geopackage where the geometries are provided in a Shape co…

[Previous page](https://discourse.julialang.org/c/domain/geo/18.md)

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