def create_simple_plan(root_type: RootType, save: Save):
parameters = {
"channel": 0,
"minimum_size": 200,
"threshold": {
"name": "Manual",
"values": {
"threshold": 13000
}
},
"noise_filtering": {
"name": "Gauss",
"values": {
"dimension_type": DimensionType.Layer,
"radius": 1.0
}
},
"side_connection": False,
}
segmentation = ROIExtractionProfile(name="test",
algorithm="Lower threshold",
values=parameters)
chosen_fields = [
MeasurementEntry(
name="Segmentation Volume",
calculation_tree=Leaf(name="Volume",
area=AreaType.ROI,
per_component=PerComponent.No),
),
]
statistic = MeasurementProfile(name="base_measure",
chosen_fields=chosen_fields,
name_prefix="")
statistic_calculate = MeasurementCalculate(
channel=-1,
units=Units.µm,
measurement_profile=statistic,
name_prefix="")
tree = CalculationTree(
root_type,
[
CalculationTree(segmentation, [
CalculationTree(statistic_calculate, []),
CalculationTree(save, [])
])
],
)
return CalculationPlan(tree=tree, name="test")
def create_calculation_plan():
parameters = {
"channel": 1,
"minimum_size": 200,
"threshold": {
"name": "Base/Core",
"values": {
"core_threshold": {
"name": "Manual",
"values": {
"threshold": 30000
}
},
"base_threshold": {
"name": "Manual",
"values": {
"threshold": 13000
}
},
},
},
"noise_filtering": {
"name": "Gauss",
"values": {
"dimension_type": DimensionType.Layer,
"radius": 1.0
}
},
"side_connection": False,
"sprawl_type": {
"name": "Euclidean",
"values": {}
},
}
segmentation = ROIExtractionProfile(
name="test",
algorithm="Lower threshold with watershed",
values=parameters)
mask_suffix = MaskSuffix(name="", suffix="_mask")
chosen_fields = [
MeasurementEntry(
name="Segmentation Volume",
calculation_tree=Leaf(name="Volume",
area=AreaType.ROI,
per_component=PerComponent.No),
),
MeasurementEntry(
name="Segmentation Volume/Mask Volume",
calculation_tree=Node(
left=Leaf(name="Volume",
area=AreaType.ROI,
per_component=PerComponent.No),
op="/",
right=Leaf(name="Volume",
area=AreaType.Mask,
per_component=PerComponent.No),
),
),
MeasurementEntry(
"Segmentation Components Number",
calculation_tree=Leaf("Components number",
area=AreaType.ROI,
per_component=PerComponent.No),
),
]
statistic = MeasurementProfile(name="base_measure",
chosen_fields=chosen_fields,
name_prefix="")
statistic_calculate = MeasurementCalculate(
channel=0,
units=Units.µm,
measurement_profile=statistic,
name_prefix="")
tree = CalculationTree(
RootType.Image,
[
CalculationTree(mask_suffix, [
CalculationTree(segmentation,
[CalculationTree(statistic_calculate, [])])
])
],
)
return CalculationPlan(tree=tree, name="test")