Dear Houdini - create points at the center of each primitive

float bound[] = primintrinsic(@OpInput1, "bounds", @primnum);

vector min = set(bound[0], bound[2], bound[4]);
vector max = set(bound[1], bound[3], bound[5]);

addpoint(geoself(),avg(min,max));
removeprim(geoself(),@primnum,1);