sc_optics.h

Go to the documentation of this file.

/*
 * sc_optics.h
  *
 * Programmer(s): Michael Paluszek
 *
 * Copyright (c) 2009, Princeton Satellite Systems. All rights reserved.
 *
 */
 
#ifndef SC_OPTICS
#define SC_OPTICS
 
#ifdef AS_OS_WINDOWS
#include "matrixlib.h"
#else
#include <MatrixLib/MatrixLib.h>
#endif
#include <vector>
#include <stdint.h>
using namespace std;
 
typedef struct camera_data_s {
        ml_matrix o;                            
        ml_matrix s;                            
        ml_matrix f;                            
        double k1;                                      
        double k2;                                      
        double alpha;                           
        ml_matrix R;                            
        ml_matrix T;                            
} camera_data_t;
 
 
typedef struct psf_data_s {
    ml_matrix r;                
    ml_matrix i;                
    int n;                      
    int m;                      
    double x_pixel;             
    double y_pixel;             
} psf_data_t;
 
typedef struct gaussian_psf_s : psf_data_t {
    double sigma;                               
} gaussian_psf_t;
 
typedef struct airy_psf_s : psf_data_t {
    double aperture;                    
    double wavelength;                  
    double f_length;                    
} airy_psf_t;
 
typedef struct airydefocused_psf_s : airy_psf_t {
    double defocus;     
} airydefocused_psf_t;
 
 
 
double          effective_focal_length( ml_matrix distance, ml_matrix focal_length );
double          diffraction_limit( double lambda, double aperture );
ml_matrix       angle_of_view( double focal_length, ml_matrix d );
double          angle_of_view( double focal_length, double d );
ml_matrix       angle_of_view( ml_matrix focal_length, double d );
double          focal_length_from_angle_of_view( double alpha, double d );
ml_matrix       separation_angles( ml_matrix u );
ml_matrix       pix_to_u( ml_matrix p, double f );
ml_matrix       u_to_pix( const ml_matrix& r, double f, double row = 0, double col = 0  );
double          laser_beam_radius( double d, double lambda, double z );
double          laser_power( double d, double lambda, double z, double power );
camera_data_t camera_calibration( const ml_matrix& pD, const ml_matrix& pW );
ml_matrix calibration_cube( double w, double r, int n );
ml_matrix camera_calibration_model( camera_data_t d, const ml_matrix& pW, double f );
ml_matrix imager_psf( ml_matrix  PSF( const ml_matrix& x, const ml_matrix& y, void* context), psf_data_t* d );
void imager_psf_uint16( ml_matrix PSF( const ml_matrix& x, const ml_matrix& y, void* context), psf_data_t* d, uint16_t* frame );
int any_less_than( const ml_matrix& x, double d );
ml_matrix rhs_gaussian_psf( const ml_matrix& x, const ml_matrix& y, void* context );
ml_matrix rhs_airy_psf(const ml_matrix& x,const ml_matrix& y,void* context);
ml_matrix rhs_airy_defocused_psf(const ml_matrix& x,const ml_matrix& y,void* context);
ml_matrix rand_pn(double f, int n, int m );
ml_matrix fit_PSF_to_data(const ml_matrix& r0, const ml_matrix& intensity, const ml_matrix& pixels, int nUse, int mUse);
double cost_psf(const ml_matrix& x, void* context);
 
 
class ImagePSF
{
                
        public:
        ml_matrix pixels;
        psf_data_t* psf_data;
    
        ImagePSF();
        ImagePSF(ml_matrix pixels);
        static double rhs_fit_psf(double[]);
};
 
 
#endif