CN117278930A - Hoa系数的获取方法和装置 - Google Patents

Hoa系数的获取方法和装置 Download PDF

Info

Publication number
CN117278930A
CN117278930A CN202310818175.2A CN202310818175A CN117278930A CN 117278930 A CN117278930 A CN 117278930A CN 202310818175 A CN202310818175 A CN 202310818175A CN 117278930 A CN117278930 A CN 117278930A
Authority
CN
China
Prior art keywords
information
reference points
angle information
trigonometric function
pitch angle
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.)
Pending
Application number
CN202310818175.2A
Other languages
English (en)
Inventor
高原
刘帅
王宾
王喆
曲天书
徐佳浩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Peking University
Huawei Technologies Co Ltd
Original Assignee
Peking University
Huawei Technologies Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Peking University, Huawei Technologies Co Ltd filed Critical Peking University
Priority to CN202310818175.2A priority Critical patent/CN117278930A/zh
Publication of CN117278930A publication Critical patent/CN117278930A/zh
Pending legal-status Critical Current

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04S—STEREOPHONIC SYSTEMS 
    • H04S7/00—Indicating arrangements; Control arrangements, e.g. balance control
    • H04S7/30—Control circuits for electronic adaptation of the sound field
    • H04S7/302—Electronic adaptation of stereophonic sound system to listener position or orientation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04S—STEREOPHONIC SYSTEMS 
    • H04S7/00—Indicating arrangements; Control arrangements, e.g. balance control
    • H04S7/30—Control circuits for electronic adaptation of the sound field
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04S—STEREOPHONIC SYSTEMS 
    • H04S5/00—Pseudo-stereo systems, e.g. in which additional channel signals are derived from monophonic signals by means of phase shifting, time delay or reverberation 
    • H04S5/005—Pseudo-stereo systems, e.g. in which additional channel signals are derived from monophonic signals by means of phase shifting, time delay or reverberation  of the pseudo five- or more-channel type, e.g. virtual surround
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
    • H04R5/00—Stereophonic arrangements
    • H04R5/02—Spatial or constructional arrangements of loudspeakers
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
    • H04R5/00—Stereophonic arrangements
    • H04R5/033—Headphones for stereophonic communication
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04S—STEREOPHONIC SYSTEMS 
    • H04S2400/00—Details of stereophonic systems covered by H04S but not provided for in its groups
    • H04S2400/11—Positioning of individual sound objects, e.g. moving airplane, within a sound field
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04S—STEREOPHONIC SYSTEMS 
    • H04S2420/00—Techniques used stereophonic systems covered by H04S but not provided for in its groups
    • H04S2420/11—Application of ambisonics in stereophonic audio systems
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Signal Processing (AREA)
  • Image Analysis (AREA)
  • Stereophonic System (AREA)
  • Complex Calculations (AREA)
  • Audible-Bandwidth Dynamoelectric Transducers Other Than Pickups (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本申请提供一种HOA系数的获取方法和装置。HOA系数的获取方法包括：获取虚拟扬声器在预设球面上的位置信息；预设球面上包含M个经线圈和N个纬线圈；根据位置信息以及预设的参考三角函数表获取位置信息对应的三角函数值，参考三角函数表包括俯仰三角函数表和/或水平三角函数表；俯仰三角函数表包括第一经线圈上的多个第一参考点的俯仰角索引对应的多个三角函数值，第一经线圈是M个经线圈的其中之一；水平三角函数表包括第一纬线圈上的多个第二参考点的水平角索引对应的多个三角函数值，第一纬线圈是N个纬线圈的其中之一；根据位置信息对应的三角函数值获取虚拟扬声器的HOA系数。本申请可以减少预先阶段的计算量，节省存储空间。

Description

HOA系数的获取方法和装置
本申请是分案申请，本申请要求于2021年3月5日提交的申请号为202110246382.6、发明名称为“HOA系数的获取方法和装置”的中国专利申请的优先权，其全部内容通过引用结合在本申请中。
技术领域
本申请涉及音频处理技术，尤其涉及一种HOA系数的获取方法和装置。
背景技术
三维音频技术是通过计算机、信号处理等方式对真实世界中声音事件和三维声场信息进行获取、处理，传输和渲染回放的音频技术。三维音频技术使声音具有强烈的空间感、包围感及沉浸感，给人以“声临其境”的听觉体验。目前主流的三维音频技术是高阶立体混响(higher order ambisonics，HOA)技术，HOA技术因其在录制和编码中与回放阶段的扬声器布局无关的性质，以及HOA格式数据的可旋转特性，使得HOA技术在三维音频回放时具有更高的灵活性，因而也得到了更为广泛的关注和研究。
HOA技术是将HOA信号转为实际扬声器信号后进行回放，或者是将HOA信号转为虚拟扬声器信号再映射为双耳信号进行回放。在上述过程中，计算HOA系数是关键步骤。相关技术中是根据P阶球谐函数展开公式计算HOA系数。
但是上述方法计算量大，且存储空间需求高。
发明内容
本申请提供一种HOA系数的获取方法和装置，以减少预先阶段的计算量，节省存储空间。
第一方面，本申请提供一种HOA系数的获取方法，包括：获取虚拟扬声器在预设球面上的位置信息，所述位置信息包括俯仰角信息和/或水平角信息；所述预设球面上包含M个经线圈和N个纬线圈，所述经线圈和所述纬线圈的交叉点称作参考点，位于同一个经线圈上的N个参考点等间距排列，位于同一个纬线圈上的M个参考点等间距排列；根据所述位置信息以及预设的参考三角函数表获取所述位置信息对应的三角函数值，所述参考三角函数表包括俯仰三角函数表和/或水平三角函数表；所述俯仰三角函数表包括多个第一参考点的俯仰角信息对应的多个三角函数值，所述多个第一参考点是第一经线圈上的参考点，所述多个第一参考点的数量不小于 表示向下取整，所述第一经线圈是所述M个经线圈的其中之一；所述水平三角函数表包括多个第二参考点的水平角信息对应的多个三角函数值，所述多个第二参考点是第一纬线圈上的参考点，所述多个第二参考点的数量不小于所述第一纬线圈是所述N个纬线圈的其中之一；根据所述位置信息对应的三角函数值获取所述虚拟扬声器的高阶立体混响HOA系数。
本申请中，预先只需要获取预设球面上的某一经线圈上的所有参考点中的四分之一个参考点的三角函数值(俯仰三角函数表)，和/或，预设球面上的某一纬线圈上的所有参考点中的四分之一个参考点的三角函数值(水平三角函数表)，可以基于前述三角函数表获取到0°～360°的任意角度的三角函数值。这样一方面可以减少预先阶段的计算量，另一方面由于三角函数表中的角度数量减少，从而需要存储的三角函数表的数据量也相应的减少，可以节省存储空间。
预设球面上包含M个经线圈和N个纬线圈，经线圈和纬线圈的交叉点称作参考点，位于同一个经线圈上的N个参考点等间距排列，位于同一个纬线圈上的M个参考点等间距排列。为了保证参考点可以尽可能覆盖到所有虚拟扬声器，因此可以满足M了保证参，其中K为虚拟扬声器总数，例如K可以取值为2048，1024，530等。本申请中，两两相邻的经线圈之间的水平角信息是相等的，即M个经线圈在球面上是均匀分布的，因此位于同一个经线圈上的N个参考点等间距排列；两两相邻的纬线圈之间的俯仰角信息是相等的，即N个纬线圈在球面上是均匀分布的，因此位于同一个纬线圈上的M个参考点等间距排列。相应的，由于位于同一个经线圈上的N个参考点等间距排列，因此位于同一个经线圈上的两两相邻的参考点的俯仰角信息差相等；由于位于同一个经线圈上的N个参考点等间距排列，因此位于同一个经线圈上的两两相邻的参考点的水平角信息差相等。
上述虚拟扬声器可以是针对声音渲染设置的任意一个虚拟扬声器。虚拟扬声器在上述球面上的位置信息包括俯仰角信息(elevation)和/或水平角信息。虚拟扬声器设置的位置和球心之间的连线与水平面(例如赤道圈所在平面)之间的夹角为虚拟扬声器的位置信息的俯仰角信息(或称作高度角信息)；虚拟扬声器设置的位置和球心之间的连线在水平面上的投影与设定初始方向的夹角为虚拟扬声器的位置信息的水平角信息。可选的，本申请中的水平角信息也可以被称作方位角(azimuth)。
参考三角函数表包括俯仰三角函数表和/或水平三角函数表。其中，
俯仰三角函数表包括多个第一参考点的俯仰角信息对应的多个三角函数值，多个第一参考点是第一经线圈上的参考点，多个第一参考点的数量不小于表示向下取整，第一经线圈是M个经线圈的其中之一。
本申请从M个经线圈中任意选取一条视作第一经线圈，用于计算俯仰三角函数表的参考点(即上述多个第一参考点)是从该第一经线圈上的N个参考点中选取的。基于三角函数的原理，可以将0°(0)～360°(2π)的角度范围划分为四个子范围，即0°(0)～90°(π/2)，90°(π/2)～180°(π)，180°(π)～270°(3π/2)以及270°(3π/2)～360°(2π)，通过该四个子范围的任意一个子范围计算得到的三角函数值可以扩展到0°(0)～360°(2π)的角度范围的三角函数值。因此本申请在选取上述多个第一参考点时，其数量不小于即从第一经线圈上的N个参考点中选取至少四分之一个参考点作为上述多个第一参考点。
可选的，上述多个第一参考点的俯仰角信息可以属于上述四个子范围的其中之一，例如，多个第一参考点的俯仰角信息在0°(0)～90°(π/2)内；或者，多个第一参考点的俯仰角信息在90°(π/2)～180°(π)内；或者，多个第一参考点的俯仰角信息在180°(π)～270°(3π/2)内；或者，多个第一参考点的俯仰角信息在270°(3π/2)～360°(2π)内。
可选的，上述多个第一参考点的俯仰角信息可以属于上述四个子范围的其中至少两个，例如，多个第一参考点中的第一部分参考点的俯仰角信息在0°(0)～90°(π/2)内，第二部分参考点的俯仰角信息在90°(π/2)～180°(π)内，这两部分参考点的俯仰角信息在计算三角函数值时不会重叠；或者，多个第一参考点中的部分参考点的俯仰角信息在90°(π/2)～180°(π)内，另一部分参考点的俯仰角信息在180°(π)～270°(3π/2)内，这两部分参考点的俯仰角信息在计算三角函数值时不会重叠；或者，多个第一参考点中的部分参考点的俯仰角信息在180°(π)～270°(3π/2)内，另一部分参考点的俯仰角信息在270°(3π/2)～360°(2π)内，这两部分参考点的俯仰角信息在计算三角函数值时不会重叠；或者，多个第一参考点中的第一部分参考点的俯仰角信息在0°(0)～90°(π/2)内，第二部分参考点的俯仰角信息在90°(π/2)～180°(π)内，第三部分参考点的俯仰角信息在180°(π)～270°(3π/2)内，这三部分参考点的俯仰角信息在计算三角函数值时不会重叠；或者，多个第一参考点中的第一部分参考点的俯仰角信息在90°(π/2)～180°(π)内，第二部分参考点的俯仰角信息在180°(π)～270°(3π/2)内，第三部分参考点的俯仰角信息在270°(3π/2)～360°(2π)内，这三部分参考点的俯仰角信息在计算三角函数值时不会重叠；或者，多个第一参考点中的第一部分参考点的俯仰角信息在0°(0)～90°(π/2)内，第二部分参考点的俯仰角信息在90°(π/2)～180°(π)内，第三部分参考点的俯仰角信息在180°(π)～270°(3π/2)内，第四部分参考点的俯仰角信息在270°(3π/2)～360°(2π)内，这四部分参考点的俯仰角信息在计算三角函数值时不会重叠。
上述不会重叠可以是指属于不同子范围的俯仰角信息计算得到的三角函数值不相等，例如，多个第一参考点共有四个(n1、n2、n3和n4)，其中，n1的俯仰角信息为0，n2的俯仰角信息为π/6，n3的俯仰角信息为π/3，n4的俯仰角信息为π/2；或者，n1的俯仰角信息为0，n2的俯仰角信息为2/3π，n3的俯仰角信息为4/3π，n4的俯仰角信息为2π。需要说明的是，本申请中对多个第一参考点的俯仰角信息所述的子范围不做具体限定，只要该多个第一参考点的俯仰角信息的三角函数值可以覆盖其中一个子范围的全部角度的三角函数值即可。
相应的，上述位置信息的俯仰角信息取值范围为0～2π。
本申请中多个第一参考点的俯仰角信息对应的多个三角函数可以是正弦函数值，也可以是余弦函数值。
可选的，第i个第一参考点的俯仰角信息对应的正弦函数值满足如下公式(1)：
其中，i＝0,1,…,N’，ri表示上述第一经线圈的半径。
可选的，第i个第一参考点的俯仰角信息对应的余弦函数值满足如下公式(2)：
其中，i＝0,1,…,N’，ri表示上述第一经线圈的半径。
水平三角函数表包括多个第二参考点的水平角信息对应的多个三角函数值，多个第二参考点是第一纬线圈上的参考点，多个第二参考点的数量不小于表示向下取整，第一纬线圈是N个纬线圈的其中之一。
本申请从N个纬线圈中任意选取一条视作第一纬线圈，用于计算水平三角函数表的参考点(即上述多个第二参考点)是从该第一纬线圈上的M个参考点中选取的。同样基于三角函数的原理，可以将0°(0)～360°(2π)的角度范围划分为四个子范围，即0°(0)～90°(π/2)，90°(π/2)～180°(π)，180°(π)～270°(3π/2)以及270°(3π/2)～360°(2π)，通过该四个子范围的任意一个子范围计算得到的三角函数值可以扩展到0°(0)～360°(2π)的角度范围的三角函数值。因此本申请在选取上述多个第二参考点时，其数量不小于即从第一纬线圈上的M个参考点中选取至少四分之一个参考点作为上述多个第二参考点。
可选的，上述多个第二参考点的水平角信息可以属于上述四个子范围的其中之一，例如，多个第二参考点的水平角信息在0°(0)～90°(π/2)内；或者，多个第二参考点的水平角信息在90°(π/2)～180°(π)内；或者，多个第二参考点的水平角信息在180°(π)～270°(3π/2)内；或者，多个第二参考点的水平角信息在270°(3π/2)～360°(2π)内。
可选的，上述多个第二参考点的水平角信息可以属于上述四个子范围的其中至少两个，例如，多个第二参考点中的第一部分参考点的水平角信息在0°(0)～90°(π/2)内，第二部分参考点的水平角信息在90°(π/2)～180°(π)内，这两部分参考点的水平角信息在计算三角函数值时不会重叠；或者，多个第二参考点中的第一部分参考点的水平角信息在90°(π/2)～180°(π)内，第二部分参考点的水平角信息在180°(π)～270°(3π/2)内，这两部分参考点的水平角信息在计算三角函数值时不会重叠；或者，多个第二参考点中的第一部分参考点的水平角信息在180°(π)～270°(3π/2)内，第二部分参考点的水平角信息在270°(3π/2)～360°(2π)内，这两部分参考点的水平角信息在计算三角函数值时不会重叠；或者，多个第二参考点中的第一部分参考点的水平角信息在0°(0)～90°(π/2)内，第二部分参考点的水平角信息在90°(π/2)～180°(π)内，第三部分参考点的水平角信息在180°(π)～270°(3π/2)内，这三部分参考点的水平角信息在计算三角函数值时不会重叠；或者，多个第二参考点中的第一部分参考点的水平角信息在90°(π/2)～180°(π)内，第二部分参考点的水平角信息在180°(π)～270°(3π/2)内，第三部分参考点的水平角信息在270°(3π/2)～360°(2π)内，这三部分参考点的水平角信息在计算三角函数值时不会重叠；或者，多个第二参考点中的第一部分参考点的水平角信息在0°(0)～90°(π/2)内，第二部分参考点的水平角信息在90°(π/2)～180°(π)内，第三部分参考点的水平角信息在180°(π)～270°(3π/2)内，第四部分参考点的水平角信息在270°(3π/2)～360°(2π)内，这四部分参考点的水平角信息在计算三角函数值时不会重叠。
上述不会重叠可以是指属于不同子范围的水平角信息计算得到的三角函数值不相等，例如，多个第二参考点共有四个(n1、n2、n3和n4)，其中，n1的水平角信息为0，n2的水平角信息为π/6，n3的水平角信息为π/3，n4的水平角信息为π/2；或者，n1的水平角信息为0，n2的水平角信息为2/3π，n3的水平角信息为4/3π，n4的水平角信息为2π。需要说明的是，本申请中对多个第二参考点的水平角信息所述的子范围不做具体限定，只要该多个第二参考点的水平角信息的三角函数值可以覆盖其中一个子范围的全部角度的三角函数值即可。
相应的，上述位置信息的水平角信息取值范围为0～2π。
本申请中多个第二参考点的水平角信息对应的多个三角函数可以是正弦函数值，也可以是余弦函数值。
可选的，第i个第二参考点的水平角信息对应的正弦函数值满足如下公式(3)：
其中，j＝0,1,…,M’，rj表示上述第一纬线圈的半径。
可选的，第i个第二参考点的水平角信息对应的余弦函数值满足如下公式(4)：
其中，j＝0,1,…,M’，rj表示上述第一纬线圈的半径。
需要说明的是，上述俯仰三角函数表和水平三角函数表可以是预先计算并存储的，也可以根据M和N的值实时计算并存储。可选的，由于编码器和解码器均采用分帧处理，因此在HOA系数计算过程中，当前音频帧的M与前一音频帧的M取值相等时，处理当前帧时可以不用重新计算水平三角函数表，直接使用前一音频帧的水平三角函数表即可；当前音频帧的N与前一音频帧的N取值相等时，处理当前帧时可以不用重新计算俯仰三角函数表，直接使用前一音频帧的俯仰三角函数表即可。
本申请中可以根据虚拟扬声器在预设球面上的位置信息的俯仰角信息以及上述俯仰三角函数表获取位置信息的俯仰角信息对应的三角函数值；和/或，根据虚拟扬声器在预设球面上的位置信息的水平角信息以及上述水平三角函数表获取位置信息的水平角信息对应的三角函数值。
可见针对虚拟扬声器在预设球面上的位置信息，可以根据其位置信息包含的具体内容，以及预先设置的参考三角函数表的具体内容来获得该位置信息对应的三角函数值，即位置信息包含俯仰角信息，可以根据该俯仰角信息和参考三角函数表(如果有俯仰三角函数表，则用俯仰三角函数表，如果没有俯仰三角函数表，则用水平三角函数表)获得位置信息的俯仰角信息对应的三角函数值；位置信息包含水平角信息，可以根据该水平角信息和参考三角函数表(如果有水平三角函数表，则用水平三角函数表，如果没有水平三角函数表，则用俯仰三角函数表)获得位置信息的水平角信息对应的三角函数值。
本申请可以采用以下几种方法获取位置信息的俯仰角信息对应的三角函数值：
一、当位置信息的俯仰角信息与多个第一参考点中的一个参考点的俯仰角信息相等时，将俯仰三角函数表中前述一个参考点的俯仰角信息对应的三角函数值作为位置信息的俯仰角信息对应的三角函数值。
位置信息的俯仰角信息与多个第一参考点中的一个参考点的俯仰角信息相等表示虚拟扬声器在预设球面上的位置对应的俯仰角信息正好与俯仰三角函数表中的某一个角度相等，即虚拟扬声器在预设球面上的位置与选取出来用于获取俯仰三角函数表的参考点的位置重合。因此可以直接从俯仰三角函数表中查询该俯仰角信息的三角函数值即可。
二、当位置信息的俯仰角信息与N个参考点中除多个第一参考点外的一个参考点的俯仰角信息相等时，将位置信息的俯仰角信息转换为多个第一参考点中的一个参考点的俯仰角信息，将俯仰三角函数表中前述一个参考点的俯仰角信息对应的三角函数值作为位置信息的俯仰角信息对应的三角函数值。
位置信息的俯仰角信息与N个参考点中除多个第一参考点外的一个参考点的俯仰角信息相等表示虚拟扬声器在预设球面上的位置对应的俯仰角信息不属于俯仰三角函数表中的角度，与第一经线圈上的某一其他参考点的俯仰角信息相等，即虚拟扬声器在预设球面上的位置与某一参考点的位置重合，该参考点没有用于获取俯仰三角函数表，但其俯仰角信息与三角函数表中的角度存在三角函数对应关系。因此可以基于三角函数的原理将位置信息的俯仰角信息转换为多个第一参考点中的一个参考点的俯仰角信息，然后查表得到多个第一参考点中的一个参考点的俯仰角信息对应的三角函数值，将多个第一参考点中的一个参考点的俯仰角信息对应的三角函数值作为位置信息的俯仰角信息对应的三角函数值。
三、当位置信息的俯仰角信息与参考点的俯仰角信息均不相等时，根据位置信息的俯仰角信息获得与位置信息对应的参考点的俯仰角信息；当与位置信息对应的参考点的俯仰角信息与多个第一参考点中的一个参考点的俯仰角信息相等时，将俯仰三角函数表中前述一个参考点的俯仰角信息对应的三角函数值作为位置信息的俯仰角信息对应的三角函数值；当与位置信息对应的参考点的俯仰角信息与N个参考点中除多个第一参考点外的一个参考点的俯仰角信息相等时，将与位置信息对应的参考点的俯仰角信息转换为多个第一参考点中的一个参考点的俯仰角信息，将俯仰三角函数表中前述一个参考点的俯仰角信息对应的三角函数值作为位置信息的俯仰角信息对应的三角函数值。
位置信息的俯仰角信息与参考点的俯仰角信息均不相等表示虚拟扬声器在预设球面上的位置既不在经线圈上，也不在纬线圈上。因此先获得与位置信息对应的参考点的俯仰角信息，再参照上述方法一或方法二的方式获取参考点的俯仰角信息对应的三角函数值，将该参考点的俯仰角信息对应的三角函数值作为位置信息的俯仰角信息对应的三角函数值。
可选的，可以通过索引的方式表示位置信息，索引可以包括俯仰角索引和水平角索引。例如，在任意一个纬线圈上，将该未分布其上的其中一个参考点的水平角度设置为0，然后根据预设的水平角度与水平角索引之间的转换公式转换获得对应的水平角索引；由于纬线圈上的任意相邻参考点之间的水平角差值是相等的，因此可以获得该纬线圈上的其他参考点的水平角度，从而根据上述转换公式获得所述其他参考点各自的水平角索引。需要说明的是，本申请对将纬线圈上的哪个参考点的水平角度设置为0不作具体限定。同理，由于在经线圈方向相邻参考点之间的俯仰角差值满足前述的要求，因此在设置了俯仰角度为0的参考点后，就可以获得其他参考点的俯仰角度，基于预设的俯仰角度和俯仰角索引之间的转换公式就可以获得经线圈上所有参考点的俯仰角索引。需要说明的是，本申请对将经线圈上哪个参考点的俯仰角度设置为0不作具体限定，例如可以是位于所述赤道圈上的参考点，或者所述位于所述南极点上的参考点，或者位于所述北极点上的参考点。
可选的，俯仰角度和该俯仰角信息的索引满足如下公式(5)(即俯仰角度和俯仰角索引的转换公式)：
其中，ri表示第一经线圈的半径，round()表示取整。
本申请可以采用以下几种方法获取位置信息的水平角信息对应的三角函数值：
一、当位置信息的水平角信息与多个第二参考点中的一个参考点的水平角信息相等时，将水平三角函数表中前述一个参考点的水平角信息对应的三角函数值作为位置信息的水平角信息对应的三角函数值。
位置信息的水平角信息与多个第二参考点中的一个参考点的水平角信息相等表示虚拟扬声器在预设球面上的位置对应的水平角信息正好与水平三角函数表中的某一个角度相等，即虚拟扬声器在预设球面上的位置与选取出来用于获取水平三角函数表的参考点的位置重合。因此可以直接从水平三角函数表中查询该水平角信息的三角函数值即可。
二、当位置信息的水平角信息与M个参考点中除多个第二参考点外的一个参考点的水平角信息相等时，将位置信息的水平角信息转换为多个第二参考点中的一个参考点的水平角信息，将水平三角函数表中前述一个参考点的水平角信息对应的三角函数值作为位置信息的水平角信息对应的三角函数值。
位置信息的水平角信息与M个参考点中除多个第二参考点外的一个参考点的水平角信息相等表示虚拟扬声器在预设球面上的位置对应的水平角信息不属于水平三角函数表中的角度，与第一经线圈上的某一其他参考点的水平角信息相等，即虚拟扬声器在预设球面上的位置与某一参考点的位置重合，该参考点没有用于获取水平三角函数表，但其水平角信息与三角函数表中的角度存在三角函数对应关系。因此可以基于三角函数的原理将位置信息的水平角信息转换为多个第二参考点中的一个参考点的水平角信息，然后查表得到多个第二参考点中的一个参考点的水平角信息对应的三角函数值，将多个第二参考点中的一个参考点的水平角信息对应的三角函数值作为位置信息的水平角信息对应的三角函数值。
三、当位置信息的水平角信息与参考点的水平角信息均不相等时，根据位置信息的水平角信息获得与位置信息对应的参考点的水平角信息；当与位置信息对应的参考点的水平角信息与多个第二参考点中的一个参考点的水平角信息相等时，将水平三角函数表中前述一个参考点的水平角信息对应的三角函数值作为位置信息的水平角信息对应的三角函数值；当与位置信息对应的参考点的水平角信息与M个参考点中除多个第二参考点外的一个参考点的水平角信息相等时，将与位置信息对应的参考点的水平角信息转换为多个第二参考点中的一个参考点的水平角信息，将水平三角函数表中前述一个参考点的水平角信息对应的三角函数值作为位置信息的水平角信息对应的三角函数值。
位置信息的水平角信息与参考点的水平角信息均不相等表示虚拟扬声器在预设球面上的位置既不在经线圈上，也不在纬线圈上。因此先获得与位置信息对应的参考点的水平角信息，再参照上述方法一或方法二的方式获取参考点的水平角信息对应的三角函数值，将该参考点的水平角信息对应的三角函数值作为位置信息的水平角信息对应的三角函数值。
可选的，水平角度θ和该水平角信息的索引θ’满足如下公式(6)(即水平角度和水平角索引的转换公式)：
其中，rj表示第一纬线圈的半径，round()表示取整。
本申请中，预先只需要获取预设球面上的某一经线圈上的所有参考点中的四分之一个参考点的三角函数值(俯仰三角函数表)，和/或，预设球面上的某一纬线圈上的所有参考点中的四分之一个参考点的三角函数值(水平三角函数表)，可以基于前述三角函数表获取到0°～360°的任意角度的三角函数值。这样一方面可以减少预先阶段的计算量，另一方面由于三角函数表中的角度数量减少，从而需要存储的三角函数表的数据量也相应的减少，可以节省存储空间。
第二方面，本申请提供一种HOA系数的获取装置，包括：获取模块，用于获取虚拟扬声器在预设球面上的位置信息，所述位置信息包括俯仰角信息和/或水平角信息；所述预设球面上包含M个经线圈和N个纬线圈，所述经线圈和所述纬线圈的交叉点称作参考点，位于同一个经线圈上的N个参考点等间距排列，位于同一个纬线圈上的M个参考点等间距排列；计算模块，用于根据所述位置信息以及预设的参考三角函数表获取所述位置信息对应的三角函数值，所述参考三角函数表包括俯仰三角函数表和/或水平三角函数表；所述俯仰三角函数表包括多个第一参考点的俯仰角索引对应的多个三角函数值，所述多个第一参考点是第一经线圈上的参考点，所述多个第一参考点的数量不小于表示向下取整，所述第一经线圈是所述M个经线圈的其中之一；所述水平三角函数表包括多个第二参考点的水平角索引对应的多个三角函数值，所述多个第二参考点是第一纬线圈上的参考点，所述多个第二参考点的数量不小于所述第一纬线圈是所述N个纬线圈的其中之一；根据所述位置信息对应的三角函数值获取所述虚拟扬声器的高阶立体混响HOA系数。
在一种可能的实现方式中，所述计算模块，具体用于根据所述位置信息的俯仰角信息以及所述俯仰三角函数表获取所述位置信息的俯仰角信息对应的三角函数值；和/或，根据所述位置信息的水平角信息以及所述水平三角函数表获取所述位置信息的水平角信息对应的三角函数值。
在一种可能的实现方式中，所述计算模块，具体用于当所述位置信息的俯仰角信息与所述多个第一参考点中的一个参考点的俯仰角索引对应上时，将所述俯仰三角函数表中所述一个参考点的俯仰角索引对应的三角函数值作为所述位置信息的俯仰角信息对应的三角函数值；或者，当所述位置信息的俯仰角信息与所述N个参考点中除所述多个第一参考点外的一个参考点的俯仰角索引对应上时，将所述N个参考点中除所述多个第一参考点外的一个参考点的俯仰角索引转换为所述多个第一参考点中的一个参考点的俯仰角索引，将所述俯仰三角函数表中所述一个参考点的俯仰角索引对应的三角函数值作为所述位置信息的俯仰角信息对应的三角函数值；或者，当所述位置信息的俯仰角信息与所述参考点的俯仰角索引均不对应时，根据所述位置信息的俯仰角信息获得与所述位置信息的俯仰角信息对应的参考点的俯仰角信息；当所述与所述位置信息的俯仰角信息对应的参考点的俯仰角信息与所述多个第一参考点中的一个参考点的俯仰角索引对应上时，将所述俯仰三角函数表中所述一个参考点的俯仰角索引对应的三角函数值作为所述位置信息的俯仰角信息对应的三角函数值；当所述与所述位置信息的俯仰角信息对应的参考点的俯仰角信息与所述N个参考点中除所述多个第一参考点外的一个参考点的俯仰角索引对应上时，将所述N个参考点中除所述多个第一参考点外的一个参考点的俯仰角索引转换为所述多个第一参考点中的一个参考点的俯仰角索引，将所述俯仰三角函数表中所述一个参考点的俯仰角索引对应的三角函数值作为所述位置信息的俯仰角信息对应的三角函数值。
在一种可能的实现方式中，所述计算模块，具体用于当所述位置信息的水平角信息与所述多个第二参考点中的一个参考点的水平角索引对应上时，将所述水平三角函数表中所述一个参考点的水平角索引对应的三角函数值作为所述位置信息的水平角信息对应的三角函数值；或者，当所述位置信息的水平角信息与所述M个参考点中除所述多个第二参考点外的一个参考点的水平角索引对应上时，将所述M个参考点中除所述多个第二参考点外的一个参考点的水平角索引转换为所述多个第二参考点中的一个参考点的水平角索引，将所述水平三角函数表中所述一个参考点的水平角索引对应的三角函数值作为所述位置信息的水平角信息对应的三角函数值；或者，当所述位置信息的水平角信息与所述参考点的水平角索引均不对应时，根据所述位置信息的水平角信息获得与所述位置信息的水平角信息对应的参考点的水平角信息；当所述与所述位置信息的水平角信息对应的参考点的水平角信息与所述多个第二参考点中的一个参考点的水平角索引对应上时，将所述水平三角函数表中所述一个参考点的水平角索引对应的三角函数值作为所述位置信息的水平角信息对应的三角函数值；当所述与所述位置信息的水平角信息对应的参考点的水平角信息与所述M个参考点中除所述多个第二参考点外的一个参考点的水平角索引对应上时，将所述M个参考点中除所述多个第二参考点外的一个参考点的水平角索引转换为所述多个第二参考点中的一个参考点的水平角索引，将所述水平三角函数表中所述一个参考点的水平角索引对应的三角函数值作为所述位置信息的水平角信息对应的三角函数值。
在一种可能的实现方式中，所述位置信息的俯仰角信息包括所述位置信息的俯仰角度或者所述位置信息的俯仰角索引；所述位置信息的水平角信息包括所述位置信息的水平角度或者所述位置信息的水平角索引。
在一种可能的实现方式中，所述三角函数值是指正弦函数值或者余弦函数值。
在一种可能的实现方式中，当所述多个第一参考点的俯仰角索引对应的多个三角函数为正弦函数值时，第i个第一参考点的俯仰角索引对应的正弦函数值满足如下公式(1)：
其中，i＝0,1,…,N’，ri表示所述第一经线圈的半径；
或，
当所述多个第一参考点的俯仰角索引对应的多个三角函数为余弦函数值时，其中，第i个第一参考点的俯仰角索引对应的余弦函数值满足如下公式(2)：
其中，i＝0,1,…,N’，ri表示所述第一经线圈的半径。
在一种可能的实现方式中，当所述多个第二参考点的水平角索引对应的多个三角函数为正弦函数值时，第j个第二参考点的水平角索引对应的正弦函数值满足如下公式(3)：
其中，j＝0,1,…,M’，rj表示所述第一纬线圈的半径；
或，
当所述多个第二参考点的水平角索引对应的多个三角函数为余弦函数值时，第j个第二参考点的水平角索引对应的余弦函数值满足如下公式(4)：
其中，j＝0,1,…,M’，rj表示所述第一纬线圈的半径。
在一种可能的实现方式中，俯仰角度和俯仰角索引满足如下公式(5)：
其中，ri表示所述第一经线圈的半径，round()表示取整。
在一种可能的实现方式中，水平角度θ和水平角索引θ’满足如下公式(6)：
其中，rj表示所述第一纬线圈的半径，round()表示取整。
在以俯仰角索引的形式表示参考点的位置的俯仰三角函数表中，如果假设N＝1024，那么俯仰三角函数表可以包含257个参考点的俯仰角索引分别对应的三角函数值，例如，sin[0]～sin[256]。
同理，在以水平角索引的形式表示参考点的位置的水平三角函数表中，如果假设M＝1024，那么水平三角函数表也可以包含257个参考点的水平角索引分别对应的三角函数值，例如，sin[0]～sin[256]。
在一种可能的实现方式中，当M≠N时，所述位置信息的俯仰角信息对应的三角函数值满足如下公式：
其中，表示所述位置信息的俯仰角信息，或者，表示所述与所述位置信息的俯仰角信息对应的参考点的俯仰角信息；sin_table_N()表示所述俯仰三角函数表；
当M≠N时，所述位置信息的水平角信息对应的三角函数值满足如下公式：
sin(θ)＝sin_table_M(θ’)
sin(θ)＝sin_table_M(M/2-θ’)
sin(θ)＝-sin_table_M(θ’-M/2)
sin(θ)＝-sin_table_M(M-θ’)
cos(θ)＝sin_table_M(M/4-θ’)
cos(θ)＝-sin_table_M(θ’-M/4)
cos(θ)＝-sin_table_M(3M/4-θ’)
cos(θ)＝sin_table_M(θ’-3M/4)
其中，θ表示所述位置信息的水平角信息，或者，θ表示所述与所述位置信息的水平角信息对应的参考点的水平角信息；sin_table_M()表示所述水平三角函数表；
或，
当M＝N时，所述位置信息的俯仰角信息对应的三角函数值满足如下公式：
其中，表示所述位置信息的俯仰角信息，或者，表示所述与所述位置信息的俯仰角信息对应的参考点的俯仰角信息；sin_table_N()表示所述俯仰三角函数表；
当M＝N时，所述位置信息的水平角信息对应的三角函数值满足如下公式：
sin(θ)＝sin_table(θ’)
sin(θ)＝sin_table(N/2-θ’)
sin(θ)＝-sin_table(θ’-N/2)
sin(θ)＝-sin_table(N-θ’)
cos(θ)＝sin_table(N/4-θ’)
cos(θ)＝-sin_table(θ’-N/4)
cos(θ)＝-sin_table(3N/4-θ’)
cos(θ)＝sin_table(θ’-3N/4)
其中，θ表示所述位置信息的水平角信息，或者，θ表示所述与所述位置信息的水平角信息对应的参考点的水平角信息；sin_table_N()表示所述俯仰三角函数表；
或，
当M≠N，M＝K1×N，K1≥2时，所述位置信息的水平角信息对应的三角函数值满足如下公式：
sin(θ)＝sin_table_M(θ’)
sin(θ)＝sin_table_M(M/2-θ’)
sin(θ)＝-sin_table_M(θ’-M/2)
sin(θ)＝-sin_table_M(M-θ’)
cos(θ)＝sin_table_M(M/4-θ’)
cos(θ)＝-sin_table_M(θ’-M/4)
cos(θ)＝-sin_table_M(3M/4-θ’)
cos(θ)＝sin_table_M(θ’-3M/4)
其中，θ表示所述位置信息的水平角信息，或者，θ表示所述与所述位置信息的水平角信息对应的参考点的水平角信息；sin_table_M()表示所述水平三角函数表；
当M≠N，M＝K1×N，K1≥2时，所述位置信息的俯仰角信息对应的三角函数值满足如下公式：
其中，表示所述位置信息的俯仰角信息，或者，表示所述与所述位置信息的俯仰角信息对应的参考点的俯仰角信息；sin_table_M()表示所述水平三角函数表；
或，
当M≠N，N＝K2×M，K2≥2时，所述位置信息的俯仰角信息对应的三角函数值满足如下公式：
其中，表示所述位置信息的俯仰角信息，或者，表示所述与所述位置信息的俯仰角信息对应的参考点的俯仰角信息；sin_table_N()表示所述俯仰三角函数表；
当M≠N，N＝K2×M，K2≥2时，所述位置信息的水平角信息对应的三角函数值满足如下公式：
sin(θ)＝sin_table_N(K2×θ’)
sin(θ)＝sin_table_N(N/2-K2×θ’)
sin(θ)＝-sin_table_N(K2×θ’-N/2)
sin(θ)＝-sin_table_N(N-K2×θ’)
cos(θ)＝sin_table_N(N/4-K2×θ’)
cos(θ)＝-sin_table_N(K2×θ’-N/4)
cos(θ)＝-sin_table_N(3N/4-K2×θ’)
cos(θ)＝sin_table_N(K2×θ’-3N/4)
其中，θ表示所述位置信息的水平角信息，或者，θ表示所述与所述位置信息的水平角信息对应的参考点的水平角信息；sin_table_N()表示所述俯仰三角函数表；
或，
当M≠N时，所述位置信息的俯仰角信息对应的三角函数值满足如下公式：
sin(θ)＝cos_table_N(N/4-θ’)
sin(θ)＝cos_table_N(θ’-N/4)
sin(θ)＝-cos_table_N(3N/4-θ’)
sin(θ)＝-cos_table_N(θ’-3N/4)
cos(θ)＝cos_table_N(θ’)
cos(θ)＝-cos_table_N(N/2-θ’)
cos(θ)＝-cos_table_N(θ’-N/2)
cos(θ)＝cos_table_N(N-θ’)
其中，表示所述位置信息的俯仰角信息，或者，表示所述与所述位置信息的俯仰角信息对应的参考点的俯仰角信息；cos_table_N()表示所述俯仰三角函数表；
当M≠N时，所述位置信息的水平角信息对应的三角函数值满足如下公式：
sin(θ)＝cos_table_M(M/4-θ’)
sin(θ)＝cos_table_M(θ’-M/4)
sin(θ)＝-cos_table_M(3M/4-θ’)
sin(θ)＝-cos_table_M(θ’-3M/4)
cos(θ)＝cos_table_M(θ’)
cos(θ)＝-cos_table_M(M/2-θ’)
cos(θ)＝-cos_table_M(θ’-M/2)
cos(θ)＝cos_table_M(M-θ’)
其中，θ表示所述位置信息的水平角信息，或者，θ表示所述与所述位置信息的水平角信息对应的参考点的水平角信息；cos_table_M()表示所述水平三角函数表。
第三方面，本申请提供一种音频处理设备，包括：一个或多个处理器；存储器，用于存储一个或多个程序；当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如上述第一方面中任一项所述的方法。
第四方面，本申请提供一种计算机可读存储介质，包括计算机程序，所述计算机程序在计算机上被执行时，使得所述计算机执行上述第一方面中任一项所述的方法。
附图说明
图1为本申请音频播放系统的一个示例性的结构图；
图2为本申请音频译码系统10的一个示例性的结构图；
图3为本申请HOA编码装置的一个示例性的结构图；
图4是本申请HOA系数的获取方法的一个示例性的流程图；
图5为本申请预设球面的一个示例性的示意图；
图6为本申请位置信息的一个示例性的示意图；
图7为本申请HOA系数的获取装置的一个示例性的结构图。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚，下面将结合本申请中的附图，对本申请中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获取的所有其他实施例，都属于本申请保护的范围。
本申请的说明书实施例和权利要求书及附图中的术语“第一”、“第二”等仅用于区分描述的目的，而不能理解为指示或暗示相对重要性，也不能理解为指示或暗示顺序。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元。方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
应当理解，在本申请中，“至少一个(项)”是指一个或者多个，“多个”是指两个或两个以上。“和/或”，用于描述关联对象的关联关系，表示可以存在三种关系，例如，“A和/或B”可以表示：只存在A，只存在B以及同时存在A和B三种情况，其中A，B可以是单数或者复数。字符“/”一般表示前后关联对象是一种“或”的关系。“以下至少一项(个)”或其类似表达，是指这些项中的任意组合，包括单项(个)或复数项(个)的任意组合。例如，a，b或c中的至少一项(个)，可以表示：a，b，c，“a和b”，“a和c”，“b和c”，或“a和b和c”，其中a，b，c可以是单个，也可以是多个。字符“～”连接的两个数值一般表示一个取值范围，该取值范围包含“～”连接的两个数值。
本申请涉及到的相关名词解释：
音频帧：音频数据是流式的，在实际应用中，为了便于音频处理和传输，通常取一时长内的音频数据量作为一帧音频，该时长被称为“采样时间”，可以根据编解码器和具体应用的需求确定其值，例如该时长为2.5ms～60ms，ms为毫秒。
音频信号：音频信号是带有语音、音乐和音效的有规律的声波的频率、幅度变化信息载体。音频是一种连续变化的模拟信号，可用一条连续的曲线来表示，称为声波。音频通过模数转换或计算机生成的数字信号即为音频信号。声波有三个重要参数：频率、幅度和相位，这也就决定了音频信号的特征。
以下是本申请所应用的系统架构。
图1为本申请音频播放系统的一个示例性的结构图，如图1所示，该音频播放系统包括：音频发送设备和音频接收设备，其中，音频发送设备包括例如手机、电脑(笔记本电脑、台式电脑等)、平板(手持平板、车载平板等)等可以进行音频编码并发送音频码流的设备；音频接收设备包括例如TWS耳机、普通无线耳机、音响、智能手表、智能眼镜等可以接收音频码流、解码音频码流并播放的设备。
音频发送设备和音频接收设备之间可以建立蓝牙连接，二者之间可以支持语音和音乐的传输。音频发送设备和音频接收设备的较为广泛的示例是手机与TWS耳机、无线头戴式耳机或者无线颈圈式耳机之间，或者手机与其他终端设备(例如智能音箱、智能手表、智能眼镜和车载音箱等)之间。可选的，音频发送设备和音频接收设备的示例也可以是平板、笔记本电脑或者台式电脑与TWS耳机、无线头戴式耳机、无线颈圈式耳机或其他终端设备(例如智能音箱、智能手表、智能眼镜和车载音箱)之间。
需要说明的是，音频发送设备和音频接收设备之间除蓝牙连接外，还可以通过其他通信方式连接，例如WiFi连接、有线连接或其他无线连接等，本申请对此不做具体限定。
图2为本申请音频译码系统10的一个示例性的结构图，如图2所示，音频译码系统10可包括源设备12和目的设备14，源设备12可以是图1的音频发送设备，目的设备14可以是图1的音频接收设备。源设备12产生经编码的码流信息，因此，源设备12也可以被称为音频编码设备。目的设备14可对由源设备12所产生的经编码的码流信息进行解码，因此，目的设备14也可以被称为音频解码设备。本申请中，源设备12、音频编码设备可以被统一称作音频发送设备，目的设备14、音频解码设备可以被统一称作音频接收设备。
源设备12包括编码器20，可选地，可包括音频源16、音频预处理器18、通信接口22。
音频源16，可以包括或可以为任何类别的音频捕获设备，例如，捕获现实世界声音，和/或任何类别的音频生成设备，例如，计算机音频处理器，或用于获取和/或提供现实世界音频、计算机动画音频(例如，屏幕内容、虚拟现实(virtual reality，VR)中的音频)的任何类别设备，和/或其任何组合(例如，增强现实(augmented reality，AR)中的音频、混合现实(mixed Reality，MR)中的音频和/或扩展现实(extended Reality，XR)中的音频)。音频源16可以为用于捕获音频的麦克风或者用于存储音频的存储器，音频源16还可以包括存储先前捕获或产生的音频和/或获取或接收音频的任何类别的(内部或外部)接口。当音频源16为麦克风时，音频源16可例如为本地的或集成在源设备中的音频采集装置；当音频源16为存储器时，音频源16可为本地的或例如集成在源设备中的集成存储器。当所述音频源16包括接口时，接口可例如为从外部音频源接收音频的外部接口，外部音频源例如为外部音频捕获设备，比如话筒、麦克风、外部存储器或外部音频生成设备，外部音频生成设备例如为外部计算机音频处理器、计算机或服务器。接口可以为根据任何专有或标准化接口协议的任何类别的接口，例如有线或无线接口、光接口。
本申请中，音频源16获取当前场景音频信号，该当前场景音频信号是指对空间中麦克风所在位置的声场进行采集得到的音频信号，当前场景音频信号也可以称为原始场景音频信号。例如，当前场景音频信号可以是通过高阶立体混响(higher order ambisonics，HOA)技术得到的音频信号。音频源16获取待编码的HOA信号，例如，可以采用实际采集设备获取HOA信号或采用人工音频对象合成HOA信号。可选的，待编码的HOA信号可以是时域HOA信号或者频域HOA信号。
音频预处理器18，用于接收原始音频信号并对原始音频信号执行预处理，以获取经预处理的音频信号。例如，音频预处理器18执行的预处理可以包括整修或去噪。
编码器20，用于接收经预处理的音频信号，对经预处理的音频信号进行处理，从而提供经编码的码流信息。
源设备12中的通信接口22可用于接收码流信息并通过通信信道13向目的设备14发送该码流。通信信道13例如为直接有线或无线连接，任何类别的网络例如为有线或无线网络或其任何组合，或任何类别的私网和公网，或其任何组合。
目的设备14包括解码器30，可选地，可包括通信接口28、音频后处理器32和播放设备34。
目的设备14中的通信接口28用于直接从源设备12接收码流信息，并将码流信息提供给解码器30。通信接口22和通信接口28可用于通过源设备12与目的设备14之间的通信信道13发送或接收码流信息。
通信接口22和通信接口28均可配置为如图2中从源设备12指向目的设备14的对应通信信道13的箭头所指示的单向通信接口，或双向通信接口，并且可用于发送和接收消息等，以建立连接，确认并交换与通信链路和/或编码音频数据等数据传输相关的任何其它信息，等等。
解码器30，用于接收码流信息，并解码码流信息得到经解码的音频数据。
音频后处理器32，用于对解码的音频数据进行后处理，得到后处理后的音频数据。音频后处理器32执行的后处理可以包括例如修剪或重采样等。
播放设备34，用于接收后处理后的音频数据，以向用户或收听者播放音频。播放设备34可以为或包括任意类型的用于播放重建后音频的播放器，例如，集成或外部扬声器。例如，扬声器可包括喇叭、音响等。
图3为本申请HOA编码装置的一个示例性的结构图，如图3所示，HOA编码装置可以应用于上述音频译码系统10的编码器20中。HOA编码装置包括：虚拟扬声器配置单元、编码分析单元、虚拟扬声器集合生成单元、虚拟扬声器选择单元、虚拟扬声器信号生成单元和核心编码器处理单元。其中，
虚拟扬声器配置单元，用于根据编码器配置信息对虚拟扬声器进行配置，以得到虚拟扬声器配置参数。编码器配置信息包括且不限于：HOA阶数，编码比特率，用户自定义信息等，虚拟扬声器配置参数包括且不限于：虚拟扬声器的个数，虚拟扬声器的HOA阶数等。
虚拟扬声器配置单元输出的虚拟扬声器配置参数作为虚拟扬声器集合生成单元的输入。
编码分析单元，用于对待编码HOA信号进行编码分析，例如分析待编码HOA信号的声场分布，包括待编码HOA信号的声源个数、方向性、弥散度等特征，作为决定如何选择目标虚拟扬声器的判断条件之一。
不限定的是，本申请中，HOA编码装置中也可以不包括编码分析单元，即HOA编码装置可以不对输入信号进行分析，则采用一种默认配置决定如何选择目标虚拟扬声器。
其中，HOA编码装置获取待编码HOA信号，例如可以将从实际采集设备记录的HOA信号或采用人工音频对象合成的HOA信号作为编码器的输入，同时编码器输入的待编码HOA信号可以是时域HOA信号也可以是频域HOA信号。
虚拟扬声器集合生成单元，用于生成虚拟扬声器集合，该虚拟扬声器集合中可以包括：多个虚拟扬声器，虚拟扬声器集合中的虚拟扬声器也可以称为“候选虚拟扬声器”。
虚拟扬声器集合生成单元生成指定的候选虚拟扬声器HOA系数。由虚拟扬声器配置单元提供的候选虚拟扬声器的坐标(即位置信息)和候选虚拟扬声器的HOA阶数用于生成候选虚拟扬声器HOA系数。候选虚拟扬声器的坐标确定方法包括且不限于按等距规则产生K个虚拟扬声器、根据听觉感知原理生成非均匀分布的K个候选虚拟扬声器。根据候选虚拟扬声器的个数生成分布均匀的候选虚拟扬声器的坐标。
接下来生成候选虚拟扬声器HOA系数。
这里采用本申请提供的HOA系数的获取方法获得每一个虚拟扬声器的HOA系数并进行表示，得到虚拟扬声器集合。
虚拟扬声器集合生成单元输出的候选虚拟扬声器的HOA系数作为虚拟扬声器选择单元的输入。
虚拟扬声器选择单元，用于根据待编码HOA信号从虚拟扬声器集合中的多个候选虚拟扬声器中选择出目标虚拟扬声器，该目标虚拟扬声器可以称为“与待编码HOA信号匹配的虚拟扬声器”，或者简称为匹配虚拟扬声器。
虚拟扬声器选择单元根据待编码HOA信号与虚拟扬声器集合生成单元输出的候选虚拟扬声器HOA系数，选择出指定的匹配虚拟扬声器。
接下来对匹配虚拟扬声器的选择方法进行举例说明：在一种可能的实现方式中，使用候选虚拟扬声器HOA系数匹配与待编码HOA信号做内积，选取内积绝对值最大的候选虚拟扬声器为目标虚拟扬声器，即匹配虚拟扬声器，并将待编码HOA信号在该候选虚拟扬声器的投影叠加到该候选虚拟扬声器HOA系数的线性组合上，然后将投影向量从待编码HOA信号中减去得到差值，对差值重复上述过程实现迭代计算，每迭代一次产生一个匹配虚拟扬声器，输出匹配虚拟扬声器坐标和匹配虚拟扬声器HOA系数。可以理解的是，匹配虚拟扬声器会选取多个，每迭代一次产生一个匹配虚拟扬声器。(除此之外，不限定其他实现方法)
虚拟扬声器选择单元输出的目标虚拟扬声器的坐标和目标虚拟扬声器的HOA系数作为虚拟扬声器信号生成单元的输入。
虚拟扬声器信号生成单元，用于根据待编码HOA信号和目标虚拟扬声器的属性信息生成虚拟扬声器信号，其中当属性信息为位置信息时，根据所述目标虚拟扬声器的位置信息确定所述目标虚拟扬声器的HOA系数，当属性信息包括HOA系数时，从所述属性信息中获取所述目标虚拟扬声器的HOA系数。
虚拟扬声器信号生成单元通过待编码HOA信号和目标虚拟扬声器的HOA系数计算虚拟扬声器信号。
虚拟扬声器的HOA系数用矩阵A表示，用矩阵A可以线性组合出待编码HOA信号，进一步的可以采用最小二乘方法求得理论的最优解w，即为虚拟扬声器信号，例如可以采用如下计算公式：
w＝A-1X，
其中，A-1代表矩阵A的逆矩阵，矩阵A的大小为(M×C)，C为目标虚拟扬声器个数，M为N阶的HOA系数的通道个数，M＝(N+1)2，a表示目标虚拟扬声器的HOA系数，例如，
X代表待编码HOA信号，矩阵X的大小为(M×L)，M为N阶的HOA系数的通道个数，L为时域或频域样点个数，x表示待编码HOA信号的系数，例如，
虚拟扬声器信号生成单元输出的虚拟扬声器信号作为核心编码器处理单元的输入。
核心编码器处理单元，用于对虚拟扬声器信号进行核心编码器处理，得到传输码流。
核心编码器处理包括且不限于变换、量化、心理声学模型、码流产生等，可以对频域传输通道进行处理也可以对时域传输通道进行处理，此处不做限定。
基于上述实施例的描述，本申请提供了一种HOA系数的获取方法。
图4是本申请HOA系数的获取方法的一个示例性的流程图。该过程400可由上述实施例中的编码器20或解码器30执行，即由音频发送设备中的编码器20实现音频编码，然后将码流信息发送给音频接收设备，由音频接收设备中的解码器30对码流信息进行解码以获得目标音频帧，进而基于该目标音频帧渲染得到对应于一个或多个虚拟扬声器的声场音频信号。过程400描述为一系列的步骤或操作，应当理解的是，过程400可以以各种顺序执行和/或同时发生，不限于图4所示的执行顺序。如图4所示，该方法包括：
步骤401、获取虚拟扬声器在预设球面上的位置信息，位置信息包括俯仰角信息和/或水平角信息。
预设球面上包含M个经线圈和N个纬线圈，M和N可以相同也可以不同，M和N均为正整数，例如，N为512，768或1024等等，M为512，768或1024等等。经线圈和纬线圈的交叉点称作参考点，位于同一个经线圈上的N个参考点等间距排列，位于同一个纬线圈上的M个参考点等间距排列。为了保证参考点可以尽可能覆盖到所有虚拟扬声器，因此可以满足M×N≥K，其中K为虚拟扬声器总数，例如K可以取值为2048，1669，1024，530等。例如，图5为本申请预设球面的一个示例性的示意图，如图5所示，将空间声场假设为一个以渲染中心为球心的球体空间，虚拟扬声器可以部署于球体空间的球面上，该球面被M个经线圈和N个纬线圈划分，其中，任意一个经线圈所在的平面均包含球心，因此M个经线圈的半径与球体空间的半径相等；纬线圈中只有赤道圈的半径与球体空间的半径相等，其他纬线圈越远离赤道圈半径越小。本申请中，两两相邻的经线圈之间的水平角度差是相等的，即M个经线圈在球面上是均匀分布的，因此位于同一个纬线圈上的M个参考点等间距排列；两两相邻的纬线圈之间的俯仰角度差也是相等的，即N个纬线圈在球面上是均匀分布的，因此位于同一个经线圈上的N个参考点等间距排列。相应的，由于位于同一个经线圈上的N个参考点等间距排列，因此位于同一个经线圈上的所有相邻的参考点的俯仰角度差相等；由于位于同一个纬线圈上的M个参考点等间距排列，因此位于同一个纬线圈上的所有相邻的参考点的水平角度差相等。
上述虚拟扬声器可以是针对声音渲染设置的任意一个虚拟扬声器。虚拟扬声器在上述球面上的位置信息包括俯仰角信息和/或水平角信息。例如，图6为本申请位置信息的一个示例性的示意图，如图6所示，虚拟扬声器设置的位置和球心之间的连线与预设水平面(例如赤道圈所在平面，或者南极点所在的平面，或者北极点所在的平面，其中，南极点所在的平面垂直于南极点和北极点之间的连线，北极点所在的平面垂直于南极点和北极点之间的连线)之间的夹角为虚拟扬声器的位置信息的俯仰角信息；虚拟扬声器设置的位置和球心之间的连线在水平面上的投影与设定初始方向的夹角为虚拟扬声器的位置信息的水平角信息。
可选的，虚拟扬声器的位置信息中的俯仰角信息可以用俯仰角度表示，也可以用俯仰角索引表示，同理，虚拟扬声器的位置信息中的水平角信息可以用水平角度表示，也可以用水平角索引表示。
例如，经线圈上有N个参考点，该N个参考点是等间距排列的，相邻两个参考点的俯仰角信息差为2π/N，因此该N个参考点各自的俯仰角信息可以由参考点的俯仰角索引表示，即该N个参考点的俯仰角信息可以表示为0～N-1。如果虚拟扬声器的位置正好选在参考点上，就可以用前述俯仰角索引表示该虚拟扬声器的俯仰角信息。同理，纬线圈上有M个参考点，该M个参考点是等间距排列的，相邻两个参考点的水平角信息差为2π/M，因此该M个参考点各自的水平角信息可以由参考点的水平角索引表示，即该M个参考点的水平角信息可以表示为0～M-1。如果虚拟扬声器的位置正好选在参考点上，就可以用前述水平角索引表示该虚拟扬声器的水平角信息。由此可见，参考点的俯仰角索引和俯仰角度之间具有一一对应关系，参考点的水平角索引和水平角度之间也具有一一对应关系。而如果虚拟扬声器的位置不在参考点的位置上，那么就无法确定是用哪一个参考点的俯仰角索引和/或水平角索引来表示其位置信息，因此可以用俯仰角度和/或水平角度的方式来表示其位置信息。应当理解的，即使虚拟扬声器的位置正好选在参考点上，也可以采用俯仰角度和/或水平角度的方式来表示其位置信息，对此不做具体限定。
步骤402、根据位置信息以及预设的参考三角函数表获取位置信息对应的三角函数值。
参考三角函数表包括俯仰三角函数表和/或水平三角函数表。其中，
俯仰三角函数表包括多个第一参考点的俯仰角索引对应的多个三角函数值，多个第一参考点是第一经线圈上的参考点，多个第一参考点的数量不小于表示向下取整，第一经线圈是M个经线圈的其中之一。
本申请从M个经线圈中任意选取一条视作第一经线圈，用于计算俯仰三角函数表的参考点(即上述多个第一参考点)是从该第一经线圈上的N个参考点中选取的。基于三角函数的原理，可以将0°(0)～360°(2π)的角度范围划分为四个子范围，即0°(0)～90°(π/2)，90°(π/2)～180°(π)，180°(π)～270°(3π/2)以及270°(3π/2)～360°(2π)，通过该四个子范围的任意一个子范围计算得到的三角函数值可以扩展到0°(0)～360°(2π)的角度范围的三角函数值。因此本申请在选取上述多个第一参考点时，其数量不小于即从第一经线圈上的N个参考点中选取至少四分之一个参考点作为上述多个第一参考点。进一步的，选取的多个第一参考点中，任意两个参考点的俯仰角度在计算三角函数值时不会重叠。
可选的，上述多个第一参考点的俯仰角度可以属于上述四个子范围的其中之一，例如，多个第一参考点的俯仰角度在0°(0)～90°(π/2)内；或者，多个第一参考点的俯仰角度在90°(π/2)～180°(π)内；或者，多个第一参考点的俯仰角度在180°(π)～270°(3π/2)内；或者，多个第一参考点的俯仰角度在270°(3π/2)～360°(2π)内。
可选的，上述多个第一参考点的俯仰角度可以属于上述四个子范围的其中至少两个，例如，多个第一参考点中的第一部分参考点的俯仰角度在0°(0)～90°(π/2)内，第二部分参考点的俯仰角度在90°(π/2)～180°(π)内，这两部分参考点的俯仰角度在计算三角函数值时不会重叠；或者，多个第一参考点中的部分参考点的俯仰角度在90°(π/2)～180°(π)内，另一部分参考点的俯仰角度在180°(π)～270°(3π/2)内，这两部分参考点的俯仰角度在计算三角函数值时不会重叠；或者，多个第一参考点中的部分参考点的俯仰角度在180°(π)～270°(3π/2)内，另一部分参考点的俯仰角度在270°(3π/2)～360°(2π)内，这两部分参考点的俯仰角度在计算三角函数值时不会重叠；或者，多个第一参考点中的第一部分参考点的俯仰角度在0°(0)～90°(π/2)内，第二部分参考点的俯仰角度在90°(π/2)～180°(π)内，第三部分参考点的俯仰角度在180°(π)～270°(3π/2)内，这三部分参考点的俯仰角度在计算三角函数值时不会重叠；或者，多个第一参考点中的第一部分参考点的俯仰角度在90°(π/2)～180°(π)内，第二部分参考点的俯仰角度在180°(π)～270°(3π/2)内，第三部分参考点的俯仰角度在270°(3π/2)～360°(2π)内，这三部分参考点的俯仰角度在计算三角函数值时不会重叠；或者，多个第一参考点中的第一部分参考点的俯仰角度在0°(0)～90°(π/2)内，第二部分参考点的俯仰角度在90°(π/2)～180°(π)内，第三部分参考点的俯仰角度在180°(π)～270°(3π/2)内，第四部分参考点的俯仰角度在270°(3π/2)～360°(2π)内，这四部分参考点的俯仰角度在计算三角函数值时不会重叠。
上述不会重叠可以是指属于不同子范围的俯仰角度计算得到的三角函数值不相等，例如，多个第一参考点共有四个(n1、n2、n3和n4)，其中，n1的俯仰角度为0，n2的俯仰角度为π/6，n3的俯仰角度为π/3，n4的俯仰角度为π/2；或者，n1的俯仰角度为0，n2的俯仰角度为2/3π，n3的俯仰角度为4/3π，n4的俯仰角度为2π。需要说明的是，本申请中对多个第一参考点的俯仰角度所述的子范围不做具体限定，只要该多个第一参考点的俯仰角度的三角函数值可以覆盖其中一个子范围的全部角度的三角函数值即可。
相应的，上述位置信息的俯仰角度取值范围为0～2π。
如上所述，上述多个第一参考点也可以用各自的俯仰角索引进行标识，基于俯仰角索引和相邻纬线圈的俯仰角度差可以在第一参考点的俯仰角度和俯仰角度差之间进行转换。
本申请中多个第一参考点的俯仰角索引对应的多个三角函数可以是正弦函数值，也可以是余弦函数值。
可选的，第i个第一参考点的俯仰角索引对应的正弦函数值满足如下公式(1)：
其中，i＝0,1,…,N’，ri表示上述第一经线圈的半径。可选的，i的取值也可以是在范围内不连续的，例如，i＝0,2,4,5,9,…,N’。
可选的，第i个第一参考点的俯仰角索引对应的余弦函数值满足如下公式(2)：
其中，i＝0,1,…,N’，ri表示上述第一经线圈的半径。
水平三角函数表包括多个第二参考点的水平角索引对应的多个三角函数值，多个第二参考点是第一纬线圈上的参考点，多个第二参考点的数量不小于表示向下取整，第一纬线圈是N个纬线圈的其中之一。
本申请从N个纬线圈中任意选取一条视作第一纬线圈，用于计算水平三角函数表的参考点(即上述多个第二参考点)是从该第一纬线圈上的M个参考点中选取的。同样基于三角函数的原理，可以将0°(0)～360°(2π)的角度范围划分为四个子范围，即0°(0)～90°(π/2)，90°(π/2)～180°(π)，180°(π)～270°(3π/2)以及270°(3π/2)～360°(2π)，通过该四个子范围的任意一个子范围计算得到的三角函数值可以扩展到0°(0)～360°(2π)的角度范围的三角函数值。因此本申请在选取上述多个第二参考点时，其数量不小于即从第一纬线圈上的M个参考点中选取至少四分之一个参考点作为上述多个第二参考点。进一步的，选取的多个第二参考点中，任意两个参考点的水平角度在计算三角函数值时不会重叠。
可选的，上述多个第二参考点的水平角度可以属于上述四个子范围的其中之一，例如，多个第二参考点的水平角度在0°(0)～90°(π/2)内；或者，多个第二参考点的水平角度在90°(π/2)～180°(π)内；或者，多个第二参考点的水平角度在180°(π)～270°(3π/2)内；或者，多个第二参考点的水平角度在270°(3π/2)～360°(2π)内。
可选的，上述多个第二参考点的水平角度可以属于上述四个子范围的其中至少两个，例如，多个第二参考点中的第一部分参考点的水平角度在0°(0)～90°(π/2)内，第二部分参考点的水平角度在90°(π/2)～180°(π)内，这两部分参考点的水平角度在计算三角函数值时不会重叠；或者，多个第二参考点中的第一部分参考点的水平角度在90°(π/2)～180°(π)内，第二部分参考点的水平角度在180°(π)～270°(3π/2)内，这两部分参考点的水平角度在计算三角函数值时不会重叠；或者，多个第二参考点中的第一部分参考点的水平角度在180°(π)～270°(3π/2)内，第二部分参考点的水平角度在270°(3π/2)～360°(2π)内，这两部分参考点的水平角度在计算三角函数值时不会重叠；或者，多个第二参考点中的第一部分参考点的水平角度在0°(0)～90°(π/2)内，第二部分参考点的水平角度在90°(π/2)～180°(π)内，第三部分参考点的水平角度在180°(π)～270°(3π/2)内，这三部分参考点的水平角度在计算三角函数值时不会重叠；或者，多个第二参考点中的第一部分参考点的水平角度在90°(π/2)～180°(π)内，第二部分参考点的水平角度在180°(π)～270°(3π/2)内，第三部分参考点的水平角度在270°(3π/2)～360°(2π)内，这三部分参考点的水平角度在计算三角函数值时不会重叠；或者，多个第二参考点中的第一部分参考点的水平角度在0°(0)～90°(π/2)内，第二部分参考点的水平角度在90°(π/2)～180°(π)内，第三部分参考点的水平角度在180°(π)～270°(3π/2)内，第四部分参考点的水平角度在270°(3π/2)～360°(2π)内，这四部分参考点的水平角度在计算三角函数值时不会重叠。
上述不会重叠可以是指属于不同子范围的水平角度计算得到的三角函数值不相等，例如，多个第二参考点共有四个(n1、n2、n3和n4)，其中，n1的水平角度为0，n2的水平角度为π/6，n3的水平角度为π/3，n4的水平角度为π/2；或者，n1的水平角度为0，n2的水平角度为2/3π，n3的水平角度为4/3π，n4的水平角度为2π。需要说明的是，本申请中对多个第二参考点的水平角度所述的子范围不做具体限定，只要该多个第二参考点的水平角度的三角函数值可以覆盖其中一个子范围的全部角度的三角函数值即可。
相应的，上述位置信息的水平角度取值范围为0～2π。
如上所述，上述多个第二参考点也可以用各自的水平角索引进行标识，基于水平角索引和相邻经线圈的水平角度差可以在第二参考点的水平角度和水平角度差之间进行转换。
本申请中多个第二参考点的水平角索引对应的多个三角函数可以是正弦函数值，也可以是余弦函数值。
可选的，第i个第二参考点的水平角索引对应的正弦函数值满足如下公式(3)：
其中，j＝0,1,…,M’，rj表示上述第一纬线圈的半径。可选的，j的取值也可以是在范围内不连续的，例如，i＝1,3,5,8,9,10,…,M’。
可选的，第i个第二参考点的水平角索引对应的余弦函数值满足如下公式(4)：
其中，j＝0,1,…,M’，rj表示上述第一纬线圈的半径。
需要说明的是，上述俯仰三角函数表和水平三角函数表可以是预先计算并存储的，也可以根据M和N的值实时计算并存储。可选的，由于编码器和解码器均采用分帧处理，因此在HOA系数计算过程中，当前音频帧的M与前一音频帧的M取值相等时，处理当前帧时可以不用重新计算水平三角函数表，直接使用前一音频帧的水平三角函数表即可；当前音频帧的N与前一音频帧的N取值相等时，处理当前帧时可以不用重新计算俯仰三角函数表，直接使用前一音频帧的俯仰三角函数表即可。
本申请中可以根据虚拟扬声器在预设球面上的位置信息的俯仰角信息以及上述俯仰三角函数表获取位置信息的俯仰角信息对应的三角函数值；和/或，根据虚拟扬声器在预设球面上的位置信息的水平角信息以及上述水平三角函数表获取位置信息的水平角信息对应的三角函数值。
可见针对虚拟扬声器在预设球面上的位置信息，可以根据其位置信息包含的具体内容，以及预先设置的参考三角函数表的具体内容来获得该位置信息对应的三角函数值，即位置信息包含俯仰角信息，可以根据该俯仰角信息和参考三角函数表(如果有俯仰三角函数表，则用俯仰三角函数表，如果没有俯仰三角函数表，则用水平三角函数表)获得位置信息的俯仰角信息对应的三角函数值；位置信息包含水平角信息，可以根据该水平角信息和参考三角函数表(如果有水平三角函数表，则用水平三角函数表，如果没有水平三角函数表，则用俯仰三角函数表)获得位置信息的水平角信息对应的三角函数值。
本申请可以采用以下几种方法获取位置信息的俯仰角信息对应的三角函数值：
一、当位置信息的俯仰角信息与多个第一参考点中的一个参考点的俯仰角索引对应上时，将俯仰三角函数表中一个参考点的俯仰角索引对应的三角函数值作为位置信息的俯仰角信息对应的三角函数值。
位置信息的俯仰角信息与多个第一参考点中的一个参考点的俯仰角索引对应上表示虚拟扬声器在预设球面上的位置对应的俯仰角度正好与俯仰三角函数表中的某一个俯仰角索引对应的俯仰角度相等，或者，虚拟扬声器在预设球面上的位置对应的俯仰角索引正好与俯仰三角函数表中的某一个俯仰角索引相等，即虚拟扬声器在预设球面上的位置与选取出来用于获取俯仰三角函数表的参考点的位置重合。因此可以直接从俯仰三角函数表中查询该俯仰角信息的三角函数值即可。
二、当位置信息的俯仰角信息与N个参考点中除多个第一参考点外的一个参考点的俯仰角索引对应上时，将N个参考点中除多个第一参考点外的一个参考点的俯仰角索引转换为多个第一参考点中的一个参考点的俯仰角索引，将俯仰三角函数表中一个参考点的俯仰角索引对应的三角函数值作为位置信息的俯仰角信息对应的三角函数值。
位置信息的俯仰角信息与N个参考点中除多个第一参考点外的一个参考点的俯仰角索引对应上表示虚拟扬声器在预设球面上的位置对应的俯仰角信息与俯仰三角函数表中的各个俯仰角索引对应的俯仰角度均不相等，但与第一经线圈上的某一其他参考点的俯仰角索引对应的角度相等，或者，虚拟扬声器在预设球面上的位置对应的俯仰角索引与第一经线圈上的某一其他参考点的俯仰角索引相等，即虚拟扬声器在预设球面上的位置与某一参考点的位置重合，该参考点没有用于获取俯仰三角函数表，但其俯仰角索引与三角函数表中的俯仰角索引存在三角函数对应关系。因此可以基于三角函数的原理将将N个参考点中除多个第一参考点外的一个参考点的俯仰角索引转换为多个第一参考点中的一个参考点的俯仰角索引，然后查表得到多个第一参考点中的一个参考点的俯仰角索引对应的三角函数值，将多个第一参考点中的一个参考点的俯仰角索引对应的三角函数值作为位置信息的俯仰角信息对应的三角函数值。
三、当位置信息的俯仰角信息与参考点的俯仰角索引均不对应时，根据位置信息的俯仰角信息获得与位置信息的俯仰角信息对应的参考点的俯仰角信息；当与位置信息的俯仰角信息对应的参考点的俯仰角信息与多个第一参考点中的一个参考点的俯仰角索引对应上时，将俯仰三角函数表中一个参考点的俯仰角索引对应的三角函数值作为位置信息的俯仰角信息对应的三角函数值；当与位置信息的俯仰角信息对应的参考点的俯仰角信息与N个参考点中除多个第一参考点外的一个参考点的俯仰角索引对应上时，将N个参考点中除多个第一参考点外的一个参考点的俯仰角索引转换为多个第一参考点中的一个参考点的俯仰角索引，将俯仰三角函数表中一个参考点的俯仰角索引对应的三角函数值作为位置信息的俯仰角信息对应的三角函数值。
位置信息的俯仰角信息与参考点的俯仰角索引均不对应表示虚拟扬声器在预设球面上的位置既不在经线圈上，也不在纬线圈上。因此先获得与位置信息的俯仰角信息对应的参考点的俯仰角信息，再参照上述方法一或方法二的方式获取参考点的俯仰角信息对应的三角函数值，将该参考点的俯仰角信息对应的三角函数值作为位置信息的俯仰角信息对应的三角函数值。
可选的，俯仰角度和俯仰角索引满足如下公式(5)：
其中，ri表示第一经线圈的半径，round()表示取整。
本申请可以采用以下几种方法获取位置信息的水平角信息对应的三角函数值：
一、当位置信息的水平角信息与多个第二参考点中的一个参考点的水平角索引对应上时，将水平三角函数表中一个参考点的水平角索引对应的三角函数值作为位置信息的水平角信息对应的三角函数值。
位置信息的水平角信息与多个第二参考点中的一个参考点的水平角索引对应上表示虚拟扬声器在预设球面上的位置对应的水平角度正好与水平三角函数表中的某一个水平角索引对应的水平角度相等，或者，虚拟扬声器在预设球面上的位置对应的水平角索引正好与水平三角函数表中的某一个水平角索引相等，即虚拟扬声器在预设球面上的位置与选取出来用于获取水平三角函数表的参考点的位置重合。因此可以直接从水平三角函数表中查询该水平角信息的三角函数值即可。
二、当位置信息的水平角信息与M个参考点中除多个第二参考点外的一个参考点的水平角索引对应上时，将M个参考点中除多个第二参考点外的一个参考点的水平角索引转换为多个第二参考点中的一个参考点的水平角索引，将水平三角函数表中一个参考点的水平角索引对应的三角函数值作为位置信息的水平角信息对应的三角函数值。
位置信息的水平角信息与M个参考点中除多个第二参考点外的一个参考点的水平角索引对应上表示虚拟扬声器在预设球面上的位置对应的水平角信息与水平三角函数表中的各个俯仰角索引对应的俯仰角度均不相等，但与第一纬线圈上的某一其他参考点的水平角索引对应的角度相等，或者，虚拟扬声器在预设球面上的位置对应的水平角索引与第一纬线圈上的某一其他参考点的水平角索引相等，即虚拟扬声器在预设球面上的位置与某一参考点的位置重合，该参考点没有用于获取水平三角函数表，但其水平角索引与三角函数表中的水平角索引存在三角函数对应关系。因此可以基于三角函数的原理将M个参考点中除多个第二参考点外的一个参考点的水平角索引转换为多个第二参考点中的一个参考点的水平角索引，然后查表得到多个第二参考点中的一个参考点的水平角索引对应的三角函数值，将多个第二参考点中的一个参考点的水平角索引对应的三角函数值作为位置信息的水平角信息对应的三角函数值。
三、当位置信息的水平角信息与参考点的水平角索引均不对应时，根据位置信息的水平角信息获得与位置信息的水平角信息对应的参考点的水平角信息；当与位置信息的水平角信息对应的参考点的水平角信息与多个第二参考点中的一个参考点的水平角索引对应上时，将水平三角函数表中一个参考点的水平角索引对应的三角函数值作为位置信息的水平角信息对应的三角函数值；当与位置信息的水平角信息对应的参考点的水平角信息与M个参考点中除多个第二参考点外的一个参考点的水平角索引对应上时，将M个参考点中除多个第二参考点外的一个参考点的水平角索引转换为多个第二参考点中的一个参考点的水平角索引，将水平三角函数表中一个参考点的水平角索引对应的三角函数值作为位置信息的水平角信息对应的三角函数值。
位置信息的水平角信息与参考点的水平角索引均不对应表示虚拟扬声器在预设球面上的位置既不在经线圈上，也不在纬线圈上。因此先获得与位置信息的水平角信息对应的参考点的水平角信息，再参照上述方法一或方法二的方式获取参考点的水平角信息对应的三角函数值，将该参考点的水平角信息对应的三角函数值作为位置信息的水平角信息对应的三角函数值。
可选的，水平角度θ和水平角索引θ’满足如下公式(6)：
其中，rj表示第一纬线圈的半径，round()表示取整。
步骤403、根据位置信息对应的三角函数值获取虚拟扬声器的HOA系数。
以下对虚拟扬声器的HOA系数的生成过程进行说明：
声波在理想介质中传播，其波速为k＝w/c，角频率w＝2πf，f表示声波频率，c表示声速。因此声压p满足如下公式(7)：
▽2p+k2p＝0 (7)
其中，▽2为拉普拉斯算子。
在球坐标下求解公式(7)，声压p可以得到如下公式(8)：
其中，r表示球半径，θ表示水平角度，表示俯仰角度，k表示波速，s表示理想平面波的幅度，m表示HOA阶数序号，表示球贝塞尔函数，亦称作径向基函数，第一个j是虚数单位， 不随角度变化，是θ和对应的球谐函数，是声源方向的球谐函数。
立体混响(Ambisonics)系数为：
因此可以得到声压p的一般展开形式(10)：
上述公式(9)可以表明声场可以在球面上按球谐函数展开，其通过Ambisonics系数进行表示。
相应的，已知Ambisonics系数则可以重建声场，将公式(9)截断到第N项，以Ambisonics系数作为对声场的近似描述，则称为N阶的HOA系数，该HOA系数亦称作Ambisonics系数。N阶Ambisonics系数共有(N+1)2个通道。可选的，HOA阶数可以为2阶～10阶，将球谐函数按照HOA信号的一个采样点对应的系数进行叠加，就能实现该采样点对应的时刻空间声场的重构。根据该原理可以生成虚拟扬声器的HOA系数。将公式(9)中的θs和分别设置为虚拟扬声器的位置信息，即水平角度和俯仰角度，根据式(9)可以获得该虚拟扬声器的HOA系数，也称作Ambisonics系数。例如，针对3阶HOA信号，假设s＝1，其对应的16通道的HOA系数可通过球谐函数得到，3阶HOA信号对应的16通道的HOA系数计算公式具体如表1所示：
表1
表1中θ表示虚拟扬声器在预设球面上的位置信息的水平角信息，表示虚拟扬声器在预设球面上的位置信息的俯仰角信息，l表示HOA阶数，l＝0,1,…,N，m表示每一阶中的方向参数，m＝-l,…,l。按照表1中的极坐标的表达式，可以根据虚拟扬声器的位置信息，获得该虚拟扬声器的3阶HOA信号对应的16个通道的HOA系数。除此之外，本发明方法也可以用于获取任意阶数的球谐函数展开HOA系数。
本申请中，预先只需要获取预设球面上的某一经线圈上的所有参考点中的四分之一个参考点的三角函数值(俯仰三角函数表)，和/或，预设球面上的某一纬线圈上的所有参考点中的四分之一个参考点的三角函数值(水平三角函数表)，可以基于前述三角函数表获取到0°～360°的任意角度的三角函数值。表1中的其他三角函数例如幂次三角函数和倍角三角函数均可以通过现有技术进行展开，得到由sin三角函数和cos三角函数组成的一般展开形式，进而可以通过本发明方法得到相应的其他三角函数值。这样一方面可以减少预先阶段的计算量，另一方面由于三角函数表中的角度数量减少，从而需要存储的三角函数表的数据量也相应的减少，可以节省存储空间。
以下通过几个具体的实施例，示例性的说明获取位置信息的俯仰角信息和水平角信息分别对应的三角函数值的方法。表示位置信息的俯仰角信息，或者表示与位置信息对应的参考点的俯仰角信息，表示俯仰角索引；θ表示位置信息的水平角信息，或者表示与位置信息对应的参考点的水平角信息，θ’表示水平角索引。
需要说明的是，本申请中俯仰角信息和水平角信息可以用索引号的方式表示，例如，经线圈上有N个参考点，该N个参考点是等间距排列的，相邻两个参考点的俯仰角信息差为2π/N，因此该N个参考点各自的俯仰角信息可以由参考点的索引号表示，即该N个参考点的俯仰角信息可以表示为0～N-1，上述多个第一参考点的数量为那么该多个第一参考点的俯仰角信息可以表示为同理，纬线圈上有M个参考点，该M个参考点是等间距排列的，相邻两个参考点的水平角信息差为2π/M，因此该M个参考点各自的水平角信息可以由参考点的索引号表示，即该M个参考点的水平角信息可以表示为0～M-1，上述多个第二参考点的数量为那么该多个第二参考点的水平角信息可以表示为本申请可以采用上述公式(5)获得俯仰角度和俯仰角索引之间的关系，采用上述公式(6)获得水平角度θ和水平角索引θ’之间的关系。
实施例一
M≠N，采用上述步骤402的方法预先设定了俯仰正弦函数表sin_table_N和水平正弦函数表sin_table_M。
可以按照以下规则获取位置信息的俯仰角信息对应的正弦函数值和余弦函数值：
可以按照以下规则获取位置信息的水平角信息对应的正弦函数值和余弦函数值：
sin(θ)＝sin_table_M(θ’)
sin(θ)＝sin_table_M(M/2-θ’)
sin(θ)＝-sin_table_M(θ’-M/2)
sin(θ)＝-sin_table_M(M-θ’)
cos(θ)＝sin_table_M(M/4-θ’)
cos(θ)＝-sin_table_M(θ’-M/4)
cos(θ)＝-sin_table_M(3M/4-θ’)
cos(θ)＝sin_table_M(θ’-3M/4)
实施例二
M＝N，采用上述步骤402的方法预先设定了俯仰正弦函数表sin_table_N，由于M和N相等，因此水平正弦函数表sin_table_M与俯仰正弦函数表sin_table_N相同，本实施例中可以只获取一个三角函数表，例如sin_table，这样可以进一步减少三角函数表的计算量和存储量，节省存储空间。
可以按照以下规则获取位置信息的俯仰角信息对应的正弦函数值和余弦函数值：
可以按照以下规则获取位置信息的水平角信息对应的正弦函数值和余弦函数值：
sin(θ)＝sin_table(θ’)
sin(θ)＝sin_table(N/2-θ’)
sin(θ)＝-sin_table(θ’-N/2)
sin(θ)＝-sin_table(N-θ’)
cos(θ)＝sin_table(N/4-θ’)
cos(θ)＝-sin_table(θ’-N/4)
cos(θ)＝-sin_table(3N/4-θ’)
cos(θ)＝sin_table(θ’-3N/4)
实施例三
M≠N，M＝K1×N，K1≥2，采用上述步骤402的方法预先设定了水平正弦函数表sin_table_M，由于M和N具有倍数关系，因此俯仰正弦函数表sin_table_N可以通过前述倍数关系结合三角函数原理获得，因此本实施例中可以只获取一个三角函数表，例如sin_table_M，这样可以进一步减少三角函数表的计算量和存储量，节省存储空间。
可以按照以下规则获取位置信息的水平角信息对应的正弦函数值和余弦函数值：
sin(θ)＝sin_table_M(θ’)
sin(θ)＝sin_table_M(M/2-θ’)
sin(θ)＝-sin_table_M(θ’-M/2)
sin(θ)＝-sin_table_M(M-θ’)
cos(θ)＝sin_table_M(M/4-θ’)
cos(θ)＝-sin_table_M(θ’-M/4)
cos(θ)＝-sin_table_M(3M/4-θ’)
cos(θ)＝sin_table_M(θ’-3M/4)
可以按照以下规则获取位置信息的俯仰角信息对应的正弦函数值和余弦函数值：
实施例四
M≠N，N＝K2×M，K2≥2，采用上述步骤402的方法预先设定了俯仰正弦函数表sin_table_N，由于M和N具有倍数关系，因此水平正弦函数表sin_table_M可以通过前述倍数关系结合三角函数原理获得，因此本实施例中可以只获取一个三角函数表，例如sin_table_N，这样可以进一步减少三角函数表的计算量和存储量，节省存储空间。
可以按照以下规则获取位置信息的俯仰角信息对应的正弦函数值和余弦函数值：
可以按照以下规则获取位置信息的水平角信息对应的正弦函数值和余弦函数值：
sin(θ)＝sin_table_N(K2×θ’)
sin(θ)＝sin_table_N(N/2-K2×θ’)
sin(θ)＝-sin_table_N(K2×θ’-N/2)
sin(θ)＝-sin_table_N(N-K2×θ’)
cos(θ)＝sin_table_N(N/4-K2×θ’)
cos(θ)＝-sin_table_N(K2×θ’-N/4)
cos(θ)＝-sin_table_N(3N/4-K2×θ’)
cos(θ)＝sin_table_N(K2×θ’-3N/4)
实施例五
M≠N，采用上述步骤402的方法预先设定了俯仰余弦函数表cos_table_N和水平余弦函数表cos_table_M。
可以按照以下规则获取位置信息的俯仰角信息对应的正弦函数值和余弦函数值：
sin(θ)＝cos_table_N(N/4-θ’)
sin(θ)＝cos_table_N(θ’-N/4)
sin(θ)＝-cos_table_N(3N/4-θ’)
sin(θ)＝-cos_table_N(θ’-3N/4)
cos(θ)＝cos_table_N(θ’)
cos(θ)＝-cos_table_N(N/2-θ’)
cos(θ)＝-cos_table_N(θ’-N/2)
cos(θ)＝cos_table_N(N-θ’)
可以按照以下规则获取位置信息的水平角信息对应的正弦函数值和余弦函数值：
sin(θ)＝cos_table_M(M/4-θ’)
sin(θ)＝cos_table_M(θ’-M/4)
sin(θ)＝-cos_table_M(3M/4-θ’)
sin(θ)＝-cos_table_M(θ’-3M/4)
cos(θ)＝cos_table_M(θ’)
cos(θ)＝-cos_table_M(M/2-θ’)
cos(θ)＝-cos_table_M(θ’-M/2)
cos(θ)＝cos_table_M(M-θ’)
图7为本申请HOA系数的获取装置的一个示例性的结构图，如图7所示，该装置可以应用于上述实施例中的编码器20或解码器30。本实施例的HOA系数的获取装置可以包括：获取模块701和计算模块702。其中，
获取模块701，用于获取虚拟扬声器在预设球面上的位置信息，所述位置信息包括俯仰角信息和/或水平角信息；所述预设球面上包含M个经线圈和N个纬线圈，所述经线圈和所述纬线圈的交叉点称作参考点，位于同一个经线圈上的N个参考点等间距排列，位于同一个纬线圈上的M个参考点等间距排列；计算模块702，用于根据所述位置信息以及预设的参考三角函数表获取所述位置信息对应的三角函数值，所述参考三角函数表包括俯仰三角函数表和/或水平三角函数表；所述俯仰三角函数表包括多个第一参考点的俯仰角索引对应的多个三角函数值，所述多个第一参考点是第一经线圈上的参考点，所述多个第一参考点的数量不小于表示向下取整，所述第一经线圈是所述M个经线圈的其中之一；所述水平三角函数表包括多个第二参考点的水平角索引对应的多个三角函数值，所述多个第二参考点是第一纬线圈上的参考点，所述多个第二参考点的数量不小于所述第一纬线圈是所述N个纬线圈的其中之一；根据所述位置信息对应的三角函数值获取所述虚拟扬声器的高阶立体混响HOA系数。
在一种可能的实现方式中，所述计算模块702，具体用于根据所述位置信息的俯仰角信息以及所述俯仰三角函数表获取所述位置信息的俯仰角信息对应的三角函数值；和/或，根据所述位置信息的水平角信息以及所述水平三角函数表获取所述位置信息的水平角信息对应的三角函数值。
在一种可能的实现方式中，所述计算模块，具体用于当所述位置信息的俯仰角信息与所述多个第一参考点中的一个参考点的俯仰角索引对应上时，将所述俯仰三角函数表中所述一个参考点的俯仰角索引对应的三角函数值作为所述位置信息的俯仰角信息对应的三角函数值；或者，当所述位置信息的俯仰角信息与所述N个参考点中除所述多个第一参考点外的一个参考点的俯仰角索引对应上时，将所述N个参考点中除所述多个第一参考点外的一个参考点的俯仰角索引转换为所述多个第一参考点中的一个参考点的俯仰角索引，将所述俯仰三角函数表中所述一个参考点的俯仰角索引对应的三角函数值作为所述位置信息的俯仰角信息对应的三角函数值；或者，当所述位置信息的俯仰角信息与所述参考点的俯仰角索引均不对应时，根据所述位置信息的俯仰角信息获得与所述位置信息的俯仰角信息对应的参考点的俯仰角信息；当所述与所述位置信息的俯仰角信息对应的参考点的俯仰角信息与所述多个第一参考点中的一个参考点的俯仰角索引对应上时，将所述俯仰三角函数表中所述一个参考点的俯仰角索引对应的三角函数值作为所述位置信息的俯仰角信息对应的三角函数值；当所述与所述位置信息的俯仰角信息对应的参考点的俯仰角信息与所述N个参考点中除所述多个第一参考点外的一个参考点的俯仰角索引对应上时，将所述N个参考点中除所述多个第一参考点外的一个参考点的俯仰角索引转换为所述多个第一参考点中的一个参考点的俯仰角索引，将所述俯仰三角函数表中所述一个参考点的俯仰角索引对应的三角函数值作为所述位置信息的俯仰角信息对应的三角函数值。
在一种可能的实现方式中，所述计算模块，具体用于当所述位置信息的水平角信息与所述多个第二参考点中的一个参考点的水平角索引对应上时，将所述水平三角函数表中所述一个参考点的水平角索引对应的三角函数值作为所述位置信息的水平角信息对应的三角函数值；或者，当所述位置信息的水平角信息与所述M个参考点中除所述多个第二参考点外的一个参考点的水平角索引对应上时，将所述M个参考点中除所述多个第二参考点外的一个参考点的水平角索引转换为所述多个第二参考点中的一个参考点的水平角索引，将所述水平三角函数表中所述一个参考点的水平角索引对应的三角函数值作为所述位置信息的水平角信息对应的三角函数值；或者，当所述位置信息的水平角信息与所述参考点的水平角索引均不对应时，根据所述位置信息的水平角信息获得与所述位置信息的水平角信息对应的参考点的水平角信息；当所述与所述位置信息的水平角信息对应的参考点的水平角信息与所述多个第二参考点中的一个参考点的水平角索引对应上时，将所述水平三角函数表中所述一个参考点的水平角索引对应的三角函数值作为所述位置信息的水平角信息对应的三角函数值；当所述与所述位置信息的水平角信息对应的参考点的水平角信息与所述M个参考点中除所述多个第二参考点外的一个参考点的水平角索引对应上时，将所述M个参考点中除所述多个第二参考点外的一个参考点的水平角索引转换为所述多个第二参考点中的一个参考点的水平角索引，将所述水平三角函数表中所述一个参考点的水平角索引对应的三角函数值作为所述位置信息的水平角信息对应的三角函数值。
在一种可能的实现方式中，所述位置信息的俯仰角信息包括所述位置信息的俯仰角度或者所述位置信息的俯仰角索引；所述位置信息的水平角信息包括所述位置信息的水平角度或者所述位置信息的水平角索引。
在一种可能的实现方式中，所述三角函数值是指正弦函数值或者余弦函数值。
在一种可能的实现方式中，当所述多个第一参考点的俯仰角索引对应的多个三角函数为正弦函数值时，第i个第一参考点的俯仰角索引对应的正弦函数值满足如下公式(1)：
其中，i＝0,1,…,N’，ri表示所述第一经线圈的半径；
或，
当所述多个第一参考点的俯仰角索引对应的多个三角函数为余弦函数值时，第i个第一参考点的俯仰角索引对应的余弦函数值满足如下公式(2)：
其中，i＝0,1,…,N’，ri表示所述第一经线圈的半径。
在一种可能的实现方式中，当所述多个第二参考点的水平角索引对应的多个三角函数为正弦函数值时，第j个第二参考点的水平角索引对应的正弦函数值满足如下公式(3)：
其中，j＝0,1,…,M’，rj表示所述第一纬线圈的半径；
或，
当所述多个第二参考点的水平角索引对应的多个三角函数为余弦函数值时，第j个第二参考点的水平角索引对应的余弦函数值满足如下公式(4)：
其中，j＝0,1,…,M’，rj表示所述第一纬线圈的半径。
在一种可能的实现方式中，俯仰角度和俯仰角索引满足如下公式(5)：
其中，ri表示所述第一经线圈的半径，round()表示取整。
在一种可能的实现方式中，水平角度θ和水平角索引θ’满足如下公式(6)：
其中，rj表示所述第一纬线圈的半径，round()表示取整。
在一种可能的实现方式中，
当M≠N时，所述位置信息的俯仰角信息对应的三角函数值满足如下公式：
其中，表示所述位置信息的俯仰角信息，或者，表示所述与所述位置信息的俯仰角信息对应的参考点的俯仰角信息；sin_table_N()表示所述俯仰三角函数表；
当M≠N时，所述位置信息的水平角信息对应的三角函数值满足如下公式：
sin(θ)＝sin_table_M(θ’)
sin(θ)＝sin_table_M(M/2-θ’)
sin(θ)＝-sin_table_M(θ’-M/2)
sin(θ)＝-sin_table_M(M-θ’)
cos(θ)＝sin_table_M(M/4-θ’)
cos(θ)＝-sin_table_M(θ’-M/4)
cos(θ)＝-sin_table_M(3M/4-θ’)
cos(θ)＝sin_table_M(θ’-3M/4)
其中，θ表示所述位置信息的水平角信息，或者，θ表示所述与所述位置信息的水平角信息对应的参考点的水平角信息；sin_table_M()表示所述水平三角函数表；
或，
当M＝N时，所述位置信息的俯仰角信息对应的三角函数值满足如下公式：
其中，表示所述位置信息的俯仰角信息，或者，表示所述与所述位置信息的俯仰角信息对应的参考点的俯仰角信息；sin_table_N()表示所述俯仰三角函数表；
当M＝N时，所述位置信息的水平角信息对应的三角函数值满足如下公式：
sin(θ)＝sin_table(θ’)
sin(θ)＝sin_table(N/2-θ’)
sin(θ)＝-sin_table(θ’-N/2)
sin(θ)＝-sin_table(N-θ’)
cos(θ)＝sin_table(N/4-θ’)
cos(θ)＝-sin_table(θ’-N/4)
cos(θ)＝-sin_table(3N/4-θ’)
cos(θ)＝sin_table(θ’-3N/4)
其中，θ表示所述位置信息的水平角信息，或者，θ表示所述与所述位置信息的水平角信息对应的参考点的水平角信息；sin_table_N()表示所述俯仰三角函数表；
或，
当M≠N，M＝K1×N，K1≥2时，所述位置信息的水平角信息对应的三角函数值满足如下公式：
sin(θ)＝sin_table_M(θ’)
sin(θ)＝sin_table_M(M/2-θ’)
sin(θ)＝-sin_table_M(θ’-M/2)
sin(θ)＝-sin_table_M(M-θ’)
cos(θ)＝sin_table_M(M/4-θ’)
cos(θ)＝-sin_table_M(θ’-M/4)
cos(θ)＝-sin_table_M(3M/4-θ’)
cos(θ)＝sin_table_M(θ’-3M/4)
其中，θ表示所述位置信息的水平角信息，或者，θ表示所述与所述位置信息的水平角信息对应的参考点的水平角信息；sin_table_M()表示所述水平三角函数表；
当M≠N，M＝K1×N，K1≥2时，所述位置信息的俯仰角信息对应的三角函数值满足如下公式：
其中，表示所述位置信息的俯仰角信息，或者，表示所述与所述位置信息的俯仰角信息对应的参考点的俯仰角信息；sin_table_M()表示所述水平三角函数表；
或，
当M≠N，N＝K2×M，K2≥2时，所述位置信息的俯仰角信息对应的三角函数值满足如下公式：
其中，表示所述位置信息的俯仰角信息，或者，表示所述与所述位置信息的俯仰角信息对应的参考点的俯仰角信息；sin_table_N()表示所述俯仰三角函数表；
当M≠N，N＝K2×M，K2≥2时，所述位置信息的水平角信息对应的三角函数值满足如下公式：
sin(θ)＝sin_table_N(K2×θ’)
sin(θ)＝sin_table_N(N/2-K2×θ’)
sin(θ)＝-sin_table_N(K2×θ’-N/2)
sin(θ)＝-sin_table_N(N-K2×θ’)
cos(θ)＝sin_table_N(N/4-K2×θ’)
cos(θ)＝-sin_table_N(K2×θ’-N/4)
cos(θ)＝-sin_table_N(3N/4-K2×θ’)
cos(θ)＝sin_table_N(K2×θ’-3N/4)
其中，θ表示所述位置信息的水平角信息，或者，θ表示所述与所述位置信息的水平角信息对应的参考点的水平角信息；sin_table_N()表示所述俯仰三角函数表；
或，
当M≠N时，所述位置信息的俯仰角信息对应的三角函数值满足如下公式：
sin(θ)＝cos_table_N(N/4-θ’)
sin(θ)＝cos_table_N(θ’-N/4)
sin(θ)＝-cos_table_N(3N/4-θ’)
sin(θ)＝-cos_table_N(θ’-3N/4)
cos(θ)＝cos_table_N(θ’)
cos(θ)＝-cos_table_N(N/2-θ’)
cos(θ)＝-cos_table_N(θ’-N/2)
cos(θ)＝cos_table_N(N-θ’)
其中，表示所述位置信息的俯仰角信息，或者，表示所述与所述位置信息的俯仰角信息对应的参考点的俯仰角信息；cos_table_N()表示所述俯仰三角函数表；
当M≠N时，所述位置信息的水平角信息对应的三角函数值满足如下公式：
sin(θ)＝cos_table_M(M/4-θ’)
sin(θ)＝cos_table_M(θ’-M/4)
sin(θ)＝-cos_table_M(3M/4-θ’)
sin(θ)＝-cos_table_M(θ’-3M/4)
cos(θ)＝cos_table_M(θ’)
cos(θ)＝-cos_table_M(M/2-θ’)
cos(θ)＝-cos_table_M(θ’-M/2)
cos(θ)＝cos_table_M(M-θ’)
其中，θ表示所述位置信息的水平角信息，或者，θ表示所述与所述位置信息的水平角信息对应的参考点的水平角信息；cos_table_M()表示所述水平三角函数表。
本实施例的装置，可以用于执行图4所示方法实施例的技术方案，其实现原理和技术效果类似，此处不再赘述。
在实现过程中，上述方法实施例的各步骤可以通过处理器中的硬件的集成逻辑电路或者软件形式的指令完成。处理器可以是通用处理器、数字信号处理器(digital signalprocessor,DSP)、特定应用集成电路(application-specific integrated circuit，ASIC)、现场可编程门阵列(field programmable gate array,FPGA)或其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。本申请公开的方法的步骤可以直接体现为硬件编码处理器执行完成，或者用编码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器，闪存、只读存储器，可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器，处理器读取存储器中的信息，结合其硬件完成上述方法的步骤。
上述各实施例中提及的存储器可以是易失性存储器或非易失性存储器，或可包括易失性和非易失性存储器两者。其中，非易失性存储器可以是只读存储器(read-onlymemory，ROM)、可编程只读存储器(programmable ROM，PROM)、可擦除可编程只读存储器(erasable PROM，EPROM)、电可擦除可编程只读存储器(electrically EPROM，EEPROM)或闪存。易失性存储器可以是随机存取存储器(random access memory，RAM)，其用作外部高速缓存。通过示例性但不是限制性说明，许多形式的RAM可用，例如静态随机存取存储器(static RAM，SRAM)、动态随机存取存储器(dynamic RAM，DRAM)、同步动态随机存取存储器(synchronous DRAM，SDRAM)、双倍数据速率同步动态随机存取存储器(doubledata rateSDRAM，DDR SDRAM)、增强型同步动态随机存取存储器(enhanced SDRAM，ESDRAM)、同步连接动态随机存取存储器(synchlink DRAM，SLDRAM)和直接内存总线随机存取存储器(directrambus RAM，DR RAM)。应注意，本文描述的系统和方法的存储器旨在包括但不限于这些和任意其它适合类型的存储器。
本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本申请的范围。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统、装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
在本申请所提供的几个实施例中，应该理解到，所揭露的系统、装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本申请各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(个人计算机，服务器，或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(read-only memory，ROM)、随机存取存储器(random access memory，RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述，仅为本申请的具体实施方式，但本申请的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本申请揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本申请的保护范围之内。因此，本申请的保护范围应以所述权利要求的保护范围为准。

Claims (26)

1.一种音频接收设备实现的高阶立体混响HOA系数的获取方法，其特征在于，包括：
所述音频接收设备接收码流，并对所述码流进行解码以获取虚拟扬声器在预设球面上的属性信息，所述属性信息包括俯仰角信息和/或水平角信息；所述预设球面上包含M个经线圈和N个纬线圈，所述经线圈和所述纬线圈的交叉点称作参考点；
根据所述属性信息以及预设的参考三角函数表获取所述属性信息对应的三角函数值，所述参考三角函数表包括俯仰三角函数表和/或水平三角函数表；所述俯仰三角函数表包括多个第一参考点的俯仰角索引对应的多个三角函数值，所述多个第一参考点是第一经线圈上的参考点，所述多个第一参考点的数量不小于 表示向下取整，所述第一经线圈是所述M个经线圈的其中之一；所述水平三角函数表包括多个第二参考点的水平角索引对应的多个三角函数值，所述多个第二参考点是第一纬线圈上的参考点，所述多个第二参考点的数量不小于所述第一纬线圈是所述N个纬线圈的其中之一；
根据所述属性信息对应的三角函数值获取所述虚拟扬声器的HOA系数。
2.根据权利要求1所述的方法，其特征在于，所述根据所述属性信息以及预设的参考三角函数表获取所述属性信息对应的三角函数值，包括：
根据所述属性信息的俯仰角信息以及所述俯仰三角函数表获取所述属性信息的俯仰角信息对应的三角函数值；和/或，
根据所述属性信息的水平角信息以及所述水平三角函数表获取所述属性信息的水平角信息对应的三角函数值。
3.根据权利要求2所述的方法，其特征在于，所述根据所述属性信息的俯仰角信息以及所述俯仰三角函数表获取所述属性信息的俯仰角信息对应的三角函数值，包括：
当所述属性信息的俯仰角信息与所述多个第一参考点中的一个参考点的俯仰角索引对应上时，将所述俯仰三角函数表中所述一个参考点的俯仰角索引对应的三角函数值作为所述属性信息的俯仰角信息对应的三角函数值；或者，
当所述属性信息的俯仰角信息与所述N个参考点中除所述多个第一参考点外的一个参考点的俯仰角索引对应上时，将所述N个参考点中除所述多个第一参考点外的一个参考点的俯仰角索引转换为所述多个第一参考点中的一个参考点的俯仰角索引，将所述俯仰三角函数表中所述一个参考点的俯仰角索引对应的三角函数值作为所述属性信息的俯仰角信息对应的三角函数值；或者，
当所述属性信息的俯仰角信息与所述参考点的俯仰角索引均不对应时，根据所述属性信息的俯仰角信息获得与所述属性信息的俯仰角信息对应的参考点的俯仰角信息；当所述与所述属性信息的俯仰角信息对应的参考点的俯仰角信息与所述多个第一参考点中的一个参考点的俯仰角索引对应上时，将所述俯仰三角函数表中所述一个参考点的俯仰角索引对应的三角函数值作为所述属性信息的俯仰角信息对应的三角函数值；当所述与所述属性信息的俯仰角信息对应的参考点的俯仰角信息与所述N个参考点中除所述多个第一参考点外的一个参考点的俯仰角索引对应上时，将所述N个参考点中除所述多个第一参考点外的一个参考点的俯仰角索引转换为所述多个第一参考点中的一个参考点的俯仰角索引，将所述俯仰三角函数表中所述一个参考点的俯仰角索引对应的三角函数值作为所述属性信息的俯仰角信息对应的三角函数值。
4.根据权利要求2所述的方法，其特征在于，所述根据所述属性信息的水平角信息以及所述水平三角函数表获取所述属性信息的水平角信息对应的三角函数值，包括：
当所述属性信息的水平角信息与所述多个第二参考点中的一个参考点的水平角索引对应上时，将所述水平三角函数表中所述一个参考点的水平角索引对应的三角函数值作为所述属性信息的水平角信息对应的三角函数值；或者，
当所述属性信息的水平角信息与所述M个参考点中除所述多个第二参考点外的一个参考点的水平角索引对应上时，将所述M个参考点中除所述多个第二参考点外的一个参考点的水平角索引转换为所述多个第二参考点中的一个参考点的水平角索引，将所述水平三角函数表中所述一个参考点的水平角索引对应的三角函数值作为所述属性信息的水平角信息对应的三角函数值；或者，
当所述属性信息的水平角信息与所述参考点的水平角索引均不对应时，根据所述属性信息的水平角信息获得与所述属性信息的水平角信息对应的参考点的水平角信息；当所述与所述属性信息的水平角信息对应的参考点的水平角信息与所述多个第二参考点中的一个参考点的水平角索引对应上时，将所述水平三角函数表中所述一个参考点的水平角索引对应的三角函数值作为所述属性信息的水平角信息对应的三角函数值；当所述与所述属性信息的水平角信息对应的参考点的水平角信息与所述M个参考点中除所述多个第二参考点外的一个参考点的水平角索引对应上时，将所述M个参考点中除所述多个第二参考点外的一个参考点的水平角索引转换为所述多个第二参考点中的一个参考点的水平角索引，将所述水平三角函数表中所述一个参考点的水平角索引对应的三角函数值作为所述属性信息的水平角信息对应的三角函数值。
5.根据权利要求2-4中任一项所述的方法，其特征在于，所述属性信息的俯仰角信息包括所述属性信息的俯仰角度或者所述属性信息的俯仰角索引；所述属性信息的水平角信息包括所述属性信息的水平角度或者所述属性信息的水平角索引。
6.根据权利要求1-4中任一项所述的方法，其特征在于，所述三角函数值是指正弦函数值或者余弦函数值。
7.根据权利要求1-4中任一项所述的方法，其特征在于，当所述多个第一参考点的俯仰角索引对应的多个三角函数为正弦函数值时，第i个第一参考点的俯仰角索引对应的正弦函数值满足如下公式(1)：
其中，i＝0,1,…,N’，ri表示所述第一经线圈的半径；
或，
当所述多个第一参考点的俯仰角索引对应的多个三角函数为余弦函数值时，第i个第一参考点的俯仰角索引对应的余弦函数值满足如下公式(2)：
其中，i＝0,1,…,N’，ri表示所述第一经线圈的半径。
8.根据权利要求1-4中任一项所述的方法，其特征在于，当所述多个第二参考点的水平角索引对应的多个三角函数为正弦函数值时，第j个第二参考点的水平角索引对应的正弦函数值满足如下公式(3)：
其中，j＝0,1,…,M’，rj表示所述第一纬线圈的半径；
或，
当所述多个第二参考点的水平角索引对应的多个三角函数为余弦函数值时，第j个第二参考点的水平角索引对应的余弦函数值满足如下公式(4)：
其中，j＝0,1,…,M’，rj表示所述第一纬线圈的半径。
9.根据权利要求2所述的方法，其特征在于，俯仰角度和俯仰角索引满足如下公式(5)：
其中，ri表示所述第一经线圈的半径，round()表示取整。
10.根据权利要求2所述的方法，其特征在于，水平角度θ和水平角索引θ’满足如下公式(6)：
其中，rj表示所述第一纬线圈的半径，round()表示取整。
11.根据权利要求9或10所述的方法，其特征在于，当M≠N时，所述属性信息的俯仰角信息对应的三角函数值满足如下公式：
其中，表示所述属性信息的俯仰角信息，或者，表示所述与所述属性信息的俯仰角信息对应的参考点的俯仰角信息；sin_table_N()表示所述俯仰三角函数表；
当M≠N时，所述属性信息的水平角信息对应的三角函数值满足如下公式：
其中，θ表示所述属性信息的水平角信息，或者，θ表示所述与所述属性信息的水平角信息对应的参考点的水平角信息；sin_table_M()表示所述水平三角函数表；
或，
当M＝N时，所述属性信息的俯仰角信息对应的三角函数值满足如下公式：
其中，表示所述属性信息的俯仰角信息，或者，表示所述与所述属性信息的俯仰角信息对应的参考点的俯仰角信息；sin_table()表示三角函数表，sin_table()为sin_table_N()、sin_table_M()中的任一个；
当M＝N时，所述属性信息的水平角信息对应的三角函数值满足如下公式：
其中，θ表示所述属性信息的水平角信息，或者，θ表示所述与所述属性信息的水平角信息对应的参考点的水平角信息；sin_table()表示三角函数表，sin_table()为sin_table_N()、sin_table_M()中的任一个；
或，
当M≠N，M＝K1×N，K1≥2时，所述属性信息的水平角信息对应的三角函数值满足如下公式：
其中，θ表示所述属性信息的水平角信息，或者，θ表示所述与所述属性信息的水平角信息对应的参考点的水平角信息；sin_table_M()表示所述水平三角函数表；
当M≠N，M＝K1×N，K1≥2时，所述属性信息的俯仰角信息对应的三角函数值满足如下公式：
其中，表示所述属性信息的俯仰角信息，或者，表示所述与所述属性信息的俯仰角信息对应的参考点的俯仰角信息；sin_table_M()表示所述水平三角函数表；
或，
当M≠N，N＝K2×M，K2≥2时，所述属性信息的俯仰角信息对应的三角函数值满足如下公式：
其中，表示所述属性信息的俯仰角信息，或者，表示所述与所述属性信息的俯仰角信息对应的参考点的俯仰角信息；sin_table_N()表示所述俯仰三角函数表；
当M≠N，N＝K2×M，K2≥2时，所述属性信息的水平角信息对应的三角函数值满足如下公式：
其中，θ表示所述属性信息的水平角信息，或者，θ表示所述与所述属性信息的水平角信息对应的参考点的水平角信息；sin_table_N()表示所述俯仰三角函数表；
或，
当M≠N时，所述属性信息的俯仰角信息对应的三角函数值满足如下公式：
其中，表示所述属性信息的俯仰角信息，或者，表示所述与所述属性信息的俯仰角信息对应的参考点的俯仰角信息；cos_table_N()表示所述俯仰三角函数表；
当M≠N时，所述属性信息的水平角信息对应的三角函数值满足如下公式：
其中，θ表示所述属性信息的水平角信息，或者，θ表示所述与所述属性信息的水平角信息对应的参考点的水平角信息；cos_table_M()表示所述水平三角函数表。
12.根据权利要求1-4、9-10中任一项所述的方法，其特征在于，位于同一个所述经线圈上的N个参考点等间距排列，位于同一个所述纬线圈上的M个参考点等间距排列。
13.一种高阶立体混响HOA系数的获取装置，其特征在于，包括：
获取模块，用于接收码流，并对所述码流进行解码以获取虚拟扬声器在预设球面上的属性信息，所述属性信息包括俯仰角信息和/或水平角信息；所述预设球面上包含M个经线圈和N个纬线圈，所述经线圈和所述纬线圈的交叉点称作参考点，位于同一个经线圈上的N个参考点等间距排列，位于同一个纬线圈上的M个参考点等间距排列；
计算模块，用于根据所述属性信息以及预设的参考三角函数表获取所述属性信息对应的三角函数值，所述参考三角函数表包括俯仰三角函数表和/或水平三角函数表；所述俯仰三角函数表包括多个第一参考点的俯仰角索引对应的多个三角函数值，所述多个第一参考点是第一经线圈上的参考点，所述多个第一参考点的数量不小于 表示向下取整，所述第一经线圈是所述M个经线圈的其中之一；所述水平三角函数表包括多个第二参考点的水平角索引对应的多个三角函数值，所述多个第二参考点是第一纬线圈上的参考点，所述多个第二参考点的数量不小于所述第一纬线圈是所述N个纬线圈的其中之一；根据所述属性信息对应的三角函数值获取所述虚拟扬声器的HOA系数。
14.根据权利要求13所述的装置，其特征在于，所述计算模块，具体用于根据所述属性信息的俯仰角信息以及所述俯仰三角函数表获取所述属性信息的俯仰角信息对应的三角函数值；和/或，根据所述属性信息的水平角信息以及所述水平三角函数表获取所述属性信息的水平角信息对应的三角函数值。
15.根据权利要求14所述的装置，其特征在于，所述计算模块，具体用于当所述属性信息的俯仰角信息与所述多个第一参考点中的一个参考点的俯仰角索引对应上时，将所述俯仰三角函数表中所述一个参考点的俯仰角索引对应的三角函数值作为所述属性信息的俯仰角信息对应的三角函数值；或者，当所述属性信息的俯仰角信息与所述N个参考点中除所述多个第一参考点外的一个参考点的俯仰角索引对应上时，将所述N个参考点中除所述多个第一参考点外的一个参考点的俯仰角索引转换为所述多个第一参考点中的一个参考点的俯仰角索引，将所述俯仰三角函数表中所述一个参考点的俯仰角索引对应的三角函数值作为所述属性信息的俯仰角信息对应的三角函数值；或者，当所述属性信息的俯仰角信息与所述参考点的俯仰角索引均不对应时，根据所述属性信息的俯仰角信息获得与所述属性信息的俯仰角信息对应的参考点的俯仰角信息；当所述与所述属性信息的俯仰角信息对应的参考点的俯仰角信息与所述多个第一参考点中的一个参考点的俯仰角索引对应上时，将所述俯仰三角函数表中所述一个参考点的俯仰角索引对应的三角函数值作为所述属性信息的俯仰角信息对应的三角函数值；当所述与所述属性信息的俯仰角信息对应的参考点的俯仰角信息与所述N个参考点中除所述多个第一参考
点外的一个参考点的俯仰角索引对应上时，将所述N个参考点中除所述多个第一参考点外的一个参考点的俯仰角索引转换为所述多个第一参考点中的一个参考点的俯仰角索引，将所述俯仰三角函数表中所述一个参考点的俯仰角索引对应的三角函数值作为所述属性信息的俯仰角信息对应的三角函数值。
16.根据权利要求14所述的装置，其特征在于，所述计算模块，具体用于当所述属性信息的水平角信息与所述多个第二参考点中的一个参考点的水平角索引对应上时，将所述水平三角函数表中所述一个参考点的水平角索引对应的三角函数值作为所述属性信息的水平角信息对应的三角函数值；或者，当所述属性信息的水平角信息与所述M个参考点中除所述多个第二参考点外的一个参考点的水平角索引对应上时，将所述M个参考点中除所述多个第二参考点外的一个参考点的水平角索引转换为所述多个第二参考点中的一个参考点的水平角索引，将所述水平三角函数表中所述一个参考点的水平角索引对应的三角函数值作为所述属性信息的水平角信息对应的三角函数值；或者，当所述属性信息的水平角信息与所述参考点的水平角索引均不对应时，根据所述属性信息的水平角信息获得与所述属性信息的水平角信息对应的参考点的水平角信息；当所述与所述属性信息的水平角信息对应的参考点的水平角信息与所述多个第二参考点中的一个参考点的水平角索引对应上时，将所述水平三角函数表中所述一个参考点的水平角索引对应的三角函数值作为所述属性信息的水平角信息对应的三角函数值；当所述与所述属性信息的水平角信息对应的参考点的水平角信息与所述M个参考点中除所述多个第二参考点外的一个参考点的水平角索引对应上时，将所述M个参考点中除所述多个第二参考点外的一个参考点的水平角索引转换为所述多个第二参考点中的一个参考点的水平角索引，将所述水平三角函数表中所述一个参考点的水平角索引对应的三角函数值作为所述属性信息的水平角信息对应的三角函数值。
17.根据权利要求14-16中任一项所述的装置，其特征在于，所述属性信息的俯仰角信息包括所述属性信息的俯仰角度或者所述属性信息的俯仰角索引；所述属性信息的水平角信息包括所述属性信息的水平角度或者所述属性信息的水平角索引。
18.根据权利要求13-16中任一项所述的装置，其特征在于，所述三角函数值是指正弦函数值或者余弦函数值。
19.根据权利要求13-16中任一项所述的装置，其特征在于，当所述多个第一参考点的俯仰角索引对应的多个三角函数为正弦函数值时，第i个第一参考点的俯仰角索引对应的正弦函数值满足如下公式(1)：
其中，i＝0,1,…,N’，ri表示所述第一经线圈的半径；
或，
当所述多个第一参考点的俯仰角信息对应的多个三角函数为余弦函数值时，第i个第一参考点的俯仰角信息对应的余弦函数值满足如下公式(2)：
其中，i＝0,1,…,N’，ri表示所述第一经线圈的半径。
20.根据权利要求13-16中任一项所述的装置，其特征在于，当所述多个第二参考点的水平角索引对应的多个三角函数为正弦函数值时，第j个第二参考点的水平角索引对应的正弦函数值满足如下公式(3)：
其中，j＝0,1,…,M’，rj表示所述第一纬线圈的半径；
或，
当所述多个第二参考点的水平角信息对应的多个三角函数为余弦函数值时，其中，第j个第二参考点的水平角信息对应的余弦函数值满足如下公式(4)：
其中，j＝0,1,…,M’，rj表示所述第一纬线圈的半径。
21.根据权利要求14所述的装置，其特征在于，俯仰角度和俯仰角索引满足如下公式(5)：
其中，ri表示所述第一经线圈的半径，round()表示取整。
22.根据权利要求14所述的装置，其特征在于，水平角度θ和水平角索引θ’满足如下公式(6)：
其中，rj表示所述第一纬线圈的半径，round()表示取整。
23.根据权利要求21或22所述的装置，其特征在于，当M≠N时，所述属性信息的俯仰角信息对应的三角函数值满足如下公式：
其中，表示所述属性信息的俯仰角信息，或者，表示所述与所述属性信息的俯仰角信息对应的参考点的俯仰角信息；sin_table_N()表示所述俯仰三角函数表；
当M≠N时，所述属性信息的水平角信息对应的三角函数值满足如下公式：
其中，θ表示所述属性信息的水平角信息，或者，θ表示所述与所述属性信息的水平角信息对应的参考点的水平角信息；sin_table_M()表示所述水平三角函数表；
或，
当M＝N时，所述属性信息的俯仰角信息对应的三角函数值满足如下公式：
其中，表示所述属性信息的俯仰角信息，或者，表示所述与所述属性信息的俯仰角信息对应的参考点的俯仰角信息；sin_table()表示三角函数表，sin_table()为sin_table_N()、sin_table_M()中的任一个；
当M＝N时，所述属性信息的水平角信息对应的三角函数值满足如下公式：
其中，θ表示所述属性信息的水平角信息，或者，θ表示所述与所述属性信息的水平角信息对应的参考点的水平角信息；sin_table()表示三角函数表，sin_table()为sin_table_N()、sin_table_M()中的任一个，所述俯仰三角函数表包括多个第一参考点的俯仰角信息对应的多个正弦函数值；
或，
当M≠N，M＝K1×N，K1≥2时，所述属性信息的水平角信息对应的三角函数值满足如下公式：
其中，θ表示所述属性信息的水平角信息，或者，θ表示所述与所述属性信息的水平角信息对应的参考点的水平角信息；sin_table_M()表示所述水平三角函数表，所述水平三角函数表包括多个第二参考点的水平角信息对应的多个正弦函数值；
当M≠N，M＝K1×N，K1≥2时，所述属性信息的俯仰角信息对应的三角函数值满足如下公式：
其中，表示所述属性信息的俯仰角信息，或者，表示所述与所述属性信息的俯仰角信息对应的参考点的俯仰角信息；sin_table_M()表示所述水平三角函数表；
或，
当M≠N，N＝K2×M，K2≥2时，所述属性信息的俯仰角信息对应的三角函数值满足如下公式：
其中，表示所述属性信息的俯仰角信息，或者，表示所述与所述属性信息的俯仰角信息对应的参考点的俯仰角信息；sin_table_N()表示所述俯仰三角函数表；
当M≠N，N＝K2×M，K2≥2时，所述属性信息的水平角信息对应的三角函数值满足如下公式：
其中，θ表示所述属性信息的水平角信息，或者，θ表示所述与所述属性信息的水平角信息对应的参考点的水平角信息；sin_table_N()表示所述俯仰三角函数表；
或，
当M≠N时，所述属性信息的俯仰角信息对应的三角函数值满足如下公式：
其中，表示所述属性信息的俯仰角信息，或者，表示所述与所述属性信息的俯仰角信息对应的参考点的俯仰角信息；cos_table_N()表示所述俯仰三角函数表；
当M≠N时，所述属性信息的水平角信息对应的三角函数值满足如下公式：
其中，θ表示所述属性信息的水平角信息，或者，θ表示所述与所述属性信息的水平角信息对应的参考点的水平角信息；cos_table_M()表示所述水平三角函数表。
24.根据权利要求13-16、21-22中任一项所述的装置，其特征在于，位于同一个所述经线圈上的N个参考点等间距排列，位于同一个所述纬线圈上的M个参考点等间距排列。
25.一种音频处理设备，其特征在于，包括：
一个或多个处理器；
存储器，用于存储一个或多个程序；
当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如权利要求1-12中任一项所述的方法。
26.一种计算机可读存储介质，其特征在于，包括计算机程序，所述计算机程序在计算机上被执行时，使得所述计算机执行权利要求1-12中任一项所述的方法。
CN202310818175.2A 2021-03-05 2021-03-05 Hoa系数的获取方法和装置 Pending CN117278930A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310818175.2A CN117278930A (zh) 2021-03-05 2021-03-05 Hoa系数的获取方法和装置

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202310818175.2A CN117278930A (zh) 2021-03-05 2021-03-05 Hoa系数的获取方法和装置
CN202110246382.6A CN115038027B (zh) 2021-03-05 2021-03-05 Hoa系数的获取方法和装置

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
CN202110246382.6A Division CN115038027B (zh) 2021-03-05 2021-03-05 Hoa系数的获取方法和装置

Publications (1)

Publication Number Publication Date
CN117278930A true CN117278930A (zh) 2023-12-22

Family

ID=83118433

Family Applications (4)

Application Number Title Priority Date Filing Date
CN202310818006.9A Pending CN117395591A (zh) 2021-03-05 2021-03-05 Hoa系数的获取方法和装置
CN202110246382.6A Active CN115038027B (zh) 2021-03-05 2021-03-05 Hoa系数的获取方法和装置
CN202310821008.3A Pending CN117041856A (zh) 2021-03-05 2021-03-05 Hoa系数的获取方法和装置
CN202310818175.2A Pending CN117278930A (zh) 2021-03-05 2021-03-05 Hoa系数的获取方法和装置

Family Applications Before (3)

Application Number Title Priority Date Filing Date
CN202310818006.9A Pending CN117395591A (zh) 2021-03-05 2021-03-05 Hoa系数的获取方法和装置
CN202110246382.6A Active CN115038027B (zh) 2021-03-05 2021-03-05 Hoa系数的获取方法和装置
CN202310821008.3A Pending CN117041856A (zh) 2021-03-05 2021-03-05 Hoa系数的获取方法和装置

Country Status (9)

Country Link
US (1) US20230421978A1 (zh)
EP (1) EP4294047A4 (zh)
JP (1) JP2024509179A (zh)
KR (1) KR20230152139A (zh)
CN (4) CN117395591A (zh)
BR (1) BR112023017357A2 (zh)
CA (1) CA3212643A1 (zh)
MX (1) MX2023010359A (zh)
WO (1) WO2022184096A1 (zh)

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3217653B1 (en) * 2009-12-24 2023-12-27 Nokia Technologies Oy An apparatus
US9412385B2 (en) * 2013-05-28 2016-08-09 Qualcomm Incorporated Performing spatial masking with respect to spherical harmonic coefficients
US9466305B2 (en) * 2013-05-29 2016-10-11 Qualcomm Incorporated Performing positional analysis to code spherical harmonic coefficients
WO2017157803A1 (en) * 2016-03-15 2017-09-21 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Apparatus, method or computer program for generating a sound field description
US10015618B1 (en) * 2017-08-01 2018-07-03 Google Llc Incoherent idempotent ambisonics rendering
US10674301B2 (en) * 2017-08-25 2020-06-02 Google Llc Fast and memory efficient encoding of sound objects using spherical harmonic symmetries
BR112020011026A2 (pt) * 2017-11-17 2020-11-17 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e. V. aparelho e método para codificar ou decodificar parâmetros de codificação de áudio direcional com o uso de quantização e codificação de entropia
CN111542877B (zh) * 2017-12-28 2023-11-24 诺基亚技术有限公司 空间音频参数编码和相关联的解码的确定

Also Published As

Publication number Publication date
CN115038027B (zh) 2023-07-07
BR112023017357A2 (pt) 2023-10-03
MX2023010359A (es) 2023-09-18
CN117395591A (zh) 2024-01-12
US20230421978A1 (en) 2023-12-28
CN115038027A (zh) 2022-09-09
WO2022184096A1 (zh) 2022-09-09
CA3212643A1 (en) 2022-09-09
EP4294047A1 (en) 2023-12-20
CN117041856A (zh) 2023-11-10
KR20230152139A (ko) 2023-11-02
EP4294047A4 (en) 2024-08-14
JP2024509179A (ja) 2024-02-29

Similar Documents

Publication Publication Date Title
CN111542877B (zh) 空间音频参数编码和相关联的解码的确定
US11558707B2 (en) Sound field adjustment
JP2023551040A (ja) オーディオの符号化及び復号方法及び装置
CN115176486A (zh) 利用空间元数据内插的音频渲染
CN114582357A (zh) 一种音频编解码方法和装置
CN114067810A (zh) 音频信号渲染方法和装置
CN116671132A (zh) 利用空间元数据内插和源位置信息的音频渲染
CN106165451B (zh) 对高阶高保真立体声信号应用动态范围压缩的方法和设备
US20230412981A1 (en) Method and apparatus for determining virtual speaker set
CN117501362A (zh) 音频渲染系统、方法和电子设备
CN115346537A (zh) 一种音频编码、解码方法及装置
CN115038027B (zh) Hoa系数的获取方法和装置
KR20220047821A (ko) 공간 오디오 방향 파라미터의 양자화
CN117546236A (zh) 音频渲染系统、方法和电子设备
WO2020089509A1 (en) Determination of spatial audio parameter encoding and associated decoding
WO2024212894A1 (zh) 场景音频信号的解码方法和装置
WO2024149557A1 (en) A method and apparatus for complexity reduction in 6dof audio rendering
Favilla et al. Acoustic sound localisation: visualisations of a 1st order ambisonic microphone array
CN118696374A (zh) 使用经优化球面量化字典对球面坐标进行编码和解码

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