Photometry
Installation
To install from the REPL, enter Pkg-mode (]
)
pkg> add Photometry
Getting Started
Here is a basic example to do some aperture photometry using CircularAperture
and CircularAnnulus
. The aperture_photometry
function performs the photometry using a given method.
data = ones(100, 100)
err = ones(100, 100)
ap1 = CircularAperture(50, 50, 3)
# partial overlap
ap2 = CircularAperture(0.5, 0.5, 5)
results = aperture_photometry([ap1, ap2], data, err)
@assert results.aperture_sum[1] ≈ 9π
@assert results.aperture_sum[2] ≈ 25π / 4
results
# output
2×4 DataFrames.DataFrame
│ Row │ xcenter │ ycenter │ aperture_sum │ aperture_sum_err │
│ │ Any │ Any │ Any │ Any │
├─────┼─────────┼─────────┼──────────────┼──────────────────┤
│ 1 │ 50 │ 50 │ 28.2743 │ 5.31736 │
│ 2 │ 0.5 │ 0.5 │ 19.635 │ 4.43113 │
Contributing
If you're interested in contributing, go ahead and check out the issues or make a pull request. If you add a new feature, please write appropriate unit tests for it and bump the package's minor version.