一种厚度测量的中值滤波方法技术

本发明专利技术属于厚度测量技术领域，公开了一种厚度测量的中值滤波方法，包括，步骤S1：首先通过测距装置对被测物体的表面进行N次预测量，获取N个测量值si；步骤S2：创建一个容量为n的先进先出队列来存储测量值，当队列里的测量值数量达到n时，队列丢弃最早进入队列的测量值，并将新的测量值放入队列；步骤S3：通过快速中值滤波的方式获取队列中的中值a；步骤S4：预测量测距装置使用中的背景距离b，被测物体的厚度值L即为b‑a。本发明专利技术可以很好地避免偶尔出现的测量误差，选择的中值即可作为准确测量值来计算被测物体的厚度，来发出合适的工作管电压kvp。

A Median Filtering Method for Thickness Measurement

The invention belongs to the technical field of thickness measurement, and discloses a median filtering method for thickness measurement, which includes: firstly, N times of prediction is made on the surface of the object to be measured by a distance measuring device to obtain N measurement values si; secondly, an advanced first-out queue with a capacity of n is created to store measurement values, and when the number of measurement values in the queue reaches n, the queue discards the earliest entry. The measured value of the queue is put into the queue, and the new measured value is put into the queue. The median a in the queue is obtained by means of fast median filtering. The prediction of background distance b in the use of the measuring distance device means that the thickness value L of the measured object is b_a. The method can well avoid the occasional measurement error, and the selected median value can be used as an accurate measurement value to calculate the thickness of the measured object and to emit the appropriate working tube voltage kvp.

【技术实现步骤摘要】
一种厚度测量的中值滤波方法
本专利技术属于厚度测量
，尤其涉及厚度测量中值滤波领域，具体涉及一种厚度测量的中值滤波方法。
技术介绍
X光是一种射线，就是我们常说的X射线，是一种有能量的电磁波或辐射。当高速移动的电子撞击任何形态的物质时，X光便有可能发生。X光具有穿透性，对不同密度的物质有不同的穿透能力。在医学上X光用来投射人体形成影像，用来辅助诊断；或照射病灶用于治疗。对于准确获取待拍摄病患所需X-射线辐射质量的需求由来已久，X-射线辐射质量主要体现在X-射线球管，X-射线源的工作管电压kvp，以及工作电流积mA·s，前者决定X-射线球管发出的X-射线的能级，更高的管电压对应更高的X-射线能级，也就对应着更强的X-射线穿透力；后者为电流与时间的乘积，共同决定了发出X-射线的量，更高的电流，更长的时间，对应更多的射线量，其中，X-射线源的工作管电压kvp决定了X-射线穿透力，过高的工作管电压kvp会产生更高能级以及更多的射线量，自然会对患者产生更大的伤害；而过低的工作管电压kvp会产生更低能级以及更少的射线量，则无法产生清晰准确的图像。同时，由于患者辐射区域的身体表面不是平面，具有一定的高度起伏，也就是每个部位的厚度不同，若采用同一大小的工作管电压kvp，则会出现一部分图像不清晰，而一部分身体部位可能遭到过大的射线损伤，因此能够实时根据身体厚度的变化来调整工作管电压kvp大小，是目前X光成像技术的尚待解决的难题。目前现有的X光成像设备，还没有类似测量厚度的系统，并且根据测量的厚度来调整工作管电压kvp，要解决上述问题，身体部位的厚度测量技术以及准确性是关键，在测量厚度时，如果只测量一次，则可能因为设备的测量误差而导致工作管电压kvp严重不符合真实厚度所对应的工作管电压kvp，对人体造成很大的损伤或者成像很不清晰，因此需要进行多次测量，并对多次的测量数据取中值，由此得到可靠的厚度测量值；同时现有的中值处理方法，若测量的数据较多，其算法所对应的处理时间也会增长，导致处理结果跟不上设备发射X射线的更新频率，或者X射线发射更新的频率变小，都会进一步导致设备的工作效率降低。
技术实现思路
本专利技术提供一种厚度测量的中值滤波方法，主要用于X光射线成像技术中，对患者不同身体部位厚度的测量，用于解决现有技术中在测量厚度时，可能因为设备的测量误差而导致工作管电压kvp严重不符合真实厚度所对应的工作管电压kvp，对人体造成很大的损伤或者成像很不清晰的技术问题；同时，本专利技术还能解决现有测量厚度设备的中值处理方法所对应的处理时间也会增长，导致处理结果跟不上设备发射X射线的更新频率，或者X射线发射更新的频率变小，进一步导致设备的工作效率降低的问题。本专利技术所采用的技术方案为：一种厚度测量的中值滤波方法，为解决现有技术中在测量厚度时，可能因为设备的测量误差而导致工作管电压kvp严重不符合真实厚度所对应的工作管电压kvp，对人体造成很大的损伤或者成像很不清晰的问题，包括以下步骤，步骤S1：首先通过测距装置对被测物体的表面进行N次预测量，获取N个测量值si，si表示测距装置到被测物体表面的距离，i表示测量序号，取值范围为i＝1，2，…，N；步骤S2：创建一个容量为n的先进先出队列来存储测量值，当队列里的测量值数量达到n时，队列丢弃最早进入队列的测量值，并将新的测量值放入队列，即丢弃旧的测量值，并加入新的测量值，这样队列就保存了最新的n个测量值；步骤S3：通过快速中值滤波的方式获取队列中的中值a，中值a即为N个预测量值的准确值；步骤S4：预测量测距装置使用中的背景距离b，被测物体的厚度值L即为b-a。通过对一个部位进行多次测量，在多个测量值中选择中值的方式，可以很好地避免偶尔出现的测量误差，选择的中值即可作为准确测量值来计算被测物体的厚度，来发出合适的工作管电压kvp，比如测距装置对被测物体的表面一点进行5次测量，5次测量的结果按照先后顺序依次为27、28、29、29、30，创建的队列容量为3，则依照先进先出的原则，则先的前两个数据将被丢弃，队列中剩余29、29、30三个测量值，然后中值滤波取中值，得到中值a＝29，然后就可以通过预先测量的背景距离得到被测物体的厚度；当然，正常的测量时，一个测量周期所测量的数据远大于5次，一般为20-50次，并且队列的容量也远大于3，一般是几十到几百不等。进一步的，为解决现有测量厚度设备的中值处理方法所对应的处理时间也会增长，导致处理结果跟不上设备发射X射线的更新频率，或者X射线发射更新的频率变小，进一步导致设备的工作效率降低的技术问题，步骤S3中，快速中值滤波算法的排序步骤具体包括，步骤S301：预先设置预测量值的取值范围[m，M]，创建一个容量为M-m+1的数组A，所有元素初始化为0，当一个新测量值mk加入队列，则将Amk增大1，当一个旧测量值mj被丢弃，则将Amj减小1；数组A记录了从m到M的各测量值出现的次数；步骤S302：对于该数组，同时从前往后累计局部和若εi-1＜n/2，且εi＞n/2，则测量值i-m+1即为中值。通过快速中值滤波，可以大大减小算法运算的时间，特别是在处理大量数据的时候，可以很快地得到中值，提高工作效率；继续假设有五个测量值27、28、29、29、30，即n＝5，数组A的取值范围为[25，30]，则在输入五个测量值后，数组A即变成[0，0，1，1，2，1]，可得到当i＝3时，εi-1＜n/2，εi＞n/2，由此得到中值为29。通常的，在进行患者身体部位测量时，数组A的取值范围一般为[20，80]，而在测量时，测量值都是以厘米为单位的整数，因此在取值范围内就会有61个取值点，而在测量时，一个测量周期的测量次数可能为几百次，如果按照传统的中值滤波方法，则需要进行不断地排序和运算，而且数组的容量也为几百个，因此算法的耗时较长，而采用本快速中值滤波后，数组的容量为不变的61个，变化的是数组中Amj的增大和减小，此快速中值滤波提供了一种时间复杂度为O(1)的方式，使得算法运算时间大大减小，大大提高了工作效率。进一步的，为了进一步提高算法处理效率，快速排序算法根据新旧测量值不断调整缩小取值范围[m，M]，其缩小取值范围的步骤包括，步骤S3011：当一个新测量值mk加入队列时，同时比较mk和m的大小，若M＞mk＞m，则将取值范围[m，M]变换为[mk，M]，此时数组A的容量变为M-mk+1；当另一个新测量值mp加入队列时，同时比较mp和mk的大小，若mk＞mp＞m，则将取值范围[mk，M]变换为[mp，mk]，此时数组A的容量变为mk-mp+1；若M＞mp＞mk，则将取值范围[mk，M]变换为[mk，mp]，此时数组A的容量变为mp-mk+1；步骤S3012：当再有新测量值mn加入队列时，重复步骤S3011对比mn和mp或mk的大小，重新确定数组的容量范围及容量大小。比如在测量时，假定数组A的取值范围为[20，80]，但是对于某些部位，其厚度范围可以确定改为50-60，则数组A中会有大量为0的Amk，因此会浪费处理器的运算时间，如果能够将范围迅速缩小到实际测量范围，则处理器在进行中值滤波的时候，会节省很多时间；因此本缩小取值范围的步骤的具体实现为，假设原本数组A的取值范围为[20本文档来自技高网...

【技术保护点】
1.一种厚度测量的中值滤波方法，其特征在于：包括以下步骤，步骤S1：首先通过测距装置对被测物体的表面进行N次预测量，获取N个测量值si，si表示测距装置到被测物体表面的距离，i表示测量序号，取值范围为i＝1，2，…，N；步骤S2：创建一个容量为n的先进先出队列来存储测量值，当队列里的测量值数量达到n时，队列丢弃最早进入队列的测量值，并将新的测量值放入队列，即丢弃旧的测量值，并加入新的测量值，这样队列就保存了最新的n个测量值；步骤S3：通过快速中值滤波的方式获取队列中的中值a，中值a即为N个预测量值的准确值；步骤S4：预测量测距装置使用中的背景距离b，被测物体的厚度值L即为b‑a。

【技术特征摘要】
1.一种厚度测量的中值滤波方法，其特征在于：包括以下步骤，步骤S1：首先通过测距装置对被测物体的表面进行N次预测量，获取N个测量值si，si表示测距装置到被测物体表面的距离，i表示测量序号，取值范围为i＝1，2，…，N；步骤S2：创建一个容量为n的先进先出队列来存储测量值，当队列里的测量值数量达到n时，队列丢弃最早进入队列的测量值，并将新的测量值放入队列，即丢弃旧的测量值，并加入新的测量值，这样队列就保存了最新的n个测量值；步骤S3：通过快速中值滤波的方式获取队列中的中值a，中值a即为N个预测量值的准确值；步骤S4：预测量测距装置使用中的背景距离b，被测物体的厚度值L即为b-a。2.根据权利要求1所述的一种厚度测量的中值滤波方法，其特征在于：所述步骤S3中，快速中值滤波算法的排序步骤具体包括，步骤S301：预先设置预测量值的取值范围[m，M]，创建一个容量为M-m+1的数组A，所有元素初始化为0，当一个新测量值mk加入队列，则将Amk增大1，当一个旧测量值mj被丢弃，则将Amj减小1；数组A记录了从m到M的各测量值出现的次数；步骤S302：对于该数组，同时从前往后累计局部和若εi-1＜n/2，且εi＞n/2，则测量值i-m+1即为中值。3.根据权利要求2所述的一种厚度测量的中值滤波方法，其特征在于：所述快速排序算法根据新旧测量值不断调整缩小取值范围[m，M]，其缩小取值范围的步骤包括，...

【专利技术属性】
技术研发人员：张军，何香颖，侯雨舟，
申请(专利权)人：晓智科技成都有限公司，
类型：发明
国别省市：四川,51