Last chance! 50% off unlimited learning
Sale ends in
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.
profoundSegimWarp(segim_in = NULL, header_in = NULL, header_out = NULL)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.
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.
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.
Integer matrix; the remapped image using the target WCS.
This function uses the interpolation='nearest' and doscale=FALSE in magwarp.
# 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)
# }
Run the code above in your browser using DataLab