Utils.c File Reference

Functions for spectra correction + Others. More...

#include <windows.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <malloc.h>
#include <math.h>
#include <io.h>
#include <fcntl.h>
#include <sys/stat.h>
#include "mgui.h"
#include "DAS_Spatram.h"
#include "dcl.h"
#include "dil.h"
#include "bil.h"
#include "DOAS.h"
#include "nrutil.h"
#include "Marq.h"

Include dependency graph for Utils.c:

Go to the source code of this file.

Macros
#define	FREERETURN   {free_vector(d,1,n);free_vector(c,1,n);return;}
#define	PI   3.1415926
#define	TINY   1.0e-25

Functions
double	CalcFluxIndex (unsigned short mis[], int dx, int dy, float te)
	calculates the Flux Index of a measurement More...
void	D_CloseVR (MOBJECT p, void *od, void *ad)
	Callback closing the status bar. More...
void	D_Print_BIL_H_ASCII (int mod, FILE *fff)
void	D_PrintBILSPH (FILE *fff)
	Print Header. More...
void	D_PrintBILSPH_BIN (int fd)
	Print on file the binary header of the measured spectrum. More...
void	D_PrintSpectrum (FILE *fff)
	Print Spectrum. More...
void	D_VRange (char *str, int val, int r, int g, int b)
	Create and display the status bar. More...
void	D_VRangeSim (int id, char *str, int val, int r, int g, int b, int x, int y)
	init progress bar during simultaneous positioning More...
int	DB_CalcFluxI (char *nomefile)
void __stdcall	DD_VRange (const char *str, int val, int Port)
int	DT_ColFormat (void)
	Callback for Column format. More...
int	DT_CorrFileEvora (void)
	Used to correct some data files. More...
int	DT_CorrHeader ()
	Used to correct the header in some data files. More...
int	DT_FluxIndex (void)
	Callback for Flux Index Calculation. More...
int	DT_LinkSpectra (void)
	For the linked spectra (..on going!!!) More...
int	Model_read (char *buff)
void	polint (float xa[], float ya[], int n, float x, float *y, float *dy)
int	polintlambda (void)
void	ShiftImg (int hp, int vp)
	Used to correct some data files. More...

Variables
bil	BIL
int	D_VOpen = 0
int	D_VRip = -1
int	DD_VOpen [16]
int	DD_VRip [16]
palette	DPAL

Detailed Description

Functions for spectra correction + Others.

Author

Daniele Bortoli

Definition in file Utils.c.

Macro Definition Documentation

FREERETURN

#define FREERETURN   {free_vector(d,1,n);free_vector(c,1,n);return;}

Definition at line 32 of file Utils.c.

PI

#define PI   3.1415926

Definition at line 33 of file Utils.c.

TINY

#define TINY   1.0e-25

Definition at line 31 of file Utils.c.

Function Documentation

CalcFluxIndex()

double CalcFluxIndex	(	unsigned short	mis[],
		int	dx,
		int	dy,
		float	te
	)

calculates the Flux Index of a measurement

Parameters

mis[]	data array
dx	horizontal pixels
dy	vertical pixels
te	exposure time

Returns

: Flux Index

Parameters

*nomefile

source file name

Definition at line 306 of file Utils.c.

Referenced by DB_CalcFluxI().

{

    int x;
    long sumdata = 0;
    double average = 0;
    double fluxi = 0;


    for(x = 200; x < dx - 200 ; x++)
    {
        sumdata += mis[ dy * dx + x];
    }
   
    average = sumdata / (double) (dx - 400);
        
    fluxi = average / (double)(te * 1000000);

    return fluxi;

}

Here is the caller graph for this function:

D_CloseVR()

void D_CloseVR	(	MOBJECT	p,
		void *	od,
		void *	ad
	)

Callback closing the status bar.

Definition at line 140 of file Utils.c.

References MShellDestroy(), and MShellUnrealize().

Referenced by D_VRange(), D_VRangeSim(), and DD_VRange().

{
    MShellUnrealize(p);
    MShellDestroy(p);
    
    
}

Here is the call graph for this function:

Here is the caller graph for this function:

D_Print_BIL_H_ASCII()

void D_Print_BIL_H_ASCII	(	int	mod,
		FILE *	fff
	)

Definition at line 621 of file Utils.c.

References spectrumheader::altitude, spectrumheader::ccdex, spectrumheader::ccdey, spectrumheader::ccdmax, spectrumheader::ccdmin, spectrumheader::ccdnr, spectrumheader::ccdsx, spectrumheader::ccdsy, spectrumheader::ccdtemp, spectrumheader::ccdtexp, spectrumheader::ccdvb, spectrumheader::ccdwx, spectrumheader::ccdwy, spectrumheader::comment, ddate::da_day, ddate::da_mon, ddate::da_year, spectrumheader::DATE, spectrumheader::dummy1, spectrumheader::dummy2, spectrumheader::dummy3, spectrumheader::dummy4, spectrumheader::filpos, spectrumheader::imgtype, spectrumheader::latitude, spectrumheader::location, spectrumheader::longitude, spectrumheader::mirpos, bil::SPH, spectrumheader::sza, dtime::ti_hour, dtime::ti_min, dtime::ti_sec, spectrumheader::TIME, and spectrumheader::wavelen.

Referenced by SaveDataOffLine(), and SaveIo().

{
    
    fprintf(fff, "LOC: %s ", BIL.SPH.location); //Name of the station
    fprintf(fff, "LAT: %.2lf ", BIL.SPH.latitude);//Station latitude
    fprintf(fff, "LON: %.2lf ", BIL.SPH.longitude);//Station latitude
    fprintf(fff, "ALT[m/asl]: %u ", BIL.SPH.altitude);//Station altitude
    fprintf(fff, "SZA: %.2lf\n", BIL.SPH.sza);//SZA

    fprintf(fff, "%d/%02d/%02d ", BIL.SPH.DATE.da_year, BIL.SPH.DATE.da_mon, BIL.SPH.DATE.da_day);
    fprintf(fff, "%02d:%02d:%02d LT ",BIL.SPH.TIME.ti_hour, BIL.SPH.TIME.ti_min,BIL.SPH.TIME.ti_sec);
    fprintf(fff, "%s ", BIL.SPH.mirpos);
    fprintf(fff, "%d ", BIL.SPH.wavelen);
    fprintf(fff, "%s ", BIL.SPH.filpos);
    fprintf(fff, "%.2lf ", BIL.SPH.ccdtemp); 
    fprintf(fff, "%.2lf ", BIL.SPH.ccdtexp);
    fprintf(fff, "%5u ", BIL.SPH.ccdmax);
    fprintf(fff, "%5u\n", BIL.SPH.ccdmin);
    
    
    // Image format
    fprintf(fff, "%d ", BIL.SPH.imgtype);
    fprintf(fff, "%d ", BIL.SPH.ccdwx);
    fprintf(fff, "%d ", BIL.SPH.ccdwy);
    fprintf(fff, "%d ", BIL.SPH.ccdsx);
    fprintf(fff, "%d ", BIL.SPH.ccdsy);
    fprintf(fff, "%d ", BIL.SPH.ccdex);
    fprintf(fff, "%d ", BIL.SPH.ccdey);
    fprintf(fff, "%d ", BIL.SPH.ccdvb);
    fprintf(fff, "%d ", BIL.SPH.ccdnr);
    fprintf(fff, "%d ", BIL.SPH.dummy1);
    fprintf(fff, "%d ", BIL.SPH.dummy2);
    fprintf(fff, "%d ", BIL.SPH.dummy3);
    fprintf(fff, "%d\n", BIL.SPH.dummy4);
    fprintf(fff, "%s\n", BIL.SPH.comment);


}

Here is the caller graph for this function:

D_PrintBILSPH()

void D_PrintBILSPH

(

FILE *

fff

)

Print Header.

Parameters

*fff	ASCII file number

Definition at line 581 of file Utils.c.

References spectrumheader::altitude, spectrumheader::ccdex, spectrumheader::ccdey, spectrumheader::ccdmax, spectrumheader::ccdmin, spectrumheader::ccdnr, spectrumheader::ccdsx, spectrumheader::ccdsy, spectrumheader::ccdtemp, spectrumheader::ccdtexp, spectrumheader::ccdvb, spectrumheader::ccdwx, spectrumheader::ccdwy, spectrumheader::comment, ddate::da_day, ddate::da_mon, ddate::da_year, spectrumheader::DATE, spectrumheader::dummy1, spectrumheader::dummy2, spectrumheader::dummy3, spectrumheader::dummy4, spectrumheader::filpos, spectrumheader::imgtype, spectrumheader::latitude, spectrumheader::location, spectrumheader::longitude, spectrumheader::mirpos, bil::SPH, spectrumheader::sza, dtime::ti_hour, dtime::ti_min, dtime::ti_sec, spectrumheader::TIME, and spectrumheader::wavelen.

Referenced by DT_CorrFileEvora(), DT_CorrHeader(), and WriteDataFile().

{
    
    fprintf(fff, "LOC: %s ", BIL.SPH.location); //Name of the station
    fprintf(fff, "LAT: %.2lf ", BIL.SPH.latitude);//Station latitude
    fprintf(fff, "LON: %.2lf ", BIL.SPH.longitude);//Station latitude
    fprintf(fff, "ALT[m/asl]: %u ", BIL.SPH.altitude);//Station altitude
    fprintf(fff, "SZA: %.2lf\n", BIL.SPH.sza);//SZA

    fprintf(fff, "%d/%02d/%02d ", BIL.SPH.DATE.da_year, BIL.SPH.DATE.da_mon, BIL.SPH.DATE.da_day);
    fprintf(fff, "%02d:%02d:%02d LT ",BIL.SPH.TIME.ti_hour, BIL.SPH.TIME.ti_min,BIL.SPH.TIME.ti_sec);
    fprintf(fff, "%s ", BIL.SPH.mirpos);
    fprintf(fff, "%d ", BIL.SPH.wavelen);
    fprintf(fff, "%s ", BIL.SPH.filpos);
    fprintf(fff, "%.2lf ", BIL.SPH.ccdtemp); 
    fprintf(fff, "%.2lf ", BIL.SPH.ccdtexp);
    fprintf(fff, "%5u ", BIL.SPH.ccdmax);
    fprintf(fff, "%5u\n", BIL.SPH.ccdmin);
    
    
    // Image format
    fprintf(fff, "%d ", BIL.SPH.imgtype);
    fprintf(fff, "%d ", BIL.SPH.ccdwx);
    fprintf(fff, "%d ", BIL.SPH.ccdwy);
    fprintf(fff, "%d ", BIL.SPH.ccdsx);
    fprintf(fff, "%d ", BIL.SPH.ccdsy);
    fprintf(fff, "%d ", BIL.SPH.ccdex);
    fprintf(fff, "%d ", BIL.SPH.ccdey);
    fprintf(fff, "%d ", BIL.SPH.ccdvb);
    fprintf(fff, "%d ", BIL.SPH.ccdnr);
    fprintf(fff, "%d ", BIL.SPH.dummy1);
    fprintf(fff, "%d ", BIL.SPH.dummy2);
    fprintf(fff, "%d ", BIL.SPH.dummy3);
    fprintf(fff, "%d\n", BIL.SPH.dummy4);
    fprintf(fff, "%s\n", BIL.SPH.comment);


}

Here is the caller graph for this function:

D_PrintBILSPH_BIN()

void D_PrintBILSPH_BIN

(

int

fd

)

Print on file the binary header of the measured spectrum.

Parameters

fd	number of destination file

Definition at line 540 of file Utils.c.

References spectrumheader::altitude, spectrumheader::ccdex, spectrumheader::ccdey, spectrumheader::ccdmax, spectrumheader::ccdnr, spectrumheader::ccdsx, spectrumheader::ccdsy, spectrumheader::ccdtemp, spectrumheader::ccdtexp, spectrumheader::ccdvb, spectrumheader::ccdwx, spectrumheader::ccdwy, spectrumheader::comment, spectrumheader::DATE, spectrumheader::dummy1, spectrumheader::dummy2, spectrumheader::dummy3, spectrumheader::dummy4, spectrumheader::filpos, spectrumheader::imgtype, spectrumheader::latitude, spectrumheader::location, spectrumheader::longitude, spectrumheader::mirpos, bil::SPH, spectrumheader::sza, spectrumheader::TIME, and spectrumheader::wavelen.

Referenced by WriteDataFile().

{


    unsigned long       info = 0;


    info =0;

    info += write(fd, BIL.SPH.location, sizeof(BIL.SPH.location));
    info += write(fd, &BIL.SPH.latitude, sizeof(BIL.SPH.latitude));
    info += write(fd, &BIL.SPH.longitude, sizeof(BIL.SPH.longitude));
    info += write(fd, &BIL.SPH.altitude, sizeof(BIL.SPH.altitude));
    info += write(fd, &BIL.SPH.sza, sizeof(BIL.SPH.sza));
    info += write(fd, &BIL.SPH.DATE, sizeof(struct ddate));
    info += write(fd, &BIL.SPH.TIME, sizeof(struct dtime));

    info += write(fd, BIL.SPH.mirpos, sizeof(BIL.SPH.mirpos));
    info += write(fd, &BIL.SPH.wavelen, sizeof(BIL.SPH.wavelen));
    info += write(fd, BIL.SPH.filpos, sizeof(BIL.SPH.filpos));

    info += write(fd, &BIL.SPH.ccdtemp, sizeof(BIL.SPH.ccdtemp));
    info += write(fd, &BIL.SPH.ccdtexp, sizeof(BIL.SPH.ccdtexp));
    info += write(fd, &BIL.SPH.ccdmax, sizeof(BIL.SPH.ccdmax));

    info += write(fd, &BIL.SPH.imgtype, sizeof(BIL.SPH.imgtype));
    info += write(fd, &BIL.SPH.ccdwx, sizeof(BIL.SPH.ccdwx));
    info += write(fd, &BIL.SPH.ccdwy, sizeof(BIL.SPH.ccdwy));
    info += write(fd, &BIL.SPH.ccdsx, sizeof(BIL.SPH.ccdsx));
    info += write(fd, &BIL.SPH.ccdsy, sizeof(BIL.SPH.ccdsy));
    info += write(fd, &BIL.SPH.ccdex, sizeof(BIL.SPH.ccdex));
    info += write(fd, &BIL.SPH.ccdey, sizeof(BIL.SPH.ccdey));
    info += write(fd, &BIL.SPH.ccdvb, sizeof(BIL.SPH.ccdvb));
    info += write(fd, &BIL.SPH.ccdnr, sizeof(BIL.SPH.ccdnr));
    info += write(fd, &BIL.SPH.dummy1, sizeof(BIL.SPH.dummy1));
    info += write(fd, &BIL.SPH.dummy2, sizeof(BIL.SPH.dummy2));
    info += write(fd, &BIL.SPH.dummy3, sizeof(BIL.SPH.dummy3));
    info += write(fd, &BIL.SPH.dummy4, sizeof(BIL.SPH.dummy4));
    info += write(fd, &BIL.SPH.comment, sizeof(BIL.SPH.comment));       

}

Here is the caller graph for this function:

D_PrintSpectrum()

void D_PrintSpectrum

(

FILE *

fff

)

Print Spectrum.

Parameters

*fff	ASCII file number

Definition at line 517 of file Utils.c.

References bil::DPLOT, IDX, IDY, and d_view::ImatGraf.

Referenced by DT_CorrFileEvora(), and DT_CorrHeader().

{

    int x, y;

    for(y = 0; y < IDY; y++)
    {
        for(x = 0; x < IDX; x++)
        {
            if((x % 10) == 0)
            {
                if(y == 0 && x == 0)
                {
                }
                else
                    fprintf(fff, "\n", NULL);
            }   
            fprintf(fff, "%5u ", BIL.DPLOT.ImatGraf[y * IDX + x]);
        }
    }
    fprintf(fff, "\n", NULL);

}

Here is the caller graph for this function:

D_VRange()

void D_VRange	(	char *	str,
		int	val,
		int	r,
		int	g,
		int	b
	)

Create and display the status bar.

Parameters

*str	string as title of the status bar
val	value of the progressive index
r	red color id (0-255)
g	green color id (0-255)
b	blu color id (0-255)

Definition at line 206 of file Utils.c.

References D_CloseVR(), D_VOpen, D_VRip, palette::DB_BG, HELV_SMALL, MCreateShell(), MCreateStatusBar(), MLoopWhileEvents(), MObjectResize(), MObjectSetBackgroundRGB(), MObjectSetCompositeHeight(), MObjectSetText(), MShellDestroy(), MShellRealize(), MShellSetWMCloseCallback(), MShellUnrealize(), MStatusBarSetPos(), MStatusBarSetSliderRGB(), SF_NO_MAXIMIZE, SF_NO_MINIMIZE, and SF_SAVEUNDER.

Referenced by AMS_DoPos(), D_Homing(), D_Positioning(), DB_CalcFluxI(), DB_Column(), DB_MirrorX(), DB_ScanMaxMin(), DemoInitMotor(), DT_CorrFileEvora(), DT_CorrHeader(), EleDatOrBinFormat(), HomeAzimFast(), HomeFiltersFast(), HomeNDFiltersFast(), and HomeTrack().

{
    char msg[16];
    static MOBJECT sts[16], shell[16];
    if(D_VRip < 0 && val == 100) 
        return;
    if(D_VRip < 0 && val > 0)
        return;
    if(val == 0)
    {
        if(D_VOpen) return;
        D_VOpen = 1;
        D_VRip++;
        shell[D_VRip] = MCreateShell(str, SF_NO_MAXIMIZE | SF_NO_MINIMIZE | SF_SAVEUNDER);
        MObjectSetBackgroundRGB(shell[D_VRip], DPAL.DB_BG[0], DPAL.DB_BG[1], DPAL.DB_BG[2]);
        sts[D_VRip] = MCreateStatusBar(shell[D_VRip], " 0 %", 100, HELV_SMALL);
        MObjectSetBackgroundRGB(sts[D_VRip], DPAL.DB_BG[0], DPAL.DB_BG[1], DPAL.DB_BG[2]);
        MStatusBarSetSliderRGB(sts[D_VRip], r, g, b);
        MObjectSetCompositeHeight(sts[D_VRip], 20);

        MShellSetWMCloseCallback(shell[D_VRip], D_CloseVR, 0L);
        //MShellRealizeXY(shell[D_VRip], ShLoc[1][0], ShLoc[1][1]);
        MShellRealize(shell[D_VRip]);
        MObjectResize(shell[D_VRip], 300, 35);
        // Chiama MGUI per modificare size shell
        MLoopWhileEvents(0);
        
    }
    else
    {
        if(val == 100)
        {
            MShellUnrealize(shell[D_VRip]);
            MShellDestroy(shell[D_VRip]);
            D_VRip--;
            D_VOpen = 0;
        }
        else
        {
            sprintf(msg, "%d %%", val);
            MObjectSetText(sts[D_VRip], msg);
            MStatusBarSetPos(sts[D_VRip], val);
            
        }
        MLoopWhileEvents(0);
        
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

D_VRangeSim()

void D_VRangeSim	(	int	id,
		char *	str,
		int	val,
		int	r,
		int	g,
		int	b,
		int	x,
		int	y
	)

init progress bar during simultaneous positioning

Definition at line 258 of file Utils.c.

References D_CloseVR(), palette::DB_BG, DD_VOpen, DD_VRip, HELV_SMALL, MCreateShell(), MCreateStatusBar(), MLoopWhileEvents(), MObjectResize(), MObjectSetBackgroundRGB(), MObjectSetCompositeHeight(), MObjectSetText(), MShellDestroy(), MShellRealizeXY(), MShellSetWMCloseCallback(), MShellUnrealize(), MStatusBarSetPos(), MStatusBarSetSliderRGB(), SF_NO_MAXIMIZE, SF_NO_MINIMIZE, and SF_SAVEUNDER.

Referenced by SD_PosAADev(), SD_Position(), and SD_PosOMUDev().

{
    char msg[16];
    static MOBJECT sts[16][16], shell[16][16];
    if(DD_VRip[id] < 0 && val == 100) 
        return;
    if(DD_VRip[id] < 0 && val > 0)
        return;
    if(val == 0)
    {
        if(DD_VOpen[id]) return;
        DD_VOpen[id] = 1;
        DD_VRip[id]++;
        shell[id][DD_VRip[id]] = MCreateShell(str, SF_NO_MAXIMIZE | SF_NO_MINIMIZE | SF_SAVEUNDER);
        MObjectSetBackgroundRGB(shell[id][DD_VRip[id]], DPAL.DB_BG[0], DPAL.DB_BG[1], DPAL.DB_BG[2]);
        sts[id][DD_VRip[id]] = MCreateStatusBar(shell[id][DD_VRip[id]], " 0 %", 100, HELV_SMALL);
        MObjectSetBackgroundRGB(sts[id][DD_VRip[id]], DPAL.DB_BG[0], DPAL.DB_BG[1], DPAL.DB_BG[2]);
        MStatusBarSetSliderRGB(sts[id][DD_VRip[id]], r, g, b);
        MObjectSetCompositeHeight(sts[id][DD_VRip[id]], 20);

        MShellSetWMCloseCallback(shell[id][DD_VRip[id]], D_CloseVR, 0L);
        MShellRealizeXY(shell[id][DD_VRip[id]], x, y);
        MObjectResize(shell[id][DD_VRip[id]], 300, 35);
        // Chiama MGUI per modificare size shell
        MLoopWhileEvents(0);        
        
    }
    else
    {
        if(val >= 100)
        {
            MShellUnrealize(shell[id][DD_VRip[id]]);
            MShellDestroy(shell[id][DD_VRip[id]]);
            DD_VRip[id]--;
            DD_VOpen[id] = 0;
        }
        else
        {
            sprintf(msg, "%d %%", val);
            MObjectSetText(sts[id][DD_VRip[id]], msg);
            MStatusBarSetPos(sts[id][DD_VRip[id]], val);
        }
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

DB_CalcFluxI()

int DB_CalcFluxI

(

char *

nomefile

)

Definition at line 328 of file Utils.c.

References AllocUSMat(), CalcFluxIndex(), spectrumheader::ccdtexp, D_VRange(), ddate::da_day, ddate::da_mon, ddate::da_year, spectrumheader::DATE, bil::DFILEM, bil::DPLOT, IDX, IDY, d_view::ImatGraf, spectrumheader::mirpos, MLoopWhileEvents(), O_BINARY, ReadHeader_BIN(), ReadIMG_BIN(), SAVECNT, bil::SPH, dtime::ti_hour, dtime::ti_min, dtime::ti_sec, spectrumheader::TIME, spectrumheader::wavelen, and d_fileman::WORKDIRECTORY.

Referenced by DT_FluxIndex(), and FluxIndexCB().

{

    FILE    *fs;
    FILE    *fd;
    char buf[512],str[512];
    int er = 0, i, k;
    char string[4] = ".FIN";
    int fb;
    double fluxindex;
    char filedatan[256];
    long numspettri;
    

// Lettura file elenco
    i=0;
    fs = fopen(nomefile, "r");
        if(fs < 0) return 1;
    while( !feof( fs ) )
    {
        fgets(str, 80, fs);
        i++;
    }
    fclose(fs);
    
    numspettri = i * SAVECNT;

    D_VRange("Flux Index Calculation",0, 0, 222, 235);


    fs = fopen(nomefile, "r");
    if(fs < 0) return 2;

    //Read the file names
    k = 0;
    while( !feof( fs ) )
    {
        fscanf(fs, "%s", filedatan);
        // build the complete file name
//      strcpy(buf, BIL.DFILEM.WORKDIRECTORY);
        strcat(buf, filedatan);

        
        fb = open(filedatan,  O_RDONLY | O_BINARY, S_IREAD );

        while( !eof( fb ) )
        {

            MLoopWhileEvents(2);
            k ++;
            if((k % 10) == 0 && k > 0)
                D_VRange(NULL, (long) k * 100 / numspettri, 0, 222, 235 );
            er = ReadHeader_BIN(fb);

            if(BIL.DPLOT.ImatGraf == 0)
                BIL.DPLOT.ImatGraf = AllocUSMat(IDX, IDY);

            ReadIMG_BIN(fb, BIL.DPLOT.ImatGraf, IDX, IDY);

            if( (strcmp(BIL.SPH.mirpos,"HorizDX") == 0) | (strcmp(BIL.SPH.mirpos,"Zenith") == 0) )
            {
                //Filter Data
                //if (BIL.SPH.ccdwx - (BIL.SPH.ccdex +1) != 0) // se i pixel del sensore sono piu' di quelli in param.ini
                //  BlindPix_1  = BIL.SPH.ccdwx - (BIL.SPH.ccdex +1) ;
            
                //FilterData(BIL.DPLOT.ImatGraf, IDX, IDY, 15, BlindPix_1);

                
                fluxindex = CalcFluxIndex( BIL.DPLOT.ImatGraf, IDX , IDY/2, (float)BIL.SPH.ccdtexp);
                
                sprintf(buf,"%d.FIN", BIL.SPH.wavelen);
                strcpy(str, BIL.DFILEM.WORKDIRECTORY);
                strcat(str,buf);
                //MMessageDialog("Full FileName:", str, "ok", NULL);

                fd = fopen(str, "a");
                    if(fd < 0) return 1;
                
                sprintf(str,"%d/%02d/%02d %02d:%02d:%02d %d %lf\n", BIL.SPH.DATE.da_year, BIL.SPH.DATE.da_mon, BIL.SPH.DATE.da_day, BIL.SPH.TIME.ti_hour, BIL.SPH.TIME.ti_min, BIL.SPH.TIME.ti_sec, BIL.SPH.wavelen, fluxindex);
                fprintf(fd, "%s", str);

                fclose(fd);
            }



        }
//      D_VRange(NULL, 100, 128, 254, 128);
        close (fb);



    }

    D_VRange(NULL, 100, 0,0,0);

    fclose(fs);

    return 0;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DD_VRange()

void __stdcall DD_VRange	(	const char *	str,
		int	val,
		int	Port
	)

Definition at line 153 of file Utils.c.

References D_CloseVR(), D_VOpen, D_VRip, palette::DB_BG, HELV_SMALL, MCreateShell(), MCreateStatusBar(), MLoopWhileEvents(), MObjectResize(), MObjectSetBackgroundRGB(), MObjectSetCompositeHeight(), MObjectSetText(), MShellDestroy(), MShellRealize(), MShellSetWMCloseCallback(), MShellUnrealize(), MStatusBarSetPos(), MStatusBarSetSliderRGB(), SF_NO_MAXIMIZE, SF_NO_MINIMIZE, and SF_SAVEUNDER.

Referenced by DB_InitCamera().

{


    char msg[16];
    static MOBJECT sts[16], shell[16];
    if(D_VRip < 0 && val == 100) 
        return;
    if(D_VRip < 0 && val > 0)
        return;
    if(val == 0)
    {
        if(D_VOpen) return;
        D_VOpen = 1;
        D_VRip++;
        shell[D_VRip] = MCreateShell(str, SF_NO_MAXIMIZE | SF_NO_MINIMIZE | SF_SAVEUNDER);
        MObjectSetBackgroundRGB(shell[D_VRip], DPAL.DB_BG[0], DPAL.DB_BG[1], DPAL.DB_BG[2]);
        sts[D_VRip] = MCreateStatusBar(shell[D_VRip], " 0 %", 100, HELV_SMALL);
        MObjectSetBackgroundRGB(sts[D_VRip], DPAL.DB_BG[0], DPAL.DB_BG[1], DPAL.DB_BG[2]);
        MStatusBarSetSliderRGB(sts[D_VRip], 0, 255, 0);
        MObjectSetCompositeHeight(sts[D_VRip], 20);

        MShellSetWMCloseCallback(shell[D_VRip], D_CloseVR, 0L);
        //MShellRealizeXY(shell[D_VRip], ShLoc[1][0], ShLoc[1][1]);
        MShellRealize(shell[D_VRip]);
        MObjectResize(shell[D_VRip], 300, 35);
        // Chiama MGUI per modificare size shell
        MLoopWhileEvents(0);        
        
    }
    else
    {
        if(val == 100)
        {
            MShellUnrealize(shell[D_VRip]);
            MShellDestroy(shell[D_VRip]);
            D_VRip--;
            D_VOpen = 0;
        }
        else
        {
            sprintf(msg, "%d %%", val);
            MObjectSetText(sts[D_VRip], msg);
            MStatusBarSetPos(sts[D_VRip], val);
        }
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

DT_ColFormat()

int DT_ColFormat

(

void

)

Callback for Column format.

Definition at line 431 of file Utils.c.

References _MAX_PATH, FB_PREFS::allow_dir_create, FB_PREFS::allow_dir_delete, FB_PREFS::allow_dir_rename, FB_PREFS::bg_image_file, FB_PREFS::button_font, FB_PREFS::change_cwd, DAS, d_fileman::DATADIRECTORY, DB_Column(), d_fileman::DBFILENAME, bil::DFILEM, das::DOption, False, FB_PREFS::font, HELV_SMALL, d_fileman::HOMEDIRECTORY, MFileBrowseExt(), optionini::savingmode, FB_PREFS::transparent_tree, and True.

Referenced by ColumnFormatCB().

{
    
    int er =0;
    char fname[128];
    char buff[_MAX_PATH];
//  int l;

    
    FB_PREFS fbp;

  fname[0] = '\0';
  buff[0] = '\0';
  
    memset (&fbp, 0, sizeof (fbp));
 
    fbp.change_cwd = True;
    fbp.allow_dir_create = True;
    fbp.allow_dir_rename = True;
    fbp.allow_dir_delete = True;
    fbp.transparent_tree = False;
    fbp.font = HELV_SMALL;
    fbp.button_font = HELV_SMALL;
    //Set Img file
    sprintf(buff, "%s\\%s", BIL.DFILEM.HOMEDIRECTORY, "sys\\Tile2.bmp" );
    fbp.bg_image_file =buff;

/*
    BIL.DFILEM.DBFILENAME[0] = 0;

    MDCreateFM_Shell(FOPEN, FDAT);  //Open File, Dat ext.

    if(BIL.DFILEM.DBFILENAME[0] == 0)
        return 1;

    DB_Column(BIL.DFILEM.DBFILENAME);
*/

                if(DAS.DOption.savingmode)
                    er = MFileBrowseExt ("Select File for Column Format", "*.bif", fname, BIL.DFILEM.DATADIRECTORY, &fbp);
                else
                    er = MFileBrowseExt ("Select File for column Format", "*.dat", fname, BIL.DFILEM.DATADIRECTORY, &fbp);
                if(er)
                {
                    sprintf(BIL.DFILEM.DBFILENAME,"%s%s",BIL.DFILEM.DATADIRECTORY, fname);
                    //l = strlen(fname);
                    //strncpy(buff, fname+l-4, 4);

                    //if( (strcmp(fname+l-4,".dat") == 0) |(strcmp(fname+l-4,".Dat") == 0) )
                    //  BIL.DFILEM.FKIND = FDAT; 
                    //else if( (strcmp(fname+l-4,".bif") == 0) | (strcmp(fname+l-4,".Bif") == 0) )
                    //  BIL.DFILEM.FKIND = FBIF;
                    //else
                    //{
                    //  MMessageDialog("DAS information", "Unsupported format!!", " Ok ", NULL);  
                    //  return 1;
                    //}
                    DB_Column(BIL.DFILEM.DBFILENAME);   
                }
                else
                    return 1;
    

    return 0;

}

Here is the call graph for this function:

Here is the caller graph for this function:

DT_CorrFileEvora()

int DT_CorrFileEvora

(

void

)

Used to correct some data files.

Definition at line 786 of file Utils.c.

References AllocUSMat(), spectrumheader::ccdex, spectrumheader::ccdey, spectrumheader::ccdsx, spectrumheader::ccdsy, spectrumheader::ccdwx, spectrumheader::ccdwy, D_PrintBILSPH(), D_PrintSpectrum(), D_VRange(), d_fileman::DBFILENAME, DeallocUSMat(), bil::DFILEM, bil::DPLOT, FELE, FOPEN, IDX, IDY, d_view::ImatGraf, d_view::ImatOrig, MDCreateFM_Shell(), MLoopWhileEvents(), ReadHeaderOnly(), ReadIMG(), SAVECNT, ShiftImg(), bil::SPH, and d_fileman::WORKDIRECTORY.

{

    FILE    *fs, *fd, *fc;
    char buf[256],str[256];
    int er = 0, i, k;

    char filedatan[32];
    long numspettri;
    

    BIL.DFILEM.DBFILENAME[0] = 0;
    MDCreateFM_Shell(FOPEN, FELE);  //Open File, Ele ext.

    if(BIL.DFILEM.DBFILENAME[0] == 0)
        return 0;

// Lettura file elenco
    i=0;
    fs = fopen(BIL.DFILEM.DBFILENAME, "r");
        if(fs < 0) return 1;
    //Prima lettura file elenco per determinare approssimativamente in numero di spettri presenti   
    while( !feof( fs ) )
    {
        fgets(str, 80, fs);
        i++;
    }
    fclose(fs);
    
    numspettri = i * SAVECNT;

    D_VRange("Header file correction",0, 0, 222, 235);


    fs = fopen(BIL.DFILEM.DBFILENAME, "r");
    if(fs < 0) return 1;

    //Read the file names
    k = 0;
    while( !feof( fs ) )
    {
        fscanf(fs, "%s", filedatan);
        // build the complete file name
        strcpy(buf, BIL.DFILEM.WORKDIRECTORY);
        strcat(buf, filedatan);
    
        fc = fopen(buf, "r");
        if(fc < 0) return 1;

        while( !feof( fc ) )
        {
            MLoopWhileEvents(0);
            k ++;
            if((k % 10) == 0 && k > 0)
                D_VRange(NULL, (long) k * 100 / numspettri, 0, 222, 235 );
            er = ReadHeaderOnly(fc);
            //Modifica alcuni parametri del Header
            BIL.SPH.ccdwx = 1024;
            BIL.SPH.ccdwy = 254;
            BIL.SPH.ccdsx = 0;
            BIL.SPH.ccdsy = 0;
            BIL.SPH.ccdex = 1023;
            BIL.SPH.ccdey = 253;


            if(BIL.DPLOT.ImatGraf == 0)
                BIL.DPLOT.ImatGraf = AllocUSMat(IDX, IDY);
            
            

            ReadIMG(fc, BIL.DPLOT.ImatGraf, IDX, IDY);
            
            
            ShiftImg(IDX, IDY);

    
                sprintf(buf,"%s\\%s", "CORRECTED",filedatan);
                strcpy(str, BIL.DFILEM.WORKDIRECTORY);
                strcat(str,buf);
                //MMessageDialog("Full FileName:", str, "ok", NULL);
                fd = fopen(str, "a");
                    if(fd < 0) return 1;
                
                //Print
                D_PrintBILSPH(fd);
                D_PrintSpectrum(fd);
                    
                fclose(fd);
            

        }
        fclose(fc);

    }

    D_VRange(NULL, 100, 0,0,0);

    fclose(fs);

    if(!BIL.DPLOT.ImatGraf == 0)
    {
        BIL.DPLOT.ImatGraf = DeallocUSMat(BIL.DPLOT.ImatGraf);
        BIL.DPLOT.ImatOrig = DeallocUSMat(BIL.DPLOT.ImatOrig);
    }
    return 0;
}

Here is the call graph for this function:

DT_CorrHeader()

int DT_CorrHeader

(

)

Used to correct the header in some data files.

Definition at line 662 of file Utils.c.

References AllocUSMat(), spectrumheader::ccdex, spectrumheader::ccdey, spectrumheader::ccdsx, spectrumheader::ccdsy, spectrumheader::ccdwx, spectrumheader::ccdwy, D_PrintBILSPH(), D_PrintSpectrum(), D_VRange(), d_fileman::DBFILENAME, DeallocUSMat(), bil::DFILEM, bil::DPLOT, FELE, FOPEN, IDX, IDY, d_view::ImatGraf, d_view::ImatOrig, MDCreateFM_Shell(), MLoopWhileEvents(), ReadHeaderOldF(), ReadIMG(), SAVECNT, bil::SPH, and d_fileman::WORKDIRECTORY.

{

    FILE    *fs, *fd, *fc;
    char buf[256],str[256];
    int er = 0, i, k;

    char filedatan[32];
    long numspettri;
    

    BIL.DFILEM.DBFILENAME[0] = 0;
    MDCreateFM_Shell(FOPEN, FELE);  //Open File, Ele ext.

    if(BIL.DFILEM.DBFILENAME[0] == 0)
        return 0;

// Lettura file elenco
    i=0;
    fs = fopen(BIL.DFILEM.DBFILENAME, "r");
        if(fs < 0) return 1;
    //Prima lettura file elenco per determinare approssimativamente in numero di spettri presenti   
    while( !feof( fs ) )
    {
        fgets(str, 80, fs);
        i++;
    }
    fclose(fs);
    
    numspettri = i * SAVECNT;

    D_VRange("Header file correction",0, 0, 222, 235);


    fs = fopen(BIL.DFILEM.DBFILENAME, "r");
    if(fs < 0) return 1;

    //Read the file names
    k = 0;
    while( !feof( fs ) )
    {
        fscanf(fs, "%s", filedatan);
        // build the complete file name
        strcpy(buf, BIL.DFILEM.WORKDIRECTORY);
        strcat(buf, filedatan);
    
        fc = fopen(buf, "r");
        if(fc < 0) return 1;

        while( !feof( fc ) )
        {
            MLoopWhileEvents(0);
            k ++;
            if((k % 10) == 0 && k > 0)
                D_VRange(NULL, (long) k * 100 / numspettri, 0, 222, 235 );
            er = ReadHeaderOldF(fc);
            //Modifica alcuni parametri del Header
            BIL.SPH.ccdwx = 1024;
            BIL.SPH.ccdwy = 254;
            BIL.SPH.ccdsx = 0;
            BIL.SPH.ccdsy = 0;
            BIL.SPH.ccdex = 1023;
            BIL.SPH.ccdey = 253;

            if(BIL.DPLOT.ImatGraf == 0)
                BIL.DPLOT.ImatGraf = AllocUSMat(IDX, IDY);
            
            ReadIMG(fc, BIL.DPLOT.ImatGraf, IDX, IDY);

                
                sprintf(buf,"%s\\%s", "CORRECTED",filedatan);
                strcpy(str, BIL.DFILEM.WORKDIRECTORY);
                strcat(str,buf);
                //MMessageDialog("Full FileName:", str, "ok", NULL);
                fd = fopen(str, "a");
                    if(fd < 0) return 1;
                
                //Print
                D_PrintBILSPH(fd);
                D_PrintSpectrum(fd);
                    
                fclose(fd);
            

        }
        fclose(fc);

    }

    D_VRange(NULL, 100, 0,0,0);

    fclose(fs);

    if(!BIL.DPLOT.ImatGraf == 0)
    {
        BIL.DPLOT.ImatGraf = DeallocUSMat(BIL.DPLOT.ImatGraf);
        BIL.DPLOT.ImatOrig = DeallocUSMat(BIL.DPLOT.ImatOrig);
    
    }
    return 0;
}

Here is the call graph for this function:

DT_FluxIndex()

int DT_FluxIndex

(

void

)

Callback for Flux Index Calculation.

Definition at line 499 of file Utils.c.

References DB_CalcFluxI(), d_fileman::DBFILENAME, bil::DFILEM, FELE, FOPEN, and MDCreateFM_Shell().

Referenced by DB_ToolsKeyB_CB().

{
    
    

    BIL.DFILEM.DBFILENAME[0] = 0;
    MDCreateFM_Shell(FOPEN, FELE);  //Open File, Ele ext.

    if(BIL.DFILEM.DBFILENAME[0] == 0)
        return 1;

    DB_CalcFluxI(BIL.DFILEM.DBFILENAME);


    
    return 0;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DT_LinkSpectra()

int DT_LinkSpectra

(

void

)

For the linked spectra (..on going!!!)

Definition at line 915 of file Utils.c.

References d_fileman::DBFILENAME, bil::DFILEM, FMOD, FOPEN, MDCreateFM_Shell(), and Model_read().

Referenced by DB_ToolsKeyB_CB().

{

    BIL.DFILEM.DBFILENAME[0] = 0;
    MDCreateFM_Shell(FOPEN, FMOD);  //Open File, mod ext.

    if(BIL.DFILEM.DBFILENAME[0] == 0)
        return 1;

    Model_read(BIL.DFILEM.DBFILENAME);
    
    return 0;


}

Here is the call graph for this function:

Here is the caller graph for this function:

Model_read()

int Model_read

(

char *

buff

)

Definition at line 893 of file Utils.c.

References d_fileman::DBFILENAME, and bil::DFILEM.

Referenced by DT_LinkSpectra().

{

    FILE    *fs; 
    

    fs = fopen(BIL.DFILEM.DBFILENAME, "r");
    if(fs < 0) return 1;

    while( !feof( fs ) )
    {
//      fscanf(fs, "%s", str); //comment


    }


    fclose(fs);
    return 0;
}

Here is the caller graph for this function:

polint()

void polint	(	float	xa[],
		float	ya[],
		int	n,
		float	x,
		float *	y,
		float *	dy
	)

Definition at line 42 of file Utils.c.

References free_vector(), nrerror(), and vector().

Referenced by polintlambda().

{
    int i,m,ns=1;
    float den,dif,dift,ho,hp,w;
    float *c,*d;
    dif=(float)fabs(x-xa[1]);
    c=vector(1,n);
    d=vector(1,n);
    for (i=1;i<=n;i++)
    {                       //Here we .nd the index ns of the closest table entry,
        if((dift=(float)fabs(x-xa[i])) < dif)
        {
            ns=i;
            dif=dift;
        }
        c[i]=ya[i]; //and initialize the tableau of c’s and d’s.
        d[i]=ya[i];
    }
    *y=ya[ns--]; //This is the initial approximation to y.
    for (m=1;m<n;m++)
    {                   //For each column of the tableau,
        for(i=1;i<=n-m;i++)
        {           //we loop over the current c’s and d’s and update them.     
            ho=xa[i]-x;
            hp=xa[i+m]-x;
            w=c[i+1]-d[i];
            if((den=ho-hp) == 0.0) //This error can occur only if two input 
                                    //xa’s are (to within roundo.) identical.
                nrerror("Error in routine polint");
            den=w/den;
            d[i]=hp*den; //Here the c’s and d’s are updated.
            c[i]=ho*den;
        }
        *y += (*dy=(2*ns < (n-m) ? c[ns+1] : d[ns--]));
        /*
        After each column in the tableau is completed, we decide which correction,
        c or d, we want to add to our accumulating value of y,
        i.e., which path to take through the tableau forking up or down.
        We do this in such a way as to take the most  straight
        line route through the tableau to its apex, updating ns accordingly ù       to keep track of where we are. 
        This route keeps the partial approximations centered (insofar as possible)
        on the target x. The last dy added is thus the error indication.
        */
    }
    free_vector(d,1,n);
    free_vector(c,1,n);
}

        to keep track of where we are. 
        This route keeps the partial approximations centered (insofar as possible)
        on the target x. The last dy added is thus the error indication.
        */
    }
    free_vector(d,1,n);
    free_vector(c,1,n);
}

Here is the call graph for this function:

Here is the caller graph for this function:

polintlambda()

int polintlambda

(

void

)

Definition at line 94 of file Utils.c.

References dvector(), nrerror(), and polint().

{

    int i, l = 1, linit;
    float dy,y;
    FILE *nf;
    float *xa,*ya;
    
    xa = (float *)dvector(1,137);
    ya = (float *)dvector(1,137) ;


    if((nf = fopen("In.dat", "r")) == NULL)
        nrerror("In.dat not found\n");
    while(!feof(nf))
    {
        fscanf(nf, "%f%f", &xa[l], &ya[l]);
        printf("%d %12.6f %12.6f \n", l, xa[l],ya[l]);

        l++;
    }
    fclose(nf);
    linit = (int) xa[1];

    nf = fopen("Is.dat", "w");

    for (i=1; i<l-2;i++)
    {
        polint( &xa[i],  &ya[i],2,(float) i+linit,&y,&dy);
        printf("\n%d %12.6f %12.6f ", i, xa[i],ya[i]);
        printf("\n%d %d %12.6f %12.6f", i, i+linit,y,dy);
        printf("\n**************");

        fprintf(nf, "%12.6lf %12.6lf %d  %12.6lf \n",xa[i],ya[i], i+linit, y);

    }

    fclose(nf);

    return 0;

}

Here is the call graph for this function:

ShiftImg()

void ShiftImg	(	int	hp,
		int	vp
	)

Used to correct some data files.

Parameters

hp	horizontal pixels
vp	vertical pixels

Definition at line 765 of file Utils.c.

References bil::DPLOT, IDX, and d_view::ImatGraf.

Referenced by DT_CorrFileEvora().

{

    int x, y;

    for(y = 0; y < vp; y++)
    {
        for(x = 0; x < hp; x++)
        {
            if((x==0) & (y!=0))
                BIL.DPLOT.ImatGraf[y * IDX + x] =BIL.DPLOT.ImatGraf[y * IDX + x+2]/2 + BIL.DPLOT.ImatGraf[y * IDX + x+3]/2;
            else if((x==1) & (y!=0))
                BIL.DPLOT.ImatGraf[y * IDX + x] =BIL.DPLOT.ImatGraf[y * IDX + x+1]/2 + BIL.DPLOT.ImatGraf[y * IDX + x+3]/2;
            else
                BIL.DPLOT.ImatGraf[y * IDX + x] =BIL.DPLOT.ImatGraf[y * IDX + x];
        }
    }


}

Here is the caller graph for this function:

Variable Documentation

BIL

bil BIL

Definition at line 27 of file Utils.c.

D_VOpen

int D_VOpen = 0

Definition at line 148 of file Utils.c.

Referenced by D_VRange(), and DD_VRange().

D_VRip

int D_VRip = -1

Definition at line 148 of file Utils.c.

Referenced by D_VRange(), and DD_VRange().

DD_VOpen

int DD_VOpen[16]

Definition at line 151 of file Utils.c.

Referenced by D_VRangeSim(), SD_PosAADev(), SD_Position(), and SD_PosOMUDev().

DD_VRip

int DD_VRip[16]

Definition at line 150 of file Utils.c.

Referenced by D_VRangeSim().

DPAL

palette DPAL

Definition at line 28 of file Utils.c.