WO2021117343A1 - 空間温度推定システム、温冷感推定システム、空間温度推定方法、温冷感推定方法、及び、プログラム - Google Patents

空間温度推定システム、温冷感推定システム、空間温度推定方法、温冷感推定方法、及び、プログラム Download PDF

Info

Publication number
WO2021117343A1
WO2021117343A1 PCT/JP2020/039188 JP2020039188W WO2021117343A1 WO 2021117343 A1 WO2021117343 A1 WO 2021117343A1 JP 2020039188 W JP2020039188 W JP 2020039188W WO 2021117343 A1 WO2021117343 A1 WO 2021117343A1
Authority
WO
WIPO (PCT)
Prior art keywords
temperature
room
estimation
person
floor plan
Prior art date
Application number
PCT/JP2020/039188
Other languages
English (en)
French (fr)
Inventor
延亮 島本
那由多 南
Original Assignee
パナソニックＩｐマネジメント株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by パナソニックＩｐマネジメント株式会社 filed Critical パナソニックＩｐマネジメント株式会社
Priority to CN202080084044.0A priority Critical patent/CN114829871A/zh
Priority to US17/783,529 priority patent/US20220398764A1/en
Priority to JP2021563767A priority patent/JPWO2021117343A1/ja
Publication of WO2021117343A1 publication Critical patent/WO2021117343A1/ja

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/70—Determining position or orientation of objects or cameras
    • 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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00—Measuring arrangements characterised by the use of optical techniques
    • G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J5/0022—Radiation pyrometry, e.g. infrared or optical thermometry for sensing the radiation of moving bodies
    • G01J5/0025—Living bodies
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J5/02—Constructional details
    • G01J5/026—Control of working procedures of a pyrometer, other than calibration; Bandwidth calculation; Gain control
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J5/02—Constructional details
    • G01J5/0275—Control or determination of height or distance or angle information for sensors or receivers
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F24—HEATING; RANGES; VENTILATING
    • F24F—AIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
    • F24F2110/00—Control inputs relating to air properties
    • F24F2110/10—Temperature
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F24—HEATING; RANGES; VENTILATING
    • F24F—AIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
    • F24F2110/00—Control inputs relating to air properties
    • F24F2110/30—Velocity
    • 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
    • F24F2120/12—Position of occupants
    • 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
    • F24F2120/14—Activity of occupants
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J2005/0077—Imaging
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10048—Infrared image

Definitions

  • the present disclosure relates to a space temperature estimation system, a temperature sensation estimation system, a space temperature estimation method, a temperature sensation estimation method, and a program. More specifically, the present disclosure describes a space temperature estimation system using a CRI (Contribution Ratio of Indoor climate) coefficient, a heating / cooling sensation estimation system using this space temperature estimation system, and a space temperature estimation method. Regarding the hot / cold feeling estimation method and the program.
  • CRI Consumer Ratio of Indoor climate
  • the air conditioner described in Patent Document 1 includes a thermal image acquisition unit (thermography), a calculation unit, and a control unit.
  • the thermal image acquisition unit acquires a thermal image showing the temperature distribution in space.
  • the calculation unit identifies a region corresponding to a person in the thermal image acquired by the thermal image acquisition unit.
  • the arithmetic unit determines the human body temperature, which is the temperature of the person in the space, based on the temperature distribution of the region corresponding to the person in the thermal image, and obtains it from the human body temperature and the temperature of the region other than the region corresponding to the person.
  • the feeling of warmth and coldness of a person in the space is estimated based on the difference value from the ambient temperature.
  • the control unit controls at least one of the air volume, air temperature, and wind direction of the air conditioner based on the feeling of warmth and coldness estimated by the calculation unit.
  • the temperature obtained from a region other than the region corresponding to the person is defined as the temperature of the space around the person.
  • the thermal image represents the temperature of the thermal factors (walls, ceilings, floors, etc.) reflected in the thermal image, and does not necessarily represent the temperature of the air in the space reflected in the thermal image. Therefore, in the air conditioner described in Patent Document 1, the temperature of the air in the space around the person may not be appropriately obtained.
  • An object of the present disclosure is to provide a space temperature estimation system, a temperature sensation estimation system, a space temperature estimation method, a temperature sensation estimation method, and a program capable of estimating the temperature of air in an indoor space by simple calculation. To do.
  • the space temperature estimation system of one aspect of the present disclosure includes a temperature detecting means and a space temperature estimating means.
  • the temperature detecting means detects the temperature of at least one surface of the ceiling, floor and a plurality of walls of a room.
  • the space temperature estimating means estimates the temperature of air in the indoor space of the room by using the CRI coefficient for at least one surface and the temperature detected by the temperature detecting means.
  • the warm / cold sensation estimation system of one aspect of the present disclosure includes the space temperature estimation system, a human position detecting means, and a warm / cold sensation index calculating means.
  • the person position detecting means detects the position of a person in the room.
  • the hot / cold feeling index calculation unit calculates an index representing the hot / cold feeling of the person by using the temperature of the air at the position of the person estimated by the space temperature estimation system.
  • the space temperature estimation method of one aspect of the present disclosure includes a temperature detection step and a space temperature estimation step.
  • the temperature detection step detects the temperature of at least one surface of the ceiling, floor and a plurality of walls of the room.
  • the space temperature estimation step estimates the temperature of air in the indoor space of the room by using the CRI coefficient for at least one surface and the temperature detected by the temperature detecting means.
  • the heating / cooling sensation estimation method of one aspect of the present disclosure is a heating / cooling sensation estimation method using the space temperature estimation method.
  • the hot / cold feeling estimation method includes a human position detection step and a hot / cold feeling index calculation step.
  • the person position detection step detects the position of a person in the room.
  • an index representing the warm / cold sensation of the person is calculated using the temperature of the air at the position of the person estimated by the space temperature estimation method.
  • the program of one aspect of the present disclosure is a program for causing a computer system to execute the space temperature estimation method.
  • the program of one aspect of the present disclosure is a program for causing a computer system to execute the above-mentioned heating / cooling sensation estimation method.
  • FIG. 1 is a block diagram of a heating / cooling sensation estimation system according to an embodiment.
  • FIG. 2 is an explanatory diagram for explaining a method of estimating the space temperature of the space temperature estimation system provided in the above-mentioned temperature sensation estimation system.
  • FIG. 3 is an explanatory diagram for explaining the installation position of the air conditioner in the room.
  • 4A to 4C are explanatory views for explaining the composition of the thermal image when the installation positions of the air conditioner of the same are the "center", "left end", and "right end", respectively.
  • FIG. 5 is a flowchart illustrating the operation of the above-mentioned heating / cooling sensation estimation system.
  • FIG. 6 is a block diagram of the heating / cooling sensation estimation system according to the first modification.
  • FIG. 7A to 7C are block diagrams of the heating / cooling sensation estimation system according to the first modification.
  • FIG. 8 is a block diagram of the heating / cooling sensation estimation system according to the second modification.
  • FIG. 9A is an explanatory diagram illustrating an example of a thermal image.
  • FIG. 9B is an explanatory diagram illustrating an example of the restored image.
  • FIG. 10 is a block diagram of the heating / cooling sensation estimation system according to the third modification.
  • FIG. 11 is a block diagram of the heating / cooling sensation estimation system according to the modified example 4.
  • FIG. 12 is a block diagram of the heating / cooling sensation estimation system according to the modified example 5.
  • FIG. 13 is a block diagram of the heating / cooling sensation estimation system according to the modified example 6.
  • the space temperature estimation system 1 is a system capable of estimating the temperature distribution of the indoor space of the target room.
  • the temperature distribution in the indoor space is the temperature distribution of air at each position in the indoor space.
  • the thermal sensation estimation system 100 is a system capable of estimating the thermal sensation of a person in a target room by using the temperature distribution of the indoor space estimated by the space temperature estimation system 1.
  • the hot / cold feeling estimation system 100 controls an air conditioner that controls the temperature in the target room based on the estimated hot / cold feeling. As a result, the temperature of the target room is automatically controlled to a temperature at which the person in the room feels comfortable.
  • the space temperature estimation system 1 and the hot / cold feeling estimation system 100 are arranged outside the air conditioner 7, but may be arranged inside the air conditioner 7.
  • the target room may be a room in a building or a room in a moving body such as a car (for example, a bus).
  • the space temperature estimation system 1 can estimate the temperature distribution of the indoor space of the target room. More specifically, the space temperature estimation system 1 can estimate the temperature (air temperature) of the air at each position in the indoor space by using the CRI (Contribution Ratio of Indoor climate) coefficient.
  • CRI Consumer Ratio of Indoor climate
  • the concept of estimating the temperature (air temperature) of the air in the indoor space in the present embodiment will be described.
  • the effect of the air flow blown out from the air outlet for example, an air conditioner
  • the air temperature at a certain position in the indoor space can be calculated by adding the degree of thermal influence from thermal factors such as the floor plan of the room and the wall floor temperature of the room under the air flow. is there.
  • the thermal factor is a factor that becomes a heat source, such as the ceiling, floor, and wall of a room.
  • the airflow distribution in the room is determined using CFD (Computational Fluid Dynamics) analysis, and based on the airflow distribution, the degree of thermal influence that contributes to a certain position in the room from each thermal factor is determined. calculate. As a result, even if the temperature of each thermal factor changes, the temperature at a certain position can be calculated simply by adding the degree of thermal influence due to the temperature change of each thermal factor.
  • CFD Computer Fluid Dynamics
  • Equation 1 a total of six surfaces, the ceiling R1 and the floor R2 of the room R10 and the surrounding (for example, four) walls R3 to R6, are thermal factors.
  • ⁇ b is a reference temperature at the position X, and is set to, for example, 0 (zero) degrees in the present embodiment.
  • ⁇ x is the change temperature ⁇ x at the position X. That is, the air temperature ⁇ x at the position X is calculated by the sum of the reference temperature ⁇ b at the position X and the change temperature ⁇ x at the position X.
  • the change temperature ⁇ x at the position X is the temperature due to the thermal influence from each thermal factor.
  • the change temperature ⁇ x at the position X is given by Equation 2.
  • Each of ⁇ 0, 1 to ⁇ 0, 6 is the change temperature of each thermal factor (ceiling R1, floor R2, and four walls R3 to R6 of the chamber R10).
  • C x, 1 to C x, 6 are CRI coefficients for each thermal factor with respect to the position X, respectively, and represent the thermal contribution rate of the change temperature from each thermal factor to the position X.
  • the CRI coefficient for each thermal factor is a value determined when the floor plan of the chamber R10 and the airflow distribution in the chamber R10 are determined.
  • the degree of thermal influence contributing to the position X is calculated from the thermal factor by the product of the change temperature of the thermal factor and the CRI coefficient related to the thermal factor. That is, the change temperature ⁇ x at the position X is calculated by the sum of the degrees of thermal influence contributing to the position X from each thermal factor.
  • the degree of thermal influence that contributes to the position X from each thermal factor is the change temperature ⁇ 0,1 to ⁇ 0,6 of the thermal factor, and the CRI coefficients C x, 1 to C x, related to the thermal factor, respectively. It is calculated by the product of 6.
  • the airflow distribution in the chamber R10 is determined by the floor plan of the chamber R10 and the installation position of the air conditioner 7. Then, the CRI coefficient for each thermal factor is determined from the floor plan of the chamber R10 and the airflow distribution in the chamber R10 (airflow distribution of the air blown by the air conditioner 7).
  • the change temperature ⁇ 0,1 ⁇ ⁇ 0,6 of each heat factors measured, CRI for each thermal factors From the coefficient, the change temperature ⁇ x at a certain position X is determined. From this change temperature ⁇ x and the reference temperature ⁇ b , the temperature at a certain position X in the room can be obtained.
  • the configuration of the thermal sensation estimation system 100 including the space temperature estimation system 1 will be described in detail. In the following description, it is assumed that the temperature of the indoor space of the chamber R10 shown in FIG. 2 is estimated.
  • the hot / cold sensation estimation system 100 includes a space temperature estimation system 1, a human position detection unit 2, an average radiation temperature calculation unit 3, a hot / cold sensation index calculation unit 4 (warm / cold sensation index calculation means), and a temperature target value. It includes a calculation unit 5, an air conditioning control unit 6, and an air conditioning device 7.
  • the space temperature estimation system 1 includes a setting input unit 21, a thermal image acquisition unit 22 (temperature detection means, thermal image acquisition means), a CRI coefficient determination unit 23, and a space temperature estimation unit 24 (spatial temperature estimation means). I have.
  • the air conditioner 7 is included in the configuration of the hot / cold sensation estimation system 100, but may not be included in the configuration of the hot / cold sensation estimation system 100.
  • the setting input unit 21 receives input of various setting information (floor plan information and installation position information) by the operation input of the operator.
  • the floor plan information is information regarding the floor plan of the target room R10.
  • the floor plan information is, for example, information on the floor surface shape (for example, vertically long rectangle, horizontally long rectangle or square) and floor area of the room R10. Since the capacity of the air conditioner 7 is determined according to the floor area of the room R10, in the present embodiment, information (rated power consumption, etc.) regarding the air conditioner capacity of the air conditioner 7 is acquired from the air conditioner 7 and automatically. Set.
  • the setting input unit 21 may be installed on the wall of the room R10, or may be configured as a remote control device. Further, it may also be used as a remote controller for the air conditioner 7.
  • the installation position information is information regarding the installation position of the air conditioner 7 in the room R10.
  • the installation position information is information on the height of the installation position of the air conditioner 7 (height from the floor R2) and the position of the wall R3 on which the air conditioner 7 is installed in the left-right direction (left end, center or right end).
  • the height of the installation position of the air conditioner 7 is approximately the same as the height of the chamber R10, and in the present embodiment, it is preset to a fixed value (for example, 2.2 m).
  • a fixed value for example, 2.2 m
  • the floor plan of the room R10, the installation position of the air conditioner 7, and the speed and direction of the air blown by the air conditioner 7 are collectively referred to as “floor plan conditions”.
  • the space temperature estimation system 1 automatically uses the thermal image acquired by the thermal image acquisition unit 22, for example, to automatically input the layout information. And installation information may be estimated. An example of this modification will be described later.
  • the CRI coefficient determination unit 23 identifies and specifies the floor plan conditions of the room R10 based on the floor plan information and the installation position information input to the setting input unit 21 and the information on the air velocity and the wind direction of the air conditioner 7. Determine the CRI coefficient corresponding to the floor plan conditions.
  • the floor plan condition of the room R10 is specified by the floor plan information and the installation position information input to the setting input unit 21, and the air velocity and the wind direction of the air conditioner 7.
  • the floor plan information includes information on the floor surface shape (vertical rectangle, horizontally long rectangle or square) of the room R10
  • the installation position information includes information on the installation position (center, left end or right end) of the air conditioner 7. ing.
  • Information on the blowing speed and the wind direction of the air conditioner 7 can be obtained from the air conditioner 7.
  • a typical wind speed and the wind direction are set in advance in the CRI coefficient determining unit 23, and the set values are used.
  • the number of possible floor plans of the chamber R10 is, for example, a total of three, which is a horizontally long rectangle, a vertically long rectangle, and a square, and the number of possible installation positions of the air conditioner 7 is, for example, a total of three at the left end, the right end, and the center. is there.
  • the CRI coefficient determination unit 23 stores a plurality of sets (9 sets in this embodiment) of CRI coefficients in a predetermined storage unit.
  • the plurality of sets of CRI coefficients have a one-to-one correspondence with the plurality of possible floor plan conditions of the chamber R10.
  • the CRI coefficient determining unit 23 specifies the floor plan condition of the room R10 based on the floor plan information and the installation position information input to the setting input unit 21 and the information on the blowing speed and the wind direction of the air conditioner 7.
  • the CRI coefficient determination unit 23 selects (determines) a set of CRI coefficients corresponding to the specified floor plan condition from the plurality of sets of CRI coefficients.
  • a set of CRI coefficients corresponding to the specified floor plan condition is selected from a plurality of sets of CRI coefficients prepared in advance. Selected (determined).
  • the processing load of estimating the space temperature of the chamber R10 can be reduced, and the CRI coefficient can be determined faster.
  • the CRI coefficient is changed according to the floor plan condition, the temperature of the air in the indoor space of the room R10 can be estimated more accurately.
  • the thermal image acquisition unit 22 acquires (that is, captures) a thermal image including one or more surfaces of each surface (ceiling R1, floor R2, and walls R3 to R6) of the chamber R10.
  • the thermal image acquisition unit 22 is installed in, for example, the air conditioner 7 or is installed around the air conditioner 7.
  • the thermal image is, for example, a temperature distribution image in which the temperature of a subject (thermal factor) reflected in the image is measured for each pixel, and each pixel is color-coded according to the temperature. From this thermal image, it is possible to measure the temperature change of each thermal factor (ceiling R1, floor R2 and walls R3 to R6 of the chamber R10).
  • the thermal image acquisition unit 22 may be, for example, an infrared array sensor.
  • the infrared array sensor is a sensor in which a plurality of infrared light receiving elements are arranged vertically and horizontally.
  • the thermal image acquisition unit 22 is a temperature detecting means for detecting the temperature of one or more of the surfaces (ceiling R1, floor R2, and walls R3 to R6) of the chamber R10.
  • each surface of the room R10 faces different directions from each other.
  • the thermal image acquisition unit 22 acquires a thermal image including four faces facing different directions, but acquires a thermal image including three or more faces or two or more faces facing different directions. May be good. It should be noted that the more surfaces included in the thermal image, the more the accuracy of estimating the temperature distribution in the indoor space of the chamber R10 can be improved.
  • the person position detection unit 2 detects the position X of a person in the room R10 based on the thermal image acquired by the thermal image acquisition unit 22. More specifically, the person position detection unit 2 identifies a region corresponding to a person in the thermal image based on the temperature corresponding to the person and the shape of the temperature distribution thereof, and based on the specified region, enters the chamber R10. Detects the position X of the person who is present. For example, a predetermined position in the specified area is used as the position X of a person.
  • the space temperature estimation unit 24 estimates the temperature (air temperature) of the air in the indoor space of the room R10. More specifically, the space temperature estimation unit 24 uses the change temperature of each thermal factor (ceiling R1, floor and R2 and four walls R3 to R6) of the room R10 and the detection result of the person position detection unit 2 (that is, the room). Using the position of the person in R10) and the CRI coefficient determined by the CRI coefficient determination unit 23, the temperature at the position of the person in the room R10 is set as the temperature of a certain position X in the indoor space of the room R10. presume.
  • the space temperature estimation unit 24 detects the change temperature of each thermal factor in the chamber R10 by using the thermal image acquired by the thermal image acquisition unit 22. Specifically, the space temperature estimation unit 24 specifies the floor plan (floor area and floor surface shape) of the room R10 based on the floor plan information. Then, based on the specified floor plan, the space temperature estimation unit 24 has a correspondence relationship between the ceiling R1, the floor R2, and the walls R3 to R6 of the room R10 and each region of the thermal image (that is, the ceiling R1, the floor R2, and the wall R3). ⁇ Which region of the thermal image R6 is reflected in) is specified.
  • the imaging direction of the thermal image is specified, if the floor plan of the room R10 is specified, the composition of the thermal image (the projection area on the thermal image of the ceiling R1, the floor R2, and the walls R3 to R6) is specified. It is possible.
  • the floor R2 of the room R10 and the walls R4 to R6 are shown in the thermal image, and the ceiling R1 and the wall R3 are not shown. Therefore, from the thermal image, the change temperature of the floor R2 and the four surfaces of the walls R4 to R6 is detected, and the change temperature of the ceiling R1 and the wall R3 is not detected. In this case, it is considered that there is no thermal contribution from the ceiling R1 and the wall R3. In terms of arithmetic processing, the temperature change of the ceiling R1 and the wall R3 may be regarded as zero.
  • the space temperature estimation unit 24 determines the change temperature of each thermal factor of the detected chamber R10, the position of the person detected by the person position detection unit 2, and the CRI coefficient determined by the CRI coefficient determination unit 23 (the above).
  • the temperature at the position X of the person in the room R10 is calculated (estimated) from the CRI coefficient for each thermal factor of the room R10 at the position X of the person.
  • the space temperature estimation unit 24 multiplies the change temperature of each thermal factor in the chamber R10 by the CRI coefficient for each thermal factor at the position X of the person, and heat from each thermal factor. Calculate the degree of thermal influence and add up the calculated degree of thermal influence.
  • the temperature change at the position X of the person is calculated.
  • the space temperature estimation unit 24 calculates (estimates) the temperature at the person's position X by adding the reference temperature at the person's position X and the change temperature at the person's position X.
  • the average radiation temperature calculation unit 3 calculates the average radiation temperature (MRT: Mean Radiant Temperature) at the position X of the person in the room R10 (that is, the position of the person detected by the person position detection unit 2).
  • the average radiant temperature is a value represented by a temperature display of the average value of heat (radiant heat) radiated from each thermal factor (ceiling R1, floor R2 and each wall R3 to R6) of the chamber R10. More specifically, the average radiation temperature calculation unit 3 has the spatial temperature estimation unit 24 based on the layout information and installation position information input to the setting input unit 21 and the thermal image acquired by the thermal image acquisition unit 22.
  • the change temperature of each thermal factor of the chamber R10 is detected in the same manner as in the case of. Then, the average radiation temperature calculation unit 3 calculates the average radiation temperature at the position X of the person in the room 10 due to the radiant heat from each thermal factor based on the detected change temperature of each thermal factor.
  • the warm / cool feeling index calculation unit 4 calculates an index representing the warm / cold feeling of a person (a person detected by the person position detection unit 2) in the room R10.
  • the warm / cool feeling index calculation unit 4 calculates, for example, PMV (predicted average warm / cold feeling index) as the above index.
  • PMV is calculated from the above-mentioned six elements of heat for a person (air temperature, radiation temperature, relative humidity, wind speed, amount of clothing, amount of activity) using a predetermined calculation formula. Since the predetermined calculation formula for calculating PMV is a well-known calculation formula, the description thereof will be omitted.
  • the PMV is, for example, a value in the range of -3 ⁇ PMV ⁇ +3.
  • a predetermined range for example, ⁇ 0.5 or more and +0.5 or less.
  • the calculated PMV is smaller than the lower limit value (for example, ⁇ 0.5) in the predetermined range, it can be estimated that the person's feeling of warmth and coldness is cold.
  • the smaller the calculated PMV is from the lower limit of the predetermined range, the colder the cold level is.
  • the calculated PMV is larger than the upper limit value (+0.5) in the predetermined range, it can be estimated that the person's feeling of warmth and coldness is hot. At this time, the larger the calculated PMV is from the upper limit value (+0.5) in the above-mentioned predetermined range, the hotter the heat level can be estimated.
  • the air temperature is the air temperature at the position X of the above person, and the estimation result of the space temperature estimation unit 24 is used as this air temperature.
  • the radiation temperature is the radiation temperature from each thermal factor of the chamber R10 at the position X of the person, and the calculation result of the average radiation temperature calculation unit 3 is used as this radiation temperature.
  • the relative humidity is the relative humidity in the room R10, and the measurement result of the humidity sensor installed in the room R10 is used as the relative humidity.
  • the humidity sensor may be, for example, a humidity sensor provided in the air conditioner 7.
  • the wind speed is the wind speed at the position X of the person, and the wind speed of the air conditioner 7 may be used as the wind speed.
  • Information on the wind speed of the air blown by the air conditioner 7 can be obtained from the air conditioner 7.
  • a typical wind speed of the wind speed of the air conditioner 7 may be used.
  • the amount of clothing is the amount of clothing of the above-mentioned person, and as the amount of clothing, a representative value of the amount of clothing preset in the warm / cold feeling index calculation unit 4 is used.
  • the amount of activity (that is, the amount of metabolism) is the amount of activity of the above-mentioned person, and as this amount of activity, a representative value of the amount of activity preset in the warm / cold sensation index calculation unit 4 is used.
  • the temperature target value calculation unit 5 calculates the temperature target value of the air conditioner of the air conditioner 7 based on the calculation result of the temperature / cold feeling index calculation unit 4. More specifically, when the PMV calculated by the temperature target value calculation unit 4 is within the above-mentioned predetermined range (for example, ⁇ 0.5 or more and +0.5 or less), the temperature target value calculation unit 5 is an air conditioner. Maintain (reset) the temperature target value of the air conditioner in 7 to the current temperature target value. When the PMV calculated by the thermal sensation index calculation unit 4 is larger than the upper limit value (for example, +0.5) in the predetermined range, the temperature target value calculation unit 5 sets the temperature target value of the air conditioning of the air conditioner 7 as the temperature target value.
  • the predetermined range for example, ⁇ 0.5 or more and +0.5 or less
  • the temperature target value calculation unit 5 sets the temperature target value of the air conditioning of the air conditioner 7 as the temperature target value. The smaller the PMV compared to the current temperature target value, the larger the value is changed (set).
  • the air conditioning control unit 6 controls the air conditioning device 7 based on the temperature target value set by the temperature target value calculation unit 5. More specifically, the air conditioner 7 controls the temperature of the air blown from the air conditioner 7 so that the air temperature in the chamber R10 approaches the set temperature set in the air conditioner 7. The air conditioning control unit 6 controls the temperature inside the room R10 to a temperature that the person in the room R10 feels comfortable by setting the above set temperature to the temperature target value set by the temperature target value calculation unit 5. To do.
  • the air conditioner 7 controls the temperature of the air blown from the air conditioner 7 so that the air temperature in the room R10 approaches the set temperature set in the air conditioner 7. More specifically, the air conditioner 7 has a temperature sensor that detects the air temperature in the chamber R10, and controls the temperature of the air blown from the air conditioner 7 so that the detected temperature of the temperature sensor approaches the set temperature. To do. Further, the air conditioner 7 controls the humidity of the air blown from the air conditioner 7 so that the humidity in the room R10 approaches the set humidity set in the air conditioner 7. More specifically, the air conditioner 7 has a humidity sensor that detects the humidity in the chamber R10, and controls the humidity of the air blown from the air conditioner 7 so that the detection temperature of the humidity sensor approaches the set humidity. To do.
  • the space temperature estimation unit 24 composes the thermal image acquired by the thermal image acquisition unit 22 based on the installation position information (information on the installation position of the air conditioner 7) input to the setting input unit 21.
  • the reflection area of the ceiling R1, the floor R2, and the walls R3 to R6) is specified.
  • the installation position P1 of the air conditioner is the “center” (center of the upper part of the wall R3) PS1
  • the composition of the thermal image G1 is as shown in FIG. 4A.
  • the floor R2 and the wall R5 are specified to be centered in the left-right direction in the thermal image G1.
  • the space temperature estimation unit 24 detects the change temperature of each thermal factor (floor R2 and each wall R4 to R6).
  • the installation position P1 of the air conditioner 7 is the “right end” (right end of the upper part of the wall R3) RR1, as shown in FIG. 4B, the composition of the thermal image G1 is the floor R2.
  • the wall R5 is specified to be located to the left of the center in the left-right direction in the thermal image G1. Then, based on this specified composition, the space temperature estimation unit 24 detects the change temperature of each thermal factor (floor R2 and each wall R4 to R6).
  • the space temperature estimation unit 24 detects the change temperature of each thermal factor (floor R2 and each wall R4 to R6).
  • the floor plan information and the set position information are input to the setting input unit 21 (S1).
  • the CRI coefficient determining unit 23 determines the room R10 based on the floor plan information and the installation position information input to the setting input unit 21 and the information on the blowing speed and the wind direction of the air conditioner 7 acquired from the air conditioner 7. Identify the floor plan conditions.
  • the CRI coefficient determination unit 23 determines a set of CRI coefficients corresponding to the specified floor plan condition from among a plurality of sets of preset CRI coefficients (S2).
  • the thermal image acquisition unit 22 includes heat including any one or more surfaces (for example, four surfaces of the floor R2 and walls R4 to R6) of each surface (ceiling R1, floor R2 and walls R3 to R6) of the chamber R10. Acquire an image (S3). Then, the person position detection unit 2 detects the position X of the person in the room R10 based on the thermal image acquired by the thermal image acquisition unit 22 (S4).
  • the space temperature estimation unit 24 is present in the change temperature of each thermal factor (ceiling R1, floor and R2 and four walls R3 to R4) of the room R10 and the detection result of the person position detection unit 2 (that is, in the room R10).
  • the temperature at the person's position X in the room R10 is estimated (S5).
  • the above-mentioned "change temperature of each thermal factor of the chamber R10" is determined by the space temperature estimation unit 24 in the thermal image acquisition unit 22 based on the above-mentioned layout information and installation position information input to the setting input unit 21.
  • the composition of the acquired thermal image (the projection area of the ceiling R1, the floor R2, and the walls R3 to R6) is specified, and the composition is detected from the specified thermal image.
  • the average radiation temperature calculation unit 3 calculates the average radiation temperature (MRT) at the position (position detected by the person position detection unit 2) X of the person in the room R10 based on the specified thermal image of the composition. Calculate (S6).
  • the temperature / cold sensation index calculation unit 4 averages the temperature of the above-mentioned person's position X estimated by the space temperature estimation unit 24 and the above-mentioned person's position X calculated by the average radiation temperature calculation unit 3. Based on the radiation temperature, the wind speed of the wind blown from the air conditioner 7, the humidity of the room R10, and the amount of clothing and activity of the above person, as a hot / cold feeling index representing the above person's feeling of warmth and coldness. The PMV is calculated (S7). Then, the temperature target value calculation unit 5 calculates the temperature target value of the air conditioning of the air conditioner 7 based on the calculation result of the hot / cold feeling index calculation unit 4 (S8), and the air conditioning control unit 6 calculates the temperature target value. The air conditioner 7 is controlled based on the temperature target value calculated by the unit 5 (S9).
  • the temperature of the air in the indoor space of the chamber R10 is estimated using the CRI coefficient for at least one surface of R3 to R6. Therefore, even if the temperature of each thermal factor (ceiling R1, floor R2 and each wall R3 to R6) changes, it is only necessary to add up the influence (thermal contribution amount) of each thermal factor due to the temperature change (that is,).
  • the temperature of the air in the indoor space of the chamber R10 can be estimated (by a simple calculation).
  • the hot / cold feeling estimation system 100 can be estimated more accurately by using the above space temperature estimation system 1.
  • the space temperature estimation method includes a detection step and a space temperature estimation step.
  • the detection step detects the temperature of the ceiling of the room R10, the floor R2, and at least one surface of each wall R3 to R6.
  • the space temperature estimation step has a CRI coefficient for at least one of the above surfaces, and estimates the temperature of the air in the indoor space of the room R10 using the temperature detected in the detection step and the above CRI coefficient.
  • the hot / cold feeling estimation method is a hot / cold feeling estimation method using the above space temperature estimation method.
  • the hot / cold sensation estimation method includes a human position detection step and a hot / cold sensation estimation step.
  • the person position detection step detects the position X of a person in the room R10.
  • an index representing the warm / cold sensation of a person is calculated using the temperature of the air at the position X of the person estimated by the above space temperature estimation method.
  • the computer program according to one aspect is a program for causing a computer system to execute the above space temperature estimation method.
  • the computer program according to one aspect is a program for causing the computer system to execute the above-mentioned heating / cooling sensation estimation method.
  • this modification includes a floor plan estimation unit 25 instead of the setting input unit 21 in the above embodiment.
  • the floor plan estimation unit 25 estimates the floor plan of the room R10 and the installation position of the air conditioner 7 by using the thermal image acquired by the thermal image acquisition unit 22. That is, in this modification, the floor plan of the room R10 and the installation position of the air conditioner 7 are estimated from the above thermal image by the floor plan estimation unit 25 instead of being input by the user to the setting input unit 21.
  • the thermal image acquisition unit 22 (temperature detecting means) and the floor plan estimation unit 25 constitute the floor plan estimation means 40.
  • the thermal image acquisition unit 22 acquires a thermal image including the floor plan of the chamber R10.
  • the thermal image acquisition unit 22 may be an infrared array sensor.
  • the floor plan estimation unit 25 estimates the floor plan of the room R10 and the installation position of the air conditioner 7 by using the history of the position of the person (person position history) reflected in the thermal image.
  • the person position history is a set of people's positions plotted on a thermal image over a certain period of time.
  • the positions of people in the thermal image are mainly distributed on the region corresponding to the floor R2 in the thermal image, and hardly distributed on the regions corresponding to the walls R3 to R6 in the thermal image.
  • the floor plan estimation unit 25 detects a region (that is, a person) corresponding to a person from the thermal image, and accumulates a history of the position of the person (person position history) in the thermal image over a certain period of time. Then, the floor plan estimation unit 25 estimates the floor plan of the room R10 and the installation position of the air conditioner 7 based on the accumulated person position history.
  • the installation position of the air conditioner 7 is estimated as follows. That is, when the installation position of the air conditioner 7 is the "right end” (right end of the upper part of the wall R3) PR1 (see FIG. 3), when the person position history is superimposed on the thermal image G1, as shown in FIG. 7A, On the right side of the thermal image G1, a person-absent region F1 in which no person exists is generated. Therefore, when the absentee region F1 is generated on the right side of the thermal image G1, the floor plan estimation unit 25 estimates this region F1 as the wall R6, and air-conditions the area F1 from the reflection position of the wall R6 on the thermal image G1. The installation position of the device 7 is estimated to be the "right end".
  • the installation position of the air conditioner 7 is “center” (center of the upper part of the wall R3) PS1 (see FIG. 3)
  • the person position history is superimposed on the thermal image G1, as shown in FIG. 7B
  • Human absent areas F2 and F3 are generated on both the left and right sides of the thermal image G1. Therefore, when the area F2 and the area F3 are generated on the left and right sides of the thermal image G1, the floor plan estimation unit 25 estimates these areas F2 and F3 as the walls R4 and R6, and the wall R4 in the thermal image G1.
  • the installation position of the air conditioner 7 is estimated to be "center” from the reflection position of R6.
  • the installation position of the air conditioner 7 is the "left end" (the left end of the upper part of the wall R3) PL1 (see FIG. 3)
  • the person position history is superimposed on the thermal image G1, as shown in FIG. 7C
  • An absentee region F4 is generated on the left side of the thermal image G1. Therefore, when the absentee region F4 is generated on the left side of the thermal image G1, the floor plan estimation unit 25 estimates this region F4 as the wall R4, and air-conditions the area F4 from the reflection position of the wall R4 on the thermal image G1.
  • the installation position of the device 7 is estimated to be the "left end".
  • the floor plan of room R10 is estimated as follows. That is, in the thermal image G1 in which the person position history is superimposed, the floor plan estimation unit 25 has a floor R2 with a distance H1 between the position XP1 of the person farthest from the lower side of the thermal image G1 in the upward direction and the lower side of the thermal image G1. It is estimated that the depth is W1 (see FIG. 3) (see FIGS. 7A to 7C).
  • the floor area of the chamber R10 is estimated from the air conditioning capacity of the air conditioner 7 as described in the above embodiment.
  • the floor plan estimation unit 25 estimates the width W2 (see FIG. 3) of the room R10 from the estimated floor area and depth W1, and from the estimated width W2 and depth W1, the floor shape (horizontal length) of the floor R2 of the room R10. Estimate a rectangle, portrait rectangle or square).
  • the floor plan estimation unit 25 estimates the floor plan of the room R10 and the installation position of the air conditioner 7 from the thermal image, so that the user sets the floor plan information of the room R10 and the installation position information of the air conditioner 7. No need to enter.
  • the space temperature estimation unit 24 changes the CRI coefficient used in the space temperature estimation according to the floor plan of the room R10 estimated by the floor plan estimation unit 25 and the installation position of the air conditioner 7.
  • Modification 2 In the first modification, the floor plan of the room R10 and the process of estimating the installation position of the air conditioner 7 are all performed in the floor plan estimation unit 25 (that is, in the space temperature estimation system 1). However, as shown in FIG. 8, the thermal image G1 is transmitted from the floor plan estimation unit 25 to the external device 27 via the communication network CR1 (cloud), and the thermal image G1 (FIG. 9A) is restored by the external device 27 (high). The restored image G2 (FIG. 9B) may be generated by performing resolution). Then, the restored image G2 may be returned from the external device 27 to the floor plan estimation unit 25, and the floor plan estimation unit 25 may estimate the floor plan of the room R10 and the installation position of the air conditioner 7 based on the restored image G2.
  • the floor plan estimation unit 25 estimates, for example, the boundary line between the floor R2 and each wall R4 to R6 based on the contour of the floor R2 of the restored image G2, and the composition of the restored image G2 (floor R2 and each wall R4 to R6). (Implementation area) is specified. Then, the floor plan estimation unit 25 estimates the installation position (right end, center or left end) of the air conditioner 7 and the floor plan (vertical rectangle, horizontally long rectangle or square) of the room R10 from the specified composition.
  • the floor plan (floor shape) of the room R10 can be estimated as follows. That is, the floor plan estimation unit 25 obtains the depth W1 (see FIG. 3) of the room R10 from the specified composition, and obtains the width W2 (see FIG. 3) of the floor R2 from the floor area and depth W1 of the floor R2. The floor shape of the floor R2 is obtained from the depth W1 and the width W2. Further, the installation position of the air conditioner 7 can be estimated as follows. That is, the floor plan estimation unit 25 estimates the installation position of the air conditioner 7 according to the position of the wall R5 (that is, the wall in front of the air conditioner 7) in the left-right direction in the specified composition, for example.
  • the thermal image G1 can be restored using the external device 27, and the floor plan estimation unit 25 can estimate the floor plan of the room R10 and the installation position of the air conditioner 7 using the restored image G2. Therefore, the floor plan of the room R10 and the installation position of the air conditioner 7 can be estimated more accurately.
  • the spatial temperature estimation system 1 further includes a scanning device 28 that scans the imaging direction (detection direction) of the thermal image acquisition unit 22 (temperature detecting means) in the above embodiment. Be prepared.
  • the scanning device 28 has a rotation shaft portion and a rotation drive source (for example, an electric motor).
  • the rotation shaft portion can be rotated by a rotation drive source.
  • a thermal image acquisition unit 22 is fixed to the rotating shaft unit.
  • the thermal image acquisition unit 22 rotates around the rotation shaft unit as the rotation shaft unit rotates, and the imaging direction (detection direction) of the thermal image acquisition unit 22 is scanned in, for example, the left-right direction of the chamber R10.
  • the rotation drive source scans the imaging direction of the thermal image acquisition unit 22 in the left-right direction by alternately rotating the rotation shaft portion in the forward rotation direction and the reverse rotation direction.
  • the imaging range (detection range) of the thermal image acquisition unit 22 can be expanded to a wider range.
  • this modification includes a plurality of thermal image acquisition units 22 (thermal image acquisition means) and an image composition unit 31 in the modification 1.
  • the plurality of thermal image acquisition units 22 image different regions in the chamber R10 (more specifically, regions having different floor plans of the chamber R10).
  • the image composition unit 31 combines the thermal images acquired by each of the plurality of thermal image acquisition units 22 to generate a single composite thermal image.
  • the floor plan estimation unit 25 of this modification estimates the floor plan of the room R10 and the installation position of the air conditioner 7 by using the synthetic thermal image generated by the image synthesis unit 31.
  • the plurality of thermal image acquisition units 22, the image composition unit 31, and the floor plan estimation unit 25 constitute the floor plan estimation means 40.
  • the person position detecting unit 2 detects the position X of the person in the room R10 by using the synthetic thermal image generated by the image synthesizing unit 31.
  • the space temperature estimation unit 24 detects the change temperature of each thermal factor in the chamber R10 by using the synthetic thermal image generated by the image synthesis unit 31.
  • the average radiation temperature calculation unit 3 uses the synthetic thermal image generated by the image synthesis unit 31 to radiate the average radiation at the position X of the person in the room R10 (that is, the position of the person detected by the person position detection unit 2). Calculate the temperature.
  • this modified example further includes an imaging device 32 that acquires an captured image including the floor plan of the chamber R10 in the modified example 1.
  • the image pickup apparatus 32 is installed in or around the air conditioner 7.
  • the floor plan estimation unit 25 of this modification uses the image captured by the image pickup device 32 instead of the thermal image acquired by the thermal image acquisition unit 22, and determines the floor plan of the room R10 and the installation position of the air conditioner 7. presume.
  • the floor plan estimation means 40 includes a floor plan estimation unit 25 and an image pickup device 32.
  • the image pickup device 32 includes an image pickup element such as a CCD and an image pickup lens.
  • the floor plan estimation unit 25 binarizes the captured image of the image pickup device 32 to extract the boundary lines of the ceiling R1, the floor R2, and the walls R3 to R6 of the chamber R10. Then, the floor plan estimation unit 25 specifies the composition of the captured image (the projection area of the ceiling R1, the floor R2, and the walls R3 to R6) based on the extracted boundary line, and the room R10 is based on the specified composition.
  • the floor plan and the installation position of the air conditioner 7 are estimated.
  • this modification includes a TOF (Time Of Flight) camera 33 that acquires a distance distribution image including the floor plan of the chamber R10 instead of the image pickup device 32 in the modification 4.
  • the TOD camera 33 is installed in or around the air conditioner 7.
  • the FOT camera 33 measures the distance to the subject for each pixel according to the reflection time of the pulsed light (for example, near-infrared light) from the subject, and colors each pixel according to the distance to represent the image (separation). It is a camera that generates a distribution image).
  • the floor plan estimation means 40 includes a TOF camera 33 and a floor plan estimation unit 25.
  • the floor plan estimation unit 25 extracts the boundary lines of the ceiling R1, the floor R2, and the walls R3 to R6 of the room R10 from the image (distance distribution image) captured by the TOF camera 33. Then, the floor plan estimation unit 25 specifies the composition of the captured image (the projection area of the ceiling R1, the floor R2, and the walls R3 to R6) based on the extracted boundary line, and the room R10 is based on the specified composition. The floor plan and the installation position of the air conditioner 7 are estimated.
  • the space temperature estimation system (1) includes a temperature detecting means (22) and a space temperature estimating means (24).
  • the temperature detecting means (22) detects the temperature of at least one surface of the ceiling (R1), the floor (R2), and the plurality of walls (R3 to R6) of the room (R10).
  • the space temperature estimating means (24) estimates the temperature of the air in the indoor space of the room (R10) using the CRI coefficient for at least one surface and the temperature detected by the temperature detecting means (22).
  • the temperature of at least one surface of the ceiling (R1), floor (R2) and each wall (R3 to R6) of the room (R10) and the CRI coefficient for at least one surface of the room are used.
  • the temperature of the air in the indoor space of (R10) is estimated. Therefore, even if the temperature of each thermal factor (ceiling (R1), floor (R2) and each wall (R3 to R6)) changes, it is only necessary to add up the effects of the temperature change of each thermal factor (that is,).
  • the temperature of the air in the indoor space of the room (R10) can be estimated (by a simple calculation).
  • the temperature detecting means (22) is among the ceiling (R1), the floor (R2), and the plurality of walls (R4 to R6). The temperature of two or more surfaces facing different directions is detected.
  • the temperature detecting means (22) is among the ceiling (R1), the floor (R2), and the plurality of walls (R4 to R6). The temperature of two or more surfaces facing different directions is detected.
  • the accuracy of estimating the temperature of the air in the indoor space of the room (R10) can be further improved.
  • the space temperature estimation means (24) has a plurality of chambers (R10) having different floor plans and one. Store a plurality of CRI coefficients corresponding to one-to-one.
  • the optimum CRI coefficient can be used according to the floor plan of the room (R10).
  • the accuracy of estimating the temperature of the air in the indoor space of the room (R10) can be further improved.
  • the temperature detecting means (22) is an infrared array sensor.
  • an infrared array sensor can be used as the temperature detecting means (22).
  • the space temperature estimation system (1) includes a scanning device (28) that scans the detection direction of the temperature detecting means (22) in the fifth aspect.
  • the detection direction of the temperature detecting means (22) can be scanned.
  • the detection range of the temperature detecting means (22) can be expanded to a wider range.
  • the space temperature estimation system (1) according to the seventh aspect further includes a floor plan estimation means (40) for estimating the floor plan of the room (R10) in the fourth aspect.
  • the space temperature estimation means (24) changes the CRI coefficient according to the estimation result of the floor plan estimation means (40).
  • the CRI coefficient is changed according to the floor plan of the room (R10), the temperature of the air in the indoor space of the room (R10) can be estimated more accurately. In addition, it is possible to save the trouble of manually setting and inputting the floor plan information of the room (R10).
  • the floor plan estimation means (40) includes an image pickup device (32) that acquires an image taken including the floor plan of the room (R10).
  • the floor plan of the room (R10) can be estimated using the image captured by the image pickup device (32).
  • the floor plan estimation means (40) includes a TOF camera (33) that acquires a distance distribution image including the floor plan of the room (R10). ..
  • the floor plan of the room (R10) can be estimated using the distance distribution image of the TOF camera (33).
  • the floor plan estimation means (40) includes a thermal image acquisition means (22) for acquiring a thermal image including the floor plan of the room (R10). Including.
  • the floor plan of the room (R10) can be estimated using the thermal image.
  • the floor plan estimation means (40) includes a plurality of thermal image acquisition means (22).
  • the floor plan estimation means (40) synthesizes the thermal images acquired by each of the plurality of thermal image acquisition means (22), and estimates the floor plan of the chamber (R10) based on the image-combined thermal images.
  • thermo image acquisition means (22) According to this configuration, a wide range of thermal images can be acquired by using a plurality of thermal image acquisition means (22).
  • the space temperature estimation system (1) according to any one of the first to eleventh aspects, the person position detecting means (2), and the temperature. It is provided with a cooling sensation index calculation means (4).
  • the person position detecting means (2) detects the position (X) of a person in the room (R10).
  • the heating / cooling sensation index calculating means (4) uses the temperature of the air at the position (X) of the person estimated by the space temperature estimation system (1) to generate an index (for example, PMV) representing the feeling of warmth / coldness of the person. calculate.
  • PMV index
  • the feeling of warmth and coldness of a person in the room (R10) can be estimated by using the above-mentioned space temperature estimation system (1).
  • the heating / cooling sensation estimation system further includes an average radiation temperature calculating means (3) for calculating the average radiation temperature at the position (X) of a person in the room (R10). ..
  • the hot / cold sensation index calculating means (4) further uses the average radiant temperature acquired by the average radiant temperature calculating means (3) to calculate an index (for example, PMV) representing a person's warm / cold sensation.
  • the average radiation temperature at the position of the person in the room (R10) is further used to calculate an index representing the feeling of warmth and coldness of the person in the room (R10), so that the person is in the room (R10).
  • the feeling of warmth and coldness of a person can be estimated more accurately.
  • the space temperature estimation method includes a temperature detection step and a space temperature estimation step.
  • the temperature detection step detects the temperature of at least one surface of the ceiling (R1), the floor (R2), and the plurality of walls (R3 to R6) of the room (R10).
  • the space temperature estimation step estimates the temperature of the air in the indoor space of the room (R10) using the CRI coefficient for at least one surface and the temperature detected by the temperature detecting means (22).
  • the temperature of at least one surface of the ceiling (R1), floor (R2) and each wall (R3 to R6) of the room (R10) and the CRI coefficient for at least one surface of the room are used.
  • the temperature of the air in the indoor space of (R10) is estimated. Therefore, even if the temperature of each thermal factor (ceiling (R1), floor (R2) and each wall (R3 to R6)) changes, it is only necessary to add up the effects of the temperature change of each thermal factor (that is,).
  • the temperature of the air in the indoor space of the room (R10) can be estimated (by a simple calculation).
  • the hot / cold feeling estimation method is a hot / cold feeling estimation method using the space temperature estimation method according to the fourteenth aspect.
  • the hot / cold feeling estimation method includes a human position detection step and a hot / cold feeling index calculation step.
  • the person position detection step detects the position (X) of a person in the room (R10).
  • an index representing the warm / cold sensation of a person is calculated using the temperature of the air at the position (X) of the person estimated by the space temperature estimation method.
  • the program according to the sixteenth aspect is a program for causing a computer system to execute the space temperature estimation method according to the fourteenth aspect.
  • the program according to the seventeenth aspect is a program for causing a computer system to execute the hot / cold feeling estimation method according to the fifteenth aspect.
  • Spatial temperature estimation system 3 Average radiation temperature calculation unit (average radiation temperature calculation means) 4 Warm / cold sensation index calculation unit (warm / cold sensation index calculation means) 22 Thermal image acquisition unit (temperature detection means, thermal image acquisition means) 24 Spatial temperature estimation unit (spatial temperature estimation means) 29 Scanning device 32 Imaging device 33 TOF camera 40 Floor plan estimation means 100 Temperature / cold feeling estimation system R1 Ceiling R2 Floor R3 to R6 Wall R10 Room X Position of people

Abstract

簡単な計算で、室内空間の空気の温度を推定できる空間温度推定システムを提供する。空間温度推定システム（１）は、熱画像取得部（２２）と、空間温度推定部（２４）と、を備える。熱画像取得部（２２）は、室の天井、床及び複数の壁の少なくとも１面の温度を検出する。空間温度推定部（２４）は、少なくとも１面に関するＣＲＩ係数、及び熱画像取得部（２２）で検出された温度を用いて、室（Ｒ１０）の室内空間の空気の温度を推定する。

Description

空間温度推定システム、温冷感推定システム、空間温度推定方法、温冷感推定方法、及び、プログラム
　本開示は、空間温度推定システム、温冷感推定システム、空間温度推定方法、温冷感推定方法、及び、プログラムに関する。より詳細には、本開示は、ＣＲＩ（Contribution Ratio of Indoor Climate：温熱環境形成寄与率）係数を用いる空間温度推定システム、この空間温度推定システムを用いた温冷感推定システム、空間温度推定方法、温冷感推定方法、及び、プログラムに関する。
　特許文献１に記載の空気調和機は、熱画像取得部（サーモグラフィ）、演算部、及び制御部を備える。熱画像取得部は、空間の温度分布を表す熱画像を取得する。演算部は、熱画像取得部が取得した熱画像内における人に該当する領域を特定する。そして、演算部は、熱画像内における人に該当する領域の温度分布に基づいて空間に居る人の温度である人体温度を定め、人体温度と、人に該当する領域以外の領域の温度から得られる周囲温度との差分値に基づいて空間に居る人の温冷感を推定する。制御部は、演算部が推定した温冷感に基づいて、空気調和機の風量、風温、風向のうち少なくとも１つを制御する。
　特許文献１記載の空気調和機では、熱画像取得部（サーモグラフィ）によって取得された熱画像において、人に該当する領域以外の領域から得られる温度を、人の周囲の空間の温度とする。しかし、熱画像は、熱画像に映った熱的要因（壁、天井及び床など）の温度を表し、熱画像に映った空間の空気の温度を必ずしも表していない。このため、特許文献１に記載の空気調和機では、人の周囲の空間の空気の温度を適切に求められていない場合がある。
国際公開第２０１５／１２２２０１号
　本開示は、簡単な計算で、室内空間の空気の温度を推定できる空間温度推定システム、温冷感推定システム、空間温度推定方法、温冷感推定方法、及び、プログラムを提供することを目的とする。
　本開示の一態様の空間温度推定システムは、温度検出手段と、空間温度推定手段と、を備える。前記温度検出手段は、室の天井、床及び複数の壁の少なくとも１面の温度を検出する。前記空間温度推定手段は、前記少なくとも１面に関するＣＲＩ係数、及び前記温度検出手段で検出された前記温度を用いて、前記室の室内空間の空気の温度を推定する。
　本開示の一態様の温冷感推定システムは、前記空間温度推定システムと、人位置検出手段と、温冷感指標算出手段と、を備える。前記人位置検出手段は、前記室に居る人の位置を検出する。前記温冷感指標算出部は、前記空間温度推定システムで推定された前記人の位置での空気の温度を用いて、前記人の温冷感を表す指標を算出する。
　本開示の一態様の空間温度推定方法は、温度検出工程と、空間温度推定工程と、を備える。前記温度検出工程は、室の天井、床及び複数の壁の少なくとも１面の温度を検出する。前記空間温度推定工程は、前記少なくとも１面に関するＣＲＩ係数、及び前記温度検出手段で検出された前記温度を用いて、前記室の室内空間の空気の温度を推定する。
　本開示の一態様の温冷感推定方法は、前記空間温度推定方法を用いた温冷感推定方法である。前記温冷感推定方法は、人位置検出工程と、温冷感指標算出工程と、を備える。前記人位置検出工程は、前記室内に居る人の位置を検出する。前記温冷感指標算出工程は、前記空間温度推定方法によって推定された前記人の位置での空気の温度を用いて、前記人の温冷感を表す指標を算出する。
　本開示の一態様のプログラムは、前記空間温度推定方法をコンピュータシステムに実行させるためのプログラムである。
　本開示の一態様のプログラムは、前記温冷感推定方法をコンピュータシステムに実行させるためのプログラムである。
図１は、実施形態に係る温冷感推定システムのブロック図である。 図２は、同上の温冷感推定システムに備えられた空間温度推定システムの空間温度推定方法を説明するための説明図である。 図３は、室内での空調装置の設置位置を説明する説明図である。 図４Ａ～図４Ｃはそれぞれ、同上の空調装置の設置位置が「中央」「左端」「右端」である場合の熱画像の構図を説明する説明図である。 図５は、同上の温冷感推定システムの動作を説明するフローチャートである。 図６は、変形例１に係る温冷感推定システムのブロック図である。 図７Ａ～図７Ｃはそれぞれ、変形例１に係る温冷感推定システムのブロック図である。 図８は、変形例２に係る温冷感推定システムのブロック図である。 図９Ａは、熱画像の一例を説明する説明図である。図９Ｂは、復元画像の一例を説明する説明図である。 図１０は、変形例３に係る温冷感推定システムのブロック図である。 図１１は、変形例４に係る温冷感推定システムのブロック図である。 図１２は、変形例５に係る温冷感推定システムのブロック図である。 図１３は、変形例６に係る温冷感推定システムのブロック図である。
　（実施形態）
　図１～図５を参照して、本実施形態に係る空間温度推定システム１及び温冷感推定システム１００について説明する。
　空間温度推定システム１（図１参照）は、対象となる室の室内空間の温度分布を推定可能なシステムである。室内空間の温度分布とは、室内空間の各位置での空気の温度分布である。温冷感推定システム１００（図１参照）は、空間温度推定システム１によって推定された室内空間の温度分布を用いて、対象となる室内に居る人の温冷感を推定可能なシステムである。本実施形態では、温冷感推定システム１００は、推定した温冷感に基づいて、対象となる室内の温度を制御する空調装置を制御する。これにより、対象となる室内の温度が、室内に居る人が快適に感じる温度に自動的に制御される。
　なお、本実施形態では、空間温度推定システム１及び温冷感推定システム１００は、空調装置７の外部に配置されるが、空調装置７の内部に配置されてもよい。対象となる室は、建物内の室でもよいし、自動車（例えばバス）などの移動体内の室でもよい。
　空間温度推定システム１は、上述のように、対象となる室の室内空間の温度分布を推定可能である。より詳細には、空間温度推定システム１は、ＣＲＩ（Contribution Ratio of Indoor Climate：温熱環境形成寄与率）係数を用いて室内空間の各位置の空気の温度（気温）を推定可能である。
　ここで、本実施形態での室内空間の空気の温度（気温）の推定の考え方を説明する。室の間取りと室内の空気の流れとが決まると、吹出口（例えば空調装置）から吹き出された空気の流れが室内の或る位置に居る人に及ぼす影響は一定であると見なせる。このため、室内空間の或る位置の気温は、決められた室の間取り及び空気の流れの下での室の壁床温度といった熱的要因からの熱的影響度を足し合わせることで計算可能である。このような熱的影響度をＣＲＩ係数を用いて算出することで、室内空間の或る位置の気温を簡易に算出可能（推定可能）である。なお、熱的要因とは、熱源となる要因であり、例えば室の天井、床及び壁などである。
　すなわち、最初にＣＦＤ（Computational Fluid Dynamics：数値流体力学）解析を用いて室内の気流分布を決定し、その気流分布に基づいて各熱的要因から室内の或る位置に寄与する熱的影響度を算出する。これにより、各熱的要因の温度が変化した場合でも、各熱的要因の温度変化による熱的影響度を足し合わせるだけで、或る位置での気温を算出可能である。
　具体的には、図２に示すように、空調装置７が設置された室Ｒ１０を想定する。この場合、室Ｒ１０内の或る位置Ｘでの気温θｘは、式１及び式２で計算可能である。なお、図２の例では、室Ｒ１０の天井Ｒ１、床Ｒ２及び周囲（例えば４つ）の壁Ｒ３～Ｒ６の計６つの面が熱的要因である。
　θｘ＝θｂ＋Δθｘ ・・・式１
　Δθｘ＝Ｃｘ,１・θ０,１＋Ｃｘ,２・θ０,２＋・・・＋Ｃｘ,６・θ０,６ ・・・式２
　ここで、θｂは、位置Ｘでの基準温度であり、本実施形態では例えば０（ゼロ）度に設定されている。Δθｘは、位置Ｘでの変化温度Δθｘである。すなわち、位置Ｘでの空気の温度θｘは、位置Ｘでの基準温度θｂと、位置Ｘでの変化温度Δθｘとの和で算出される。
　位置Ｘでの変化温度Δθｘは、各熱的要因からの熱的影響による温度である。位置Ｘでの変化温度Δθｘは、式２で与えられる。θ０,１～θ０,６はそれぞれ、各熱的要因（室Ｒ１０の天井Ｒ１、床Ｒ２及び４つの壁Ｒ３～Ｒ６）の変化温度である。Ｃｘ,１～Ｃｘ,６はそれぞれ、位置Ｘに対する各熱的要因に関するＣＲＩ係数であり、各熱的要因から位置Ｘへの変化温度の熱的寄与率を表す。なお、各熱的要因に関するＣＲＩ係数は、室Ｒ１０の間取りと室Ｒ１０内の気流分布が決まると、決まる値である。熱的要因の変化温度と、熱的要因に関するＣＲＩ係数との積によって、熱的要因から位置Ｘに寄与する熱的影響度が算出される。すなわち、位置Ｘでの変化温度Δθｘは、各熱的要因から位置Ｘに寄与する熱的影響度の和で算出される。そして、各熱的要因から位置Ｘに寄与する熱的影響度はそれぞれ、熱的要因の変化温度θ０,１～θ０,６と、熱的要因に関するＣＲＩ係数Ｃｘ,１～Ｃｘ,６との積で算出される。
　室Ｒ１０の間取りと空調装置７の設置位置とから、室Ｒ１０内の気流分布が決まる。そして、室Ｒ１０の間取りと室Ｒ１０内の気流分布（空調装置７の送風の気流分布）とから、各熱的要因に関するＣＲＩ係数が決まる。よって、各熱的要因の変化温度θ０,１～θ０,６を測定することで、測定した各熱的要因の変化温度θ０,１～θ０,６と、各熱的要因に関するＣＲＩ係数とから、或る位置Ｘでの変化温度Δθｘが決まる。この変化温度Δθｘと基準温度θｂとから、室内の或る位置Ｘでの温度が求まる。
　次に図１を参照して、空間温度推定システム１を含む温冷感推定システム１００の構成について詳しく説明する。以下の説明では、図２に示す室Ｒ１０の室内空間の温度を推定する場合を想定する。
　温冷感推定システム１００は、空間温度推定システム１と、人位置検出部２と、平均放射温度算出部３と、温冷感指標算出部４（温冷感指標算出手段）と、温度目標値算出部５と、空調制御部６と、空調装置７とを備えている。空間温度推定システム１は、設定入力部２１と、熱画像取得部２２（温度検出手段、熱画像取得手段）と、ＣＲＩ係数決定部２３と、空間温度推定部２４（空間温度推定手段）とを備えている。なお、本実施形態では、空調装置７は、温冷感推定システム１００の構成に含まれるが、温冷感推定システム１００の構成に含まれなくてもよい。
　設定入力部２１は、操作者の操作入力によって、各種の設定情報（間取り情報及び設置位置情報）の入力を受け付ける。間取り情報は、対象となる室Ｒ１０の間取りに関する情報である。間取り情報は、例えば室Ｒ１０の床面形状（例えば縦長矩形、横長矩形又は正方形）及び床面積に関する情報である。なお、空調装置７の能力は、室Ｒ１０の床面積に応じて決めるため、本実施形態では、空調装置７の空調能力に関する情報（定格消費電力など）を空調装置７から取得して自動的に設定される。設定入力部２１は、室Ｒ１０の壁に設置されていてもよく、リモコン装置として構成されていてもよい。また、空調装置７のリモコン装置と兼用されていてもよい。
　設置位置情報は、室Ｒ１０内での空調装置７の設置位置に関する情報である。設置位置情報は、空調装置７の設置位置の高さ（床Ｒ２からの高さ）、及び空調装置７が設置された壁Ｒ３の左右方向の位置（左端、中央又は右端）に関する情報である。空調装置７の設置位置の高さは、室Ｒ１０の高さとおおよそ同じであり、本実施形態では、固定値（例えば２．２ｍ）に予め設定されている。以下、室Ｒ１０の間取り、空調装置７の設置位置、及び空調装置７の送風の速度及び風向を、まとめて「間取り条件」と記載する。
　なお、操作者が設定入力部２１から間取り情報及び設置情報を設定入力する代わりに、空間温度推定システム１が、例えば、熱画像取得部２２で取得された熱画像を用いて自動的に間取り情報及び設置情報を推定してもよい。この変形例については、後述する。
　ＣＲＩ係数決定部２３は、設定入力部２１に入力された間取り情報及び設置位置情報と、空調装置７の送風の速度及び風向の情報とに基づいて、室Ｒ１０の間取り条件を特定し、特定した間取り条件に対応したＣＲＩ係数を決定する。
　より詳細には、室Ｒ１０の間取り条件は、設定入力部２１に入力された間取り情報及び設置位置情報と、空調装置７の送風の速度及び風向とで特定される。間取り情報には、室Ｒ１０の床面形状（縦長矩形、横長矩形又は正方形）の情報が含まれ、設置位置情報には、空調装置７の設置位置（中央、左端又は右端）の情報が含まれている。空調装置７の送風の速度及び風向の情報は、空調装置７から取得可能である。本実施形態では、空調装置７の送風の速度及び風向の情報は、代表的な風速及び風向をＣＲＩ係数決定部２３に予め設定されており、その設定値が用いられる。室Ｒ１０の取り得る間取りの数は、例えば、横長矩形、縦長矩形及び正方形の計３つであり、空調装置７の取り得る設置位置の数は、例えば、左端、右端及び中央の計３つである。空調装置７の送風の取り得る速度は、例えば予め設定された設置値の１つであり、空調装置７の送風の取り得る風向は、例えば予め設定された設定値の１つである。よって、本実施形態では、室Ｒ１０の取り得る間取り条件の数は、例えば９個（＝３×３×１×１）である。ＣＲＩ係数は、間取り条件ごとに決まるため、本実施形態では、例えば９組のＣＲＩ係数が存在する。
　ＣＲＩ係数決定部２３は、複数組（本実施形態では９組）のＣＲＩ係数を所定の記憶部に保存している。複数組のＣＲＩ係数は、室Ｒ１０の取り得る複数の間取り条件に１対１に対応している。そして、ＣＲＩ係数決定部２３は、設定入力部２１に入力された間取り情報及び設置位置情報と、空調装置７の送風の速度及び風向の情報とに基づいて室Ｒ１０の間取り条件を特定する。そして、ＣＲＩ係数決定部２３は、複数組のＣＲＩ係数の中から、特定した間取り条件に対応したＣＲＩ係数の組を選択（決定）する。
　すなわち、本実施形態では、特定された間取り条件に基づいてＣＲＩ係数が算出される代わりに、予め用意された複数組のＣＲＩ係数の中から、特定された間取り条件に対応したＣＲＩ係数の組が選択（決定）される。これにより、室Ｒ１０の空間温度推定の処理負担を軽減できると共にＣＲＩ係数をより速く決定できる。また、間取り条件に応じてＣＲＩ係数が変更されるため、より精度よく、室Ｒ１０の室内空間の空気の温度を推定できる。
　熱画像取得部２２は、室Ｒ１０の各面（天井Ｒ１、床Ｒ２及び壁Ｒ３～Ｒ６）の何れか１面以上の面を含む熱画像を取得（すなわち撮像）する。熱画像取得部２２は、例えば、空調装置７に設置されるか、又は空調装置７の周囲に設置されている。熱画像は、例えば、画像に映る被写体（熱的要因）の温度を画素毎に測定して、各画素を温度に応じて色分けして表した温度分布画像である。この熱画像から、各熱的要因（室Ｒ１０の天井Ｒ１、床Ｒ２及び壁Ｒ３～Ｒ６）の変化温度を測定可能である。本実施形態では、例えば、４面（床Ｒ２と、空調装置７が設置された壁Ｒ３以外の３つの壁Ｒ４～Ｒ６）を含む画像を取得する。熱画像取得部２２は、例えば、赤外線アレイセンサであってもよい。赤外線アレイセンサは、複数の赤外線受光素子を縦横に配置させたセンサである。熱画像取得部２２は、室Ｒ１０の各面（天井Ｒ１、床Ｒ２及び壁Ｒ３～Ｒ６）の何れか１面以上の面の温度を検出する温度検出手段である。
　なお、室Ｒ１０の各面（天井Ｒ１、床Ｒ２及び壁Ｒ３～Ｒ６）は、互いに異なる向きを向いている。本実施形態では、熱画像取得部２２、互いに異なる向きを向いた４面を含む熱画像を取得するが、互いに異なる向きを向いた３面以上、又は２面以上を含む熱画像を取得してもよい。なお、熱画像が含む面が多いほど室Ｒ１０の室内空間の温度分布の推定精度をより向上させることができる。
　人位置検出部２は、熱画像取得部２２で取得された熱画像に基づいて、室Ｒ１０内に居る人の位置Ｘを検出する。より詳細には、人位置検出部２は、人に該当する温度及びその温度分布の形状などに基づいて、熱画像において人に該当する領域を特定し、特定した領域に基づいて、室Ｒ１０に居る人の位置Ｘを検出する。例えば、特定した領域内の所定位置が人の位置Ｘとして用いられる。
　空間温度推定部２４は、室Ｒ１０の室内空間の空気の温度（気温）を推定する。より詳細には、空間温度推定部２４は、室Ｒ１０の各熱的要因（天井Ｒ１、床及Ｒ２び４つの壁Ｒ３～Ｒ６）の変化温度と、人位置検出部２の検出結果（すなわち室Ｒ１０に居る人の位置）と、ＣＲＩ係数決定部２３で決定されたＣＲＩ係数とを用いて、室Ｒ１０の室内空間の或る位置Ｘの気温として、室Ｒ１０に居る人の位置での気温を推定する。
　更に詳細には、空間温度推定部２４は、熱画像取得部２２で取得された熱画像を用いて、室Ｒ１０の各熱的要因の変化温度を検出する。具体的には、空間温度推定部２４は、間取り情報に基づいて室Ｒ１０の間取り（床面積及び床面形状）を特定する。そして、空間温度推定部２４は、特定した間取りに基づいて、室Ｒ１０の天井Ｒ１、床Ｒ２及び各壁Ｒ３～Ｒ６と熱画像の各領域との対応関係（すなわち天井Ｒ１、床Ｒ２及び壁Ｒ３～Ｒ６が熱画像のどの領域に映っているか）を特定する。すなわち、熱画像の撮像方向は特定されているため、室Ｒ１０の間取りが特定されれば、熱画像の構図（天井Ｒ１、床Ｒ２及び壁Ｒ３～Ｒ６の熱画像上の映込領域）が特定可能である。
　本実施形態では、熱画像には、室Ｒ１０の床Ｒ２及び各壁Ｒ４～Ｒ６が映り、天井Ｒ１及び壁Ｒ３は映らない。このため、熱画像からは、床Ｒ２及び各壁Ｒ４～Ｒ６の４面の変化温度が検出され、天井Ｒ１及び壁Ｒ３の変化温度は検出されない。この場合は、天井Ｒ１及び壁Ｒ３からの熱的寄与量は無いと見なされる。なお、演算処理上は、天井Ｒ１及び壁Ｒ３の変化温度はゼロと見なされてもよい。
　そして、空間温度推定部２４は、検出した室Ｒ１０の各熱的要因の変化温度と、人位置検出部２で検出された人の位置と、ＣＲＩ係数決定部２３で決定されたＣＲＩ係数（上記人の位置Ｘにおける室Ｒ１０の各熱的要因に関するＣＲＩ係数）とから、室Ｒ１０に居る人の位置Ｘの気温を算出（推定）する。具体的には、空間温度推定部２４は、室Ｒ１０の各熱的要因の変化温度と、上記の人の位置Ｘにおける各熱的要因に関するＣＲＩ係数とを掛け合わせて各熱的要因からの熱的影響度を算出し、算出した各熱的影響度を足し合わせる。この結果、上記の人の位置Ｘでの変化温度が算出される。そして、空間温度推定部２４は、人の位置Ｘでの基準温度と人の位置Ｘでの変化温度とを足し合わせて、人の位置Ｘでの気温を算出（推定）する。
　平均放射温度算出部３は、室Ｒ１０に居る人の位置（すなわち人位置検出部２で検出された人の位置）Ｘでの平均放射温度（ＭＲＴ：Mean Radiant Temperature）を算出する。平均放射温度とは、室Ｒ１０の各熱的要因（天井Ｒ１、床Ｒ２及び各壁Ｒ３～Ｒ６）から放射される熱（放射熱）の平均値を温度表示で表した値である。より詳細には、平均放射温度算出部３は、設定入力部２１に入力された間取り情報及び設置位置情報と、熱画像取得部２２で取得された熱画像とに基づいて、空間温度推定部２４の場合と同様に、室Ｒ１０の各熱的要因の変化温度を検出する。そして、平均放射温度算出部３は、検出した各熱的要因の変化温度に基づいて、各熱的要因からの放射熱による室１０に居る人の位置Ｘでの平均放射温度を算出する。
　温冷感指標算出部４は、室Ｒ１０に居る人（人位置検出部２で検出された人）の温冷感を表す指標を算出する。本実施形態では、温冷感指標算出部４は、上記の指標として、例えばＰＭＶ（予測平均温冷感指標）を算出する。ＰＭＶは、上記の人に対する温熱６要素（気温、放射温度、相対湿度、風速、着衣量、活動量）から所定の演算式を用いて算出される。なお、ＰＭＶを算出するための所定の演算式は、周知の演算式であるため、その説明は省略する。
　なお、ＰＭＶは、例えば－３≦ＰＭＶ≦＋３の範囲の値である。このＰＭＶを用いた場合、算出されたＰＭＶが所定範囲（例えば－０．５以上且つ＋０．５以下）のとき、上記の人の温冷感は快適であると推定可能である。また、算出されたＰＭＶが上記所定範囲の下限値（例えば－０．５）よりも小さいとき、上記の人の温冷感は寒いと推定可能である。このとき、算出されたＰＭＶが上記所定範囲の下限値から小さいほど、寒さのレベルはより寒いと推定可能である。また、算出されたＰＭＶが上記所定範囲の上限値（＋０．５）よりも大きいとき、上記の人の温冷感は暑いと推定可能である。このとき、算出されたＰＭＶが上記所定範囲の上限値（＋０．５）から大きいほど、暑さのレベルはより暑く推定可能である。
　上記の温熱６要素のうち、気温は、上記の人の位置Ｘでの気温であり、この気温として空間温度推定部２４の推定結果が用いられる。放射温度は、上記の人の位置Ｘでの室Ｒ１０の各熱的要因からの放射温度であり、この放射温度として平均放射温度算出部３の算出結果が用いられる。相対湿度は、室Ｒ１０内の相対湿度であり、この相対湿度として、室Ｒ１０に設置された湿度センサの測定結果が用いられる。この湿度センサは、例えば空調装置７に備えられた湿度センサであってもよい。風速は、上記の人の位置Ｘでの風速であり、この風速として、空調装置７の送風の風速が用いられてもよい。空調装置７の送風の風速の情報は、空調装置７から取得可能である。空調装置７の送風の風速は、空調装置７の風速の代表的な風速を用いてもよい。着衣量は、上記の人の着衣量であり、この着衣量として、温冷感指標算出部４に予め設定された着衣量の代表値が用いられる。活動量（すなわち代謝量）は、上記の人の活動量であり、この活動量として、温冷感指標算出部４に予め設定された活動量の代表値が用いられる。
　温度目標値算出部５は、温冷感指標算出部４の算出結果に基づいて、空調装置７の空調の温度目標値を算出する。より詳細には、温度目標値算出部５は、温冷感指標算出部４で算出されたＰＭＶが上記所定範囲（例えば－０．５以上且つ＋０．５以下）内である場合は、空調装置７の空調の温度目標値を現在の温度目標値に維持（再設定）する。温度目標値算出部５は、温冷感指標算出部４で算出されたＰＭＶが上記所定範囲の上限値（例えば＋０．５）よりも大きい場合は、空調装置７の空調の温度目標値を、現在の温度目標値と比べてＰＭＶが大きいほど小さい値に変更（設定）する。温度目標値算出部５は、温冷感指標算出部４で算出されたＰＭＶが上記所定範囲の下限値（－０．５）よりも小さい場合は、空調装置７の空調の温度目標値を、現在の温度目標値と比べてＰＭＶが小さいほど大きい値に変更（設定）する。
　空調制御部６は、温度目標値算出部５で設定された温度目標値に基づいて、空調装置７を制御する。より詳細には、空調装置７は、室Ｒ１０内の気温が空調装置７に設定された設定温度に近づくように、空調装置７から送風される風の温度を制御する。空調制御部６は、上記の設定温度を、温度目標値算出部５で設定された温度目標値に設定することで、室Ｒ１０内の気温を、室Ｒ１０に居る人が快適に感じる温度に制御する。
　空調装置７は、上述のように、室Ｒ１０内の気温を、空調装置７に設定された設定温度に近づくように、空調装置７から送風される風の温度を制御する。より詳細には、空調装置７は、室Ｒ１０内の気温を検出する温度センサを有し、その温度センサの検出温度が設定温度に近づくように、空調装置７から送風される風の温度を制御する。また、空調装置７は、室Ｒ１０内の湿度を、空調装置７に設定された設定湿度に近づくように、空調装置７から送風される風の湿度を制御する。より詳細には、空調装置７は、室Ｒ１０内の湿度を検出する湿度センサを有し、その湿度センサの検出温度が設定湿度に近づくように、空調装置７から送風される風の湿度を制御する。
　次に図３及び図４を参照して、空間温度推定部２４による熱画像の構図の特定の仕方について説明を補足する。
　本実施形態では、空間温度推定部２４は、設定入力部２１に入力された設置位置情報（空調装置７の設置位置の情報）に基づいて、熱画像取得部２２で取得された熱画像の構図（天井Ｒ１、床Ｒ２及び壁Ｒ３～Ｒ６の映込領域）を特定する。具体的には、図３に示すように、空調装置の設置位置Ｐ１が「中央」（壁Ｒ３の上部の中央）ＰＳ１である場合は、図４Ａに示すように、熱画像Ｇ１の構図は、床Ｒ２及び壁Ｒ５が熱画像Ｇ１における左右方向の中央に配置するように特定される。そして、この特定された構図に基づいて、空間温度推定部２４は、各熱的要因（床Ｒ２及び各壁Ｒ４～Ｒ６）の変化温度を検出する。
　また、図３に示すように、空調装置７の設置位置Ｐ１が「右端」（壁Ｒ３の上部の右端）ＲＲ１である場合は、図４Ｂに示すように、熱画像Ｇ１の構図は、床Ｒ２及び壁Ｒ５が熱画像Ｇ１における左右方向の中央よりも左寄りに配置するように特定される。そして、この特定された構図に基づいて、空間温度推定部２４は、各熱的要因（床Ｒ２及び各壁Ｒ４～Ｒ６）の変化温度を検出する。
　また、図３に示すように、空調装置７の設置位置Ｐ１が「左端」（壁Ｒ３の上部の左端）ＲＬ１である場合は、図４Ｃに示すように、熱画像Ｇ１の構図は、床Ｒ２及び壁Ｒ５が熱画像Ｇ１の中央よりも右寄りに配置するように特定される。そして、この特定された構図に基づいて、空間温度推定部２４は、各熱的要因（床Ｒ２及び各壁Ｒ４～Ｒ６）の変化温度を検出する。
　次に図５を参照して、温冷感推定システム１００の動作を説明する。
　まず、設定入力部２１に間取り情報及び設定位置情報が入力される（Ｓ１）。そして、ＣＲＩ係数決定部２３が、設定入力部２１に入力された間取り情報及び設置位置情報と、空調装置７から取得された空調装置７の送風の速度及び風向の情報とに基づいて、室Ｒ１０の間取り条件を特定する。そして、ＣＲＩ係数決定部２３は、予め設定された複数組のＣＲＩ係数の内から、特定した間取り条件に対応したＣＲＩ係数の組を決定する（Ｓ２）。
　そして、熱画像取得部２２が、室Ｒ１０の各面（天井Ｒ１、床Ｒ２及び壁Ｒ３～Ｒ６）の何れか１面以上の面（例えば床Ｒ２及び壁Ｒ４～Ｒ６の４面）を含む熱画像を取得する（Ｓ３）。そして、人位置検出部２が、熱画像取得部２２で取得された熱画像に基づいて、室Ｒ１０内に居る人の位置Ｘを検出する（Ｓ４）。
　そして、空間温度推定部２４が、室Ｒ１０の各熱的要因（天井Ｒ１、床及Ｒ２び４つの壁Ｒ３～Ｒ４）の変化温度と、人位置検出部２の検出結果（すなわち室Ｒ１０に居る人の位置Ｘ）と、ＣＲＩ係数決定部２３で決定されたＣＲＩ係数とを用いて、室Ｒ１０に居る人の位置Ｘでの気温を推定する（Ｓ５）。
　なお、上記の「室Ｒ１０の各熱的要因の変化温度」は、空間温度推定部２４が、設定入力部２１に入力された上記の間取り情報及び設置位置情報に基づいて熱画像取得部２２で取得された熱画像の構図（天井Ｒ１、床Ｒ２、壁Ｒ３～Ｒ６の映込領域）を特定し、その構図が特定された熱画像から検出する。そして、平均放射温度算出部３が、構図の特定された熱画像に基づいて、室Ｒ１０に居る人の位置（人位置検出部２で検出された位置）Ｘでの平均放射温度（ＭＲＴ）を算出する（Ｓ６）。
　そして、温冷感指標算出部４が、空間温度推定部２４で推定された上記の人の居る位置Ｘの気温と、平均放射温度算出部３で算出された上記の人の居る位置Ｘの平均放射温度と、空調装置７から送風される風の風速と、室Ｒ１０の湿度と、上記の人の着衣量及び活動量とに基づいて、上記の人の温冷感を表す温冷感指標としてＰＭＶを算出する（Ｓ７）。そして、温度目標値算出部５が、温冷感指標算出部４の算出結果に基づいて、空調装置７の空調の温度目標値を算出し（Ｓ８）、空調制御部６が、温度目標値算出部５で算出された温度目標値に基づいて、空調装置７を制御する（Ｓ９）。
　以上のように、本実施形態に空間温度推定システム１によれば、室Ｒ１０の天井Ｒ１、床Ｒ２及び各壁Ｒ３～Ｒ６の少なくとも１面の温度と、上記の天井Ｒ１、床Ｒ２及び各壁Ｒ３～Ｒ６の少なくとも１面に関するＣＲＩ係数とを用いて、室Ｒ１０の室内空間の空気の温度を推定する。このため、各熱的要因（天井Ｒ１、床Ｒ２及び各壁Ｒ３～Ｒ６）の温度が変化した場合でも、各熱的要因の温度変化による影響（熱的寄与量）を足し合わせるだけで（すなわち簡単な計算で）、室Ｒ１０の室内空間の空気の温度を推定できる。
　また、本実施形態に係る温冷感推定システム１００によれば、上記の空間温度推定システム１を用いて、室Ｒ１０に居る人の温冷感をより精度よく推定することができる。
　（変形例）
　以下に説明する変形例は、適宜組み合わせて適用可能である。また、空間温度推定システム１と同様の機能は、空間温度推定方法又はコンピュータプログラムで具現化されてもよい。また、温冷感推定システム１００と同様の機能は、温冷感推定方法又はコンピュータプログラムで具現化されてもよい。
　一態様に係る空間温度推定方法は、検出工程と、空間温度推定工程と、を備える。検出工程は、室Ｒ１０の天井、床Ｒ２及び各壁Ｒ３～Ｒ６の少なくとも１面の温度を検出する。空間温度推定工程は、上記の少なくとも１面に関するＣＲＩ係数を有し、検出工程で検出された温度と上記のＣＲＩ係数とを用いて、室Ｒ１０の室内空間の空気の温度を推定する。
　一態様に係る温冷感推定方法は、上記の空間温度推定方法を用いた温冷感推定方法である。温冷感推定方法は、人位置検出工程と、温冷感推定工程と、を備える。人位置検出工程は、室Ｒ１０内に居る人の位置Ｘを検出する。温冷感推定工程は、上記の空間温度推定方法で推定された人の位置Ｘでの空気の温度を用いて、人の温冷感を表す指標を算出する。
　一態様に係るコンピュータプログラムは、上記の空間温度推定方法をコンピュータシステムに実行させるためのプログラムである。
　一態様に係るコンピュータプログラムは、上記の温冷感推定方法をコンピュータシステムに実行させるためのプログラムである。
　以下に説明する変形例では、上記の実施形態と異なる点を中心に説明する。また、以下の変形例では、上記の実施形態と同じ部分については、同じ符号を付して説明を省略する場合がある。
　（変形例１）
　本変形例は、図６に示すように、上記の実施形態において、設定入力部２１の代わりに、間取り推定部２５を備えている。間取り推定部２５は、熱画像取得部２２で取得された熱画像を用いて、室Ｒ１０の間取り及び空調装置７の設置位置を推定する。すなわち、本変形例では、室Ｒ１０の間取り及び空調装置７の設置位置は、使用者が設定入力部２１に入力するのではなく、間取り推定部２５によって上記の熱画像から推定される。なお、熱画像取得部２２（温度検出手段）と間取り推定部２５は、間取り推定手段４０を構成している。本変形例では、熱画像取得部２２は、室Ｒ１０の間取りを含む熱画像を取得する。本変形例では、熱画像取得部２２は、赤外線アレイセンサであってもよい。
　間取り推定部２５は、熱画像に映る人の位置の履歴（人位置履歴）を用いて、室Ｒ１０の間取り及び空調装置７の設置位置を推定する。なお、人位置履歴とは、一定期間にわたっての熱画像上にプロットした人の位置の集合である。熱画像における人の位置は、熱画像における床Ｒ２に該当する領域上に主に分布し、熱画像における壁Ｒ３～Ｒ６に該当する領域上には殆ど分布しない。この特性を利用すれば、熱画像から、室Ｒ１０の間取り及び空調装置７の設置位置を推定可能である。
　より詳細には、間取り推定部２５は、熱画像から人に該当する領域（すなわち人）を検出し、一定期間にわたって、熱画像において人の位置の履歴（人位置履歴）を蓄積する。そして、間取り推定部２５は、蓄積した人位置履歴に基づいて、室Ｒ１０の間取り及び空調装置７の設置位置を推定する。
　具体的には、空調装置７の設置位置の推定は、下記のように行われる。すなわち、空調装置７の設置位置が「右端」（壁Ｒ３の上部の右端）ＰＲ１である場合（図３参照）は、人位置履歴が熱画像Ｇ１に重ねられると、図７Ａに示すように、熱画像Ｇ１の右側に人が全く存在しない人不在領域Ｆ１が発生する。したがって、間取り推定部２５は、熱画像Ｇ１の右側に人不在領域Ｆ１が発生している場合は、この領域Ｆ１を壁Ｒ６と推定し、熱画像Ｇ１での壁Ｒ６の映込位置から、空調装置７の設置位置を「右端」と推定する。また、空調装置７の設置位置が「中央」（壁Ｒ３の上部の中央）ＰＳ１である場合（図３参照）は、人位置履歴が熱画像Ｇ１に重ねられると、図７Ｂに示すように、熱画像Ｇ１の左右両側に人不在領域Ｆ２，Ｆ３が発生する。したがって、間取り推定部２５は、熱画像Ｇ１の左右両側にそれぞれ領域Ｆ２及び領域Ｆ３が発生している場合は、これら領域Ｆ２，Ｆ３を壁Ｒ４，Ｒ６と推定し、熱画像Ｇ１での壁Ｒ４，Ｒ６の映込位置から、空調装置７の設置位置は、「中央」と推定する。また、空調装置７の設置位置が「左端」（壁Ｒ３の上部の左端）ＰＬ１である場合（図３参照）は、人位置履歴が熱画像Ｇ１に重ねられると、図７Ｃに示すように、熱画像Ｇ１の左側に人不在領域Ｆ４が発生する。したがって、間取り推定部２５は、熱画像Ｇ１の左側に人不在領域Ｆ４が発生している場合は、この領域Ｆ４を壁Ｒ４と推定し、熱画像Ｇ１での壁Ｒ４の映込位置から、空調装置７の設置位置を「左端」と推定する。
　また、室Ｒ１０の間取りの推定は、下記のように行われる。すなわち、間取り推定部２５は、人位置履歴が重ねられた熱画像Ｇ１において、熱画像Ｇ１の下辺から上方向に最も離れた人の位置ＸＰ１と熱画像Ｇ１の下辺との間隔Ｈ１が床Ｒ２の奥行きＷ１（図３参照）であると推定する（図７Ａ～図７Ｃ参照）。室Ｒ１０の床面積は、上記の実施形態で述べた通り、空調装置７の空調能力から推定される。間取り推定部２５は、推定された床面積と奥行きＷ１とから、室Ｒ１０の横幅Ｗ２（図３参照）を推定し、推定した横幅Ｗ２と奥行きＷ１とから室Ｒ１０の床Ｒ２の床形状（横長矩形、縦長矩形又は正方形）を推定する。
　この変形例によれば、間取り推定部２５によって室Ｒ１０の間取り及び空調装置７の設置位置が熱画像から推定されるため、使用者が室Ｒ１０の間取り情報及び空調装置７の設置位置情報を設定入力する必要がない。
　本変形例では、空間温度推定部２４は、間取り推定部２５で推定された室Ｒ１０の間取り及び空調装置７の設置位置に応じて、空間温度推定で用いるＣＲＩ係数を変更する。
　（変形例２）
　変形例１では、室Ｒ１０の間取り及び空調装置７の設置位置の推定の処理は全て、間取り推定部２５内（すなわち空間温度推定システム１内）で行われる。ただし、図８に示すように、熱画像Ｇ１を間取り推定部２５から通信回線網ＣＲ１（クラウド）を介して外部装置２７に送信し、外部装置２７で熱画像Ｇ１（図９Ａ）を復元（高解像度化）して復元画像Ｇ２（図９Ｂ）を生成してもよい。そして、その復元画像Ｇ２を外部装置２７から間取り推定部２５に返送し、間取り推定部２５でその復元画像Ｇ２に基づいて室Ｒ１０の間取り及び空調装置７の設置位置を推定してもよい。
　この場合、復元画像Ｇ２では、復元によって例えば床Ｒ２の輪郭が明確に可視化される。間取り推定部２５は、例えば、復元画像Ｇ２の床Ｒ２の輪郭に基づいて、床Ｒ２と各壁Ｒ４～Ｒ６の境界線を推定して、復元画像Ｇ２の構図（床Ｒ２及び各壁Ｒ４～Ｒ６の映込領域）を特定する。そして、間取り推定部２５は、特定した構図から、空調装置７の設置位置（右端、中央又は左端）及び室Ｒ１０の間取り（縦長矩形、横長矩形又は正方形）を推定する。
　より詳細には、室Ｒ１０の間取り（床形状）は下記のように推定可能である。すなわち、間取り推定部２５は、特定した構図から室Ｒ１０の奥行きＷ１（図３参照）を求め、床Ｒ２の床面積と奥行きＷ１とから床Ｒ２の横幅Ｗ２（図３参照）を求め、求めた奥行きＷ１及び横幅Ｗ２から、床Ｒ２の床形状を求める。また、空調装置７の設置位置は下記のように推定可能である。すなわち、間取り推定部２５は、例えば、特定した構図における壁Ｒ５（すなわち空調装置７の正面の壁）の左右方向の位置に応じて、空調装置７の設置位置を推定する。
　この構成によれば、外部装置２７を用いて熱画像Ｇ１を復元し、間取り推定部２５は、復元画像Ｇ２を用いて室Ｒ１０の間取り及び空調装置７の設置位置を推定することができる。このため、より精度よく、室Ｒ１０の間取り及び空調装置７の設置位置を推定できる。
　（変形例３）
　本変形例に係る空間温度推定システム１は、図１０に示すように、上記の実施形態において、熱画像取得部２２（温度検出手段）の撮像方向（検出方向）を走査する走査装置２８を更に備える。走査装置２８は、回転軸部と、回転駆動源（例えば電動機）とを有する。回転軸部は、回転駆動源によって回転可能である。回転軸部には、熱画像取得部２２が固定されている。回転軸部の回転に伴って熱画像取得部２２が回転軸部の回りに回転し、熱画像取得部２２の撮像方向（検出方向）が例えば室Ｒ１０の左右方向に走査される。回転駆動源は、回転軸部を正転方向及び逆転方向に交互に回転させることで、熱画像取得部２２の撮像方向を左右方向に走査する。これにより、熱画像取得部２２の撮像範囲（検出範囲）をより広範囲に広げることができる。
　（変形例４）
　本変形例は、図１１に示すように、変形例１において、複数の熱画像取得部２２（熱画像取得手段）と、画像合成部３１とを備えている。複数の熱画像取得部２２は、室Ｒ１０内の異なる領域（より詳細には室Ｒ１０の間取りの異なる領域）を撮像する。画像合成部３１は、複数の熱画像取得部２２の各々で取得された熱画像を画像合成して単一の合成熱画像を生成する。本変形例の間取り推定部２５は、画像合成部３１で生成された合成熱画像を用いて、室Ｒ１０の間取り及び空調装置７の設置位置を推定する。本変形例では、複数の熱画像取得部２２、画像合成部３１、及び間取り推定部２５は、間取り推定手段４０を構成している。
　本変形例では、人位置検出部２は、画像合成部３１で生成された合成熱画像を用いて、室Ｒ１０に居る人の位置Ｘを検出する。空間温度推定部２４は、画像合成部３１で生成された合成熱画像を用いて、室Ｒ１０の各熱的要因の変化温度を検出する。平均放射温度算出部３は、画像合成部３１で生成された合成熱画像を用いて、室Ｒ１０に居る人の位置（すなわち人位置検出部２で検出された人の位置）Ｘでの平均放射温度を算出する。
　（変形例５）
　本変形例は、図１２に示すように、変形例１において、室Ｒ１０の間取りを含む撮像画像を取得する撮像装置３２を更に備えている。撮像装置３２は、空調装置７に設置されるか、空調装置７の周囲に設置されている。本変形例の間取り推定部２５は、熱画像取得部２２で取得された熱画像の代わりに、撮像装置３２で撮像された撮像画像を用いて、室Ｒ１０の間取り及び空調装置７の設置位置を推定する。本変形例では、間取り推定手段４０は、間取り推定部２５と、撮像装置３２とで構成される。
　撮像装置３２は、例えばＣＣＤなどの撮像素子と撮像レンズとを備えている。本変形例では、間取り推定部２５は、撮像装置３２の撮像画像を二値化処理して室Ｒ１０の天井Ｒ１、床Ｒ２及び各壁Ｒ３～Ｒ６の境界線を抽出する。そして、間取り推定部２５は、抽出した境界線に基づいて、撮像画像の構図（天井Ｒ１、床Ｒ２及び各壁Ｒ３～Ｒ６の映込領域）を特定し、特定した構図に基づいて室Ｒ１０の間取り及び空調装置７の設置位置を推定する。
　（変形例６）
　本変形例は、図１３に示すように、変形例４において、撮像装置３２の代わりに、室Ｒ１０の間取りを含む距離分布画像を取得するＴＯＦ（Time Of Flight）カメラ３３を備えている。ＴＯＤカメラ３３は、空調装置７に設置されるか、空調装置７の周囲に設置されている。ＦＯＴカメラ３３は、パルス発光された光（例えば近赤外光）の被写体からの反射時間によって被写体までの距離を画素毎に計測し、各画素を距離に応じて色分けして表した画像（離分布画像）を生成するカメラである。本変形例では、間取り推定手段４０は、ＴＯＦカメラ３３と、間取り推定部２５とで構成される。
　本変形例では、間取り推定部２５は、ＴＯＦカメラ３３の撮像画像（距離分布画像）から室Ｒ１０の天井Ｒ１、床Ｒ２及び各壁Ｒ３～Ｒ６の境界線を抽出する。そして、間取り推定部２５は、抽出した境界線に基づいて、撮像画像の構図（天井Ｒ１、床Ｒ２及び各壁Ｒ３～Ｒ６の映込領域）を特定し、特定した構図に基づいて室Ｒ１０の間取り及び空調装置７の設置位置を推定する。
　（まとめ）
　第１の態様に係る空間温度推定システム（１）は、温度検出手段（２２）と、空間温度推定手段（２４）と、を備える。温度検出手段（２２）は、室（Ｒ１０）の天井（Ｒ１）、床（Ｒ２）及び複数の壁（Ｒ３～Ｒ６）の少なくとも１面の温度を検出する。空間温度推定手段（２４）は、上記の少なくとも１面に関するＣＲＩ係数、及び温度検出手段（２２）で検出された温度を用いて、室（Ｒ１０）の室内空間の空気の温度を推定する。
　この構成によれば、室（Ｒ１０）の天井（Ｒ１）、床（Ｒ２）及び各壁（Ｒ３～Ｒ６）の少なくとも１面の温度と、上記の少なくとも１面に関するＣＲＩ係数とを用いて、室（Ｒ１０）の室内空間の空気の温度を推定する。このため、各熱的要因（天井（Ｒ１）、床（Ｒ２）及び各壁（Ｒ３～Ｒ６））の温度が変化した場合でも、各熱的要因の温度変化による影響を足し合わせるだけで（すなわち簡単な計算で）、室（Ｒ１０）の室内空間の空気の温度を推定できる。
　第２の態様に係る空間温度推定システム（１）では、第１の態様において、温度検出手段（２２）は、天井（Ｒ１）、床（Ｒ２）及び複数の壁（Ｒ４～Ｒ６）のうち、互いに異なる向きを向いた２つ以上の面の温度をそれぞれ検出する。
　この構成によれば、室（Ｒ１０）の室内空間の空気の温度の推定精度を向上できる。
　第３の態様に係る空間温度推定システム（１）では、第１の態様において、温度検出手段（２２）は、天井（Ｒ１）、床（Ｒ２）及び複数の壁（Ｒ４～Ｒ６）のうち、互いに異なる向きを向いた２つ以上の面の温度をそれぞれ検出する。
　この構成によれば、室（Ｒ１０）の室内空間の空気の温度の推定精度を更に向上できる。
　第４の態様に係る空間温度推定システム（１）では、第１～第３の態様の何れか１つの態様において、空間温度推定手段（２４）は、間取りが異なる複数の室（Ｒ１０）と１対１に対応する複数のＣＲＩ係数を記憶する。
　この構成によれば、室（Ｒ１０）の間取りに応じて最適なＣＲＩ係数を用いることができる。この結果、室（Ｒ１０）の室内空間の空気の温度の推定精度を更に向上できる。
　第５の態様に係る空間温度推定システム（１）では、第１～第４の態様の何れか１つの態様において、温度検出手段（２２）は、赤外線アレイセンサである。
　この構成によれば、温度検出手段（２２）として赤外線アレイセンサを用いることができる。
　第６の態様に係る空間温度推定システム（１）は、第５の態様において、温度検出手段（２２）の検出方向を走査する走査装置（２８）を備える。
　この構成によれば、温度検出手段（２２）の検出方向を走査することができる。この結果、温度検出手段（２２）の検出範囲をより広範囲に広げることができる。
　第７の態様に係る空間温度推定システム（１）では、第４の態様において、室（Ｒ１０）の間取りを推定する間取り推定手段（４０）を更に備える。空間温度推定手段（２４）は、間取り推定手段（４０）の推定結果に応じてＣＲＩ係数を変更する。
　この構成によれば、室（Ｒ１０）の間取りに応じてＣＲＩ係数が変更されるため、より精度よく、室（Ｒ１０）の室内空間の空気の温度を推定できる。また、室（Ｒ１０）の間取り情報を手動で設定入力する手間が省ける。
　第８の態様に係る空間温度推定システム（１）では、第７の態様において、間取り推定手段（４０）は、室（Ｒ１０）の間取りを含む撮像画像を取得する撮像装置（３２）を含む。
　この構成によれば、撮像装置（３２）の撮像画像を用いて室（Ｒ１０）の間取りを推定できる。
　第９の態様に係る空間温度推定システム（１）では、第７の態様において、間取り推定手段（４０）は、室（Ｒ１０）の間取りを含む距離分布画像を取得するＴＯＦカメラ（３３）を含む。
　この構成によれば、ＴＯＦカメラ（３３）の距離分布画像を用いて室（Ｒ１０）の間取りを推定できる。
　第１０の態様に係る空間温度推定システム（１）では、第７の態様において、間取り推定手段（４０）は、室（Ｒ１０）の間取りを含む熱画像を取得する熱画像取得手段（２２）を含む。
　この構成によれば、熱画像を用いて室（Ｒ１０）の間取りを推定できる。
　第１１の態様に係る空間温度推定システム（１）では、第１０の態様において、間取り推定手段（４０）は、熱画像取得手段（２２）を複数含む。間取り推定手段（４０）は、複数の熱画像取得手段（２２）の各々で取得された熱画像を画像合成し、画像合成した熱画像に基づいて室（Ｒ１０）の間取りを推定する。
　この構成によれば、複数の熱画像取得手段（２２）を用いて広範囲の熱画像を取得することができる。
　第１２の態様に係る温冷感推定システム（１００）では、第１～第１１の態様の何れか１つの態様に係る空間温度推定システム（１）と、人位置検出手段（２）と、温冷感指標算出手段（４）と、を備える。人位置検出手段（２）は、室（Ｒ１０）に居る人の位置（Ｘ）を検出する。温冷感指標算出手段（４）は、空間温度推定システム（１）で推定された人の位置（Ｘ）での空気の温度を用いて、人の温冷感を表す指標（例えばＰＭＶ）を算出する。
　この構成によれば、上記の空間温度推定システム（１）を用いて、室（Ｒ１０）に居る人の温冷感を推定することができる。
　第１３の態様に係る温冷感推定システムは、第１２の態様において、室（Ｒ１０）に居る人の位置（Ｘ）での平均放射温度を算出する平均放射温度算出手段（３）を更に備える。温冷感指標算出手段（４）は、平均放射温度算出手段（３）によって取得された平均放射温度を更に用いて、人の温冷感を表す指標（例えばＰＭＶ）を算出する。
　この構成によれば、室（Ｒ１０）に居る人の位置での平均放射温度を更に用いて、室（Ｒ１０）に居る人の温冷感を表す指標を算出するため、室（Ｒ１０）に居る人の温冷感をより精度よく推定できる。
　第１４の態様に係る空間温度推定方法は、温度検出工程と、空間温度推定工程と、を備える。温度検出工程は、室（Ｒ１０）の天井（Ｒ１）、床（Ｒ２）及び複数の壁（Ｒ３～Ｒ６）の少なくとも１面の温度を検出する。空間温度推定工程は、少なくとも１面に関するＣＲＩ係数、及び温度検出手段（２２）で検出された温度を用いて、室（Ｒ１０）の室内空間の空気の温度を推定する。
　この構成によれば、室（Ｒ１０）の天井（Ｒ１）、床（Ｒ２）及び各壁（Ｒ３～Ｒ６）の少なくとも１面の温度と、上記の少なくとも１面に関するＣＲＩ係数とを用いて、室（Ｒ１０）の室内空間の空気の温度を推定する。このため、各熱的要因（天井（Ｒ１）、床（Ｒ２）及び各壁（Ｒ３～Ｒ６））の温度が変化した場合でも、各熱的要因の温度変化による影響を足し合わせるだけで（すなわち簡単な計算で）、室（Ｒ１０）の室内空間の空気の温度を推定できる。
　第１５の態様に係る温冷感推定方法は、第１４の態様に係る空間温度推定方法を用いた温冷感推定方法である。温冷感推定方法は、人位置検出工程と、温冷感指標算出工程と、を備える。人位置検出工程は、室（Ｒ１０）内に居る人の位置（Ｘ）を検出する。温冷感指標算出工程は、空間温度推定方法によって推定された人の位置（Ｘ）での空気の温度を用いて、人の温冷感を表す指標を算出する。
　この構成によれば、上記の空間温度推定方法を用いて、室（Ｒ１０）内に居る人の温冷感を表す指標を算出することができる。
　第１６の態様に係るプログラムは、第１４の態様に係る空間温度推定方法をコンピュータシステムに実行させるためのプログラムである。
　この構成によれば、上記の空間温度推定方法をコンピュータシステムに実行させるためのプログラムを提供できる。
　第１７の態様に係るプログラムは、第１５の態様に係る温冷感推定方法をコンピュータシステムに実行させるためのプログラムである。
　この構成によれば、上記の温冷感推定方法をコンピュータシステムに実行させるためのプログラムを提供できる。
　１　空間温度推定システム
　３　平均放射温度算出部（平均放射温度算出手段）
　４　温冷感指標算出部（温冷感指標算出手段）
　２２　熱画像取得部（温度検出手段、熱画像取得手段）
　２４　空間温度推定部（空間温度推定手段）
　２９　走査装置
　３２　撮像装置
　３３　ＴＯＦカメラ
　４０　間取り推定手段
　１００　温冷感推定システム
　Ｒ１　天井
　Ｒ２　床
　Ｒ３～Ｒ６　壁
　Ｒ１０　室
　Ｘ　人の位置

Claims (17)

  1. 　室の天井、床及び複数の壁の少なくとも１面の温度を検出する温度検出手段と、
    　前記少なくとも１面に関するＣＲＩ係数、及び前記温度検出手段で検出された前記温度を用いて、前記室の室内空間の空気の温度を推定する空間温度推定手段と、を備える、
    空間温度推定システム。
  2. 　前記温度検出手段は、前記天井、前記床及び前記複数の壁のうち、互いに異なる向きを向いた２つ以上の面の温度をそれぞれ検出する、
    請求項１に記載の空間温度推定システム。
  3. 　前記温度検出手段は、前記天井、前記床及び前記複数の壁のうち、互いに異なる向きを向いた３つ以上の面の温度をそれぞれ検出する、
    請求項１に記載の空間温度推定システム。
  4. 　前記空間温度推定手段は、間取りが異なる複数の前記室と一対一に対応する複数の前記ＣＲＩ係数を記憶する、
    請求項１～３の何れか１項に記載の空間温度推定システム。
  5. 　前記温度検出手段は、赤外線アレイセンサである、
    請求項１～４の何れか１項に記載の空間温度推定システム。
  6. 　前記温度検出手段の検出方向を走査する走査装置を更に備える、
    請求項５に記載の空間温度推定システム。
  7. 　前記室の間取りを推定する間取り推定手段を更に備え、
    　前記空間温度推定手段は、前記間取り推定手段の推定結果に応じて前記ＣＲＩ係数を変更する、
    請求項４に記載の空間温度推定システム。
  8. 　前記間取り推定手段は、前記室の間取りを含む撮像画像を取得する撮像装置を含む、
    請求項７に記載の空間温度推定システム。
  9. 　前記間取り推定手段は、前記室の間取りを含む距離分布画像を取得するＴＯＦカメラを含む、
    請求項７に記載の空間温度推定システム。
  10. 　前記間取り推定手段は、前記室の間取りを含む熱画像を取得する熱画像取得手段を含む、
    請求項７に記載の空間温度推定システム。
  11. 　前記間取り推定手段は、前記熱画像取得手段を複数含み、
    　前記間取り推定手段は、前記複数の熱画像取得手段の各々で取得された熱画像を画像合成し、画像合成した熱画像に基づいて前記室の間取りを推定する、
    請求項１０に記載の空間温度推定システム。
  12. 　請求項１～１１の何れか１項に空間温度推定システムと、
    　前記室に居る人の位置を検出する人位置検出手段と、
    　前記空間温度推定システムで推定された前記人の位置での空気の温度を用いて、前記人の温冷感を表す指標を算出する温冷感指標算出手段と、を備える、
    温冷感推定システム。
  13. 　前記室に居る前記人の位置での平均放射温度を算出する平均放射温度算出手段を更に備え、
    　前記温冷感指標算出手段は、前記平均放射温度算出手段によって算出された前記平均放射温度を更に用いて、前記人の温冷感を表す指標を算出する、
    請求項１２に記載の温冷感推定システム。
  14. 　前記室の天井、床及び複数の壁の少なくとも１面の温度を検出する温度検出工程と、
    　前記少なくとも１面に関するＣＲＩ係数、前記温度検出工程で検出された前記温度を用いて、前記室の室内空間の空気の温度を推定する空間温度推定工程と、を備える、
    空間温度推定方法。
  15. 　請求項１４に記載の空間温度推定方法を用いた温冷感推定方法であって、
    　前記室内に居る人の位置を検出する人位置検出工程と、
    　前記空間温度推定方法によって推定された前記人の位置での空気の温度を用いて、前記人の温冷感を表す指標を算出する温冷感指標算出工程と、を備える、
    温冷感推定方法。
  16. 　請求項１４に記載の空間温度推定方法をコンピュータシステムに実行させるための、
    プログラム。
  17. 　請求項１５に記載の温冷感推定方法をコンピュータシステムに実行させるための、
    プログラム。
PCT/JP2020/039188 2019-12-10 2020-10-16 空間温度推定システム、温冷感推定システム、空間温度推定方法、温冷感推定方法、及び、プログラム WO2021117343A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN202080084044.0A CN114829871A (zh) 2019-12-10 2020-10-16 空间温度估计系统、暖/冷感估计系统、空间温度估计方法、暖/冷感估计方法和程序
US17/783,529 US20220398764A1 (en) 2019-12-10 2020-10-16 Spatial temperature estimation system, warm/cold sensation estimation system, spatial temperature estimation method, warm/cold sensation estimation method, and program
JP2021563767A JPWO2021117343A1 (ja) 2019-12-10 2020-10-16

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2019-223242 2019-12-10
JP2019223242 2019-12-10

Publications (1)

Publication Number Publication Date
WO2021117343A1 true WO2021117343A1 (ja) 2021-06-17

Family

ID=76328876

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2020/039188 WO2021117343A1 (ja) 2019-12-10 2020-10-16 空間温度推定システム、温冷感推定システム、空間温度推定方法、温冷感推定方法、及び、プログラム

Country Status (4)

Country Link
US (1) US20220398764A1 (ja)
JP (1) JPWO2021117343A1 (ja)
CN (1) CN114829871A (ja)
WO (1) WO2021117343A1 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114235892A (zh) * 2021-12-17 2022-03-25 南京林业大学 一种家居木制品表面冷暖感的测试方法
WO2022269781A1 (ja) * 2021-06-23 2022-12-29 三菱電機株式会社 温度測定装置、温度測定方法、及び電気機器

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03291445A (ja) * 1990-04-06 1991-12-20 Shimizu Corp 温熱環境評価支援システム
JP2008261567A (ja) * 2007-04-12 2008-10-30 Mitsubishi Electric Corp 空気調和機
JP2010091253A (ja) * 2008-09-10 2010-04-22 Mitsubishi Electric Corp 空気調和機
WO2015122201A1 (ja) * 2014-02-17 2015-08-20 パナソニック株式会社 空気調和機及び熱画像センサシステム
JP2019028006A (ja) * 2017-08-03 2019-02-21 アズビル株式会社 放射環境評価システムおよび放射環境評価方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03291445A (ja) * 1990-04-06 1991-12-20 Shimizu Corp 温熱環境評価支援システム
JP2008261567A (ja) * 2007-04-12 2008-10-30 Mitsubishi Electric Corp 空気調和機
JP2010091253A (ja) * 2008-09-10 2010-04-22 Mitsubishi Electric Corp 空気調和機
WO2015122201A1 (ja) * 2014-02-17 2015-08-20 パナソニック株式会社 空気調和機及び熱画像センサシステム
JP2019028006A (ja) * 2017-08-03 2019-02-21 アズビル株式会社 放射環境評価システムおよび放射環境評価方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
SASAMOTO, TARO: "Control of indoor thermal environment based on the concept of contribution ratio of indoor climate", JOURNAL OF ENVIRONMENTAL ENGINEERING (TRANSACTIONS OF AIJ), 2004, pages 33 - 38 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022269781A1 (ja) * 2021-06-23 2022-12-29 三菱電機株式会社 温度測定装置、温度測定方法、及び電気機器
JP7442743B2 (ja) 2021-06-23 2024-03-04 三菱電機株式会社 温度測定装置、温度測定方法、及び電気機器
CN114235892A (zh) * 2021-12-17 2022-03-25 南京林业大学 一种家居木制品表面冷暖感的测试方法

Also Published As

Publication number Publication date
US20220398764A1 (en) 2022-12-15
JPWO2021117343A1 (ja) 2021-06-17
CN114829871A (zh) 2022-07-29

Similar Documents

Publication Publication Date Title
JP5238679B2 (ja) 空調制御装置、空調制御方法及び輻射温度計測装置
WO2021117343A1 (ja) 空間温度推定システム、温冷感推定システム、空間温度推定方法、温冷感推定方法、及び、プログラム
CN101900400B (zh) 空气调节器
JP5085670B2 (ja) 空調制御システムおよび空調制御方法
CN101672498B (zh) 空调机
KR102130402B1 (ko) 열화상 카메라를 이용한 실내 온열 환경 평가 장치 및 그 방법
CN106247525A (zh) 空气调节机、传感系统及其温冷感推定方法
US20130162836A1 (en) Calculating energy dissipation in an ir image
KR20120050989A (ko) 공기 조화기
JP2017015384A (ja) 空調制御装置
JP2021152416A (ja) 空調制御装置、空調制御方法、及び空調制御システム
JP2004078297A (ja) 情報処理装置、情報処理方法及び環境制御装置
JP2016205633A (ja) 空調制御装置
JP2010216713A (ja) 空気調和機
KR102246575B1 (ko) 열화상 카메라를 이용한 실내 표면 온도 측정 장치 및 그 방법
JP7372201B2 (ja) 空調制御方法及びそのシステム
CN106196490B (zh) 空调器的控制方法及空调器
JP2014142099A (ja) 空気調和機
US20210396585A1 (en) Thermal quality mappings
JP2022149774A (ja) 情報処理装置およびプログラム
JP7389361B2 (ja) 情報処理装置、空気調和システムおよびプログラム
JP7199593B2 (ja) 人検知装置、人検知システム、設備機器システム、人検知方法及びプログラム
JPWO2018069961A1 (ja) 空気調和装置
WO2021210168A1 (ja) 空調制御システム、空調システム、空調制御方法及び空調制御プログラム
JPWO2021038696A1 (ja) 空気調和装置

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: 20900454

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2021563767

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 20900454

Country of ref document: EP

Kind code of ref document: A1