<div dir="ltr">Hi, all--<div><br></div><div>I have a clump that I've gotten with the clump finder.  I'd like to paint those zones of that clump with a flag, so I can exclude that region from other analysis.  I'm a little stuck on how to get this done-- I feel like it needs to involve clump.data._flush_data_to_grids and a cut_mask, but I'm a little stuck on how to define the field that I flush_to_grids.  </div>

<div><br></div><div>Specifically, how to get the cut_mask information to the field, and how to avoid clobbering the paint from one clump with another (or, from one grid with another)</div>
<div><br></div><div style>I'm trying to do something like:</div><div style><br></div><div style><code></div><div><div>for grid in clump.data._grids:</div><div>    cut_mask = clump.data._get_cut_mask(grid)</div>

<div>    def masked_region(field,data):<br></div><div>        cut_mask = data.get_field_parameter('cut_mask', None)</div><div>        try:<br></div><div>            base = data['masked_region']</div><div>
        except:</div>
<div>            base = na.zeros(data['Density'].shape)</div><div>        if cut_mask is not None:</div><div>            base = na.maximum(base,cut_mask)</div><div>        return base    </div><div>     add_field('masked_region',function=masked_region,validators=[ValidateGridType()], not_in_all=True)</div>

<div>    grid.set_field_parameter('cut_mask',cut_mask)</div><div>    print grid, cut_mask.sum(), grid['masked_region'].sum()</div></div><div style>clump.data._flush_to_grids('masked_region',0)</div>

<div style></code></div><div style><br></div><div style>but clump['masked_region'].sum() only returns the flagged zones in the first of the grids.</div><div style><br></div><div>I feel like this is something I've done before, but can't find any reference to it in my email or old scripts.  Apologies if this has been answered elsewhere, or it's obvious and I'm being thick.</div>


<div><br></div><div>Thanks!</div><div>d.</div><div><br></div><div><br></div></div>