sc_star_image.h

Go to the documentation of this file.

//
//
//  SCControl/Optics
//
//  Created by Stephanie Thomas on 7/26/17.
//
 
#ifndef sc_star_image_hpp
#define sc_star_image_hpp
 
#ifdef AS_OS_WINDOWS
#include "matrixlib.h"
#else
#include <MatrixLib/MatrixLib.h>
#endif
#include <list>
#include <stdint.h>
#include <time.h>
using namespace std;
 
#include <stdio.h>
 
typedef struct gaussian_fit_s {
    bool status;        
    double center;      
    double std_dev;     
    double height;      
    int iter;           
} gaussian_fit_t;
 
typedef struct gaussian_xyfit_s {
    bool status[2];        
    double center[2];      
    double std_dev[2];     
    double height[2];      
    int iter[2];           
} gaussian_xyfit_t;
 
ml_matrix       vm_to_intensity( ml_matrix v );
ml_matrix       intensity_to_vm( ml_matrix i );
double          vm_to_intensity( double v );
double          intensity_to_vm( double i );
void index_to_sub(int r, int c, int k , int* kR, int* kC );
ml_matrix centroid_COM( const ml_matrix& p, double threshold, int n );
ml_matrix coarse_centroid_COM( const ml_matrix& p, double threshold=0.0, int n=5 );
int edge_correction( int k0, int kC, int m );
ml_matrix zodiacal_light( const ml_matrix& r, const ml_matrix& beta );
double zodiacal_light( double r );
void apply_threshold(ml_matrix& x, double b=-1.0, double T=0.0);
ml_matrix fit_gaussian(const ml_matrix& pixels, const ml_matrix& r0, gaussian_xyfit_t& data);
 
struct star_blob {
    // bounding box
    double r0;
    double r1;
    double c0;
    double c1;
    // location of nonzero elements
    ml_matrix iPixels;
    int nPix;
    // maximum pixel in blob
    double pMax;
    // total counts in star
    double counts;
    ml_matrix iMax;
    double radius;
};
 
std::list<star_blob> Blobify(const ml_matrix& pixels, int minPix=4);
 
#endif /* sc_star_image_h */