evetion / geoarrays.jl Goto Github PK

ERROR: UndefVarError: coords not defined
Stacktrace:
 [1] macro expansion
   @ ~/.julia/packages/GeoArrays/oqARi/src/plot.jl:10 [inlined]
 [2] apply_recipe(plotattributes::AbstractDict{Symbol, Any}, ga::GeoArray)
   @ GeoArrays ~/.julia/packages/RecipesBase/qpxEX/src/RecipesBase.jl:289
 [3] _process_userrecipes!(plt::Any, plotattributes::Any, args::Any)
   @ RecipesPipeline ~/.julia/packages/RecipesPipeline/7ijBv/src/user_recipe.jl:36
 [4] recipe_pipeline!(plt::Any, plotattributes::Any, args::Any)
   @ RecipesPipeline ~/.julia/packages/RecipesPipeline/7ijBv/src/RecipesPipeline.jl:70
 [5] _plot!(plt::Plots.Plot, plotattributes::Any, args::Any)
   @ Plots ~/.julia/packages/Plots/9C6z9/src/plot.jl:208
 [6] plot(args::Any; kw::Base.Pairs{Symbol, V, Tuple{Vararg{Symbol, N}}, NamedTuple{names, T}} where {V, N, names, T<:Tuple{Vararg{Any, N}}})
   @ Plots ~/.julia/packages/Plots/9C6z9/src/plot.jl:91
 [7] top-level scope

julia> coords
WARNING: both Plots and GeoArrays export "coords"; uses of it in module Main must be qualified
ERROR: UndefVarError: coords not defined

fn = download("https://bitbucket.org/earthobservation/julia-remote-sensing/raw/ccea113bf43ba71d4d1c6f02e060fd2923c982b2/data/KR_S2B_MSIL2A_20210614_20.tif")
geoarray = GeoArrays.read(fn)

ga_band = geoarray[:,:,7] # Matrix{UInt16}

function ga_subset(ga, bands, ind)
	afm = ga.f
	r_min = ind[1]
	r_max = ind[3]
	c_min = ind[2]
	c_max = ind[4]
	sub_afm_l = afm.linear
	sub_a = ga.A[r_min:r_max,c_min:c_max,bands]
	sub_afm_t = GeoArrays.coords(ga, [r_min,c_min])
	sub_f = CoordinateTransformations.AffineMap(sub_afm_l, sub_afm_t)
	sub_crs = ga.crs
	sub_ga = GeoArray(sub_a, sub_f, sub_crs)
    return sub_ga
end

ga_sub = ga_subset(geoarray, [1,7], [200,200,400,300]);
typeof(ga_sub) # 201x101x2 Array{UInt16, 3} with AffineMap([20.0 0.0; 0.0 -20.0], [448980.0, 5.12002e6]) ...

julia> using GeoArrays
[ Info: Precompiling GeoArrays [2fb1d81b-e6a0-5fc5-82e6-8e06903437ab]
WARNING: could not import TableTransforms.ColSpec into GeoStatsBase
WARNING: could not import TableTransforms.Col into GeoStatsBase
WARNING: could not import TableTransforms.colspec into GeoStatsBase
WARNING: could not import TableTransforms.choose into GeoStatsBase
ERROR: LoadError: UndefVarError: `ColSpec` not defined
Stacktrace:
  [1] top-level scope
    @ ~/.julia/packages/GeoStatsBase/hN0RT/src/transforms/detrend.jl:5
  [2] include(mod::Module, _path::String)
    @ Base ./Base.jl:457
  [3] include(x::String)
    @ GeoStatsBase ~/.julia/packages/GeoStatsBase/hN0RT/src/GeoStatsBase.jl:5
  [4] top-level scope
    @ ~/.julia/packages/GeoStatsBase/hN0RT/src/transforms.jl:14
  [5] include(mod::Module, _path::String)
    @ Base ./Base.jl:457
  [6] include(x::String)
    @ GeoStatsBase ~/.julia/packages/GeoStatsBase/hN0RT/src/GeoStatsBase.jl:5
  [7] top-level scope
    @ ~/.julia/packages/GeoStatsBase/hN0RT/src/GeoStatsBase.jl:56
  [8] include
    @ ./Base.jl:457 [inlined]
  [9] include_package_for_output(pkg::Base.PkgId, input::String, depot_path::Vector{String}, dl_load_path::Vector{String}, load_path::Vector{String}, concrete_deps::Vector{Pair{Base.PkgId, UInt128}}, source::String)
    @ Base ./loading.jl:2049
 [10] top-level scope
    @ stdin:3
in expression starting at ~/.julia/packages/GeoStatsBase/hN0RT/src/transforms/detrend.jl:5
in expression starting at ~/.julia/packages/GeoStatsBase/hN0RT/src/transforms.jl:14
in expression starting at ~/.julia/packages/GeoStatsBase/hN0RT/src/GeoStatsBase.jl:5
in expression starting at stdin:3
ERROR: LoadError: Failed to precompile GeoStatsBase [323cb8eb-fbf6-51c0-afd0-f8fba70507b2] to "~/.julia/compiled/v1.9/GeoStatsBase/jl_OK8pRv".
Stacktrace:
  [1] error(s::String)
    @ Base ./error.jl:35
  [2] compilecache(pkg::Base.PkgId, path::String, internal_stderr::IO, internal_stdout::IO, keep_loaded_modules::Bool)
    @ Base ./loading.jl:2300
  [3] compilecache
    @ ./loading.jl:2167 [inlined]
  [4] _require(pkg::Base.PkgId, env::String)
    @ Base ./loading.jl:1805
  [5] _require_prelocked(uuidkey::Base.PkgId, env::String)
    @ Base ./loading.jl:1660
  [6] macro expansion
    @ ./loading.jl:1648 [inlined]
  [7] macro expansion
    @ ./lock.jl:267 [inlined]
  [8] require(into::Module, mod::Symbol)
    @ Base ./loading.jl:1611
  [9] include(mod::Module, _path::String)
    @ Base ./Base.jl:457
 [10] include(x::String)
    @ GeoArrays ~/.julia/packages/GeoArrays/Is81p/src/GeoArrays.jl:1
 [11] top-level scope
    @ ~/.julia/packages/GeoArrays/Is81p/src/GeoArrays.jl:17
 [12] include
    @ ./Base.jl:457 [inlined]
 [13] include_package_for_output(pkg::Base.PkgId, input::String, depot_path::Vector{String}, dl_load_path::Vector{String}, load_path::Vector{String}, concrete_deps::Vector{Pair{Base.PkgId, UInt128}}, source::Nothing)
    @ Base ./loading.jl:2049
 [14] top-level scope
    @ stdin:3
in expression starting at ~/.julia/packages/GeoArrays/Is81p/src/interpolate.jl:1
in expression starting at ~/.julia/packages/GeoArrays/Is81p/src/GeoArrays.jl:1
in expression starting at stdin:3
ERROR: Failed to precompile GeoArrays [2fb1d81b-e6a0-5fc5-82e6-8e06903437ab] to "~/.julia/compiled/v1.9/GeoArrays/jl_eJEmzF".
Stacktrace:
 [1] error(s::String)
   @ Base ./error.jl:35
 [2] compilecache(pkg::Base.PkgId, path::String, internal_stderr::IO, internal_stdout::IO, keep_loaded_modules::Bool)
   @ Base ./loading.jl:2300
 [3] compilecache
   @ ./loading.jl:2167 [inlined]
 [4] _require(pkg::Base.PkgId, env::String)
   @ Base ./loading.jl:1805
 [5] _require_prelocked(uuidkey::Base.PkgId, env::String)
   @ Base ./loading.jl:1660
 [6] macro expansion
   @ ./loading.jl:1648 [inlined]
 [7] macro expansion
   @ ./lock.jl:267 [inlined]
 [8] require(into::Module, mod::Symbol)
   @ Base ./loading.jl:1611

julia> ga = GeoArray(rand(Int, (100,100)))
100x100x1 Array{Int64,3} with AffineMap([1.0 0.0; 0.0 1.0], [0.0, 0.0]) and undefined CRS

julia> GeoArrays.write!("test.tif", ga)
ERROR: KeyError: key Int64 not found
Stacktrace:
 [1] getindex at ./dict.jl:477 [inlined]
 [2] #unsafe_create#30 at /Users/evetion/.julia/packages/ArchGDAL/FiLcr/src/dataset.jl:200 [inlined]
 [3] create(::GeoArrays.var"#27#28"{GeoArray{Int64},Int64}, ::String; kwargs::Base.Iterators.Pairs{Symbol,Any,NTuple{6,Symbol},NamedTuple{(:driver, :width, :height, :nbands, :dtype, :options),Tuple{ArchGDAL.Driver,Int64,Int64,Int64,DataType,Array{String,1}}}}) at /Users/evetion/.julia/packages/ArchGDAL/FiLcr/src/context.jl:213
 [4] write!(::String, ::GeoArray{Int64}, ::Nothing, ::String) at /Users/evetion/.julia/packages/GeoArrays/95hSL/src/io.jl:80
 [5] write!(::String, ::GeoArray{Int64}) at /Users/evetion/.julia/packages/GeoArrays/95hSL/src/io.jl:61
 [6] top-level scope at REPL[35]:1

A::AbstractArray{T,3}
f::AffineMap
crs::WellKnownText{GeoFormatTypes.CRS,<:String}

(v1.3) pkg> add GeoArrays
  Updating registry at `~/.julia/registries/General`
  Updating git-repo `https://github.com/JuliaRegistries/General.git`
 Resolving package versions...
ERROR: Unsatisfiable requirements detected for package GeoStatsBase [323cb8eb]:
 GeoStatsBase [323cb8eb] log:
 ├─possible versions are: [0.2.0-0.2.6, 0.3.0-0.3.6, 0.4.0-0.4.3, 0.5.0, 0.6.0-0.6.4, 0.7.0-0.7.1, 0.8.0-0.8.7, 0.9.0-0.9.5, 0.10.0-0.10.1] or uninstalled
 ├─restricted by compatibility requirements with StaticArrays [90137ffa] to versions: [0.2.0-0.2.6, 0.3.0-0.3.6, 0.4.0-0.4.3, 0.6.4, 0.7.0-0.7.1, 0.8.0-0.8.7, 0.9.0-0.9.5, 0.10.0-0.10.1] or uninstalled
 │ └─StaticArrays [90137ffa] log:
 │   ├─possible versions are: [0.8.0-0.8.3, 0.9.0-0.9.2, 0.10.0, 0.10.2-0.10.3, 0.11.0-0.11.1, 0.12.0-0.12.3] or uninstalled
 │   ├─restricted to versions 0.12 by GeometryBasics [5c1252a2], leaving only versions 0.12.0-0.12.3
 │   │ └─GeometryBasics [5c1252a2] log:
 │   │   ├─possible versions are: 0.2.13 or uninstalled
 │   │   ├─restricted to versions * by GeoJSONTables [b51d8a2e], leaving only versions 0.2.13
 │   │   │ └─GeoJSONTables [b51d8a2e] log:
 │   │   │   ├─possible versions are: 0.1.0 or uninstalled
 │   │   │   └─GeoJSONTables [b51d8a2e] is fixed to version 0.1.0
 │   │   └─GeometryBasics [5c1252a2] is fixed to version 0.2.13
 │   └─restricted to versions 0.12.3 by an explicit requirement, leaving only versions 0.12.3
 ├─restricted by compatibility requirements with DataFrames [a93c6f00] to versions: [0.2.0-0.2.6, 0.3.0-0.3.6, 0.4.0-0.4.3, 0.10.0-0.10.1] or uninstalled
 │ └─DataFrames [a93c6f00] log:
 │   ├─possible versions are: [0.11.7, 0.12.0, 0.13.0-0.13.1, 0.14.0-0.14.1, 0.15.0-0.15.2, 0.16.0, 0.17.0-0.17.1, 0.18.0-0.18.4, 0.19.0-0.19.4, 0.20.0-0.20.2, 0.21.0-0.21.2] or uninstalled
 │   ├─restricted to versions * by GeoDataFrames [fbc9b913], leaving only versions [0.11.7, 0.12.0, 0.13.0-0.13.1, 0.14.0-0.14.1, 0.15.0-0.15.2, 0.16.0, 0.17.0-0.17.1, 0.18.0-0.18.4, 0.19.0-0.19.4, 0.20.0-0.20.2, 0.21.0-0.21.2]
 │   │ └─GeoDataFrames [fbc9b913] log:
 │   │   ├─possible versions are: 0.0.0 or uninstalled
 │   │   └─GeoDataFrames [fbc9b913] is fixed to version 0.0.0
 │   └─restricted to versions 0.21.0 by an explicit requirement, leaving only versions 0.21.0
 └─restricted by compatibility requirements with GeoArrays [2fb1d81b] to versions: [0.5.0, 0.6.0-0.6.4, 0.7.0-0.7.1, 0.8.3-0.8.7] — no versions left
   └─GeoArrays [2fb1d81b] log:
     ├─possible versions are: [0.1.0-0.1.7, 0.2.0-0.2.2, 0.3.0-0.3.1] or uninstalled
     ├─restricted to versions * by an explicit requirement, leaving only versions [0.1.0-0.1.7, 0.2.0-0.2.2, 0.3.0-0.3.1]
     ├─restricted by compatibility requirements with RecipesBase [3cdcf5f2] to versions: [0.1.0-0.1.4, 0.2.2, 0.3.0-0.3.1] or uninstalled, leaving only versions: [0.1.0-0.1.4, 0.2.2, 0.3.0-0.3.1]
     │ └─RecipesBase [3cdcf5f2] log:
     │   ├─possible versions are: [0.4.0, 0.5.0, 0.6.0, 0.7.0, 0.8.0, 1.0.0-1.0.1] or uninstalled
     │   ├─restricted to versions * by GeoDataFrames [fbc9b913], leaving only versions [0.4.0, 0.5.0, 0.6.0, 0.7.0, 0.8.0, 1.0.0-1.0.1]
     │   │ └─GeoDataFrames [fbc9b913] log: see above
     │   ├─restricted by compatibility requirements with GeoInterface [cf35fbd7] to versions: [0.6.0, 0.7.0, 0.8.0, 1.0.0-1.0.1]
     │   │ └─GeoInterface [cf35fbd7] log:
     │   │   ├─possible versions are: [0.4.0-0.4.1, 0.5.0-0.5.4] or uninstalled
     │   │   ├─restricted to versions 0.4-0.5 by LibGEOS [a90b1aa1], leaving only versions [0.4.0-0.4.1, 0.5.0-0.5.4]
     │   │   │ └─LibGEOS [a90b1aa1] log:
     │   │   │   ├─possible versions are: 0.6.4 or uninstalled
     │   │   │   └─LibGEOS [a90b1aa1] is fixed to version 0.6.4
     │   │   ├─restricted to versions 0.5 by GeoMakie [db073c08], leaving only versions 0.5.0-0.5.4
     │   │   │ └─GeoMakie [db073c08] log:
     │   │   │   ├─possible versions are: 0.1.14 or uninstalled
     │   │   │   └─GeoMakie [db073c08] is fixed to version 0.1.14-DEV
     │   │   └─restricted to versions 0.5.4 by an explicit requirement, leaving only versions 0.5.4
     │   └─restricted by compatibility requirements with Plots [91a5bcdd] to versions: 1.0.0-1.0.1
     │     └─Plots [91a5bcdd] log:
     │       ├─possible versions are: [0.12.1-0.12.4, 0.13.0-0.13.1, 0.14.0-0.14.2, 0.15.0-0.15.1, 0.16.0, 0.17.0-0.17.4, 0.18.0, 0.19.0-0.19.3, 0.20.0-0.20.6, 0.21.0, 0.22.0-0.22.5, 0.23.0-0.23.2, 0.24.0, 0.25.0-0.25.3, 0.26.0-0.26.3, 0.27.0-0.27.1, 0.28.0-0.28.4, 0.29.0-0.29.9, 1.0.0-1.0.14, 1.1.0-1.1.4, 1.2.0-1.2.6, 1.3.0-1.3.7, 1.4.0-1.4.3] or uninstalled
     │       └─restricted to versions 1.0.14 by an explicit requirement, leaving only versions 1.0.14
     └─restricted by compatibility requirements with ArchGDAL [c9ce4bd3] to versions: [0.1.4-0.1.7, 0.2.0-0.2.2, 0.3.0] or uninstalled, leaving only versions: [0.1.4, 0.2.2, 0.3.0]
       └─ArchGDAL [c9ce4bd3] log:
         ├─possible versions are: [0.1.0, 0.2.0-0.2.2, 0.3.0-0.3.3, 0.4.0-0.4.1] or uninstalled
         └─restricted to versions 0.3.2 by an explicit requirement, leaving only versions 0.3.2