i.feotio2 - GRASS GIS manual

NAME

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

KEYWORDS

imagery, Moon, Clementine, Chandrayaan1, UVVIS, M3, FeO, TiO2, reflectance

SYNOPSIS

i.feotio2

i.feotio2 --help

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

Flags:

-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
--ui: Force launching GUI dialog

Parameters:

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; Options: 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)

DESCRIPTION
NOTES
REFERENCES
AUTHORS

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));
}

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.

AUTHORS

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

SOURCE CODE

Available at: i.feotio2 source code (history)