<tkinter>ボタンを押したとき最小二乗法の線の交点を求めたい

Question

前提・実現したいこと

　tkinter上に読み込んだcsvファイルを表示させ、グラフ上で、マウスの右クリックで最小二乗法したい範囲を選択し、左ダブルクリックで最小二乗法の線を引く。これを二回行い、二本の最小二乗法の線を引いた後、tkinter上に表示される"交点"というボタンを押すと二本の最小二乗法の交点の部分に点を打つというプログラムを作成したいと考えております。
tkinter上にグラフを表示させ、マウス操作を行い、ボタンを設置するところはできたのですが、ボタンを押しても交点に点を打てず、下記のようなエラーメッセージが出ます。このエラーメッセージの解決方法がわからず、悩んでおります。

発生している問題・エラーメッセージ

can't invoke "update" command: application has been destroyed

該当のソースコード

from tkinter import *
from matplotlib import pyplot as plt
from matplotlib.backends.backend_tkagg import FigureCanvasTkAgg
import numpy as np
from scipy import optimize
import pylab
import turtle
import tkinter as tk

x_data = []
y_data = []
katamuki = []
seppen = []
state = 0
f = open('xrr_data.csv', 'r')

datalist = f.readlines()
for data in datalist:
    if state == 1:
        x_data.append(float(data.split(',')[0]))
        y_data.append(float(data.split(',')[1]))
    if "Angle,Intensity" in data:
        state =1
    
f.close()
xdata = np.array(x_data)
ydata = np.array(y_data)


def drawlsm(lind,rind):
    
    global xdata2
    global ydata2
    
    xdata2 = xdata[lind:rind]
    ydata2 = ydata[lind:rind]
    
    parameter0 = [0.,0.]
    result = optimize.leastsq(fit_func,parameter0,args=(xdata2,ydata2))
    print(result)
    a_fit=result[0][0]
    b_fit=result[0][1]
    print(a_fit,b_fit)
    plt.plot(xdata,a_fit*xdata+b_fit,'k-', label='fitted line', linewidth=3, alpha=0.3)
    
    
#Least squares method with scipy.optimize
def fit_func(parameter,x,y):
    a = parameter[0]
    b = parameter[1]
    residual = y-(a*x+b)
    return residual

def oncmask(event):

    global stat
    global leftind, rightind
    
    ind=np.searchsorted(xdata,event.xdata)
    plt.title("You clicked index="+str(ind))
    if event.button==3 and stat==1:
        leftind=ind
        ax.plot([xdata[ind]],[ydata[ind]],".",color="red")
        stat=2
    elif event.button==3 and stat==2:
        rightind=ind
        ax.plot(xdata[leftind:rightind],ydata[leftind:rightind],color="red")
        stat=3
        print (leftind, rightind)
    elif event.button==1 and event.dblclick==1 and stat==3:
        plt.title("Approved")
        mask[leftind:rightind]=False
        drawlsm(leftind,rightind)
        stat=1
    elif event.button==2 and stat==3:  
        plt.title("Canceled")
        ax.plot(xdata[leftind:rightind],ydata[leftind:rightind],color="blue")
        ax.plot([xdata[leftind]],[ydata[leftind]],".",color="green")
        stat=1
    fig.canvas.draw()
    
katamuki2 = np.array(katamuki)
seppen2 = np.array(seppen)


def getintersection():#一番最後とその1つ前に引いた線の交点をget
    global xa #交点x座標
    indent = len(katamuki2)#傾きリストの要素数
    indent2= len(seppen2)#切片のリストの要素数
    print("傾きリストの要素数={}, 切片リストの要素数={}".format(indent,indent2))
    #xa求める
    xa=(seppen2[indent2-1]-seppen2[indent2-2])/(katamuki2[indent-2]-katamuki2[indent-1])#(b1-b0)/(a0-a1)
    print("交点座標={}".format(xa))
    ya=katamuki2[indent-1]*xa+seppen2[indent2-1]#x座標からy座標を求める
    plt.plot(xa,ya,".",color="m")#交点プロット
    plt.show()
    fig.canvas.draw()

#x=np.arange(0, 10, 0.1)
#y= np.random.rand(100) 

mask=np.ones(len(xdata),dtype=bool)

stat = 1
fig=plt.figure()
ax=fig.add_subplot(111)
ax.plot(xdata,ydata)

# Windowの設定
root = Tk()
root.title("Plot window")
root.geometry()

# Frameの設定
frame_1 = Frame(root, bd=4, relief=GROOVE)
frame_2 = Frame(root, bd=4, relief=GROOVE)

# widgetの設定
btn2 = Button(frame_1,text='交点',font=("",18),command=getintersection)
canvas = FigureCanvasTkAgg(fig, frame_2)#figをcanvasにのっける

# widgetの配置
frame_1.grid(row=0, column=0, sticky=W + E)
frame_2.grid(row=1, column=0)
btn2.pack(fill=X)
canvas.get_tk_widget().grid(row=0, column=0)

cid = fig.canvas.mpl_connect('button_press_event', oncmask)

canvas.start_event_loop()

試したこと

ボタンを押したときにtkinterを終了するプログラムを作成して試してみたところ、tkinterの終了はできましたが同じエラーメッセージが出ました。

補足情報（FW/ツールのバージョンなど）

Python初心者で、説明不足や至らない点があると思いますが、アドバイス頂けたら嬉しいです。
よろしくお願いします。

"追記"
excelデータはこのようになっております。

列1	列2
Angle	Intensity
0.027	879036
0.029	838848
0.031	766020
0.033	745008
0.035	714204
0.037	695844
0.039	677280
0.041	656268
0.043	628116
0.045	661164
0.047	649332
0.049	651984
0.051	640560
0.053	631788
0.055	647496
0.057	643620
0.059	656064
0.061	641988
0.063	668916
0.065	663000
0.067	671160
0.069	661164
0.071	680136
0.073	666876
0.075	691560
0.077	681564
0.079	712368
0.081	706656
0.083	744600
0.085	720732
0.087	740520
0.089	750108
0.091	745008
0.093	754800
0.095	776832
0.097	777240
0.099	778260
0.101	778260
0.103	795396
0.105	810492
0.107	812124
0.109	813960
0.111	812532
0.113	807636
0.115	838440

Comments

[teamikl]

コードの実行に必要な最低限のデータで良いので csv データの掲載をお願いします。

もしくは、テスト用に問題のエラーの再現が出来るデータを、
スクリプト内に埋め込んでコードを掲載して頂けると、
環境構築の手間が省けるので助かります。

実行はしてませんが、tkinter を一緒に使う場合は、
plt.show() の併用が問題が起こりそうなポイントです。

確認: エラーが出るのは「交点」のボタンから呼ばれる、
getintersection関数内ですか？エラーログに行数等は表示されてませんか？

[A_7]

ご指摘ありがとうございます。説明不足で申し訳ありません。
excelデータを質問の方に表で追記しました。目を通して頂けたら幸いです。
「交点」のボタンのgetintersection関数だけで実行してみたのですが、このようなエラーメッセージは出てきませんでした。
エラーメッセージの前にこのような表示がありました。
File "C:ファイル名", line 1305, in update
self.tk.call('update')

TclError: can't invoke "update" command: application has been destroyed

[teamikl]

ボタンを押したときには同様のエラーを再現できませんでした。
エラーの再現に、他に必要な操作はありますか？（グラフをクリックしてから交点？）

TypeError や IndexError がでます。（画面の枠外をクリックした時・交点ボタンをクリックした時）
また、ウィンドウを閉じてプログラムを終了したときに
「交点」のボタンを押していなくても 上記のTclError が起こります。

恐らく、３つくらいの問題があって、
報告されてるエラーは別の問題のもののようです。
終了時の問題については心当たりがあるので回答します。

[A_7]

TclErrorはボタンを押したときに起こるのではなく、tkinterのウィンドウを閉じた時に起こるものでした。
「交点」ボタンを押したときにはteamikl様の仰る通りIndexErrorが発生してしまいます。
Error動作の誤解や説明不足でお手数をお掛けしてしまい、申し訳ございません。
TclErrorについての回答をして頂き、本当にありがとうございます。

Answer 1

プログラム終了時の問題について

(※追記: 交点の問題に対する回答ではありません）

原因: 破棄(destroy)されたウィジェットに対して、
稼働中のイベントループがアクセスしようとしています。

対策: tkinter/matplotlib を合わせて使う場合は、
予めイベントループを止めてから destroy するようにします。

def _quit(): # 全てのウィンドウが閉じられた時に呼び出される
    canvas.stop_event_loop()
    root.destroy()
root.protocol("WM_DELETE_WINDOW", _quit)

因みに、canvas.start_event_loop ではなく root.mainloop を使う場合は、
canvas.stop_event_loop の代わりに root.quit と、
イベントループに対応した停止メソッドを使います。

---

通常の tkinter のプログラムでは、
ウィンドウを閉じるとプログラムの終了なので、
これらのメソッドはどちらもプログラムの終了として扱われる事が多いのですが、
イベントループが絡む場合使い分けが必要になることも有ります。

• root.destroy リソースの破棄
• root.quit イベントループの停止

Comments

[A_7]

回答をして頂き、ありがとうございます。
無事にTclErrorを解決する事ができました。
tkinterとmatplotlibを合わせて使用する時のイベントループの止め方がわからなかったので非常に勉強になりました。
本当にありがとうございました。