Click here to download this notebook.

Editings

[2]:

from casys.editing import Clip, ClipCondition, EditingComponent
from casys.readers import CLSTableReader

from casys import CasysPlot, DateHandler, Field, NadirData

NadirData.enable_loginfo()

Dataset definition

[3]:

# Reader definition
table_name = "TABLE_C_J3_B_GDRD"
orf_name = "C_J3_GDRD"
cycle_number = 122
start = DateHandler.from_orf(orf_name, cycle_number, 1, pos="first")
end = DateHandler.from_orf(orf_name, cycle_number, 154, pos="last")

reader = CLSTableReader(
    name=table_name,
    date_start=start,
    date_end=end,
    orf=orf_name,
    time="time",
    longitude="LONGITUDE",
    latitude="LATITUDE",
)

# Data container definition
ad = NadirData(source=reader)

Definition of the statistic

Using the add_editing method:

[4]:

condition_nan = ClipCondition(
    name="Wind speed is NaN", clip=Clip("EQ_DV(WIND_SPEED.ALTI)")
)
ed_c_nan = EditingComponent(
    name="Invalid values nan", invalidity_conditions=[condition_nan], value=1
)

condition = ClipCondition(
    name="Wind speed out of range", clip=Clip("WIND_SPEED.ALTI :> 5")
)
ed_c = EditingComponent(
    name="Invalid values", invalidity_conditions=[condition], value=2
)

var_edt = Field(name="EDITING.CUSTOM", description="My personal editing")

ad.add_editing(name="My editing", field=var_edt, components=[ed_c_nan, ed_c])

Compute

[5]:

ad.compute()

2025-05-14 11:03:27 INFO    Reading ['WIND_SPEED.ALTI', 'LONGITUDE', 'LATITUDE', 'time']
2025-05-14 11:03:29 INFO    Computing My editing
2025-05-14 11:03:29 INFO    Computing done.

Plot

Editings graphic

[6]:

cp = CasysPlot(data=ad, data_name="My editing", plot="editing")

cp.show()

[6]:

../_images/examples_example_basic_Editings_12_0.png

Editings map

[7]:

cp = CasysPlot(data=ad, data_name="My editing", plot="map")

cp.show()

[7]:

../_images/examples_example_basic_Editings_14_0.png

To learn more about editings definition, please visit this documentation page.