CR2RE Pipeline Reference Manual 1.6.2
hdrl_cat_areals-test.c
1/*
2 * This file is part of the HDRL
3 * Copyright (C) 2017 European Southern Observatory
4 *
5 * This program is free software; you can redistribute it and/or modify
6 * it under the terms of the GNU General Public License as published by
7 * the Free Software Foundation; either version 2 of the License, or
8 * (at your option) any later version.
9 *
10 * This program is distributed in the hope that it will be useful,
11 * but WITHOUT ANY WARRANTY; without even the implied warranty of
12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13 * GNU General Public License for more details.
14 *
15 * You should have received a copy of the GNU General Public License
16 * along with this program; if not, write to the Free Software
17 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
18 */
19
20#include <cpl_test.h>
21
22#include "../hdrl_cat_apio.h"
23#include "../hdrl_cat_areals.h"
24
25
26#define NT 117
27
28
29int main(void)
30{
31 cpl_test_init(PACKAGE_BUGREPORT,CPL_MSG_WARNING);
32
33 cpl_size x[] = {398,399,400,397,398,399,400,401,402,403,396,397,398,399,400,
34 401,402,403,404,395,396,397,398,399,400,401,402,403,404,405,
35 395,396,397,398,399,400,401,402,403,404,405,395,396,397,398,
36 399,400,401,402,403,404,405,395,396,397,398,399,400,401,402,
37 403,404,405,406,395,396,397,398,399,400,401,402,403,404,405,
38 395,396,397,398,399,400,401,402,403,404,405,395,396,397,398,
39 399,400,401,402,403,404,405,396,397,398,399,400,401,402,403,
40 404,397,398,399,400,401,402,403,398,399,400,401};
41
42 cpl_size y[] = {394,394,394,395,395,395,395,395,395,395,396,396,396,396,396,
43 396,396,396,396,397,397,397,397,397,397,397,397,397,397,397,
44 398,398,398,398,398,398,398,398,398,398,398,399,399,399,399,
45 399,399,399,399,399,399,399,400,400,400,400,400,400,400,400,
46 400,400,400,400,401,401,401,401,401,401,401,401,401,401,401,
47 402,402,402,402,402,402,402,402,402,402,402,403,403,403,403,
48 403,403,403,403,403,403,403,404,404,404,404,404,404,404,404,
49 404,405,405,405,405,405,405,405,406,406,406,406};
50
51 double z[] = {8.87152,12.515,7.69699,10.8527,22.2509,21.7368,13.0388,
52 12.1853,17.1976,7.43948,15.2245,29.1964,37.9117,57.9371,
53 71.5542,57.1288,34.7726,15.5934,11.5374,15.995,21.3606,
54 60.4006,103.46,147.55,168.274,147.476,98.9157,51.7186,20.188,
55 3.04248,5.77832,49.3103,98.2057,187.557,268.353,310.638,
56 274.295,183.969,94.6933,47.9889,20.245,26.3758,59.1781,
57 152.389,275.916,395.107,450.251,397.53,272.322,147.053,54.767,
58 11.8971,13.3888,73.3689,165.899,298.455,449.707,493.25,441.585,
59 299.31,157.474,70.1224,15.5313,8.76074,20.7188,54.5798,141.249,
60 264.87,382.736,435.452,393.871,268.175,138.485,65.9307,28.7812,
61 19.379,36.6449,93.5458,186.823,270.95,305.093,260.879,183.683,
62 100.676,32.6281,16.6497,5.94965,17.8105,57.256,106.32,145.264,
63 164.271,137.093,88.9384,60.7841,31.8582,10.0435,4.69162,
64 15.2187,32.5385,61.0381,74.5399,67.3727,43.3964,25.0956,
65 16.7595,-0.37323,21.3832,19.2497,18.5883,9.37448,19.6048,
66 11.5006,13.0159,14.5852,13.66,-1.04889};
67
68 /* Set up apm structure */
69 ap_t ap;
70 ap.lsiz = 2048;
71 ap.csiz = 2048;
72 ap.inframe = NULL;
73 ap.conframe = NULL;
74
75 /* Initialize */
76 hdrl_apinit(&ap);
77
78 ap.npl_pix = NT;
79 ap.plarray = cpl_realloc(ap.plarray, NT * sizeof(*ap.plarray));
80 for (cpl_size i = 0; i < NT; i++) {
81 ap.plarray[i].x = x[i];
82 ap.plarray[i].y = y[i];
83 ap.plarray[i].z = z[i];
84 ap.plarray[i].zsm = z[i];
85 }
86
87 ap.xintmin = 0.;
88 ap.areal_offset = 3.47165;
89 ap.thresh = 11.0936;
90 ap.fconst = 1.4427;
91
92 /* Run the test */
93 cpl_size iareal[NAREAL];
94 hdrl_areals(&ap, iareal);
95
96 cpl_test_eq(iareal[0], 104);
97 cpl_test_eq(iareal[1], 75);
98 cpl_test_eq(iareal[2], 63);
99 cpl_test_eq(iareal[3], 45);
100 cpl_test_eq(iareal[4], 25);
101 cpl_test_eq(iareal[5], 9);
102 cpl_test_eq(iareal[6], 0);
103 cpl_test_eq(iareal[7], 0);
104
105 /* Clean up */
106 hdrl_apclose(&ap);
107
108
109 return cpl_test_end(0);
110}
void hdrl_apinit(ap_t *ap)
Initialize the ap structure.
Definition: hdrl_cat_apio.c:54
void hdrl_apclose(ap_t *ap)
Close ap structure.
void hdrl_areals(ap_t *ap, cpl_size iareal[NAREAL])
Work out the areal profiles for an object.