python使用shapely.ops nearest_points做最邻近分析-CSDN博客

本文链接：https://blog.csdn.net/cuisidong1997/article/details/140876853

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)

        selectisotemp['cibetween'] = selectisotemp['Ci'] + ';' + mc['Ci']
        selectisotemp[['A', 'B']] = selectisotemp.apply(lambda x: (x.Ci, mc['Ci']), axis=1,
                                                        result_type='expand')
        selectisotemp['cibetweensame'] = selectisotemp['cibetween'].apply(lambda x: ';'.join(sorted(x.split(';'))))
        print(selectisotemp)
        # selectisotemp.drop(['nearpointA', 'nearpointB', 'bandA', 'bandB'], axis=1, inplace=True)
        print(selectisotemp)
        selectiso.append(selectisotemp)