code: def gt_corners(gt, nx, ny):
      
      ul = [gt[0], gt[3]]
      ll = [gt[0], gt[3] + (gt[5] * ny)]
      ur = [gt[0] + (gt[1] * nx), gt[3]]
      lr = [gt[0] + (gt[1] * nx), gt[3] + (gt[5] * ny)]
      return ul, ll, ur, lr
docs: Get corner coordinates based on input geotransform and raster dimensions