code: def get_area_def(self, key, info=None):
          
          # TODO Get projection information from input file
          a = 6378169.
          h = 35785831.
          b = 6356583.8
          lon_0 = 0.
          # area_extent = (-5432229.9317116784, -5429229.5285458621,
          #                5429229.5285458621, 5432229.9317116784)
          area_extent = (-5570248.4773392612, -5567248.074173444,
                         5567248.074173444, 5570248.4773392612)
          proj_dict = {'a': float(a),
                       'b': float(b),
                       'lon_0': float(lon_0),
                       'h': float(h),
                       'proj': 'geos',
                       'units': 'm'}
  
          area = geometry.AreaDefinition(
              'LI_area_name',
              "LI area",
              'geosli',
              proj_dict,
              self.ncols,
              self.nlines,
              area_extent)
          self.area = area
          logger.debug("Dataset area definition: \n {}".format(area))
          return area
docs: Create AreaDefinition for specified product.
  
          Projection information are hard coded for 0 degree geos projection
          Test dataset doesn't provide the values in the file container.
          Only fill values are inserted.