] update
in your environment or ] update GeoStatsBase
to only update the package.
sol = solve(prob, solver);
Hello everyone, I have a question about exporting/writing spatial data. After doing simulations with GeoStats.jl, I got data with the data structure "GeoData{CartesianGrid{2,Float64},TypedTables.Table{NamedTuple{(:X1,),Tuple{Float64}},1,NamedTuple{(:X1,),Tuple{Array{Float64,1}}}}}". How to export it as raster data?
I tried NCDatasets.jl, but seems like they don't provide a method to convert GeoData to netCDF? Do you happen to know the best package to export the simulation results of GeoStats.jl? Thanks!
r = ensemble[1]
will give you the first realization of the process. You can then retrieve the domain of this realization with d = domain(r)
. Let's assume that d
is a CartesianGrid{2}
. This domain has finite elements, in this case quadrangles. For example, d[1]
is the first quadrangle. You can get the centroid of all quadrangles with centroid.(d)
and the result will be a vector of Point
. Finally, you can get the coordinates of these points in a global reference system using coordinates.(centroid.(d))
.