RFsig.h

Go to the documentation of this file.

/* © 1998 JRM Technologies, Inc.*/



/***********************************
Copyright JRM Enterprises, Inc. 1998
All rights reserved.

This code is the intellectual property of JRM Enterprises, Inc.
It may not be used or released as source or compiled binary form
without the prior written consent of JRM Enterprises, Inc.
***********************************/


#include "SigSimStructs.h"
#include "Objects.h"
#include "SigSimAPI.h"

#ifdef __cplusplus
extern "C"
{
#endif

#define NUMBER_RCS_HEADER_FIELDS                         12
#define NUMBER_RCS_LINE_FIELDS                              15

   /*// In SigSimAPI.h:
   double UnityBeamwidth2Gain(double beamwidth_rad);
   double UnityGain2Beamwidth(double gain);
   double UnityDirectivity2Gain(double directivity);
   double UnityGain2Directivity(double gain);


   double LobeBeamwidth2Gain(double beamwidth_rad);
   double LobeGain2Beamwidth(double gain);
   double LobeDirectivity2Gain(double directivity);
   double LobeGain2Directivity(double gain);

   double RectBeamwidth2Gain(double beamwidth_rad);
   double RectGain2Beamwidth(double gain);
   double RectDirectivity2Gain(double directivity);
   double RectGain2Directivity(double gain);
   */

    int FillAtObjectRFSourceFromSensor(SENSOR *sens, METIN *metin, OBJECT *object, INCIDENT_RADIANCE *incident);
    int GetAtTerrainSignatureRF(VECTOR to_view_NED, double *weight, MATL_PROP_TABLE *mp_table, VECTOR normal_NED, double surface_area_m2, ENVIRONMENTAL_INFLUENCE_MAP eimap, double *exitent_h, double *exitent_v);
    int Get_HV_Axes_NED(VECTOR NED, VECTOR *horizontal, VECTOR *vertical);
    int GetAtObjectSignatureRF(VECTOR to_view_NED, double *weight, OBJECT *object, ENVIRONMENTAL_INFLUENCE_MAP eimap, double *exitent_h, double *exitent_v);

    int FreeOLDRCSTable(OLD_RCS_TABLE *rcs);

    double Sigma0(double gamma, double A, double grazing_angle_rad);
    double UnpolarizedAbsoluteRCSatOrientationFRD(double GHz, VECTOR FRD, OLD_RCS_TABLE rcs);
    //int ReadOLDRCSFile(char *in_file, OLD_RCS_TABLE **rcs_table);
    //double BAND_UnpolarizedAbsoluteRCSatOrientationFRD(BAND band, VECTOR FRD, OLD_RCS_TABLE rcs);

   int BandAveragedArealRCS(MATL_PROP_TABLE mp_table, BAND band, double Zenith, enum Polarization PL, double *RCS, double *STD);
   int BandAveragedUlabyParams(MATL_PROP_TABLE mp_table, BAND band, enum Polarization PL, double *P1, double *P2, double *P3, double *P4, double *P5, double *P6, double *M1, double *M2, double *M3);
   int BAND_UlabyParams(MATL_PROP_TABLE mp_table, BAND band, enum Polarization PL, double *P1, double *P2, double *P3, double *P4, double *P5, double *P6, double *M1, double *M2, double *M3);
   //int BAND_EstimateRCSShininess(MATL_PROP_TABLE mp_table, BAND band);


    // BAND_TerrainRCS() provides Ulaby-Dobson terrain RCS returns for a given zenith angle, wavelength, and polarization: 
    // Inputs: band = desired waveband
    //         Zenith angle [rad] 
   //       Polarization [any]
    //         mp_table = relevant MATL_PROP_TABLE whose top layer contains spectral RCS data
    // Uses :  mp_table.{nspecrecs, reflectflag[i]=rf, polarflag[i]==PL, mp_table.lambda_um[i], param{1-9}[i]}.
    // Outputs:   *RCS = mean unpolarized terrain radar cross section (sigma_0) [dB] for the specified waveband and angle.
    //         *STD = standard deviation of the RCS value [dB] for the specified waveband and angle.
    // Assumes:   mp_table.lambda_um values are at beginning of bins, and are in increasing order for a given polarization.
    // Returns:   SIGSIM_DATA_NOT_FOUND if there exist no spectral records in mp_table for which (mp_table.reflectflag[i]==rf)&&(mp_table.polarflag[i]==PL).
    //         SIGSIM_SUCCESS upon successful completion.
    //int BAND_TerrainRCS(MATL_PROP_TABLE mp_table, BAND band, double Zenith, enum Polarization PL, double *RCS, double *STD);

    // TerrainRCS() provides Ulaby-Dobson terrain RCS returns for a given zenith angle, wavelength, and polarization: 
    // Inputs: lambda [um] = desired wavelength 
    //         Zenith angle [rad] 
   //       Polarization [any]
    //         mp_table = relevant MATL_PROP_TABLE whose top layer contains spectral RCS data
    // Uses :  mp_table.{nspecrecs, reflectflag[i]=rf, polarflag[i]==PL, mp_table.lambda_um[i], param{1-9}[i]}.
    // Outputs:   *RCS = mean unpolarized terrain radar cross section (sigma_0) [dB] for the specified wavelength and angle.
    //         *STD = standard deviation of the RCS value [dB] for the specified wavelength and angle.
    // Assumes:   mp_table.lambda_um values are at beginning of bins, and are in increasing order for a given polarization.
    // Returns:   SIGSIM_DATA_NOT_FOUND if there exist no spectral records in mp_table for which (mp_table.reflectflag[i]==rf)&&(mp_table.polarflag[i]==PL).
    //         SIGSIM_SUCCESS upon successful completion.
    //int TerrainRCS(MATL_PROP_TABLE mp_table, double lambda, double Zenith, enum Polarization PL, double *RCS, double *STD);

    double Tsigma0(double Zenith, double P1, double P2, double P3, double P4, double P5, double P6);
    double TSTDDev(double Zenith, double M1, double M2, double M3);

    // TerrainRCS_Polynomial() provides Ulaby-Dobson terrain RCS parameters for a given angle, wavelength, and polarization:
    // Inputs: lambda [um] = desired wavelength 
    //         mp_table = relevant MATL_PROP_TABLE whose top layer contains spectral RCS data
    //          PL = polarization state (use "un" for average of all polarization states).
    // Uses :  mp_table.{nspecrecs, reflectflag[i]=rf, polarflag[i]==PL, mp_table.lambda_um[i], param{1-9}[i]}.
    // Outputs:   Coefficients of polynomial RCS = RCS0 + RCS1 * sin(Zenith) + RCS2 * sin^2(Zenith) + RCS3 * sin^3(Zenith)
    //      Coefficients of polynomial STD = STD0 + STD1 * sin(Zenith) + STD2 * sin^2(Zenith)   
    // Assumes:   mp_table.lambda_um values are at beginning of bins, and are in increasing order for a given polarization.
    // Returns:   SIGSIM_DATA_NOT_FOUND if there exist no spectral records in mp_table for which (mp_table.reflectflag[i]==rf)&&(mp_table.polarflag[i]==PL)&&(lambda>mp_table.lambda_um[i]).
    //      SIGSIM_SUCCESS upon successful completion.
    //int TerrainRCS_Polynomial(MATL_PROP_TABLE mp_table, double lambda, enum Polarization PL, double *RCS0, double *RCS1, double *RCS2, double *RCS3, double *STD0, double *STD1, double *STD2);
   int BAND_TerrainRCS_Polynomial(MATL_PROP_TABLE mp_table, BAND band, enum Polarization PL, double *RCS0, double *RCS1, double *RCS2, double *RCS3, double *STD0, double *STD1, double *STD2);


   // ---------------- CW DATALINKS and EW JAMMERS ------------------
   //int AddDATALINK_PARTICIPANT(DATALINK_PARTICIPANT **c, int *nc, DATALINK_PARTICIPANT toadd);
   //void FreeDATALINK_PARTICIPANTArray(DATALINK_PARTICIPANT **c, int *nc, jrmbool init);
   //int LoadDATALINK_PARTICIPANTs(char *filepath, DATALINK_PARTICIPANT **datalink, int *ndatalinks);
   //int FindDATALINK_PARTICIPANT(DATALINK_PARTICIPANT *datalink, int ndatalinks, char *name, DATALINK_PARTICIPANT *found);
   //int AddRF_JAMMER(RF_JAMMER **c, int *nc, RF_JAMMER toadd);
   //void FreeRF_JAMMERArray(RF_JAMMER **c, int *nc, jrmbool init);
   //int LoadRF_JAMMERs(char *filepath, RF_JAMMER **rfjammer, int *nrfjammers);
   //int FindRF_JAMMER(RF_JAMMER *rfjammer, int nrfjammers, char *name, RF_JAMMER *found);

#ifdef __cplusplus
}
#endif