/usr/include/astrometry/sip-utils.h is in astrometry.net 0.46-0ubuntu2.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 | /*
This file is part of the Astrometry.net suite.
Copyright 2007, 2010 Dustin Lang.
The Astrometry.net suite is free software; you can redistribute
it and/or modify it under the terms of the GNU General Public License
as published by the Free Software Foundation, version 2.
The Astrometry.net suite is distributed in the hope that it will be
useful, but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
General Public License for more details.
You should have received a copy of the GNU General Public License
along with the Astrometry.net suite ; if not, write to the Free Software
Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
*/
#ifndef SIP_UTILS_H
#define SIP_UTILS_H
#include "sip.h"
int sip_ensure_inverse_polynomials(sip_t* sip);
/**
sip->a_order, sip->b_order, sip->a, and sip->b must be set.
sip->ap_order and sip->bp_order must be set.
Computes sip->ap and sip->bp by inverting the forward polynomial evaluated
on a regular grid of points, NX x NY, between xlo,xhi and ylo,yhi.
If NX,NY are 0, a reasonable default will be chosen.
If xlo=xhi=0 or ylo=yhi=0, the bounds of the image (from
sip->wcstan.imagew/h) will be used.
*/
int sip_compute_inverse_polynomials(sip_t* sip, int NX, int NY,
double xlo, double xhi,
double ylo, double yhi);
/*
Finds stars that are inside the bounds of a given field (wcs).
One of "sip" or "tan" must be non-NULL; if "sip" is non-NULL it is used.
One of "xyz" or "radec" must be non-NULL. If both are non-NULL, xyz is used.
"N" indicates how many elements are in these arrays. "radec" are in degrees.
If "inds" is non-NULL, the indices of stars that are inside the field are
put there; otherwise a new int array is allocated and returned; it should
be free()'d.
The pixel (xy) positions are placed into a newly-allocated array at "xy",
unless "xy" is NULL.
The number of good stars is placed in Ngood, which must be non-NULL.
*/
int* sip_filter_stars_in_field(const sip_t* sip, const tan_t* tan,
const double* xyz, const double* radec,
int N,
double** xy, int* inds, int* Ngood);
/**
Returns the bounds of the image in RA,Dec space, approximately, by
walking the boundary of the image in steps of size "stepsize".
Return values satisfy:
-90 <= decmin <= decmax <= 90
ramin <= ramax
ramin may be < 0, or ramax > 360, if the image straddles RA=0.
*/
void sip_get_radec_bounds(const sip_t* wcs, int stepsize,
double* pramin, double* pramax,
double* pdecmin, double* pdecmax);
void sip_walk_image_boundary(const sip_t* wcs, double stepsize,
void (*callback)(const sip_t* wcs, double x, double y, double ra, double dec, void* token),
void* token);
// sets RA,Dec in degrees.
void sip_get_radec_center(const sip_t* wcs,
double* p_ra, double* p_dec);
void tan_get_radec_center(const tan_t* wcs,
double* p_ra, double* p_dec);
double sip_get_radius_deg(const sip_t* wcs);
double tan_get_radius_deg(const tan_t* wcs);
// RA hours:minutes:seconds, Dec degrees:minutes:seconds
void sip_get_radec_center_hms(const sip_t* wcs,
int* rah, int* ram, double* ras,
int* decsign, int* decd, int* decm, double* decs);
anbool sip_pixel_is_inside_image(const sip_t* wcs, double px, double py);
anbool sip_is_inside_image(const sip_t* wcs, double ra, double dec);
// Writes RA, Dec H:M:S and D:M:S strings.
void sip_get_radec_center_hms_string(const sip_t* wcs,
char* rastr, char* decstr);
void sip_get_field_size(const sip_t* wcs,
double* pw, double* ph,
char** units);
void sip_shift(const sip_t* sipin, sip_t* sipout,
double xlo, double xhi, double ylo, double yhi);
anbool tan_pixel_is_inside_image(const tan_t* wcs, double px, double py);
anbool tan_is_inside_image(const tan_t* wcs, double ra, double dec);
void tan_transform(const tan_t* tanin, tan_t* tanout,
double xlo, double xhi, double ylo, double yhi,
double scale);
void tan_scale(const tan_t* tanin, tan_t* tanout, double scale);
void tan_rotate(const tan_t* tanin, tan_t* tanout, double rotate);
#endif
|