v@size = getbbox_size(0);

i[]@pts = primpoints(0, @primnum);
i@cluster = -1;
int remove = 1;
for(int i=0; i < len(@pts); i++){
    int pt = @pts[i];
    i@tc = point(0, "cluster", pt);
    
    if(@cluster == -1){   
    @cluster = @tc;
    }else if(@cluster != @tc){
        remove = 0;
        @Cd = set(0, 0,0);
    }
}
if(remove == 1){
    removeprim(0, @primnum, 1);
}

int neighs[] = polyneighbours(0, @primnum);
int pt = addpoint(0, @P);

for(int i=0; i<len(neighs); i++){
    int neigh = neighs[i];
    vector pos = prim(0, "P", neigh);
    int npt = addpoint(0, pos);
    int line = addprim(0, "polyline", pt, npt);
}

removeprim(0, @primnum, 1);