JP2021124439A - Collected sound data display program - Google Patents

Abstract

To provide, as means for solving the problem of prior art, a collected sound data display program, with which it is possible to easily and effectively visualize the arrival direction of the sound obtained by a 3D microphone in order to evaluate sound environments.SOLUTION: A collected sound display program pertaining to the present invention analyzes collected sound data having been collected for a prescribed period by a 3D microphone, the collected sound display program comprising a function that causes an intensity calculation process, a unit angular strength calculation process and intensity distribution map creation process to be executed by a computer. In the intensity distribution map creation process out of these, the polar coordinates of arrival point of a unit angular strength are converted into a coordinate system composed of altitudes and longitudes set to a spherical surface, and an "intensity distribution map" is created in which the intensity of unit angular strength is distributed to a coordinate system composed of altitudes and longitudes.SELECTED DRAWING: Figure 1

Description

The present invention is a technique related to sound environment evaluation, and more specifically, in order to evaluate the sound environment, it is possible to analyze data collected by a microphone having a plurality of channels and effectively visualize the result. It is related to the sound collection data display program.

The conventional sound environment evaluation has been mainly performed by using the measurement result by the noise level meter. This noise level meter is equipped with one omnidirectional microphone, and although it can measure the sound pressure level existing in the environment, the sound source direction (arrival direction) of the sound and the sound source type of the sound. (Noise from construction, noise from driving a car, noise from home appliances, etc.) cannot be measured. Therefore, when trying to implement noise countermeasures, it was not possible to appropriately plan the location where noise countermeasures should be implemented and the degree of noise countermeasures based only on the measurement results of the noise level meter.

By the way, in recent years, attention has been paid to acoustic technology for three-dimensionally reproducing a sound field. Ambisonics, which was developed by an audio engineer, is one of them. For example, the sound field can be changed as appropriate according to the rotation of the listener's head, and by combining it with video, it is possible to provide realistic content. It is possible. Therefore, technologies using Ambisonics have been proposed in various fields. For example, in Patent Document 1, the Ambisonics technology is used to reproduce a sound field corresponding to each of a plurality of persons in a vehicle. We are proposing a method to do this.

Japanese Unexamined Patent Publication No. 2016-220032

Generally, for sound collection in Ambisonics, a 3D microphone (also referred to as an Ambisonic microphone) having a plurality of individual microphones having known directions is used. For example, four individual microphones (that is, from the inside of the regular tetrahedron toward each wall surface), an individual microphone facing forward and upper left from the center, an individual microphone pointing forward to the lower right, an individual microphone pointing to the lower left rear, and an individual microphone pointing to the upper right rear. A 3D microphone can be constructed by means of four individual microphones arranged in a row.

By analyzing the recording (sound collection data) collected by the 3D microphone, the direction of arrival of the sound can be obtained. That is, the direction of arrival of sound, which could not be measured by the noise level meter commonly used in the conventional sound environment evaluation, can be measured by the 3D microphone. Therefore, it is conceivable to use the sound collection data by the 3D microphone in order to evaluate the sound environment.

However, the sound environment evaluation using a 3D microphone has not been known so far. In addition, the direction of arrival of sound can be obtained by analyzing the sound collection data by a 3D microphone. Although there is an example of visualizing the direction of arrival of sound, the conventional method uses an expensive system and intuitively obtains sound. It was not possible to grasp the direction of arrival. Therefore, it has been desired to develop a method for easily and effectively visualizing the sound environment evaluation.

The object of the present invention is to solve the problem of the prior art, that is, to evaluate the sound environment, it is possible to easily and effectively visualize the arrival direction of the sound obtained by the 3D microphone. To provide a program.

The present invention has been made by paying attention to the point that a reference sphere centered on a sound receiving point is set and an intensity distribution map expressing intensity is created by using this reference sphere. It is an invention made based on an unprecedented idea.

The sound collection data display program of the present invention is a program that causes a computer to execute a function of analyzing sound collection data collected for a predetermined period by a 3D microphone and displaying the analysis result using a reference spherical surface. It has a function to make a computer execute city calculation processing, unit angle strength calculation processing, and intensity distribution map creation processing. The "reference spherical surface" here is a spherical surface set with a predetermined radius around the sound receiving point. In the intensity calculation process, the "instantaneous intensity" for each sound collection time is calculated from the sound collection data, and in the unit angle intensity calculation process, it is instantaneously divided into a divided area (a small area obtained by dividing the reference sphere into a plurality of areas). In addition to assigning the intensity of, the "unit angle intensity" is obtained by dividing the sound intensity of the instantaneous intensity by the solid angle of this divided region. In the intensity distribution map creation process, an "intensity distribution map" in which the unit angle intensity is distributed on the reference sphere is created.

The sound collection data display program of the present invention converts the divided region into a Cartesian coordinate system consisting of latitude and longitude set on the reference sphere, and also creates an "intensity distribution map" in which the unit angle intensity is distributed in this Cartesian coordinate system. It can also be created.

The sound collection data display program of the present invention may further include a function of causing a computer to execute an intensity calculation process for each arrival angle. In the strength calculation process for each arrival angle, the latitude and longitude are set on the reference sphere, and the "strength for each arrival angle", which is the unit angle strength for each unit longitude width, is calculated. In this case, an "intensity distribution map" is created in which the intensities for each arrival angle are distributed for each direction on the reference cross section for cutting the reference spherical surface. Further, in the intensity calculation process for each arrival angle, the intensity for each arrival angle can be calculated based on the instantaneous intensity within a designated latitude range.

The sound collection data display program of the present invention may further include a function of causing a computer to execute a sound source type estimation process. In the sound source type estimation process, the sound source type of intensity is estimated based on the time change of the intensity and the sound intensity of the intensity and the sound source type information stored in advance. The sound source type information includes the sound source type, the direction of arrival of the sound, the strength of the sound, and the time change of the sound. In this case, in the intensity distribution map creation process, an “intensity distribution map” is created in which the unit angle intensity based on the instantaneous intensity at which the sound source type is estimated by the sound source type estimation process is distributed. In the sound source type estimation process, in addition to (or instead of) the time change of the sound intensity of the intensity, the sound source type of the intensity can be estimated based on the time change in the direction of arrival of the intensity.

The sound collection data display program of the present invention may further include a function of causing a computer to execute a specific sound detection process. In the specific sound detection process, the instantaneous intensity in a specific arrival direction range, which is the instantaneous intensity in which the sound intensity exceeds a predetermined intensity threshold value, is detected. In this case, in the intensity distribution map creation process, an "intensity distribution map" in which the unit angle intensity based on the instantaneous intensity detected by the specific sound detection is distributed is created. Further, in the specific sound detection process, it is also possible to detect a predetermined instantaneous intensity on the condition of a specific frequency range.

The sound collection data display program of the present invention may further include a function of causing a computer to execute an intensity extraction process. In the intensity extraction process, of the two or more instantaneous intensities obtained from the sound collection data collected for a predetermined period, the instantaneous intensity is in a specific arrival direction range and the sound intensity exceeds a predetermined intensity threshold value. Extract the city. In this case, in the intensity distribution map creation process, the intensity distribution map at the pick-up time of the instantaneous intensity extracted by the intensity extraction process is created.

The sound collection data display program of the present invention calculates the "unit time intensity", which is the sound intensity of the instantaneous intensity in a predetermined time width, and executes a process of obtaining the unit angle intensity based on the unit time intensity on the computer. It can also be made to.

The sound collection data display program of the present invention may also cause a computer to execute a process of calculating an instantaneous intensity based on sound collection data extracted according to a preset frequency condition.

The sound collection data display program of the present invention may further include a function of causing a computer to execute a time series display process. In the time-series display process, the intensity distribution map at a plurality of sound-collecting times is continuously displayed in the order of the sound-collecting times.

The sound collecting data display program of the present invention has the following effects.
(1) The sound arrival direction can be grasped by using the sound collection data by the 3D microphone, and the sound arrival direction and the sound intensity can be intuitively grasped by the intensity distribution map. As a result. Appropriate sound environment evaluation can be performed, and effective noise countermeasures can be planned.
(2) If the arrival direction and frequency to be focused on are determined in advance, a specific sound (for example, a harmful sound requiring countermeasures) can be detected by comparing with the intensity threshold value.
(3) If the characteristics of sound such as ceiling ventilation and printers in the office (direction of arrival, sound intensity, and changes over time) are known in advance, the sound source type of sound collection data can be estimated. As a result, the sound to be dealt with can be efficiently detected.

The block diagram which shows the relationship of the main processing of the sound collection data display program of this invention. The flow chart which shows the main processing flow of the sound collection data display program of this invention. (A) is a mathematical expression diagram showing sound collection data in A format, and (b) is a mathematical expression diagram showing sound collection data in B format format. A spectrogram diagram showing the instantaneous intensity. A mathematical diagram for obtaining the "sound intensity" and "sound arrival direction" of the intensity based on the sound collection data in the B format. (A) is a model diagram showing an arbitrary vector on the reference sphere, (b) is a plan view in which the vector is projected on the equatorial plane (z = 0) of the reference sphere, and (c) is a vertical plane (z axis) including the vector. A cross-sectional view in which a vector is projected onto a cross section (including a cross section). A model diagram in which a large number of divided region MSs are set on the reference spherical surface BF. A model diagram explaining the arrival point. A model diagram showing an example of a Cartesian coordinate system intensity distribution diagram. A mathematical diagram for obtaining the sound intensity of the intensity in _{a predetermined range (azimuth θ 1 to θ 2} , elevation angle φ _{1 to φ 2} ) when creating a Cartesian coordinate system intensity distribution map. (A) is a model diagram showing an example of a polar coordinate system intensity distribution map, and (b) is a model diagram explaining a procedure for creating a polar coordinate system intensity distribution map. This is a mathematical formula for calculating the sound intensity for each arrival angle when creating a polar coordinate system intensity distribution map. When the XY plane is the equatorial plane, when the XY plane is the equatorial plane, the YZ A mathematical diagram showing each of the cases where the plane is the equatorial plane. The flow diagram which shows the flow of the process which the sound source type estimation process estimates a sound source type. The flow diagram which shows the flow of the process which the specific sound detection process detects a specific sound.

An example of the embodiment of the sound collection data display program of the present invention will be described with reference to the drawings.

1. 1. Overall Overview The sound collection data display program of the present invention analyzes a record of sound collection for a predetermined period (hereinafter referred to as "sound collection data") using a 3D microphone, and the sound collection position (hereinafter, "sound reception point"). ”), So to speak, the process of creating a sound distribution map is executed by a computer. By showing the distribution map, it is possible to intuitively grasp from which direction and how loud the sound is coming from the sound receiving point. For example, soundproofing that is effective when taking noise countermeasures. The facility can be deployed.

Here, the 3D microphone is a multi-channel microphone having a plurality of individual microphones whose directions are known. A 3D microphone can be composed of four individual microphones (that is, four individual microphones arranged from the inside of a regular tetrahedron toward each wall surface) of individual microphones facing rear upper right. The sound pick-up data picked up by the 3D microphone is composed of information for each individual microphone, such as the identifier of the individual microphone (which individual microphone picked up the sound), the time when the sound was picked up, and the frequency and pressure of the sound. Contains information.

The main processing of the sound collecting data display program of the present invention will be described with reference to FIG. FIG. 1 is a block diagram showing a relationship between main processes of the sound collecting data display program 100 of the present invention. As shown in this figure, the sound collection data display program 100 includes a function of causing a computer to execute an intensity calculation process 103, a unit angle intensity calculation process 104, and an intensity distribution map creation process 111, and further detects specific sound. It can also be configured to include a function of causing a computer to execute the process 107, the sound source type estimation process 109, the time series display process 112, and the intensity extraction process 113. Further, a sound collecting data display including a 3D microphone 101, a sound collecting data storage means 102, a frequency condition setting means 105, a unit angle intensity storage means 106, a focus range storage means 108, a sound characteristic storage means 110, and an output means such as a display or a printer. It can also be built as a system.

The computer device that executes each function of the sound collection data display program 100 of the present invention is a personal computer (PC), a tablet PC such as an iPad (registered trademark), a mobile terminal including a smartphone, or a PDA (Personal Data Assistance). Etc. can be used. The computer device includes a processor such as a CPU and a memory such as a ROM and a RAM, and some computer devices further include an input means such as a mouse and a keyboard and a display. Further, the sound collecting data storage means 102, the unit angle intensity storage means 106, the focus range storage means 108, and the sound characteristic storage means 110 can be constructed on a database server, for example, and can be built on a local network (LAN: Local Area Network). It can be placed, or it can be a cloud server that stores via the Internet (that is, wireless communication or wired communication).

Subsequently, the main processing of the sound collecting data display program 100 of the present invention will be described with reference to FIG. FIG. 2 is a flow chart showing the main processing flow of the sound collection data display program 100. The central column shows the processing to be performed, the left column shows the input information required for the processing, and the right column shows the input information required for the processing. The output information generated from the processing is shown.

As shown in FIG. 2, first, the "instantaneous intensity" is calculated (Step 10). Specifically, by causing the computer to execute the intensity calculation process 103, the sound collection data is read from the sound collection data storage means 102, and the instantaneous intensity is calculated based on the sound collection data. The instantaneous intensity calculated here is the result for each sound collection time, and consists of the sound arrival direction and the sound intensity.

When the instantaneous intensity is calculated, the sound source type can be estimated by causing the computer to execute the sound source type estimation process 109 (Step 20). For example, if the characteristics of sound sources (noise from construction, noise from driving a car, sound from home appliances, etc.) existing around the sound receiving point are known in advance (hereinafter, simply referred to as "sound characteristics"), it will be instantaneous. It is possible to estimate the sound source type by extracting the equipment corresponding to the known equipment from the intensities of. As a result, when taking measures against noise, construction noise and car running noise are certified as sounds that should be dealt with, and the sounds of ceiling ventilation and printers installed in the office are certified as sounds that do not need to be dealt with. be able to. Note that the sound characteristics referred to here include information such as the arrival direction, sound intensity, time change in the arrival direction, and time change in sound intensity.

Furthermore, if a sound source that is likely to require countermeasures such as construction noise or automobile running noise is known in advance, a sound that requires countermeasures by having a computer execute the specific sound detection process 107 (hereinafter referred to as "specific sound"). ) Can also be detected (Step 30). When taking noise countermeasures in this way, it is sufficient to implement soundproofing measures only when the construction noise or automobile running noise that has been paid attention to in advance exceeds the allowable sound intensity, that is, efficient soundproofing measures can be taken. ..

When the instantaneous intensity is calculated, the unit angle intensity calculation process 104 is executed by the computer to calculate the "unit angle intensity" based on the instantaneous intensity (Step 40). The unit angle intensity calculated here is a value obtained by assigning an instantaneous intensity to the "division area" set in the "reference sphere" and dividing the sound intensity of the instantaneous intensity by the solid angle of the division area. Is. The "reference spherical surface" referred to here is a spherical surface having a predetermined radius set around the sound receiving point (the size of this predetermined radius can be appropriately designed), and the "divided region" is defined as a "divided area". , A divided region (so-called mesh) obtained by dividing a reference spherical surface into a plurality of parts.

When the unit angle intensity is calculated, the intensity distribution map creation process 111 is executed by a computer to create an "intensity distribution map" based on the unit angle intensity (Step 50). When the sound source type is estimated by the sound source type estimation process 109 (Step 20), the intensity distribution map can be created based on the unit angle intensity related to the sound for which the sound source type is estimated, or the specific sound detection process 107 can be used. When detecting a specific sound (Step 30), an intensity distribution map can be created based on the unit angle intensity related to the detected sound. The intensity distribution map can be created for each sound collection time, and is a three-dimensional intensity distribution map (hereinafter referred to as "spherical intensity distribution map") represented on a spherical surface and a diagram showing the intensity distribution. Intensity distribution map represented on the plane of the Cartesian coordinate system (hereinafter referred to as "Cartesian coordinate system intensity distribution map"), intensity distribution map represented on the plane of the polar coordinate system (hereinafter referred to as "curved coordinate system intensity"). It is classified into three types: "city distribution map").

When the intensity distribution map for each sound collection time is created, the intensity distribution map can be continuously displayed by moving images or frame-by-frame by causing a computer to execute the time-series display process 112 (Step 60). Specifically, the intensity distribution map created at a plurality of sound collection times is continuously displayed in the order of the sound collection time for each preset display time. This display time can be set at relatively short intervals such as a few seconds interval or a few minutes interval, or can be set at a relatively long period such as an hour unit, a week unit, or a monthly unit.

Hereinafter, each process to be executed by the computer by the sound collecting data display program 100 will be described in detail.

2. Intensity calculation process As described above, the intensity calculation process 103 reads the sound collection data from the sound collection data storage means 102, and based on the sound collection data, the intensity for each sound collection time (hereinafter, for each sound collection time). Intensity is called "instantaneous intensity"). The sound collection data stored by the sound collection data storage means 102 is the sound collected by the 3D microphone 101. Normally, the sound collection data collected by the 3D microphone 101 is data in a format called A format, and when the computer executes the intensity calculation process 103, the sound collection data in A format is converted to B format by the encoder. Convert. 3A and 3B are mathematical formulas representing sound collection data of each format, in which FIG. 3A shows an A format format and FIG. 3B shows a B format format. However, this figure shows a case of a 3D microphone 101 composed of four individual microphones arranged from the inside of a regular tetrahedron toward each wall surface.

The instantaneous intensity obtained as a result of having the computer execute the intensity calculation process 103 using the sound intensity data converted into the B format as the input data is information for each sound intensity time and for each frequency. It is a vector consisting of the direction of arrival of sound and the intensity of sound. Since the sound intensity is obtained for each sound pick-up time and each frequency, the so-called spectrogram as shown in FIG. 4 can represent the instantaneous intensity. In this figure, the sound intensity of each sound collection data is represented by color in a coordinate system in which the time is on the horizontal axis and the frequency is on the vertical axis. The "sound arrival direction" and "sound intensity" of the intensity can be calculated from the sound collection data in the B format format as shown in FIG. 5, respectively. In FIG. 5, a mathematical formula for obtaining the “sound intensity” of the intensity based on the sound intensity data in the B format format and a mathematical formula for obtaining the “sound arrival direction” of the intensity based on the sound intensity data in the B format format. Is shown.

Here, the azimuth angle θ and the elevation angle φ used in the mathematical formula of FIG. 5 will be described with reference to FIG. FIG. 6A is a model diagram showing an arbitrary vector VC on a reference spherical surface BF set with a radius r centered on the sound receiving point SP, and FIG. 6B is an equatorial plane (Z = 0) of the reference spherical surface BF. ) Is a plan view in which the vector VC is projected, and FIG. 6 (c) is a cross-sectional view in which the vector VC is projected on a vertical surface (a cross section including the Z axis) including the vector VC. The azimuth angle θ is a horizontal angle formed by a predetermined reference direction (X-axis in this figure) and a vector VC as shown in FIG. 6 (b). One elevation angle φ is the vertical perpendicularity formed by the XY plane (or Z axis) and the vector VC as shown in FIG. 6 (c).

3. 3. Unit angle intensity calculation process As described above, in the unit angle intensity calculation process 104, the unit angle intensity is calculated based on the instantaneous intensity, and this unit angle intensity divides the sound intensity of the instantaneous intensity into the solid angle of the divided region. It is the value divided by the angle. Hereinafter, the procedure for calculating the “unit angle strength” will be described in detail with reference to FIG. 7.

FIG. 7 is a model diagram in which a large number of divided region MSs are set on the reference spherical surface BF. For convenience, FIG. 7 shows a partially cut reference spherical surface BF. The divided region MS shown in this figure is a small region surrounded by a plurality of horizontal lines parallel to the XY plane (equatorial plane) and a plurality of vertical lines connecting both poles (upper and lower single ends of the Z axis), that is, a plurality of regions. It is an area surrounded by the plumb bob. The division region MS can be set by various division methods, not limited to the mode shown in FIG. 7, as long as a plurality of small regions are set on the reference spherical surface BF.

When the divided area MS is set, the instantaneous intensity obtained by causing the computer to execute the intensity calculation process 103 is assigned to the divided area MS. Specifically, as shown in FIG. 7, the instantaneous intensity IT is assigned to the divided region MS in which the instantaneous intensity IT (that is, the vector starting from the sound source) passes through the reference spherical surface BF. Note that FIG. 7 shows only one instantaneous intensity IT for convenience, but of course, instantaneous intensity IT is required for the sound collection time, and all (or a desired part) of these are instantaneous. Intensity IT is assigned to the divided area MS. Then, the "unit angle intensity" is calculated by dividing the sound intensity of the instantaneous intensity IT by the solid angle of the divided region.

When allocating the instantaneous intensity to the divided region MS, it is also possible to obtain and allocate the “arrival point AP” shown in FIG. FIG. 8 is a model diagram for explaining the arrival point AP. As shown in this figure, the arrival point A is the intersection of the instantaneous intensity and the reference sphere BF. In this figure, the front-back direction of the sound receiving point SP is the X-axis (front is positive), the left-right direction of the sound-receiving point SP is the Y-axis (the left is positive), and the vertical direction of the sound-receiving point SP is the Z-axis (upper). Since the reference sphere BF is set in the 3-axis coordinate system in which is positive), the coordinates of the arrival point AP can be represented by the 3-axis coordinate system (X, Y, Z), but the present invention is not limited to this. It can also be expressed in a polar coordinate system based on the radius and angle (azimuth angle θ and elevation angle φ).

The unit angle intensity is obtained based on the instantaneous intensity, and the instantaneous intensity is calculated for each sound intensity and frequency. In other words, the instantaneous intensity is called the sound intensity and frequency. It also has information. The unit angle intensity calculated by causing the computer to execute the unit angle intensity calculation process 104 is stored in the unit angle intensity storage means 106 together with this information (FIG. 1).

The unit angle intensity can be obtained for each sound collection time, that is, based on the instantaneous intensity obtained for each sound collection time, or can be obtained within a predetermined time width (hereinafter, referred to as "unit time"). It can also be calculated based on the instantaneous intensity. In this case, the unit time is set by specifying a desired time zone (start time to end time) or by specifying a fixed interval (for example, 1 hour or 1 day) in advance, and the unit time is obtained within the unit time. It is preferable to calculate the unit angle intensity by synthesizing the sound intensity of the instantaneous intensity (integrated value or total value) and dividing the combined sound intensity by the division region MS.

By the way, the intensity distribution map is created based on the unit angle intensity, but in some cases, the intensity distribution map is limited to a specific frequency range or excluded from the specific frequency range. Sometimes I want to create it. in this case. The operator sets a desired frequency range using the frequency condition setting means 105 (FIG. 1). Then, the intensity calculation process 103 extracts the sound collection data in the frequency range set by the frequency condition setting means 105 from the sound collection data, and creates an instantaneous intensity based on the sound collection data. As a result, the unit angle intensity corresponding to the desired frequency range is obtained, and the intensity distribution map corresponding to the desired frequency range is created.

4. Intensity distribution map creation process As described above, the intensity distribution map creation process 111 creates an intensity distribution map based on the unit angle intensity. This intensity distribution map is, so to speak, a diagram showing the distribution of intensity intensity, and is classified into three types: a spherical intensity distribution map, a Cartesian coordinate system intensity distribution map, and a polar coordinate system intensity distribution map. Therefore, each of the various intensity distribution maps will be described in detail.

(Spherical intensity distribution map)
The spherical intensity distribution map is a distribution map obtained by imparting a unit angle intensity to the divided region MS set on the reference spherical surface BF, and is a distribution map that three-dimensionally expresses the distribution of the unit angle intensity. be. The unit angle intensity can be expressed by a color model (for example, RGB or CMYK) consisting of three attributes of hue, saturation, and lightness according to the intensity (large or small value), and is also expressed by shades such as gray scale. It can also be represented by a combination of color models and shades.

(Cartesian coordinate system intensity distribution map)
FIG. 9 is a model diagram showing an example of a Cartesian coordinate system intensity distribution diagram. The Cartesian coordinate system intensity distribution map is a distribution map showing the unit angle strength on the plane of the Cartesian coordinate system. Specifically, the divided region MS is set in the Cartesian coordinate system (from the latitude and longitude set in the reference sphere BF). It is a distribution map obtained by giving a unit angle strength to the divided region MS which has been converted into the coordinate system (coordinate system) and which has been subjected to the coordinate conversion. In FIG. 9, the unit angle intensity is represented by color in the Cartesian coordinate system with the azimuth angle θ as the horizontal axis and the elevation angle φ as the vertical axis. Can be represented by a color model, shades, or a combination of these.

Further, the unit angle intensity shown in the Cartesian coordinate system intensity distribution map can be expressed as a point (that is, for each division region MS) or can be expressed by a unit angle intensity in an arbitrary range. FIG. 10 is a mathematical formula for obtaining the sound intensity (level in this case) of the intensity in a predetermined range (azimuth angles θ _{1 to θ 2} , elevation angles φ _{1 to φ 2).} The sound intensity of an arbitrary range of intensity is obtained using the formula shown in this figure, and the unit angle strength is obtained by dividing by the solid angle corresponding to the range, and this unit angle strength is used in the Cartesian coordinate system. The intensity distribution map is created.

(Polar coordinate system intensity distribution map)
FIG. 11A is a model diagram showing an example of a polar coordinate system intensity distribution map. As shown in this figure, the polar coordinate system intensity distribution map is a distribution map showing the distribution of unit angle intensities on the plane of the polar coordinate system consisting of a radius and a central angle. Hereinafter, the procedure for creating the polar coordinate system intensity distribution map will be described in detail with reference to FIG.

The unit angle strength is calculated by causing the computer to execute the unit angle strength calculation process 104, that is, when the unit angle strength is obtained for each division region MS, the arrival angle is calculated by causing the computer to execute the strength calculation process for each arrival angle. "Another strength" is calculated. The intensity for each arrival angle is, for example, a divided region MS within a range surrounded by a vertical line (that is, a longitude line) shown in FIG. 7, in other words, a predetermined azimuth interval Δθ (hereinafter, referred to as “unit longitude width”). It is a value obtained by synthesizing (integrated value or total value) the unit angle intensity of a certain divided region MS. At this time, a desired elevation angle interval Δφ (hereinafter referred to as “latitude range”) surrounded by the horizontal line (that is, the latitude line) shown in FIG. 7 is specified, and the split area MS having a unit longitude width and within the latitude range is specified. It is also possible to obtain the strength for each arrival angle by synthesizing the unit angle strength.

A polar coordinate system intensity distribution map is created by distributing the intensity for each arrival angle calculated by causing a computer to execute the intensity calculation process for each arrival angle on the radar chart shown in FIG. 11B. Specifically, the intensity for each arrival angle is arranged on the radar chart according to the azimuth information (direction and unit longitude width) of the intensity for each arrival angle, and the region extends in the radial direction according to the intensity (large or small value). Is set and a color or the like is applied to the area to create a polar coordinate system intensity distribution map. In this case as well, the unit angle intensity can be expressed by a color model, a shade, or a combination thereof, as in the spherical intensity distribution map and the Cartesian coordinate system intensity distribution map.

The unit longitude width (or the latitude range orthogonal to the unit longitude width) when determining the strength for each arrival angle can be appropriately changed by changing the direction of the latitude and longitude set in the reference spherical surface BF. For example, when the XYZ coordinate systems shown in FIGS. 6 and 8 are fixed, the strength for each arrival angle can be obtained after the ZZ plane is the equator plane and the X-axis direction is the latitude direction. , The arrival angle can be obtained by setting the ZZ plane as the equatorial plane and the Y-axis direction as the latitude direction, or setting the equator plane inclined with respect to the XYZ and then the arrival angle. It is also possible to obtain another strength. In FIG. 12, when the XY plane is the equatorial plane, the XY plane is the equatorial plane, and the YY plane is the equatorial plane, the sound intensity according to the arrival angle is determined in each case. The formula to be obtained is shown, and the strength of each arrival angle can be obtained by dividing the sound intensity by the solid angle.

5. Sound source type estimation process As described above, the sound source type estimation process 109 estimates the sound source type (noise caused by construction work, sound caused by driving a car, sound of home appliances, etc.) based on the instantaneous intensity. FIG. 13 is a flow chart showing a flow of a process of estimating a sound source type by causing a computer to execute a sound source type estimation process 109. When the instantaneous intensity is calculated as shown in this figure, in the sound source type estimation process 109, the sound characteristics of this instantaneous intensity (arrival direction, sound intensity, time change of arrival direction and sound intensity). And the sound characteristics of known equipment are compared. More specifically, the instantaneous calculated by having the computer execute the sound source type estimation process 109 shown in FIG. 1, reading the sound characteristics of the known equipment from the sound characteristic storage means 110, and causing the computer to execute the intensity calculation process 103. Compare with the sound characteristics of the intensity of. At this time, it is possible to compare the sound characteristics with the sound characteristics of known equipment based on the direction of arrival of the instantaneous intensity, the sound intensity, and the time change of the sound intensity, or the sound characteristics can be compared with the sound characteristics of the known equipment. It can be compared with the sound characteristics of known equipment based on the direction of arrival of intensity, sound intensity, and time change in the direction of arrival, and among the sound characteristics, the direction of arrival of intensity, sound intensity, etc. It is also possible to compare the sound characteristics of known equipment with the time change of the arrival direction and the time change of the sound intensity. Then, if the sound characteristics of the instantaneous intensity and the sound characteristics of the known equipment are similar (or the same), the instantaneous intensity is recognized as the sound source of the known equipment, and the sound source type is selected. If it is estimated and the sound characteristics of the instantaneous intensity and the sound characteristics of the known equipment are different, the instantaneous intensity is recognized as an unknown sound source.

For example, the ceiling ventilation installed in the office has a known direction (azimuth and elevation) with respect to the sound receiving point SP, and the strength and frequency of the sound can also be known, and the sound intensity can be further known. The temporal change (so to speak, the periodic change) can also be known. Then, if the instantaneous intensity of the sound characteristic equal to or similar to the sound characteristic of the ceiling ventilation is detected, the instantaneous intensity is recognized as the sound source of the ceiling ventilation. Of course, not limited to ceiling ventilation, various sound source types can be estimated by storing the sound characteristics of various known facilities in the sound characteristic storage means 110 and causing the computer to execute the sound source type estimation process 109.

When the sound source type is estimated by causing the computer to execute the sound source type estimation process 109, it is possible to create an intensity distribution map of only the instantaneous intensity at which the sound source type is estimated, or to make the instantaneous intensity. It is also possible to create an intensity distribution map in which the divided region MS is represented by a color or shade different from the others.

6. Specific sound detection process As described above, the specific sound detection process 107 detects a specific sound based on the instantaneous intensity. In this case, the focus range and the intensity threshold value for detecting the specific sound are stored by the focus range storage means 108 (FIG. 1). The range of interest is a range including a "set arrival direction range" set in advance for detecting a specific sound, and may be a range including a "set frequency range".

FIG. 14 is a flow chart showing a flow of processing for detecting a specific sound by causing a computer to execute the specific sound detection process 107. When the instantaneous intensity is calculated as shown in this figure, it is determined whether or not the instantaneous intensity is within the set arrival direction range. Specifically, by causing the computer to execute the specific sound detection process 107, the focus range is read from the focus range storage means 108, and the set arrival direction range of the focus range and the arrival direction (or arrival point) of the instantaneous intensity are obtained. Compare AP). Then, when the arrival direction is within the set arrival direction range (Yes), the process proceeds to the next process, and when the arrival direction is not within the set arrival direction range (No), it is determined that the sound is not a specific sound. However, when the set frequency range is not included in the focus range and the arrival direction is within the set arrival direction range (Yes), the process proceeds to the step of comparing with the intensity threshold value as it is.

When it is determined that the arrival direction is within the set arrival direction range, the next step is to compare the set frequency range in the range of interest with the frequency of the instantaneous intensity. Then, when the instantaneous intensity frequency is within the set frequency range (Yes), the process proceeds to the next process, and when the instantaneous intensity frequency is outside the set frequency range (No), it is determined that the sound is not a specific sound.

When it is determined that the frequency of the instantaneous intensity is within the set frequency range, the intensity threshold is then compared with the sound intensity (or unit angle intensity) of the instantaneous intensity. When the intensity of the instantaneous intensity sound exceeds the intensity threshold (Yes), it is recognized as a specific sound, and when the intensity of the instantaneous intensity sound is below the intensity threshold (No), it is recognized as not a specific sound. do.

7. Intensity extraction process If an intensity distribution map is created at all sound intensity collection times of the sound intensity data, the calculation cost to create it will be high and the time and effort to check it will increase. As a result, it may take a lot of time to discover the problem when trying to take noise countermeasures. Therefore, it is preferable to create an intensity distribution map after extracting only the sounds of interest in advance.

In the intensity extraction process 113, the instantaneous intensity of interest is extracted from the two or more instantaneous intensities obtained from the sound collection data collected for a certain period of time. Specifically, a specific arrival direction range (hereinafter referred to as “extraction arrival direction range”) and intensity threshold value (hereinafter referred to as “extraction intensity threshold”) are determined in advance, and the instantaneous arrival of intensity is achieved. When the direction is within the range of the arrival direction for extraction and the intensity of the sound exceeds the intensity threshold for extraction, the instantaneous intensity is extracted. Then, the sound intensity pick-up time of the instantaneous intensity extracted by causing the computer to execute the intensity extraction process 113 is selected, and the intensity is determined by the instantaneous intensity of the sound pick-up data related to the sound pick-up time selected here. Create a distribution map. It is also possible to create an intensity distribution map for all instantaneous intensities related to the selected sound collection time (including those that are not extracted by causing the computer to execute the intensity extraction process 113). It is also possible to create an intensity distribution map for the instantaneous intensity related to the selected sound collection time extracted by the intensity extraction process 113.

8. Other Examples Up to this point, an example of creating a “intensity distribution map” using the “unit angle intensity” obtained by dividing the sound intensity by the solid angle of the divided region MS has been described. The display program 100 sets the "unit angle intensity" as it is without dividing the sound intensity by the stereoscopic angle, and then assigns this "unit angle intensity" to the corresponding divided region MS to obtain the "intensity distribution map". Can also be created.

The sound collection data display program of the present invention can be effectively used not only for facilities on the noise receiving side such as houses, schools and hospitals, but also for facilities on the noise receiving side such as event venues and construction sites. Considering that the invention of the present application contributes to effective noise control, provides a good sound environment, and eventually provides a comfortable living environment, it can be expected to contribute not only industrially but also socially. It can be said to be an invention.

100 Sound collection data display program 101 (Sound collection data display program) 3D microphone 102 (Sound collection data display program) Sound collection data storage means 103 (Sound collection data display program) Intensity calculation processing 104 (Sound collection data display) Unit angle intensity calculation process (of program) 105 (of sound collection data display program) Frequency condition setting means 106 (of sound collection data display program) Unit angle intensity storage means 107 (of sound collection data display program) Specific sound detection process 108 ( Range of interest storage means (of sound collection data display program) 109 Sound type estimation processing (of sound collection data display program) 110 Sound characteristic storage means (of sound collection data display program) 111 Intensity distribution map (of sound collection data display program) Creation processing 112 (Sound collection data display program) Time-series display processing 113 (Sound collection data display program) Intensity extraction processing AP arrival point BF Reference spherical SP Sound reception point MS Division area IT Instantaneous intensity

Claims (12)

It is a sound collection data display program that causes a computer to execute a function of analyzing sound collection data collected for a predetermined period by a 3D microphone and displaying the analysis result using a reference spherical surface.
The reference spherical surface is a spherical surface set with a predetermined radius around the sound receiving point.
Intensity calculation processing that calculates the instantaneous intensity for each sound collection time from the sound collection data, and
The instantaneous intensity is assigned to the divided region obtained by dividing the reference spherical surface into a plurality of parts, and the "unit angle intensity" obtained by dividing the sound intensity of the instantaneous intensity by the solid angle of the divided region is obtained. Unit angle strength calculation processing and
It has a function of causing the computer to execute an intensity distribution map creation process for creating an "intensity distribution map" in which the unit angle intensity is distributed on the reference sphere.
A sound collection data display program characterized by this. It is a sound collection data display program that causes a computer to execute a function of analyzing sound collection data collected for a predetermined period by a 3D microphone and displaying the analysis result using a reference spherical surface.
The reference spherical surface is a spherical surface set with a predetermined radius around the sound receiving point.
Intensity calculation processing that calculates the instantaneous intensity for each sound collection time from the sound collection data, and
The instantaneous intensity is assigned to the divided region obtained by dividing the reference sphere into a plurality of parts, and the "unit angle intensity" obtained by dividing the sound intensity of the instantaneous intensity by the solid angle of the divided region is obtained. Unit angle strength calculation processing and
An intensity distribution map that transforms the divided region into a Cartesian coordinate system consisting of latitude and longitude set on the reference sphere and creates an "intensity distribution map" in which the unit angle intensity is distributed in the Cartesian coordinate system. It has a function to make the computer execute the creation process.
A sound collection data display program characterized by this. It is a sound collection data display program that causes a computer to execute a function of analyzing sound collection data collected for a predetermined period by a 3D microphone and displaying the analysis result using a reference spherical surface.
The reference spherical surface is a spherical surface set with a predetermined radius around the sound receiving point.
Intensity calculation processing that calculates the instantaneous intensity for each sound collection time from the sound collection data, and
The instantaneous intensity is assigned to the divided region obtained by dividing the reference spherical surface into a plurality of parts, and the "unit angle intensity" obtained by dividing the sound intensity of the instantaneous intensity by the solid angle of the divided region is obtained. Unit angle strength calculation processing and
In addition to setting the latitude and longitude on the reference sphere, the strength calculation process for each arrival angle, which calculates the "strength for each arrival angle", which is the unit angle strength for each unit longitude width,
It is provided with a function of causing the computer to execute an intensity distribution map creation process for creating an "intensity distribution map" in which the intensities for each arrival angle are distributed for each direction on a reference cross section for cutting the reference spherical surface.
A sound collection data display program characterized by this. In the intensity calculation process for each arrival angle, the intensity for each arrival angle is calculated based on the instantaneous intensity within a designated latitude range.
3. The sound collection data display program according to claim 3. The computer is made to execute a sound source type estimation process for estimating the sound source type of the intensity based on the intensity, the time change of the sound intensity of the intensity, and the sound source type information stored in advance. With more features,
The sound source type information includes the sound source type, the direction of arrival of the sound, the strength of the sound, and the time change of the sound.
In the intensity distribution map creation process, an "intensity distribution map" in which the unit angle intensity is distributed based on the instantaneous intensity at which the sound source type is estimated by the sound source type estimation process is created.
The sound collection data display program according to any one of claims 1 to 4, wherein the sound collection data display program is characterized. A function of causing the computer to execute a sound source type estimation process for estimating a sound source type of the intensity based on the intensity, a time change in the direction of arrival of the intensity, and sound source type information stored in advance. Further prepare
The sound source type information includes a sound source type, a sound arrival direction, a sound intensity, and a time change in the arrival direction.
In the intensity distribution map creation process, an "intensity distribution map" in which the unit angle intensity is distributed based on the instantaneous intensity at which the sound source type is estimated by the sound source type estimation process is created.
The sound collection data display program according to any one of claims 1 to 5, wherein the sound collection data display program is characterized. Further, a function of causing the computer to execute a specific sound detection process for detecting the instantaneous intensity in a specific arrival direction range and the sound intensity exceeds a predetermined intensity threshold value. Prepare,
In the intensity distribution map creation process, an "intensity distribution map" in which the unit angle intensity based on the instantaneous intensity detected by the specific sound detection process is distributed is created.
The sound collection data display program according to any one of claims 1 to 6, wherein the sound collection data display program is characterized. In the specific sound detection process, the instantaneous intensity in a specific frequency range is further detected.
7. The sound collection data display program according to claim 7. Of the two or more instantaneous intensities obtained from the sound collection data collected for a predetermined period, the instantaneous intensity that is in a specific arrival direction range and the sound intensity exceeds a predetermined intensity threshold is extracted. Further provided with a function of causing the computer to execute the intensity extraction process to be performed.
In the intensity distribution map creation process, an intensity distribution map at the sound intensity pick-up time of the instantaneous intensity extracted by the intensity extraction process is created.
The sound collection data display program according to any one of claims 1 to 8, wherein the sound collection data display program is characterized in that. In the unit angle intensity calculation process, the "unit time intensity", which is the sound intensity of the instantaneous intensity in a predetermined time width, is obtained, and the unit angle intensity is obtained based on the unit time intensity.
The sound collection data display program according to any one of claims 1 to 9, wherein the sound collection data display program is characterized. In the intensity calculation process, the instantaneous intensity is calculated based on the sound collection data extracted according to the preset frequency conditions.
The sound collection data display program according to any one of claims 1 to 10. 1. The sound collection data display program described in any of.

Images