def bad_delta(): axis_1d = wcs.CoordAxis1D(wcs.Axis("UTC", "d")) temporal_wcs = chunk.TemporalWCS(axis_1d) # delta == 0.0 is bad ref_coord = wcs.RefCoord(float(1.0), float(2.0)) temporal_wcs.axis.function = CoordFunction1D(int(100), 0.0, ref_coord) return temporal_wcs
def get_temporal_wcs(self): axis = self.get_coord_axis1d() time = chunk.TemporalWCS(axis) if self.complete: time.exposure = 1.0 time.resolution = 2.0 time.timesys = "UTC" time.trefpos = "TOPOCENTER" time.mjdref = 3.0 return time
def get_test_function(complete, px, sx, nx, ds): axis_1d = wcs.CoordAxis1D(wcs.Axis("UTC", "d")) if complete: wcs.exposure = 300.0 wcs.resolution = 0.1 temporal_wcs = chunk.TemporalWCS(axis_1d) ref_coord = wcs.RefCoord(px, sx) temporal_wcs.axis.function = wcs.CoordFunction1D(nx, ds, ref_coord) return temporal_wcs
def bad_range_wcs(): px = float(0.5) sx = float(54321.0) nx = 200 ds = float(0.01) axis_1d = wcs.CoordAxis1D(wcs.Axis("UTC", "d")) temporal_wcs = chunk.TemporalWCS(axis_1d) temporal_wcs.exposure = 300.0 temporal_wcs.resolution = 0.1 # divide into 2 samples with a gap between c1 = wcs.RefCoord(px, sx) c2 = wcs.RefCoord(0.0, 0.0) c3 = wcs.RefCoord(px + nx * 0.66, sx + nx * ds * 0.66) c4 = wcs.RefCoord(px + nx, sx + nx * ds) temporal_wcs.axis.bounds = wcs.CoordBounds1D() temporal_wcs.axis.bounds.samples.append(wcs.CoordRange1D(c1, c3)) temporal_wcs.axis.bounds.samples.append(wcs.CoordRange1D(c4, c2)) return temporal_wcs