Prepare data for time domainΒΆ
In order to mimic the realistic situation where only a limited number of observed epochs is available at each
day, it is necessary to prepare our simulate data resemble this scenario. In actsnclass
this is done in
5 steps:
- Determine minimum and maximum MJD for the entire SNPCC sample;
- For each day of the survey, run through the entire data sample and select only the observed epochs which were obtained prior to it;
- Perform the feature extraction process considering only the photometric points which survived item 2.
- Check if in the MJD in question the object is available for querying.
- Join all information in a standard features file.
You can perform the entire analysis for one day of the survey using the actsnclass.time_domain
module:
1 2 3 4 5 6 7 8 9 10 | >>> from actsnclass.time_domain import SNPCCPhotometry >>> path_to_data = 'data/SIMGEN_PUBLIC_DES/' >>> output_dir = 'results/time_domain/' >>> day = 20 >>> data = SNPCCPhotometry() >>> data.create_daily_file(output_dir=output_dir, day=day) >>> data.build_one_epoch(raw_data_dir=path_to_data, day_of_survey=day, time_domain_dir=output_dir) |
Alternatively you can use the command line to prepare a sequence of days in one batch:
>>> build_time_domain.py -d 20 21 22 23 -p <path to raw data dir> -o <path to output time domain dir>