i.feotio2

Calculates the FeO or TiO2 contents from the Clementine project Moon data.

Command linePython (grass.script)

i.feotio2 [-t] band0=name band1=name output=name [param0=float] [param1=float] eqname=type [--overwrite] [--verbose] [--quiet] [--qq] [--ui]

Example:

i.feotio2 band0=name band1=name output=name eqname=feolucey2000

grass.script.run_command("i.feotio2", band0, band1, output, param0=None, param1=None, eqname="feolucey2000", flags=None, overwrite=False, verbose=False, quiet=False, superquiet=False)

Example:

gs.run_command("i.feotio2", band0="name", band1="name", output="name", eqname="feolucey2000")

Parameters

Command linePython (grass.script)

band0=name [required]
    reflectance band at 750 nm (UVVIS 2, M3 9)
band1=name [required]
    reflectance band at 950 nm (UVVIS 4 for FeO/OMAT, M3 19) or at 415 nm (UVVIS 1 for TiO2)
output=name [required]
    Name for output raster map
param0=float
    Value of theta (FeO) or y0Ti (TiO2)
param1=float
    Value of s0Ti (TiO2)
eqname=type [required]
    Type of equation
    Allowed values: feolucey2000, feolawrence2002, feowilcox2005_setparam, feowilcox2005, feozhang2013, omatlucey2000, omatwilcox2005, tio2lucey2000_setparam, tio2lucey2000
    Default: feolucey2000
    feolucey2000: FeO (Lucey et al., 2000)
    feolawrence2002: FeO (lawrence et al., 2002)
    feowilcox2005_setparam: FeO (Wilcox et al., 2005) set parameter theta
    feowilcox2005: FeO (Wilcox et al., 2005)
    feozhang2013: FeO (Zhang et al., 2013)
    omatlucey2000: OMAT (Lucey et al., 2000)
    omatwilcox2005: OMAT (Wilcox et al., 2005)
    tio2lucey2000_setparam: TiO2 (Lucey et al., 2000) set parameters
    tio2lucey2000: TiO2 (Lucey et al., 2000)
-t
    TiO2 (refl@415, refl@750 & y0Ti, s0Ti)
--overwrite
    Allow output files to overwrite existing files
--help
    Print usage summary
--verbose
    Verbose module output
--quiet
    Quiet module output
--qq
    Very quiet module output
--ui
    Force launching GUI dialog

band0 : str, required
    reflectance band at 750 nm (UVVIS 2, M3 9)
    Used as: input, raster, name
band1 : str, required
    reflectance band at 950 nm (UVVIS 4 for FeO/OMAT, M3 19) or at 415 nm (UVVIS 1 for TiO2)
    Used as: input, raster, name
output : str, required
    Name for output raster map
    Used as: output, raster, name
param0 : float, optional
    Value of theta (FeO) or y0Ti (TiO2)
param1 : float, optional
    Value of s0Ti (TiO2)
eqname : str, required
    Type of equation
    Used as: type
    Allowed values: feolucey2000, feolawrence2002, feowilcox2005_setparam, feowilcox2005, feozhang2013, omatlucey2000, omatwilcox2005, tio2lucey2000_setparam, tio2lucey2000
    feolucey2000: FeO (Lucey et al., 2000)
    feolawrence2002: FeO (lawrence et al., 2002)
    feowilcox2005_setparam: FeO (Wilcox et al., 2005) set parameter theta
    feowilcox2005: FeO (Wilcox et al., 2005)
    feozhang2013: FeO (Zhang et al., 2013)
    omatlucey2000: OMAT (Lucey et al., 2000)
    omatwilcox2005: OMAT (Wilcox et al., 2005)
    tio2lucey2000_setparam: TiO2 (Lucey et al., 2000) set parameters
    tio2lucey2000: TiO2 (Lucey et al., 2000)
    Default: feolucey2000
flags : str, optional
    Allowed values: t
    t
        TiO2 (refl@415, refl@750 & y0Ti, s0Ti)
overwrite: bool, optional
    Allow output files to overwrite existing files
    Default: False
verbose: bool, optional
    Verbose module output
    Default: False
quiet: bool, optional
    Quiet module output
    Default: False
superquiet: bool, optional
    Very quiet module output
    Default: False

DESCRIPTION

i.feo_tio2 calculates the FeO or TiO2 contents from Clementine data. Clementine UVVIS sensor has five bands, UVVIS1 = 415, UVVIS2 = 750, UVVIS3 = 900, UVVIS4 = 950 and UVVIS5 = 1000 nm.

wt%TiO2 = 3.708 arctan ((R415/R750)-y0Ti)/(R750-s0Ti) (1)

wt%FeO = -137.97 ( R750 sin(theta) + (R950/R750) cos(theta) ) + 57.46   (2)

with theta = 1.3885 rad, the 'average slope of the trends in the mare' from Wilcox et al (2005). Map-a-Planet Explorer mentions this set of equations:

double feolucey2000(double uvvis2, double uvvis4){
     //\cite{lucey2000lunar}
     return(17.427*(-atan2f(((uvvis4/uvvis2)-1.19)/(uvvis2-0.08)))-7.565);
}

double feolawrence2002(double uvvis2, double uvvis4){
     //\cite{lawrence2002iron}
     return(5.7*((-0.147+0.372*(-(uvvis4/uvvis2-1.22)/(uvvis2-0.04))+(-0.036)*pow((-(uvvis4/uvvis2-1.22)/(uvvis2-0.04)),2)))+2.15);
}

double feowilcox2005(double uvvis2, double uvvis4){
      //\cite{wilcox2005mapping}
      return(-137.97*((uvvis2*0.9834)+(uvvis4/uvvis2*0.1813))+57.46);
}

double omatlucey2000(double uvvis2, double uvvis4){
     //\cite{lucey2000lunar}
     return(sqrtf(pow((uvvis2-0.08),2)+pow(((uvvis4/uvvis2)-1.19),2)));
}

double omatwilcox2005(double uvvis2, double uvvis4){
     //\cite{wilcox2005mapping}
     return((uvvis2*0.1813)-((uvvis4/uvvis2)*0.9834));
}

double tio2lucey2000(double uvvis1, double uvvis2){
     //\cite{lucey2000lunar}
     return(3.708*pow((atan2f(((uvvis1/uvvis2)-0.42)/(uvvis2-0.0))),5.979));
}

NOTES

Initially created for Clementine data.

REFERENCES

(1) Lucey, P.G., Blewett, D.T., Jolliff, B.L., 2000. Lunar iron and titanium abundance algorithms based on final processing of Clementine ultraviolet-visible images. J. Geophys. Res. 105(E8): 20297-20305. (2) Wilcox, B.B., Lucey, P.G., Gillis, J.J., 2005. Mapping iron in the lunar mare: An improved approach. J. Geophys. Res. 110(E11):2156-2202. (3) Lawrence, DJ and Feldman, WC and Elphic, RC and Little, RC and Prettyman, TH and Maurice, S and Lucey, PG and Binder, AB, 2002. Iron abundances on the lunar surface as measured by the Lunar Prospector gamma-ray and neutron spectrometers. Journal of Geophysical Research: Planets (1991--2012), 107(E12):13-1.

AUTHOR

Yann Chemin (B.Sc. student), Birkbeck, University of London.

SOURCE CODE

Available at: i.feotio2 source code (history)
Latest change: Thursday Mar 20 21:36:57 2025 in commit 7286ecf