Skip to content

Validation

correct_goes16_bands(ds)

Corrects the band coordinates in a GOES-16 dataset.

Parameters:

Name Type Description Default
ds Dataset

The input dataset containing GOES-16 bands.

 required

Returns:

Type Description
Dataset

xr.Dataset: The corrected dataset with updated band coordinates.
Source code in rs_tools/_src/geoprocessing/goes/validation.py 
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
def correct_goes16_bands(ds: xr.Dataset) -> xr.Dataset:
    """
    Corrects the band coordinates in a GOES-16 dataset.

    Parameters:
        ds (xr.Dataset): The input dataset containing GOES-16 bands.

    Returns:
        xr.Dataset: The corrected dataset with updated band coordinates.

    """
    # reassign coordinate
    band_id_attrs = ds.band_id.attrs
    ds = ds.assign_coords(band=ds.band_id.values)
    ds.band.attrs = band_id_attrs

    # drop bandid dims
    ds = ds.drop_vars("band_id")
    return ds

correct_goes16_satheight(ds)

Convert measurement angle of GOES-16 satellite data to horizontal distance (in meters).

Parameters: - ds (xr.Dataset): The input dataset containing the GOES-16 satellite data.

Returns: - xr.Dataset: The dataset with corrected perspective height.

Source code in rs_tools/_src/geoprocessing/goes/validation.py 
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
def correct_goes16_satheight(ds: xr.Dataset) -> xr.Dataset:
    """
    Convert measurement angle of GOES-16 satellite data to horizontal distance (in meters).

    Parameters:
    - ds (xr.Dataset): The input dataset containing the GOES-16 satellite data.

    Returns:
    - xr.Dataset: The dataset with corrected perspective height.
    """

    # get perspective height
    sat_height = ds.goes_imager_projection.attrs["perspective_point_height"]

    # reassign coordinates to correct height
    x_attrs = ds.x.attrs
    ds = ds.assign_coords({"x": ds.x.values * sat_height})
    ds["x"].attrs = x_attrs
    ds["x"].attrs["units"] = "meters"

    y_attrs = ds.y.attrs
    ds = ds.assign_coords({"y": ds.y.values * sat_height})
    ds["y"].attrs = y_attrs
    ds["y"].attrs["units"] = "meters"

    return ds