from shapely.ops import nearest_points
selectisotemp[['nearpointA', 'nearpointB']] = selectisotemp.apply(
lambda x: nearest_points(x.mrband, mc['mrband']), axis=1, result_type='expand')
selectisotemp[['nearpointA', 'nearpointB']] = selectisotemp.apply(
lambda x: (x.nearpointA.coords[0], x.nearpointB.coords[0]), axis=1, result_type='expand')
selectisotemp[['bandA', 'bandB']] = selectisotemp.apply(
lambda x: (list(x.mrband.exterior.coords), list(mc['mrband'].exterior.coords)), axis=1,
result_type='expand')
selectisotemp['pointATRUE'] = selectisotemp.apply(lambda x: 1 if x.nearpointA in x.bandA else 0, axis=1)
selectisotemp['pointBTRUE'] = selectisotemp.apply(lambda x: 1 if x.nearpointB in x.bandB else 0, axis=1)
selectis