Method EPToctree

LASlib/src/lascopc.cpp:119–176

}

EPToctree::EPToctree(const LASheader& header)
{
  if (header.vlr_copc_info)
  {
    xmin = header.vlr_copc_info->center_x - header.vlr_copc_info->halfsize;
    ymin = header.vlr_copc_info->center_y - header.vlr_copc_info->halfsize;
    zmin = header.vlr_copc_info->center_z - header.vlr_copc_info->halfsize;
    xmax = header.vlr_copc_info->center_x + header.vlr_copc_info->halfsize;
    ymax = header.vlr_copc_info->center_y + header.vlr_copc_info->halfsize;
    zmax = header.vlr_copc_info->center_z + header.vlr_copc_info->halfsize;
    point_spacing = header.vlr_copc_info->spacing;
    max_depth = 0;
    grid_size = 0;
  }
  else
  {
    F64 center_x = (header.min_x + header.max_x)/2;
    F64 center_y = (header.min_y + header.max_y)/2;
    F64 center_z = (header.min_z + header.max_z)/2;
    F64 halfsize = MAX3(header.max_x - header.min_x, header.max_y - header.min_y, header.max_z - header.min_z)/2;

    xmin = center_x - halfsize;
    ymin = center_y - halfsize;
    zmin = center_z - halfsize;
    xmax = center_x + halfsize;
    ymax = center_y + halfsize;
    zmax = center_z + halfsize;
    point_spacing = 0;
    max_depth = 0;
    grid_size = 0;
  }

  if (header.vlr_copc_info && header.vlr_copc_entries)
  {
    registry.reserve(header.number_of_copc_entries);

    U64 ni = 0;
    U64 nf = 0;
    for (U32 i = 0 ; i < header.number_of_copc_entries ; i++)
    {
      if (header.vlr_copc_entries[i].point_count > 0)
      {
        ni = nf;
        nf = ni + header.vlr_copc_entries[i].point_count-1;
        EPToctant octant(header.vlr_copc_entries[i], xmin, ymin, zmin, xmax, ymax, zmax, ni, nf);
        registry[(EPTkey)octant] = octant;
        if (octant.d > max_depth) max_depth = octant.d;
        nf++;
      }
      else
      {
        // Octants with 0 point must be added to be able to reccurse the octree
        EPToctant octant(header.vlr_copc_entries[i], xmin, ymin, zmin, xmax, ymax, zmax, 0, 0);
        registry[(EPTkey)octant] = octant;
      }
    }
  }
}