3*3のマス目代入方法について

前提・実現したいこと

センサーの座標（X[,]）とその取得温度(t[])、観測領域の広さ（a,b）が分かっています。そこで、ある地点(x,y)を指定した際に、その地点を中心として、33のマス目を形成します。各マス目の大きさはNとし、その各マス目(dot[][])に含まれるセンサーの温度を代入していきたいです。
・マス内に複数のセンサーが含まれた場合、温度の合計を算出
・マス内にセンサーが無い場合、０を代入
・中心のマス目に他のセンサーが含まれる場合、値は代入しない（＝中心セルの値は０）
・33のマス目を形成した際に、観測領域の範囲を超えてしまった場合、エラーを表示
・各マス目に含まれるセンサー番号（配列の要素番号）を保持

以上を実現したく、苦戦しております。考え方だけでも助かりますので、ご教授頂けると幸いです。

該当のソースコード

python
1#X：センサーの座標
2X = np.array( [
3    [21.5,23.0],
4    [24.5,20.0],
5    [19.5,19.0],
6    [22.5,15.0],
7    [24.5,12.0],
8    [19.5,12.0],
9    [22.5, 8.0],
10    [24.5, 4.0],
11    [21.5, 2.0],
12    [19.5, 5.0],
13    [16.5, 3.0],
14    [13.5, 1.0],
15    [12.5, 5.0],
16    [ 8.5, 6.0],
17    [ 5.5, 3.0],
18    [ 1.5, 2.0],
19    [ 1.5, 8.0],
20    ])
21
22#t：各センサーの温度
23t=np.array([
24　19.9884 ,
25　19.567 , 
26　20.204 , 
27　19.8806 ,
28　20.4196 ,
29　19.4298 ,
30　19.1848 ,
31　19.0672 ,
32　19.3612 ,
33　18.6556 ,
34　17.803 , 
35　17.8422 ,
36　17.852 , 
37　23.0754 ,
38　17.2444 ,
39　16.9308 ,
40　19.2534 ,
41])
42
43#観測領域の広さa*b
44a=40
45b=35
46
47#各マス目の辺の長さN*N
48N=3
49
50#中心地点(指定地点)
51x= 10
52y= 10
53
54"""
55処理
56"""
57
58