CN120446830B - 一种针对医用回旋加速器磁铁的磁场分布探测方法 - Google Patents

Abstract

本申请涉及磁场分布探测技术领域，具体涉及一种针对医用回旋加速器磁铁的磁场分布探测方法，该方法包括：针对待探测的医用回旋加速器磁铁，将其均分为各磁场区域；采集各磁场区域内各时刻的磁场强度，将所有磁场区域划分为各类别；确定当前时刻任一磁场区域的磁场电流受扰值；获取当前时刻所述任一磁场区域的线圈电流受扰度；结合各磁场区域的所述磁场电流受扰值与所述线圈电流受扰度，确定各磁场区域所测磁场强度的测量误差度，对各磁场区域的磁场强度数据进行判断筛选。本申请提高了医用回旋加速器磁铁的磁场分布探测的精度。

Description

一种针对医用回旋加速器磁铁的磁场分布探测方法

技术领域

本申请涉及磁场分布探测技术领域，具体涉及一种针对医用回旋加速器磁铁的磁场分布探测方法。

背景技术

医用回旋加速器是用于生产PET/CT检查当中放射性药物的重要设备，保障其正常运行至关重要。医用回旋加速器中，磁场的均匀性及稳定性直接影响粒子的加速效率及聚焦度，磁场测量可以及时发现并纠正磁场中的偏差和波动，如因磁极表面缺陷及外界干扰对磁场的影响，以确保加速器稳定高效地运行，提高粒子的加速效率和束流质量；故高精度的磁场测量至关重要。

现有技术对于磁场分布的探测主要利用霍尔探头来扫描磁场，使探头移动到各类磁铁所需测试位置处，获取任意空间磁场分布情况，而传统探测方法往往专注于判断磁铁本身的缺陷所带来的磁场影响，容易忽略磁铁内外部电流对监测的影响，若测量的电磁铁磁场受到电流的影响较大，在磁场均匀性测试时，则会带来较大的偏差，影响磁场探测精度。

发明内容

为了解决上述技术问题，本申请提供一种针对医用回旋加速器磁铁的磁场分布探测方法，以解决现有的问题。

本申请的一种针对医用回旋加速器磁铁的磁场分布探测方法采用如下技术方案：

本申请一个实施例提供了一种针对医用回旋加速器磁铁的磁场分布探测方法，该方法包括以下步骤：

针对待探测的医用回旋加速器磁铁，将其均分为各磁场区域；采集各磁场区域内各时刻的磁场强度，将各磁场区域当前时刻前所有时刻的磁场强度组成磁场强度序列；

利用各磁场区域的磁场强度序列中峰值的所占比例，确定各磁场区域的波动频率；基于任意两磁场区域之间磁场强度序列的相似性，将所有磁场区域划分为各类别；

根据所述划分的类别数量和任一磁场区域所属类别中所有磁场区域之间的度量距离，结合所述任一磁场区域与除其所属类别中所有磁场区域之外，剩余各磁场区域之间的所述波动频率的差异，以及磁场强度序列之间峰值分布的差异，确定当前时刻所述任一磁场区域的磁场电流受扰值；

分析当前时刻所述任一磁场区域的径向方向上所有磁场区域与除所述任一磁场区域外其径向方向上剩余所有磁场区域的磁场强度的趋势性强度的差异，结合所述任一磁场区域与其余各磁场区域的磁场强度序列中相同位序的峰值的时间差异，确定当前时刻所述任一磁场区域的线圈电流受扰度；

结合各磁场区域的所述磁场电流受扰值与所述线圈电流受扰度，确定各磁场区域所测磁场强度的测量误差度，对各磁场区域的磁场强度数据进行判断筛选。

在其中一种实施例中，所述波动频率为各磁场区域的磁场强度序列中峰值的总个数与采样时刻的总个数的比值。

在其中一种实施例中，所述将所有磁场区域划分为各类别，包括：

将所述相似性大于等于预设阈值的两磁场区域作为同类别的磁场区域。

在其中一种实施例中，所述确定当前时刻所述任一磁场区域的磁场电流受扰值，包括：

计算所述任一磁场区域的所属类别中所有任意两磁场区域之间的度量距离的均值，计算所述均值与所述划分的类别数量的乘积；

将除所述任一磁场区域的所属类别中所有磁场区域之外，所述任一磁场区域与剩余各磁场区域的波动频率的差异，记为第一差异；

计算各磁场区域的磁场强度序列中所有峰值的平均值，将除所述任一磁场区域的所属类别中所有磁场区域之外，所述任一磁场区域与剩余各磁场区域的所述平均值的差异，记为第二差异；

融合所述第一差异与所述第二差异，结合所述乘积，确定所述磁场电流受扰值。

在其中一种实施例中，所述磁场电流受扰值与所述融合的结果成正相关，与所述乘积成负相关。

在其中一种实施例中，所述磁场电流受扰值的进一步确定包括：

计算所述第一差异与所述第二差异的和值，计算所述任一磁场区域与剩余所有磁场区域的所述和值的累加和，所述磁场电流受扰值为所述累加和与所述乘积的比值。

在其中一种实施例中，所述确定当前时刻所述任一磁场区域的线圈电流受扰度，包括：

针对所述任一磁场区域，获取磁铁中心指向所述任一磁场区域并延长至磁铁边缘的径向方向，将径向方向上当前时刻所有磁场区域的磁场强度组成磁场径向序列，将磁场径向序列中剔除所述任一磁场区域的磁场强度后的序列，记为磁场径向序列；

计算磁场径向序列与磁场径向序列的趋势性强度的差异，记为第三差异，所述线圈电流受扰度与所述第三差异成正相关，与所述时间差异成负相关。

在其中一种实施例中，所述线圈电流受扰度的表达式为：

；式中，为当前时刻第i个磁场区域的线圈电流受扰度，为当前时刻第i个磁场区域的磁场径向序列的趋势性强度，为当前时刻第i个磁场区域的磁场径向序列的趋势性强度，为当前时刻第i个磁场区域的磁场强度序列中第k个峰值的对应时刻，为当前时刻第n个磁场区域的磁场强度序列中第k个峰值的对应时刻，N为医用回旋加速器磁铁划分的磁场区域的数量，为当前时刻第i个磁场区域的磁场强度序列中峰值的数量，为预设大于0的数值。

在其中一种实施例中，所述测量误差度为所述磁场电流受扰值与所述线圈电流受扰度的乘积的归一化值。

在其中一种实施例中，所述对各磁场区域的磁场强度数据进行判断筛选，包括：

获取当前所有磁场区域的测量误差度的分割阈值，若所述任一磁场区域的测量误差度大于等于所述分割阈值，判定所述任一磁场区域的磁场强度测量存在偏差，将该磁场强度进行剔除，否则，判定所述任一磁场区域的磁场强度测量精准，将该磁场强度进行保留。

本申请至少具有如下有益效果：

本申请针对待探测的医用回旋加速器磁铁，将其均分为各磁场区域；采集各磁场区域内各时刻的磁场强度，将各磁场区域当前时刻前所有时刻的磁场强度组成磁场强度序列；提高了局部磁场动态演化的捕捉能力，增强了异常磁场区域的早期识别能力；基于任意两磁场区域之间磁场强度序列的相似性，将所有磁场区域划分为各类别；优化了磁场分布的空间关联建模，避免了孤立分析忽视磁场协同畸变现象；进一步，确定当前时刻各磁场区域的磁场电流受扰值，提升了电流扰动溯源精度，量化了磁场强度受到的电流干扰程度，提升了磁场强度测量精度分析的准确性与可靠性；确定当前时刻各磁场区域的线圈电流受扰度，实现了径向磁场失衡的精准定位，增强了磁场强度时域-空域的耦合解析能力，反映了磁场强度受到的线圈电流的干扰程度；结合各磁场区域的所述磁场电流受扰值与所述线圈电流受扰度，确定各磁场区域所测磁场强度的测量误差度，对各磁场区域的磁场强度数据进行判断筛选，提高了各磁场区域的磁场强度的异常识别能力，提升了磁场强度测量的可信度，并且提高了磁场分布探测的精度。

附图说明

为了更清楚地说明本申请实施例或现有技术中的技术方案和优点，下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其它附图。

图1为本申请提供的一种针对医用回旋加速器磁铁的磁场分布探测方法的步骤流程图；

图2为磁场区域所测磁场强度的测量误差度确定流程图。

具体实施方式

为了更进一步阐述本申请为达成预定发明目的所采取的技术手段及功效，以下结合附图及较佳实施例，对依据本申请提出的一种针对医用回旋加速器磁铁的磁场分布探测方法，其具体实施方式、结构、特征及其功效，详细说明如下。在下述说明中，不同的“一个实施例”或“另一个实施例”指的不一定是同一实施例。此外，一或多个实施例中的特定特征、结构或特点可由任何合适形式组合。

除非另有定义，本文所使用的所有的技术和科学术语与属于本申请的技术领域的技术人员通常理解的含义相同。

下面结合附图具体的说明本申请所提供的一种针对医用回旋加速器磁铁的磁场分布探测方法的具体方案。

本申请一个实施例提供的一种针对医用回旋加速器磁铁的磁场分布探测方法，具体的，提供了如下的一种针对医用回旋加速器磁铁的磁场分布探测方法，请参阅图1，该方法包括以下步骤：

步骤S001，针对待探测的医用回旋加速器磁铁，将其均分为各磁场区域；采集各磁场区域内各时刻的磁场强度，将各磁场区域当前时刻前所有时刻的磁场强度组成磁场强度序列。

针对待探测的医用回旋加速器磁铁，本实施例将磁铁分为相等大小的N个磁场区域，并利用霍尔探头，每间隔时间T对各个磁场区域的中心点进行磁场强度数据的测量，最后根据磁场强度的采样时刻顺序构建医用回旋加速器磁铁各个磁场区域当前时刻的磁场强度序列，即将各磁场区域当前时刻前所有时刻的磁场强度按时间先后顺序组成各磁场区域的磁场强度序列。

本实施例中，N=64，T=0.5s，实施者可根据实际医用回旋加速器的型号规模对N和T的取值自行调整，本实施例对此不做限制。

步骤S002，利用各磁场区域的磁场强度序列中峰值的所占比例，确定各磁场区域的波动频率；基于任意两磁场区域之间磁场强度序列的相似性，将所有磁场区域划分为各类别。

通常情况下，医用回旋加速器为径向扇型等时性回旋加速器，具体的，这种回旋加速器的磁场分布沿径向和方位角方向变化，形成“峰区”和“谷区”。在不受干扰的情况下，磁场强度沿径向方向随着半径的增大而增大，与磁铁中心距离相近的磁铁区域磁场强度越接近，整体的分布较为均匀。

在实际应用时，医用回旋加速器的磁铁在浇铸过程中，盖板内部可能会出现缩孔、夹杂等缺陷，这些缺陷会影响材料的磁性能，进而影响磁场的分布，造成局部区域的磁场强度均匀度降低，与其余磁场区域产生一定的磁场强度差异，且差异程度不均匀，磁场的整体分布规律发生改变。具体表现为，本实施例采集到的磁场强度数据从磁铁中心至磁铁外侧边缘的分布规律发生改变，磁场强度数据与其所在区域位置的相关程度下降，存在部分磁场区域与其余磁场区域的磁场强度差异较大，且分布较为离散。

而在回旋加速器运行过程中，射频系统、束流提取系统的电流变化也可能会导致回旋加速器的磁场受到干扰，同样会导致磁场的整体分布规律发生变化，但此类干扰的范围较为集中，且由于干扰来源一致，因而，磁场强度受到干扰的程度和频率相似，如束流提取系统的影响区域通常更接近回旋加速器的离子出口。具体表现为，单个磁场区域的磁场强度发生波动情况的频率与其距离较远的磁场区域相比差异较大，但与其周围的磁场区域发生波动情况相近，磁场强度变化情况近似的磁场区域分布更为密集。

基于上述分析，针对当前时刻的各磁场区域，本实施例采用峰值查找算法获取磁场强度序列中的各峰值及其对应时刻，将各磁场区域的磁场强度序列中峰值的总个数与采样时刻的总个数的比值，作为当前时刻各磁场区域的波动频率。

进一步地，计算当前时刻医用回旋加速器磁铁划分的N个磁场区域中，任意两磁场区域之间的磁场强度序列的相似性，相似性越大，表明对应的两个磁场区域的磁场强度分布近似程度越高，因此，本实施例预设阈值Q，将所述相似性大于等于阈值Q的两个磁场区域作为相同类别的磁场区域，基于此，可将N个磁场区域划分为多个类别。

需要说明的是，若存在磁场区域与剩余所有磁场区域的所述相似性均小于阈值Q，则将该磁场区域单独作为一个类别；本实施例中所述相似性采用余弦相似度进行计算，实施者可自行选择现有其他可行的相似性计算方法，例如皮尔逊相关系数等；本实施例中阈值Q设定为0.97，实施者可根据实际情况自行设定，本实施例对此不做限制。

步骤S003，根据所述划分的类别数量和任一磁场区域所属类别中所有磁场区域之间的度量距离，结合所述任一磁场区域与除其所属类别中所有磁场区域之外，剩余各磁场区域之间的所述波动频率的差异，以及磁场强度序列之间峰值分布的差异，确定当前时刻所述任一磁场区域的磁场电流受扰值。

基于上述分析，本实施例计算当前时刻各磁场区域的磁场电流受扰值，具体为：

计算任一磁场区域的所属类别中所有任意两磁场区域之间的度量距离的均值，计算所述均值与所述划分的类别数量的乘积；

将除所述任一磁场区域的所属类别中所有磁场区域之外，所述任一磁场区域与剩余各磁场区域的波动频率的差异，记为第一差异；

计算各磁场区域的磁场强度序列中所有峰值的平均值，将除所述任一磁场区域的所属类别中所有磁场区域之外，所述任一磁场区域与剩余各磁场区域的所述平均值的差异，记为第二差异；

融合所述第一差异与所述第二差异，结合所述乘积，确定所述磁场电流受扰值。

需要说明的是，差异表示两个变量之间的差别程度，具体可以采用差值绝对值、差值的平方、比值等方式进行计算，本实施例对此不做限制；融合表示将多个变量进行结合，具体可以采用相加、相乘、加乘混合等方式进行计算。

本实施例中所述磁场电流受扰值的计算方式为：

；式中，为当前时刻第i个磁场区域的磁场电流受扰值，为医用回旋加速器磁铁的所有磁场区域划分的类别数量，为当前时刻第i个磁场区域的所属类别中所有任意两磁场区域之间的度量距离的均值，分别为第i个、第j个磁场区域的波动频率，分别为第i个、第j个磁场区域的磁场强度序列中所有峰值的均值，J为除第i个磁场区域所属类别的所有磁场区域以外，剩余所有磁场区域的数量。将记为第一差异，记为第二差异。

需要说明的是，本实施例中所述度量距离采用欧氏距离进行计算，实施者可自行选择现有其他可行的度量距离计算方式；磁场区域之间的欧氏距离计算方式为，获取所有磁场区域的中心点，计算两磁场区域中心点之间的欧氏距离作为磁场区域之间的欧氏距离。当第i个磁场区域所属类别中仅存在第i个磁场区域时，令。

应当理解的是，当第i个磁场区域所属类别的所有磁场区域之间欧氏距离越小，与其余类别中的磁场区域的波动频率、峰值均值差异越大，且整体磁场区域进行分类后的类别越少，则说明第i个磁场区域所测量出的磁场强度越可能受到了电流的干扰，测量的准确程度较低。

步骤S004，分析当前时刻所述任一磁场区域的径向方向所有磁场区域与除所述任一磁场区域外其径向方向上剩余所有磁场区域的磁场强度的趋势性强度的差异，结合所述任一磁场区域与其余各磁场区域的磁场强度序列中相同位序的峰值的时间差异，确定当前时刻所述任一磁场区域的线圈电流受扰度。

医用回旋加速器在运行过程中，磁场主要通过加速器的主线圈产生，其电流的稳定性直接影响磁场整体的均匀性和强度，当主线圈的电流出现不稳定或波动时，则会导致磁场强度的变化，具体表现为整体磁场强度都会受到较为近似的干扰，进而若单个磁场区域为磁铁本身缺陷的影响，磁场电流受扰值也较大，进而容易将该磁场区域判断为受电流干扰的区域，从而误判该磁场区域的真实情况，因此对磁场区域的磁场强度进一步分析。

当磁场受到主线圈的电流不稳定波动影响时，磁场整体都会受到相应的干扰，从而所有磁场区域都会受到近似的干扰变化，且出现干扰的时刻也近似，而相较于存在磁铁缺陷的区域，受到主线圈电流影响的区域磁场强度虽然会整体发生改变，但由于受影响程度近似，因而磁场分布的规律受影响程度较低，即磁场强度沿径向方向随着半径的增大而增大，相反，存在磁铁缺陷的区域由于其本身就对磁场产生了一定的影响，因而，其所在径向方向上的磁场分布规律则存在较大影响。

为表征磁铁中磁场强度分布规律，针对单个磁场区域，以第i个磁场区域为例，由磁铁中心指向第i个磁场区域的方向，将磁铁中心与第i个磁场区域中心相连，并将连接线延长至磁铁边缘，记为第i个磁场区域的径向方向，将当前时刻第i个磁场区域的径向方向上所有磁场区域的磁场强度，以距离磁铁中心点由近至远的顺序，组成磁场径向序列，进一步地，将磁场径向序列中第i个磁场区域的磁场强度数据剔除后，重新构建序列，记为磁场径向序列。

结合上述分析，本实施例计算当前时刻各磁场区域的线圈电流受扰度，具体表达式为：

；式中，为当前时刻第i个磁场区域的线圈电流受扰度，为当前时刻第i个磁场区域的磁场径向序列的趋势性强度，为当前时刻第i个磁场区域的磁场径向序列的趋势性强度，为当前时刻第i个磁场区域的磁场强度序列中第k个峰值的对应时刻，为当前时刻第n个磁场区域的磁场强度序列中第k个峰值的对应时刻，N为医用回旋加速器磁铁划分的磁场区域的数量，为当前时刻第i个磁场区域的磁场强度序列中峰值的数量，为预设大于0的数值，避免分母为0，本实施例中，实施者可根据实际情况自行设定，本实施例对此不做限制。将记为第三差异。本实施例中磁场径向序列的趋势性强度与磁场径向序列的趋势性强度，采用STL（Seasonal andTrend decomposition using Loess）分解算法中的计算方式进行获取，STL分解算法为现有公知技术，具体过程不做赘述。

应当理解的是，当第i个磁场区域的磁场径向序列的趋势性强度越强，且磁场径向序列与磁场径向序列的趋势性强度差异越大，以及第i个磁场区域与其余磁场区域的磁场强度序列中出现波峰的采样时间越接近，则说明第i个磁场区域中采集到的磁场强度数据越可能受到了线圈电流不稳定波动的影响，测量值准确度越低。

步骤S005，结合各磁场区域的所述磁场电流受扰值与所述线圈电流受扰度，确定各磁场区域所测磁场强度的测量误差度，对各磁场区域的磁场强度数据进行判断筛选。

进一步，本实施例融合各磁场区域的所述磁场电流受扰值与所述线圈电流受扰度，确定各磁场区域所测磁场强度的测量误差度，表征单个磁场区域的所测磁场强度值存在误差的可能性，表达式为：

；式中，为当前时刻第i个磁场区域所测磁场强度的测量误差度，为当前时刻第i个磁场区域的磁场电流受扰值，为当前时刻第i个磁场区域的线圈电流受扰度，norm（）为归一化函数。磁场区域所测磁场强度的测量误差度确定流程图如图2所示。

当第i个磁场区域所测得的磁场强度数据受到磁场的电流干扰以及线圈电流干扰程度越高，则说明该磁场强度数据存在偏差的可能性越强，越不应该用于判断医用回旋加速器的磁场实际情况。

通过上述方式，可计算出当前所有磁场区域的测量误差度，将所有测量误差度作为输入，利用交叉验证的方式，获得测量误差度的分割阈值，当磁场区域的测量误差度大于等于分割阈值时，则说明该磁场区域的测量值越可能受到了回旋加速器内外部的电流干扰，因此该测量值存在较大的偏差，反之，则判定该磁场区域的测量值为精准，可用于判断磁场的分布均匀程度。

基于此，可判断每次磁场强度测量时，各个磁场区域的测量误差情况，当出现某磁场区域的测量误差度大于等于分割阈值时，则对该磁场区域的磁场强度数据进行剔除，重置测量工况后，重新进行磁场强度测量，直至医用回旋加速器磁铁的所有磁场区域测量误差度小于分割阈值，实现医用回旋加速器磁铁的磁场分布探测。

需要说明的是：上述本申请实施例先后顺序仅仅为了描述，不代表实施例的优劣。且上述对本说明书特定实施例进行了描述。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。

本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同或相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。

以上所述实施例仅用以说明本申请的技术方案，而非对其限制；对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换，并不使相应技术方案的本质脱离本申请各实施例技术方案的范围，均应包含在本申请的保护范围之内。

Claims (8)

1.一种针对医用回旋加速器磁铁的磁场分布探测方法，其特征在于，该方法包括以下步骤：

针对待探测的医用回旋加速器磁铁，将其均分为各磁场区域；采集各磁场区域内各时刻的磁场强度，将各磁场区域当前时刻前所有时刻的磁场强度组成磁场强度序列；

利用各磁场区域的磁场强度序列中峰值的所占比例，确定各磁场区域的波动频率；基于任意两磁场区域之间磁场强度序列的相似性，将所有磁场区域划分为各类别；

计算任一磁场区域的所属类别中所有任意两磁场区域之间的度量距离的均值，计算所述均值与划分的类别数量的乘积；将除所述任一磁场区域的所属类别中所有磁场区域之外，所述任一磁场区域与剩余各磁场区域的波动频率的差异，记为第一差异；计算各磁场区域的磁场强度序列中所有峰值的平均值，将除所述任一磁场区域的所属类别中所有磁场区域之外，所述任一磁场区域与剩余各磁场区域的所述平均值的差异，记为第二差异；融合所述第一差异与所述第二差异，结合所述乘积，确定当前时刻所述任一磁场区域的磁场电流受扰值；

分析当前时刻所述任一磁场区域的径向方向上所有磁场区域与除所述任一磁场区域外其径向方向上剩余所有磁场区域的磁场强度的趋势性强度的差异，结合所述任一磁场区域与其余各磁场区域的磁场强度序列中相同位序的峰值的时间差异，确定当前时刻所述任一磁场区域的线圈电流受扰度；

结合各磁场区域的所述磁场电流受扰值与所述线圈电流受扰度，确定各磁场区域所测磁场强度的测量误差度，对各磁场区域的磁场强度数据进行判断筛选；

针对所述任一磁场区域，获取磁铁中心指向所述任一磁场区域并延长至磁铁边缘的径向方向，将径向方向上当前时刻所有磁场区域的磁场强度组成磁场径向序列，将磁场径向序列中剔除所述任一磁场区域的磁场强度后的序列，记为磁场径向序列；

计算磁场径向序列与磁场径向序列的趋势性强度的差异，记为第三差异，所述线圈电流受扰度与所述第三差异成正相关，与所述时间差异成负相关。

2.如权利要求1所述的一种针对医用回旋加速器磁铁的磁场分布探测方法，其特征在于，所述波动频率为各磁场区域的磁场强度序列中峰值的总个数与采样时刻的总个数的比值。

3.如权利要求1所述的一种针对医用回旋加速器磁铁的磁场分布探测方法，其特征在于，所述将所有磁场区域划分为各类别，包括：

将所述相似性大于等于预设阈值的两磁场区域作为同类别的磁场区域。

4.如权利要求1所述的一种针对医用回旋加速器磁铁的磁场分布探测方法，其特征在于，所述磁场电流受扰值与所述第一差异、所述第二差异均成正相关，与所述乘积成负相关。

5.如权利要求1所述的一种针对医用回旋加速器磁铁的磁场分布探测方法，其特征在于，所述磁场电流受扰值的进一步确定包括：

计算所述第一差异与所述第二差异的和值，计算所述任一磁场区域与剩余所有磁场区域的所述和值的累加和，所述磁场电流受扰值为所述累加和与所述乘积的比值。

6.如权利要求1所述的一种针对医用回旋加速器磁铁的磁场分布探测方法，其特征在于，所述线圈电流受扰度的表达式为：

；式中，为当前时刻第i个磁场区域的线圈电流受扰度，为当前时刻第i个磁场区域的磁场径向序列的趋势性强度，为当前时刻第i个磁场区域的磁场径向序列的趋势性强度，为当前时刻第i个磁场区域的磁场强度序列中第k个峰值的对应时刻，为当前时刻第n个磁场区域的磁场强度序列中第k个峰值的对应时刻，N为医用回旋加速器磁铁划分的磁场区域的数量，为当前时刻第i个磁场区域的磁场强度序列中峰值的数量，为预设大于0的数值。

7.如权利要求1所述的一种针对医用回旋加速器磁铁的磁场分布探测方法，其特征在于，所述测量误差度为所述磁场电流受扰值与所述线圈电流受扰度的乘积的归一化值。

8.如权利要求1所述的一种针对医用回旋加速器磁铁的磁场分布探测方法，其特征在于，所述对各磁场区域的磁场强度数据进行判断筛选，包括：

获取当前所有磁场区域的测量误差度的分割阈值，若所述任一磁场区域的测量误差度大于等于所述分割阈值，判定所述任一磁场区域的磁场强度测量存在偏差，将该磁场强度进行剔除，否则，判定所述任一磁场区域的磁场强度测量精准，将该磁场强度进行保留。