Skip to content

Rescale

convert_integers2radiances(da)

Function to convert scaled integers to radiances. Scaled integers are unitless, radiances are given in W^2/m^2/um/sr

Source code in rs_tools/_src/geoprocessing/modis/rescale.py 
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
def convert_integers2radiances(
        da: xr.DataArray, 
    ) -> xr.DataArray:
    """
    Function to convert scaled integers to radiances.
    Scaled integers are unitless, radiances are given in W^2/m^2/um/sr
    """

    radiance_scale = da.radiance_scales
    radiance_offsets = da.radiance_offsets 

    radiance_scale = np.expand_dims(radiance_scale, axis=(1,2))
    radiance_offsets = np.expand_dims(radiance_offsets, axis=(1,2))

    assert radiance_offsets.shape == radiance_scale.shape

    corrected_data = (da - radiance_offsets)*radiance_scale

    #TODO - change attributes to have units
    #TODO - change attributes to change name of variable
    #TODO - change attributes to change rescaling

    return corrected_data

convert_integers2reflectances(da)

Function to convert scaled integers to reflectances. Scaled integers and reflectances are both unitless.

Source code in rs_tools/_src/geoprocessing/modis/rescale.py 
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
def convert_integers2reflectances(
        da: xr.DataArray, 
    ) -> xr.DataArray:
    """
    Function to convert scaled integers to reflectances.
    Scaled integers and reflectances are both unitless.
    """

    reflectance_scale = da.reflectance_scales
    reflectance_offsets = da.reflectance_offsets 

    reflectance_scale = np.expand_dims(reflectance_scale, axis=(1,2))
    reflectance_offsets = np.expand_dims(reflectance_offsets, axis=(1,2))

    assert reflectance_offsets.shape == reflectance_scale.shape

    corrected_data = (da - reflectance_offsets)*reflectance_scale

    return corrected_data