# Geo **URL:** https://discourse.julialang.org/c/domain/geo/18.md?page=4 [Latest](https://discourse.julialang.org/latest.md) · [Categories](https://discourse.julialang.org/categories.md) · [Tags](https://discourse.julialang.org/tags.md) **Page:** 5 --- ## [Plot Raster depth vs value](https://discourse.julialang.org/t/plot-raster-depth-vs-value/89150)
**Author:** [@ryofurue](https://discourse.julialang.org/u/ryofurue)\ **Replies:** 8\ **Last updated:** [November 5, 2022, 10:44am UTC](https://discourse.julialang.org/t/plot-raster-depth-vs-value/89150 "2022-11-05T10:44:27Z")
I’m trying to plot a vertical profile, a 1D plot, of a variable, a multidimensional Raster on a depth-time grid: using Rasters using Plots infile = "sample.nc" varname = "myvar" myvar = Raster(infile; name=myvar) # 2D… --- ## [Plotting data with 2D lat/lons in geomakie](https://discourse.julialang.org/t/plotting-data-with-2d-lat-lons-in-geomakie/88648)
**Author:** [@k945](https://discourse.julialang.org/u/k945)\ **Replies:** 8\ **Last updated:** [November 2, 2022, 9:37am UTC](https://discourse.julialang.org/t/plotting-data-with-2d-lat-lons-in-geomakie/88648 "2022-11-02T09:37:24Z")
Is it possible to plot data with 2d lat/lons using GeoMakie? I have so far been unable to properly project WRF simulation data onto a map. Using test data, this is what I have tried so far: using GeoMakie, CairoMa… --- ## [LightOSM and OpenStreetMapX read different nodes from one and the same osm file](https://discourse.julialang.org/t/lightosm-and-openstreetmapx-read-different-nodes-from-one-and-the-same-osm-file/89628)
**Author:** [@ctrebbau](https://discourse.julialang.org/u/ctrebbau)\ **Replies:** 0\ **Last updated:** [November 1, 2022, 7:13pm UTC](https://discourse.julialang.org/t/lightosm-and-openstreetmapx-read-different-nodes-from-one-and-the-same-osm-file/89628 "2022-11-01T19:13:50Z")
Anyone have a quick idea why this fails? g1 = LightOSM.graph\_from\_file( data\_file, ... largest\_connected\_component=true) g2 = OpenStreetMapX.get\_map\_data( … --- ## [LightOSM: precomputing dijkstra\_states throws an UndefVarError](https://discourse.julialang.org/t/lightosm-precomputing-dijkstra-states-throws-an-undefvarerror/89551)
**Author:** [@ctrebbau](https://discourse.julialang.org/u/ctrebbau)\ **Replies:** 1\ **Last updated:** [October 31, 2022, 10:26am UTC](https://discourse.julialang.org/t/lightosm-precomputing-dijkstra-states-throws-an-undefvarerror/89551 "2022-10-31T10:26:39Z")
From LightOSM.jl/benchmarks.jl at master · DeloitteOptimalReality/LightOSM.jl · GitHub point = GeoLocation(-37.8142176, 144.9631608) # Melbourne, Australia radius = 5 # km data\_file = joinpath(@\_\_DIR\_\_, "benchmark\_map.o… --- ## [Repartition data back and forward from one administrative division to another](https://discourse.julialang.org/t/repartition-data-back-and-forward-from-one-administrative-division-to-another/88389)
**Author:** [@sylvaticus](https://discourse.julialang.org/u/sylvaticus)\ **Replies:** 1\ **Last updated:** [October 7, 2022, 11:55am UTC](https://discourse.julialang.org/t/repartition-data-back-and-forward-from-one-administrative-division-to-another/88389 "2022-10-07T11:55:17Z")
I have some data relative to a certain administrative division, let’s say forest volumes per region. I need to translate this data in a second definition of “regions”, and sometimes I need to do the opposite. I have th… --- ## [Converting UTM to EPSG:5389](https://discourse.julialang.org/t/converting-utm-to-epsg-5389/87280)
**Author:** [@ameresv](https://discourse.julialang.org/u/ameresv)\ **Replies:** 5\ **Last updated:** [September 14, 2022, 10:20pm UTC](https://discourse.julialang.org/t/converting-utm-to-epsg-5389/87280 "2022-09-14T22:20:54Z")
It could sound like a silly question, but I am having trouble converting points from UTM to long - lat (researching a bit, I found that I need is EPSG:5389 because is on the Southern side of Peru). So, the data looks li… --- ## [NCDatasets: copy the properties of a variable](https://discourse.julialang.org/t/ncdatasets-copy-the-properties-of-a-variable/85721)
**Author:** [@ryofurue](https://discourse.julialang.org/u/ryofurue)\ **Replies:** 11\ **Last updated:** [September 3, 2022, 2:13pm UTC](https://discourse.julialang.org/t/ncdatasets-copy-the-properties-of-a-variable/85721 "2022-09-03T14:13:13Z")
Is there a facility to copy a netCDF variable including its axes (dimensions) and attributes at once? A scenario: ds\_org = Dataset("original.nc", "r") org = ds\["orgvar"\] Dataset("new.nc", "c") do ds\_new defVar(ds\_new… --- ## [DimensionMismatch for a Raster](https://discourse.julialang.org/t/dimensionmismatch-for-a-raster/86677)
**Author:** [@ryofurue](https://discourse.julialang.org/u/ryofurue)\ **Replies:** 1\ **Last updated:** [September 2, 2022, 9:33am UTC](https://discourse.julialang.org/t/dimensionmismatch-for-a-raster/86677 "2022-09-02T09:33:45Z")
A Raster behaves just like an Array in suitable contexts, but I found the following difference in broadcast behavior. Is it 1) an intended behavior, 2) an unavoidable inconsistency, or 3) a bug? Basically, avr = sum(v… --- ## [TimeZone from latitude/longtitude](https://discourse.julialang.org/t/timezone-from-latitude-longtitude/55129)
**Author:** [@lubo](https://discourse.julialang.org/u/lubo)\ **Replies:** 7\ **Last updated:** [August 22, 2022, 3:50pm UTC](https://discourse.julialang.org/t/timezone-from-latitude-longtitude/55129 "2022-08-22T15:50:44Z")
Hi, maybe very silly question, but can not find solution so far. I have gps coordinates of the city, ie place in Stockholm/Sweden: \[59.3725,18.0\]. I would like to get TimeZone for this place: ‘Europe/Stockholm’ In py… --- ## [Visualising NetCDF using Raster.jl](https://discourse.julialang.org/t/visualising-netcdf-using-raster-jl/85761)
**Author:** [@Michaela\_Vollmer](https://discourse.julialang.org/u/Michaela_Vollmer)\ **Replies:** 10\ **Last updated:** [August 16, 2022, 9:29pm UTC](https://discourse.julialang.org/t/visualising-netcdf-using-raster-jl/85761 "2022-08-16T21:29:25Z")
Hello! I am trying to plot spatial data of a NetCDF file from a WFlow example. I have followed the steps outlined in the following question Raster.jl discourse post where a similar example is shown where the data can be … --- ## [Error while reading and writing raster dataset using ArchGDAL](https://discourse.julialang.org/t/error-while-reading-and-writing-raster-dataset-using-archgdal/85268)
**Author:** [@anu21](https://discourse.julialang.org/u/anu21)\ **Replies:** 2\ **Last updated:** [August 5, 2022, 6:20pm UTC](https://discourse.julialang.org/t/error-while-reading-and-writing-raster-dataset-using-archgdal/85268 "2022-08-05T18:20:36Z")
I am new to the Julia. I am using ArchGDAL to open and write raster data, but I keep getting this error. GDALError (CE\_Failure, code 1): GDALOpen() called on “/path/to/file” recursively Is anyone familiar with this err… --- ## [Kriging on surface of sphere](https://discourse.julialang.org/t/kriging-on-surface-of-sphere/85058)
**Author:** [@yakir12](https://discourse.julialang.org/u/yakir12)\ **Replies:** 5\ **Last updated:** [August 1, 2022, 1:47pm UTC](https://discourse.julialang.org/t/kriging-on-surface-of-sphere/85058 "2022-08-01T13:47:39Z")
I have a few random points on a sphere’s surface and I’d like to interpolate and extrapolate those using kriging. But I’m getting artefacts in the result that look like an issue with the fact that the coordinates are sph… --- ## [Splitting polygons into narrow and wide parts](https://discourse.julialang.org/t/splitting-polygons-into-narrow-and-wide-parts/84307)
**Author:** [@Wikunia](https://discourse.julialang.org/u/Wikunia)\ **Replies:** 3\ **Last updated:** [July 16, 2022, 10:59pm UTC](https://discourse.julialang.org/t/splitting-polygons-into-narrow-and-wide-parts/84307 "2022-07-16T22:59:59Z")
Hey everyone, I’m interested in a general polygon splitting algorithm that is able to identify very narrow regions and can split up narrow from wide regions. I found this via google: and was wondering whether this is… --- ## [Display ArchGDAL geotiff raster wants unsigned type?](https://discourse.julialang.org/t/display-archgdal-geotiff-raster-wants-unsigned-type/82813)
**Author:** [@floswald](https://discourse.julialang.org/u/floswald)\ **Replies:** 2\ **Last updated:** [June 15, 2022, 8:02pm UTC](https://discourse.julialang.org/t/display-archgdal-geotiff-raster-wants-unsigned-type/82813 "2022-06-15T20:02:06Z")
I wanted to look at this raster with ArchGDAL: using ArchGDAL (loaded raster from disk) julia\> e GDAL Dataset (Driver: GTiff/GeoTIFF) File(s): /Users/florian.oswald/Downloads/population\_fra\_2019-07-01\_geotiff/popula… --- ## [NCDatasets.jl: define values of dimensions](https://discourse.julialang.org/t/ncdatasets-jl-define-values-of-dimensions/65892)
**Author:** [@Inna](https://discourse.julialang.org/u/Inna)\ **Replies:** 4\ **Last updated:** [June 15, 2022, 9:48am UTC](https://discourse.julialang.org/t/ncdatasets-jl-define-values-of-dimensions/65892 "2022-06-15T09:48:03Z")
Dear all, Using NCDatasets, To write a data to netcdf, i can define the dimensions with this way: ds = NCDataset("/tmp/test.nc","c") defDim(ds,"lon",100) However, is there any way to define the values of my dimension… --- ## [ArchGDAL error after recent update](https://discourse.julialang.org/t/archgdal-error-after-recent-update/80987)
**Author:** [@scotts](https://discourse.julialang.org/u/scotts)\ **Replies:** 2\ **Last updated:** [May 13, 2022, 8:03am UTC](https://discourse.julialang.org/t/archgdal-error-after-recent-update/80987 "2022-05-13T08:03:21Z")
I had some very basic geocode working for plotting a tif file that has recently stopped working. It’s very basic, from here: https://www.acgeospatial.co.uk/julia-prt3/ using ArchGDAL, Plots, Colors, Images, GeoData con… --- ## [Mapping polygon or mesh data from array input for satellite swaths](https://discourse.julialang.org/t/mapping-polygon-or-mesh-data-from-array-input-for-satellite-swaths/80568)
**Author:** [@haffner](https://discourse.julialang.org/u/haffner)\ **Replies:** 10\ **Last updated:** [May 7, 2022, 6:21pm UTC](https://discourse.julialang.org/t/mapping-polygon-or-mesh-data-from-array-input-for-satellite-swaths/80568 "2022-05-07T18:21:44Z")
I am new to Julia and work on satellite remote sensing of the atmosphere and ocean. I am looking for some advice on mapping point and polygon (vector) data in geographic coordinates using Julia. I produce relatively co… --- ## [How can I get Proj4 to transform multiple locations?](https://discourse.julialang.org/t/how-can-i-get-proj4-to-transform-multiple-locations/80098)
**Author:** [@chadagreene](https://discourse.julialang.org/u/chadagreene)\ **Replies:** 11\ **Last updated:** [April 29, 2022, 12:32am UTC](https://discourse.julialang.org/t/how-can-i-get-proj4-to-transform-multiple-locations/80098 "2022-04-29T00:32:20Z")
This is probably more of a Julia question than a Proj4 question, but how can I perform coordinate transformations on multiple locations at once? I’m writing a function that would need to handle the following three cases: … --- ## [WebMercatorfromLLA not defined](https://discourse.julialang.org/t/webmercatorfromlla-not-defined/78358)
**Author:** [@julio.ca2020](https://discourse.julialang.org/u/julio.ca2020)\ **Replies:** 9\ **Last updated:** [April 26, 2022, 7:59pm UTC](https://discourse.julialang.org/t/webmercatorfromlla-not-defined/78358 "2022-04-26T19:59:33Z")
I am trying to use WebMercatorfromLLA… Julia tells me that is not defined. Any other transformation, function work ERROR: UndefVarError: WebMercatorfromLLA not defined Stacktrace: \[1\] top-level scope @ ~/Documents/Hardw… --- ## [Reading ENVI style multispectral images with archGDAL](https://discourse.julialang.org/t/reading-envi-style-multispectral-images-with-archgdal/79702)
**Author:** [@Arrigo\_Benedetti](https://discourse.julialang.org/u/Arrigo_Benedetti)\ **Replies:** 4\ **Last updated:** [April 20, 2022, 2:09pm UTC](https://discourse.julialang.org/t/reading-envi-style-multispectral-images-with-archgdal/79702 "2022-04-20T14:09:13Z")
I have several multispectral images with separate ENVI style header text files and I would like to load them with archGDAL.jl. I looked at the documentation of archGDAL.read() but I could not find how to do it. Maybe arc… --- ## [Calculate area of a raster](https://discourse.julialang.org/t/calculate-area-of-a-raster/79223)
**Author:** [@ahmoreira](https://discourse.julialang.org/u/ahmoreira)\ **Replies:** 4\ **Last updated:** [April 20, 2022, 8:30am UTC](https://discourse.julialang.org/t/calculate-area-of-a-raster/79223 "2022-04-20T08:30:43Z")
Hi all, I am trying to calculate the area of urban land cover only for cells in which 100% of the cell is urban cover. I am calculating the area as the product of the resolution times the number of cells, but it is a bi… --- ## [What is the state of Rasterization and Rasterstats in Julia](https://discourse.julialang.org/t/what-is-the-state-of-rasterization-and-rasterstats-in-julia/72366)
**Author:** [@Ordens\_Ritter](https://discourse.julialang.org/u/Ordens_Ritter)\ **Replies:** 26\ **Last updated:** [April 8, 2022, 11:48am UTC](https://discourse.julialang.org/t/what-is-the-state-of-rasterization-and-rasterstats-in-julia/72366 "2022-04-08T11:48:56Z")
I have followed the progress of the geo community in Julia every know and than. To me there seams to be a slow but steady progress. That both for raster and vector files. Although, I for instead didn’t found out thinks,… --- ## [Speed up ArchGDAL.contains for thousands of data](https://discourse.julialang.org/t/speed-up-archgdal-contains-for-thousands-of-data/77847)
**Author:** [@Frankiewaang](https://discourse.julialang.org/u/Frankiewaang)\ **Replies:** 8\ **Last updated:** [March 25, 2022, 12:39pm UTC](https://discourse.julialang.org/t/speed-up-archgdal-contains-for-thousands-of-data/77847 "2022-03-25T12:39:52Z")
Hi Guys, I have geopoints like this import ArchGDAL as AG lon = LinRange(-180,180,7500) lat = LinRange(-90,90,5500) function meshgrid(x, y) X = \[i for i in x, j in 1:length(y)\] Y = \[j for i in 1:length(x), j i… --- ## [WMS - Web Map Services displays](https://discourse.julialang.org/t/wms-web-map-services-displays/14483)
**Author:** [@johnh](https://discourse.julialang.org/u/johnh)\ **Replies:** 12\ **Last updated:** [February 23, 2022, 8:30pm UTC](https://discourse.julialang.org/t/wms-web-map-services-displays/14483 "2022-02-23T20:30:45Z")
May I ask if anyone is familiar with WMS displays of data Web Map Service - MATLAB & Simulink - MathWorks United Kingdom Web Map Service - Wikipedia I ask for a good reason, as I would like to do some example plots usi… --- ## [Calculating a line of bearing between two LLAs](https://discourse.julialang.org/t/calculating-a-line-of-bearing-between-two-llas/75244)
**Author:** [@devin](https://discourse.julialang.org/u/devin)\ **Replies:** 4\ **Last updated:** [January 27, 2022, 7:56pm UTC](https://discourse.julialang.org/t/calculating-a-line-of-bearing-between-two-llas/75244 "2022-01-27T19:56:33Z")
Hi, I’ve got two lat/lon/alt locations and I’d like to compute the line of bearing between them (az/el). Is there a way to do this using geodesy.jl? I’ve looked around and found GreatCircle.jl, geodesics.jl and LibGeo… --- ## [So now I convert my OSGB36 data to LLA ... oh ... that looks ... hard](https://discourse.julialang.org/t/so-now-i-convert-my-osgb36-data-to-lla-oh-that-looks-hard/74400)
**Author:** [@lawless-m](https://discourse.julialang.org/u/lawless-m)\ **Replies:** 10\ **Last updated:** [January 13, 2022, 2:15pm UTC](https://discourse.julialang.org/t/so-now-i-convert-my-osgb36-data-to-lla-oh-that-looks-hard/74400 "2022-01-13T14:15:23Z")
Hi, The UK’s Office of National Statistics publishes data which uses OSGB36 East / North coordinates. Of which I have 39 million addresses from the ONSUD database e.g. East/North “426642,380231” Using an online con… --- ## [GeoStats.jl workshop at FOSS4G 2021](https://discourse.julialang.org/t/geostats-jl-workshop-at-foss4g-2021/65889)
**Author:** [@juliohm](https://discourse.julialang.org/u/juliohm)\ **Replies:** 0\ **Last updated:** [August 5, 2021, 1:42pm UTC](https://discourse.julialang.org/t/geostats-jl-workshop-at-foss4g-2021/65889 "2021-08-05T13:42:05Z")
FOSS4G is the largest international conference on geospatial open source software. I invite everyone to participate in the GeoStats.jl workshop therein. We will cover the basics of the Julia programming language and the… --- ## [Delaunay Triangularization](https://discourse.julialang.org/t/delaunay-triangularization/36172)
**Author:** [@Donald\_Lacombe](https://discourse.julialang.org/u/Donald_Lacombe)\ **Replies:** 20\ **Last updated:** [January 4, 2022, 1:55am UTC](https://discourse.julialang.org/t/delaunay-triangularization/36172 "2022-01-04T01:55:26Z")
Hello! I’d like to be able to obtain the following output from MATLAB but in Julia. Here is a toy example in MATLAB: \>\> x = rand(5,1) x = 0.8147 0.9058 0.1270 0.9134 0.6324 \>\> y = rand(5,1) y = … --- ## [Read multi .nc files along new dimension](https://discourse.julialang.org/t/read-multi-nc-files-along-new-dimension/70696)
**Author:** [@Inna](https://discourse.julialang.org/u/Inna)\ **Replies:** 10\ **Last updated:** [December 22, 2021, 7:52am UTC](https://discourse.julialang.org/t/read-multi-nc-files-along-new-dimension/70696 "2021-12-22T07:52:37Z")
Dear all, I have some netcdf files containg temperature and precipitation, with size 360 x 180. There isn’t inside the .nc file the dimension of time, but it is on the name of each file (2012.01.nc, 2012.02.nc, … ). I w… --- ## [Error in resampling with Rasters.jl](https://discourse.julialang.org/t/error-in-resampling-with-rasters-jl/73261)
**Author:** [@MJulia](https://discourse.julialang.org/u/MJulia)\ **Replies:** 6\ **Last updated:** [December 22, 2021, 2:01am UTC](https://discourse.julialang.org/t/error-in-resampling-with-rasters-jl/73261 "2021-12-22T02:01:02Z")
I’m trying to resample a data of lambert\_azimuthal\_equal\_area projection to normal Lat/Lon coordinates by using Rasters.jl as following: using NCDatasets using ArchGDAL using DimendionalData using Rasters ds = Data… [Previous page](https://discourse.julialang.org/c/domain/geo/18.md?page=3) [Next page](https://discourse.julialang.org/c/domain/geo/18.md?page=5)