simple.py

from shapely.geometry import Polygon, Point

def blockInZone(block, zone):
  boundaryPoints = []
  for pair in block:
    # reverse the order because from the file 
    #  coordinates data was read in lat, long
    #  and for Shapely we need to put the data
    #  as long, lat
    boundaryPoints.append((pair[1], pair[0]))
  blockPolygon = Polygon(boundaryPoints)

  boundaryPoints = []
  for pair in zone:
    boundaryPoints.append((pair[1], pair[0]))
  zonePolygon = Polygon(boundaryPoints)

  return blockPolygon.intersects(zonePolygon)