convert_to_rect2.py

def create_geo(LATRES, LONRES, ORIGIN_LAT, ORIGIN_LON, rowno, startcol, endcol):
  # represent each rectangle by a polygon of its corners
  top = ORIGIN_LAT - rowno * LATRES
  bot = ORIGIN_LAT - (rowno+1) * LATRES
  left = ORIGIN_LON + startcol * LONRES
  right = ORIGIN_LON + endcol * LONRES

  poly = json.dumps({
      'type': 'Polygon',
      'coordinates': [
          [ # polygon 1
            [left, top],  # topleft
            [left, bot],  # botleft
            [right, bot], # botright
            [right, top], # topright
            [left, top]   # same as first point
          ]
      ]
  })

  center = 'POINT({:2f} {:2f})'.format( (left+right)/2, (top+bot)/2 )
  return poly, center, endcol