NAME

htrdr‑planeto — simulate radiative transfer in 3D planetary atmosphere

SYNOPSIS

DESCRIPTION

htrdr‑planeto simulates the radiative transfer of a terrestrial planet in the visible or the infrared part of the spectrum. The planet's ground (option -G) can be any set of triangles with BRDFs and temperatures defined per triangle. The atmosphere is composed of a gas mixture (option -g) and a potentially empty set of aerosols (option -a). Both can have arbitrary tetrahedral meshes with per-node radiative properties. Rayleigh is used as a gas phase function. The temperature of the gas is defined on the mesh nodes. Aerosol phase functions (Henyey and Greenstein or user defined) are also defined per node.

htrdr‑planeto is mainly a renderer that calculates an image (option -i) for a given observation position (option -C). Its internal rendering algorithm is based on Monte Carlo integration, which consists for each pixel of simulating a given number of optical paths from the sensor, taking into account the phenomena of light absorption and scattering.

htrdr‑planeto offers three ways to perform spectral integration (option -s). By default, it calculates an image for the visible part of the spectrum between 380 and 780 nanometers, for the three components of the CIE 1931 XYZ color space which are then recombined to obtain the final color for each pixel. The other two methods are to explicitly define the longwave or shortwave spectral range to be integrated and continuously sample a wavelength in this range. In fact, longwave and shortwave are keywords that mean that the source of radiation is either internal or external to the medium, respectively. In shortwave, only radiance is evaluated and stored in the output image. For longwave rendering, this estimated radiance is then converted to brightness temperature and both are recorded in the image.

In htrdr‑planeto, the spatial unit 1.0 corresponds to one meter and temperatures are expressed in Kelvin. The estimated radiances are given in W/sr/m^2, except for monochromatic calculations where the calculated spectral radiance is defined in W/sr/m^2/nm.

htrdr‑planeto implements mixed parallelism. On a single computer (that is, a node), it uses shared memory parallelism while it relies on Message Passing Interface (MPI) to parallelize calculations between multiple nodes. htrdr‑planeto can therefore be launched either directly or via a process launcher such as mpirun(1) to distribute the rendering on several computers.

The options are as follows:

-a aerosol_opt[:aerosol_opt ...]

Define an aerosol. Use this option once per aerosol, and duplicate it as many times as necessary.

The aerosol options are as follows:

mesh=volumetric_mesh

Aerosol tetrahedral mesh saved in smsh(5) format.

name=string

Name assined to the aerosol.

radprop=radiative_properties

Radiatve properties of the aerosol saved in sars(5) format. Radiative properties are defined per volumetric mesh node. This file and the tetrahedral mesh (option mesh) must therefore be consistent with each other, i.e. the nodes must be listed in the same order.

phasefn=phase_functions_list

List in rnsl(5) format of phase functions to be loaded. Each phase function is saved in rnsf(5) format. The correspondence between these phase functions and the nodes of the volumetric mesh is defined in another file (option phaseids).

phaseids=per_node_phase_function

Path to the rnpfi(5) file that stores the index of the phase function to be used per volumetric mesh node. The list of phase function is defined in another file (option phasefn). Note that this file and the tetrahedral mesh (option mesh) must be consistent with each other, i.e. the nodes must be listed in the same order.

-C persp_camera_opt[:persp_camera_opt ...]

Set up a pinhole or thin-lens perspective camera.

The options for a perspective camera are as follows:

focal-dst=distance

Distance to focus on with a thin lens camera, that is, a camera whose lens-radius is not zero. The default focal distance is 1 meters.

focal-length=length

Focal length of a camera lens. It is another way to control the field of view of a thin lens camera. By default, the field of view is set through the fov parameter.

fov=angle

Vertical field of view of the camera in ]0.0, 180.0[ degrees. The default field of view is 70 degrees.

lens-radius= radius

Radius of the camera lens. A non-zero radius means that the camera is a thin lens camera while a zero radius defines a pinhole camera. The default lens radius is 0.

pos=x,y,z

Camera position. Default is 0,0,0.

tgt=x,y,z

Targeted position Default is 0,1,0.

up=x,y,z

Upward vector that the top of the camera is pointing towards. Default is 0,0,1.

-d

Write atmospheric acceleration structures to output. Each structure is saved in legacy VTK format. To divide the resulting output into N files (N > 1), each storing an acceleration structure, one can use the csplit(1) command as below:

csplit -f octree -k output %^#\ vtk% /^#\ vtk/ \
       {$(($(grep -ce "^# vtk" output)-2))}

-f

Force overwriting of output file.

-G ground_opt[:ground_opt ...]

The planet's ground.

The ground options are as follows:

brdf=brdfs_list

List in rnsl(5) format of the BRDFs to be loaded. Each BRDF is saved in mrumtl(5) format. The correspondence between these BRDFs and the triangles of the surface mesh is defined in another file (option prop).

mesh=surface_mesh

Ground triangular mesh saved in smsh(5) format.

name=string

Name assigned to the ground.

prop=surface_properties

Ground surface properties, i.e. BRDF index and temperature, both defined by triangle. The list of BRDF is defined in another file (option brdf). Note that this file and the surface mesh must be consistent (option mesh), i.e. the triangles must be listed in the same order.

-g gas_opt[:gas_opt ...]

Gas mixture.

The gas options are as follows:

mesh=volumetric_mesh

Gas tetrahedral mesh saved in smsh(5) format.

ck=correlated_k

Correlated K fof the gas saved in sck(5) format. The correlated K are defined per volumetric mesh node. This file and the tetrahedral mesh (option mesh) must therefore be consistent with each other, i.e. the nodes must be listed in the same order.

temp=temperature

Gas temperatures saved in rngt(5) format. The temperature is defined per volumetric mesh node. This file and the tetrahedral mesh (option mesh) must therefore be consistent with each other, i.e. the nodes must be listed in the same order.

-h

Display short help and exit.

-i image_opt[:image_opt ...]

Configure sensor image.

The image options are as follows:

def=widthxheight

Image definition. Default is 320x240.

spp=samples_per_pixel

Number of samples to solve the Monte Carlo estimation of each pixel. Default is 1.

-N

Precalculate tetrahedron normals. This speeds up runtime performance by calculating normals once and for all rather than re-evaluating them every time a tetrahedron is queried at a given position. In return, the memory space used to store normals increases the memory footprint.

-O accel_struct_storage

File where atmospheric acceleration structures are stored/loaded.

If accel_struct_storage does not exist, it is created and is used to store the built acceleration structures.

If accel_struct_storage exists, acceleration structures are loaded from it rather than built from scratch, resulting in significant acceleration of the preprocessing step. Note that if the data structures stored in accel_struct_storage are not as expected (that is, the input atmospheric data or construction parameters are different), an error is notified and execution is stopped, thus avoiding the use of incorrect acceleration structures.

-o output

Output file. If not defined, data is written to standard output.

-S source_opt[:source_opt ...]

Define the external source.

The source options are as follows:

lat=latitude

The latitude of the source, i.e. its angle between [-90, 90] degrees about the Y axis. The default latitude of 0 is that of the X axis.

lon=longitude

The longitude of the source, i.e. its angle between [-180, 180] degrees about the Z axis. The default longitude of 0 is that of the X axis.

dst=distance

Distance in kilometers from source to origin. Default is 0.

rad=radiance_distribution

Source radiance distribution saved in rnrl(5) format. This option is not compatible with the temperature setting of the source (option temp) which also defines its radiance distribution.

radius=real

Source radius in kilometers.

temp=temperature

Source temperature in Kelvin. When this option is set, the radiance distribution of the source is Planck, at the specified temperature. This option is not compatible with the rad option that explicitly defines the source radiance distribution.

-s spectral_opt[:spectral_opt ...]

Configure spectral integration.

The spectral integration options are as follows:

cie_xyz

Calculate the radiance for the visible part of the spectrum between [380, 780] nanometers using the XYZ CIE 1931 color matching functions. This is the default behavior.

lw=wlen_min,wlen_max

Calculate the radiance using the internal source of radiation, i.e. the radiance emitted by the medium and its boundaries (ground and space).

Calculations are performed between [wlen_min, wlen_max] nanometers according to Planck's function for a reference temperature defined as the maximum ground temperature.

If wlen_min and wlen_max are equal, the calculation is monochromatic.

sw=wlen_min,wlen_max

Calculate the radiance using the external source of radiance (option -S).

Calculations are performed between [wlen_min, wlen_max] nanometers according to the radiance distribution of the external source (see -S option)

If wlen_min and wlen_max are equal, the calculation is monochromatic.

-T optical_thickness

Optical thickness used as threshold criterion for building the acceleration structures. Default is 1.

-t

Advice on the number of threads to use. By default, htrdr‑planeto uses many threads as processor cores.

-V accel_struct_definition

Advice on the definition of the atmospheric acceleration structures. Default is 512.

-v

Make htrdr‑planeto verbose.

OUTPUT IMAGE

Images calculated by htrdr‑planeto are saved in htrdr-image(5) format. This section describes the nature and arrangement of image data depending on the type of calculation performed.

EXIT STATUS

The htrdr‑planeto utility exits 0 on success, and >0 if an error occurs.

EXAMPLES

The following command line runs htrdr‑planeto in a verbose way (option -v) to calculate an 800 by 600 pixel image by sampling 64 radiative paths per pixel (option -i) for the 3 components of the CIE 1931 XYZ color space (option -s). The external source is positioned at -45 degrees longitude and 50 degrees latitude relative to the absolute referential (option -S). The camera (option -C) looks at the origin (tgt=0,0,0) and is positioned at 1.5e7 meters along the Y axis with an image plane aligned along the Z axis (up=0,0,1). Its vertical field of view is 70 degrees. The gas of the planetary atmosphere is described by the tetrahedral mesh recorded in the gas.smsh file, while its spectral data and temperature are given by the files gas.sck and gas.rngt, respectively. Two aerosols complete the planetary atmosphere: one for clouds and one for haze. Their respective meshes are stored in the a.smsh and b.smsh files while their radiative properties are given by the a.sars and b.sars files. Finally, their phase functions are described by a set of 2 files: the a.rnsf and b.rnsf files which list the aerosol phase functions, and the a.rnpfi and b.rnpfi files which reference them by volumetric mesh node. To speed up rendering time, the normals of the tetrahedral meshes of the gas and aerosols are precalculated once and for all (option -N). The ground geometry is stored in the ground.smsh file with its triangle properties (temperature and BRDF) defined in the ground.rnsp file. The referenced BRDFs are listed in the ground.rnsl file. The definition of acceleration structures cannot exceed 512^3 and its voxels can be merged until their optical thickness is greater than 10. These structures are either reloaded from storage_cie.bin or built from scratch and stored in storage_cie.bin depending on whether that file exists or not. Finally, the calculated images are stored in the image_CIE_XYZ.ht file even if the file already exists (options -fo):

htrdr-planeto -v -N \
-i def=800x600:spp=64 \
-s cie_xyz \
-S lon=-45:lat=50:dst=1.5e8:radius=6.9e5:temp=5778 \
-C pos=0,1.5e7,0:tgt=0,0,0:up=0,0,1:fov=70 \
-g mesh=gas.smsh:ck=gas.sck:temp=gas.rngt \
-a mesh=a.smsh:radprop=a.sars:phasefn=a.rnsf:phaseids=a.rnpfi:name=clouds \
-a mesh=b.smsh:radprop=b.sars:phasefn=b.rnsf:phaseids=b.rnpfi:name=haze \
-G mesh=ground.smsh:prop=ground.rnsp:brdf=ground.rnsl:name=namek \
-V 512 -T 10 -O storage_cie.bin \
-fo image_CIE_XYZ.ht

The next command line is the same as the previous one, except that it calculates an infrared image between 10,000 nm and 20,000 nm (option -s). Note that the acceleration structure storage file is no longer the same (storage_lw.bin rather than storage_cie.bin). Indeed, the previous one records the acceleration structures for the spectral range of the CIE 1931 XYZ color space (i.e. between [380, 780] nm), while one wants to store/reload the acceleration structures for a spectral range between 10 and 20 microns. In any case, if the previous storage had been submitted, the command would have stopped with an error message, thus avoiding the use of the wrong acceleration structures (see -O option).

htrdr-planeto -v -N \
-i def=800x600:spp=64 \
-s lw=10000,20000 \
-C pos=0,1.5e7,0:tgt=0,0,0:up=0,0,1:fov=70 \
-g mesh=gas.smsh:ck=gas.sck:temp=gas.rngt \
-a mesh=a.smsh:radprop=a.sars:phasefn=a.rnsf:phaseids=a.rnpfi:name=clouds \
-a mesh=b.smsh:radprop=b.sars:phasefn=b.rnsf:phaseids=b.rnpfi:name=haze \
-G mesh=ground.smsh:prop=ground.rnsp:brdf=ground.rnsl:name=namek \
-V 512 -T 10 -O storage_lw.bin \
-fo image_infrared.ht

SEE ALSO

htrdr-image(5), mrumtl(5), rngt(5), rnpfi(5), rnrl(5), rnsf(5), sars(5), sck(5), smsh(5)

STANDARDS

International Organization for Standardization / CIE, Colorimetry - Part 1: CIE standard colorimetric observers, ISO/CIE 11664-1:2019, June 2019.

OpenMP Architecture Review Board, OpenMP C and C++ Application Interface, March 2002, version 2.0.

Message Passing Interface Forum, MPI-2: Extensions to The Message-Passing Interface, July 1997.

HISTORY

htrdr‑planeto has been developed as part of ANR-21-CE49-0020 RaD-net project.