Skip to content

Commit

Permalink
Cleanup, add docs.
Browse files Browse the repository at this point in the history
  • Loading branch information
DinoBektesevic committed Apr 12, 2024
1 parent f25f0de commit d198e5c
Show file tree
Hide file tree
Showing 3 changed files with 430 additions and 102 deletions.
23 changes: 18 additions & 5 deletions src/kbmod/mocking/fits.py
Original file line number Diff line number Diff line change
Expand Up @@ -4,12 +4,25 @@

import numpy as np
from astropy.wcs import WCS
from astropy.io.fits import HDUList, PrimaryHDU, CompImageHDU, BinTableHDU, Header
from astropy.modeling import models
from astropy.io.fits import (
HDUList,
PrimaryHDU,
CompImageHDU,
BinTableHDU,
Header
)

from .headers import HeaderFactory, ArchivedHeader
from .catalogs import gen_catalog, SimpleSourceCatalog, SimpleObjectCatalog
from .fits_data import ZeroedData, SimpleImage, SimpleVariance, SimpleMask, add_model_objects
from .fits_data import (
DataFactory,
ZeroedData,
SimpleImage,
SimpleVariance,
SimpleMask,
add_model_objects
)


__all__ = [
Expand Down Expand Up @@ -175,10 +188,10 @@ def __init__(self, shape, noise=0, noise_std=1, dt=0.001, source_catalog=None, o
shape=shape,
noise=noise,
noise_std=noise_std,
catalog=self.src_cat
src_cat=self.src_cat,
)
self.var_data = SimpleVariance(self.img_data.base, read_noise=noise, gain=1.0)
self.mask_data = ZeroedData(np.zeros(shape))
self.mask_data = DataFactory(np.zeros(shape))

# Now we can build the HDU map and the HDUList layout
layout = [
Expand Down Expand Up @@ -207,7 +220,7 @@ def mock(self, n=1):
if self.obj_cat is not None:
obj_cats = self.obj_cat.mock(n, dt=self.dt)

images = self.img_data.mock(n, obj_cats)
images = self.img_data.mock(n, obj_cats=obj_cats)
variances = self.var_data.mock(images=images)
mask = self.mask_data.mock()
imghdr, varhdr, maskhdr = self.image.mock(), self.variance.mock(), self.mask.mock()
Expand Down
Loading

0 comments on commit d198e5c

Please sign in to comment.