text: #!/usr/bin/env pytest
  ###############################################################################
  # $Id$
  #
  # Project:  GDAL/OGR Test Suite
  # Purpose:  Test CTG driver
  # Author:   Even Rouault, <even dot rouault at spatialys.com>
  #
  ###############################################################################
  # Copyright (c) 2011, Even Rouault <even dot rouault at spatialys.com>
  #
  # Permission is hereby granted, free of charge, to any person obtaining a
  # copy of this software and associated documentation files (the "Software"),
  # to deal in the Software without restriction, including without limitation
  # the rights to use, copy, modify, merge, publish, distribute, sublicense,
  # and/or sell copies of the Software, and to permit persons to whom the
  # Software is furnished to do so, subject to the following conditions:
  #
  # The above copyright notice and this permission notice shall be included
  # in all copies or substantial portions of the Software.
  #
  # THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
  # OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
  # FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
  # THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
  # LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
  # FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
  # DEALINGS IN THE SOFTWARE.
  ###############################################################################
  
  from osgeo import gdal
  
  
  import gdaltest
  
  ###############################################################################
  # Test a fake CTG dataset
  
  
  def test_ctg_1():
  
      tst = gdaltest.GDALTest('CTG', 'ctg/fake_grid_cell', 1, 21)
      expected_gt = [421000.0, 200.0, 0.0, 5094400.0, 0.0, -200.0]
      expected_srs = """PROJCS["WGS 84 / UTM zone 14N",
      GEOGCS["WGS 84",
          DATUM["WGS_1984",
              SPHEROID["WGS 84",6378137,298.257223563,
                  AUTHORITY["EPSG","7030"]],
              AUTHORITY["EPSG","6326"]],
          PRIMEM["Greenwich",0,
              AUTHORITY["EPSG","8901"]],
          UNIT["degree",0.0174532925199433,
              AUTHORITY["EPSG","9122"]],
          AUTHORITY["EPSG","4326"]],
      UNIT["metre",1,
          AUTHORITY["EPSG","9001"]],
      PROJECTION["Transverse_Mercator"],
      PARAMETER["latitude_of_origin",0],
      PARAMETER["central_meridian",-99],
      PARAMETER["scale_factor",0.9996],
      PARAMETER["false_easting",500000],
      PARAMETER["false_northing",0],
      AUTHORITY["EPSG","32614"],
      AXIS["Easting",EAST],
      AXIS["Northing",NORTH]]"""
      ret = tst.testOpen(check_gt=expected_gt, check_prj=expected_srs)
  
      if ret == 'success':
          ds = gdal.Open('data/ctg/fake_grid_cell')
          lst = ds.GetRasterBand(1).GetCategoryNames()
          assert lst is not None and lst, 'expected non empty category names for band 1'
          lst = ds.GetRasterBand(2).GetCategoryNames()
          assert lst is None, 'expected empty category names for band 2'
          assert ds.GetRasterBand(1).GetNoDataValue() == 0, \
              'did not get expected nodata value'
  
      return ret
  
  
  
  