#-*- coding:utf-8 -*-
# 计算点到标准轮廓距离的极性

import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
import cv2

data = pd.read_csv("D:\\result.csv")
standard = pd.read_csv("D:\\standard.csv")
dataArray = np.array(data)[:,0:2:1]*1000
standardArray = np.array(standard)*1000

polarity = np.array([0])

for index in range(0,dataArray.shape[0],1):
    temp = cv2.pointPolygonTest(standardArray.astype(int),tuple(dataArray[index,:]),False)
    polarity =np.r_[polarity,temp]

polarity = np.delete(polarity.reshape(1,-1).transpose(),0,0)
polarityDistance = np.array(data)[:,2] * polarity.reshape(1,-1)

newData = np.array(data)[:,0:2:1]
newData = np.c_[newData,polarityDistance.transpose()]

polarityData = pd.DataFrame(newData).sort_values(by=0)
polarityData.to_csv("D:\\polarity.csv",header = False, index = False)