imProc.cat - EranOfek/AstroPack GitHub Wiki

Description

High-level functions for manipulation of astronomical catalogs in form of AstroCatalog objects, or such objects that reside in AstroImage objects.

Additional relevant functions are available in the AstroTable, AstroImage, and AstroWCS classes.

Functions

  • imProc.cat.allFunList - Functions and Classes list for the imProc.cat package
  • imProc.cat.affine2d - Apply affine transformation to AstroCatalog, or catalog in AstroImage
  • imProc.cat.applyProperMotion - Apply proper motion and parallax to sources in AstroCatalog object
  • imProc.cat.filterForAstrometry - Given two catalogs, match their surface density and filter sources.
  • imProc.cat.findSourceInCat - find a single source by coordinates in a list of AstroCatalog of FITS tables.
  • imProc.cat.fitPeakMultipleColumns - Given N columns with some property (e.g., S/N) fit a parabola
  • imProc.cat.getAstrometricCatalog - Get Astrometric catalog from local/external database
  • imProc.cat.insertAzAlt - Calculate and insert Az, Alt, AirMass, ParAng columns to AstroCatalog object
  • imProc.cat.insertImageVal - Insert Image values at specific positions into an AstroCatalog
  • imProc.cat.unitTest - unitTest for imProc.cat

Selected examples

To apply an affine transformation to an AstroCatalog object or AstroCatalog in AstroImage object:

% Create a random catalog
AC=AstroCatalog({rand(10,2)},'ColNames',{'X','Y'})
% apply inverse affine transformation on catalog
AC=imProc.cat.affine2d(AC, [0 1 10 5],1,'AllCols',false);

To apply proper motion to catalog - i.e., convert coordinates to a new epoch (at the same equinox):

% Get a catalog
C = catsHTM.cone_search('GAIADR2',1,1,1000,'OutType','astrocatalog');
% apply proper motion, but not parallax. From epoch 2015 to 2021
Result = imProc.cat.applyProperMotion(C, 2015,2021,'EpochInUnits','J','EpochOutUnits','J','ApplyPlx',0)

To insert Az, Alt, airmass, and parallactic angle to an AstroCatalog object:

% generate a random catalog 
AC = AstroCatalog({rand(100,2)},'ColNames',{'RA','Dec'},'ColUnits',{'rad','rad'});
% Insert Az/Alt/...
imProc.cat.insertAzAlt(AC, 'JD',2451545, 'ObsCoo',[35 32]);