CN113589208A - 射频系统的频率确定方法、装置、磁共振设备及存储介质 - Google Patents
射频系统的频率确定方法、装置、磁共振设备及存储介质 Download PDFInfo
- Publication number
- CN113589208A CN113589208A CN202110835685.1A CN202110835685A CN113589208A CN 113589208 A CN113589208 A CN 113589208A CN 202110835685 A CN202110835685 A CN 202110835685A CN 113589208 A CN113589208 A CN 113589208A
- Authority
- CN
- China
- Prior art keywords
- result
- frequency
- determining
- line information
- spectral line
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 108
- 230000003595 spectral effect Effects 0.000 claims abstract description 147
- 230000009466 transformation Effects 0.000 claims abstract description 135
- 238000013519 translation Methods 0.000 claims abstract description 126
- 238000012545 processing Methods 0.000 claims abstract description 101
- 238000003384 imaging method Methods 0.000 claims description 16
- 238000010606 normalization Methods 0.000 claims description 15
- 238000001228 spectrum Methods 0.000 claims description 13
- 238000004590 computer program Methods 0.000 claims description 6
- 230000001131 transforming effect Effects 0.000 claims 1
- 238000012937 correction Methods 0.000 description 46
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 11
- 230000008569 process Effects 0.000 description 10
- 239000000126 substance Substances 0.000 description 10
- 238000010586 diagram Methods 0.000 description 6
- 238000011430 maximum method Methods 0.000 description 6
- 230000003287 optical effect Effects 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 230000007423 decrease Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000011158 quantitative evaluation Methods 0.000 description 4
- 230000035945 sensitivity Effects 0.000 description 4
- 241000894007 species Species 0.000 description 4
- 125000004435 hydrogen atom Chemical group [H]* 0.000 description 3
- 238000012805 post-processing Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 210000003205 muscle Anatomy 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 239000004519 grease Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 230000008707 rearrangement Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000000741 silica gel Substances 0.000 description 1
- 229910002027 silica gel Inorganic materials 0.000 description 1
- 238000000264 spin echo pulse sequence Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
- G01R33/32—Excitation or detection systems, e.g. using radio frequency signals
- G01R33/34—Constructional details, e.g. resonators, specially adapted to MR
- G01R33/34007—Manufacture of RF coils, e.g. using printed circuit board technology; additional hardware for providing mechanical support to the RF coil assembly or to part thereof, e.g. a support for moving the coil assembly relative to the remainder of the MR system
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
- G01R33/38—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/148—Wavelet transforms
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Algebra (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明实施例公开了一种射频系统的频率确定方法、装置、磁共振设备及存储介质,该方法通过预设变换尺度对频率谱线信息进行小波变换,得到频率谱线信息的小波变换的第一结果,并基于平移频率间距以及第一结果,确定出频率谱线信息的小波变换平移的第二结果,进而对第一结果和第二结果进行加权处理,得到小波变换平移相加的第三结果,实现了通过小波变换以及平移相加的方式,综合了频率谱线信息中各特征峰以及与各特征峰相邻的特征峰信息,确定出射频系统的中心频率,提高了中心频率的准确性。并且,该射频系统的频率确定方法适用于各种场强,在不同的场强下均可基于该方法准确提取出射频系统的中心频率。
Description
技术领域
本发明实施例涉及磁共振技术领域,尤其涉及一种射频系统的频率确定方法、装置、磁共振设备及存储介质。
背景技术
在磁共振扫描时,为了得到被扫描者身体组织中特定物质的成像核素的分布信息,如水,脂肪等特定物质,通常需要将磁共振设备射频系统的中心频率对准特定物质中的成像核素的拉莫频率,以保证选层正确、无伪影,达到最高的扫描效率。
然而,由于不同的特定物质中所包含的成像核素的拉莫频率不同,如,脂肪中的氢原子与水中氢原子的拉莫频率存在差异,不同成像核素之间的差异,以及人体组织对成像所需的主磁场产生的影响,主磁场的强度随着时间变化等因素,导致磁共振信号的频率谱线发生变化。因此,需要在进行临床图像采集前调整射频系统的中心频率,以达到最好的成像效果。
现有技术通常采用极值分析、拟合分析等方法,从频率谱线中提取出中心频率。然而,现有技术存在提取到的中心频率不准确,且无法适应于各种场强的技术问题。
发明内容
本发明实施例提供了一种射频系统的频率确定方法、装置、磁共振设备及存储介质,以提高提取射频系统的中心频率的准确性,适用于各种场强下的磁共振系统。
第一方面,本发明实施例提供了一种射频系统的频率确定方法,可用于确定磁共振设备扫描时射频系统的频率确定,包括:
获取扫描对象的频率谱线信息;
基于预设变换尺度对所述频率谱线信息进行小波变换,得到所述频率谱线信息的第一结果;
基于平移频率间距以及所述第一结果确定第二结果,对所述第一结果和所述第二结果进行加权处理,得到第三结果;
基于所述第三结果确定射频系统对应的目标中心频率。
可选的,在所述基于平移频率间距以及所述第一结果确定第二结果之前,还包括:
对所述第一结果中的正值部分和/或负值部分进行修正处理;
基于修正处理结果更新所述第一结果。
可选的,所述基于预设变换尺度对所述频率谱线信息进行小波变换,得到所述频率谱线信息的第一结果,包括:
对进行小波变换的小波基中的正值部分和/或负值部分进行修正处理,基于修正处理结果更新所述小波基;
基于所述预设变换尺度和所述小波基对所述频率谱线信息进行小波变换,得到所述频率谱线信息的第一结果。
可选的,所述基于所述第三结果确定射频系统对应的目标中心频率,包括:
若所述第三结果的数量为至少两个,则确定各所述第三结果对应的候选中心频率;
基于各所述候选中心频率确定射频系统对应的目标中心频率。
可选的,所述基于平移频率间距以及所述第一结果确定第二结果,包括:
基于所述第一结果以及至少两个平移频率间距确定至少两个第二结果;
相应的,所述对所述第一结果和所述第二结果进行加权处理,得到第三结果,包括:
对所述第一结果和至少两个所述第二结果进行加权处理,得到至少两个第三结果。
可选的,所述基于预设变换尺度对所述频率谱线信息进行小波变换,得到所述频率谱线信息的第一结果,包括:
基于至少两个预设变换尺度对所述频率谱线信息进行小波变换,得到所述频率谱线信息的至少两个第一结果;
相应的,所述基于平移频率间距以及所述第一结果确定第二结果,对所述第一结果和所述第二结果进行加权处理,得到第三结果,包括:
基于平移频率间距以及至少两个所述第一结果确定至少两个第二结果,对所述第一结果和至少两个所述第二结果进行加权处理,得到至少两个第三结果。
所述基于所述第三结果确定射频系统对应的目标中心频率,包括:
若所述第三结果包括多个特征峰,则通过全局最大法或全局频率比对法,在第三结果中确定目标特征峰;
基于所述目标特征峰对应的频率确定射频系统对应的目标中心频率。
可选的,在所述基于预设变换尺度对所述频率谱线信息进行小波变换之前,还包括:
对所述频率谱线信息进行归一化处理,基于归一化处理结果更新频率谱线信息;
相应的,所述基于所述第三结果确定射频系统对应的目标中心频率,包括:
基于第三结果确定射频系统对应的归一化中心频率,基于归一化中心频率确定射频系统对应的目标中心频率。
第二方面,本发明实施例还提供了一种射频系统的频率确定方法,包括:
获取扫描对象的频率谱线信息;
基于平移频率间距对第一小波基进行平移处理,得到第二小波基,对所述第一小波基和所述第二小波基进行加权处理,得到第三小波基;
基于预设变换尺度和所述第三小波基对所述频率谱线信息进行小波变换,得到所述频率谱线信息的第三结果;
基于所述第三结果确定射频系统对应的目标中心频率。
第三方面,本发明实施例还提供了一种射频系统的频率确定装置,包括:
频率谱线获取模块,用于获取扫描对象的频率谱线信息;
第一结果确定模块,用于基于预设变换尺度对所述频率谱线信息进行小波变换,得到所述频率谱线信息的第一结果;
第三结果确定模块,用于基于平移频率间距以及所述第一结果确定第二结果,对所述第一结果和所述第二结果进行加权处理,得到第三结果;
中心频率确定模块,用于基于所述第三结果确定射频系统对应的目标中心频率。
第四方面,本发明实施例还提供了一种磁共振设备,包括:
处理器,所述处理器用于在所述磁共振设备对扫描对象进行成像扫描前,基于本发明任意实施例提供的射频系统的频率确定方法,确定所述磁共振设备的射频系统的目标中心频率。
第五方面,本发明实施例还提供了一种电子设备,所述电子设备包括:
一个或多个处理器;
存储装置,用于存储一个或多个程序,
当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现如本发明任意实施例提供的射频系统的频率确定方法。
第六方面,本发明实施例还提供了一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现如本发明任意实施例提供的射频系统的频率确定方法。
上述发明中的实施例具有如下优点或有益效果:
通过预设变换尺度对频率谱线信息进行小波变换,得到频率谱线信息的小波变换的第一结果,并基于平移频率间距以及第一结果,确定出频率谱线信息的小波变换平移的第二结果,进而对第一结果和第二结果进行加权处理,得到小波变换平移相加的第三结果,实现了通过小波变换以及平移相加的方式,综合了频率谱线信息中各特征峰以及与各特征峰相邻的特征峰信息,确定出射频系统的中心频率,提高了中心频率的准确性。并且,该射频系统的频率确定方法适用于各种场强,在不同的场强下均可基于该方法准确提取出射频系统的中心频率。
附图说明
为了更加清楚地说明本发明示例性实施例的技术方案,下面对描述实施例中所需要用到的附图做一简单介绍。显然,所介绍的附图只是本发明所要描述的一部分实施例的附图,而不是全部的附图,对于本领域普通技术人员,在不付出创造性劳动的前提下,还可以根据这些附图得到其他的附图。
图1A为本发明实施例所提供的一种射频系统的频率确定方法的流程示意图;
图1B为本发明实施例所提供的一种归一化处理结果;
图2为本发明实施例所提供的另一种射频系统的频率确定方法的流程示意图;
图3为本发明实施例所提供的另一种射频系统的频率确定方法的流程示意图;
图4为本发明实施例所提供的另一种射频系统的频率确定方法的流程示意图;
图5A为本发明实施例所提供的另一种射频系统的频率确定方法的流程示意图;
图5B为本发明实施例所提供的两种可选的射频系统的频率确定方法的流程示意图;
图6为本发明另一实施例所提供的一种射频系统的频率确定装置的结构示意图;
图7为本发明另一实施例所提供的一种电子设备的结构示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是,此处所描述的具体实施例仅仅用于解释本发明,而非对本发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本发明相关的部分而非全部结构。
图1A为本发明实施例提供的一种射频系统的频率确定方法的流程示意图,本实施例可适用于基于扫描对象的频率谱线信息提取射频系统的中心频率的情况,该方法可以由射频系统的频率确定装置来执行,该装置可以由硬件和/或软件来实现,该方法具体包括如下步骤:
在对本实施例提供的射频系统的频率确定方法进行详细介绍之前,对本实施例提供的射频系统的频率确定方法的应用场景进行示例性说明。示例性的,本方法可以适用于在不同的主磁场强度下,或在磁共振系统的主磁场的强度发生变化时,如,主磁场强度由6T变为7T;或者,扫描对象发生变化时,如,头部变为腿部;又或者,成像核素发生变化时,如氢变为其它成像核素,等,均可以根据本申请提供的射频系统的频率确定方法,重新确定出当前扫描对象对应的射频系统的中心频率。所述射频系统可以是和磁共振设备集成一体的,也可以是单独的射频系统;所述射频系统可以实现射频信号的发射和接收,也可以仅仅用于射频信号的发射。
S110、获取扫描对象的频率谱线信息。
其中,频率谱线信息的横轴为频率,纵轴为信号强度;频率谱线信息可以反映出不同频率处的磁共振强度信息。具体的,可以使用预设采集序列对扫描对象进行扫描,从扫描得到的磁共振信号中提取扫描对象的频率谱线信息。扫描对象可以是被扫描部位。可选的,预设采集序列可以是自旋回波序列。扫描可以是定位像扫描。
示例性的,可以通过预设采集序列对扫描对象进行预扫描,得到磁共振数据,通过对磁共振数据进行傅里叶变换,得到频率谱线信息。频率谱线信息中可以包含随着频率的变化,幅度也发生变化所形成的峰和谷。需要说明的是,频率谱线信息中峰的数量可以是一个,也可以是多个,本申请对此不作限定。当然,若频率谱线信息中峰的数量为一个,则频率谱线信息中不包含谷。
S120、基于预设变换尺度对频率谱线信息进行小波变换,得到频率谱线信息的第一结果。
其中,预设变换尺度可以是预先设置的对频率谱线信息进行小波变换的尺度。具体的,预设变换尺度可以是预先设置的用于控制小波变换伸缩量的参数。预设变换尺度的大小与频率谱线信息进行小波变换后所得到的第一结果中的峰的展宽相关。具体的,不同第一结果中峰的展宽在不同的预设变化尺度下有不同的效果。例如,若第一结果中的峰的展宽越大,则使用越大的预设变换尺度可以得到更大的第一结果值。
示例性的,预设变换尺度可以是0.1ppm、0.2ppm、0.4ppm、0.6ppm等,其具体可以根据实际扫描情况进行调整,本申请对此不作限定。
当然,预设变换尺度的数量可以是一个,也可以是多个;当预设变换尺度为至少两个时,可以基于至少两个预设变换尺度对频率谱线信息进行小波变换,得到频率谱线信息的至少两个第一结果。
需要说明的是,本申请对频率谱线信息进行小波变换的目的在于:对频率谱线信息中的峰和谷的特征进行量化评价,并综合考虑频率谱线信息中峰的展宽、偏态、间距等特征,从而实现正确得到目标中心频率的技术效果。
S130、基于平移频率间距以及第一结果确定第二结果,对第一结果和第二结果进行加权处理,得到第三结果。
其中,平移频率间距可以是预先设置的对第一结果进行平移的间距。在一种实施方式中,基于平移频率间距以及第一结果确定第二结果,可以是:基于平移频率间距对第一结果进行平移处理,得到第二结果。即,直接通过平移频率间距对第一结果进行平移,得到第二结果。
在另一种实施方式中,基于平移频率间距以及第一结果确定第二结果,还可以是:基于平移频率间距对第一结果对应的第一小波基进行平移处理,得到第二小波基,并基于第二小波基以及预设变换尺度对频率谱线信息进行小波变换,得到第二结果。其中,第一结果对应的第一小波基可以是进行小波变换的初始小波基,如,母小波基。换言之,该可选的方法不对第一结果直接进行平移,而是利用第一结果对应的第一小波基,将小波变换后的平移处理放置在第一小波基的处理中,对进行小波变换时所应用的第一小波基进行平移,基于平移后得到的第二小波基对频率谱线信息进行小波变换,生成第二结果。即,通过平移第一小波基,再基于平移后的小波基进行小波变换的方式,代替直接对小波变换后得到的第一结果进行平移的方式。
在本实施例中,平移频率间距可以根据成像核素进行调整。在一种实施方式中,平移频率间距可以基于成像核素各特征峰间的理论化学位移差确定,如,平移频率间距可以是将理论化学位移差扩大预设范围后得到的区间内的任意值。示例性的,成像核素为氢质子时,水脂对的理论化学位移差一般为3.4ppm,则平移频率间距可以是3.4ppm;或者,预设平移间距还可以是3.4ppm扩大预设范围后的区间内的任一数值。当然,理论化学位移差扩大的范围其可以根据实际成像需求进行调整。在另一种实施方式中,还可以根据实际测试效果确定平移频率间距,如,可以在模体上、扫描对象或动物上进行模拟测试,根据模拟测试数据确定平移频率间隔;又或者,还可以根据历史数据得到预设平移间距。
需要说明的是,平移频率间距的数量可以是一个,也可以是多个,本申请对此不作限定。例如,若要识别的第一结果中特征峰对的数量大于一个,则可以在特征峰中确定一个主特征峰,各平移频率间距可以分别为除主特征峰之外的其余特征峰相对于主特征峰的偏移。示例性的,特征峰对的数量为两个,平移频率间距可以分别为a-b峰对间的偏移、a-c峰对间的偏移。具体的,若平移频率间距的数量为至少两个,则针对每一个平移频率间距,均可以基于平移频率间距以及第一结果确定出第二结果,从而得到至少两个第二结果。
具体的,在得到第二结果后,对第一结果以及第二结果进行加权计算,得到第三结果。其中,在加权处理时,可以对第一结果对应的第一权重以及第二结果对应的第二权重进行预先设置。示例性的,针对H谱,水峰的频率比脂主峰高,因此,可以保证正向平移时第一权重大于第二权重,即,若平移频率间距大于零,则第一权重大于第二权重,如,第一权重为1.1,第二权重为1,第一权重略大于第二权重;若平移频率间距小于零(反向平移时),则第一权重小于第二权重,如,第一权重为1,第二权重为1.1,第二权重略大于第一权重。
可选的,在设置第一权重以及第二权重时,保证加权处理时的主特征峰的权重高于其余特征峰的权重。
本实施例在进行小波变换以及平移后,对小波变换得到的第一结果和平移处理后的第二结果进行加权计算的目的在于:在对频率谱线信息进行量化评估之后,综合考虑各特征峰之间的位置关系,以增强符合平移频率间距的特征峰对得到的量化评估结果,进而提高提取出的目标中心频率的准确性。
比如,水脂的化学位移反映在频率谱线信息的偏移上,将小波变换的结果进行平移相加,有水脂对存在时,水与脂的峰都会对小波变换结果产生贡献,从而可提取出潜在的水脂对;而只有单峰存在时,仅有单峰处有较大的正值。通过此方法也能将水、脂、硅胶的化学位移纳入考虑,从而应用于各种临床场景。例如,在对A峰(如,水峰)进行分析时,通过平移、加权一并将其相邻的B峰(如,脂肪峰)纳入考虑,平移、加权处理后的A峰所在频率处的第三结果包含原A峰信息以及与其相邻的B峰信息,此时,不同频率处的量化评估结果将是多个特征峰评价结果的综合,从而可以提高分析出的中心频率的准确性。
示例性的,对第一结果和第二结果进行加权处理,得到第三结果,满足如下公式:
wt(fn)=d0×wx(fn)+d1×wx(fn-x),
其中,wt(fn)为第三结果,wx(fn)为第一结果,fn为频率谱线信息中的频率或归一化频率,wx(fn-x)为第二结果,x表示平移频率间距,d0、d1分别为第一结果对应的加权权重值、第二结果对应的加权权重值。
当然,若要识别的第一结果中特征峰对的数量大于一个,相应的,平移频率间距的数量大于一个,基于平移频率间距所得到的第二结果的数量也大于一个;此时,对第一结果和第二结果进行加权处理,得到第三结果,满足如下公式:
wt(fn)=d0×wx(fn)+d1×wx(fn-x1)+…+dn×wx(fn-xn)
其中,wt(fn)为第三结果,wx(fn)为第一结果,fn为频率谱线信息中的频率或归一化频率,x1~xn为各特征峰相对主峰的偏移,即各平移频率间距,wx(fn-x1)~wx(fn-xn)分别为各平移频率间距对应的第二结果,d0为第一结果对应的加权权重值,d1~dn为各个第二结果对应的加权权重值。
S140、基于第三结果确定射频系统对应的目标中心频率。
其中,第三结果为频率谱线信息经过小波变换、平移、平移后相加的结果。第三结果的数量可以为一个,也可以为多个。
具体的,当只有一个第三结果时,该结果包含的量化评价将直接用于确定目标中心频率。可选的,基于第三结果确定射频系统对应的目标中心频率,包括:若第三结果包括多个特征峰,则通过全局最大法或全局频率比对法,在第三结果中确定目标特征峰;基于目标特征峰对应的频率确定射频系统对应的目标中心频率。
其中,全局最大法可以是将最大值特征峰作为目标特征峰;全局频率比对法可以是将最大值特征峰和次大值特征峰中频率较低特征峰作为目标特征峰;,最大值特征峰为幅值最大的特征峰,次大值特征峰为幅值次大值的特征峰。目标特征峰的筛选规则还可以根据实际需求进行设置,本申请对此不作限定。
具体的,在对扫描对象进行压脂相关序列的扫描时,如果射频系统的频率校准的第三结果中包括两个特征峰,此时,选择频率低的特征峰作为目标特征峰,可能会出现脂肪抑制效果降低但肌肉水信号不会损失的情况;选频率高的特征峰作为目标特征峰,可能会出现脂肪抑制效果好但肌肉水信号丢失的情况。因此,综合考虑上述两种情况,本方法可以在第三结果中包括多个特征峰时,可以在最大值特征峰和次大值特征峰中选择频率较低的特征峰作为目标特征峰。
又或者,基于第三结果确定射频系统对应的目标中心频率,还可以是:对第三结果进行后处理,基于后处理的结果确定射频系统对应的目标中心频率。其中,后处理可以是对第三结果进行降噪处理,滤除第三结果中的孤立噪声点。
需要说明的是,在第三结果的数量为多个时,可以将各第三结果中的最大值特征峰作为各第三结果对应的目标特征峰;或者,将各第三结果中最大值特征峰和次大值特征峰中频率较低的特征峰对应的频率作为该第三结果的代表值,对各个第三结果进行横向比较,将各第三结果对应的代表值中最大的代表值对应的特征峰作为目标特征峰。又或者,还可以综合考虑所有第三结果中量化评价结果值最高的频率,作为目标中心频率,得到该结果所用的小波变换尺度及预设平移间距即是描述该谱线的峰的展宽和特征峰对间距的相关指标。
本实施例的技术方案,通过预设变换尺度对频率谱线信息进行小波变换,得到频率谱线信息的小波变换的第一结果,并基于平移频率间距以及第一结果,确定出频率谱线信息的小波变换平移的第二结果,进而对第一结果和第二结果进行加权处理,得到小波变换平移相加的第三结果,实现了通过小波变换以及平移相加的方式,综合了频率谱线信息中各特征峰以及与各特征峰相邻的特征峰信息,确定出射频系统的中心频率,提高了中心频率的准确性。并且,该射频系统的频率确定方法适用于各种场强,在不同的场强下均可基于该方法准确提取出射频系统的中心频率。
在一种可选的实施方式中,还可以在进行小波变换、平移、加权之前,对频率谱线信息进行频率归一化处理,以进一步适用于各种主磁场场强下提取射频系统的中心频率。即,在对频率谱线信息进行小波变换之前,还包括:对频率谱线信息进行归一化处理,基于归一化处理结果更新频率谱线信息;相应的,基于第三结果确定射频系统对应的目标中心频率,包括:基于第三结果确定射频系统对应的归一化中心频率;基于归一化中心频率确定射频系统对应的目标中心频率。
其中,对频率谱线信息进行归一化处理,可以是:基于预设中心频率对频率谱线信息进行归一化处理。具体的,可以是将频率谱线信息除以预设中心频率。其中,预设中心频率可以是射频系统默认频率,也可以是预先估计的中心频率。对频率谱线信息进行归一化处理,还可以是:频率谱线信息的频率除以主场强与旋磁比的积。示例性的,如图1B所示,展示了一种归一化处理结果。在归一化处理结果中,纵轴不变,横轴的单位变为ppm,即,从强度随频率变化的频率谱线信息变为强度随归一化频率变化的归一化频率谱线信息。
具体的,在对频率谱线信息进行归一化处理之后,第一结果、第二结果以及第三结果的横轴均为归一化频率,此时,基于第三结果确定射频系统对应的目标中心频率,可以是先基于第三结果确定射频系统对应的归一化中心频率;再基于归一化中心频率确定射频系统对应的目标中心频率。其中,归一化中心频率的确定方式可参见上述方法,对此不再进行赘述。
在该可选的实施方式中,通过对频率谱线信息进行归一化处理,进一步地减少了场强对提取到的中心频率的影响,从而进一步地提高了本实施例提供的射频系统的频率确定方法的适用范围,使得方法可被应用于各种主场强下。示例性的,在扫描场强发生变化时,如,1.5T到5T时,仍可以采用本方法提取射频系统的中心频率。
可选的,本实施例提供的射频系统的频率确定方法还包括:基于目标中心频率,对射频系统的中心频率进行调整,以基于射频系统调整后的中心频率对扫描对象进行扫描成像。
本实施例的技术方案,通过预设变换尺度对频率谱线信息进行小波变换,得到频率谱线信息的小波变换的第一结果,并基于平移频率间距以及第一结果,确定出频率谱线信息的小波变换平移的第二结果,进而对第一结果和第二结果进行加权处理,得到小波变换平移相加的第三结果,实现了通过小波变换以及平移相加的方式,综合了频率谱线信息中各特征峰以及与各特征峰相邻的特征峰信息,确定出射频系统的中心频率,提高了中心频率的准确性。并且,该射频系统的频率确定方法适用于各种场强,在不同的场强下均可基于该方法准确提取出射频系统的中心频率。
图2为本发明实施例提供的另一种射频系统的频率确定方法的流程示意图,可选的,在所述对所述第一结果和所述第二结果进行加权处理,得到第三结果之前,还包括:对所述第一结果中的正值部分和/或负值部分进行修正处理;基于修正处理结果更新所述第一结果。
其中与上述相同或相应的术语的解释在此不再赘述。参见图2,本实施例提供的另一种射频系统的频率确定方法包括如下步骤:
S210、获取扫描对象的频率谱线信息。
S220、基于预设变换尺度对频率谱线信息进行小波变换,得到频率谱线信息的第一结果。
S230、对第一结果中的正值部分和/或负值部分进行修正处理,基于修正处理结果更新第一结果,并基于平移频率间距以及第一结果确定第二结果。
其中,频率谱线信息中峰的小波变换结果为正值部分,谷的小波变换结果为负值部分。或者,频率谱线信息仅包括单峰时,对单峰进行小波变换后,也可以得到正值部分和负值部分。在本实施例中,考虑到在平移后加权时负值部分会降低正值部分的大小,即,平移后加权时特征谷会降低特征峰的峰值,因此,可以在进行平移、加权处理之前,对第一结果的正值部分和/或负值部分进行修正处理,以降低负值部分的强度的绝对值,和/或,增加正值部分的强度的绝对值。
需要说明的是,本实施例可以对第一结果中的正值部分和负值部分同时进行修正处理,在增加正值部分的强度的绝对值的同时,降低负值部分的强度的绝对值,以减少负值部分在加权时对正值部分的影响;还可以仅对第一结果中的正值部分或负值部分进行修正处理。
可选的,对第一结果中的正值部分和/或负值部分进行修正处理,可以是直接对第一结果中的正值部分和/或负值部分的幅值进行调整;也可以是基于预设的正值修正系数和负值修正系数中的至少一种,对第一结果中的正值部分和/或负值部分进行修正处理。如,将正值修正系数与第一结果中的正值部分相乘;和/或,将负值修正系数与第一结果中的负值部分相乘。
其中,预设的正值修正系数可以显著大于负值修正系数。例如,预设的正值修正系数可以为10,预设的负值修正系数可以为1;或者,预设的正值修正系数可以为1,预设的负值修正系数可以为0.1,等。
具体的,可以在第一结果中的正值部分和/或负值部分进行修正处理后,将修正处理结果作为第一结果,进而基于修正更新后的第一结果确定第二结果。其中,第二结果为经过小波变换、修正、平移的结果。
当然,若第一结果中仅包括一个特征峰,即第一结果为单特征峰结果时,仍可以对该单特征峰进行修正处理。本实施例提供的射频系统的频率确定方法无需判断第一结果是否为单特征峰结果,换言之,若第一结果中仅包含一个特征峰,此时仍可对该特征峰进行修正处理。
S240、对第一结果和第二结果进行加权处理,得到第三结果。
在本实施例中,第一结果为经过小波变换、修正处理后的结果。因此,将修正更新后的第一结果与第二结果进行加权处理,可以实现在对特征峰进行量化评估时,减少特征谷对特征峰的影响,调整对特征峰和特征谷的敏感度,提高第三结果的准确性。
S250、基于第三结果确定射频系统对应的目标中心频率。
本实施例的技术方案,在对第一结果和第二结果进行加权处理之前,通过对第一结果中的正值部分和/或负值部分进行修正处理,降低平移加权时特征谷对特征峰的影响,提高第三结果的准确性,从而提高提取到的中心频率的准确性。
图3为本发明实施例提供的另一种射频系统的频率确定方法的流程示意图,可选的,所述基于预设变换尺度对所述频率谱线信息进行小波变换,得到所述频率谱线信息的第一结果,包括:对进行小波变换的小波基中的正值部分和/或负值部分进行修正处理,基于修正处理结果更新所述小波基;基于所述预设变换尺度和所述小波基对所述频率谱线信息进行小波变换,得到所述频率谱线信息的第一结果。
其中与上述相同或相应的术语的解释在此不再赘述。参见图3,本实施例提供的另一种射频系统的频率确定方法包括如下步骤:
S310、获取扫描对象的频率谱线信息。
S320、对进行小波变换的小波基中的正值部分和/或负值部分进行修正处理,基于修正处理结果更新小波基。
本实施例可以在基于预设变换尺度对频率谱线信息进行小波变换之前,对小波变换的小波基的正值部分和/或负值部分进行修正处理,降低负值部分的强度的绝对值,和/或,增加正值部分的强度的绝对值。
需要说明的是,本实施例也可以同时对小波基中的正值部分和负值部分进行修正处理,在增加正值部分的强度的绝对值的同时,降低负值部分的强度的绝对值;或者,还可以仅对第一结果中的正值部分或负值部分进行修正处理。
可选的,对小波基中的正值部分和/或负值部分进行修正处理,可以是直接对小波基中的正值部分和/或负值部分的幅值进行调整;还可以是基于预设的正值修正系数和负值修正系数中的至少一种,对小波基中的正值部分和/或负值部分进行修正处理。其中,预设的正值修正系数可以显著大于预设的负值修正系数。示例性的,正值修正系数为1,负值修正系数为0.1。
S330、基于预设变换尺度和小波基对频率谱线信息进行小波变换,得到频率谱线信息的第一结果。
具体的,在对小波基的正值部分和/或负值部分进行修正处理后,通过该修正后的小波基对频率谱线信息进行小波变换,可以得到已修正的小波变换结果。
S340、基于平移频率间距以及第一结果确定第二结果,对第一结果和第二结果进行加权处理,得到第三结果。
其中,第一结果为已修正的小波变换结果,第二结果为修正、小波变换、平移的结果。通过对小波基进行修正,并基于修正的小波基对频率谱线信息进行小波变换,可以直接得到已修正的小波变换结果,减少了平移相加时特征谷对特征峰的影响,提高了第三结果的准确性。
S350、基于第三结果确定射频系统对应的目标中心频率。
在本实施例中,通过对小波基中的正值部分和/或负值部分进行修正处理,以基于修正后的小波基对频率谱线信息进行小波变换,直接得到已修正的小波变换结果,降低了平移相加时特征谷对特征峰的影响,提高了第三结果的准确性,从而提高了提取到的中心频率的准确性。
图4为本发明实施例提供的另一种射频系统的频率确定方法的流程示意图,可选的,所述基于所述第三结果确定射频系统对应的目标中心频率,包括:若所述第三结果的数量为至少两个,则确定各所述第三结果对应的候选中心频率;基于各所述候选中心频率确定射频系统对应的目标中心频率。
其中与上述相同或相应的术语的解释在此不再赘述。参见图4,本实施例提供的另一种射频系统的频率确定方法包括如下步骤:
S410、获取扫描对象的频率谱线信息,基于预设变换尺度对所述频率谱线信息进行小波变换,得到所述频率谱线信息的第一结果。
S420、基于平移频率间距以及所述第一结果确定第二结果,对所述第一结果和所述第二结果进行加权处理,得到第三结果。
S430、若第三结果的数量为至少两个,则确定各第三结果对应的候选中心频率,基于各候选中心频率确定射频系统对应的目标中心频率。
其中,第三结果基于第一结果和第二结果加权得到,当第一结果和/或第二结果的数量为多个时,第三结果的数量可以是多个。
在一种实施方式中,可以是通过多个预设变换尺度得到多个第三结果。即,可选的,基于预设变换尺度对频率谱线信息进行小波变换,得到频率谱线信息的第一结果,包括:基于至少两个预设变换尺度对频率谱线信息进行小波变换,得到频率谱线信息的至少两个第一结果;相应的,基于平移频率间距以及第一结果确定第二结果,对第一结果和第二结果进行加权处理,得到第三结果,包括:基于平移频率间距以及至少两个第一结果确定至少两个第二结果,对第一结果和至少两个第二结果进行加权处理,得到至少两个第三结果。
在该可选的实施方式中,预设变换尺度的数量为至少两个时,小波变换时所生成的第一结果的数量也为至少两个。相应的,基于平移频率间距也可以得到每一个第一结果所对应的第二结果。通过将每个第一结果与其对应的第二结果进行加权处理,得到至少两个第三结果。
示例性的,预设变换尺度的数量为4个,如,0.1ppm、0.2ppm、0.4ppm、0.6ppm,根据4个预设变换尺度可以得到4个第一结果(如,A1、B1、C1、D1),并通过一个平移频率间距生成每个第一结果对应的第二结果(如,A2、B2、C2、D2),此时,分别将每个第一结果与其对应的第二结果进行加权处理(如,A1+A2、B1+B2、C1+C2、D1+D2),得到4个第三结果。
在另一种实施方式中,还可以是通过平移频率间距得到多个第三结果。即,可选的,基于平移频率间距以及第一结果确定第二结果,包括:基于第一结果以及至少两个平移频率间距确定至少两个第二结果;相应的,对第一结果和第二结果进行加权处理,得到第三结果,包括:对第一结果和至少两个第二结果进行加权处理,得到至少两个第三结果。
在该可选的实施方式中,平移频率间距的数量为至少两个时,生成的第二结果的数量也为至少两个。每个第二结果经过加权处理后均可得到其对应的第三结果,从而得到至少两个第三结果。
示例性的,平移频率间距的数量为4个,如,0.1ppm、0.2ppm、0.4ppm、0.6ppm,根据4个平移频率间距可以得到4个第二结果(如,A2、B2、C2、D2),此时,分别将第一结果与各第二结果分别进行加权处理(如,A1+A2、A1+B2、A1+C2、A1+D2),得到4个第三结果。
需要说明的是,在另一种可选的方式中,可以同时通过多个预设变换尺度和多个平移频率间距得到多个第三结果。具体的,可以是基于各预设变换尺度对频率谱线信息进行小波变换,得到频率谱线信息的各第一结果;基于各平移频率间距以及各第一结果确定各第一结果对应的各第二结果,对各第一结果和各第一结果对应的各第二结果分别进行加权处理,得到各第三结果。即,若预设变换尺度的数量为p,平移频率间距的数量为q,则针对每一个预设变换尺度,均可得到q个第三结果,第三结果的数量为p×q。
例如,预设变换尺度的数量为3个,平移频率间距的数量为4个。通过3个预设变换尺度对频率谱线信息进行小波变换,得到A1、B1、C1;基于4个平移频率间距,确定出A1对应的A21、A22、A23、A24,B1对应的B21、B22、B23、B24,C1对应的C21、C22、C23、C24;对A1以及A1对应的A21、A22、A23、A24分别进行加权处理,得到A1+A21、A1+A22、A1+A23、A1+A24,对B1以及B1对应的B21、B22、B23、B24分别进行加权处理,得到B1+B21、B1+B22、B1+B23、B1+B24,对C1以及C1对应的C21、C22、C23、C24分别进行加权处理,得到C1+C21、C1+C22、C1+C23、C1+C24。即,得到3×4个第三结果。
在本实施例中,当第三结果为多个时,可以针对每一个第三结果,均确定出其对应的候选中心频率,从而在各个候选中心频率中确定出射频系统对应的目标中心频率。其中,各第三结果对应的候选中心频率可以是最大特征峰对应的频率,或者,根据全局频率比对法确定出的特征峰对应的频率,或者,滤除噪声后最大特征峰对应的频率,等。
在一种实施方式中,基于各候选中心频率确定射频系统对应的目标中心频率,可以是:在各候选中心频率对应的特征峰中确定最大值特征峰,将该最大值特征峰对应的中心频率确定为目标中心频率;或者,将各候选中心频率进行横向比较,确定出最大值候选中心频率,将该最大值候选中心频率作为目标中心频率。
本实施例的技术方案,考虑到不同小波变换尺度的小波对不同展宽的峰的敏感度不同,不同的平移频率间距的小波对不同间距的峰的敏感度不同,因此,采用不同的小波变换尺度和/或平移频率间距,可以增强本实施例提供的射频系统的频率确定方法的适用性,进而可以综合各预设变换尺度和/或平移频率间距下计算得到的候选中心频率,确定出最优的目标中心频率,进一步提高了提取出的射频系统的中心频率的准确性。
图5A为本发明实施例提供的另一种射频系统的频率确定方法的流程示意图。由于小波变换属于线性变换,因此,对小波变换后的结果进行平移、平移后相加,与对小波基进行平移、平移后相加是等价的,因此,本实施例还可以将平移以及平移后相加的处理放在小波基的处理中,可以适用于基于扫描对象的频率谱线信息提取射频系统的中心频率的情况,该方法具体包括如下步骤,其中与上述相同或相应的术语的解释在此不再赘述:
S510、获取扫描对象的频率谱线信息。
S520、基于平移频率间距对第一小波基进行平移处理,得到第二小波基,对第一小波基和第二小波基进行加权处理,得到第三小波基。
其中,第一小波基可以是进行小波变换的初始小波基;第二小波基可以是平移后的小波基。在一种实施方式中,基于平移频率间距对第一小波基进行平移处理,得到第二小波基,可以是:基于平移频率间距对第一结果进行平移处理,得到第二小波基。
在得到平移的第二小波基后,可以对第一小波基和第二小波基进行加权处理,得到第三小波基。其中,第三小波基可以理解为平移、平移后相加的小波基。在加权处理时第一小波基对应的第一权重和第二小波基对应的第二权重的设置可以参考上述实施例。
S530、基于预设变换尺度和第三小波基对频率谱线信息进行小波变换,得到频率谱线信息的第三结果。
其中,第三结果为频率谱线信息经过平移、加权处理以及小波变换的结果。本实施例通过对小波基进行平移以及平移后相加,并基于平移以及平移后相加的小波基进行小波变换,可以直接得到频率谱线信息的小波变换、平移以及平移后相加的结果。
具体的,本实施例可以通过将上述实施例中小波变换后的平移、加权的处理,变换为对小波基进行平移、加权的处理。即,通过在进行小波变换之前,对小波基进行平移、加权的处理,代替上述实施例中在小波变换后的平移、加权处理,得到频率谱线信息的小波变换与平移加权结合的结果。
当然,还可以是在对第一小波基和第二小波基进行加权处理得到第三小波基之前,对第一小波基进行平移处理之后,对第一小波基的正值部分和/或负值部分进行修正,基于修正的结果更新第一小波基,以减少第一小波基和第二小波基相加时特征谷对特征峰的影响。
S540、基于第三结果确定射频系统对应的目标中心频率。
可选的,基于第三结果确定射频系统对应的目标中心频率,包括:若第三结果的数量为至少两个,则确定各第三结果对应的候选中心频率;基于各候选中心频率确定射频系统对应的目标中心频率。
其中,可以通过至少两个平移频率间距和/或至少两个预设变换尺度得到至少两个第三结果。
如,基于至少两个平移频率间距对第一小波基进行平移处理,得到至少两个第二小波基,对第一小波基和至少两个第二小波基进行加权处理,得到至少两个第三小波基;基于预设变换尺度和至少两个第三小波基对频率谱线信息进行小波变换,得到频率谱线信息的至少两个第三结果。
或者,基于至少两个预设变换尺度和第三小波基对频率谱线信息进行小波变换,得到频率谱线信息的至少两个第三结果。
或者,基于至少两个预设变换尺度和至少两个平移频率间距,得到各平移频率间距对应的第三小波基,基于各预设变换尺度分别对各平移频率间距对应的第三小波基进行小波变换,得到多个第三结果。
本实施例的技术方案,通过对小波基进行平移以及平移后相加,并对平移以及平移后相加的小波基进行小波变换,可以直接得到频率谱线信息小波变换、平移以及平移后相加的结果,实现了通过小波变换以及平移相加的方式,综合了频率谱线信息中各特征峰以及与各特征峰相邻的特征峰信息,确定出射频系统的中心频率,提高了中心频率的准确性。并且,该射频系统的频率确定方法适用于各种场强,在不同的场强下均可基于该方法准确提取出射频系统的中心频率。
本实施例还提供两种可选的射频系统的频率确定方法,如图5B所示,展示了两种可选的射频系统的频率确定方法的流程示意图。通过图中的任一种方法,均可确定出射频系统的中心频率。
结合图5B所示内容,其中一种射频系统的频率确定方法包括如下步骤:获取磁共振信号;对磁共振信号进行傅里叶变换,得到频率谱线信息;对频率谱线信息进行归一化处理,基于归一化处理结果更新频率谱线信息;对频率谱线信息进行小波变换,得到第一结果;对第一结果的正值部分和/或负值部分进行修正处理,基于修正处理的结果更新第一结果;对第一结果进行平移处理,得到第二结果;对第一结果和第二结果进行加权相加,得到第三结果;在第三结果中确定出目标特征峰,从而基于目标特征峰对应的频率确定目标中心频率,其中,确定方法可以是全局最大法或全局频率比对法;并通过目标中心频率对射频系统的中心频率进行调整。
结合图5B所示内容,另一种射频系统的频率确定方法如下步骤:获取磁共振信号;对磁共振信号进行傅里叶变换,得到频率谱线信息;对频率谱线信息进行归一化处理,基于归一化处理结果更新频率谱线信息;对第一小波基进行修正处理,基于修正处理结果更新第一小波基;对第一小波基进行平移处理,得到第二小波基;对第二小波基与第一小波基进行加权相加,得到第三小波基;基于第三小波基对频率谱线信息进行小波变换,得到第三结果;在第三结果中确定出目标特征峰,确定方法可以是全局最大法或全局频率比对法,从而基于目标特征峰对应的频率确定目标中心频率;并通过目标中心频率对射频系统的中心频率进行调整。
其中,在上述两种射频系统的频率确定方法中,均可以在不同的预设变换尺度和/或不同的平移频率间隔下重复计算第三结果,从而在多个第三结果中确定目标中心频率。
在上述可选的两种射频系统的频率确定方法中,均实现了通过小波变换、平移、平移后相加的方式,综合了频率谱线信息中各特征峰以及与各特征峰相邻的特征峰信息,确定出射频系统的中心频率,提高了中心频率的准确性。并且,该射频系统的频率确定方法适用于各种场强,在不同的场强下均可基于该方法准确提取出射频系统的中心频率。
图6为本发明另一实施例提供的一种射频系统的频率确定装置的结构示意图,本实施例可适用于基于扫描对象的频率谱线信息提取射频系统的中心频率的情况,该装置具体包括:频率谱线获取模块610、第一结果确定模块620、第三结果确定模块630以及中心频率确定模块640。
频率谱线获取模块610,用于获取扫描对象的频率谱线信息;
第一结果确定模块620,用于基于预设变换尺度对所述频率谱线信息进行小波变换,得到所述频率谱线信息的第一结果;
第三结果确定模块630,用于基于平移频率间距以及所述第一结果确定第二结果,对所述第一结果和所述第二结果进行加权处理,得到第三结果;
中心频率确定模块640,用于基于所述第三结果确定射频系统对应的目标中心频率。
可选的,所述射频系统的频率确定装置还包括修正处理模块,用于在所述基于平移频率间距以及所述第一结果确定第二结果之前,对所述第一结果中的正值部分和/或负值部分进行修正处理;基于修正处理结果更新所述第一结果。
可选的,第一结果确定模块具体用于:
对进行小波变换的小波基中的正值部分和/或负值部分进行修正处理,基于修正处理结果更新所述小波基;
基于所述预设变换尺度和所述小波基对所述频率谱线信息进行小波变换,得到所述频率谱线信息的第一结果。
可选的,中心频率确定模块包括中心频率筛选单元,用于若所述第三结果的数量为至少两个,则确定各所述第三结果对应的候选中心频率;基于各所述候选中心频率确定射频系统对应的目标中心频率。
可选的,第三结果确定模块630具体用于基于所述第一结果以及至少两个平移频率间距确定至少两个第二结果,对所述第一结果和至少两个所述第二结果进行加权处理,得到至少两个第三结果。
可选的,第一结果确定模块具体用于基于至少两个预设变换尺度对所述频率谱线信息进行小波变换,得到所述频率谱线信息的至少两个第一结果;相应的,第三结果确定模块630具体用于基于平移频率间距以及至少两个所述第一结果确定至少两个第二结果,对所述第一结果和至少两个所述第二结果进行加权处理,得到至少两个第三结果。
可选的,所述中心频率确定模块640具体用于若所述第三结果包括多个特征峰,则通过全局最大法或全局频率比对法,在第三结果中确定目标特征峰;基于所述目标特征峰对应的频率确定射频系统对应的目标中心频率。
可选的,所述射频系统的频率确定装置还包括归一化模块,用于在所述基于预设变换尺度对所述频率谱线信息进行小波变换之前,对所述频率谱线信息进行归一化处理,基于归一化处理结果更新频率谱线信息;相应的,中心频率确定模块640具体用于基于第三结果确定射频系统对应的归一化中心频率;基于归一化中心频率确定射频系统对应的目标中心频率。
在本实施例中,通过第一结果确定模块,基于预设变换尺度对频率谱线信息进行小波变换,得到频率谱线信息的小波变换的第一结果,并通过第三结果确定模块,基于平移频率间距以及第一结果,确定出频率谱线信息的小波变换平移的第二结果,进而对第一结果和第二结果进行加权处理,得到小波变换平移相加的第三结果,实现了通过小波变换以及平移相加的方式,综合了频率谱线信息中各特征峰以及与各特征峰相邻的特征峰信息,确定出射频系统的中心频率,提高了中心频率的准确性。并且,该射频系统的频率确定方法适用于各种场强,在不同的场强下均可基于该方法准确提取出射频系统的中心频率。
本发明实施例所提供的射频系统的频率确定装置可执行本发明任意实施例所提供的射频系统的频率确定方法,具备执行方法相应的功能模块和有益效果。
值得注意的是,上述系统所包括的各个单元和模块只是按照功能逻辑进行划分的,但并不局限于上述的划分,只要能够实现相应的功能即可;另外,各功能单元的具体名称也只是为了便于相互区分,并不用于限制本发明实施例的保护范围。
在另一个实施例中,提供了一种磁共振设备,包括处理器,所述处理器用于在所述磁共振设备对扫描对象进行成像扫描前,基于本发明任意实施例所提供的射频系统的频率确定方法,确定所述磁共振设备的射频系统的目标中心频率。
本实施例提供的磁共振设备,通过处理器,基于预设变换尺度对频率谱线信息进行小波变换,得到频率谱线信息的小波变换的第一结果,并基于平移频率间距以及第一结果,确定出频率谱线信息的小波变换平移的第二结果,进而对第一结果和第二结果进行加权处理,得到小波变换平移相加的第三结果,实现了通过小波变换以及平移相加的方式,综合了频率谱线信息中各特征峰以及与各特征峰相邻的特征峰信息,确定出射频系统的中心频率,提高了中心频率的准确性。并且,该射频系统的频率确定方法适用于各种场强,在不同的场强下均可基于该方法准确提取出射频系统的中心频率。
图7是本发明另一实施例提供的一种电子设备的结构示意图。图7示出了适于用来实现本发明实施方式的示例性电子设备12的框图。图7显示的电子设备12仅仅是一个示例,不应对本发明实施例的功能和使用范围带来任何限制。设备12典型的是承担确定射频系统的中心频率的功能的电子设备。
如图7所示,电子设备12以通用计算设备的形式表现。电子设备12的组件可以包括但不限于:一个或者多个处理器或者处理单元16,存储器28,连接不同组件(包括存储器28和处理单元16)的总线18。
总线18表示几类总线结构中的一种或多种,包括存储器总线或者存储器控制器,外围总线,图形加速端口,处理器或者使用多种总线结构中的任意总线结构的局域总线。举例来说,这些体系结构包括但不限于工业标准体系结构(Industry StandardArchitecture,ISA)总线,微通道体系结构(Micro Channel Architecture,MCA)总线,增强型ISA总线、视频电子标准协会(Video Electronics Standards Association,VESA)局域总线以及外围组件互连(Peripheral Component Interconnect,PCI)总线。
电子设备12典型地包括多种计算机可读介质。这些介质可以是任何能够被电子设备12访问的可用介质,包括易失性和非易失性介质,可移动的和不可移动的介质。
存储器28可以包括易失性存储器形式的计算机装置可读介质,例如随机存取存储器(Random Access Memory,RAM)30和/或高速缓存存储器32。电子设备12可以进一步包括其它可移动/不可移动的、易失性/非易失性计算机存储介质。仅作为举例,存储装置34可以用于读写不可移动的、非易失性磁介质(图7未显示,通常称为“硬盘驱动器”)。尽管图7中未示出,可以提供用于对可移动非易失性磁盘(例如“软盘”)读写的磁盘驱动器,以及对可移动非易失性光盘(例如只读光盘(Compact Disc-Read Only Memory,CD-ROM)、数字视盘(Digital Video Disc-Read Only Memory,DVD-ROM)或者其它光介质)读写的光盘驱动器。在这些情况下,每个驱动器可以通过一个或者多个数据介质接口与总线18相连。存储器28可以包括至少一个程序产品40,该程序产品40具有一组程序模块42,这些程序模块被配置以执行本发明各实施例的功能。程序产品40,可以存储在例如存储器28中,这样的程序模块42包括但不限于一个或者多个应用程序、其它程序模块以及程序数据,这些示例中的每一个或某种组合中可能包括网络环境的实现。程序模块42通常执行本发明所描述的实施例中的功能和/或方法。
电子设备12也可以与一个或多个外部设备14(例如键盘、鼠标、摄像头等和显示器)通信,还可与一个或者多个使得用户能与该电子设备12交互的设备通信,和/或与使得该电子设备12能与一个或多个其它计算设备进行通信的任何设备(例如网卡,调制解调器等等)通信。这种通信可以通过输入/输出(I/O)接口22进行。并且,电子设备12还可以通过网络适配器20与一个或者多个网络(例如局域网(Local Area Network,LAN),广域网WideArea Network,WAN)和/或公共网络,例如因特网)通信。如图所示,网络适配器20通过总线18与电子设备12的其它模块通信。应当明白,尽管图中未示出,可以结合电子设备12使用其它硬件和/或软件模块,包括但不限于:微代码、设备驱动器、冗余处理单元、外部磁盘驱动阵列、磁盘阵列(Redundant Arrays of Independent Disks,RAID)装置、磁带驱动器以及数据备份存储装置等。
处理器16通过运行存储在存储器28中的程序,从而执行各种功能应用以及数据处理,例如实现本发明上述实施例所提供的射频系统的频率确定方法,包括:
获取扫描对象的频率谱线信息;
基于预设变换尺度对所述频率谱线信息进行小波变换,得到所述频率谱线信息的第一结果;
基于平移频率间距以及所述第一结果确定第二结果,对所述第一结果和所述第二结果进行加权处理,得到第三结果;
基于所述第三结果确定射频系统对应的目标中心频率。
当然,本领域技术人员可以理解,处理器还可以实现本发明任意实施例所提供的射频系统的频率确定方法的技术方案。
本发明实施例还提供一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现如本发明任意实施例所提供的射频系统的频率确定方法步骤,该方法包括:
获取扫描对象的频率谱线信息;
基于预设变换尺度对所述频率谱线信息进行小波变换,得到所述频率谱线信息的第一结果;
基于平移频率间距以及所述第一结果确定第二结果,对所述第一结果和所述第二结果进行加权处理,得到第三结果;
基于所述第三结果确定射频系统对应的目标中心频率。
本发明实施例的计算机存储介质,可以采用一个或多个计算机可读的介质的任意组合。计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件,或者任意以上的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括:具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑磁盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本文件中,计算机可读存储介质可以是任何包含或存储程序的有形介质,该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。
计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号,其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式,包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质,该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。
计算机可读介质上包含的程序代码可以用任何适当的介质传输,包括——但不限于无线、电线、光缆、RF等等,或者上述的任意合适的组合。
可以以一种或多种程序设计语言或其组合来编写用于执行本发明实施例操作的计算机程序代码,所述程序设计语言包括面向对象的程序设计语言—诸如Java、Smalltalk、C++,还包括常规的过程式程序设计语言——诸如“C”语言或类似的程序设计语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中,远程计算机可以通过任意种类的网络——包括局域网(LAN)或广域网(WAN)—连接到用户计算机,或者,可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
注意,上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解,本发明不限于这里所述的特定实施例,对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此,虽然通过以上实施例对本发明进行了较为详细的说明,但是本发明不仅仅限于以上实施例,在不脱离本发明构思的情况下,还可以包括更多其他等效实施例,而本发明的范围由所附的权利要求范围决定。
Claims (10)
1.一种射频系统的频率确定方法,其特征在于,包括:
获取扫描对象的频率谱线信息;
基于预设变换尺度对所述频率谱线信息进行小波变换,得到所述频率谱线信息的第一结果;
基于平移频率间距以及所述第一结果确定第二结果,对所述第一结果和所述第二结果进行加权处理,得到第三结果;
基于所述第三结果确定射频系统对应的目标中心频率。
2.根据权利要求1所述的方法,其特征在于,在所述基于平移频率间距以及所述第一结果确定第二结果之前,还包括:
对所述第一结果中的正值部分和/或负值部分进行修正处理;
基于修正处理结果更新所述第一结果。
3.根据权利要求1所述的方法,其特征在于,所述基于所述第三结果确定射频系统对应的目标中心频率,包括:
若所述第三结果的数量为至少两个,则确定各所述第三结果对应的候选中心频率;
基于各所述候选中心频率确定射频系统对应的目标中心频率。
4.根据权利要求3所述的方法,其特征在于,所述基于平移频率间距以及所述第一结果确定第二结果,包括:
基于所述第一结果以及至少两个平移频率间距确定至少两个第二结果;
相应的,所述对所述第一结果和所述第二结果进行加权处理,得到第三结果,包括:
对所述第一结果和至少两个所述第二结果进行加权处理,得到至少两个第三结果。
5.根据权利要求3所述的方法,其特征在于,所述基于预设变换尺度对所述频率谱线信息进行小波变换,得到所述频率谱线信息的第一结果,包括:
基于至少两个预设变换尺度对所述频率谱线信息进行小波变换,得到所述频率谱线信息的至少两个第一结果;
相应的,所述基于平移频率间距以及所述第一结果确定第二结果,对所述第一结果和所述第二结果进行加权处理,得到第三结果,包括:
基于平移频率间距以及至少两个所述第一结果确定至少两个第二结果,对所述第一结果和至少两个所述第二结果进行加权处理,得到至少两个第三结果。
6.根据权利要求1所述的方法,其特征在于,在所述基于预设变换尺度对所述频率谱线信息进行小波变换之前,还包括:
对所述频率谱线信息进行归一化处理,基于归一化处理结果更新频率谱线信息;
相应的,所述基于所述第三结果确定射频系统对应的目标中心频率,包括:
基于第三结果确定射频系统对应的归一化中心频率,基于归一化中心频率确定射频系统对应的目标中心频率。
7.一种射频系统的频率确定方法,其特征在于,包括:
获取扫描对象的频率谱线信息;
基于平移频率间距对第一小波基进行平移处理,得到第二小波基,对所述第一小波基和所述第二小波基进行加权处理,得到第三小波基;
基于预设变换尺度和所述第三小波基对所述频率谱线信息进行小波变换,得到所述频率谱线信息的第三结果;
基于所述第三结果确定射频系统对应的目标中心频率。
8.一种射频系统的频率确定装置,其特征在于,包括:
频率谱线获取模块,用于获取扫描对象的频率谱线信息;
第一结果确定模块,用于基于预设变换尺度对所述频率谱线信息进行小波变换,得到所述频率谱线信息的第一结果;
第三结果确定模块,用于基于平移频率间距以及所述第一结果确定第二结果,对所述第一结果和所述第二结果进行加权处理,得到第三结果;
中心频率确定模块,用于基于所述第三结果确定射频系统对应的目标中心频率。
9.一种磁共振设备,其特征在于,包括:
处理器,所述处理器用于在所述磁共振设备对扫描对象进行成像扫描前,基于权利要求1-7中任一所述的射频系统的频率确定方法,确定所述磁共振设备的射频系统的目标中心频率。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现如权利要求1-7中任一所述的射频系统的频率确定方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110835685.1A CN113589208B (zh) | 2021-07-23 | 射频系统的频率确定方法、装置、磁共振设备及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110835685.1A CN113589208B (zh) | 2021-07-23 | 射频系统的频率确定方法、装置、磁共振设备及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113589208A true CN113589208A (zh) | 2021-11-02 |
CN113589208B CN113589208B (zh) | 2024-05-31 |
Family
ID=
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10258974A (ja) * | 1996-09-13 | 1998-09-29 | Toshiba Corp | 非定常信号解析装置及び非定常信号解析プログラムを記録した媒体 |
JPH1139287A (ja) * | 1997-07-22 | 1999-02-12 | Chiyouko Saitou | 離散ウェーブレット変換を用いたシステムパラメータ 同定方式 |
JP2000285104A (ja) * | 1999-01-28 | 2000-10-13 | Atr Ningen Joho Tsushin Kenkyusho:Kk | 信号処理方法および装置 |
US20100016693A1 (en) * | 2008-07-15 | 2010-01-21 | Nellcor Puritan Bennett Ireland | Systems And Methods For Identifying Pulse Rates |
CN102706555A (zh) * | 2012-01-05 | 2012-10-03 | 重庆大学 | 一种复解析最优小波解调法 |
CN107229597A (zh) * | 2017-05-31 | 2017-10-03 | 成都理工大学 | 同步挤压广义s变换信号时频分解与重构方法 |
CN109725274A (zh) * | 2018-12-30 | 2019-05-07 | 上海联影医疗科技有限公司 | 磁共振波谱扫描以及其扫描调整方法、装置、设备和存储介质 |
KR102079034B1 (ko) * | 2019-01-30 | 2020-02-19 | 윤상흠 | 전류 공간 벡터와 푸리에 또는 웨이블릿 변환을 이용하여 유도 전동기의 회전자 결함 및 부하 결함을 진단하는 방법 및 장치 |
CN111141520A (zh) * | 2020-02-24 | 2020-05-12 | 江南大学 | 一种基于改进经验小波变换的滚动轴承故障诊断方法 |
CN111259323A (zh) * | 2020-02-22 | 2020-06-09 | 西安科技大学 | 一种采煤机摇臂机械传动系统故障精确定位方法 |
CN112861740A (zh) * | 2021-02-18 | 2021-05-28 | 哈尔滨工程大学 | 基于复合评价指标和小波熵的小波阈值去噪参数选取方法 |
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10258974A (ja) * | 1996-09-13 | 1998-09-29 | Toshiba Corp | 非定常信号解析装置及び非定常信号解析プログラムを記録した媒体 |
JPH1139287A (ja) * | 1997-07-22 | 1999-02-12 | Chiyouko Saitou | 離散ウェーブレット変換を用いたシステムパラメータ 同定方式 |
JP2000285104A (ja) * | 1999-01-28 | 2000-10-13 | Atr Ningen Joho Tsushin Kenkyusho:Kk | 信号処理方法および装置 |
US20100016693A1 (en) * | 2008-07-15 | 2010-01-21 | Nellcor Puritan Bennett Ireland | Systems And Methods For Identifying Pulse Rates |
CN102706555A (zh) * | 2012-01-05 | 2012-10-03 | 重庆大学 | 一种复解析最优小波解调法 |
CN107229597A (zh) * | 2017-05-31 | 2017-10-03 | 成都理工大学 | 同步挤压广义s变换信号时频分解与重构方法 |
CN109725274A (zh) * | 2018-12-30 | 2019-05-07 | 上海联影医疗科技有限公司 | 磁共振波谱扫描以及其扫描调整方法、装置、设备和存储介质 |
KR102079034B1 (ko) * | 2019-01-30 | 2020-02-19 | 윤상흠 | 전류 공간 벡터와 푸리에 또는 웨이블릿 변환을 이용하여 유도 전동기의 회전자 결함 및 부하 결함을 진단하는 방법 및 장치 |
CN111259323A (zh) * | 2020-02-22 | 2020-06-09 | 西安科技大学 | 一种采煤机摇臂机械传动系统故障精确定位方法 |
CN111141520A (zh) * | 2020-02-24 | 2020-05-12 | 江南大学 | 一种基于改进经验小波变换的滚动轴承故障诊断方法 |
CN112861740A (zh) * | 2021-02-18 | 2021-05-28 | 哈尔滨工程大学 | 基于复合评价指标和小波熵的小波阈值去噪参数选取方法 |
Non-Patent Citations (2)
Title |
---|
D. D. ARIANANDA 等: "A Wavelet Packet Transceiver for Spectral Analysis and Dynamic Spectrum Access", 《2011 IEEE 73RD VEHICULAR TECHNOLOGY CONFERENCE (VTC SPRING)》, 18 July 2011 (2011-07-18) * |
蔡笑风 等: "基于经验小波变换干耦合超声检测Lamb波信号分析", 《中国测试》, vol. 45, no. 1 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109709503B (zh) | 一种磁共振系统控制方法、磁共振系统、计算机设备 | |
CN107656224B (zh) | 一种磁共振成像方法、装置及系统 | |
CN113298831A (zh) | 图像分割方法、装置、电子设备及存储介质 | |
CN114167334B (zh) | 磁共振图像的重建方法、装置及电子设备 | |
CN114399483A (zh) | 皮层靶点确定方法、装置、电子设备及存储介质 | |
CN114944229A (zh) | 基于深度学习和磁共振结构脑影像的脑龄预测方法 | |
US10620287B2 (en) | Method and apparatus for denoising magnetic resonance diffusion tensor, and computer program product | |
CN111445546A (zh) | 图像重建方法、装置、电子设备、存储介质 | |
US20110101978A1 (en) | Motion information capture and automatic motion correction for imaging systems | |
US9590711B2 (en) | Managing beamformed signals to optimize transmission rates of sensor arrays | |
CN113589208B (zh) | 射频系统的频率确定方法、装置、磁共振设备及存储介质 | |
CN113589208A (zh) | 射频系统的频率确定方法、装置、磁共振设备及存储介质 | |
CN112581385A (zh) | 基于多先验约束扩散峰度成像张量估计方法、介质和设备 | |
CN116484184A (zh) | 一种电力设备局部放电缺陷样本增强方法及装置 | |
CN112101396A (zh) | 一种分类方法、装置、设备和存储介质 | |
Pieciak et al. | Variance stabilization of noncentral-chi data: Application to noise estimation in MRI | |
CN111239657B (zh) | 谱图的相位校正方法、装置及设备 | |
Chang et al. | Instrument variables for reducing noise in parallel MRI reconstruction | |
CN114494014A (zh) | 一种磁共振图像超分辨率重建方法及装置 | |
CN114897032A (zh) | 一种基于宽度学习的电流互感器故障诊断方法、装置 | |
CN115153476A (zh) | 基于多维数据的睡眠评估方法、装置、电子设备及介质 | |
CN114972565A (zh) | 图像的降噪方法、装置、电子设备及介质 | |
CN111612865A (zh) | 一种基于条件生成对抗网络的mri成像方法及装置 | |
CN113702666A (zh) | 一种光纤陀螺惯性测量单元信号联合降噪方法 | |
CN112800831A (zh) | 一种用于时变重力场的emd滤波方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant |