| 62 | } |
| 63 | |
| 64 | void VertexIsotropicOffsetParameter::process_roi() { |
| 65 | const size_t dim = m_wire_network->get_dim(); |
| 66 | MatrixFr vertices = m_wire_network->get_vertices(); |
| 67 | VectorF center = m_wire_network->center(); |
| 68 | |
| 69 | m_transforms.clear(); |
| 70 | IsotropicTransforms iso_trans(dim); |
| 71 | size_t roi_size = m_roi.size(); |
| 72 | size_t seed_vertex_index = m_roi.minCoeff(); |
| 73 | VectorF seed_dir = vertices.row(seed_vertex_index).transpose() - center; |
| 74 | |
| 75 | for (size_t i=0; i<roi_size; i++) { |
| 76 | VectorF v_dir = vertices.row(m_roi[i]).transpose() - center; |
| 77 | MatrixF trans = iso_trans.fit(seed_dir, v_dir); |
| 78 | m_transforms.push_back(trans); |
| 79 | } |
| 80 | } |