Tools for interstellar dust extinction in astronomy
Documentation: DustExtinction
From package manager (] key)
pkg> add DustExtinctionColor laws:
julia> using DustExtinction
julia> CCM89(Rv=3.1)(4000)
1.4645557029425842
julia> CCM89(Rv=3.1).([4000, 5000])
2-element Vector{Float64}:
1.46456
1.12225Dust maps:
julia> dustmap = SFD98Map()
SFD98Map("[...]")
julia> dustmap(0.1, 0.1)
0.793093095733043
julia> dustmap.([0.1, 0.2], [0.1, 0.2])
2-element Vector{Float64}:
0.793093
0.539507Reddening:
julia> wave = [4000., 5000.]
2-element Vector{Float64}:
4000.0
5000.0
julia> flux = [1.0, 1.5]
2-element Vector{Float64}:
1.0
1.5
julia> red = redden.(CCM89, wave, flux; Av=0.3, Rv=3.1)
2-element Vector{Float64}:
0.6671958182723856
1.1000733242882896
julia> deredden.(CCM89(Rv=3.1), wave, red; Av=0.3)
2-element Vector{Float64}:
1.0
1.5We provide first-class support for Unitful.jl and Measurements.jl packages, too! Check out the documentation for more examples.
Feel free to open an issue or a pull-request for any discussion, suggestions, new features, or patches!