DIRSIG4 Interactive Mode Guide

Interactive Mode Ready

# Simulation Parameters #
Scene Bounding Box:
    minimum: (-521.713,-543.484,-953.792)
    maximum: (1654.42,1632.64,1222.34)

# Solution Results #
View Ray:
    origin: (800,800,1000)
    direction: (0,0,-1)
Solid Angle: 1
Surface Count: 1
Path Radiance:
    [0.47 microns] 0.00279327
    [0.48 microns] 0.00281209
    [0.49 microns] 0.00275102
    [0.5 microns] 0.0028744

# Surface 1 #
Entry Distance: 859.254
Entry Angle: 0.485803
Material ID: 25

Interactive Mode Finished

<?xml version="1.0" standalone="no"?>
<!DOCTYPE interactivesession
  SYSTEM "http://www.dirsig.org/dtds/interactive_out.dtd">
<interactivesession>
  <simulationinfo>
    <simscenebbox>
	<box>
          <lowerextent>
            <point><x>-521.713</x><y>-543.484</y><z>-953.792</z></point>
          </lowerextent>
          <upperextent>
            <point><x>1654.42</x><y>1632.64</y><z>1222.34</z></point>
          </upperextent>
        </box>
    </simscenebbox>
  </simulationinfo>
  <solutioninfo>
    <pathviewray>
      <ray>
        <origin><point><x>800</x><y>800</y><z>1000</z></point></origin>
        <direction><vector><x>0</x><y>0</y><z>-1</z></vector></direction>
      </ray>
    </pathviewray>
    <pathsolidangle>1</pathsolidangle>
    <pathsurfacecount>1</pathsurfacecount>
    <pathradiance>
      <spectralvector>
        <stokesvector location="0.47"> 0.00279327 </stokesvector>
        <stokesvector location="0.48"> 0.00281209 </stokesvector>
        <stokesvector location="0.49"> 0.00275102 </stokesvector>
        <stokesvector location="0.5"> 0.0028744 </stokesvector>
      </spectralvector>
    </pathradiance>
    <surfaceinfo>
      <surfaceentrydistance>859.254</surfaceentrydistance>
      <surfaceentryangle>0.485803</surfaceentryangle>
      <surfacematid>25</surfacematid>
    </surfaceinfo>
  </solutioninfo>
</interactivesession>

<?xml version="1.0" standalone="no"?>
<!DOCTYPE interactive SYSTEM "http://www.dirsig.org/dtds/interactive.dtd">
<interactive format="binary" pausepolicy="never"
             startupoutputredirect="startup.log">

  <simulationcollectors>
  </simulationcollectors>

  <solutioncollectors>
    <pathradiance/>
  </solutioncollectors>

  <surfacecollectors>
  </surfacecollectors>
</interactive>

<bandpasslist>
  <bandpass spectralunits="microns">
    <minimum>0.3</minimum>
    <maximum>0.5</maximum>
    <delta>0.1</delta>
  </bandpass>
</bandpasslist>

0 1E-6 0 0 800 0 0 -1
0 1E-6 0 1 800 0 0 -1

$ dirsig --mode=interactive \
        --scene_filename=demo.scene \
        --atm_filename=uniform.atm \
        --collector_filename=collectors.xml \
        --bandpass_filename=bandpasses.xml \
        --datetime=2008-07-01T12:00:00.0000-05:00 < query_input.txt > query_output.txt

#include <fstream>
#include <iostream>
#include <istream>
#include <utility>
#include <vector>

typedef std::vector<double> StokesVector;
typedef std::pair<double,StokesVector> SpectralData;
typedef std::vector<SpectralData> SpectralVector;

/* Parse a stokes vector.  */
void parseSV( std::istream& is, StokesVector& stokes ) {

    // read polarization status
    char polStatus = -1;
    is.read( &polStatus, 1 );

    // number of data elements in stokes depends on polarization state
    const int numElements = ( polStatus == 1 ) ? 4 : 1;

    // read all stokes elements
    for ( int i = 0; i < numElements; i++ ) {

        // read current stokes vector element
        double stokesData = -1.0;
        is.read( reinterpret_cast<char*>( &stokesData ), sizeof( double ) );

        // stash parsed data value
        stokes.push_back( stokesData );
    }

}

/* Parse spectral vector of stokes vectors */
void parseSpecVecSV( std::istream& is, SpectralVector& specvec ) {

    // read number of stokes vectors in spectral vector
    int vectorCount = -1;
    is.read( reinterpret_cast<char*>( &vectorCount ), sizeof( int ) );

    // read each wavelength / stokes vector pair
    for ( int i = 0; i < vectorCount; i++ ) {

        double wavelength = -1.0;
        is.read( reinterpret_cast<char*>( &wavelength ), sizeof( double ) );

        StokesVector stokes;
        parseSV( is, stokes );

        specvec.push_back( SpectralData( wavelength, stokes ) );
    }
}

int main( int argc, char* argv ) {

    // Name of query output data file
    std::string dataFilename = "queryoutput";

    // Number of queries we are parsing data for
    const int numQueries = 2;

    // Open binary data file
    std::ifstream is( dataFilename.c_str(), std::ios::in | std::ios::binary );
    if ( !is ) {
        std::cerr << "Error: unable to open query result file" << std::endl;
        return 1;
    }

    // Drive through each query result
    for ( int i = 0; i < numQueries; i++ ) {

        // Parse query info
        SpectralVector pathRadiance;
        parseSpecVecSV( is, pathRadiance );

        // Pretty print query info
        std::cout << "Query # " << (i+1) << std::endl;
        for ( int l = 0; l < pathRadiance.size(); l++ ) {
            // dump wavelength
            std::cout << "  [" << pathRadiance.at(l).first << "] ";

            // dump stokes
            for ( int stokesIdx = 0;
                  stokesIdx < pathRadiance.at(l).second.size();
                  stokesIdx++ ) {
                std::cout << pathRadiance.at(l).second.at(stokesIdx) << " ";
            }
            std::cout << std::endl;
        }
        std::cout << std::endl;
    }

    return 0;
}

Query # 1
  [0.3] 0.0236427
  [0.4] 0.0235044
  [0.5] 0.00549655

Query # 2
  [0.3] 0.0236427
  [0.4] 0.0234982
  [0.5] 0.0051655

% open the input file as read-only and big endian order
fid = fopen(fname,'r','ieee-be');

% read the number of bands (values) in the spectral radiance vector
numBands = fread(fid,1,'int32');

% create storage for the spectral locations (wavelengths) and data values
specVec = zeros(numBands,1);
dataVec = zeros(numBands,1);

% read in each spectral value one at a time
for f = 1:numBands

    % read the spectral location (wavelength)
    specVec(f) = fread(fid,1,'double');

    % read the polarization flag
    polarInfo = fread(fid,1,'int8');

    % we are expecting unpolarized data, so make sure that is true
    if(polarInfo~=0)
        error('Polarized Vector')
    else
        % read in the spectral value
        dataVec(f) = fread(fid,1,'double');
    end
end

% close the file
fclose(fid);