WO2011027452A1 - 存在人数推定装置 - Google Patents
存在人数推定装置 Download PDFInfo
- Publication number
- WO2011027452A1 WO2011027452A1 PCT/JP2009/065432 JP2009065432W WO2011027452A1 WO 2011027452 A1 WO2011027452 A1 WO 2011027452A1 JP 2009065432 W JP2009065432 W JP 2009065432W WO 2011027452 A1 WO2011027452 A1 WO 2011027452A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- room
- value
- unit
- information
- absence
- Prior art date
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04Q—SELECTING
- H04Q9/00—Arrangements in telecontrol or telemetry systems for selectively calling a substation from a main station, in which substation desired apparatus is selected for applying a control signal thereto or for obtaining measured values therefrom
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F24—HEATING; RANGES; VENTILATING
- F24F—AIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
- F24F11/00—Control or safety arrangements
- F24F11/30—Control or safety arrangements for purposes related to the operation of the system, e.g. for safety or monitoring
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F24—HEATING; RANGES; VENTILATING
- F24F—AIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
- F24F11/00—Control or safety arrangements
- F24F11/62—Control or safety arrangements characterised by the type of control or by internal processing, e.g. using fuzzy logic, adaptive control or estimation of values
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F24—HEATING; RANGES; VENTILATING
- F24F—AIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
- F24F11/00—Control or safety arrangements
- F24F11/62—Control or safety arrangements characterised by the type of control or by internal processing, e.g. using fuzzy logic, adaptive control or estimation of values
- F24F11/63—Electronic processing
- F24F11/64—Electronic processing using pre-stored data
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F24—HEATING; RANGES; VENTILATING
- F24F—AIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
- F24F2120/00—Control inputs relating to users or occupants
- F24F2120/10—Occupancy
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04Q—SELECTING
- H04Q2209/00—Arrangements in telecontrol or telemetry systems
- H04Q2209/40—Arrangements in telecontrol or telemetry systems using a wireless architecture
- H04Q2209/43—Arrangements in telecontrol or telemetry systems using a wireless architecture using wireless personal area networks [WPAN], e.g. 802.15, 802.15.1, 802.15.4, Bluetooth or ZigBee
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04Q—SELECTING
- H04Q2209/00—Arrangements in telecontrol or telemetry systems
- H04Q2209/80—Arrangements in the sub-station, i.e. sensing device
- H04Q2209/82—Arrangements in the sub-station, i.e. sensing device where the sensing device takes the initiative of sending data
- H04Q2209/823—Arrangements in the sub-station, i.e. sensing device where the sensing device takes the initiative of sending data where the data is sent when the measured values exceed a threshold, e.g. sending an alarm
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04Q—SELECTING
- H04Q2209/00—Arrangements in telecontrol or telemetry systems
- H04Q2209/80—Arrangements in the sub-station, i.e. sensing device
- H04Q2209/84—Measuring functions
- H04Q2209/845—Measuring functions where the measuring is synchronized between sensing devices
Definitions
- the present invention relates to an existence number estimating device for estimating the number of persons existing in each area such as a residence and an office.
- Patent Document 1 discloses a monitoring system that uses a combination of a pyroelectric sensor and a number detection sensor to estimate the occupancy status of a resident in each area.
- the existence number estimating device is required to accurately estimate the number of existence persons.
- Highly accurate estimation of the number of people in each area enables device control according to the number of people in each area, for example, lighting on / off, operation mode of air conditioning equipment, etc., reducing unnecessary energy consumption can do.
- the present invention has been made to solve the above problems, and an object of the present invention is to provide a presence person estimation device capable of accurately estimating the number of people for each area at a low cost.
- a plurality of devices are installed in any one of a plurality of regions within a predetermined range to be estimated, and the presence or absence of humans in the region or the presence or absence of human movements in the region is detected.
- a human body detection sensor that outputs a detection signal; and a human body that collects the detection signal and includes time information in which the time at which the presence / absence of the human being or the presence / absence of the human motion was last detected in the region is associated with each region.
- a collection unit that outputs detection information, a plurality of realized values indicating the estimated number of persons for each region, a connection relationship between the regions, and a connection relationship between the region and the outside of the region from which the number of people is not estimated.
- a storage unit that stores connection relation information including the overlapping relation information indicating the overlapping area where the areas overlap, and a first transition that indicates a probability that a person will transition between the areas using the connection relation information. From probability and region Calculating a second transition probability indicating the probability that a person will transition out of the region, creating a transition matrix having the first and second transition probabilities as matrix elements, and using the transition matrix
- An actual value predicting unit that predicts an actual value at the current time from each of the actual values, and a likelihood calculation unit that calculates the likelihood of each actual value at the current time using the time information and the overlapping relationship information
- an actual value selection unit that selects an actual value having a likelihood greater than a preset threshold among the actual values at the current time, and calculates the number of identical actual values from the selected actual value.
- an output unit that generates output information in which the number of realization values and the number of people in each region are associated with each other, and outputs the output information.
- the existence number estimating device of the present invention the number of existence persons for each area can be estimated with high accuracy.
- FIG. 1 shows a schematic configuration of an existence number estimating device 100 according to the first embodiment of the present invention.
- This existence number estimating device 100 can estimate the number of persons existing indoors such as a residence or office as will be described later with reference to FIG. 3 for each room (area).
- the room described in this specification is not limited to a room such as a living room, a bedroom, and a kitchen in a house, and may mean an arbitrary area (area) for which the number of people is estimated.
- a human body detection sensor group 101 including a plurality of human body detection sensors 101A to 101C, and each of the human body detection sensors 101A to 101C is assigned a specific area as a detection area. .
- Each of the human body detection sensors 101A to 101C detects the presence or absence of a human in the detection area or the presence or absence of a human motion, and transmits a detection signal to the detection information collecting unit 102.
- the detection information collection unit 102 is connected to the human body detection sensors 101A to 101C through a wired or wireless route.
- Examples of the human body detection sensors 101A to 101C include sensors that detect the presence of a moving human such as a pyroelectric infrared sensor and a thermopile infrared sensor.
- a moving human such as a pyroelectric infrared sensor and a thermopile infrared sensor.
- the human body detection sensors 101A to 101C are arranged on the ceiling or the like, and detect the presence or absence of a human by detecting the movement of a human in the detection area. More specifically, the human body detection sensors 101A to 101C including the pyroelectric infrared sensor detect a change in infrared energy emitted from the detection region as an electrical change by a pyroelectric element, so that a human who emits infrared rays.
- the human body detection sensors 101A to 101C are not limited to the above example, and may be any camera as long as they can detect the presence or absence of humans or the presence or absence of human movements. It may be configured with household electrical appliances that can be output as The human body detection sensor group 101 may be configured by combining various sensors as described above.
- a detection information collection unit (also simply referred to as a collection unit) 102 collects detection signals from the human body detection sensors 101A to 101C and transmits the collected detection signals to the control unit 103 and the actual value evaluation unit 105 as human body detection information. .
- the detection information collecting unit 102 can transmit the human body detection information including time information related to the time when the detection signal is received, that is, the time when the presence or absence of a human being or the presence or absence of a human motion is detected. More specifically, the detection information collection unit 102 stores the time (final detection time) at which the presence or absence of a person or the presence or absence of human movement was last detected in each room, and includes a human body including the final detection time for each room. Send detection information periodically. Alternatively, the detection information collection unit 102 may calculate the undetected time indicating the elapsed time from the last detection time, and transmit the human body detection information including the undetected time for each room.
- the control unit 103 sequentially operates the one-year ahead prediction unit 104, the actual value evaluation unit 105, and the occupancy number output unit 106 at a predetermined timing, for example, immediately after receiving the human body detection information from the detection information collection unit 102.
- the control unit 103 can initialize a set 111 of realized values storing a plurality of realized values.
- the actual value refers to data having an estimated value of the number of people in each area (for example, a room) for which the number of people is to be estimated.
- These realization values have a data structure as shown in FIG. 11, as will be described later.
- the first-term prediction unit 104 includes a transition matrix creation unit 104A and an actual value prediction unit 104B, and the transition matrix creation unit 104A includes connection relation information 109 between rooms and connection relation information 110 with the outdoors.
- the realization value prediction unit 104B is initialized by the control unit 103 using the transition matrix, or updated by the one-time ahead prediction unit 104 and the realization value evaluation unit 105 before using the transition matrix.
- the actual value at the current time is predicted from each realized value (also referred to as the actual value one hour before), and the set 111 of actual values is updated with the predicted actual value.
- Realized value evaluation unit 105 uses the human body detection information received from detection information collection unit 102 to evaluate the actual value predicted by first-term ahead prediction unit 104. More specifically, the actual value evaluation unit 105 deletes the actual value that deviates from the occupancy status of each room based on the human body detection information, and duplicates the actual value that matches this occupancy status to The set 111 is updated.
- the realization value evaluation unit 105 calculates the likelihood for each realization value by using the human body detection information and the overlapping relationship information 112 regarding the area where the detection areas of the human body detection sensors 101A to 101C overlap (also referred to as an overlap area).
- a likelihood calculation unit 105A is provided, and a realization value selection unit 105B that compares the realization value of each realization value and deletes the realization value and a realization value to be duplicated.
- the occupant number output unit (also referred to as an output unit) 106 generates output information related to the occupancy number for each room from the actual values in the actual value set 111 updated by the evaluation unit 105, and displays the output information on the display device 107. And output to the control target device 108.
- the display device 107 includes a display such as a display that can display the output information of the occupant number output unit 106, a printer, and the like, and can include the display unit 204 shown in FIG.
- the control target device 108 is configured by a general-purpose computer, a home appliance, a residential device, or the like that is controlled according to the output information from the occupancy number output unit 106.
- the occupancy output unit 106 is connected to the display device 107 and the control target device 108 via a communication network such as a wired LAN (Local Area Network), a wireless LAN, Bluetooth (registered trademark), serial communication, or parallel communication. .
- FIG. 2 shows a hardware configuration example of the presence number estimating device 100.
- the existence number estimating device 100 can be realized by using the general-purpose computer device 200 as basic hardware.
- the general-purpose computer device 200 includes a CPU 202, an input unit 203, a display unit 204, a communication unit 205, a main storage unit 206, and an external storage unit 207, which communicate with each other via a bus 201. Connected as possible.
- the input unit 203 includes input devices such as a keyboard and a mouse, and outputs an operation signal generated by operating the input device to the CPU 202.
- the display unit 204 includes a display such as an LCD (Liquid Crystal Display) and a CRT (Cathode Ray Tube).
- a display such as an LCD (Liquid Crystal Display) and a CRT (Cathode Ray Tube).
- the communication unit 205 includes communication means such as Ethernet (registered trademark), wireless LAN (Local Area Network), and Bluetooth (registered trademark), and communicates with the human body detection sensor 210.
- communication means such as Ethernet (registered trademark), wireless LAN (Local Area Network), and Bluetooth (registered trademark), and communicates with the human body detection sensor 210.
- the external storage unit 207 includes, for example, a hard disk or a storage medium such as a CD-R, CD-RW, DVD-RAM, DVD-R, and the like, and includes the above-described detection information collection unit 102, control unit 103, and one-time destination.
- a control program for causing the CPU 202 to execute processing by the prediction unit 104, the actual value evaluation unit 105, and the occupant number output unit 106 is stored.
- the main storage unit 206 is configured by a memory or the like, expands a control program stored in the external storage unit 207 under the control of the CPU 202, data necessary for execution of the program, and data generated by execution of the program Memorize etc.
- the existence estimation device 100 may be realized by installing the control program in a computer device in advance, or may be stored in a storage medium such as a CD-ROM or distributed through the network. Then, this program may be realized by appropriately installing it in a computer device. Also, the inter-room connection relationship information 109, the outdoor connection relationship information 110, the realization value set 111, and the overlapping relationship information 112 shown in FIG. It can be realized by appropriately using a memory, a hard disk, or a storage medium such as a CD-R, CD-RW, DVD-RAM, or DVD-R, such as the external storage unit 207.
- connection relation information 109 between the rooms in order to print the connection relation information 109 between the rooms, the outdoor connection relation information 110, the realization value set 111, the information stored in the duplication relation information 112, an abnormality notification, and the like.
- the printer may be provided. Further, the hardware configuration of the existence number estimating device 100 shown in FIG. 2 may be changed according to the target device for collecting the usage status.
- FIG. 3 shows a floor plan of a house as an example to which the existence number estimating device 100 shown in FIG. 1 is applied.
- the human body detection sensor 101 ⁇ / b> A is disposed in an entrance (also referred to as a room A) 301, and an area 311 ⁇ / b> A in the entrance 301 is assigned as a detection area.
- the human body detection sensors 101B to 101E include a kitchen (also referred to as a room B) 302, a north side of the western room 303 (also referred to as a room C), a south side of the western room 303 (also referred to as a room D), and a washroom (also referred to as a room E).
- a human body detection sensor is not installed in the bathroom (also referred to as room F) 305.
- the arrows shown in FIG. 3 indicate the connection relationship between movable rooms.
- a person in the room A can move to either the room B or the outdoor (room X, also referred to as outside the area) that is not subject to estimation of the number of people.
- Examples of the connection relation information 109 between rooms and the connection relation information 110 with the outdoors in the floor plan of FIG. 3 are shown in FIGS. 4 and 5, respectively.
- Each field of the connection relationship information 109 between rooms shown in FIG. 4 can be moved between rooms, that is, when the rooms are adjacent to each other, 1 is stored, and when the rooms are not adjacent to each other, 0 is stored.
- the connection relationship information 109 between rooms shown in FIG. 4 indicates that the room A can move to the room B and can stay in the room A.
- each field of the outdoor connection relationship information 110 shown in FIG. 5 stores 1 when the room and the outdoor are adjacent to each other, and 0 when the room and the outdoor are not adjacent. Stored.
- the connection relation information 110 with the outdoors shown in FIG. 5 indicates that it is possible to move from the room A to the outdoors (room X) and from the outdoors to the room A.
- connection relationship between the rooms and the connection relationship between the room and the outdoors may be prepared as separate information, or may be prepared together in the same connection relationship information.
- FIG. 6 shows an example of the overlapping relationship information 112 regarding the overlapping region of the human body detection sensor in the floor plan shown in FIG.
- the diagonal element that is, the element having the same row and column index
- the human body detection sensor is installed in the room specified by the index.
- the diagonal element is 0, it indicates that the human body detection sensor is not installed in the room specified by the index.
- the off-diagonal element is 1, it indicates that the detection areas of the two rooms specified by the index overlap each other.
- the human body detection sensors are installed in the rooms A to E, the human body detection sensors are not installed in the rooms F and G, and the detection areas of the rooms C and D overlap.
- the case where the overlapping area is a virtual room G is shown.
- FIG. 7 shows an example of the processing procedure of the existence number estimating device 100.
- the control unit 103 initializes a set 111 of realized values in step S701.
- the realization value set 111 is initialized by the following method. That is, the control unit 103 generates a plurality of realization values with random numbers so as to satisfy preset constraint conditions such as the total number of people in the residence and the maximum number of people in each room, or realizes all states satisfying the constraint conditions. Generate as a value.
- step S702 the control unit 103 temporarily stops processing until it receives detection state information from the detection information collection unit 102 or is periodically called by a built-in timer. After returning from the temporary stop, in step S703, the one-year ahead prediction process shown in FIG. 8 is executed, in step S704, the actual value evaluation process shown in FIG. 14 is executed, and in step S705, the occupancy number output process is executed. Is executed. The series of processing shown in steps S702 to S705 is repeatedly executed until the control unit 103 instructs to end (step S706), and the set of realized values 111 is updated one after another.
- FIG. 8 shows an example of the procedure of the one-year ahead prediction process shown in step S703 of FIG.
- the transition matrix creation unit 104A creates a transition matrix from the connection relation information 109 between rooms and the connection relation information 110 with the outdoors.
- the actual value prediction unit 104B uses the created transition matrix to predict the actual value at the current time from each actual value one time before, and updates the set 111 of actual values with the predicted actual value. To do.
- FIG. 9 shows an example of the procedure of the transition matrix creation process shown in step S801 of FIG.
- the transition matrix creation unit 104A determines the transition probability P ij when a person moves from one room (also referred to as the source room) i to another room (also referred to as the destination room) j. calculate.
- the transition probability P ij is a probability of moving to the adjacent room, a being a probability of moving to the outside, b, a connection relation between the room i and the room j, d ij , and a connection relation between the outdoor and the room i (room j), g.
- the movement probability a represents that the movement to another room is likely to occur as the value is large
- the movement probability b represents that the movement to the outdoors and the movement from the outdoors are likely to occur as the value is large.
- Transition probabilities P ij is calculated for all the rooms j other than the source room i (step S902).
- step S ⁇ b > 903 the first-term prediction unit 104 calculates a transition probability P ii when staying in the source room i.
- the transition probability P ii is calculated from the following formula 2 using the transition probability P ij calculated in steps S901 and S902.
- steps S901 to S903 are executed for all the movement source rooms i (step S904), and a transition matrix P having the transition probabilities calculated in steps S901 to S904 as matrix elements is created.
- FIG. 11 shows an example of the data structure of realized values in the set 111 of realized values.
- the realization value set 111 stores four realization values with realization values 1 to 4 and indexes.
- the number N of actual values stored in the actual value set 111 can be appropriately changed according to the situation such as the number of rooms.
- FIG. 12 shows the procedure of the actual value prediction process shown in step S802 of FIG.
- the realization value prediction unit 104B uses the transition matrix P calculated in step S801 to realize a realization value k one time ago (where k satisfies 1 ⁇ k ⁇ N). It is determined to which room j a person in room i in (integer) moves to the current time or to move outdoors. More specifically, since the horizontal vector in the i-th row of the transition matrix has a probability of moving from room i to each room, random numbers are generated according to this probability to sample one destination room j. For example, for a person in room G, as shown in FIG. 10, the probability of staying in room G is 0.8, and the probability of moving to room C and room D is 0.1. Whether to stay in room G or move to room C or room D is determined by a random number according to this probability.
- step S1202 the actual value prediction unit 104B decreases the number of persons in the source room i by 1 according to the destination room j determined in step S1201 (the room j may be the same as the room i). , Increase the number of people in the destination room j by one.
- the destination room determined in step S1201 is the room X
- the number of people in the source room i is decreased by one.
- the processing in steps S1201 and S1202 is repeated as many times as the number of people in the room i of the movement source at the realization value k one time before (step S1203).
- step S1204 the processing shown in steps S1201 to S1203 is repeatedly executed for all rooms, for example, room A to room F.
- step S1205 the processing shown in steps S1201 to S1204 is repeatedly executed for all realized values.
- FIG. 13 shows an example of the real value set 111 after the real value prediction process is performed on the real value set 111 shown in FIG.
- realization value 1 one person moves from room G to room D, and one person moves from room E to room F.
- realization value 2 one person moves from room E to room F, and in realization value 3, one person moves from room D to room G.
- realization value 4 one person moves from the room G to the room D, and one person moves from the room E to the room F.
- the one-period ahead prediction unit 103 predicts the actual value at the current time from the actual value one time before, and updates the set 111 of actual values.
- transition matrix creation process shown in step S801 of FIG. 8 is not necessarily executed every time, and the transition matrix P may be calculated in advance and stored in the external storage unit 207 or the like.
- FIG. 14 shows an example of the procedure of the actual value evaluation process shown in step S704 of FIG.
- the realization value evaluation unit 105 executes the likelihood calculation unit 105A and the realization value evaluation unit 105, and each realization value in the realization value set 111 updated by the one-prediction prediction unit 104 is based on human body detection information. Evaluate whether it is consistent with the occupancy status of each room.
- the likelihood calculation unit 105A executes a likelihood calculation process for calculating the likelihood w of each realization value in the realization value set 111 using the overlap relationship information 112 and the human body detection information. To do.
- the likelihood w of the realization value k is calculated by the following mathematical formula 3.
- a value indicating the presence / absence of a person in the room i calculated from the number of people in each room of the actual value k (also referred to as a first presence / absence value) s i and each specified by the human body detection information.
- H jj indicates the overlapping relationship information 112. Therefore, a virtual room (for example, room G) due to overlapping detection areas and a room (for example, room F) in which no human body detection sensor is installed are excluded.
- the room i becomes a calculation target of the likelihood w.
- the value s i is calculated by Equation 5 below.
- This equation 5 converts the number of people x i in the room i with the realization value k into presence (1) absence (0).
- the value r i is calculated by the following formula 6.
- j ⁇ adj (i) represents the adjacent room j of the room i
- y i represents the last detection time in the room i.
- the elapsed time t-y i from the last detection time y i to the current time t is consistent with the non-detection time in the room i.
- the function f is given as, for example, a function that monotonically decreases a value such as the following Equation 7 from 1 to 0 so as to reduce the probability that a person is present in the room i according to the undetected time t-yi. It is done.
- the parameter ⁇ i represents a non-detection time until no one is present in the room i
- the parameter ⁇ i is a room where no sensor is installed from the room i, outdoors, etc. (hereinafter referred to as a transition destination). It represents the undetected time until it is considered that the user has moved to the room.
- the parameters ⁇ i and ⁇ i can be set to different values for each room. For example, in the hallway (also referred to as the room S), it is unlikely that the human detection sensor will not react for a long time despite being in the room. There is a person inside but there is almost no movement, and there is a high possibility that there is no response for a long time despite being in the room. Therefore, a small value is set for ⁇ S and a large value is set for ⁇ T.
- FIG. 15 shows an example of a calculation result obtained by executing the likelihood calculation process in step S1401 on the realized value shown in FIG.
- the likelihood w of each realization value is normalized so that the sum of the likelihoods of all realization values is 1, and is represented as a ratio of likelihoods.
- step S1401 in FIG. 14 the likelihood w of the actual value k may be calculated by the following formula 8.
- x i new represents the number of people in the room i having the real value k calculated in the one-year ahead prediction process
- x i old is the room i having the real value k before being calculated in the one-year ahead prediction process.
- r is a uniform random number from 0 to 1.
- step S1402 of FIG. 14 the actual value evaluation unit 105 compares the likelihood w of each actual value calculated in step S1401, erases the actual value with a small likelihood w, and realizes the actual value with a large likelihood w.
- An actual value selection process for replicating is executed.
- the realized value selection process in step S1402 is illustrated in FIG. 16 as an example.
- step S1601 shown in FIG. 16 the actual value selection unit 105B calculates the cumulative likelihood that is the sum of the ratios of the likelihood w from the actual value 1 to the actual value k from the likelihood w of each actual value calculated in step S1401.
- the degree ratio v k is calculated as in Equation 9 below.
- N is the number of actual values stored in the actual value set 111
- w m represents the likelihood w of the actual value m.
- the actual value selection unit 105B duplicates and stores the actual value k by the number of integers n satisfying v k ⁇ 1 ⁇ n ⁇ ⁇ v k .
- the parameter ⁇ is a value exceeding 0 and less than 1, and is set in advance. Deletes the realization value k when there is no integer n satisfying v k ⁇ 1 ⁇ n ⁇ ⁇ v k .
- the realization degree selection unit 105B deletes the realization value with the low likelihood w and duplicates the realization value with the high likelihood w.
- a threshold value is set in advance for the likelihood ratio, and when the likelihood ratio calculated for each realization value in step S1402 is equal to or greater than the threshold value, the realization value is copied, and the likelihood ratio is less than the threshold value In such a case, the actual value may be deleted.
- FIG. 17 shows an example of the result of executing the actual value selection process on the actual value shown in FIG.
- the realization value 3 shown in FIG. 13 is deleted and the realization value 2 shown in FIG. 2 and realization value 3.
- step S 703 the one-time prediction process in step S703 shown in FIG. 7 and the actual value evaluation process in step S704 are sequentially executed, and the real value set 111 is updated.
- step S ⁇ b> 705 the occupancy number output unit 106 outputs output information related to the occupancy number for each room from the updated real value set 111 to the display device 107 and the control target device 108.
- the occupancy number output unit 106 calculates the number of realized values with the same number of occupants in all the rooms, and outputs the realized values rearranged in the order of the number as output information. be able to.
- the control target device 108 uses various realization values, or averages the top ten realization values. Can be executed. Further, as shown in FIG. 19, the occupancy number output unit 106 calculates an average value of the occupancy number for each room from the N actual values, and outputs the average occupancy number for each room as a representative value. May be. By outputting one data like the representative value, it becomes easy to replace the output information from the occupancy output unit 106 with the output information related to the occupancy for each room estimated by other methods. In addition, there is an advantage that the output information is easy to understand for humans.
- Examples of the control target device 108 controlled according to the number of people in the room include lighting, an air conditioner, a fan, an air cleaner, a television, and a personal computer.
- the lighting installed in the room 303 shown in FIG. 3 is controlled to be turned off when the total number of people in the rooms C, D, and G is estimated to be 0.2 or less.
- the air conditioner, the fan, the air cleaner, and the like can be controlled in strength or operation mode according to the number of people in the room. Thus, unnecessary power consumption can be suppressed by controlling the operation of various devices in accordance with the number of people in the room.
- the occupant number output unit 106 may output output information to a control device (not shown), and the control device may control the operation of the control target device 108 according to the number of occupants in the room.
- the room that is, the area
- the room for which the number of persons is estimated is not limited to the example shown in FIG. 3, and may include a living room, a toilet, a bedroom, a meeting room in an office, and the like.
- the existence number estimation device prepares a plurality of actual values indicating the number of persons for each room, predicts the actual value of the current time from these actual values, and predicts the current An actual value that matches the occupancy status based on the human body detection information is selected from the actual value of the time by likelihood calculation.
- This existence number estimation device can accurately estimate the number of persons for each room (area) by selecting an actual value that is highly consistent with the actual occupancy status from the predicted actual value.
- FIG. 20 to FIG. 22 the existence number estimating device according to the second embodiment will be described.
- An existence number estimating device 100 according to the second embodiment is shown in a simplified manner in FIG. 20, but has the same configuration as that of the first embodiment shown in FIG. In the second embodiment, the arrangement of the human body detection sensors is changed as compared with the first embodiment.
- FIG. 21 shows an example of detection areas assigned to the human body detection sensors 101A and 101B in the second embodiment.
- an area 2111B in the room 2102 is set as a detection area
- an area 2111A including an entrance (room A) between the room 2101 and the room 2102 is set as a detection area.
- the human body detection sensor 101 ⁇ / b> A is disposed near the entrance / exit between the room 2101 and the room 2102.
- FIG. 21 shows an example of detection areas assigned to the human body detection sensors 101A and 101B in the second embodiment.
- human body detection sensors 101A and 101B are provided by being incorporated in lighting installed on the ceiling of the room 2102, and the sensor 101A is set with the entrance of the room as a detection area, The sensor 101B may be set as a detection area in the room. As described above, one or more human body detection sensors may be provided integrally with another device such as an air conditioner or a lighting device.
- the detection area of the human body detection sensor is assigned to the entrance / exit, it is possible to reliably detect movement between human rooms.
- the number of people in each (area) can be estimated with higher accuracy.
- the existence number estimating device can accurately estimate the number of persons existing indoors, such as a house and an office, for each area, and device control of lighting, air conditioner, television, etc. can be performed according to the estimated number of existence. Since it becomes possible, the industrial applicability is high.
- DESCRIPTION OF SYMBOLS 100 Presence number estimation apparatus, 101 ... Human body detection sensor group, 101A-101E ... Human body detection sensor, 102 ... Detection state collection part, 103 ... Control part, 104 ... First term prediction part, 104A ... Transition matrix production
Landscapes
- Engineering & Computer Science (AREA)
- Combustion & Propulsion (AREA)
- Signal Processing (AREA)
- General Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Chemical & Material Sciences (AREA)
- Fuzzy Systems (AREA)
- Mathematical Physics (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Air Conditioning Control Device (AREA)
- Alarm Systems (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
存在人数推定装置(100)においては、検知領域内の人間の有無を検知して検知信号を出力する人体検知センサ(101A~101C)と、前記検知信号を収集し、最終検知時刻を部屋毎に関連づけた時刻情報を含む人体検知情報を出力する収集部(102)と、部屋毎の在室人数を示す複数の実現値、前記部屋間の接続関係(109)及び部屋と屋外との接続関係(110)を含む接続関係情報、並びに検知領域の重複領域を示す重複関係情報(112)を格納する記憶部と、接続関係情報を使用して、一時刻前の実現値から現在時刻の実現値を予測する一期先予測部(104)と、時刻情報及び重複関係情報(112)を使用して1以上の実現値を選択する実現値選択部(105)と、選択された1以上の実現値に基づいた出力情報を出力する出力部(106)と、を具備する。
Description
本発明は、住居、オフィス等の各領域に存在する人数を推定する存在人数推定装置に関する。
生活者の見守り、及び状況に合わせた機器制御を行う場合には、生活者の位置等の状態を監視する必要がある。位置情報の取得方法としては、従来、監視カメラ、赤外線イメージセンサ、フロア圧力センサ、超音波センサ、無線タグ及びリーダ等を利用して行う方法が知られている。しかしながら、これら高度なセンサを使用する場合には、プライバシ及びコストに関する問題がある。
特許文献1には、焦電センサと人数検知センサを組み合わせて使用し、領域毎の生活者の在室状況を推定する監視システムが開示されている。
特許文献1に開示される監視システムでは、出入り口の通過人数を検知して特定の領域内に存在する人数を推定するが、この出入り口においては同時に1人の通過しか考慮しておらず、正確な存在人数を領域毎に推定することが困難とされる。また、人数検知センサが別途必要なことから、コスト面の問題がある。
従って、存在人数推定装置においては、存在人数を精度良く推定することが求められている。高精度な領域毎の存在人数の推定により、領域毎の存在人数に応じた機器制御、例えば、照明のオン・オフ、空調設備の運転モード等の制御が可能となり、不必要なエネルギー消費を低減することができる。
本発明は、上記問題点を解決するためになされてものであり、その目的は、低コストで、人数を領域毎に精度良く推定することができる存在人数推定装置を提供することにある。
本発明によれば、人数を推定する対象となる所定の範囲の複数の領域のいずれかに複数個設置され、前記領域内の人間の有無又は前記領域内の人間の動作の有無を検知して検知信号を出力する人体検知センサと、前記検知信号を収集し、前記領域で前記人間の有無又は前記人間の動作の有無が最後に検知された時刻を前記領域毎に関連づけた時刻情報を含む人体検知情報を出力する収集部と、前記領域毎の推定存在人数を示す複数の実現値と、前記領域間の接続関係及び前記領域と前記人数を推定する対象外となる領域外との接続関係を含む接続関係情報と、前記領域が重なり合う重複領域を示す重複関係情報と、を格納する記憶部と、前記接続関係情報を使用して、前記領域間を人間が遷移する確率を示す第1の遷移確率、及び領域から領域外へ人間が遷移する確率を示す第2の遷移確率を計算し、前記第1及び第2の遷移確率を行列要素に有する遷移行列を作成する遷移行列作成部と、前記遷移行列を使用して前記実現値の各々から現在時刻の実現値を予測する実現値予測部と、前記時刻情報及び前記重複関係情報を使用して前記現在時刻の実現値の各々の尤度を計算する尤度計算部と、前記現在時刻の各実現値のうち前記尤度が予め設定される閾値よりも大きい実現値を選択する実現値選択部と、前記選択された実現値から同一の実現値の数を計算し、該実現値の数と前記領域毎の存在人数とを関連づけた出力情報を生成して、この出力情報を出力する出力部と、を具備することを特徴とする存在人数推定装置が提供される。
本発明の存在人数推定装置によれば、領域毎の存在人数を精度良く推定することができる。
以下、必要に応じて図面を参照しながら、本発明の一実施の形態に係る存在人数推定装置を説明する。
[第1の実施の形態]
図1は、本発明の第1の実施の形態に係る存在人数推定装置100の概略構成を示している。この存在人数推定装置100は、図3を参照して後に説明されるような住居、オフィス等の屋内に存在する人数を部屋(領域)毎に推定することができる。本明細書に記述される部屋は、例えば、住宅における居間、寝室、台所等の部屋を指すように限定されず、人数を推定する対象となる任意の領域(エリア)を意味してもよい。
図1は、本発明の第1の実施の形態に係る存在人数推定装置100の概略構成を示している。この存在人数推定装置100は、図3を参照して後に説明されるような住居、オフィス等の屋内に存在する人数を部屋(領域)毎に推定することができる。本明細書に記述される部屋は、例えば、住宅における居間、寝室、台所等の部屋を指すように限定されず、人数を推定する対象となる任意の領域(エリア)を意味してもよい。
図1に示される存在人数推定装置100は、複数の人体検知センサ101A~101Cを含む人体検知センサ群101を備え、人体検知センサ101A~101Cは、夫々特定の領域を検知領域として割り当てられている。各人体検知センサ101A~101Cは、検知領域内の人間の有無又は人間の動作の有無を検知し、検知信号を検知情報収集部102に送信する。検知情報収集部102は、有線又は無線経路で各人体検知センサ101A~101Cに接続されている。
人体検知センサ101A~101Cとしては、例えば、焦電型赤外線センサ、サーモパイル型赤外線センサ等の動いている人間の存在を検知するセンサがある。人体検知センサ101A~101Cに焦電型赤外線センサが利用される場合、各人体検知センサ101A~101Cは、天井等に配置され、検知領域における人間の動きを検出して人間の有無を検知する。より詳細には、焦電型赤外線センサを備える人体検知センサ101A~101Cは、検知領域から発せられる赤外線エネルギーの変化を焦電素子によって電気的な変化として検出することで、赤外線を発している人間の動きを検出する。人体検知センサ101A~101Cは、上記の例に限定されず、人間の有無又は人間の動作の有無を検知することができればよく、カメラ、フロア圧力センサ、超音波センサ、操作されたことを検知信号として出力することができる家電機器等で構成されてもよい。人体検知センサ群101は、上記のような各種センサを組み合わせて構成されてもよい。
検知情報収集部(単に、収集部とも言う)102は、各人体検知センサ101A~101Cから検知信号を収集し、収集した検知信号を人体検知情報として制御部103及び実現値評価部105に送信する。また、検知情報収集部102は、検知信号の受信時の時刻、即ち、人間の有無又は人間の動作の有無が検知された時刻に関する時刻情報を人体検知情報に含めて送信することができる。より具体的には、検知情報収集部102は、各部屋において人間の有無又は人間の動作の有無が最後に検知された時刻(最終検知時刻)を記憶し、部屋毎の最終検知時刻を含む人体検知情報を定期的に送信する。或いは、検知情報収集部102は、最終検知時刻からの経過時間を示す未検知時間を計算し、部屋毎の未検知時間を含む人体検知情報を送信してもよい。
制御部103は、所定のタイミングで、例えば、検知情報収集部102から人体検知情報を受信した直後に一期先予測部104、実現値評価部105及び在室人数出力部106を順次動作させる。また、制御部103は、複数の実現値を格納している実現値の集合111を初期化することができる。ここで、実現値とは、人数を推定する対象となる領域(例えば、部屋)毎の存在人数の推定値を有するデータを指す。これらの実現値は、後に説明されるように、図11に示されるようなデータ構造を有している。
図1に示される一期先予測部104は、部屋間の接続関係を示す部屋間の接続関係情報109、及び部屋と屋外との接続関係を示す屋外との接続関係情報110を使用して実現値の集合111を更新する。より具体的には、一期先予測部104は、遷移行列作成部104A及び実現値予測部104Bを備え、遷移行列作成部104Aは、部屋間の接続関係情報109及び屋外との接続関係情報110を使用して遷移行列を作成し、実現値予測部104Bは、遷移行列を使用して、制御部103により初期化された、或いは一期先予測部104及び実現値評価部105により以前に更新された各実現値(一時刻前の実現値とも言う)から現在時刻の実現値を予測し、予測した実現値で実現値の集合111を更新する。
実現値評価部105は、検知情報収集部102から受信した人体検知情報を使用して、一期先予測部104により予測された実現値を評価する。より具体的には、実現値評価部105は、人体検知情報に基づいた各部屋の在室状況と乖離する実現値を消去し、この在室状況と整合する実現値を複製して実現値の集合111を更新する。この実現値評価部105は、人体検知情報、及び人体検知センサ101A~101Cの検知領域が重なり合う領域(重複領域とも言う)に関する重複関係情報112を使用して、実現値毎に尤度を計算する尤度計算部105Aと、各実現値の尤度を比較して消去する実現値と複製する実現値を選択する実現値選択部105Bとを備えている。
在室人数出力部(出力部とも言う)106は、評価部105により更新された実現値の集合111内の実現値から部屋毎の在室人数に関する出力情報を生成し、出力情報を表示装置107及び制御対象機器108に出力する。表示装置107は、在室人数出力部106の出力情報を表示することができる表示ディスプレイ、プリンタ等の機器で構成され、図2に示される表示部204を含むことができる。制御対象機器108は、在室人数出力部106の出力情報に応じて制御される汎用計算機、家電機器、住設機器等から構成される。在室人数出力部106は、有線LAN(Local Area Network)、無線LAN、Bluetooth(登録商標)、シリアル通信、又はパラレル通信等の通信ネットワークを介して表示装置107及び制御対象機器108に接続される。
図2は、存在人数推定装置100のハードウェア構成例を示している。存在人数推定装置100は、汎用コンピュータ装置200を基本ハードウェアとして使用することで実現することができる。汎用コンピュータ装置200は、図2に示されるように、CPU202、入力部203、表示部204、通信部205、主記憶部206及び外部記憶部207により構成され、これらは、バス201により相互に通信可能に接続される。
入力部203は、キーボード、マウス等の入力デバイスを備え、入力デバイスの操作による操作信号をCPU202に出力する。
表示部204は、LCD(Liquid Crystal Display)、CRT(Cathode Ray Tube)等の表示ディスプレイからなる。
通信部205は、Ethernet(登録商標)、無線LAN(Local Area Network)、Bluetooth(登録商標)等の通信手段を有し、人体検知センサ210との間で通信を行う。
外部記憶部207は、例えば、ハードディスク、若しくはCD-R、CD-RW、DVD-RAM、DVD-R等の記憶媒体等から構成され、上記した検知情報収集部102、制御部103、一期先予測部104、実現値評価部105及び在室人数出力部106による処理をCPU202に実行させるための制御プログラムが記憶されている。
主記憶部206は、メモリ等により構成され、CPU202による制御の下で、外部記憶部207に記憶された制御プログラムを展開し、当該プログラムの実行時に必要なデータ、当該プログラムの実行により生じたデータ等を記憶する。
尚、存在人数推定装置100は、上記制御プログラムをコンピュータ装置にあらかじめインストールすることで実現してもよいし、CD-ROM等の記憶媒体に記憶して、或いはネットワークを介して上記のプログラムを配布して、このプログラムをコンピュータ装置に適宜インストールすることで実現してもよい。また、図1に示した部屋間の接続関係情報109、屋外との接続関係情報110、実現値の集合111、重複関係情報112は、上記のコンピュータ装置に内蔵又は外付けされた主記憶部206、外部記憶部207といった、メモリ、ハードディスク、若しくはCD-R、CD-RW、DVD-RAM、DVD-R等の記憶媒体等を適宜利用して実現することができる。
尚、上述の構成要素の他に、部屋間の接続関係情報109、屋外との接続関係情報110、実現値の集合111、重複関係情報112に記憶されている情報及び異常通知等を印刷するためのプリンタを備えるようにしてもよい。また、使用状況を収集する対象の機器に応じて、図2に示す存在人数推定装置100のハードウェア構成を変更してもよい。
図3は、図1に示される存在人数推定装置100が適用される一例としての住居の間取りを示している。人体検知センサ101Aは、図3に示されるように、玄関(部屋Aとも言う)301に配置され、玄関301内の領域311Aを検知領域として割り当てられている。同様に、人体検知センサ101B~101Eは、夫々台所(部屋Bとも言う)302、洋室303の北側(部屋Cとも言う)、洋室303の南側(部屋Dとも言う)及び洗面所(部屋Eとも言う)304に配置され、夫々領域311B~311Eを検知領域として割り当てられている。洋室303では、検知領域311C及び311Dが部分的に重なり合っている。検知領域311C及び311Dが重なり合う領域は、仮想的な部屋Gに設定される。また、浴室(部屋Fとも言う)305には、人体検知センサが設置されていない。
図3に示される矢印は、移動可能な部屋間の接続関係を示している。例えば、部屋Aにいる人は、部屋B、又は人数を推定する対象外となる屋外(部屋X、領域外とも言う)のいずれかに移動することができる。図3の間取りにおける部屋間の接続関係情報109及び屋外との接続関係情報110の例が夫々図4及び図5に示されている。図4に示される部屋間の接続関係情報109の各フィールドには、部屋間の移動が可能、即ち、部屋同士が隣接している場合、1が格納され、部屋同士が非隣接である場合、0が格納される。図4に示される部屋間の接続関係情報109は、部屋Aからは、部屋Bに移動することができ、且つ部屋Aに滞在することができることを示している。同様に、図5に示される屋外との接続関係情報110の各フィールドには、部屋と屋外とが隣接している場合、1が格納され、部屋と屋外とが非隣接である場合、0が格納される。図5に示される屋外との接続関係情報110は、部屋Aから屋外(部屋X)へ、及び屋外から部屋Aへ移動することができることを示している。
尚、部屋間の接続関係及び部屋と屋外との接続関係は、別々の情報として用意されてもよく、或いは同一の接続関係情報にまとめて用意されてもよい。
図6は、図3に示される間取りにおける人体検知センサの重複領域に関する重複関係情報112の一例を示している。図6に示されるように、重複関係情報112は、対角要素(即ち、行と列のインデックスが同じ要素)が1の場合、そのインデックスで特定される部屋に人体検知センサが設置されていることを示し、対角要素が0の場合、そのインデックスにより特定される部屋に人体検知センサが非設置であることを示す。また、非対角要素が1の場合、そのインデックスで特定される2つの部屋の検知エリアが互いに重なり合っていることを示す。図6に示される重複関係情報112は、部屋A~Eには人体検知センサが設置され、部屋F及びGには人体検知センサが非設置であり、部屋Cと部屋Dと検知エリアが重複し、重複エリアを仮想的な部屋Gとした場合を示している。
以下では、図3に示される住居を具体例として適宜参照しながら、存在人数推定装置100の具体的な処理手順を説明する。
図7は、存在人数推定装置100の処理手順の一例を示している。図7に示されるように、存在人数推定装置100が動作開始されると、制御部103は、ステップS701において、実現値の集合111を初期化する。初期化時に、実際の在室状況、即ち、人体検知情報から得られる在室状況と大きく異なって実現値が設定された場合、正しい推定値を計算する処理に時間がかかるため、制御部103は、例えば、次に示す方法で実現値の集合111を初期化する。即ち、制御部103は、住居全体の人数及び各部屋の最大人数等の予め設定された制約条件を満たすように乱数で複数の実現値を生成し、或いはこの制約条件を満たす全ての状態を実現値として生成する。
ステップS702では、制御部103は、検知情報収集部102から検知状態情報を受信するか、或いは内蔵されたタイマーによって定期的に呼び出されるまで処理を一時停止する。一時停止復帰後、ステップS703において、図8に示される一期先予測処理が実行され、ステップS704において、図14に示される実現値評価処理が実行され、さらに、ステップS705において在室人数出力処理が実行される。ステップS702~S705に示される一連の処理は、制御部103が終了を指示するまで、繰り返し実行され(ステップS706)、実現値の集合111が次々に更新される。
図8は、図7のステップS703に示される一期先予測処理の手順の一例を示している。図8のステップS801において、遷移行列作成部104Aは、部屋間の接続関係情報109及び屋外との接続関係情報110から遷移行列を作成する。ステップS802において、実現値予測部104Bは、作成された遷移行列を使用して一時刻前の各実現値から現在時刻の実現値を予測して、予測した実現値で実現値の集合111を更新する。
図9は、図8のステップS801に示される遷移行列作成処理の手順の一例を示している。図9のステップS901において、遷移行列作成部104Aは、ある部屋(移動元の部屋とも言う)iから他の部屋(移動先の部屋とも言う)jに人が移動する場合の遷移確率Pijを計算する。遷移確率Pijは、隣室へ移動する確率をa、屋外へ移動する確率をbとし、部屋iと部屋jとの接続関係をdij、屋外と部屋i(部屋j)との接続関係をgi(gj)とすると、下記数式1から計算される。移動確率aは、その値が大きいほど他の部屋への移動が起こり易いことを表わし、移動確率bは、その値が大きいほど屋外への移動及び屋外からの移動が起こり易いことを表わす。遷移確率Pijは、移動元の部屋i以外の全ての部屋jについて計算される(ステップS902)。
ここで、i≠jである。
ステップS903において、一期先予測部104は、移動元の部屋iに滞在する場合の遷移確率Piiを計算する。遷移確率Piiは、ステップS901及びS902で算出された遷移確率Pijを使用して下記数式2から計算される。
ステップS901~S903に示される計算は、全ての移動元の部屋iに対して実行され(ステップS904)、ステップS901~S904で計算された遷移確率を行列要素に有する遷移行列Pが作成される。図10には、図4及び図5に示される接続関係情報109及び110を使用して、a=0.1、b=0.01として計算された遷移行列Pが示されている。
図11は、実現値の集合111内の実現値のデータ構造例を示している。図11に示される例では、実現値の集合111には、4つの実現値が夫々実現値1から実現値4とインデックスを付されて格納されている。実現値の集合111に格納される実現値の数Nは、部屋の数等の状況に応じて適宜変更されることができる。
図12は、図8のステップS802に示される実現値予測処理の手順を示している。図12のステップ1201において、実現値予測部104Bは、ステップS801で計算された遷移行列Pを使用して、一時刻前のある実現値k(ここで、kは、1≦k≦Nを満たす整数)における部屋iに在室する人間が現在時刻においてどの部屋jに移動するか、或いは屋外に移動するかを決定する。より具体的には、遷移行列のi行目の横ベクトルが部屋iから各部屋に移動する確率になっているため、この確率に従って乱数を発生させて1つの移動先の部屋jをサンプリングする。例えば、部屋Gに在室する人に関しては、図10に示されるように、部屋Gに滞在する確率が0.8、部屋C及び部屋Dに移動する確率が各々0.1であり、この人が部屋Gに滞在するか、部屋C又は部屋Dに移動するかをこの確率に従って乱数で決定する。
ステップS1202において、実現値予測部104Bは、ステップS1201で決定された移動先の部屋j(部屋jは、部屋iと同一であってもよい)に従って、移動元の部屋iの人数を1減少させ、移動先の部屋jの人数を1増加させる。ここで、ステップS1201で決定された移動先の部屋が部屋Xである場合は、移動元の部屋iの人数を1減少させる。ステップS1201及びS1202における処理は、一時刻前の実現値kにおける移動元の部屋iの在室人数の回数だけ繰り返される(ステップS1203)。ステップS1204では、ステップS1201~S1203に示される処理が全ての部屋、例えば、部屋Aから部屋Fに対して繰り返し実行される。次に、ステップS1205に示されるように、全ての実現値に対してステップS1201~S1204に示される処理が繰り返し実行される。
図13は、図11に示した実現値の集合111に対して実現値予測処理を実行した後の実現値の集合111の例を示している。図11及び図13に示されるように、実現値1においては、部屋Gから部屋Dに1人移動し、部屋Eから部屋Fに1人移動している。また、実現値2においては、部屋Eから部屋Fに1人移動し、実現値3においては、部屋Dから部屋Gに1人移動している。さらに、実現値4においては、部屋Gから部屋Dに1人移動し、部屋Eから部屋Fに1人移動している。
上述のような処理に従って、一期先予測部103は、一時刻前の実現値から現在時刻の実現値を予測し、実現値の集合111を更新する。
尚、図8のステップS801に示される遷移行列作成処理は、必ずしも毎回実行される必要はなく、遷移行列Pは、予め計算されて外部記憶部207等に格納されてもよい。
図14は、図7のステップS704に示される実現値評価処理の手順の一例を示している。実現値評価部105は、尤度計算部105A及び実現値評価部105を実行し、一期先予測部104によって更新された実現値の集合111内の各実現値が人体検知情報に基づいた部屋毎の在室状況と整合するか否かを評価する。図14のステップS1401において、尤度計算部105Aは、重複関係情報112及び人体検知情報を使用して、実現値の集合111内の各実現値の尤度wを計算する尤度計算処理を実行する。実現値kの尤度wは、下記数式3により計算される。
ここで、実現値kの各部屋の在室人数から計算される部屋i内の人間の在不在を示す値(第1の在不在値とも言う)siと、人体検知情報で特定される各部屋の在室状況から計算される部屋i内の人間の在不在を示す値(第2の在不在値とも言う)riとの差が、平均μが0、標準偏差がσの正規分布N(x;μ=0,σ)に従うとして尤度wを計算している。正規分布N(x;μ=0,σ)は、下記数式4のように定義されている。
また、Hjjは、重複関係情報112を示しており、従って、検知領域の重複による仮想的な部屋(例えば、部屋G)及び人体検知センサが未設置の部屋(例えば、部屋F)を除いた部屋iが尤度wの計算対象になる。値siは、下記数式5により計算される。
この数式5は、実現値kの部屋iの人数xiを在(1)不在(0)に変換している。例えば、図3に示される実現値1では、sD=sF=1、sA=sB=sC=sE=0となる。また、値riは、下記数式6により計算される。
ここで、j∈adj(i)は、部屋iの隣室jを表わし、yiは、部屋iにおける最終検知時刻を表わしている。最終検知時刻yiから現在時刻tまでの経過時間t-yiは、部屋iでの未検知時間に一致する。関数fは、未検知時間t-yiに応じて部屋iに人間が在室している確率を小さくするように、例えば、下記数式7のような値が1から0に単調減少する関数として与えられる。
ここで、パラメータαiは、部屋iに1人も在室していないとみなすまでの未検知時間を表わし、パラメータβiは、部屋iからセンサ未設置の部屋及び屋外等(以下、遷移先の部屋と言う)に移動したとみなすまでの未検知時間を表わしている。パラメータαi及びβiは、部屋毎に異なる値を設定することができる。例えば、廊下(部屋Sとも言う)では、在室しているにもかかわらず長時間人間検知センサに反応がない可能性は低いが、寝室(部屋Tとも言う)では、睡眠時等、検知領域内に人が存在するが動きがほとんどないことがあり、在室しているにもかかわらず長時間反応がない可能性は高い。従って、αSには小さい値が設定され、αTには大きい値が設定される。
数式6では、3つの場合に場合分けされている。部屋i及びその隣室jのうち部屋iで最後に動いている人間の有無が検知され、且つ隣室に遷移先の部屋が含まれていない場合に、値riが数式6の上式で計算される。この場合、人が隣室に移動しているとは考えられないので、必ず部屋iに在室し、即ち、常にri=1になる。また、部屋i及び隣室のうち部屋iで最後に動いている人間の有無が検知され、且つ隣室に1以上の遷移先の部屋が存在する場合に、値riが数式6の中式で計算される。この場合、人間が部屋iから遷移先の部屋に移動している可能性があるため、0≦ri≦1になる。さらに、部屋iより後に隣室で動いている人間の有無が検知される場合、値riが数式6の下式で計算される。この場合、人間が部屋iから隣室に移動している可能性が高いため、0≦ri≦0.5になる。
図13に示される実現値に対してステップS1401の尤度計算処理を実行して得られた計算結果の一例が図15に示されている。図15において、各実現値の尤度wは、全ての実現値の尤度の合計が1となるように規格化され、尤度の比として示されている。
ここで、xi
newは、一期先予測処理で計算された実現値kの部屋iの人数を表し、xi
oldは、一期先予測処理で計算される前の実現値kの部屋iの人数を表している。また、rは、0から1の一様乱数とする。
図14のステップS1402において、実現値評価部105は、ステップS1401で計算された各実現値の尤度wを比較して、尤度wの小さい実現値を消去し、尤度wの大きい実現値を複製する実現値選択処理を実行する。ステップS1402の実現値選択処理は、一例として図16に示されている。図16に示されるステップS1601において、実現値選択部105Bは、ステップS1401で計算した各実現値の尤度wから、実現値1から実現値kまでの尤度wの比の合計である累積尤度比vkを下記数式9のように計算する。
ここで、Nは、実現値の集合111内に格納されている実現値の数であり、wmは、実現値mの尤度wを表わす。ステップS1502では、実現値選択部105Bは、vk-1<n-ε≦vkを満たす整数nの数だけ実現値kを複製して保存する。ここで、パラメータεは、0を超え、1未満の値であり、予め設定される。は、vk-1<n-ε≦vkを満たす整数nが存在しない場合は、この実現値kを消去する。これらの処理によって、実現度選択部105Bは、尤度wの低い実現値を消去し、尤度wの高い実現値を複製している。図16に示される実現値選択処理は、一例であって、他の方法で消去及び複製する実現値を選択してもよい。例えば、尤度の比に関して予め閾値が設定され、ステップS1402で各実現値に対して計算される尤度の比が閾値以上となる場合にその実現値を複製し、尤度の比が閾値未満となる場合にその実現値を消去するようにしてもよい。
図17は、図13に示された実現値に対して実現値選択処理を実行した結果の一例を示している。図13に示される各実現値の累積尤度比は、図15に示される値を使用すると、v1=0.96、v2=3.00、v3=3.04、V4=4.00と求まる。この結果、例えばパラメータεが0.1と設定される場合、図17に示されるように、図13に示される実現値3が削除され、図13に示される実現値2が複製されて実現値2及び実現値3とされる。
上述のように、図7に示されるステップS703の一期先予測処理、及びステップS704の実現値評価処理が順次実行されて、実現値の集合111が更新される。ステップS705において、在室人数出力部106は、更新された実現値の集合111から部屋毎の在室人数に関する出力情報を表示装置107及び制御対象装置108に出力する。在室人数出力部106は、図18に示されるように、全ての部屋の在室人数が同一の実現値の数を計算し、数の多い順序に並び替えた実現値を出力情報として出力することができる。このように、同一の実現値の総数とともに実現値を出力することにより、制御対象装置108が、総数の多い実現値を使用したり、上位10個の実現値の平均を取る等の種々の処理を実行することが可能となる。また、在室人数出力部106は、図19に示されるように、N個の実現値から部屋毎に在室人数の平均値を計算し、部屋毎の平均在室人数を代表値として出力してもよい。代表値のように1つのデータを出力することにより、在室人数出力部106からの出力情報を他の手法で推定した部屋毎の在室人数に関する出力情報との置き換えが容易になる。また、出力情報が人間にとって理解しやすくなるという利点がある。
在室人数に応じて制御される制御対象機器108としては、例えば、照明、エアコン、扇風機、空気清浄機、テレビ、パーソナルコンピュータ等がある。例えば、図3に示される部屋303に設置される照明においては、部屋C、D及びGの合計の在室人数が0.2人以下と推定された場合にオフするように制御される。また、エアコン、扇風機、空気清浄機等は、在室人数に応じて強弱、又は運転モード等が制御されることができる。このように在室人数に応じて種々の機器の動作を制御することにより、不必要な電力消費を抑えることができる。
尚、在室人数出力部106が図示しない制御装置に出力情報を出力し、この制御装置が部屋の在室人数に応じて上記の制御対象機器108の動作を制御してもよい。
また、人数を推定する対象となる部屋(即ち、領域)は、図3に示した例に限定されず、リビング、トイレ、寝室、オフィスにおける会議室等を含むことができる。
以上のように、第1の実施の形態に係る存在人数推定装置は、部屋毎の存在人数を示す複数の実現値を用意し、これら実現値から現在時刻の実現値を予測し、予測した現在時刻の実現値から人体検知情報に基づいた在室状況と整合する実現値を尤度計算により選択している。この存在人数推定装置は、予測した実現値から実際の在室状況と整合性の高い実現値を選択することにより、部屋(領域)毎の存在人数を精度良く推定することができる。
[第2の実施の形態]
図20から図22を参照して、第2の実施の形態に係る存在人数推定装置を説明する。第2の実施の形態に係る存在人数推定装置100は、図20に簡略化して示されているが、図1に示される第1の実施の形態と同一の構成を有している。第2の実施の形態では、第1の実施の形態と比較して、人体検知センサの配置が変更されている。
図20から図22を参照して、第2の実施の形態に係る存在人数推定装置を説明する。第2の実施の形態に係る存在人数推定装置100は、図20に簡略化して示されているが、図1に示される第1の実施の形態と同一の構成を有している。第2の実施の形態では、第1の実施の形態と比較して、人体検知センサの配置が変更されている。
図21は、第2の実施の形態における人体検知センサ101A及び101Bに割り当てられる検知領域の一例を示している。人体検知センサ101Bは、部屋2102内の領域2111Bを検知領域として設定され、人体検知センサ101Aは、部屋2101及び部屋2102間の出入り口(部屋A)を含む領域2111Aを検知領域として設定されている。この場合、人体検知センサ101Aは、部屋2101及び部屋2102間の出入り口付近に配置される。他の例では、人体検知センサ101A及び101Bは、図22に示されるように、部屋2102の天井に設置された照明に組み込まれて設けられ、センサ101Aが部屋の出入り口を検知領域に設定され、センサ101Bが部屋内を検知領域に設定されてもよい。このように、1以上の人体検知センサがエアコン、照明等の他の1つの機器に一体化して設けられてもよい。
以上のように、第2の実施の形態に係る存在人数推定装置においては、人体検知センサの検知領域が出入り口に割り当てられることから、人間の部屋間の移動を確実に検知することができ、部屋(領域)毎の存在人数をより精度良く推定することができる。
本発明に係る存在人数推定装置は、住宅、オフィス等の屋内に存在する人数を領域毎に精度良く推定することができ、推定された存在人数に応じて照明、エアコン、テレビ等の機器制御が可能となることから産業上の利用可能性が高い。
100…存在人数推定装置、101…人体検知センサ群、101A~101E…人体検知センサ、102…検知状態収集部、103…制御部、104…一期先予測部、104A…遷移行列生成部、104B…実現値予測部、105…実現値評価部、105A…計算部、105B…実現値選択部、106…在室人数出力部、107…表示装置、108…制御対象機器、109…部屋間の接続関係情報、110…屋外との接続関係情報、111…実現値の集合
Claims (5)
- 人数を推定する対象となる所定の範囲の複数の領域のいずれかに複数個設置され、前記領域内の人間の有無又は前記領域内の人間の動作の有無を検知して検知信号を出力する人体検知センサと、
前記検知信号を収集し、前記領域で前記人間の有無又は前記人間の動作の有無が最後に検知された時刻を前記領域毎に関連づけた時刻情報を含む人体検知情報を出力する収集部と、
前記領域毎の推定存在人数を示す複数の実現値と、前記領域間の接続関係及び前記領域と前記人数を推定する対象外となる領域外との接続関係を含む接続関係情報と、前記領域が重なり合う重複領域を示す重複関係情報と、を格納する記憶部と、
前記接続関係情報を使用して、前記領域間を人間が遷移する確率を示す第1の遷移確率、及び領域から領域外へ人間が遷移する確率を示す第2の遷移確率を計算し、前記第1及び第2の遷移確率を行列要素に有する遷移行列を作成する遷移行列作成部と、
前記遷移行列を使用して前記実現値の各々から現在時刻の実現値を予測する実現値予測部と、
前記時刻情報及び前記重複関係情報を使用して前記現在時刻の実現値の各々の尤度を計算する尤度計算部と、
前記現在時刻の各実現値のうち前記尤度が予め設定される閾値よりも大きい実現値を選択する実現値選択部と、
前記選択された実現値から同一の実現値の数を計算し、該実現値の数と前記領域毎の存在人数とを関連づけた出力情報を生成して、この出力情報を出力する出力部と、
を具備することを特徴とする存在人数推定装置。 - 前記遷移行列作成部は、前記第1及び第2の遷移確率から、同一の領域に滞在する確率を示す滞在確率をさらに計算し、
前記遷移行列は、前記滞在確率を行列要素に有することを特徴とする請求項1に記載の存在人数推定装置。 - 前記尤度計算部は、前記現在時刻の実現値から前記領域毎の第1の在不在値を計算し、前記人体検知信号に基づいて前記領域毎の第2の在不在値を計算して、これら第1及び第2の在不在値を使用して尤度を計算することを特徴とする請求項1に記載の存在人数推定装置。
- 前記出力情報は、前記領域毎の平均の推定存在人数を含むことを特徴とする請求項1に記載の存在人数推定装置。
- 前記出力情報は、通信可能に接続される機器の電源、運転モード、又は運転の強弱の制御に使用されることを特徴とする請求項1に記載の存在人数推定装置。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2009/065432 WO2011027452A1 (ja) | 2009-09-03 | 2009-09-03 | 存在人数推定装置 |
JP2011529742A JP5162030B2 (ja) | 2009-09-03 | 2009-09-03 | 存在人数推定装置 |
US13/410,332 US9445164B2 (en) | 2009-09-03 | 2012-03-02 | Existent person count estimation apparatus |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2009/065432 WO2011027452A1 (ja) | 2009-09-03 | 2009-09-03 | 存在人数推定装置 |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/410,332 Continuation US9445164B2 (en) | 2009-09-03 | 2012-03-02 | Existent person count estimation apparatus |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2011027452A1 true WO2011027452A1 (ja) | 2011-03-10 |
Family
ID=43649014
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2009/065432 WO2011027452A1 (ja) | 2009-09-03 | 2009-09-03 | 存在人数推定装置 |
Country Status (3)
Country | Link |
---|---|
US (1) | US9445164B2 (ja) |
JP (1) | JP5162030B2 (ja) |
WO (1) | WO2011027452A1 (ja) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014183553A (ja) * | 2013-03-21 | 2014-09-29 | Toshiba Corp | 在室確率推定装置およびその方法、ならびにプログラム |
JP2015507263A (ja) * | 2011-12-19 | 2015-03-05 | マイクロソフト コーポレーション | 複数センサ入力のためのセンサフュージョンインタフェース |
JP2015153156A (ja) * | 2014-02-14 | 2015-08-24 | Kddi株式会社 | 消費電力量に基づく整合性のある生活状況の推定が可能な装置、プログラム及び方法 |
JP7550522B2 (ja) | 2020-03-13 | 2024-09-13 | 三菱電機ビルソリューションズ株式会社 | 人監視システム、及び人監視方法 |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014120180A1 (en) * | 2013-01-31 | 2014-08-07 | Hewlett-Packard Development Company, L.P. | Area occupancy information extraction |
US9689583B2 (en) * | 2013-09-10 | 2017-06-27 | Honeywell International Inc. | Occupancy based energy optimization systems and methods |
US20150088272A1 (en) * | 2013-09-23 | 2015-03-26 | Emerson Electric Co. | Energy Management Based on Occupancy and Occupant Activity Level |
JP6428144B2 (ja) * | 2014-10-17 | 2018-11-28 | オムロン株式会社 | エリア情報推定装置、エリア情報推定方法、および空気調和装置 |
US10270609B2 (en) * | 2015-02-24 | 2019-04-23 | BrainofT Inc. | Automatically learning and controlling connected devices |
KR102051855B1 (ko) | 2018-01-02 | 2019-12-04 | 한양대학교 산학협력단 | Uwb 레이더를 이용하는 인원 계수 방법 및 장치 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0948565A (ja) * | 1995-08-08 | 1997-02-18 | Toshiba Corp | 居住空間監視制御装置 |
JP2000035243A (ja) * | 1998-07-15 | 2000-02-02 | Yamatake Corp | ビル管理システム |
JP2008077361A (ja) * | 2006-09-20 | 2008-04-03 | Kanazawa Inst Of Technology | 監視方法および監視システム |
JP2008217351A (ja) * | 2007-03-02 | 2008-09-18 | Mitsubishi Electric Corp | 移動量算出装置 |
JP2008298353A (ja) * | 2007-05-31 | 2008-12-11 | Mitsubishi Electric Building Techno Service Co Ltd | 建屋内空調管理システム |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7139409B2 (en) * | 2000-09-06 | 2006-11-21 | Siemens Corporate Research, Inc. | Real-time crowd density estimation from video |
US7596241B2 (en) * | 2005-06-30 | 2009-09-29 | General Electric Company | System and method for automatic person counting and detection of specific events |
US7868812B2 (en) * | 2005-08-12 | 2011-01-11 | Patriot Technologies, Llp | Surveillance and warning system |
US8098485B2 (en) * | 2007-07-03 | 2012-01-17 | 3M Innovative Properties Company | Wireless network sensors for detecting events occurring proximate the sensors |
CA2700342A1 (en) * | 2007-09-19 | 2009-03-26 | United Technologies Corporation | System and method for occupancy estimation |
JP4908541B2 (ja) | 2009-03-31 | 2012-04-04 | 株式会社東芝 | 在室人数推定装置及び在室人数推定方法 |
-
2009
- 2009-09-03 WO PCT/JP2009/065432 patent/WO2011027452A1/ja active Application Filing
- 2009-09-03 JP JP2011529742A patent/JP5162030B2/ja active Active
-
2012
- 2012-03-02 US US13/410,332 patent/US9445164B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0948565A (ja) * | 1995-08-08 | 1997-02-18 | Toshiba Corp | 居住空間監視制御装置 |
JP2000035243A (ja) * | 1998-07-15 | 2000-02-02 | Yamatake Corp | ビル管理システム |
JP2008077361A (ja) * | 2006-09-20 | 2008-04-03 | Kanazawa Inst Of Technology | 監視方法および監視システム |
JP2008217351A (ja) * | 2007-03-02 | 2008-09-18 | Mitsubishi Electric Corp | 移動量算出装置 |
JP2008298353A (ja) * | 2007-05-31 | 2008-12-11 | Mitsubishi Electric Building Techno Service Co Ltd | 建屋内空調管理システム |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015507263A (ja) * | 2011-12-19 | 2015-03-05 | マイクロソフト コーポレーション | 複数センサ入力のためのセンサフュージョンインタフェース |
JP2017215970A (ja) * | 2011-12-19 | 2017-12-07 | マイクロソフト テクノロジー ライセンシング,エルエルシー | 複数センサ入力のためのセンサフュージョンインタフェース |
JP2014183553A (ja) * | 2013-03-21 | 2014-09-29 | Toshiba Corp | 在室確率推定装置およびその方法、ならびにプログラム |
JP2015153156A (ja) * | 2014-02-14 | 2015-08-24 | Kddi株式会社 | 消費電力量に基づく整合性のある生活状況の推定が可能な装置、プログラム及び方法 |
JP7550522B2 (ja) | 2020-03-13 | 2024-09-13 | 三菱電機ビルソリューションズ株式会社 | 人監視システム、及び人監視方法 |
Also Published As
Publication number | Publication date |
---|---|
US20120209567A1 (en) | 2012-08-16 |
JPWO2011027452A1 (ja) | 2013-01-31 |
US9445164B2 (en) | 2016-09-13 |
JP5162030B2 (ja) | 2013-03-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5162030B2 (ja) | 存在人数推定装置 | |
EP3832556B1 (en) | Action optimization device, method and program | |
Tekler et al. | A scalable Bluetooth Low Energy approach to identify occupancy patterns and profiles in office spaces | |
Candanedo et al. | Accurate occupancy detection of an office room from light, temperature, humidity and CO2 measurements using statistical learning models | |
Jung et al. | Human-in-the-loop HVAC operations: A quantitative review on occupancy, comfort, and energy-efficiency dimensions | |
WO2017216511A1 (en) | Estimating a number of occupants in a region | |
JP5406006B2 (ja) | 人数推定装置及び人数推定方法 | |
KR102472214B1 (ko) | 에어 컨디셔닝 시스템에서 복합 제어 장치 및 방법 | |
Dong et al. | Sensing and data acquisition | |
US9116515B2 (en) | In-room probability estimating apparatus, method therefor and program | |
JP6654091B2 (ja) | 監視装置、監視方法、およびプログラム | |
JP2008269215A (ja) | 特異パターン検出システム、モデル学習装置、特異パターン検出方法、及び、コンピュータプログラム | |
WO2009038561A1 (en) | System and method for threat propagation estimation | |
US20130006899A1 (en) | Activity Recognition in Multi-Entity Environments | |
WO2019199919A1 (en) | A system for monitoring occupancy and activity in a space | |
Gunay et al. | Detecting occupants’ presence in office spaces: a case study | |
JP2007249647A (ja) | 環境制御システム及び環境制御方法 | |
Rajabi et al. | MODES: M ulti-sensor o ccupancy d ata-driven e stimation s ystem for smart buildings | |
Chu et al. | Influential variables impacting the reliability of building occupancy sensor systems: A systematic review and expert survey | |
Mashuk et al. | Demonstrating the potential of indoor positioning for monitoring building occupancy through ecologically valid trials | |
Andersen et al. | Exploring occupant detection model generalizability for residential buildings using supervised learning with IEQ sensors | |
Carvalho et al. | Adaptive tracking model in the framework of medical nursing home using infrared sensors | |
WO2023010220A1 (en) | Determining spatial maps based on user input and motion-sensing data derived from wireless signals | |
JP4908541B2 (ja) | 在室人数推定装置及び在室人数推定方法 | |
De Simone et al. | Occupant behavior: a †œnew†factor in energy performance of buildings. Methods for its detection in houses and in offices. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 09848983 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2011529742 Country of ref document: JP |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 09848983 Country of ref document: EP Kind code of ref document: A1 |