Member clipm_align_points (const cpl_matrix *ref_points, const cpl_matrix *in_points, const cpl_matrix *ref_variances, const cpl_matrix *in_variances, clipm_align_opt align_mode_bitmask, cpl_matrix **transform_matrix, cpl_matrix **shift, cpl_matrix **residuals)

API change to accept uncertainties of points

Member clipm_aperture_characterise_circular (const cpl_image *input, const cpl_size window_xxyy[4], double *out_centre_x, double *out_centre_y, double *out_radius, double *out_r_sigma, cpl_matrix **statistics)

Export statistics computation into other function
check bad pixel handling

Member clipm_aperture_characterise_rectangular (const cpl_image *input, const cpl_size window_xxyy[4], double *centre_x, double *centre_y, double *angle, double *length, double *width, cpl_matrix **statistics)

If input is integer, convert rotated to double inside clipm_priv_image_lineartrans_bilinear()
Improve object detection step: just taking the barycentre is poor...:-(
Export statistics computation into other function
create a fitting routine to optimize the result (or alternatively, measure the angle again in the region around the rectangle, then correct)
document bad pixel handling

Member clipm_aperture_detect_circles (const cpl_image *input, const cpl_size window_xxyy[4], cpl_matrix **out_circles)

try clipm_image_signal_estimate_bg_in_region() for the Background Sigma Estimation
split up the function into more atomic parts
check bad pixel handling

Member clipm_centroiding_gauss (const cpl_image *image, const cpl_size window_xxyy[4], double *xy_centre, double *xy_centre_err, double *xy_sigma, double *xy_sigma_err, double *xy_fwhm, double *xy_fwhm_err, double *centre_intensity, int robustness)

flag to allow enlargement of window

Member clipm_centroiding_moment (const cpl_image *image, const cpl_size window_xxyy[4], int allow_wdw_enlarge, double gain, double *out_xy_centre, double *out_xy_centre_err, double *out_xy_sigma, double *out_xy_fwhm, double *centre_intensity)

CLIPM_CENTROIDING_MIN_WINDOWSIZE
replace +/- 2 fwhm by _clipm_priv_image_get_optimal_box_width() (or better, use the output from clipm_priv_image_estimate_fwhm_xy())

Member clipm_image_signal_get_barycentre (const cpl_image *img, const cpl_size window_xxyy[4], double bg_level, double lower_cutlevel, double *out_xy_centre, double *out_weight, cpl_size *out_nused)

:

test bad pixel handling in unit test

Member clipm_image_sim_circle (cpl_image *img_modified, double centre_x, double centre_y, double radius, double brightness)

implement unit test

Member clipm_image_sim_noise_gaussian (cpl_image *img_modified, const cpl_size window_xxyy[4], double sigma)

implement unit test

Member clipm_image_sim_rectangle (cpl_image *img_modified, double centre_x, double centre_y, double angle, double length, double width, double brightness)

implement unit test

Member clipm_math_rng_poisson_pointpattern_2d (double xmin, double ymin, double xmax, double ymax, cpl_size homogeneity, cpl_size N, cpl_type type)

default values?
error handling
optimize computation, and test using same random initialization
Correct eclipse bugs:
- first entry always (0,0): done
- round error for negative numbers: done
- correct rectangle limits for integers: done

Member clipm_priv_array_estimate_fwhm (const cpl_array *input, double centre, double bg_level, cpl_size *out_maxindex, double *out_middlepos, double *out_edgeslope)

Force _clipm_priv_vector_find_limits() to fail in unit test
Check out_edgeslope in unit test, also with flat object

Member clipm_priv_image_collapse (const cpl_image *image, const cpl_size window_xxyy[4], cpl_array **horizontal, cpl_array **vertical, cpl_array **x_weight_map, cpl_array **y_weight_map)

implement unit test

Member clipm_priv_image_compute_gradients (const cpl_image *input, const cpl_size window_xxyy[4], cpl_image **out_magnitudes, cpl_image **out_angles)

Implement care for bad pixels map

Member clipm_priv_image_compute_mean_angle (const cpl_image *magnitudes, const cpl_image *angles, const cpl_size window_xxyy[4], double periodic_factor, double *out_norm)

adapt cpl_error handling
ignore bad pixels: done
return also $\sqrt{|\sum/N|}$ as measurement for the quality

Member clipm_priv_image_conv_matrix (const cpl_image *input, const cpl_matrix *kernel, const cpl_size window_xxyy[4], int extend_bpm, int int2double)

:

implement unit test

Member clipm_priv_image_estimate_fwhm_xy (const cpl_image *img, const double *xy_peakpos, double bg_level, double *out_xy_fwhm, double *out_xy_middle, double *out_xy_edge_sigma)

test with flat object in unit test

Member clipm_priv_image_extract_cast (const cpl_image *input, const cpl_size *window_xxyy, cpl_type outtype)

implement unit test

Member clipm_priv_image_fill (cpl_image *image, const cpl_size *window_xxyy, double value)

implement unit test

Member clipm_priv_image_fill_rectangle (cpl_image *image, const double centre_xy[2], double angle, const double size_lw[2], double brightness, int anti_alias, int additive)

implement unit test

Member clipm_priv_image_filter_lowpass (const cpl_image *image, const cpl_size window_xxyy[4], double sigma)

implement unit test

Member clipm_priv_image_get_kappa_sigma (const cpl_image *img, const cpl_size window_xxyy[4], double kappa, double initial_limits[], int nmax_iterations, double *out_mean, double *out_sigma, double *out_kappasigma, cpl_size *out_nused, int *out_niterations)

test initial_limits in unit test, and that they can be left

Member clipm_priv_image_get_psf_sigma (const cpl_image *img, const cpl_size window_xxyy[4], const double *centre_xy, double bg_level, double cut_lower, double cut_upper, double *sigma_xy, double gain, double *centre_err_xy)

:

implement unit test

Member clipm_priv_image_lineartrans_bilinear (const cpl_image *input, const cpl_matrix *transform, cpl_matrix *shift, const cpl_size window_xxyy[4], int adjust_size, int use_outside_wdw, int norm_flag, int set_bpm_flag, cpl_image **contrib_map)

Support transforming integer into double

Member clipm_priv_image_resample_bilinear (const cpl_image *input, double scale)

Implement bad pixel handling

Member clipm_priv_image_save_debug (const cpl_image *img, const char *filename)

implement unit test

Member clipm_priv_matrix_create_corners_rectangle (const double centre_xy[2], double angle, const double size_lw[2])

implement unit test

Member clipm_priv_matrix_new_from_image_window (const cpl_image *image, const cpl_size window_xxyy[4], const int vflip)

Remove when cpl_matrix_new_from_image_window() is available.

Member clipm_priv_matrix_round (cpl_matrix *matrix)

: Remove when cpl_matrix_round() is available.

Member clipm_priv_matrix_transform_invert (const cpl_matrix *transf, const cpl_matrix *shift, cpl_matrix **inv_transf, cpl_matrix **inv_shift)

implement unit test

Member clipm_priv_optimize_gaussian (const cpl_vector *xvalues, const cpl_vector *xsigmas, const cpl_vector *yvalues, const cpl_vector *ysigmas, cpl_fit_mode fit_opts, int max_iterations, double *x0, double *x0_uncertainty, double *sigma, double *sigma_uncertainty, double *fwhm, double *fwhm_uncertainty, double *area, double *offset, double *mse, double *red_chisq, cpl_matrix **covariance, cpl_size *lower_ndx, cpl_size *upper_ndx)

Enable cpl_vector_fit_gaussian() error history again when fixed [DFS05410].
Return area and offset uncertainties.
Fix iteration window for non-equidistant xvalues.

Member clipm_priv_optimize_gaussian (const cpl_vector *xvalues, const cpl_vector *xsigmas, const cpl_vector *yvalues, const cpl_vector *ysigmas, cpl_fit_mode fit_opts, int max_iterations, double *x0, double *x0_uncertainty, double *sigma, double *sigma_uncertainty, double *fwhm, double *fwhm_uncertainty, double *area, double *offset, double *mse, double *red_chisq, cpl_matrix **covariance, cpl_size *lower_ndx, cpl_size *upper_ndx)

Enable cpl_vector_fit_gaussian() error history again when fixed [DFS05410].
Return area and offset uncertainties.
Fix iteration window for non-equidistant xvalues.

Member clipm_test_aperture_characterise_rectangular (void)

add clean failure tests

Member clipm_test_image_display_insert_gaps (void)

improve test

Member clipm_test_priv_image_estimate_fwhm_xy (void)

test object shape: step function smoothed (for out_edge_sigma)

Member clipm_test_priv_matrix_transform_invert (void)

complete test by clean failure testing

Member clipm_testlib_ctype_2_cpltype (const char *ctype_str)

Modify unit tests to drop this function.

Member component

rewrite whole unit test file