docs: def _load_GeoTransform(self):
          """Calculate latitude and longitude variable calculated from the
          gdal.Open.GetGeoTransform method"""
          def load_lon():
              return arange(ds.RasterXSize)*b[1]+b[0]
  
          def load_lat():
              return arange(ds.RasterYSize)*b[5]+b[3]
          ds = self.ds
          b = self.ds.GetGeoTransform()  # bbox, interval
          if with_dask:
              lat = Array(
                  {('lat', 0): (load_lat,)}, 'lat', (self.ds.RasterYSize,),
                  shape=(self.ds.RasterYSize,), dtype=float)
              lon = Array(
                  {('lon', 0): (load_lon,)}, 'lon', (self.ds.RasterXSize,),
                  shape=(self.ds.RasterXSize,), dtype=float)
          else:
              lat = load_lat()
              lon = load_lon()
          return Variable(('lat',), lat), Variable(('lon',), lon)