core/vil/file_formats/vil_j2k_nitf2_pyramid_image_resource.cxx

Go to the documentation of this file.

#include "vil_j2k_nitf2_pyramid_image_resource.h"
//:
// \file
// Do not remove the following notice
// Approved for public Release, distribution unlimited
// DISTAR Case 14074

#include <vcl_cmath.h>

// By definition, each level is a factor of 2 reduced in scale
static float scale_at_level(unsigned level)
{
  if (level == 0)
    return 1.0f;
  float s = vcl_pow(2.0f, -static_cast<float>(level));
  return s;
}

vil_j2k_nitf2_pyramid_image_resource::
vil_j2k_nitf2_pyramid_image_resource(vil_image_resource_sptr const& nitf2): nitf2_sptr_(nitf2)
{
  ptr_ = 0;
  if (!nitf2_sptr_)
    return;
  vcl_string fmt = nitf2_sptr_->file_format();
  if (fmt=="nitf21")
    ptr_=static_cast<vil_nitf2_image*>(nitf2_sptr_.ptr());
}

unsigned vil_j2k_nitf2_pyramid_image_resource::nplanes()const
{
  unsigned ret = 0;
  if (ptr_)
    ret =  ptr_->nplanes();
  return ret;
}

//: The number of pixels in each row.
// Dimensions:  Planes x ni x nj.
// This method refers to the base (max resolution) image
unsigned vil_j2k_nitf2_pyramid_image_resource::ni() const
{
  unsigned ret = 0;
  if (ptr_)
    ret = ptr_->ni();
  return ret;
}

//: The number of pixels in each column.
// Dimensions:  Planes x ni x nj.
// This method refers to the base (max resolution) image
unsigned vil_j2k_nitf2_pyramid_image_resource::nj() const
{
  unsigned ret = 0;
  if (ptr_)
    ret =  ptr_->nj();
  return ret;
}

//: Pixel Format.
vil_pixel_format vil_j2k_nitf2_pyramid_image_resource::pixel_format() const
{
  if (ptr_)
    return  ptr_->pixel_format();
  return VIL_PIXEL_FORMAT_UNKNOWN;
}

//: Return a string describing the file format.
// Only file images have a format, others return 0
char const* vil_j2k_nitf2_pyramid_image_resource::file_format() const
{
  return "j2k_nitf2";
}


  // === Methods particular to pyramid resource ===

  //: Number of pyramid levels.
  //  Defined to assign 1000 pixels to the largest
  //  dimension of the least resolution level
unsigned vil_j2k_nitf2_pyramid_image_resource::nlevels() const
{
  if (!ptr_)
    return 0;
  unsigned max_dim = this->ni(), nj = this->nj();
  if (nj>max_dim)
    max_dim = nj;
  if (max_dim<1000)
    return 1;
  // Assume top level of the pyramid has maximum dimension of 1K pixels
  double scale = max_dim/1000.0;
  if (scale<=1.0)
    return 1;
  double lscale = vcl_log(scale);
  unsigned nlev = static_cast<unsigned>(lscale/vcl_log(2.0));
  return nlev;
}

  //: Get a partial view from the image from a specified pyramid level
vil_image_view_base_sptr
vil_j2k_nitf2_pyramid_image_resource::get_copy_view(unsigned i0, unsigned n_i,
                                                    unsigned j0, unsigned n_j,
                                                    unsigned level) const
{
  if (!ptr_||!(ptr_->is_jpeg_2000_compressed()))
     return 0;
  if (level>this->nlevels())
    level = this->nlevels()-1;
  double s = scale_at_level(level);
  double factor = static_cast<unsigned>(1/s);
  return ptr_->get_copy_view_decimated_j2k(i0, n_i, j0, n_j, factor, factor);
}

  //: Get a partial view from the image in the pyramid closest to scale.
  // The origin and size parameters are in the coordinate system of the base image.
  // The scale factor is with respect to the base image (base scale = 1.0).
vil_image_view_base_sptr
vil_j2k_nitf2_pyramid_image_resource::get_copy_view(unsigned i0, unsigned n_i,
                                                    unsigned j0, unsigned n_j,
                                                    const float scale,
                                                    float& actual_scale) const
{
  if (scale>=1.0f)
  {
    actual_scale = 1.0f;
    return this->get_copy_view(i0, n_i, j0, n_j, 0);
  }
  float f_lev = -vcl_log(scale)/vcl_log(2.0f);
  unsigned level = static_cast<unsigned>(f_lev);
  if (level>this->nlevels())
    level = this->nlevels()-1;
  actual_scale = scale_at_level(level);
  return this->get_copy_view(i0, n_i, j0, n_j, level);
}

//: Get an image resource from the pyramid at the specified level
vil_image_resource_sptr
vil_j2k_nitf2_pyramid_image_resource::get_resource(const unsigned level) const
{
  if (level==0)
    return nitf2_sptr_;
  return 0;
}

//: for debug purposes
void vil_j2k_nitf2_pyramid_image_resource::print(const unsigned level)
{
}