Method SetValuesInRadius

src/readybase/ImageRD.cpp:1617–1678

// --------------------------------------------------------------------------------

void ImageRD::SetValuesInRadius(float x,float y,float z,float r,float val,const Properties& render_settings)
{
    const int X = this->GetX();
    const int Y = this->GetY();
    const int Z = this->GetZ();

    // which chemical was clicked-on?
    float offset_x = 0.0f;
    bool show_multiple_chemicals = render_settings.GetProperty("show_multiple_chemicals").GetBool();
    int iChemical;
    if(show_multiple_chemicals && this->GetArenaDimensionality()==1)
    {
        // detect which chemical was drawn on from the click position
        const double image_height = X / this->image_ratio1D;
        iChemical = int(floor((- y + this->image_top1D + image_height)/(image_height*2)));
        iChemical = min(this->GetNumberOfChemicals()-1,max(0,iChemical)); // clamp to allowed range (just in case)
    }
    else if(show_multiple_chemicals && this->GetArenaDimensionality()>=2)
    {
        // detect which chemical was drawn on from the click position
        const float x_gap = this->x_spacing_proportion * this->GetX();
        iChemical = int(floor((x + x_gap  / 2) / (X + x_gap)));
        iChemical = min(this->GetNumberOfChemicals()-1,max(0,iChemical)); // clamp to allowed range (just in case)
        offset_x = iChemical * (X + x_gap);
    }
    else
    {
        // only one chemical is shown, must be that one
        iChemical = IndexFromChemicalName(render_settings.GetProperty("active_chemical").GetChemical());
    }

    double *dataset_bbox = this->images.front()->GetBounds();
    r *= hypot3(dataset_bbox[1]-dataset_bbox[0],dataset_bbox[3]-dataset_bbox[2],dataset_bbox[5]-dataset_bbox[4]);

    int ix,iy,iz;
    ix = int(floor(x-offset_x));
    iy = int(floor(y));
    iz = int(floor(z));
    ix = min(X-1,max(0,ix));
    iy = min(Y-1,max(0,iy));
    iz = min(Z-1,max(0,iz));

    for(int tz=max(0,int(iz-r));tz<=min(Z-1,int(iz+r));tz++)
    {
        for(int ty=max(0,int(iy-r));ty<=min(Y-1,int(iy+r));ty++)
        {
            for(int tx=max(0,int(ix-r));tx<=min(X-1,int(ix+r));tx++)
            {
                if(hypot3(ix-tx,iy-ty,iz-tz)<r)
                {
                    float old_val = this->GetImage(iChemical)->GetScalarComponentAsFloat(tx,ty,tz,0);
                    int ijk[3] = { tx, ty, tz };
                    vtkIdType iCell = this->GetImage(iChemical)->ComputeCellId(ijk);
                    this->StorePaintAction(iChemical,iCell,old_val);
                    this->GetImage(iChemical)->SetScalarComponentFromFloat(tx,ty,tz,0,val);
                }
            }
        }