GRASS logo

NAME

i.nightlights.intercalibration - Performs inter-satellite calibration on DMSP-OLS Nighttime Lights Time Series

KEYWORDS

imagery, inter-satellite, calibration, nighttime lights, time series, DMSP-OLS

SYNOPSIS

i.nightlights.intercalibration
i.nightlights.intercalibration --help
i.nightlights.intercalibration [-iegkznt] image=name[,name,...] suffix=suffix model=author [--overwrite] [--help] [--verbose] [--quiet] [--ui]

Flags:

-i
Print out calibration equations
-e
Evaluation based on the Normalised Difference Index
-g
Print in shell script style (currently only NDI via -e)
-k
Keep current computational region settings
-z
Exclude zero values from the analysis (retain zero cells in output)
-n
Exclude zero values from the analysis (set zero cells to NULL in output)
-t
Do not try to transfer timestamps (for input without timestamp)
--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:

image=name[,name,...] [required]
Clean average DMSP-OLS visible band digital number image(s)
suffix=suffix [required]
output file(s) suffix
Suffix for calibrated average digital number output image(s)
Default: c
model=author [required]
Calibration model
Inter-satellite calibration model for average DMSP-OLS nighttime lights time series
Options: elvidge2009, elvidge2014, liu2012, wu2013
Default: elvidge2014

Table of contents

DESCRIPTION

i.nightlights.intercalibration is a GRASS GIS module performing inter-satellite calibration on DMSP-OLS nighttime lights time series. Based on "well known" emprirical regression models, it calibrates average visible band Digital Number values.

Note, the module is still under testing. Eventually minor, but important, changes might be applied to the intercalibration process.

Overview


+----------------------------------------------------------------------+
|                                                                      |
|          +-----------------+                                         |
| DNi +--> |Calibration Model| +--> Calibrated DN                      |
|          +---^-------------+            ^                            |
|              |                          |                            |
|              |             +--Evaluation+Methods-------------------+ |
|              |             |                                       | |
|              |             | ?                                     | |
|              |             |                                       | |
|              |             +---------------------------------------+ |
|                                                                    | |
| +--Regression+Models-----------------------------------------------+ |
| |                                                                  | |
| | Elvidge, 2009/2014: DNc = C0 + C1xDN + C2xDNv2                   | |
| |                                                                  | |
| | Liu,2012: based on Elvidge's model + optimal threshold method    | |
| |                                                                  | |
| | Wu, 2014:           DNc + 1 = ax(DN + 1)^b                       | |
| |                                                                  | |
| | Others?                                                          | |
| |                                                                  | |
| +------------------------------------------------------------------+ |
|                                                                      |
+----------------------------------------------------------------------+

Details

From a review paper:

"Several methods were proposed to overcome the lack of inter-satellite calibration. These include the invariant region and the quadratic regression method proposed by Elvidge et al. [23], the second-order regression and optimal threshold method proposed by Liu et al. [24], and a power-law regression method proposed by Wu et al. [25]. Although studies based on these calibration methods showed performance improvement after the rectification [24,25], the assumption that the nighttime light remains stableover time in a particular area requires a careful choice of the invariant region manually." [Huang 2014]

References above are: [23]: [Elvidge 2009] [24]: [Liu 2012], [25]: [Wu 2013]

NOTES

?

EXAMPLES

Given all maps are imported in GRASS' data base, which are:

g.list rast pattern="F*" sep=,

F101992,F101993,F101994,F121994,F121995,F121996,F121997,F121998,F121999,F141997, F141998,F141999,F142000,F142001,F142002,F142003,F152000,F152001,F152002,F152003, F152004,F152005,F152006,F152007,F162004,F162005,F162006,F162007,F162008, F162009,F182010,F182011,F182012

the default inter-calibration, based on [Elvidge 2014], can be performed as:

i.nightlights.intercalibration image=g.list rast pattern="F*" sep=, model=wu2013 output=l

Remarks

The calibration models do not include regression coefficients for all of the yearly products. In which case, the module will fail and inform with an error message like:

i.nightlights.intercalibration image=g.list rast pattern="F??????" sep=, model=liu2012 --v

... ValueError: The selected model does not know about this combination of Satellite + Year!

Example figures

To add...

TODO

in general:

in i.nightlights.intercalibration.py:

in calibration_models.py:

another module?

REFERENCES

[Review paper(s)]

[Empirical second order regression model by Elvidge, 2009 | Y = C0 + C1*X + C2*X^2 ]

[Second order regression model & optimal threshold method by Liu, 2012]

[Non-linear, power regression model]

AUTHORS

Nikos Alexandris

SOURCE CODE

Available at: i.nightlights.intercalibration source code (history)


Main index | Imagery index | Topics index | Keywords index | Graphical index | Full index

© 2003-2018 GRASS Development Team, GRASS GIS 7.4.1svn Reference Manual