profoundSegimWarp: Remap Segmentation Map via Warping

Description

Remaps an input segmentation map WCS Tan Gnomonic projection system to a different target WCS. This uses magwarp with sensible settings, but magwarp can be used more directly if the other lower level options are required. This interface should cover most practical use cases though. Using profoundProFound with a remapped segmentation map is likely to be more sensible than remapping image flux since it will not produce flux interpolation errors.

Usage

profoundSegimWarp(segim_in = NULL, header_in = NULL, header_out = NULL)

Arguments

segim_in

Integer matrix; required, the segmentation map we want to remap. If segim_in is a list as created by readFITS, read.fits of magcutoutWCS then the image part of the list is parsed to segim_in and the correct header part is passed to header_in.

header_in

Full FITS header in table or vector format. This should be the header WCS that matches segim_in. Legal table format headers are provided by the read.fitshdr function or the hdr list output of read.fits in the astro package; the hdr output of readFITS in the FITSio package or the header output of magcutoutWCS. If header_in is provided then key words will be taken from here as a priority. Missing header keywords are printed out and other header option arguments are used in these cases.

header_out

Full FITS header in table or vector format. This is the target WCS projection that segim_in will be mapped onto. Legal table format headers are provided by the read.fitshdr function or the hdr list output of read.fits in the astro package; the hdr output of readFITS in the FITSio package or the header output of magcutoutWCS. If header_out is provided then key words will be taken from here as a priority. Missing header keywords are printed out and other header option arguments are used in these cases.

Value

Integer matrix; the remapped image using the target WCS.

Details

This function uses the interpolation='nearest' and doscale=FALSE in magwarp.

Examples

Run this code

# NOT RUN {
VST_r=readFITS(system.file("extdata", 'VST_r.fits', package="magicaxis"))
GALEX_NUV=readFITS(system.file("extdata", 'GALEX_NUV.fits', package="magicaxis"))

profound_KiDS=profoundProFound(VST_r, sky=0, skycut=1, sigma=2, tolerance=8, plot=TRUE)

segimFUV=profoundSegimWarp(profound_KiDS$segim, profound_KiDS$header, GALEX_NUV$hdr)

profoundSegimPlot(GALEX_NUV, segim = segimFUV)

profound_GALEX=profoundProFound(GALEX_NUV, segim=segimFUV, plot=TRUE)
# }