JP2503965B2 - End-tidal curve automatic recognition device - Google Patents

End-tidal curve automatic recognition device

Info

Publication number
JP2503965B2
JP2503965B2 JP61100231A JP10023186A JP2503965B2 JP 2503965 B2 JP2503965 B2 JP 2503965B2 JP 61100231 A JP61100231 A JP 61100231A JP 10023186 A JP10023186 A JP 10023186A JP 2503965 B2 JP2503965 B2 JP 2503965B2
Authority
JP
Japan
Prior art keywords
curve
gas
exhalation
monitor
exhalation monitor
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP61100231A
Other languages
Japanese (ja)
Other versions
JPS62254731A (en
Inventor
信行 山岡
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shimadzu Corp
Original Assignee
Shimadzu Corp
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 Shimadzu Corp filed Critical Shimadzu Corp
Priority to JP61100231A priority Critical patent/JP2503965B2/en
Publication of JPS62254731A publication Critical patent/JPS62254731A/en
Application granted granted Critical
Publication of JP2503965B2 publication Critical patent/JP2503965B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/507Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for determination of haemodynamic parameters, e.g. perfusion CT

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Engineering & Computer Science (AREA)
  • Radiology & Medical Imaging (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Dentistry (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Investigating Or Analysing Biological Materials (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Nuclear Medicine (AREA)

Description

【発明の詳細な説明】 産業上の利用分野 この発明は、核医学の分野で、133Xeガスなどの放射
性ガスを用いて行なわれるCBF測定(脳血流量測定)に
おいて動脈血中の放射性物質の量を示すエンドタルダル
カーブ(Endtidal Curve)を呼気モニタカーブから自動
的に認識する装置に関する。
TECHNICAL FIELD The present invention is used in the field of nuclear medicine to measure the amount of radioactive substances in arterial blood in CBF measurement (cerebral blood flow measurement) performed using radioactive gas such as 133 Xe gas. The present invention relates to a device for automatically recognizing an Endtidal Curve showing from the expiration monitor curve.

従来の技術 核医学の分野で、CBF測定は脳梗塞などの能の血管障
害を早期に発見するものとして期待されている。放射線
ガスを用いて行なうCBF測定は、133Xeガスなどの放射性
ガスを被検者の呼吸器官に送り込み、呼吸器官によって
RI(放射性同位元素)を血液中に取り入れさせ、動脈血
中のRI(放射性同位元素)量を測定するとともに、他方
で、リング型ECT装置(エミッションコンピュータ断層
撮影装置)により脳の特定の部位の血液中のRI量を測定
して、上記の動脈血中のRI量と比較することにより行な
われる。
2. Description of the Related Art In the field of nuclear medicine, CBF measurement is expected as an early detection of functional vascular disorders such as cerebral infarction. CBF measurement performed using radiation gas involves sending radioactive gas such as 133 Xe gas to the respiratory organs of a subject,
RI (radioisotope) is introduced into blood to measure the amount of RI (radioisotope) in arterial blood, while on the other hand, a ring-type ECT device (emission computed tomography device) is used to measure blood in a specific region of the brain. It is performed by measuring the amount of RI in the blood and comparing it with the amount of RI in the arterial blood.

この場合、動脈血中のRI量のカーブ(エンドタイダル
カーブ)は呼気モニタカーブから推測することができ
る。そのため、呼気中のRI量を呼気モニタ検出器で検出
して呼気モニタカーブを作成し、これからエンドタイダ
ルカーブを作成する作業が必要である。
In this case, the curve of the RI amount in arterial blood (end tidal curve) can be estimated from the exhalation monitor curve. Therefore, it is necessary to create an exhalation monitor curve by detecting the RI amount in exhalation with an exhalation monitor detector and then create an end tidal curve.

従来では被検者の呼吸が一定間隔であるとして一定間
隔で呼気モニタカーブをサンプリングし、エンドタイダ
ルカーブを作ったり、呼気モニタカーブを微分して変曲
点を求めてエンドタイダルカーブを作ったりすることが
試みられている。
Conventionally, an exhalation monitor curve is sampled at regular intervals assuming that the subject's breathing is at regular intervals to create an end tidal curve, or the exhalation monitor curve is differentiated to obtain an inflection point to create an end tidal curve. Is being attempted.

発明が解決しようとする問題点 しかし、被験者の呼吸(肺機能)は人によってかなり
異なり、また同一人でも一定であるとは限らない。ま
た、吸気モニタカーブには雑音成分が多く含まれてお
り、これを微分することは雑音を強調することになる。
Problems to be Solved by the Invention However, the breathing (lung function) of a subject varies considerably from person to person, and is not always constant in the same person. Also, the intake monitor curve contains a lot of noise components, and differentiating this component emphasizes noise.

そのため、従来では、呼気モニタカーブから適切なエ
ンドタイダルカーブを認識することが困難であった。
Therefore, conventionally, it was difficult to recognize an appropriate end tidal curve from the exhalation monitor curve.

この発明は、CBF測定における呼気モニタカーブから
適切なエンドタイダルカーブを自動的に認識できる装置
を提供することを目的とする。
An object of the present invention is to provide a device capable of automatically recognizing an appropriate end tidal curve from an exhalation monitor curve in CBF measurement.

問題点を解決するための手段 この発明のエンドタイダルカーブ自動認識装置は、放
射性ガスを用いて行なわれるCBF測定において呼気モニ
タカーブを作成する手段と、該呼気モニタカーブから平
均カーブを算出する手段と、上記呼気モニタカーブから
上記平均カーブを減算する手段と、放射性ガスを投与し
ている期間について上記の減算結果として得られるカー
ブにおける極小位置を求めるとともに該ガスの投与後の
期間について上記カーブにおける極大位置を求める手段
と、上記呼気モニタカーブにおける上記極小位置および
極大位置の値を求めてこれらを結ぶカーブを作成する手
段とからなる。
Means for Solving Problems The automatic end tidal curve recognition device of the present invention comprises means for creating an exhalation monitor curve in CBF measurement performed using radioactive gas, and means for calculating an average curve from the exhalation monitor curve. A means for subtracting the average curve from the exhalation monitor curve and a minimum position in the curve obtained as a result of the subtraction for a period during which radioactive gas is being administered, and a maximum in the curve for a period after administration of the gas. It comprises means for obtaining the position and means for obtaining the values of the minimum position and the maximum position in the exhalation monitor curve and creating a curve connecting them.

作用 呼気モニタカーブから平均カーブを算出し、呼気モニ
タカーブより平均カーブを減算して、この減算結果とし
て得られるカーブにおける、放射性ガスを投与している
期間については極小位置を求め、放射性ガスの投与後の
期間については極大位置を求め、これら極小位置、極大
位置での呼気モニタカーブの値(極小値、極大値)を求
めて、これら極小値、極大値を結んでカーブを作成して
いるので、133Xeガス投与中および投与後の呼気時の133
Xe濃度をプロットしてカーブを作成したこととなって、
こうして得たカーブは適切なエンドタイダルカーブとい
うことになる。
Action Calculate the average curve from the exhalation monitor curve, subtract the average curve from the exhalation monitor curve, find the minimum position for the period during which radioactive gas is administered in the curve obtained as the result of this subtraction, and administer the radioactive gas. For the subsequent period, the maximum position is calculated, the values of the exhalation monitor curve (minimum value, maximum value) at these minimum position and maximum position are calculated, and the curve is created by connecting these minimum value and maximum value. 133 133 During and after administration of Xe gas 133
It means that the curve was created by plotting the Xe concentration.
The curve thus obtained is a proper end tidal curve.

実 施 例 第1図において、患者1に対してガスコントローラ
(たとえば安西総業株式会社製「キセノンガスコントロ
ールシステムAZ−701−NTS)2で133Xeガスを投与し、
ガスコントローラ2に内蔵の呼気モニタ検出器24で呼気
中に含まれる133Xe濃度を検出する。ガスコントローラ
2内のスパイロメータ21にあらかじめ133Xeガスが入れ
られており、三方バルブ22を操作することによって空気
または133Xeガスが患者1にセットされたマスク23に送
られ、呼吸器官を介して患者1に対して133Xeガスの投
与が行なわれる。
Example 1 In FIG. 1, a patient 1 was administered 133 Xe gas with a gas controller (for example, “Xenon gas control system AZ-701-NTS” manufactured by Anzai Sogyo Co., Ltd.),
The exhalation monitor detector 24 incorporated in the gas controller 2 detects the 133 Xe concentration contained in the exhaled breath. The spirometer 21 in the gas controller 2 is prefilled with 133 Xe gas, and air or 133 Xe gas is sent to the mask 23 set on the patient 1 by operating the three-way valve 22, and is transmitted through the respiratory organs. Patient 1 is administered 133 Xe gas.

呼気モニタ検出器24から出力される133Xe濃度信号は
呼気モニタカーブ作成器3に送られ、上記濃度信号を、
時間軸を横軸にとったグラフ上にプロットすることによ
り、たとえば第2図のような呼気モニタカーブAが作成
される。この第2図では、T1の時点で三方バルブ22が操
作され、133Xeガスの患者1に対する投与が開始される
ものとしている。この133Xeガス投与中、まず、患者1
の呼吸運動の吸気動作で濃いガスが検出器24および肺に
入るため、T2の時点ではガス濃度が高くなる。次にT3の
時点は呼気過程で肺からガスが吐き出されるが、肺にあ
る空気により薄められるので呼気中のガス濃度は低いも
のとなる。次の吸気過程で再びガスが肺に吸い込まれる
が、この時点T4ではある量が人体中に取り込まれるので
先の時点T2より濃度が下がる。
The 133 Xe concentration signal output from the exhalation monitor detector 24 is sent to the exhalation monitor curve generator 3 and the concentration signal
By plotting the graph on the horizontal axis of the time axis, for example, the exhalation monitor curve A as shown in FIG. 2 is created. In FIG. 2, it is assumed that the three-way valve 22 is operated at time T1 to start administration of 133 Xe gas to the patient 1. During this 133 Xe gas administration, first, patient 1
At the time of T2, the gas concentration becomes high because the dense gas enters the detector 24 and the lungs by the inspiratory action of the respiratory motion of. Next, at the time of T3, gas is exhaled from the lung during the exhalation process, but since it is diluted by the air in the lung, the gas concentration in the exhalation becomes low. In the next inhalation process, gas is again inhaled into the lungs, but at this time T4, a certain amount is taken into the human body, so the concentration becomes lower than at the previous time T2.

こうして133Xeガス投与中は吸気時には濃度が高く、
呼気時には濃度が低いというカーブが繰り返され、ある
時間(約1分間)経過した時点T6で三方バルブ22が操作
されて133Xeガスの投与が停止され、空気が送られるよ
うになる。すると今度は空気が吸われた時に133Xeガス
が下がるので、吸気時に極小値を示し呼気時に極大値を
示す。時点T7は吸気時で、空気が吸われることによって
ガス濃度が下がる。時点T8は吸気時で、肺中の133Xeガ
スが吐き出されるのでガス濃度が高くなる。このような
呼吸運動にともなう山・谷が連続していく。
Thus, the concentration is high during inhalation during 133 Xe gas administration,
During exhalation, the curve of low concentration is repeated, and at a time T6 when a certain time (about 1 minute) has elapsed, the three-way valve 22 is operated to stop the administration of 133 Xe gas and to send air. Then, when the air is sucked in, 133 Xe gas goes down this time, so it shows a minimum value during inspiration and a maximum value during expiration. At time T7, the gas concentration is lowered because the air is sucked in at the time of intake. At time point T8, the gas concentration is high because 133 Xe gas in the lungs is exhaled during inspiration. The mountains and valleys that accompany such breathing exercises continue.

平均カーブ算出器4は、この呼気モニタカーブAをス
ムージングして平均カーブBを求める。スムージングす
る幅は1呼吸間隔である。この呼吸間隔を求めるには、
まず、標準的な呼吸間隔でスムージングし、呼気モニタ
カーブから引き算し、極大値(もしくは極小値)間の間
隔を計測し、それらの平均値をその患者1の呼吸間隔と
する。こうして第2図の点線Bで示すような平均カーブ
を得る。
The average curve calculator 4 smoothes this exhalation monitor curve A to obtain an average curve B. The width of smoothing is one breath interval. To find this breathing interval,
First, the breathing interval of the patient 1 is smoothed at a standard breathing interval, subtracted from the exhalation monitor curve, the interval between maximum values (or minimum values) is measured, and the average value thereof is taken as the breathing interval of the patient 1. In this way, an average curve as shown by the dotted line B in FIG. 2 is obtained.

次に減算器5によって呼気モニタカーブAから平均カ
ーブBを減算して第4図、第6図のような減算カーブD
を得る。第4図の減算カーブDは133Xeガス投与中のも
のであり、第6図の減算カーブDは投与後のものであ
る。
Next, the average curve B is subtracted from the exhalation monitor curve A by the subtractor 5 to obtain a subtraction curve D as shown in FIGS.
Get. The subtraction curve D in FIG. 4 is during administration of 133 Xe gas, and the subtraction curve D in FIG. 6 is after administration.

この減算カーブDに基づき、極大・極小算出器6が極
大点と極小点とを求める。133Xeガス投与中では、第4
図のように減算カーブDの極小となる時刻T3、T5、…が
求められる、133Xeガス投与後では、第6図のように減
算カーブDの極大となる時刻T8、T10、…が求められ
る。
Based on the subtraction curve D, the maximum / minimum calculator 6 obtains the maximum point and the minimum point. No. 4 during 133 Xe gas administration
As shown in the figure, the times T3, T5, ... At which the subtraction curve D becomes the minimum are obtained. After administration of 133 Xe gas, the times T8, T10, ... At which the subtraction curve D becomes the maximum are obtained as shown in FIG. .

次に、第3図に示すような133Xeガス投与中の呼気モ
ニタカーブAの極小時刻T3、T5、…での濃度値f
(3),f(5),…を求め、第5図に示す133Xeガス投
与後の呼気モニタカーブAの極大時刻の濃度値f
(8),f(10),…を求める。これらの点を結ぶことに
よってカーブC(第2図の1点鎖線)が作成される。こ
の動作はエンドタイダルカーブ作成器7により行なわれ
る。
Next, the concentration value f at the minimum times T3, T5, ... Of the exhalation monitor curve A during 133 Xe gas administration as shown in FIG.
(3), f (5), ... Is obtained, and the concentration value f at the maximum time of the exhalation monitor curve A after 133 Xe gas administration shown in FIG.
Find (8), f (10), ... A curve C (one-dot chain line in FIG. 2) is created by connecting these points. This operation is performed by the end tidal curve generator 7.

こうして求めたカーブCは、133Xeガス投与中および
投与後の呼気時の133Xe濃度をプロットしたカーブであ
るから、エンドタイダルカーブを表わすことになる。こ
のように、呼気モニタカーブAからエンドタイダルカー
ブCを自動的に認識できる。そして、上記のように、ス
ムージングと減算の操作のみでよいため、呼気モニタカ
ーブAの雑音成分に影響されることなく最適なエンドタ
イダルカーブCを自動作成できる。
The curve C thus obtained is a curve plotting 133 Xe concentration during and after administration of 133 Xe gas, and thus represents an end tidal curve. In this way, the end tidal curve C can be automatically recognized from the exhalation monitor curve A. As described above, since only smoothing and subtraction operations are required, the optimum end tidal curve C can be automatically created without being affected by the noise component of the exhalation monitor curve A.

なお、この呼気モニタカーブ作成器3、平均カーブ算
出器4、減算器5、極大・極小算出器6、エンドタイダ
ルカーブ作成器7は、コンピュータによりソフト的に構
成することが可能である。
The expiration monitor curve generator 3, the average curve calculator 4, the subtractor 5, the maximum / minimum calculator 6, and the end tidal curve generator 7 can be configured by a computer in software.

発明の効果 この発明によれば、呼気モニタカーブに多くの雑音成
分が含まれていても、また被験者の呼吸間隔にばらつき
があっても、呼気モニタカーブから適切なエンドタイダ
ルカーブを認識することが可能となる。
EFFECTS OF THE INVENTION According to the present invention, it is possible to recognize an appropriate end tidal curve from the exhalation monitor curve even if the exhalation monitor curve contains many noise components and the breathing interval of the subject varies. It will be possible.

【図面の簡単な説明】[Brief description of drawings]

第1図はこの発明の一実施例の模式図、第2図、第3
図、第4図、第5図および第6図は動作説明のための波
形図である。 1……患者、2……ガスコントローラ 21……スパイロメータ、22……三方バルブ 23……マスク、24……呼気モニタ検出器 3……呼気モニタカーブ作成器 4……平均カーブ算出器、5……減算器 6……極大・極小算出器 7……エンドタイダルカーブ作成器
FIG. 1 is a schematic view of an embodiment of the present invention, FIG. 2 and FIG.
FIG. 4, FIG. 5, FIG. 5 and FIG. 6 are waveform diagrams for explaining the operation. 1 ... Patient, 2 ... Gas controller 21 ... Spirometer, 22 ... Three-way valve 23 ... Mask, 24 ... Exhalation monitor detector 3 ... Exhalation monitor curve generator 4 ... Average curve calculator, 5 …… Subtractor 6 …… Maximum / minimum calculator 7 …… End tidal curve generator

───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.6 識別記号 庁内整理番号 FI 技術表示箇所 // G01N 33/497 7638−2J A61B 6/00 350A ─────────────────────────────────────────────────── ─── Continuation of the front page (51) Int.Cl. 6 Identification code Office reference number FI technical display location // G01N 33/497 7638-2J A61B 6/00 350A

Claims (1)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】放射性ガスを用いて行なわれる脳血流量測
定において呼気モニタカーブを作成する手段と、該呼気
モニタカーブから平均カーブを算出する手段と、上記呼
気モニタカーブから上記平均カーブを減算する手段と、
放射性ガスを投与している期間について上記の減算結果
として得られるカーブにおける極小位置を求めるととも
に該ガスの投与後の期間について上記カーブにおける極
大位置を求める手段と、上記呼気モニタカーブにおける
上記極小位置および極大位置の値を求めてこれらを結ぶ
カーブを作成する手段とからなるエンドタイダルカーブ
自動認識装置。
1. A means for creating an exhalation monitor curve in cerebral blood flow measurement performed using radioactive gas, a means for calculating an average curve from the exhalation monitor curve, and a subtraction of the average curve from the exhalation monitor curve. Means and
Means for determining the minimum position in the curve obtained as the result of the subtraction for the period during which the radioactive gas is administered, and means for determining the maximum position in the curve for the period after the administration of the gas, and the minimum position in the exhalation monitor curve and An automatic end-tidal curve recognition device comprising means for obtaining the maximum position value and creating a curve connecting these.
JP61100231A 1986-04-30 1986-04-30 End-tidal curve automatic recognition device Expired - Fee Related JP2503965B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP61100231A JP2503965B2 (en) 1986-04-30 1986-04-30 End-tidal curve automatic recognition device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP61100231A JP2503965B2 (en) 1986-04-30 1986-04-30 End-tidal curve automatic recognition device

Publications (2)

Publication Number Publication Date
JPS62254731A JPS62254731A (en) 1987-11-06
JP2503965B2 true JP2503965B2 (en) 1996-06-05

Family

ID=14268500

Family Applications (1)

Application Number Title Priority Date Filing Date
JP61100231A Expired - Fee Related JP2503965B2 (en) 1986-04-30 1986-04-30 End-tidal curve automatic recognition device

Country Status (1)

Country Link
JP (1) JP2503965B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2829389B2 (en) * 1996-10-16 1998-11-25 工業技術院長 Respiration detector
JP6155004B2 (en) * 2012-10-01 2017-06-28 国立研究開発法人国立循環器病研究センター Face mask mechanism and PET inspection method using the face mask mechanism

Also Published As

Publication number Publication date
JPS62254731A (en) 1987-11-06

Similar Documents

Publication Publication Date Title
JP3761465B2 (en) Non-invasive determination of cardiac output, pulmonary blood flow and blood gas volume
US6648833B2 (en) Respiratory analysis with capnography
JP3521915B2 (en) Apparatus and method for non-invasive measurement of cardiac output
JP3477200B2 (en) In particular, a respiratory assist device for treating sleep apnea
US20070123792A1 (en) System and method for determining airway obstruction
US6200271B1 (en) Bi-directional partial re-breathing method
CN108024758B (en) Simplified display of end-tidal CO2
US6210342B1 (en) Bi-directional partial re-breathing method
EP1229826A2 (en) Method and apparatus for determining cardiac output
EP1238631B1 (en) A method for non-invasively determining conditions in the circulatory system of a subject
US6302851B1 (en) Method and apparatus for determining a pulmonary function parameter for gas exchange
US20130324872A1 (en) Real-time airway check status indicator
US20090308393A1 (en) Medical diagnostic cart and method of use
JP6615317B2 (en) Capnography with decision support system architecture
JPH08243165A (en) Respiratory organ ventilation / anesthesia device
GB2055046A (en) Determining hypersensitivity of the respiratory system
JP2503965B2 (en) End-tidal curve automatic recognition device
CN112754465A (en) Method for estimating quasi-static compliance of lung under pressure-controlled mechanical ventilation
EP0180876A2 (en) Apparatus for measuring xenon concentration in xenon cerebral blood-flow studies
JP2802655B2 (en) Enhanced gas concentration measuring device for measuring enhanced gas concentration in end-tidal gas and medical diagnostic imaging apparatus equipped with such enhanced gas concentration measuring device
Kuhlen et al. Validation and clinical application of a continuous P0. 1 measurement using standard respiratory equipment
JP2001276057A (en) Method for deciding constant in relational expression between end expiratory velocity constant and arterial blood velocity constant and xenon ct device
KR102578530B1 (en) APPARATUS AND METHOD FOR PERFORMING BILEVEL HIGH FLOW THERAPY BASED ON SpO2
CN115714011A (en) Intelligent management system for household remote control oxygen therapy of children suffering from bronchopulmonary dysplasia
Badgwell Respiratory gas monitoring in the pediatric patient

Legal Events

Date Code Title Description
LAPS Cancellation because of no payment of annual fees