Optical determination of white blood cell concentration in fluids

Abstract

The concentration of white blood cells in a fluid is determined by an optical detection apparatus. The apparatus comprises a light emitting arrangement for generating a light beam (300) with a wavelength in the range of 350-575 nm, and a light detection arrangement (31) for detecting scattered light from a thus-illuminated region in the fluid. A calculation arrangement (2) determines a plurality of properties of the scattered light, and operates a calculation function on the plurality of properties to estimate the concentration of white blood cells in the fluid. The calculation arrangement may additionally determine and use a property of the light of the light beam that is transmitted by the fluid and / or a property of light that is scattered by the fluid from a second light beam with a wavelength in the range of 600-1000 nm.

Description

[0001] OPTICAL DETERMINATION OF WHITE BLOOD CELL CONCENTRATION IN FLUIDS

[0002] Technical Field

[0003] The present disclosure relates generally to optical techniques for measuring density of particles in fluids, and in particular the concentration of white blood cells in effluent from peritoneal dialysis.

[0004] Background Art

[0005] Inflammation in the peritoneum is common among patients undergoing peritoneal dialysis (PD). Early detection of infection or inflammation ("peritonitis") is essential to avoid suffering and therapy drop-out. Basically, there are two modalities for carrying out PD: automated peritoneal dialysis (APD) and a manual non-automated procedure denoted continuous ambulatory peritoneal dialysis (CAPD). In CAPD, infection may be detected by visual inspection of effluent bags in which spent dialysis fluid ("effluent") is collected. A cloudy effluent bag is a sign of peritonitis. The cloudiness is caused by increased presence of white blood cells (WBCs) caused by the infection. In APD, the effluent is often passed through an effluent line directly to the drain, and no visual inspection is possible. The infection is therefore detected late, when other signs such as stomach pain appear, and the peritoneum may be damaged.

[0006] It may also be relevant to detect presence of red blood cells (RBCs) in the effluent, since RBCs is a sign of bleeding within the peritoneal cavity (hemoperitoneum). The RBCs may originate from the peritoneal membrane, from the intraperitoneal organs (organs that are fully encapsulated by the visceral peritoneal membrane), or from partially or completely extraperitoneal structures. For women, the RBCs may be caused by normal ovulation or menstruation. If this origin can be ruled out, the RBCs may be a sign of issues with the PD catheter, or a more acute medical condition such as splenic laceration, liver rupture, Ever or renal cyst rupture, erosion of mesenteric vessel, bleeding from a malignant tumor, etc.

[0007] The prior art comprises W02022 / 008213, which discloses a technique of illuminating a fluid with a light beam and detecting scattered and / or transmitted light, where the particle density in the fluid is given by the temporal variability of the scattered and / or transmitted light. The technique is useful for determining the concentration of WBCs in PD effluent, presuming that the particles are WBCs.

[0008] WO2019 / 118929 proposes a technique of illuminating PD effluent with a light beam and detecting scattered and / or transmitted light, where the particle density in the PD effluent is given by the magnitude of the scattered and / or transmitted light. It is also proposed to use two light beams for the illumination to discriminate between RBCs and WBCs; a first light beam in the infrared (IR), and a second light beam in the range of 260-550 nm. A magnitude value of the scattered and / or transmitted light is determined for the first beam and second beam, respectively. The proposed evaluation technique results in a concentration value of WBCs or RBCs in the PD effluent and uses a first correlation plot that relates magnitude values for the first beam to concentration of RBCs and WBCs, respectively, and a second correlation plot that relates the ratio of magnitude values for the first and second beams to the concentration of RBCs and WBCs, respectively. As far can be understood, the proposed evaluation technique presumes that the effluent contains either WBCs or RBCs.

[0009] Summary

[0010] It is an objective to at least partly overcome one or more limitations of the prior art.

[0011] One objective is to provide an optical technique for estimating the concentration of white blood cells (WBCs) in a fluid.

[0012] Another objective is to provide such a technique that is capable of estimating the concentration of WBCs in the presence of red blood cells (RBCs).

[0013] A further objective is to provide such a technique that is capable of estimating the concentration of RBCs in the presence of WBCs.

[0014] One or more of these objectives, as well as further objectives that may appear from the description below, are at least partly achieved by optical detection apparatuses, a calculation arrangement, an apparatus for automated peritoneal dialysis, computer- implemented methods, and a computer-readable medium according to the independent claims, embodiments thereof being defined by the dependent claims.

[0015] The present disclosure emanates from a significant experimental effort by the Applicant to understand the dependence of scattered and transmitted light from a fluid, when illuminated by light beams at different wavelengths, on the concentration of WBCs and RBCs in the fluid. Surprisingly, the Applicant has found that the concentration of WBCs in a fluid can be estimated, even if the fluid contains RBCs, by use of a first time -dependent signal from a detector arranged to receive scattered light from the fluid when illuminated by a light beam with a wavelength in the range of 350- 575 nm. According to the Applicant's findings, it is possible to correctly estimate the concentration of WBCs based on two or more properties that represent the first timedependent signal. In this context, a correct estimation of WBC concentration lies within ±20% of the ground truth or ±100 cells / pL, whichever is the largest. The Applicant has also identified further time -dependent signals that may be used to improve the estimation.

[0016] The Applicant has also found that it is possible to estimate the concentration of WBCs in a fluid as a function of at least one property that represents scattered light from the fluid when illuminated by the light beam with a wavelength in the range of 350-575 nm, and at least one property that represents light transmitted through the fluid when illuminated by the light beam with a wavelength in the range of 350-575 nm.

[0017] Further, the Applicant has found that the combined concentration of WBCs and RBCs in a fluid can be estimated, by use of the first time -dependent signal in combination with a second time-dependent signal from a detector arranged to receive scattered light from the fluid when illuminated by a light beam with a wavelength in the range of 600-1000 nm. When the combined concentration has been estimated, the concentration of RBCs in the fluid can be estimated by accounting for the estimated concentration of WBCs.

[0018] Still other objectives as well as embodiments, features, advantages and technical effects may appear from the following detailed description, from the attached claims as well as from the drawings.

[0019] Brief Description of the Drawings

[0020] FIG. 1A is a side view of an example system for peritoneal dialysis comprising an optical detection apparatus, and FIG. IB is a side view, partly in section, of an example optical detection apparatus.

[0021] FIGS 2A-2C are side views of example arrangements for detection of transmitted and scattered light in relation to a fluid.

[0022] FIGS 3A-3B are graphs of an example signal representing scattered light after normalization.

[0023] FIGS 4A-4B are graphs of an example signal representing transmitted light after normalization.

[0024] FIGS 5A-5B are histograms corresponding to the signals in FIG. 3A and FIG. 4A.

[0025] FIG. 6 is a graph of optical absorption spectra for hemoglobin.

[0026] FIG. 7 is a block diagram of an example calculation device for determining densities of white blood cells (WBCs) and red blood cells (RBCs) in a fluid.

[0027] FIG. 8A is a flow chart of an example method of determining the concentration of WBCs in a fluid, and FIG. 8B is a flow chart of an example method of determining the concentration of RBCs in the fluid.

[0028] FIGS 9A-9B are block diagrams of example calculation devices for determining the concentration of WBCs and the total concentration of cells, respectively, in a fluid. FIG. 10A shows the average (top) and the standard deviation (bottom) of detected light from a WBC-containing fluid upon illumination by a laser beam at 650 nm as a function of detection angle, and FIG. 10B shows the average (top) and the standard deviation (bottom) of detected light from an RBC-containing fluid upon illumination by a laser beam at 650 nm as a function of detection angle.

[0029] FIG. 11A shows the average (top) and the standard deviation (bottom) of detected light from a WBC-containing fluid upon illumination by a laser beam at 405 nm as a function of detection angle, and FIG. 1 IB shows the average (top) and the standard deviation (bottom) of detected light from an RBC-containing fluid upon illumination by a laser beam at 405 nm as a function of detection angle.

[0030] FIG. 12 is a graph of an estimation error as a function of the number of properties used to estimate total cell concentration in accordance with an example.

[0031] FIGS 13-14 are graphs of estimated WBC concentration versus actual WBC concentration for a first estimation example and a second estimation example, respectively.

[0032] FIGS 15-16 are graphs of estimated total cell concentration versus actual total cell concentration for a first estimation example and a second estimation example, respectively.

[0033] FIGS 17-18 are graphs of estimated WBC cell concentration versus actual WBC cell concentration for further estimation examples.

[0034] FIGS 19A-19B are flow charts of alternative methods of determining the concentration of WBCs in a fluid.

[0035] Detailed Description of Example Embodiments

[0036] Embodiments will now be described more fully hereinafter with reference to the accompanying drawings, in which some, but not all, embodiments are shown. Indeed, the subject of the present disclosure may be embodied in many different forms and should not be construed as limited to the embodiments set forth herein; rather, these embodiments are provided so that this disclosure may satisfy applicable legal requirements.

[0037] Also, it will be understood that, where possible, any of the advantages, features, functions, devices, and / or operational aspects of any of the embodiments described and / or contemplated herein may be included in any of the other embodiments described and / or contemplated herein, and / or vice versa. In addition, where possible, any terms expressed in the singular form herein are meant to also include the plural form and / or vice versa, unless explicitly stated otherwise. As used herein, "at least one" shall mean "one or more" and these phrases are intended to be interchangeable. Accordingly, the terms "a" and / or "an" shall mean "at least one" or "one or more", even though the phrase "one or more" or "at least one" is also used herein. As used herein, except where the context requires otherwise owing to express language or necessary implication, the word "comprise" or variations such as "comprises" or "comprising" is used in an inclusive sense, that is, to specify the presence of the stated features but not to preclude the presence or addition of further features in various embodiments. The term "and / or" includes any and all combinations of one or more of the associated listed items.

[0038] It will furthermore be understood that, although the terms first, second, etc. may be used herein to describe various elements, these elements should not be limited by these terms. These terms are only used to distinguish one element from another. For example, a first element could be termed a second element, and, similarly, a second element could be termed a first element, without departing the scope of the present disclosure.

[0039] Well-known functions or constructions may not be described in detail for brevity and / or clarity. Unless otherwise defined, all terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this disclosure belongs.

[0040] Like reference signs refer to like elements throughout.

[0041] Embodiments will be described with reference to deployment in conjunction with automated peritoneal dialysis ("APD") therapy. In peritoneal dialysis (PD), dialysis fluid is infused into a patient's peritoneal cavity. This cavity is lined by the peritoneal membrane ("peritoneum") which is highly vascularized. Substances are removed from the patient's blood by diffusion across the peritoneum into the dialysis fluid. Excess fluid (water) is also removed by osmosis induced by the dialysis fluid being hypertonic. Automated peritoneal dialysis ("APD") is performed by an APD machine, commonly known as a "cycler". The cycler is operable to automatically perform one or more treatment cycles including fill, dwell and drain phases, for example while the patient sleeps. The cycler is fluidly connected to an implanted catheter, to a source of dialysis fluid and to a fluid drain. The cycler is operated to pump fresh dialysis fluid from the source, through the catheter, into the patient's peritoneal cavity and to allow the dialysis fluid to dwell within the cavity for the transfer of waste, toxins and excess water to take place. The cycler is then operated to pump spent dialysis fluid from the peritoneal cavity, through the catheter, to the drain. Spent dialysis fluid is commonly known as "effluent". The drain may be a fixed installation such as a sink, bathtub, toilet, floor drain, etc. Alternatively, the drain may be a reservoir or container for collecting fluid.

[0042] FIG. 1 is a schematic elevated side view of an arrangement 11 for peritoneal dialysis. The arrangement 11 comprises an APD machine Ila ("cycler"). The arrangement further comprises a disposable unit 11b mounted onto the cycler Ila. The cycler Ila comprises a combination of a control system, sensors and actuators to properly move fluid inside a hydraulic circuit of the disposable unit 11b. Although not shown in FIG. 1, the cycler Ila may also comprise a user interface for input / output of data. The disposable unit 11b comprises a cassette 12, as well as a set of tubes ("tubing set") connected to the cassette 12. In the illustrated example, the tubing set includes a container line 14 terminating in a container 13 configured for holding a treatment fluid ("dialysis fluid"). The container 13 may be in the form of a collapsible bag to be positioned, for example, in a dedicated tray on the cycler Ila. The container 13 may be delivered as a ready-made bag of dialysis fluid, or the container 13 may be filled by dialysis fluid prepared on-line by the cycler Ila or a separate machine (not shown), for example by mixing one or more concentrates with water. The cycler Ila may comprise a heater (not shown) for heating the treatment fluid before it is supplied to the patient. The tubing set further comprises a patient line 15 for connection to a catheter (not shown) implanted in the patient, and a drain line 16 for dispensing spent treatment fluid ("effluent") to a drain (not shown), for example a container or a sink. During operation, the arrangement 11 performs a fill phase in which a pumping mechanism in the cycler Ila actuates the cassette 12 to pump treatment fluid from the container 13 to the patient through the lines 14, 15, a dwell phase in which the treatment fluid is left within the peritoneal cavity of the patient, and a drain phase in which spent treatment fluid (effluent) is pumped to the drain through the patient and drain lines 15, 16.

[0043] Patients on PD are exposed to an elevated risk of attracting infection or inflammation in the peritoneal cavity, caused by bacteria admitted via the indwelling catheter. Often such infection or inflammation is located at the peritoneum and is denoted "peritonitis". Developed peritonitis may be manifested by the patient experiencing fever, diffuse abdominal pain, and nausea. Peritonitis represents a medical emergency, and early detection and treatment is essential to reduce morbidity and mortality in PD patients. In addition, repeated episodes of peritonitis may contribute to vascular proliferation and interstitial fibrosis, with ensuing loss of ultrafiltration capacity and therapy failure. In PD, peritonitis may be detected by extracting and analyzing the concentration of white blood cells (WBCs) in the effluent. According to established practice, a WBC concentration above 100 cells / pL in the effluent of PD is regarded as a sign of peritonitis and may result in the patient being given antibiotics.

[0044] The present disclosure relates to techniques that enable early detection of peritonitis in patients undergoing PD without the need to extract and analyze samples of the effluent. This is achieved by use of an optical detection apparatus (OPA) for mounting on the drain line. The OPA is thereby operable to produce a signal representative of the WBC concentration in the effluent that flows through the drain line.

[0045] In FIG. 1A, the OPA is represented by reference numeral 20 and comprises a measurement device 21, which is mounted onto the drain line 16, and a calculation device or arrangement 22 which is configured to receive and process one or more output signals of the measurement device 21 for determination of a value representative of the WBC concentration. It is conceivable that the measurement and calculation devices 21, 22 are integrated in a single unit. In a variant, the measurement device 21 is configured to wirelessly transmit its output signal(s) to the calculation device 22, which thus may be located remotely from the measurement device 21. In a further variant, the measurement device 21 and / or the calculation device 22 is included in the cycler I la. It is also conceivable that the measurement device 21 is mounted onto the patient line 15.

[0046] In some embodiments, the measurement device 21 is mounted onto a "tubing portion" which is transparent or at least translucent. The tubing portion may be included in the drain line 16, as shown in FIG. 1, or in the patient line 15. As used herein, a "transparent" material has the property of transmitting light with little scattering so that objects are clearly visible through the material, and a "translucent" material has the property of permitting passage of light while diffusing it so that objects are not clearly visible through the material. Commonly, the tubings in the disposable 11b are made of transparent or translucent material. The tubing portion may thus be an integral part of the drain line 16 or patient line 15. However, should the line 16 / 15 be made of opaque material, the tubing portion may be a separate tubing of transparent or translucent material that is spliced into the line 16 / 15. However, it is also conceivable that the tubing portion is a specialized cavity or chamber that is installed in the line 16 / 15 to facilitate the measurements by measurement device 21. Such a specialized cavity may be or comprise a cuvette with planar optical surfaces of transparent material.

[0047] In some embodiments, the OPA 20 is arranged to communicate with the cycler Ila. For example, the OPA 20 may receive, from the cycler Ila, a signal that indicates start of a drain phase and thereby presence of effluent in the drain line 16. Alternatively or additionally, the signal from the cycler I la may trigger the OPA 20 to measure the WBC concentration. It is also conceivable that the OPA 20 transmits a signal indicative of the measurement result to the cycler I la, which may operate a feedback unit to inform a user thereof.

[0048] In some embodiments, the OPA 20 is configured to illuminate a sample of effluent by light at one or more wavelengths, detect light that is scattered by the sample and / or light that is transmitted through the sample as a result of the illumination, and analyze the detected light for calculation of an output value that represents the concentration of WBCs in the effluent. For practical use, it is desirable for the OPA 20 to properly estimate the concentration of WBCs even if the effluent comprises other types of particles. One type of particle that is seen as problematic in this context is red blood cells (RBCs), which are known to interact with light by both absorption and scattering. As described in the Background section, RBCs may be present in the effluent for natural reasons (menstruation) or as a result of a complication or illness.

[0049] After significant experimentation and testing, the present Applicant has developed a technique of detecting an optical response to illumination that is specific to WBCs, irrespective of RBCs. The technique may be applied to determine the WBC concentration in effluent from PD therapy, in accordance with a method to be described below with reference to FIG. 8A. The technique enables early detection of peritonitis in patients undergoing PD and reduces the risk for serious complications.

[0050] The technique may be extended to estimate the RBC concentration in the effluent, even if the effluent contains WBCs. The extended technique enables early detection of blood in the effluent, allowing the caretaker to investigate the origin of the blood. If the origin is the catheter, an early remedy reduces the risk for complications. If the origin is a serious health issue, early detection may be lifesaving. The extended technique is described below with reference to FIG. 8B .

[0051] To facilitate understanding of the developed technique, a description will first be given of example detection arrangements that have been used by the Applicant when developing the technique. A description is also given of available measurement signals, as well as experimental results on which the technique is at least partly based.

[0052] FIG. IB is a side view, partly in section, of an example measurement device 21 arranged on a tubing portion 17, which may be part of the patient line 15 or the drain line 16. The measurement device 21 includes an illumination system ("light emitting arrangement") 30 and a detection system ("light detection arrangement") 31. In the illustrated example, the measurement device 21 also includes a holder 33, which is configured to engage with the tubing portion 17 to achieve a proper alignment of the tubing portion 17 with the illumination system 30 and the detection system 31. In FIG. IB, the holder 33 comprises first and second walls 37, 38 which define a slot for receiving the tubing portion 17. The spacing between the walls 37, 38 is equal to or smaller than the outer diameter of the tubing portion 17, so that the tubing portion 17 is squeezed, Fictionally held or otherwise fixed in the holder 33. At least part of the respective wall 37, 38 is transparent or translucent. In the example of FIG. IB, the respective wall 37, 38 defines a through-hole or opening 37a, 38a. In some embodiments, transparent / translucent windows may be arranged in the openings 37a, 38a. The illumination system 30 comprises at least one light source 34, which is configured to emit a light beam 300. The light source 34 may comprise a light-emitting diode (LED) or a miniaturized laser device such as a diode laser. In some embodiments, the illumination system 30 further comprises beam-forming optics 35, which may be arranged to focus the light beam 300 inside the holder 30, for example at a nominal location halfway between the walls 37, 38. Alternatively or additionally, the beamforming optics 35 may be configured to achieve a predefined transverse beam profile of the light beam 300. The beam-forming optics 35 may comprise one or more lenses. The light emitting device 34 is aligned with the walls 37, 38 so that the light beam 300 passes the openings 37a, 38a. The light beam 300 thereby defines a target volume inside the tubing portion 17. The light source 34 may be configured to generate time- continuous light or pulsed light.

[0053] The detection system 31 comprises light-detection devices ("detectors") 36a, 36b, which are responsive to the emitted light and provide a respective output signal OS1, OS2 that represents the amount of incident light ("intensity") on the respective detector 36a, 36b. The detector 36a is arranged to detect scattered light and is offset transversely to the light beam 300. The detector 36b is arranged to detect transmitted light and is thus aligned with the light beam 300. In some embodiments (not shown), the detection system 31 comprises detection optics to direct incoming light onto the respective detector 36a, 36b. The output signals OS1, OS2 are time -varying and comprises signal values that represent the momentary amount of incident light at different times. The detectors 36a, 36b may be separate units, such as photodiodes, photoresistors, phototransistors, etc. Alternatively, the detectors 36a, 36b may be formed by different portions of a light-sensing device, for example a photodiode, an array sensor, etc.

[0054] As indicated in FIG. IB, the light beam 300 interacts with particles 100 contained in the effluent within the tubing portion 17. In this interaction, the respective particle 100 may absorb and / or scatter photons in the light beam 300. Scattering is represented by dashed arrows in FIG. IB.

[0055] The effluent may be pumped through the tubing portion 17 at a flow rate F. In some embodiments, the detection system 31 is operated to detect light during a sequence of detection periods, where each detection period results in a signal value in the output signal. The minimum time between starts of detection periods may be set in relation to the expected or actual flow rate of the fluid through the tubing portion 17. A detection period may be achieved by selectively activating the respective detector 36a, 36b to be responsive to light and / or by selectively opening a shutter (not shown) in front of the detector 36a, 36b. Additionally or alternatively, detection periods may be achieved by pulsing the light beam 300. In a variant, the effluent is stationary in the tubing portion 17 during the measurement.

[0056] The use of pulsed light allows for the impact of ambient light on the measurement to be suppressed, if the detection system 31 is operated to detect light during and between light pulses, respectively. Thereby, the light detected between light pulses represents ambient light and may be subtracted from the light detected during light pulses to substantially remove the influence of ambient light. Such subtraction may be performed by the detection system 31 or the calculation device 22.

[0057] The overall operation of the measurement device 21 is controlled by a control unit 40, which may be configured to generate control signals Cl, C2 for the illumination system 30 and the detection system 31. The control signals Cl, C2 may control activation of the light source(s) 34 and the detectors 36a, 36b, respectively, as well as any shutter, if present.

[0058] FIG. 2A is a side view, partly in section, of an example arrangement on the detection side of an OPA. Although the tubing portion 17 is shown with a circular cross-section, it could have shape, for example rectangular. If the tubing portion 17 is deformable, the tubing portion may be deformed to define approximately planar and parallel wall portions facing the illumination system 30 and the detection system 31, for example as described in aforesaid W02022 / 008213.

[0059] As seen in FIG. 2A, the light beam 300 has a main direction 310 towards the tubing portion 17 and defines a target volume 320 within the effluent in the tubing portion 17. The target volume 320 is the volume from which the detection system 31 receives the main portion of the scattered light as the light beam 300 interacts with particles in the effluent. In the illustrated example, the light beam 300 is focused at the target volume 320 to increase the radiance of the light that interacts with the particles, to thereby increase the amount of scattered light. In alternative configurations, the light beam 300 may be focused elsewhere along the main direction 310 or be non-focused, for example collimated. The detector 36b is aligned with the main direction 310 to receive at least part of the light that passes the effluent ("transmitted light"). The detector 36a is arranged to receive scattered light at a predefined detection angle a to the main direction 310 and provide an output signal OS1. As used herein, a detection angle is given by the angle between the main direction 310 and the center line 311 of the scattered light that reaches the detector 36a, with the detection angle being measured at the center of the target volume 320. Since the light-sensing area of the detector 36a has an extent, the detector 36a will typically receive scattered light within a detection cone 301 ("angular range") that has an angular width Aa. The angular width of the detection cone 301 may differ in different directions depending on the shape of the light-sensing area. As used herein, the angular width refers to the largest dimension of the detection cone. Generally, the scattered light is symmetrically distributed around the main direction 310. In some embodiments, the detector 36a is ring-shaped to maximize the amount of scattered light that is detected. The width of the ring defines the angular width Aa. For example, the detector 36a may be shaped as a circular ring with uniform width. As shown, the detector 36b may be configured to detect transmitted light within a detection cone 300' with angular width Ap and provide an output signal OS2. As will be described further below, angular specificity in light detection may serve to improve the accuracy of the concentration values that are calculated by the calculation device 22.

[0060] The calculation device 22 is coupled to the detectors 36a, 36b to receive the output signals OS1, OS2 and output a result signal RS with concentration values. The concentration values may designate the concentration of WBCs, RBCs or both in the effluent.

[0061] The calculation device 22 may be implemented by hardware or a combination of software and hardware. In the example of FIG. 2A, the calculation device 22 comprises processor circuitry 22a and memory 22b. The processor circuitry 22a may e.g., include one or more of a CPU ("Central Processing Unit"), a DSP ("Digital Signal Processor"), a microprocessor, a microcontroller, an ASIC ("Application-Specific Integrated Circuit"), a combination of discrete analog and / or digital components, or some other programmable logical device, such as an FPGA ("Field Programmable Gate Array"). The memory 22b may include any form of conventional computer memory. The calculation device 22 may be at least partly operated in accordance with a control program comprising computer instructions. The control program is stored in the memory 22b and executed by the processor circuitry 22a. The control program may be supplied to the calculation device on a computer-readable medium, which may be a tangible (non-transitory) product (e.g., magnetic medium, optical disk, read-only memory, flash memory, etc.) or a propagating signal. Calculation devices in other embodiments described herein may be configured by analogy with the calculation device 22 in FIG. 2A.

[0062] FIG. 2B depicts an alternative configuration in which scattered light is detected at two spatially separate detection angles al and a2 by a first detector 36al and a second detector 36a2, which generate output signals OS 11, OS 12. Scattered light is thus detected by a detector 36a, which comprises detectors 36al, 36a2, within a first detection cone 301a ("first angular range") that has an angular width Aal and within a second detection cone 301b ("second angular range") that has an angular width Aa2. The detection angle al is defined between the main direction 310 and the center line 311a of the first detection cone 301a. The detection angle a2 is defined between the main direction 310 and the center line 311b of the second detection cone 301b. The first detection cone 301a is located closer to the main direction 310 than the second detection cone 301b. Thus, the output signal OS 11 represents "inner scattered light", and the output signal OS 12 represents "outer scattered light". As indicated, the output signals OS 11, OS 12 collectively correspond to the output signal OS1 in FIG. 2A.

[0063] In FIG. 2B, Eke in FIG. 2A, transmitted light is received by a detector 36b, which is aligned with the light beam 300 and configured to detect transmitted light within a detection cone 300' with angular width A . The output signal OS2 of the detector 36b represents the transmitted light.

[0064] As will be explained further below, it may be beneficial to illuminate the effluent by light beams at two different wavelengths. The illustration in FIG. 2B may be seen to represent illumination of the effluent by a first light beam 300 at a shorter wavelength located within a first wavelength band. FIG. 2C is identical to FIG. 2B and shows illumination of the effluent by a second light beam 400 at a longer wavelength located within a second wavelength band.

[0065] Each of the light beams 300, 400 have a confined spectral width. In the experiments presented herein, narrowband light is used for the illumination of the target regions 320, 420. In the context of the present disclosure, narrowband light has a spectral width of less than 20 nm, 10 nm or 5 nm, given as FWHM (full width at half maximum). To generate narrowband light, the light source 36 typically includes a narrowband LED or a laser, for example a semiconductor-based laser comprising one or more laser diodes. However, it is currently believed that the results and embodiments presented herein are equally applicable also when the respective light beam 300, 400 does not qualify as narrowband.

[0066] The Applicant has identified two wavelength bands with differing interaction between light and RBCs, on the one hand, and between light and WBCs, on the other hand. A first (lower) wavelength band extends from about 350 nm to about 575 nm. The first wavelength band may extend further into the ultraviolet (UV), for example to 200 nm, but such shorter wavelengths are currently not believed to be of practical use because of a lack of commercially available small-size light sources. Further, the tubing portion 17 may be made of plastics, which typically exhibit significant absorption of UV light. There is also significant light absorption by oxygen molecules below 200 nm. A second (higher) wavelength band extends from about 600 nm to about 1000 nm. The second wavelength band may extend further into the infrared, for example to 1500 nm.

[0067] Reverting to FIGS 2A-2C, the wavelength of the first light beam 300 is confined to the first wavelength band, and the wavelength of the second light beam 400 is confined to the second wavelength band. The selection of wavelength bands may be understood based on the absorption spectrum of RBCs. FIG. 6 shows curves 150, 151 of the molar extinction coefficient (MEC) as a function of wavelength for hemoglobin (Hb) and oxyhemoglobin (HbO2), respectively. MEC is a measure of how strongly a chemical species absorbs, and thereby attenuates, light at a given wavelength. Hemoglobin and oxyhemoglobin are carried by RBCs. As seen, absorption is significant at wavelengths below about 600 nm. Above about 600 nm, the absorption decreases rapidly, given that the scale is logarithmic in FIG. 6. Thus, in the first wavelength band Wl, the light is primarily absorbed by RBCs, if present, rather than being scattered. Conversely, in the second wavelength band W2, light is primarily scattered by RBCs, if present, rather than being absorbed. For WBCs, the absorption is less dependent on wavelength. Assuming light interacts with WBCs through Mie scattering, scattering will increase with decreasing wavelength. Thus, the interaction of RBCs and WBCs with light will differ between the first and second wavelength bands, which is seen in experiments conducted by the Applicant (below).

[0068] As seen in FIG. 6, absorption of RBCs is larger below about 450 nm, with a maximum at 410-430 nm. It may thus be advantageous for the wavelength of the first light beam to be below 450 nm, and in particular in the range of 410-430 nm. In some embodiments, the first light beam has a peak intensity in the wavelength range of 400- 425 nm, and the second light beam has a peak intensity in the wavelength range of 630- 655 nm. The experimental results presented herein are obtained for illuminating light at 405 nm and 650 nm, respectively.

[0069] Reverting to FIG. 2C, the second light beam 400 has a main direction 410 towards the tubing portion 17 and defines a target volume 420 within the effluent in the tubing portion 17. Scattered light is detected at two spatially separate detection angles al' and a2' by a first detector 36al' and a second detector 36a2', which generate output signals OSH', OS 12'. As indicated, the output signals OSH', OS 12' may been seen to collectively correspond to the output signal OS1'. Scattered light is thus detected by a detector 36a', which comprises detectors 36al', 36a2', within a first detection cone 401a ("first angular range") that has an angular width Aal' and within a second detection cone 401b ("second angular range") that has an angular width Aa2'. The detection angle al' is defined between the main direction 410 of the second light beam 400, and the center line 411a of the first detection cone 401a. The detection angle a2 is defined between the main direction 410 and the center line 411b of the second detection cone 401b . The first and second detection cones 401a, 401b may be dimensioned in correspondence with the first and second detection cones 301a, 301b. The output signal OSH' of the detector 36al' represents "inner scattered light", and the output signal OS 12' of the detector 36a2' represents "outer scattered light". Transmitted light is received by a detector 36b', which is aligned with the light beam 400 and configured to detect transmitted light within a detection cone 400' with angular width A0'. The output signal OS2' of the detector 36b' thus represents the transmitted light.

[0070] FIGS 2B-2C show the experimental set-up that has been used to explore, develop and refine the technique presented herein. In a practical implementation of the measurement device 21, the configuration may be different depending on the measurement data that is used for estimating the WBC concentration and / or the RBC concentration. For example, for the first light beam 300 and / or the second light beam 400, a single detection cone may be used, or detection of transmitted light may be omitted. Thus, the number of detectors may be reduced. It is also conceivable to use more than two detection cones for the first light beam 300 and / or the second light beam 400. Likewise, it is conceivable to use more than two light beams at different wavelengths.

[0071] It may also be noted that the first and second light beams 300, 400 may illuminate different portions of the effluent in the tubing portion 17. In other words, the target volumes 320, 420 may be spatially separated, for example shifted along the extent of the tubing portion 17. Such a configuration may be simpler to implement, at the cost of more equipment.

[0072] In an alternative configuration, the illumination system 30 is configured to direct the first and second light beams to illuminate approximately the same portion of the effluent, but at different time points. This means that the target volumes 320, 420 are very close to each other or effectively overlap. This makes it possible to use the same detectors for detecting scattered and / or transmitted light from the first and second light beams 300, 400.

[0073] FIG. 3A shows a pre-processed version of an example signal 130 that represents scattered light received by the detector 36al in FIG. 2B. The signal 130 is given by the signal OS 11 after processing for removal of a baseline signal by subtraction, to yield normalized intensity. The signal 130 is shown for a time period of 90 seconds. In this example, signal values are provided every 2 ms, which corresponds to the above- mentioned detection period. FIG. 3B is an enlarged view of FIG. 3 A and shows individual signal values more clearly.

[0074] One reason for the normalization of signals representing scattered light is to reduce the impact of light that originates from scattering of the light beam 300 by the tubing portion 17. The baseline signal may be given by the lower temporal envelope of the signal OS 11. The lower temporal envelope may be determined by any conventional signal processing technique, as readily available to the person skilled in the art, for example by extraction of selected values from the signal OS 11 or by operating a Hilbert transformer on the signal OS 11. In a non-limiting example, the lower temporal envelope is given by determining the minimum for a sliding window. Alternatively, the baseline signal may be given by a single value, which is calculated from the signal OS 11 and subtracted from all signal values in the signal OS 11. In a variant, the baseline signal is given by the upper temporal envelope of the signal OS2 (below), which is subtracted from the signal OS 11. These examples of normalization are applicable to any signal that represents scattered light.

[0075] In a further variant, the normalization is performed by dividing the signal representing scattered light (e.g., OS 11) by a baseline signal given by the signal representing transmitted light (e.g., OS2). For example, the baseline signal may represent the upper temporal envelope of the signal representing transmitted light. This specific type of normalization is denoted "cross-normalization" herein.

[0076] FIG. 4A shows a pre-processed version of an example signal 140 that represents transmitted light received by the detector 36b in FIG. 2B. The signal 140 is given by the signal OS2 after processing for removal of a baseline signal by division, to yield normalized intensity. The signal 140 is given for a time period of 90 seconds, with signal values being provided every 2 ms. FIG. 4B is an enlarged view of FIG. 4A and shows individual signal values more clearly.

[0077] One reason for the normalization of signals representing transmitted light is to reduce the impact of changes in the intensity of the light beam 300 over time. The baseline signal may be given by the upper temporal envelope of the signal OS2. The upper temporal envelope may be determined in correspondence with the lower temporal envelope. In a non-limiting example, the upper temporal envelope is given by determining the maximum for a sliding window. As an alternative or supplement to using a baseline signal for normalizing the signal OS2, the energy of the laser beam 300 may be measured by a light detector (not shown), for example in the illumination system 30 (FIG. IB), and used for normalizing the signal OS2, for example by division. These examples of normalization are applicable to any signal that represents transmitted light.

[0078] Generally, pre-processing by normalization of the signals from the detectors in the measurement device 21 may improve the accuracy and robustness of the calculated concentration values. However, the data analysis presented further below indicates that acceptable accuracy and robustness is possible without normalization. Further, the data analysis indicates that it may be preferable to use both normalized signals and nonnormalized signals in the calculation of concentration values. FIG. 5A is a histogram showing the distribution of normalized signal values in the signal 130 in FIG. 3A. Similarly, FIG. 5B is a histogram showing the distribution of normalized signal values in the signal 140 in FIG. 4A. The Applicant has found that, by proper choice of wavelength of the first light beam 300 and the second light beam 400, the distribution of signal values that are acquired during a measurement time period (MTP) changes in dependence of the WBC concentration and / or RBC concentration. This is especially significant for scattered light, but may also be visible in transmitted light.

[0079] The MTP may be set to provide a sufficient number of signal values to represent the distribution. In the examples given herein, the distribution is analyzed based on approximately 45 000 signal values, corresponding to MTP being 90 seconds. It is currently believed that the MTP should result in at least 5 000 signal values, and preferably at least 10000 signal values. Alternatively or additionally, the MTP may be set to be less than about 200, 150, or 100 seconds, and larger than about 5, 10 or 15 seconds. Alternatively or additionally, the MTP may be set in view of the flow rate of the effluent (cf. F in FIG. IB) and / or a desired accuracy of the calculated concentration value(s).

[0080] The Applicant has chosen to characterize the distribution of signal values by two main categories: magnitude and variability. Both of these categories are thus estimated based on an ensemble of signal values obtained during the MTP. The variability represents the variation over time ("temporal variability"). In the following, the distribution of normalized signal values is denoted "normalized distribution", by contrast to the distribution of non-normalized signal values, which is denoted "original distribution".

[0081] For normalized distributions, for example shown in FIGS 5A-5B, the magnitude may be given by the average or mean of the included signal values ("normalized average"), or by any equivalent measure, such as median or norm. By correlation analysis, the Applicant has found that the k:th percentiles of the original distribution correlates to some degree with the normalized average, at least for k > 5. The term k:th percentile is used in its ordinary meaning to denote a score at or below which a k:th percentage of the signal values falls. The 50:th percentile is equal to the median. Likewise, various quantiles of the original distribution correlate to some degree with the normalized average, including the inter-quartile range (IQR). Further, the average or mean, or any equivalent measure, of the signal values in the original distribution correlates with the normalized average. Thus, any one of the foregoing measures of the original distribution corresponds to the normalized average and may thus be used to represent the magnitude of the signal values obtained during an MTP. For normalized distributions, the variability may be given by the variance of the included signal values ("normalized variability"), or by any equivalent measure, such as any of the variability measures described in aforesaid W02022 / 008213, including but not limited to energy, standard deviation, coefficient of variation, variance-to-mean, or Median Absolute Deviation or Mean Absolute Deviation (MAD). By correlation analysis, the Applicant has found that the k:th percentiles of the normalized distribution correlates to some degree with the normalized variability for k > 1. Likewise, various quantiles of the normalized distribution correlate to some degree with the normalized variability, including the inter-quartile range (IQR). Further, the variance, or any equivalent measure, of the signal values in the original distribution correlates with the normalized variability. Thus, any of the foregoing measures of the original distribution corresponds to the normalized variability and may thus be used to represent the variability of the signal values obtained during an MTP.

[0082] As used herein, the term "magnitude" designates an average parameter that is calculated for the ensemble of signal values within the MTF after removal of a baseline signal, typically by subtraction for scattered light or by division for transmitted light, or an alternative parameter that is correlated with the average parameter and calculated for the signal values within the MTF without removal of the baseline signal. Further, as used herein, the term "variability" designates a variability parameter that is calculated for the ensemble of signal values within the MTF after removal of a baseline signal, typically by subtraction for scattered light or by division for transmitted light, or an alternative parameter that is correlated with the variability parameter and calculated for the signal values within the MTF without removal of the baseline signal.

[0083] The present Applicant has conducted experiments to analyze the impact of presence of WBCs and RBCs, respectively, on the magnitude and the variability of transmitted and scattered light for a beam of narrowband laser light with a wavelength in the first wavelength band and the second wavelength band, respectively. In these experiments, the detectors 36al, 36a2, 36b in FIG. 2B and the detectors 36al', 36a2' and 36b' in FIG. 2C were replaced with a two-dimensional CCD detector, and the light received by the CCD detector was separated by detection angle to the main direction 310. A few examples of the experimental results are shown in FIGS 10-11.

[0084] For the experiments, samples with different concentrations of RBCs were prepared by adding RBCs to phosphate buffered saline (PBS). Similarly, samples with different concentrations of WBCs were prepared by adding WBCs to PBS. FIGS 10-11 are given for RBC concentrations of 0 cells / pL, 63 cells / pL and 1000 cells / pL, and WBC concentrations of 0 cells / pL, 63 cells / pL and 1000 cells / pL. In FIGS 10-11, results for concentrations of 0, 63 and 1000 cells / pL are indicated by solid lines, dashed lines and dot-dashed lines, respectively. The results are presented for a laser beam at 405 nm (first wavelength band, Wl) and a laser beam at 650 nm (second wavelength band). The laser beam at 405 nm will be referred to as "blue light", and the laser beam at 650 nm will be referred to as "red light". Each of FIGS 10A, 10B, 11 A, 11B includes a top graph that shows magnitude as a function of detection angle, and a bottom graph that shows variability as a function of detection angle. The magnitude is given by the average of the detected light at the respective angle within the MTP, and the variability is given by the standard deviation of the detected light at the respective angle within the MTP.

[0085] FIG. 10A shows experimental results for illumination of WBC samples with the laser beam at 650 nm (red light). The top graph, in encircled region 160, indicates that the magnitude of scattered light changes with WBC concentration, at least for angles in the range of 7°-30°. The bottom graph, in encircled region 161, indicates that the variability of the transmitted light changes with WBC concentration.

[0086] FIG. 10B shows experimental results for illumination of RBC samples with the laser beam at 650 nm (red light). The top graph, in encircled region 162, indicates that the magnitude of scattered light changes with RBC. The bottom graph, in encircled region 163, indicates that the variability of the transmitted light changes with RBC concentration. For red light, the response to RBCs (regions 162, 163) is thus similar to the response to WBCs (regions 160, 161).

[0087] FIG. 11A shows experimental results for illumination of WBC samples with the laser beam at 405 nm (blue light). The top graph, in encircled region 164, indicates that the magnitude of scattered light changes with WBC concentration, at least for angles in the range of 7°-30°. The bottom graph, in encircled region 165, indicates that the variability of scattered light changes with WBC concentration, at least for angles in the range of 7°-30°. The bottom graph further indicates, in encircled region 166, that the variability of transmitted light changes with WBC concentration.

[0088] FIG. 11B shows experimental results for illumination of RBC samples with the laser beam at 405 nm (blue light). In the top graph, the magnitude of the detected light is effectively independent of RBC concentration. The bottom graph, in encircled region 167, indicates that the variability of transmitted light changes with RBC concentration.

[0089] The results in FIGS 10-11 indicate that blue light is the best candidate for determining the concentration of WBCs in the presence of RBCs since the both the magnitude and the variability of the scattered light seem to be dependent on WBC concentration (FIG. 11 A) and almost independent of, or at least much less dependent, on the RBC concentration (FIG. 11B). The present Applicant has found that it might be easier to estimate, based on measurement data for scattered and / or absorbed light, the total concentration of particles and then determine the RBC concentration from the total concentration and WBC concentration. In some embodiments, the RBC concentration is given by the difference between the total concentration and the WBC concentration, assuming the contribution to the measurement data from other particles is small. Experiments indicate that the total cell concentration can be determined with relatively high accuracy based on at partly the same signals that are used for determining the WBC concentration.

[0090] HG. 7 is a block diagram of an example calculation device 22 in accordance with an embodiment. In the illustrated example, the calculation device 22 is configured to receive output signals OS1, OST and OS2 from three detectors in the measurement device 21 (not shown). With reference to FIG. 2 A, signal OS1 represents scattered light of wavelength in band W1 ("blue light") and is obtained from detector 36a, signal OST represents scattered light of wavelength in band W2 ("red light") and is obtained from a detector not shown in HG. 2A (cf. detector 36a' in HG. 2C), and signal OS2 represents transmitted light of wavelength W1 ("blue light") and is obtained from detector 36b. The calculation device 22 comprises a first parameter calculation unit 23, which is configured to calculate values of a first set of predefined parameters ("properties") based on the signals OS1, OST, OS2. Each parameter or property represents one of the above-mentioned categories, namely magnitude or variability, or a combination thereof. The resulting set of property values, [P], is provided to a WBC calculation unit 24, which is configured to estimate the WBC concentration, C_WBC, based on [P], by use of a first calculation function Fl. The function Fl defines a predefined relation between C_WBC and [P] and may be given as a mathematical expression or a look-up table.

[0091] In the illustrated example, the calculation device 22 comprises a second parameter calculation unit 25, which is configured to calculate values of an additional set of predefined parameters ("properties") based on the signals OS1, OS1', OS2. Each parameter or property represents one of the above-mentioned categories, namely magnitude or variability, or a combination thereof. The resulting set of property values, [P*], is provided to a TC calculation unit 26, which is configured to estimate the total concentration of particles, C_TC, based on [P*], by use of a second calculation function F2. The function F2 defines a predefined relation between C_TC and [P*] and may be given as a mathematical expression or a look-up table.

[0092] The calculation device 22 further comprises an RBC calculation unit 27, which is configured to calculate the RBC concentration, C_RBC, based on C_TC and C_WBC. As noted above, C_RBC may be obtained by subtracting C_WBC from C_TC. In the illustrated example, the calculation device 22 is configured to output both C_WBC and C_RBC. With reference to FIG. 2A, C_WBC and C_RBC would be included in the result signal, RS. In a variant, only C_RBC is output. In a further variant, only C_WBC is output. If only C_WBC is calculated and output, units 25-27 are superfluous and may be omitted.

[0093] In a further variant, the total concentration of particles, C_TC, is obtained from another measurement device, which may determine C_TC by use of the optical technique described in aforesaid W02022 / 008213, turbidimetry, optical coherence tomography (OCT), direct imaging, a Coulter counter, or any other commercially available particle counter.

[0094] The calculation device 22 may further include a pre-processing unit, which is configured to pre-process one or more of the signals OS1, OST, OS2, for example to remove outlier data, to perform the above-mentioned normalization, to perform noise reduction, to perform a high-pass filtration, etc.

[0095] The separation into units 23-27 in FIG. 7 is done for illustrative purposes. In practice, the functionality of the calculation device 22 may be implemented by any number of units. It is also to be understood that the calculation device 22 may be modified to only use a sub-set of the signals OS1, OST and OS2 for the calculations, for example in accordance with the methods in FIGS 8A-8B.

[0096] FIG. 8 A is a flowchart of an example method Ml for operating an OPA 20 to determine the concentration of WBCs in a fluid. The fluid may be effluent from an APD cycler. The method Ml will be described with reference to the OPA 20 in FIGS 1-2. Optional steps are indicated by dashed fines.

[0097] In step S10, the illumination system 30 is operated to illuminate the fluid by light in the first wavelength band, W1 (FIG. 6). The illumination system 30 thereby generates the first light beam 300, which is directed to a first region in the fluid. The first region corresponds to the target volume 320 in FIG. 2B. In step S12, the detection system 31 is operated to detect scattered light from the first region and provide an output signal. Here, the scattered light that is detected in step S12 is denoted "first scattered light". Depending on configuration of the detection system 31, the output signal may be OS1 (FIG. 2A), one of OS 11, OS 12, or a combination of OS 11 and OS 12 (FIG. 2B). Steps S 10 and Sil may be performed by the control device 40, shown in FIG. IB.

[0098] In step SI 2, the calculation device 22 obtains the output signal from step Sil and determines, based on this signal, a plurality of first properties, [Pl], that represent the first scattered light. In step S 18, the calculation device 22 operates a first calculation function on [Pl] to determine the WBC concentration. The first calculation function corresponds to Fl in FIG. 7, in which the set of property values [P] corresponds to [Pl], optionally combined with [P2] and / or [P3] (below).

[0099] As used herein, a "property" is a characteristic of light received by a detector and is typically a characteristic of the histogram of signal values within a time window in the output signal from the detector (cf. FIGS 5A-5B). In other words, the property is determined to represent an ensemble or time series of signal values. The time window corresponds to the above-mentioned measurement time period, MTP. Typically, the property falls into one of two main categories: magnitude or variability. As noted above, both magnitude and variability may be calculated in many different ways. After significant experimentation, the present Applicant has made the surprising finding that the WBC concentration (C_WBC) may be estimated with adequate accuracy based on a plurality of first properties, [Pl], for example two first properties. As shown by step S12a, the properties in [Pl] may be determined to represent magnitude and / or variability. In some embodiments, [Pl] only represents magnitude and thus contains different measures of magnitude for an ensemble of signal values. In some embodiments, [Pl] only represents variability and thus contains different measures of variability for an ensemble of signal values. To be useful, these different measures should not be perfectly correlated with each other. It is currently believed that different measures of magnitude are useful if they have a correlation coefficient to each other in the range of 0-0.8. The same applies to different measures of variability. In some embodiments, [Pl] includes one or more properties that represent magnitude and one or more properties that represent variability. Generally, the use of both magnitude and variability has been found to improve the accuracy of the estimated C_WBC.

[0100] Based on the experiments described with reference to FIGS 10-11, the Applicant has achieved good results when the first scattered light is received at a detection angle a, which is in a range from about 6° to about 35° and defined in accordance with FIG. 2A. The detection cone of the detector 36a may have any suitable angular width, Aa (FIG. 2A), for example l°-10°. With reference to FIG. 2A, the detector 36a may be arranged to not receive scattered light at angles outside the range 6°-35° to maximize the relevance of the detected light.

[0101] The Applicant has identified a potential for further improvement by through angularly specific detection of the first scattered light. Surprisingly, it has been found that the first scattered light that is received at different detection angles a have different dependencies on the WBC concentration. Thus, in some embodiments, the first scattered light is detected in first and second angular ranges, which differ from each other. The first and second angular ranges corresponds to the detection cones 301a, 301b in FIG. 2B. The first and second angular ranges may partially overlap. However, to maximize the collected information about WBCs, it is believed that the first and second angular ranges should be non-overlapping. Overlapping or not, one angular range is located closer to the main direction 301 than the other. In the following, the first and second angular ranges are therefore denoted "inner angular range", and "outer angular range", respectively. In FIG. 2B, the inner angular range corresponds to detection cone 301a, and the outer angular range corresponds to detection cone 301b. The Applicant achieved good results with the inner angular range 301a being fully located within 7°-16° to the main direction 310, and the outer angular range 301b being fully located within 16-35° to the main direction 310. In the notation of FIG. 2B, the inner angular range 301a fulfils: al - al / 2 > 7° al + Acrl / 2 < 16° and the outer angular range 301b fulfils: a2 - a2 / 2 > 16° a2 + Aa2 / 2 < 35°

[0102] Returning to FIG. 8 A, the method Ml may include a step SI 2b of determining a first magnitude of the first scattered light in the inner angular range 301a, and a step S 12c of determining a second magnitude of the first scattered light in the outer angular range 301b. Data analysis indicates that the accuracy of step S18 may be improved if [Pl] represents the first and second magnitudes from steps S12b-S12c.

[0103] Further analysis of the experimental results indicates that it may be beneficial to include a property in [Pl] that represents the variability of the first scattered light detected in the outer angular range 301b. This property may or may not be combined with the first and / or second magnitudes from steps S12b-S 12c, or a magnitude determined by step S 12a. Thus, the method Ml may include a step S12d of determining the temporal variability of the first scattered light in the outer angular range 301b.

[0104] Experiments also indicate that the accuracy of step S 18 may be improved by steps S15-S17. In step S 15, the illumination system 30 is operated to illuminate the fluid by light in the second wavelength band, W2 (FIG. 6). The illumination system 30 thereby generates the second light beam 400, which is directed to a second region in the fluid. The second region corresponds to the target volume 420 in FIG. 2C. In step S 16, the detection system 31 is operated to detect scattered light from the second region and provide an output signal. Here, the scattered light that is detected in step S 16 is denoted "second scattered light". Depending on configuration of the detection system 31, the output signal may be the signal OS 1' , or one of OSH', OS 12', or a combination of OSH' and OS 12' (FIG. 2C). Steps S15 and S 16 may be performed by the control device 40, shown in FIG. IB.

[0105] In step S17, the calculation device 22 obtains the output signal from step S 16 and determines, based on this signal, at least one second property, [P2], that represents the second scattered light. In step S I 8, the calculation device 22 operates the first calculation function on [Pl] and [P2] to determine the WBC concentration. Data analysis indicates an improvement when [P2] includes a property that represents the variability of the second scattered light. Thus, the method Ml may include a step S 17a of determining the temporal variability of the second scattered light. It is currently believed, backed by experimental data, to be beneficial if [P2] includes the variability of the second scattered light in the inner angular range. Thus, step S17a may be replaced by a step S 17b of determining the temporal variability of the second scattered light in the inner angular range.

[0106] Experiments also indicate that the accuracy of step S 18 may be improved by steps S13-S14. In step S 13, the detection system 31 is operated to detect transmitted light from the first region and provide an output signal. Step S 13 is performed based on the first light beam (in the first wavelength band, Wl) that is generated by step S10. Step S13 may or may not be performed concurrently with step Si l. The output signal corresponds to OS 2 in FIGS 2A-2B. Step S 13 may be performed by the control device 40, shown in FIG. IB.

[0107] Based on the experiments described with reference to FIGS 10-11, the Applicant has found that good results are achieved when the transmitted light is detected within a detection cone that has a width A (FIGS 2A-2B) of less than 6° and preferably less than 5° or 4°, to limit the impact of the first scattered light on the measurement.

[0108] In step SI 4, the calculation device 22 obtains the output signal from step S 13 and determines, based on this signal, at least one third property, [P3], that represents the transmitted light. The transmitted light that is detected in step S 14 may also be denoted "first transmitted light". In step S 18, the calculation device 22 operates the first calculation function on [Pl] and [P3] to determine the WBC concentration. Data analysis indicates an improvement when [P3] includes a property that represents the magnitude of the first transmitted light. Thus, the method Ml may include a step SI 4a of determining the magnitude of the first transmitted light.

[0109] Data analysis also indicates that it may be beneficial to operate the first calculation function on [Pl], [P2] and [P3] to determine the WBC concentration in step S18. Data analysis further indicates that the estimate of the WBC concentration may be improved by using, in step SI 8, at least one additional property [P'J that represents the transmitted light from the second region ("second transmitted light"). Thus, although not shown in FIG. 8 A, the method Ml may include a step of operating the detection system 31 to detect second transmitted light from the second region and provide an output signal. This step is performed based on the second light beam (in the second wavelength band, W2) that is generated by step S15. The output signal corresponds to OS2' in FIG. 2C. Data analysis indicates a potential improvement when [P'] represents the variability of the second transmitted light. Thus, the method Ml may include a step of determining the variability of the second transmitted light from the second region. In step S I 8, [P'] may replace or supplement [P2] and / or [P3].

[0110] As shown by step S19, the method Ml may involve an evaluation of the WBC concentration from step SI 8, for example for detection of potential peritonitis. The evaluation may comprise comparing the WBC concentration to a threshold value. If an elevated WBC concentration is detected, the caretaker may be alerted thereof. The evaluation may be performed by the control system of the cycler 1 la (FIG. 1 A), and the caretaker may be alerted through a feedback unit. The feedback unit may comprise one or more of a display, an indicator lamp, a speaker, a buzzer, a beeper, a vibrator, etc. The feedback unit may be included in the cycler Ila, the OPA 20, or in a separate device. Alternatively or additionally, step S19 may involve presenting the WBC concentration to the caretaker through the feedback device. Alternatively or additionally, step S 19 may involve storing the WBC concentration for the patient, for example on a memory of the cycler 1 la or on a remote storage device, such as a server or a cloud service.

[0111] FIG. 9A is a block diagram of an example calculation device 22 which is configured to implement embodiments of the method Ml of FIG. 8A. The first parameter calculation unit 23 is configured to receive at least the signals OS 11 and OS 12, which represent the scattered light of the first light beam 300 in the inner and outer angular ranges 301a, 301b (FIG. 2B). As indicated by dashed arrows, the unit 23 may optionally be configured to receive the signal OSH', which represents the second scattered light in the inner angular range 401a (FIG. 2C), and / or the signal OS2, which represents the first transmitted light of the first light beam 300 (FIG. 2B). The unit 23 is configured to determine the plurality of first properties, [Pl], based on OS 11 and OS 12, in accordance with step S12 (FIG. 8A). As indicated by italic characters, the unit 23 may optionally be configured to determine the at least one second property, [P2], based on OSH', in accordance with step S17. Alternatively or additionally, the unit 23 may be configured to determine the at least one third property, [P3], based on OS2, in accordance with step S14. The WBC calculation unit 24 is configured to determine C_WBC by operating the first calculation function, Fl, on [Pl], optionally in combination with [P2] and / or [P3], in accordance with step S17. The calculation device 22 is configured to output C_WBC for further processing, for example in accordance with step S19.

[0112] The number of properties that are extracted and used in the first and second calculation functions Fl, F2 has been found to influence the accuracy of the calculated concentrations. FIG. 12 is a graph of prediction error (RMSE) as a function of the number of properties included in the second calculation function for a linear dependence on the properties (cf. Equation 3, below). The prediction error is determined based on the experiments and results described below with reference to FIGS 15-16. As seen, the prediction error decreases significantly when going from one to two properties, which is believed to be the minimum number of properties to be used. It is also seen that the prediction error levels out when more than three properties are included. Thus, based on FIG. 12, it seems as if using four or more properties does not add much in terms of prediction accuracy of the second calculation function. The same dependence of prediction error on the number of properties is seen for the first calculation function.

[0113] Experiments have been performed to determine the first calculation function, Fl. In these experiments, reference fluids were prepared with known concentrations of WBCs and RBCs. Specifically, the reference fluids had WBC concentrations in the range of 0-15000 cells / pL, and RBCs were included at concentration ratios of 0, 0.2, 0.5 and 0.8. The method Ml was performed for each of the reference fluids for different combinations of extracted properties, resulting in predicted values of WBC concentration. Some example results are presented in FIGS 13-14.

[0114] FIG. 13 is a graph of predicted WBC concentration versus known WBC concentration for a linear dependence on three specific properties. The solid line represents a one-to-one relation between predicted and known WBC concentrations. The high degree of correlation shows that WBC concentration in a fluid can be determined by use of the method Ml irrespective of RBC concentration in the fluid. In the example of FIG. 13, the predicted WBC concentration is determined based on [Pl] only, given by step S12 in FIG. 8A. In the illustrated example, [Pl] consists of a property Pl a representing the magnitude of first scattered light in the outer angular range 301b (FIG. 2B), a property Plb representing the variability of first scattered light in the outer angular range 301b, and a property Pic representing the magnitude of first scattered light in the inner angular range 301a. A linear correlation was used, represented as: y = a + b ■ Pla + c ■ Plb + d ■ Pic (1) with y being the predicted WBC concentration, and a-d being weights or coefficients given by the linear regression. The first calculation function Fl, used by step S 18, may be given by Equation (1).

[0115] In the specific example of FIG. 13, Pla is given by the 90:th percentile of the signal values within a time window (MTP = 90 seconds) in the signal OS 12 (FIG. 2B), Plb is given by the 75:th percentile of the signal values within the MTP in the signal OS 12 after normalization (by subtraction of a baseline signal), and Pic is given by the 90: th percentile of the signal values within the MTP in the signal OS 11 (FIG. 2B). The correlation of data in FIG. 13 yields R = 0.985, MAE = 143.3 and RMSE = 268.7, with R being the correlation coefficient, MAE being the mean absolute error, and RMSE being the root means square error.

[0116] FIG. 14 is a graph of predicted WBC concentration versus known WBC concentration when including a second order dependence on specific properties. FIG. 14 shows that WBC concentration in a fluid can be determined irrespective of RBC concentration. In the example of FIG. 14, the predicted WBC concentration is determined based on [Pl] and [P2], given by steps S12 and S17 in FIG. 8A. In the illustrated example, [Pl] consists of a property Pla representing the magnitude of first scattered light in the outer angular range 301b (FIG. 2B), a property Plb representing the variability of first scattered light in the outer angular range 301b, and a property P2a representing the variability of second scattered light in the inner angular range 401a (FIG. 2C). A linear correlation was used, represented as: with y being the predicted WBC concentration, and a-j being weights or coefficients given by the linear regression. The first calculation function Fl, used by step S 18, may be given by Equation (2). If one or more of the coefficients a-j is small, it may be set to zero (0). For example, in the analysis underlying FIG. 14, coefficients d and j are about 1 / 1000 and 1 / 100, respectively, of the size of the other coefficients.

[0117] In the specific example of FIG. 14, Pla is given by the 90:th percentile of the signal values within a time window (MTP = 90 seconds) in the signal OS 12 (FIG. 2B), Plb is given by the 75:th percentile of the signal values within the MTP in the signal OS 12 after normalization (by subtraction of a baseline signal), and P2a is given by the MAD of the signal values within the MTP in the signal OSH' after normalization (by subtraction of a baseline signal). The correlation of data in FIG. 14 yields R = 0.989, MAE = 62.9 and RMSE = 232.1.

[0118] FIG. 17 is another graph of predicted WBC concentration versus known WBC concentration when including a second order dependence on specific properties. In the example of FIG. 17, the calculation function Fl includes two properties Pla, Plb given by steps S 12 and S 14 in FIG. 8A. The calculation function is defined by y = (a + b ■ Pla + c ■ Plb~)2, where Pla represents the magnitude of first scattered light in the inner angular range 301a (FIG. 2B), and Plb represents the variability of first scattered light in the outer angular range 301b. In the specific example of FIG. 17, Pla is given by the mean of the signal values within a time window (MTP) in the signal OS 11 after normalization (by division by OS2), and Plb is given by the 99:th percentile of the signal values within MTP in the signal OS 12 after normalization (by division by OS2). Thus, in this specific example, Pla and Plb are properties obtained by crossnormalization.

[0119] Generally, the first calculation function Fl may be defined as a linear or nonlinear function of the included properties. The non-linear function may be a polynomial function of any order, for example a second-order function. Thus, the calculation function Fl may include "power terms", in which a property is multiplied by itself to a given power, for example by being squared as in Equation 2. Further, the calculation function Fl may include "interaction terms", which comprise a multiplication of two or more different properties. It should be noted that the extracted properties and the functional dependences used in Equations 1 and 2 are merely given as non-limiting examples.

[0120] FIG. 8B is a flowchart of an example method M2 for operating an OPA 20 to determine the concentration of RBCs in a fluid. The method M2 will be described with reference to the OPA 20 in FIGS 1-2. Optional steps are indicated by dashed lines. Steps S20-S24 are performed by the calculation device 22 and presumes that at least steps S10-S11 and S 15 -S 16 of the method Ml have been performed.

[0121] In step S20, the calculation device 22 obtains the output signal from step Sil (FIG. 8 A) and determines, based on this signal, at least one fourth property, [P4], that represents the first scattered light. As shown by step S20a, [P4] may be determined to represent magnitude and / or variability. As shown, the method M2 may include a step S20b of determining a property for [P4] that represents the magnitude of the first scattered light in the inner angular range 301a, and / or a step S20c of determining a property for [P4] that represents the variability of the first scattered light in the inner angular range 301a. Data analysis indicates that the accuracy of step S20 may be improved if [P4] represents at least one of the properties from steps S20b-S20c.

[0122] In step S21, the calculation device 22 obtains the output signal from step S 16 (FIG. 8 A) and determines, based on this signal, at least one fifth property, [P5], that represents the second scattered light. Data analysis indicates that the accuracy of step S21 may be improved if [P5] represents magnitude. Thus, step S21 may include a step S21a of determining the magnitude of the second scattered light. It is currently believed, backed by experimental data, to be beneficial if [P5] includes the magnitude of the second scattered light in the inner angular range. Thus, step S21a may be replaced by a step S 2 lb of determining the magnitude of the second scattered light in the inner angular range.

[0123] In step S23, the calculation device 22 operates a second calculation function on [P4] and [P5] to determine the TC concentration. The second calculation function corresponds to F2 in FIG. 7.

[0124] In step S24, the calculation device 22 determines the RBC concentration based on the TC concentration from step S23 and the WBC concentration from step SI 8, for example by subtracting the WBC concentration from the TC concentration.

[0125] As shown by step S25, the method M2 may involve an evaluation of the RBC concentration from step S24, for example for detection of a risk for complications with the catheter or illness of the patient. The evaluation may comprise comparing the RBC concentration to a threshold value and, if deemed necessary, alert the caretaker. The evaluation may be performed by analogy with step S 19. Alternatively or additionally, step S25 may involve presenting the RBC concentration to the caretaker. Alternatively or additionally, step S25 may involve storing the RBC concentration for the patient.

[0126] Experiments indicate that the accuracy of step S23 may be improved by step S22.

[0127] In step S22, the calculation device 22 obtains the output signal from step S 14 and determines, based on this signal, at least one sixth property, [P6], that represents the first transmitted light of the first light beam. In step S23, the calculation device 22 operates the second calculation function on [P4], [P5] and [P6] to determine the TC concentration. Data analysis indicates an improvement when [P6] includes a property that represents the variability of the first transmitted light. Thus, step S22 may include a step S22a of determining the temporal variability of the first transmitted light.

[0128] FIG. 9B is a block diagram of an example calculation device 22 which is configured to implement embodiments of the method M2 of FIG. 8B. The second parameter calculation unit 25 is configured to receive at least the signals OS 11 and OSH', which represent the scattered light of the first light beam 300 in the inner angular ranges 301a (FIG. 2B), and the scattered light of the second light beam 400 in the inner angular range 401a (FIG. 2C). As indicated by dashed arrows, the unit 25 may optionally be configured to receive the signal OS2, which represents the transmitted light of the first light beam 300 (FIG. 2B). The unit 25 is configured to determine at least one fourth property, [P4], based on OS 11, in accordance with step S20 (FIG. 8B), and at least one fifth property [P5], based on OSH', in accordance with step S21. As indicated by italic characters, the unit 25 may optionally be configured to determine at least one sixth property, [P6], based on OS2, in accordance with step S22. The RBC calculation unit 26 is configured to determine C_TC by operating the second calculation function, F2, on [P4] and [P5], optionally in combination with [P6], in accordance with step S23. The calculation device 22 is configured to output C_TC for further processing.

[0129] Experiments have been performed to determine the second calculation function F2. The same reference fluids were used as in the experiments presented with reference to FIGS 13-14. The method M2 was performed for each of the reference fluids for different combinations of extracted properties, resulting in predicted values of TC concentration. Some example results are presented in FIGS 15-16.

[0130] FIG. 15 is a graph of predicted TC concentration versus known TC concentration for a linear dependence on three specific properties. The solid line represents a one-to- one relation between predicted and known TC concentrations. The correlation shows that the TC concentration in a fluid can be determined by use of the method M2. In the example of FIG. 15, the predicted TC concentration is determined based on [P4], [P5] and [P6], given by steps S20, S21 and S22 in FIG. 8B. In the illustrated example, [P4] consists of a property P4a representing the magnitude of first scattered light in the inner angular range 301a, [P5] consists of a property P5a representing the magnitude of second scattered light in the inner angular range 401a, and [P6] consists of a property P6a representing the variability of first transmitted light. A linear correlation was used, represented as: y' - a' + b' ■ P4a + c' ■ P5a + d' ■ P6a (3) with y' being the predicted TC concentration, and a'-d' being weights or coefficients given by the linear regression. The second calculation function F2, used by step S23, may be given by Equation (3).

[0131] In the specific example of FIG. 15, P4a is given by the 5 :th percentile of the signal values within a time window (MTP = 90 seconds) in the signal OS 11 (FIG. 2B), P5a is given by the 25:th percentile of the signal values within the MTP in the signal OSH' (FIG. 2C), and P6a is given by the IQR of the signal values within the MTP in the signal OS2 (FIG. 2B) after normalization (by division with a baseline signal). The correlation of data in FIG. 15 yields R = 0.959, MAE = 595.1 and RMSE = 873.3.

[0132] FIG. 16 is a graph of predicted TC concentration versus known TC concentration when including a second order dependence on specific properties. FIG. 16 is based on the same properties P4a, P5a and P6a as FIG. 15. A linear correlation was used, represented as: with y' being the predicted TC concentration, and a'-j' being weights or coefficients given by the linear regression. The second calculation function F2, used by step S23, may be given by Equation (4). If one or more of the coefficients a'-j' is small, it may be set to zero (0).

[0133] The correlation of data in FIG. 16 yields R = 0.990, MAE = 440.9 and RMSE = 253.4.

[0134] It is to be understood that the second calculation function F2 may be defined as a linear or non-linear function of the included properties. The non-linear function may be a polynomial function of any order, for example a second-order function. Thus, the calculation function F2 may include "power terms", in which a property multiplied by itself to a given power, for example by being squared as in Equation 2. Further, the calculation function F2 may include "interaction terms", which comprise a multiplication of two or more different properties. It should be noted that the extracted properties and the functional dependences used in Equations 3 and 4 are merely given as non-limiting examples.

[0135] To give a further non-limiting example, a result comparable to the one in FIG. 15 may be obtained by using only [P4] and [P5] to determine the predicted TC concentration. For example, this may be achieved when [P4] comprises the above- mentioned properties P4a and P5a, and a property P4b that represents the variability of first scattered light in the inner angular range. In a specific example, P4a is given by the median of the signal values within the MTP in the signal OS 11 (FIG. 2B), P4b is given by the 1 :st percentile of the signal values within the MTP in the signal OS 11 after normalization (by subtraction of a baseline signal), and P5a is given by the 25:th percentile of the signal values within the MTP in the signal OSH' (PIG. 2C).

[0136] The available data indicates that calculation functions that have a second order dependence on extracted properties may give a higher accuracy of the predicted concentration compared to calculation functions that have a linear dependence on extracted properties. In the results presented herein, this is particularly noticeable for the second calculation function.

[0137] It is conceivable to use calculation functions that have a third or higher order dependence on extracted properties, or another type of non-linear dependence.

[0138] In some embodiments, the respective calculation function Fl, F2 comprises at least one trained machine learning-based model, which is operated on the included properties or a subset thereof.

[0139] It is also conceivable to employ the concept of stacking models, also known as stacked generalization, to determine the WBC concentration, the TC concentration or the RBC concentration. In stacking models, multiple models are combined to improve accuracy. Thus, one or both of the calculation functions Fl, F2 may comprise stacked models. In one example, the stacked models are configured to determine the same parameter, whereupon the outcomes of the models are combined, typically by another model, to produce a final determination. In another example, a first model is operated on a first plurality of properties, and if the outcome of the first model fulfils a predefined condition, a second model is operated on a second plurality of properties. The first and second pluralities of properties may or may not overlap. The models used as stacking models may be of any type, for example a predefined mathematical function or a machine learning-based model.

[0140] Experiments also indicate that the calculation function Fl may include a single property Pl a representing first scattered light if it also includes a property P3a presenting transmitted light of the first light beam (cf. 300' in FIGS 2A-2B). Such a calculation function Fl is denoted "single scatter function" herein. The property Pla may or may not be obtained by cross-normalization (by division with a baseline signal in OS2). The single scatter function Fl has been found to at least be useful in detecting high WBC concentrations, for example above 500-3000 cells / pL.

[0141] FIG. 18 is a graph of predicted WBC concentration versus known WBC concentration for a single scatter function Fl with a second order dependence on two properties Pla, P3a given by steps S 12 and S 14 in FIG. 8A. Here, the single scatter function Fl is given by y — (a + b ■ Pla + c ■ P3a)z, where Pla represents the variability of first scattered light in the outer angular range 301b, and P3a represents the magnitude of first transmitted light in the angular range 300 (FIG. 2B). In the specific example of FIG. 17, Pla is given by MAD squared (MAD2) for the signal values within a time window (MTP) in the signal OS 12 after normalization (by division with the baseline signal in OS2), and P3a is given by the 50:th percentile of the signal values within MTP in the signal OS2 after normalization (by division with a baseline signal). Thus, in this specific example, Pla is a property obtained by cross-normalization. The single scatter function Fl is not limited to this specific example but generally includes a single property Pla representing first scattered light and at least one property [P3] representing first transmitted light.

[0142] FIG. 19A is a flow chart of an example method M3 that utilizes a single scatter function. The method M3 is performed by the calculation device 22. Steps S10-S11 are the same as in the method Ml of FIG. 8A and results in a signal OS 1, OS 11 or OS 12 that represents first scattered light. Step S 12' corresponds to step S12 in FIG. 8A but only results in a single property that represents the first scattered light. The single property may be a magnitude or a temporal variability and may be determined for either the inner angular range or the outer angular range, if the detection system 31 is configured to define such a confined angular range. Steps S13-S14 are the same as in the method Ml of FIG. 8 A and result in at least one property of first transmitted light, for example its magnitude. Step S 18' corresponds to step S18 in FIG. 8A but operates a single scatter function Fl on the single property from step S 12' and the at least one property from step S14. Step S 18' results in a WBC concentration value. Like in FIG. 8A, the method M3 may also include a step S19 of evaluating the WBC concentration value for detection of peritonitis, presenting the WBC concentration value, or storing the WBC concentration value.

[0143] In some embodiments, the single scatter function Fl may be used as a first model in a stacked model as described hereinabove. As noted, the single scatter function Fl may be particularly useful to detect high WBC concentrations in the effluent. If the first model detects a WBC concentration below a predefined concentration value, for example in the range of 500-3000 cells / pL, a second model may be operated to determine the WBC concentration. Typically, the second model is a further calculation function Fl that yields a higher accuracy for WBC concentrations below the predefined concentration value. The second model may also have an increased robustness to presence of interfering particles (non-WBCs) in the effluent, such as RBCs, fibrine, fat particles, triglycerides, etc. For example, the second model may operate on two or more properties of the first scattered light, such as any of the calculation functions Fl that have been described above with reference to FIGS 13, 14 or 17.

[0144] FIG. 19B is a flow chart of an example method M4 that implements a stacked model comprising a first model and a second model. The first model is the single scatter function of the method M3. The method M4 is performed by the calculation device 22. In the illustrated example, the method M4 starts by performing the method M3, using the single scatter function, to determine a first WBC value. In step S40, the first WBC value is evaluated for certainty (or equivalently, accuracy). If the certainty is deemed insufficient, the method M4 proceeds to operate a second model, which is configured to determine a second WBC value with higher certainty than the first model. In the illustrated embodiment, the second model is the calculation function of the method Ml ("further first calculation function"). If the certainty is deemed sufficient in step S40, step S19 is performed on the first WBC value. In some embodiments, step S40 involves evaluating the first WBC value in relation to a predefined threshold. In some embodiments, the certainty is deemed insufficient when the first WBC value is below the predefined threshold. The second model need not be implemented by the method Ml but may be implemented by any method that provides a second WBC value with higher certainty than the first model.

[0145] The method M4 in FIG. 19B is provided as a basic example of a stacked model. It is conceivable that the stacked model includes further models that are operated in sequence and / or in parallel with the first and second models. Such further models may include calculation functions operating on any combination of properties described herein, as well as other properties or parameters not derived from the output signals of the optical measurement device 31 (FIG. IB).

[0146] Reverting to FIG. 19A, it is realized that the method M3 is converted into an embodiment of the method Ml if steps S 12' and S 18' are replaced by steps S12 and S18. Thus, generally, the present disclosure comprises a method that estimates the WBC concentration as a function of at least one property [Pl] of first scattered light and at least one property [P3] of first transmitted light. This method may be used on its own to determine the WBC concentration or be included as part of a stacked model.

[0147] While the subject of the present disclosure has been described in connection with what is presently considered to be the most practical embodiments, it is to be understood that the subject of the present disclosure is not to be limited to the disclosed embodiments, but on the contrary, is intended to cover various modifications and equivalent arrangements included within the scope of the appended claims.

[0148] For example, the techniques described in the foregoing are not limited to APD but are equally applicable to other types of PD therapy such as CAPD. The techniques are not limited to PD effluent but are equally applicable to other medical fluids that may contain WBCs and RBCs.

[0149] Further, while operations are depicted in the drawings in a particular order, this should not be understood as requiring that such operations be performed in the particular order shown or in sequential order, or that all illustrated operations be performed, to achieve desirable results.

[0150] In the following, clauses are recited to summarize some aspects and embodiments as disclosed in the foregoing. Cl. An optical detection apparatus for determining a concentration of white blood cells in a fluid, said apparatus comprising: a light emitting arrangement (30), which is configured to generate a first light beam with a wavelength within a range of 350-575 nm and arranged to direct the first light beam (300) along a first main direction (310) into a first region (320) in the fluid; a light detection arrangement (31), which is arranged to detect first scattered light originating from the first light beam (300) in the first region (320) and provide a first signal (OS1; OS 11, OS 12) that represents the first scattered light as a function of time; and a calculation arrangement (22), which is configured to: determine, based on the first signal (OS1; OS 11, OS 12), a plurality of first properties ([Pl]) representing the first scattered light; and operate a first calculation function (Fl) on the plurality of first properties ([Pl]) to estimate the concentration of white blood cells in the fluid.

[0151] C2. The apparatus of Cl, wherein the calculation arrangement (22) is configured to estimate the concentration of white blood cells in the presence of red blood cells within the fluid.

[0152] C3. The apparatus of Cl or C2, wherein the light detection arrangement (31) is arranged to detect the first scattered light at an angle to the first main direction (310) of the first light beam (300) through the first region (320), said angle being in a range extending from about 6° to about 35°.

[0153] C4. The apparatus of any preceding clause, wherein the calculation arrangement (22) is configured to determine a respective first property to represent an ensemble of signal values within a respective time window of the first signal (OS1; OS 11, OS 12) or a signal derived from the first signal.

[0154] C5. The apparatus of any preceding clause, wherein the plurality of first properties ([Pl]) comprises a magnitude of the first scattered light.

[0155] C6. The apparatus of any preceding clause, wherein the plurality of first properties ([Pl]) comprises a temporal variability of the first scattered light.

[0156] C7. The apparatus of any preceding clause, wherein the light detection arrangement (31) is configured to detect the first scattered light in a first angular range (301a) and a second angular range (301b), wherein the first angular range (301a) differs from the second angular detection range (301b), and wherein the plurality of first properties ([Pl]) comprise a first magnitude of the first scattered light within the first angular range (301a), and a second magnitude of the first scattered light within the second angular range (301b).

[0157] C8. The apparatus of C7, wherein the first angular range (301a) is nonoverlapping with the second angular range (301b). C9. The apparatus of C7 or C8, wherein the first angular range (301a) is closer to the first main direction (310) than the second angular range (301b).

[0158] CIO. The apparatus of any one of C7-C9, wherein the first angular range (301a) is located within 7°-16° to the first main direction (310), and the second angular range (301b) is located within 16°-35° to the first main direction (310).

[0159] Cll. The apparatus of C6 in combination with any one of C7-C10, wherein calculation arrangement (22) is configured to determine the temporal variability of the first scattered light for the second angular range (301b).

[0160] C12. The apparatus of any preceding clause, wherein the light emitting arrangement (30) is configured to generate a second light beam (400) with a wavelength within a range of 600-1000 nm and arranged to direct the second light beam (400) along a second main direction (410) into a second region (420) in the fluid, wherein the light detection arrangement (31) is arranged to detect second scattered light originating from the second light beam (400) in the second region (420) and provide a second signal (OST) that represents the second scattered light as a function of time, and wherein the calculation arrangement (22) is configured to determine, based on the second signal (OST), at least one second property ([P2]) representing the second scattered light, and operate the first calculation function (Fl) on the plurality of first properties ([Pl]) and the at least one second property ([P2]) to estimate the concentration of white blood cells.

[0161] C13. The apparatus of Cl 2, wherein the calculation arrangement (22) is configured to determine the at least one second property ([P2]) to represent an ensemble of signal values within a second time window of the second signal (OS1') or a signal derived from the second signal.

[0162] C14. The apparatus of C12 or C13, wherein the at least one second property ([P2]) comprises a temporal variability of the second scattered light.

[0163] C15. The apparatus of C14, wherein the light detection arrangement (31) is configured to detect the second scattered light in a third angular range (401a), which is located within 7°-16° to the second main direction (410), and wherein the calculation arrangement (22) is configured to determine the temporal variability of the second scattered light for the third angular range (401a).

[0164] C16. The apparatus of any preceding clause, wherein the light detection arrangement (31) is further configured to detect transmitted light of the first light beam (300) by the first region (320) and provide a third signal (OS2) representing the transmitted light, and wherein the calculation arrangement (22) is configured to normalize the first signal (OS 1) by dividing the first signal (OS1) with a baseline signal given by the third signal (OS2), and determine at least one of the plurality of first properties ([Pl]) based on the thus-normalized first signal (OS1).

[0165] C17. The apparatus of any one of Cl -C 15, wherein the light detection arrangement (31) is further configured to detect transmitted light of the first light beam (300) by the first region (320) and provide a third signal (OS2) representing the transmitted light, and wherein the calculation arrangement (22) is configured to determine, based on the third signal (OS2), at least one third property ([P3]) representing the transmitted light, and operate the first calculation function (Fl) on the plurality of first properties ([Pl]) and the at least one third property ([P3]) to estimate the concentration of white blood cells in the fluid.

[0166] C18. The apparatus of C17, wherein the at least one third property ([P3]) comprises a magnitude of the transmitted light.

[0167] C19. The apparatus of C17 or C18, wherein the light detection arrangement (31) is arranged to detect the transmitted light within an angular range extending to less than 6° from the first main direction (310).

[0168] C20. The apparatus of any preceding clause, wherein the calculation arrangement (22) is further configured to estimate a total particle concentration in the fluid, and estimate a concentration of red blood cells in the fluid as a function of the total particle concentration and the concentration of white blood cells.

[0169] C21. The apparatus of any one of C12-C15, wherein the calculation arrangement (22), to estimate a total particle concentration in the fluid, is further configured to: determine, based on the first signal (OS1; OS 11, OS 12), at least one fourth property ([P4]) representing the first scattered light; determine, based on the second signal (OST), at least one fifth property ([P5]) representing the second scattered light; and operate a second calculation function (F2) on the at least one fourth property ([P4]) and the at least one fifth property ([P5 ]) to estimate the total particle concentration in the fluid.

[0170] C22. The apparatus of C21, wherein the at least one fourth property ([P4]) comprises a magnitude of the first scattered light.

[0171] C23. The apparatus of C22 in combination with any one of C7-C11, wherein the calculation arrangement (22) is configured to determine the magnitude of the first scattered light for the first angular range (301a).

[0172] C24. The apparatus of any one of C21-C23, wherein the at least one fourth property ([P4]) comprises a temporal variability of the first scattered light.

[0173] C25. The apparatus of C24 in combination with any one of C7-C11, wherein the calculation arrangement (22) is configured to determine the temporal variability of the first scattered light for the first angular range (301a). C26. The apparatus of any one of C21-C25, wherein the at least one fifth property ([P5]) comprises a magnitude of the second scattered light.

[0174] C27. The apparatus of C26 in combination with C15, wherein the calculation arrangement (22) is configured to determine the magnitude of the second scattered light for the third angular range (401a).

[0175] C28. The apparatus of any one of C21-C27 in combination with any one of Cl 7- C19, wherein the calculation arrangement (22) is further configured to determine, based on the third signal (OS2), at least one sixth property ([P6]) representing the transmitted light, and operate the second calculation function (F2) on the at least one fourth property ([P4]), the at least one fifth property ([P5]) and the at least one sixth property ( [P6]) to estimate the total particle concentration in the fluid.

[0176] C29. The apparatus of C28, wherein the at least one sixth property ([P6]) comprises a temporal variability of the transmitted light.

[0177] C30. The apparatus of any preceding clause, wherein the first light beam (300) has a spectral width below 20 nm.

[0178] C31. The apparatus of any preceding clause, wherein the first calculation function (Fl) comprises a weighted combination of the plurality of first properties ([Pl]), as well as any further properties determined by the calculation arrangement (22).

[0179] C32. The apparatus of C31, wherein the first calculation function (Fl) is a linear function.

[0180] C33. The apparatus of C31, wherein the first calculation function (Fl) is a nonlinear function.

[0181] C34. An optical detection apparatus for determining a concentration of white blood cells in a fluid, said apparatus comprising: a light emitting arrangement (30), which is configured to generate a first light beam with a wavelength within a range of 350-575 nm and arranged to direct the first light beam (300) along a first main direction (310) into a first region (320) in the fluid; a light detection arrangement (31), which is arranged to detect first scattered light originating from the first light beam (300) in the first region (320) and provide a first signal (OS1; OS 11, OS 12) that represents the first scattered light as a function of time, and to detect transmitted light of the first light beam (300) by the first region (320) and provide a third signal (OS2) representing the transmitted light; and a calculation arrangement (22), which is configured to: determine, based on the first signal (OS1; OS 11, OS 12), at least one first property ([Pl]) representing the first scattered light; determine, based on the third signal (OS2), at least one at least one third property ([P3]) representing the transmitted light; and operate a first calculation function (Fl) on the at least one first property ([Pl]) and the at least one third property ( [P3]) to estimate the concentration of white blood cells in the fluid. C35. The optical detection apparatus of C34, wherein the at least one first property ([Pl]) comprises a magnitude or a temporal variability of the first scattered light.

[0182] C36. The optical detection apparatus of C34 or C35, wherein the at least one third property ([P3]) comprises a magnitude of the transmitted light.

[0183] C37. The optical detection apparatus of any one of C34-C37, wherein the calculation arrangement (22) is configured to normalize the first signal (OS1) by dividing the first signal (OS1) with a baseline signal given by the third signal (OS2), and determine the at least one first property ([Pl]) based on the thus-normalized first signal.

[0184] C38. The optical detection apparatus of any one of C34-C37, wherein the calculation arrangement (22) is configured to perform an evaluation of a first concentration value provided by the first calculation function (Fl) in relation to a predefined concentration value and selectively, based on the evaluation, operate a further first calculation function (Fl) to determine a second concentration value for the white blood cells in the fluid.

[0185] C39. The optical detection apparatus of C38, wherein the calculation arrangement (22) is configured to operate the further first calculation function (Fl) when the first concentration value is below a predefined concentration value.

[0186] C40. The optical detection apparatus of C38 or C39, wherein the calculation arrangement (22) is configured to determine, based on the first signal (OS1; OS 11, OS 12), a plurality of first properties ([Pl]) representing the first scattered light, and operate the further first calculation function (Fl) on the plurality of first properties ([Pl]) to determine the second concentration value for the white blood cells in the fluid.

[0187] C41. The optical detection apparatus of C40, wherein the calculation arrangement (22) is configured to normalize the first signal (OS1) by dividing the first signal (OS1) with a baseline signal given by the third signal (OS2), and determine one or more of the plurality of first properties ([Pl]) based on the thus-normalized first signal (OS1).

[0188] C42. A calculation arrangement for use in the optical detection apparatus of any one of C1-C41.

[0189] C43. An apparatus for automated peritoneal dialysis comprising the optical detection apparatus of any one of C1-C41.

[0190] C44. A computer-implemented method for determining a concentration of white blood cells in a fluid, said method comprising: obtaining (S 10, S 11) a first signal that represents first scattered light received as a function of time by a light detection arrangement from a first region in the fluid when the first region is illuminated by a first light beam, said first light having a wavelength within a range of 350-575 nm and being directed along a first main direction into the first region; determining (S12), based on the first signal, a plurality of first properties representing the first scattered light; and operating (S 18) a first calculation function on the plurality of first properties to estimate the concentration of white blood cells in the fluid.

[0191] C45. The computer-implemented method of C44, further comprising: obtaining (S15, S16) a second signal that represents second scattered light received as a function of time by the light detection arrangement from a second region in the fluid when the second region is illuminated by a second light beam, said second light having a wavelength within a range of 600-1000 nm and being directed along a second main direction into the second region; and determining (S17), based on the second signal, at least one second property representing the second scattered light, wherein the first calculation function is operated on the plurality of first properties and the at least one second property to estimate the concentration of white blood cells in the fluid.

[0192] C46. The computer-implemented method of C44 or C45, further comprising: obtaining (SI 3) a third signal representing transmitted light of the first light beam by the first region; and determining (S14), based on the third signal, at least one third property representing the transmitted light, wherein the first calculation function is operated on the plurality of first properties and the at least one third property to estimate the concentration of white blood cells in the fluid.

[0193] C47. The computer-implemented method of C45, further comprising estimating a total particle concentration in the fluid, wherein said estimating the total particle concentration comprises: determining (S20), based on the first signal, at least one fourth property representing the first scattered light; determining (S21), based on the second signal, at least one fifth property representing the second scattered light; and operating (S23) a second calculation function on the at least one fourth property and the at least one fifth property to estimate the total particle concentration in the fluid.

[0194] C48. The computer-implemented method of C47 in combination with C46, wherein said estimating the total particle concentration comprises: determining (S22), based on the third signal, at least one sixth property representing the transmitted light, wherein the second calculation function is operated on the at least one fourth property, the at least one fifth property, and the at least one sixth property to estimate the total particle concentration in the fluid.

[0195] C49. A computer-implemented method for determining a concentration of white blood cells in a fluid, said method comprising: obtaining (S 10, S 11) a first signal that represents first scattered light received as a function of time by a light detection arrangement from a first region in the fluid when the first region is illuminated by a first light beam, said first light having a wavelength within a range of 350-575 nm and being directed along a first main direction into the first region; obtaining (S13) a third signal that represents transmitted light of the first light beam by the first region; determining (S 12'), based on the first signal, at least one first property representing the first scattered light; determining (S14), based on the third signal, at least one third property representing the transmitted light; and operating (S 18') a first calculation function on the at least one first property and the at least one third property to estimate the concentration of white blood cells in the fluid.

[0196] C50. The method of C49, wherein said determining (S 12') the at least one first property comprises: normalizing the first signal by dividing the first signal with a baseline signal given by the third signal, and determining the at least one first property based on the thus-normalized first signal.

[0197] C51. The method of C49 or C50, further comprising: evaluating (S40) a first concentration value, which is provided by the first calculation function, in relation to a predefined concentration value, and selectively operating (Ml), based on said evaluating (S40), a further first calculation function to estimate a second concentration value for the white blood cells in the fluid.

[0198] C52. The method of C51, wherein said selectively operating (Ml) the further first calculation function comprises: determining (S12), based on the first signal, a plurality of first properties representing the first scattered light, and operating (S 18) the further first calculation function on the plurality of first properties to estimate the second concentration value for the white blood cells in the fluid.

[0199] C53. A computer-readable medium comprising instructions which when executed by processor circuitry causes the processor circuitry to perform the method of any one of C44-C52.

[0200] C54. An optical detection apparatus for determining a total concentration of white blood cells and red blood cells in a fluid, said apparatus comprising: a light emitting arrangement (30), which is configured to generate a first light beam with a wavelength within a range of 350-575 nm and arranged to direct the first light beam (300) along a first main direction (310) into a first region (320) in the fluid and to generate a second light beam (400) with a wavelength within a range of 600-1000 nm and is arranged to direct the second light beam (400) along a second main direction (410) into a second region (420) in the fluid; a light detection arrangement (31), which is arranged to detect first scattered light originating from the first light beam (300) in the first region (320) and provide a first signal (OS1; OS 11, OS 12) that represents the first scattered light as a function of time and to detect second scattered light originating from the second light beam (400) in the second region (420) and provide a second signal (OST) that represents the second scattered light as a function of time; and a calculation arrangement (22), which is configured to: determine, based on the first signal (OS1; OS 11, OS 12), at least one fourth property ([P4]) representing the first scattered light; determine, based on the second signal (OS2'), at least one fifth property ([P5]) representing the second scattered light; and operate a second calculation function (F2) on the at least one fourth property ([P4]) and the at least one fifth property ([P5]) to estimate the total concentration of white blood cells and red blood cells in the fluid.

[0201] C55. A computer-implemented method for determining a total concentration of white blood cells and red blood cells in a fluid, said method comprising: obtaining (S 10, Si l) a first signal that represents first scattered light received as a function of time by a light detection arrangement from a first region in the fluid when the first region is illuminated by a first light beam, said first light having a wavelength within a range of 350-575 nm and being directed along a first main direction into the first region; obtaining (S15, S16) a second signal that represents second scattered light received as a function of time by the light detection arrangement from a second region in the fluid when the second region is illuminated by a second light beam, said second light having a wavelength within a range of 600-1000 nm and being directed along a second main direction into the second region; determining (S20), based on the first signal, at least one fourth property representing the first scattered light; determining (S21), based on the second signal, at least one fifth property representing the second scattered light; and operating (S23) a second calculation function on the at least one fourth property and the at least one fifth property to estimate the total concentration of white blood cells and red blood cells in the fluid.

Claims

42CLAIMS1. An optical detection apparatus for determining a concentration of white blood cells in a fluid, said apparatus comprising: a light emitting arrangement (30), which is configured to generate a first light beam with a wavelength within a range of 350-575 nm and arranged to direct the first light beam (300) along a first main direction (310) into a first region (320) in the fluid, a light detection arrangement (31), which is arranged to detect first scattered light originating from the first light beam (300) in the first region (320) and provide a first signal (OS1; OS 11, OS 12) that represents the first scattered light as a function of time, and a calculation arrangement (22), which is configured to: determine, based on the first signal (OS1; OS 11, OS 12), a plurality of first properties ([Pl]) representing the first scattered light, and operate a first calculation function (Fl) on the plurality of first properties ([Pl]) to estimate the concentration of white blood cells in the fluid.

2. The optical detection apparatus of claim 1, wherein the calculation arrangement (22) is configured to estimate the concentration of white blood cells in the presence of red blood cells within the fluid.

3. The optical detection apparatus of claim 1 or 2, wherein the light detection arrangement (31) is arranged to detect the first scattered light at an angle to the first main direction (310) of the first light beam (300) through the first region (320), said angle being in a range extending from about 6° to about 35°.

4. The optical detection apparatus of any preceding claim, wherein the calculation arrangement (22) is configured to determine a respective first property to represent an ensemble of signal values within a respective time window of the first signal (OS1;OS 11, OS 12) or a signal derived from the first signal (OS1; OS 11, OS 12).

5. The optical detection apparatus of any preceding claim, wherein the plurality of first properties ([Pl]) comprises a magnitude of the first scattered light.

6. The optical detection apparatus of any preceding claim, wherein the plurality of first properties ([Pl]) comprises a temporal variability of the first scattered light.

437. The optical detection apparatus of any preceding claim, wherein the light detection arrangement (31) is configured to detect the first scattered light in a first angular range (301a) and a second angular range (301b), wherein the first angular range (301a) differs from the second angular detection range (301b), and wherein the plurality of first properties ([Pl]) comprise a first magnitude of the first scattered light within the first angular range (301a), and a second magnitude of the first scattered light within the second angular range (301b).

8. The optical detection apparatus of claim 7, wherein the first angular range (301a) is non-overlapping with the second angular range (301b).

9. The optical detection apparatus of claim 7 or 8, wherein the first angular range (301a) is closer to the first main direction (310) than the second angular range (301b).

10. The optical detection apparatus of any one of claims 7-9, wherein the first angular range (301a) is located within 7°-16° to the first main direction (310), and the second angular range (301b) is located within 16°-35° to the first main direction (310).

11. The optical detection apparatus of claim 6 in combination with any one of claims 7-10, wherein calculation arrangement (22) is configured to determine the temporal variability of the first scattered light for the second angular range (301b).

12. The optical detection apparatus of any preceding claim, wherein the light emitting arrangement (30) is configured to generate a second light beam (400) with a wavelength within a range of 600-1000 nm and arranged to direct the second light beam (400) along a second main direction (410) into a second region (420) in the fluid, wherein the light detection arrangement (31) is arranged to detect second scattered light originating from the second light beam (400) in the second region (420) and provide a second signal (OST) that represents the second scattered light as a function of time, and wherein the calculation arrangement (22) is configured to determine, based on the second signal (OST), at least one second property ([P2]) representing the second scattered light, and operate the first calculation function (Fl) on the plurality of first properties ([Pl]) and the at least one second property ([P2]) to estimate the concentration of white blood cells.4413. The optical detection apparatus of claim 12, wherein the calculation arrangement (22) is configured to determine the at least one second property ([P2]) to represent an ensemble of signal values within a second time window of the second signal (OST) or a signal derived from the second signal (OST).

14. The optical detection apparatus of claim 12 or 13, wherein the at least one second property ([P2]) comprises a temporal variability of the second scattered light.

15. The optical detection apparatus of claim 14, wherein the light detection arrangement (31) is configured to detect the second scattered light in a third angular range (401a), which is located within 7°-16° to the second main direction (410), and wherein the calculation arrangement (22) is configured to determine the temporal variability of the second scattered light for the third angular range (401a).

16. The optical detection apparatus of any preceding claim, wherein the light detection arrangement (31) is further configured to detect transmitted light of the first light beam (300) by the first region (320) and provide a third signal (OS2) representing the transmitted light, and wherein the calculation arrangement (22) is configured to normalize the first signal (OS 1) by dividing the first signal (OS1) with a baseline signal given by the third signal (OS2), and determine at least one of the plurality of first properties ([Pl]) based on the thus-normalized first signal (OS1).

17. The optical detection apparatus of any one or claims 1-15, wherein the light detection arrangement (31) is further configured to detect transmitted light of the first light beam (300) by the first region (320) and provide a third signal (OS2) representing the transmitted light, and wherein the calculation arrangement (22) is configured to determine, based on the third signal (OS2), at least one third property ([P3]) representing the transmitted light, and operate the first calculation function (Fl) on the plurality of first properties ([Pl]) and the at least one third property ([P3]) to estimate the concentration of white blood cells in the fluid.

18. The optical detection apparatus of claim 17, wherein the at least one third property ([P3]) comprises a magnitude of the transmitted light.

19. The optical detection apparatus of claim 17 or 18, wherein the light detection arrangement (31) is arranged to detect the transmitted light within an angular range extending to less than 6° from the first main direction (310).

20. The optical detection apparatus of any preceding claim, wherein the calculation arrangement (22) is further configured to estimate a total particle concentration in the fluid, and estimate a concentration of red blood cells in the fluid as a function of the total particle concentration and the concentration of white blood cells.

21. The optical detection apparatus of any one of claims 12-15, wherein the calculation arrangement (22), to estimate a total particle concentration in the fluid, is further configured to: determine, based on the first signal (OS1; OS 11, OS 12), at least one fourth property ([P4]) representing the first scattered light; determine, based on the second signal (OST), at least one fifth property ([P5]) representing the second scattered light; and operate a second calculation function (F2) on the at least one fourth property ([P4]) and the at least one fifth property ([P5 ]) to estimate the total particle concentration in the fluid.

22. The optical detection apparatus of claim 21, wherein the at least one fourth property ([P4]) comprises a magnitude of the first scattered light.

23. The optical detection apparatus of claim 22 in combination with any one of claims 7-11, wherein the calculation arrangement (22) is configured to determine the magnitude of the first scattered light for the first angular range (301a).

24. The optical detection apparatus of any one of claims 21-23, wherein the at least one fourth property ([P4]) comprises a temporal variability of the first scattered light.

25. The optical detection apparatus of claim 24 in combination with any one of claims 7-11, wherein the calculation arrangement (22) is configured to determine the temporal variability of the first scattered light for the first angular range (301a).

26. The optical detection apparatus of any one of claims 21-25, wherein the at least one fifth property ([P5]) comprises a magnitude of the second scattered light.

27. The optical detection apparatus of claim 26 in combination with claim 15, wherein the calculation arrangement (22) is configured to determine the magnitude of the second scattered light for the third angular range (401a).

28. The optical detection apparatus of any one of claims 21-27 in combination with any one of claims 17-19, wherein the calculation arrangement (22) is further configured to determine, based on the third signal (OS2), at least one sixth property ([P6]) representing the transmitted light, and operate the second calculation function (F2) on the at least one fourth property ([P4]), the at least one fifth property ([P5]) and the at least one sixth property ([P6]) to estimate the total particle concentration in the fluid.

29. The optical detection apparatus of claim 28, wherein the at least one sixth property ([P6]) comprises a temporal variability of the transmitted light.

30. The optical detection apparatus of any preceding claim, wherein the first light beam (300) has a spectral width below 20 nm.

31. The optical detection apparatus of any preceding claim, wherein the first calculation function (Fl) comprises a weighted combination of the plurality of first properties ([Pl]), as well as any further properties determined by the calculation arrangement (22).

32. The optical detection apparatus of claim 31, wherein the first calculation function (Fl) is a linear function.

33. The optical detection apparatus of claim 31, wherein the first calculation function (Fl) is a non-linear function.

34. An optical detection apparatus for determining a concentration of white blood cells in a fluid, said apparatus comprising: a light emitting arrangement (30), which is configured to generate a first light beam with a wavelength within a range of 350-575 nm and arranged to direct the first light beam (300) along a first main direction (310) into a first region (320) in the fluid; a light detection arrangement (31), which is arranged to detect first scattered light originating from the first light beam (300) in the first region (320) and provide a first signal (OS1; OS 11, OS 12) that represents the first scattered light as a function of time, and to detect transmitted light of the first light beam (300) by the first region (320) and provide a third signal (OS2) representing the transmitted light, and a calculation arrangement (22), which is configured to:47 determine, based on the first signal (OS1; OS 11, OS 12), at least one first property ([Pl]) representing the first scattered light, determine, based on the third signal (OS2), at least one at least one third property ([P3]) representing the transmitted light, and operate a first calculation function (Fl) on the at least one first property ([Pl]) and the at least one third property ([P3]) to estimate the concentration of white blood cells in the fluid.

35. The optical detection apparatus of claim 34, wherein the at least one first property ([Pl]) comprises a magnitude or a temporal variability of the first scattered light.

36. The optical detection apparatus of claim 34 or 35, wherein the at least one third property ( [P3]) comprises a magnitude of the transmitted light.

37. The optical detection apparatus of any one of claims 34-37, wherein the calculation arrangement (22) is configured to normalize the first signal (OS1) by dividing the first signal (OS1) with a baseline signal given by the third signal (OS2), and determine the at least one first property ([Pl]) based on the thus-normalized first signal.

38. The optical detection apparatus of any one of claims 34-37, wherein the calculation arrangement (22) is configured to perform an evaluation of a first concentration value provided by the first calculation function (Fl) in relation to a predefined concentration value and selectively, based on the evaluation, operate a further first calculation function (Fl) to determine a second concentration value for the white blood cells in the fluid.

39. The optical detection apparatus of claim 38, wherein the calculation arrangement (22) is configured to operate the further first calculation function (Fl) when the first concentration value is below a predefined concentration value.

40. The optical detection apparatus of claim 38 or 39, wherein the calculation arrangement (22) is configured to determine, based on the first signal (OS1; OS 11, OS 12), a plurality of first properties ([Pl]) representing the first scattered light, and operate the further first calculation function (Fl) on the plurality of first properties ([Pl]) to determine the second concentration value for the white blood cells in the fluid.4841. The optical detection apparatus of claim 40, wherein the calculation arrangement (22) is configured to normalize the first signal (OS 1) by dividing the first signal (OS1) with a baseline signal given by the third signal (OS2), and determine one or more of the plurality of first properties ([Pl]) based on the thus-normalized first signal (OS 1).

42. A calculation arrangement for use in the optical detection apparatus of any one of claims 1-41.

43. An apparatus for automated peritoneal dialysis comprising the optical detection apparatus of any one of claims 1-41.

44. A computer-implemented method for determining a concentration of white blood cells in a fluid, said method comprising: obtaining (S10, Sil) a first signal that represents first scattered light received as a function of time by a light detection arrangement from a first region in the fluid when the first region is illuminated by a first light beam, said first light having a wavelength within a range of 350-575 nm and being directed along a first main direction into the first region; determining (S12), based on the first signal, a plurality of first properties representing the first scattered light; and operating (S 18) a first calculation function on the plurality of first properties to estimate the concentration of white blood cells in the fluid.

45. The method of claim 44, further comprising: obtaining (S15, S16) a second signal that represents second scattered light received as a function of time by the light detection arrangement from a second region in the fluid when the second region is illuminated by a second light beam, said second light having a wavelength within a range of 600-1000 nm and being directed along a second main direction into the second region; and determining (S 17), based on the second signal, at least one second property representing the second scattered light, wherein the first calculation function is operated on the plurality of first properties and the at least one second property to estimate the concentration of white blood cells in the fluid.

46. The method of claim 44 or 45, further comprising: obtaining (S13) a third signal representing transmitted light of the first light beam by the first region; and49 determining (S14), based on the third signal, at least one third property representing the transmitted light, wherein the first calculation function is operated on the plurality of first properties and the at least one third property to estimate the concentration of white blood cells in the fluid.

47. The method of claim 45, further comprising estimating a total particle concentration in the fluid, wherein said estimating the total particle concentration comprises: determining (S20), based on the first signal, at least one fourth property representing the first scattered light; determining (S21), based on the second signal, at least one fifth property representing the second scattered light; and operating (S23) a second calculation function on the at least one fourth property and the at least one fifth property to estimate the total particle concentration in the fluid.

48. The method of claim 47 in combination with claim 46, wherein said estimating the total particle concentration comprises: determining (S22), based on the third signal, at least one sixth property representing the transmitted light, wherein the second calculation function is operated on the at least one fourth property, the at least one fifth property, and the at least one sixth property to estimate the total particle concentration in the fluid.

49. A computer-implemented method for determining a concentration of white blood cells in a fluid, said method comprising: obtaining (S10, Sil) a first signal that represents first scattered light received as a function of time by a light detection arrangement from a first region in the fluid when the first region is illuminated by a first light beam, said first light having a wavelength within a range of 350-575 nm and being directed along a first main direction into the first region; obtaining (SI 3) a third signal that represents transmitted light of the first light beam by the first region; determining (S 12'), based on the first signal, at least one first property representing the first scattered light; determining (S14), based on the third signal, at least one third property representing the transmitted light; and operating (S 18') a first calculation function on the at least one first property and the at least one third property to estimate the concentration of white blood cells in the fluid.5050. The method of claim 49, wherein said determining (S 12') the at least one first property comprises: normalizing the first signal by dividing the first signal with a baseline signal given by the third signal, and determining the at least one first property based on the thus-normalized first signal.

51. The method of claim 49 or 50, further comprising: evaluating (S40) a first concentration value, which is provided by the first calculation function, in relation to a predefined concentration value, and selectively operating (Ml), based on said evaluating (S40), a further first calculation function to estimate a second concentration value for the white blood cells in the fluid.

52. The method of claim 51, wherein said selectively operating (Ml) the further first calculation function comprises: determining (S12), based on the first signal, a plurality of first properties representing the first scattered light, and operating (SI 8) the further first calculation function on the plurality of first properties to estimate the second concentration value for the white blood cells in the fluid.

53. A computer-readable medium comprising instructions which when executed by processor circuitry (22a) causes the processor circuitry (22a) to perform the method of any one of claims 44-52.

54. An optical detection apparatus for determining a total concentration of white blood cells and red blood cells in a fluid, said apparatus comprising: a light emitting arrangement (30), which is configured to generate a first light beam with a wavelength within a range of 350-575 nm and arranged to direct the first light beam (300) along a first main direction (310) into a first region (320) in the fluid and to generate a second light beam (400) with a wavelength within a range of 600- 1000 nm and is arranged to direct the second light beam (400) along a second main direction (410) into a second region (420) in the fluid; a light detection arrangement (31), which is arranged to detect first scattered light originating from the first light beam (300) in the first region (320) and provide a first signal (OS1; OS 11, OS 12) that represents the first scattered light as a function of time and to detect second scattered light originating from the second light beam (400) in the second region (420) and provide a second signal (OST) that represents the second scattered light as a function of time; and a calculation arrangement (22), which is configured to:51 determine, based on the first signal (OS1; OS 11, OS 12), at least one fourth property ([P4]) representing the first scattered light; determine, based on the second signal (OS2'), at least one fifth property ([P5]) representing the second scattered light; and operate a second calculation function (F2) on the at least one fourth property ([P4]) and the at least one fifth property ([P5 ]) to estimate the total concentration of white blood cells and red blood cells in the fluid.

55. A computer-implemented method for determining a total concentration of white blood cells and red blood cells in a fluid, said method comprising: obtaining (S10, Si l) a first signal that represents first scattered light received as a function of time by a light detection arrangement from a first region in the fluid when the first region is illuminated by a first light beam, said first light having a wavelength within a range of 350-575 nm and being directed along a first main direction into the first region; obtaining (S15, S16) a second signal that represents second scattered light received as a function of time by the light detection arrangement from a second region in the fluid when the second region is illuminated by a second light beam, said second light having a wavelength within a range of 600-1000 nm and being directed along a second main direction into the second region; determining (S20), based on the first signal, at least one fourth property representing the first scattered light; determining (S21), based on the second signal, at least one fifth property representing the second scattered light; and operating (S23) a second calculation function on the at least one fourth property and the at least one fifth property to estimate the total concentration of white blood cells and red blood cells in the fluid.

Images