NL7905282A

NL7905282A - METHOD AND APPARATUS FOR DETERMINING A RADIATION ABSORPTION DISTRIBUTION IN A SURFACE OF A BODY.

Info

Publication number: NL7905282A
Application number: NL7905282A
Authority: NL
Original assignee: Philips Nv
Priority date: 1979-07-06
Filing date: 1979-07-06
Publication date: 1981-01-08
Also published as: CA1150858A; DE3070773D1; EP0022299B1; US4374419A; JPS5611041A; JPS6411298B2; EP0022299A1; IL60520A; IL60520A0

Description

a Λ « ΡΗΝ 9535 ί N.V. Philips' Gloeilampenfabrieken te Eindhoven.a ΡΗΝ «ΡΗΝ 9535 ί N.V. Philips' Incandescent light factories in Eindhoven.

Werkwijze en inrichting voor het bepalen van een stralings-absorptieverdeling in een vlak van een lichaam.Method and device for determining a radiation absorption distribution in a plane of a body.

De uitvinding heeft betrekking op een werkwijze voor het bepalen van een stralingsabsorptieverdeling in een vlak van een lichaam, waarbij het vlak met een in het vlak liggende platte, waaiervormige, het lichaam doordrin-5 gende en omvattende bundel straling in een veelvoud van equiangulaire richtingen wordt doorstraald en het lichaam gepasseerde straling wordt gedetecteerd voor het genereren van groepen van meetsignalen, die de verzwakking van straling in het lichaam weergeven langs groepen van meetwegen, 10 waarbij elke groep van meetwegen vanuit een knooppunt waaiervormig uitloopt en één waaiervormige bundel straling opdeelt, waarna voor het bepalen van de stralingsabsorptieverdeling elke groep van meetsignalen na elkaar met een reeks van getallen wordt geconvolueerd voor het verkrijgen 15 van een groep convolutiegetallen en de groep convolutiege-tallen op een matrix van elementen worden teruggeprojecteerd langs wegen, die door een bijbehorende groep van meetwegen zijn bepaald, waarbij de matrix van elementen een vaste relatie heeft tot het vlak van het lichaam en een 20 bijdrage in een absorptiewaarde in een element wordt bepaald door het produkt van een convolutiegetal behorende bij de door het element gaande meetweg en van een weegge-tal, dat omgekeerd evenredig is met het kwadraat van de afstand tussen het element en het knooppunt van meetwegen, 7 α n n 9 Λ4) i 'i PHN 9535 . .................................The invention relates to a method for determining a radiation absorption distribution in a plane of a body, wherein the plane with an in-plane flat, fan-shaped, penetrating and enveloping beam of radiation in a plurality of equiangular directions irradiated and radiation passed through the body is detected to generate groups of measurement signals representing the attenuation of radiation in the body along groups of measurement paths, each group of measurement paths extending from a node fan-shaped and dividing one fan-shaped beam of radiation, after which determining the radiation absorption distribution, each group of measurement signals is consecutively sequenced with a series of numbers to obtain a group of convolution numbers and the group of convolution numbers on a matrix of elements are projected back along roads, which are by an associated group of measurement paths determined, incl The matrix of elements has a fixed relation to the plane of the body and a contribution to an absorption value in an element is determined by the product of a convolution number corresponding to the measuring path passing through the element and of a weight that is reversed is proportional to the square of the distance between the element and the junction of measuring paths, 7 α nn 9 Λ4) i 'i PHN 9535. .................................

2 waarbij een absorptieverdeling wordt gevormd door de in de elementen op te bouwen absorptiewaarden door sommatie van de bijdragen.2 wherein an absorption distribution is formed by the absorption values to be built up in the elements by summing the contributions.

De uitvinding heeft verder betrekking op een in-5 richting voor het bepalen van een stralingsabsorptieverde- t •r ling in een vlak van een lichaam, welke inrichting bevat: een strelingsbron voor het opwekken van een platte, in het vlak liggende waaiervormige bundel het lichaam doordringende straling, een detectorinrichting voor het detecteren van de 10 het lichaam gepasseerde straling en het leveren van groepen van meetsignalen, een draaggestel voor de stralingsbron en de detectorinrichting, aandrijfmiddelen voor het bewegen van ten minste de stralingsbron in het draaggestel voor het uit een veelvoud van richtingen doorstralen van het 15 lichaam, waarbij de detectorinrichting voor elke richting een groep van meetsignalen levert, die een stralingsver-zwakking weergeven langs een groep van, van de stralingsbron waaiervormig uitgaande meetwegen, een bewerkings-inrichting voor het convolueren van een groep van meet-20 signalen met een reeks getallen en voor het terugprojecteren van de aldus geconvolueerde meetsignalen voor verkrijgen van absorptiewaarden, een geheugeninrichting voor het opslaan van de absorptiewaarden en een weergeefinrichting voor het weergeven van de absorptiewaarden in elemen-25 ten van een matrix, waarin de stralingsabsorptieverdeling wordt afgebeeld.The invention further relates to a device for determining a radiation absorption distribution in a plane of a body, comprising apparatus: a source of radiation for generating a flat, in-plane fan-shaped beam penetrating radiation, a detector device for detecting the radiation passed through the body and supplying groups of measuring signals, a support frame for the radiation source and the detector device, driving means for moving at least the radiation source in the support frame for a plurality of radiating directions from the body, the detector device providing for each direction a group of measuring signals representing a radiation attenuation along a group of measuring paths radiating fan-shaped from the radiation source, a processing device for convolving a group of measuring 20 signals with a series of numbers and for backprojecting the aldu s convolved measurement signals for obtaining absorption values, a memory device for storing the absorption values and a display device for displaying the absorption values in elements of a matrix, in which the radiation absorption distribution is displayed.

Een dergelijke werkwijze en inrichting is beschreven in het ter inzage gelegde Nederlandse octrooi-schrift 7614230* In de genoemde aanvrage is duidelijk aan-30 gegeven hoe de meetsignalen nodig voor het reconstrueren van een absorptieverdeling worden gemeten met behulp van een waaiervormige bundel röntgenstraling en een rij van detectoren, die voor elke positie van de röntgenbron een groep van meetsignalen levert. De reeks van getallen waar-35 790 5 2 82 * * 3 ΡΗΝ 9535 mee elk van deze groepen van meetsignalen geconvolueerd moet worden voor het berekenen van de groepen van convolu-tiegetallen is eveneens in de genoemde aanvrage afgeleid.Such a method and device is described in Netherlands Patent Specification 7614230 laid open to public inspection. In the said application, it is clearly indicated how the measuring signals required for reconstructing an absorption distribution are measured using a fan-shaped beam of X-rays and a row of detectors, which provide a group of measurement signals for each position of the X-ray source. The series of numbers with which to convolve each of these groups of measurement signals to calculate the groups of conversion numbers has also been derived in said application.

Gok is aangeduid dat een absorptiewaarde in een bepaald 5 punt wordt bepaald door sommatie van de convolutiegetallen 4 behorende bij wegen, die door dat punt gaan, waarbij elk convolutiegetal voor sommatie met een weeggetal woxdt vermenigvuldigd, dat omgekeerd evenredig is met de afstand van het punt tot de bij het convolutiegetal horende positie van de rontgenbron. De beschreven wijze van terugprojectie is een bewerkelijke en tijdrovende zaak hetgeen van nadeel is.It has also been suggested that an absorbance value in a given 5 point is determined by summing the convolution numbers 4 associated with roads passing through that point, where each convolution number is multiplied by a weight number which is inversely proportional to the distance from the point to the position of the X-ray source corresponding to the convolution number. The described manner of back projection is a laborious and time-consuming matter, which is of disadvantage.

De uitvinding beoogt een werkwijze en inrichting aan te geven, waarmee en waarin een terugprojectie op relatief eenvoudige wijze wordt uitgevoerd.The object of the invention is to indicate a method and apparatus with which and in which a back projection is carried out in a relatively simple manner.

15 De werkwijze volgens de uitvinding heeft daartoe tot kenmerk, dat de matrix van elementen een polaire matrix is, waarin een element wordt geadresseerd met een afstands-coördinaat en een hoekcoördinaat, waarbij de weeggetallen in een polaire matrix zijn opgeslagen en met dezelfde af- 20 standscoördinaat en eenzelfde hoekcoördinaat worden geadresseerd, welke hoekcoördinaat na convolutie respectievelijk terugprojectie van een groep meetsignalen respectievelijk van een groep convolutiegetallen met een inkrement wordt vergroot, dat gelijk is aan het hoekinkrement tussen de bij 25 de op elkaar volgend te verwerken groepen horende richtingen.The method according to the invention is therefore characterized in that the matrix of elements is a polar matrix in which an element is addressed with a distance coordinate and an angle coordinate, the weight numbers being stored in a polar matrix and with the same dimensions. The position coordinate and the same angle coordinate are addressed, which angle coordinate is increased after convolution or backprojection of a group of measuring signals or of a group of convolution numbers with an increment equal to the angular increment between the directions corresponding to the groups to be processed.

Een zodanige werkwijze heeft het voordeel, dat bij terugprojectie van een groep van convolutiegetallen steeds weer dezelfde weeggetallen kunnen worden benut. De 30 geometrie van de terugprojectie verandert ten opzichte van de geometrie van de elementen van de reconstructiematrix niet uitgezonderd de hoekverdraaiing, die de stralingsbron ondergaat tussen het meten van twee groepen van meetsignalen c.q. convolutiegetallen. Door daarbij gebruik te maken 35 van een tweede matrix, waarin de weeggetallen zijn opgeslagen die eveneens een adressering heeft in poolcoördinaten, zijn na verdraaiing van de weeggetallen in de tweede matrix ~över een hoek die gelijk is aan de hoek waarover de rönt-~ i PHN 9535 _____________________________________ _____ ________...........Such a method has the advantage that the same weight numbers can be used again and again when projecting a group of convolution numbers. The geometry of the backprojection does not change with respect to the geometry of the elements of the reconstruction matrix, except for the angular rotation, which the radiation source undergoes between the measurement of two groups of measurement signals or convolution numbers. By using a second matrix, in which the weighing numbers are stored, which also has an address in polar coordinates, after rotation of the weighing numbers in the second matrix, an angle is equal to the angle over which the X-ray PHN 9535 _____________________________________ _____ ________...........

4 genbron is verdraaid, de weeggetallen met eenzelfde adressering op te roepen als de mogelijk eerder berekende bijdrage. Na vermenigvuldiging van beide met elkaar wordt de bijdrage bij een op het gegeven adres opgeslagen waarde op-5 geteld. De verdraaiing van de weeggetallen in de tweede matrix is in feite niets anders dan een vast inkrement in de hoekcoördinaat in de adressering van een weeggetal ten opzichte van de hoekcoördinaat in de adressering van de plaats, waaraan de te berekenen bijdrage wordt toegewezen.4 source is rotated, to call up the weighing numbers with the same addressing as the contribution previously calculated. After multiplying the two together, the contribution is added to -5 at a value stored at the given address. The rotation of the weighing numbers in the second matrix is in fact nothing more than a fixed increment in the angular coordinate in the addressing of a weighing number relative to the angular coordinate in the addressing of the place to which the contribution to be calculated is assigned.

^ Een inrichting volgens de uitvinding heeft tot kenmerk, dat de geheugeninrichting bevat een eerste geheugen met polaire adressering, waarin de absorptiewaarden door adressering met een afstandscoördinaat en een hoekco-ordinaat zijn op te slaan en een tweede geheugen met een ^ polaire adressering, waarin weeggetallen voor gebruik bij het terugprojecteren van de convolutiegetallen zijn opgeslagen. Een dergelijke inrichting heeft het voordeel, dat in tegenstelling tot inrichtingen met een geheugen, waarinA device according to the invention is characterized in that the memory device contains a first memory with polar addressing, in which the absorption values can be stored by addressing with a distance coordinate and an angular coordinate and a second memory with a polar addressing, in which weighing numbers stored for use in backprojecting the convolution numbers. Such a device has the advantage that, unlike devices with a memory, in which

de elementen met carthesische coördinaten worden geadres-2Dthe elements with Cartesian coordinates are addressed-2D

seerd, de geheugenruimte voor de weeggetallen relatief klein is en dus goedkoop. Wordt namelijk een geheugen in carthesische coördinaten toegepast, dan moet voor elkebron- positie (andere doorstralingsrichting door het lichaam en in feite door de matrix) een afzonderlijke verzameling van 25 weeggetallen worden toegepast, die wel van te voren berekend kunnen worden hetgeen dan wel een enorme geheugenruimte vraagt en dus duur is.The memory space for the weight numbers is relatively small and therefore cheap. If a memory in Cartesian coordinates is used, then for each source position (different irradiation direction through the body and in fact through the matrix) a separate set of 25 weighing numbers must be applied, which can be calculated in advance, which is an enormous requires memory space and is therefore expensive.

De werkwijze en inrichting volgens de uitvinding zal worden toegelicht aan de hand van een in een tekening 30 weergegeven voorbeeld, in welke tekening fig. 1 een schematisch overzicht geeft van een computertomografie-inrichting volgens de uitvinding, fig. 2 een matrix van elementen, aan de hand waarvan het bepalen van bijdragen van absorptiewaarden vol-35 gens de uitvinding zal worden toegelicht, toont, fig. 3 een blokschema van een uitvoeringsvorm van een bewerkingsinrichting en een geheugeninrichting voor de 790 5 2 82 « a* PHN 9535 5 computertomografie-inrichting ult fig. 1 weergeeft, fig. 4 een blokschema van een voorkeursuitvoeringsvorm van een deel van een bewerkingsinrichting en een geheugeninrichting voor de computertomografie-inricbting 5 uit fig. 1 toont.The method and device according to the invention will be elucidated on the basis of an example shown in a drawing, in which figure 1 shows a schematic overview of a computed tomography device according to the invention, figure 2 shows a matrix of elements. With the aid of which the determination of contributions of absorbance values according to the invention will be explained, Fig. 3 shows a block diagram of an embodiment of a processing device and a memory device for the 790 5 2 82 a PHN 9535 5 computer tomography device. Ult, Fig. 1, Fig. 4 shows a block diagram of a preferred embodiment of a part of a processing device and a memory device for the computed tomography device 5 of Fig. 1.

Een computertomografie-inricbting, zoals in fig.A computed tomography arrangement, as shown in FIG.

1 schematisch is weergegeven . bevat een stralenbron 1 die bij voorkeur een röntgenstralenbron is, die echter ook uit een radioactieve isotoop bijvoorbeeld Am 241 kan bestaan.1 is shown schematically. contains a ray source 1 which is preferably an X-ray source, but which may also consist of a radioactive isotope, for example Am 241.

10 Met behulp van een diafragma 2 wordt de door de stralenbron 1 uitgézonden straling tot een in een vlak liggende divergente stralingsbundel 3 gecollimeerd, waarbij de dikte van de stralingsbundel 3 loodrecht op het vlak bijvoorbeeld tussen 3 en 25 mm ligt en de divergentie ervan in het vlak 15 door de hoek oC ±s bepaald. De stralingsbundel 3 valt op een detectorrij 4 die uit afzonderlijke, de straling metende detectoren 5 bestaat, die stralenbundels 3a- definiëren, waarbij de breedte en de afstand van de afzonderlijke detectoren 5 ten opzichte van elkaar de ruimtelijke nauwkeurig-20 heid bepalen, waarmee een op een objecttafel 6 liggend object 7 wordt afgetast. De detectorrij 4, die symmetrisch ten opzichte van een centrale straal 8 is gepositioneerd, bevat bijvoorbeeld 300 detectoren 5, waarbij de afstand van de middelpunten van twee detectoren 5 ten opzichte van el-25 kaar enkele millimeters bedraagt. Als detector kan ook een lange, met gas gevulde ionisatiekamer worden gebruikt, waarin op een rij geplaatste afzonderlijke gebieden detecterende electroden zijn geplaatst. Het object 7 staat loodrecht op het vlak van de stralingsbundel 3 in de lengte-30 richting van de as 9» die binnen het object 7 ligt en de middelpuntsas van het cirkelvormige draaggestel 10 weergeeft, verschuifbaar, zodat verschillende lagen van het object 7 kunnen worden doorstraald.With the aid of a diaphragm 2, the radiation emitted by the radiation source 1 is collimated into an in-plane divergent radiation beam 3, the thickness of the radiation beam 3 perpendicular to the plane being for example between 3 and 25 mm and its divergence in the plane 15 determined by the angle oC ± s. The radiation beam 3 falls on a detector row 4 consisting of separate radiation-measuring detectors 5, which define beams 3a, the width and the distance of the individual detectors 5 from each other determining the spatial accuracy, with which a object 7 lying on an object table 6 is scanned. The detector row 4, which is positioned symmetrically with respect to a central beam 8, contains, for example, 300 detectors 5, the distance of the centers of two detectors 5 from each other being a few millimeters. As a detector, a long gas-filled ionization chamber can also be used, in which disposed individual areas of detecting electrodes are arranged. The object 7 is perpendicular to the plane of the radiation beam 3 in the longitudinal direction of the axis 9, which lies within the object 7 and represents the center axis of the circular supporting frame 10, so that different layers of the object 7 can be irradiated.

Het systeem stralingsbron 1 -detectorrij 4 is om de as 9 draaibaar geplaatst, zodat een laag van het object 7 in verschillende, in de laag liggende richtingen met behulp van de stralenbundel 3 kan worden doorstraald. De draaiing van het draaggestel 10, dat met behulp van “Tagers ' 7Qn K 9 A9 ΡΗΝ 9535 6 ί -4 11 wordt geleid, geschiedt door aandrijfmiddelen, zoals een tandwiel 12 en een motor 13· De draaiing van het draagge-stel 10 kan zowel continu alsook in stappen geschieden, waarbij in het laatste geval na elke stap het object 7 met 5 de stralenbron 1 in een flits wordt doorstraald.The system radiation source 1 detector row 4 is rotatably arranged about the axis 9, so that a layer of the object 7 can be irradiated in different, layer-lying directions with the aid of the beam 3. The rotation of the support frame 10, which is guided by means of "Tagers" 7Qn K 9 A9 ΡΗΝ 9535 6 ί -4 11, is effected by drive means, such as a gear wheel 12 and a motor 13 · The rotation of the support frame 10 both continuously and in steps, in the latter case after each step the object 7 with 5 the radiation source 1 is irradiated in a flash.

De meetsignalen van de detectoren 5 worden met behulp van de versterker 14 versterkt en aan een signaalom-zetter 15 toegevoerd, waarin de meetsignalen worden gecorrigeerd op bekende wijze voor "offset", gerefereerd aan een 10 referentiewaarde, gedigitaliseerd, gelogaritme erd en gecali-breerd aan de hand van de in de signaalomzetter opgenomen logaritme- en calibratietabellen. Aan de uitgang van de omzetter 15 worden digitale meetwaarden aan de bewerkingsin-richting 16 aangeboden. De in digitale vorm gezette meet-15 waarden worden met behulp van de bewerkingsinrichting 16 in een reconstructiebeeld weergevende absorptiewaarden omgezet, die in de geheugeninrichting 17 worden opgeslagen. De berekende absorptiewaarden kunnen op een weergeefinrichting zoals monitor 18 worden afgebeeld. Een teller 19 telt het aan-20 tal meetsignalen, dat per groep van meetsignalen aan de bewerkingsinrichting 16 wordt toegevoerd. Zodra het aantal meetsignalen overeenkomt met het aantal detectoren 5» wordt een stuurschakeling 20 geactiveerd, die de motor 13 kort aandrijft en daarmee een rotatie van het draaggestel 10 be-25 werkstelligt. Daarna wordt de volgende groep van meetsignalen gemeten enzovoort. Met een optische opnemer* 30 wordt door tellen van de tanden van het tandwiel 12 de hoekver, -draaiing Θ tussen de opeenvolgende meetreeksen bepaald. De door de opnemer 30 opgewekte pulsen worden aan de bewer-3° kingsinrichting 16 toegevoerd, zodat in combinatie met de in de bewerkingsinrichting vastgelegde gegevens over de geometrische opbouw van het draaggestel 10, en de bron 1 met de detectorinrichting k_ de coördinaten van alle meetwegen van elke groep van meetsignalen bepaalbaar zijn, 55 Het is voordelig gebleken, dat de afstand tussen de stralenbron 1 en het object 7 kan worden aangepast aan de diameter van het object 7· Hiertoe is het systeem stra-“---lenbron 1 -detectorrij k op een drager 21 gemonteerd^ cTïe 7905282 PHN 9535 -........The measuring signals of the detectors 5 are amplified with the aid of the amplifier 14 and applied to a signal converter 15, in which the measuring signals are corrected in known manner for "offset", referred to a reference value, digitized, logarithmic and calibrated. using the logarithm and calibration tables included in the signal converter. Digital measured values are presented to the processing device 16 at the output of the converter 15. The measured values converted into digital form are converted by means of the processing device 16 into a reconstruction image displaying reconstruction image, which are stored in the memory device 17. The calculated absorbance values can be displayed on a display device such as monitor 18. A counter 19 counts the number of measuring signals which are applied to the processing device 16 per group of measuring signals. As soon as the number of measuring signals corresponds to the number of detectors 5, a control circuit 20 is activated, which briefly drives the motor 13 and thereby effects rotation of the support frame 10. Then the next group of measurement signals is measured and so on. With an optical sensor * 30, the angular rotation Θ between the successive measuring series is determined by counting the teeth of the gear wheel 12. The pulses generated by the sensor 30 are applied to the processing device 16, so that, in combination with the data on the geometric structure of the support frame 10 and the source 1 with the detector device k_, the coordinates of all measuring paths are recorded in the processing device. of each group of measuring signals are determinable, 55 It has proved advantageous that the distance between the beam source 1 and the object 7 can be adapted to the diameter of the object 7. k mounted on a support 21 ^ cTie 7905282 PHN 9535 -........

7 » * langs de geleiderails 22 op lagers 23 en door middel van een met een motor 24 gekoppelde tandwielaandrijving 25 kan worden verplaatst. Een stuurschakeling 26 is bijvoorbeeld door middel van een handschakelaar 27 bedienbaar, de scha-5 keling 26 kan echter ook automatisch worden bediend. Vóór de aanvang van de meting worden de meetsignalen van twee detec oren 5’ en 5JI via de signaalomzetter 15 aan de stuurschakeling 26 toegevoerd. De drager 21 wordt zodanig verplaatst, dat het meetsignaal van de detector 511 maxi-10 maal wordt, terwijl het meetsignaal van de detector 5' een iets kleinere waarde heeft. In dit geval ontvangt de detector 5 * * straling, die niet door het object 7» maar volledig door de ruimte gaat, die het object 7 omgeeft, terwijl de straling, die de detector 5’ meet, door het object ^ 7 is afgezwakt. Daarna wordt de stuurschakeling 26 vergren deld, om tijdens de opname de afstand tussen de stralenbron 1 en de draaiingsas 9 constant te houden.7 »* can be moved along the guide rails 22 on bearings 23 and by means of a gear drive 25 coupled to a motor 24. A control circuit 26 can be operated, for example, by means of a hand switch 27, but the circuit 26 can also be operated automatically. Before the start of the measurement, the measuring signals of two detectors 5 and 5JI are applied to the control circuit 26 via the signal converter 15. The carrier 21 is displaced such that the measuring signal of the detector 511 becomes maximally 10 times, while the measuring signal of the detector 5 'has a slightly smaller value. In this case, the detector 5 * * receives radiation, which does not pass through the object 7 »but completely through the space surrounding the object 7, while the radiation, which the detector 5 measures, is attenuated by the object 7. The control circuit 26 is then locked in order to keep the distance between the beam source 1 and the axis of rotation 9 constant during the recording.

Aan de hand van fig. 2 zal de werkwijze en de computertomografie-inrichting volgens de uitvinding worden on toegelicht. De bron 1 neemt tijdens de meetprocedure successievelijk de positie θ ^, θ^, θ^, etc. op een cirkelomtrek 32 met middelpunt 9 in, waarbij het lichaam 7 en de een vaste relatie daarmee hebbende matrix met de elementen € stationair blijven. De elementen ζ van de matrix zijn 25 ingesloten door concentrische cirkels 31 en door door het centrum 9 gaande stralen 33» waarvan slechts een gedeelte is weergegeven. De hoek die de stralen 33 paarsgewijs insluiten, is gelijk aan de hoekA© > die twee naburige ver- bindingsli jnen tussen de posities θ , Θ en het centrum 9 30 1 ^ insluiten. Van een groep van meetsignalen, die met de bron in positie langs een groep van meetwegen (slechts een gedeelte is in de figuur weergegeven) zijn bepaald, worden op op zichzelf bekende wijze convolutiegetallen bepaald, die via de meetwegen 37 .» 37 » 37 Λ in de elementen 21™ 1 Π Π4" 1 ) · t £ moeten worden teruggeprojecteerd om bijdragen in de ab-sorptiewaarde in dat element £ te bepalen. Hiertoe wordt het convolutiegetal gewogen door vermenigvuldiging met een weeggetal, die bijvoorbeeld evenredig is met het gemeen? PHN 9535 , _____________________________________ ___________ 8The method and the computed tomography device according to the invention will be elucidated with reference to Fig. 2. During the measuring procedure, the source 1 successively assumes the position θ ^, θ ^, θ ^, etc. on a circumference 32 with center 9, the body 7 and the matrix having a fixed relationship with the elements remaining stationary. The elements ζ of the matrix are enclosed by concentric circles 31 and by rays 33 »passing through the center 9, of which only a part is shown. The angle which the rays 33 enclose in pairs is equal to the angle A ©> which two adjacent connecting lines enclose between positions θ, Θ and the center 9 30 1 ^. From a group of measuring signals, which are determined with the source in position along a group of measuring paths (only a part is shown in the figure), convolution numbers are determined in a manner known per se, which are measured via the measuring paths 37. » 37 »37 Λ in elements 21 ™ 1 Π Π4" 1) · t £ must be projected back to determine contributions in the absorbance value in that element £. To do this, the convolution number is weighted by multiplying it by a weighting number, which is for example proportional in common? PHN 9535, _____________________________________ ___________ 8

4 V4 V.

schappelijk oppervlak van het element £ en de bijbehorende meetweg 37 en omgekeerd evenredig is met de afstand R in het kwadraat tussen de bron 1 en het betreffende element £. De weeggetallen zijn opgeslagen in een matrix dat eenzelf-B de organisatie kent als de matrix van elementen £ voor het weergeven van de stralingsabsorptieverdeling. In elke ge-heugencel van de matrix liggen (ligt) twee (een) weegfactoren) opgeslagen, waarmee een convolutiegetal dat bij een meetweg hoort, die door het element £ loopt, wordt verme-10 nigvuldigd. Verder liggen in de geheugencel bij de weegfactoren) het (de) meetwegnummer(s) (..n-1, n, n+1,n+2,..) opgeslagen, die de meetweg 37n 37n» ^^n+1’ e^c* aan“ geeft, waarbij het convolutiegetal behoort,waarmee het in dezelfde geheugencel opgeslagen weegfactor moet worden ver-10 menigvuldigd. In de meeste geheugencellen van de weegfactor-matrix zijn twee weegfactoren en twee meetwegnummers opgeslagen. Voor sommige elementen zal voor het berekenen van een bijdrage slechts (zoals uit de figuur blijkt) één weeggetal en één meetwegnummer nodig zijn.surface area of the element £ and the associated measuring path 37 and inversely proportional to the distance R squared between the source 1 and the respective element £. The weighing numbers are stored in a matrix which self-B knows the organization as the matrix of elements £ to represent the radiation absorption distribution. Two (one) weighting factors (stored) are stored in each memory cell of the matrix, thereby multiplying a convolution number associated with a measuring path passing through the element £. Furthermore, in the memory cell at the weighting factors) the measuring path number (s) (..n-1, n, n + 1, n + 2, ..) are stored, which measures the measuring path 37n 37n »^^ n + 1 "e ^ c * indicates", which includes the convolution number, by which to multiply the weighting factor stored in the same memory cell. In most memory cells of the weighting factor matrix, two weighting factors and two measuring path numbers are stored. For some elements, only one weighing number and one measuring path number will be required to calculate a contribution (as can be seen from the figure).

20 Nu is eenvoudig in te zien dat, indien alle con- volutiewaarden behorende bij de bronpositie zijn verwerkt,. de convolutiewaarden behorende bij de bronpositie kunnen worden verwerkt, nadat de matrix, waarin de weeggetallen en meetwegnummers zijn opgeslagen ten opzichte van 20 de matrix van elementen over een hoek Δ9 ( = Δψΐ) is verdraaid (met de verdraaiingszin van de bron 1 meel). Nu is de configuratie van beide matrices ten opzichte van de bron 1 niet gewijzigd, uitgezonderd de hoekcoördinaat van de matrix van elementen. De matrix met de weegfactor en de meet-20 wegnummers is echter identiek gebleven, zodat elk weeggetal en bijbehorend meetwegnummer onveranderd kan worden benut. Hetgeen het mogelijk maakt het weeggetal van te voren uit te rekenen en samen met het meetwegnummer op te slaan in een geheugen. Verder is in te zien dat de matrix ten op-35 zichte van de verbindingslijn tussen bron 1 en middelpunt 9 symmetrisch is, waardoor in feite het voldoende is om de weegfactoren en bijbehorende meetwegnummers van een helft -van"de matrix op te slaan.20 Now it is easy to see that, if all conversion values associated with the source position have been processed ,. the convolution values associated with the source position can be processed after the matrix, in which the weighing numbers and measuring path numbers are stored, relative to the matrix of elements is rotated through an angle Δ9 (= Δψΐ) (with the rotation sense of the source 1 flour). Now the configuration of both matrices relative to source 1 has not changed, except for the angular coordinate of the array of elements. However, the matrix with the weighting factor and the measuring 20 road numbers has remained identical, so that every weighing number and associated measuring road number can be used unchanged. Which makes it possible to calculate the weighing number in advance and to store it in a memory together with the measuring path number. It can further be seen that the matrix is symmetrical with respect to the connecting line between source 1 and center 9, so that it is in fact sufficient to store the weighting factors and associated measuring path numbers of half of the matrix.

79Ö5282 * *.79Ö5282 * *.

PHN 9535 9PHN 9535 9

Zoals uit de figuur blijkt zijn niet alle elementen £ even groot. Om een nagenoeg uniforme informatiedicht-heid te krijgen worden de in naburige elementen (bijvoorbeeld tussen twee concentrische cirkels 31 bepaalde absorp-5 tiewaarden opgeteld, waarbij de som van de oppervlakte van de afzonderlijke elementen een benadering van een 'Tstan-daardoppervlak" dient te zijn. Heeft het centraal element 40 een straal r* (element 40 behoeft niet opgedeeld te worden, omdat dit toch op de plaats blijft ongeacht de rotatie van 10 de matrix) dan is het oppervlak lt»r2 . Het oppervlak van een element £ tussen de concentrische cirkels 31 niet een straal n%r en (n-l)*r is gelijk aan 5“ . ^(n.r)2 - (n-1)2 . r2?As can be seen from the figure, not all elements are equal in size. In order to obtain a nearly uniform information density, the absorption values determined in neighboring elements (for example between two concentric circles 31) are added, whereby the sum of the area of the individual elements must be an approximation of a "standard standard area" If the central element 40 has a radius r * (element 40 does not have to be divided, since it remains in place regardless of the rotation of the matrix), the surface area lt »r2. concentric circles 31 not a radius n% r and (nl) * r equals 5 ”. ^ (nr) 2 - (n-1) 2. r2?

15 M15 M

waarbij M het aantal elementen tussen twee cirkels 31 is.where M is the number of elements between two circles 31.

Een aantal malen P van dit element is dan gelijk aan 7¾ .r2. Hieruit volgt dat P = M/2n, waarbij P een geheel getal moet zijn. Is het aantal bronposities O·: 3^0 (een groep meet-signalen per graadrotatie). Dan ontstaan er in de tweede ring na sommatie over 90 elementen 4 absorptiewaarden met een standaardoppervlak. In de derde ring ontstaan 6 absorptiewaarden na sommatie over steeds 60 elementen, in de vierde ring ontstaan 8 waarden door sommatie over 45 elementen 25 enzovoort. De in de geheugeninrichting 17 (zie fig. 1) opgeslagen groepen van meetsignalen worden elk door de bewer-kingsinrichting 16 geconvolueerd met een op zichzelf bekende reeks getallen en daarna weer in het daarvoor bestemd deel 17a (fig. 3) van de geheugeninrichting 17 opgeslagen.A number of times P of this element equals 7¾ .r2. It follows that P = M / 2n, where P must be an integer. Is the number of source positions O: 3 ^ 0 (a group of measurement signals per degree rotation). Then in the second ring, after summation over 90 elements, 4 absorption values with a standard surface are created. In the third ring 6 absorption values arise after summation over 60 elements in each case, in the fourth ring 8 values arise through summation over 45 elements 25 and so on. The groups of measuring signals stored in the memory device 17 (see Fig. 1) are each convolved by the processing device 16 with a series of numbers known per se and then stored again in the designated part 17a (Fig. 3) of the memory device 17 .

30 .30.

Met een klok 41 worden pulsen opgewekt, die via een vierde-ler 43 aan een ψ-teller 45 worden aangeboden. De uitgang van de teller 45 vormt een eerste deel van een adres, waarop weeggetallen en meetwegnummers in een deel 17b van het geheugen 17 zijn opgeslagen. Via een M-deler 47 worden de 35 pulsen verder toegevoerd aan een ringtelier 49, waarvan de uitgang een verder deel van het adres vormt, waarop de weeggetallen en meetwegnummers zijn opgeslagen. Op de adres- 79ÖÏ2I12 * ^ PHN 9535 .. ......... _ ....................Pulses are generated with a clock 41, which are applied to a ψ-counter 45 via a fourth-timer 43. The output of the counter 45 forms a first part of an address on which weighing numbers and measuring path numbers are stored in a part 17b of the memory 17. The 35 pulses are further supplied via an M-divider 47 to a ring winch 49, the output of which forms a further part of the address, on which the weighing numbers and measuring path numbers are stored. At address 79ÖÏ2I12 * ^ PHN 9535 .. ......... _ ....................

10 seerde geh.eugen.cel in geheugen 17h zijn maximaal twee weeg-getallen en vier meetwegnummers opgeslagen. De weeggetallen zijn voor de coördinaten (n,^f) en (n,~jp) identiek, maar moeten wel met convolutiegetallen vermenigvuldigd worden, 5 die elk bij een andere meetweg horen. De vier meetwegnummers worden via pulsen, waartoe het geheugendeel 17b rechtstreeks met de klok 41 is verbonden, een voor een via de verbinding 50 doorgegeven aan het geheugendeel 17a. In het geheugendeel 17a zijn de bij de meetwegnummers behorende M convolutiegetallen opgeslagen en worden dus een voor een toegevoerd aan een vermenigvuldiger 51» Verder worden de weeggetallen simultaan met de convolutiegetallen aan de vermenigvuldiger 51 toegevoerd. Het produkt van een weegge-tal en een convolutiegetal wordt aan een opteller 53 toege-^ voerd. De opteller 53 sommeert de van de vermenigvuldiger 51 afkomstige produkten en ontvangt daartoe rechtstreeks pulsen van de klok 41 om de optellingen synchroon met de toevoer van de weeggetallen en convolutiegetallen aan de vermenigvuldiger 51 te laten verlopen. Verder worden de 30 pulsen van klok 41 aan een tweedeler 55 toegevoerd, die telkens na twee klokpulsen, een puls afgeeft naar een deel van het geheugen 17c. Door toevoeren van een puls aan het geheugendeel 17c wordt de inhoud van de opteller 53 in. het geheugendeel 17c opgeslagen. Het adres, waarop de inhoud oe van teller 53 wordt opgeslagen, wordt gevormd door de uitgangen van de ^-teller 45 en van de ringteller 49, waarbij bij een eerste puls van de tweedeler 55 de inhoud op het adres (n,p) wordt opgeslagen, en bij een tweede puls van de tweedeler 55 op het adres (n,-^). Van het adres (n,^) 30 wordt de eerste term (n) door de tellerstand van de ringteller 49 en de tweede term door de tellerstand van deψ-teller 45 bepaald. Het zal nu duidelijk zijn met welk doel de onderscheidenjke delers (43, 47, 55) worden toegepast.A maximum of two weighing numbers and four measuring path numbers are stored in memory 17h. The weighing numbers are identical for the coordinates (n, ^ f) and (n, ~ jp), but must be multiplied by convolution numbers, 5 which each belong to a different measuring path. The four measuring path numbers are transmitted one by one via the connection 50 to the memory part 17a via pulses, for which the memory part 17b is directly connected to the clock 41. In the memory part 17a, the M convolution numbers associated with the measuring path numbers are stored and are thus applied one by one to a multiplier 51. Furthermore, the weight numbers are simultaneously supplied to the multiplier 51 with the convolution numbers. The product of a weighing number and a convolution number is fed to an adder 53. The adder 53 sums the products from the multiplier 51 and receives pulses directly from the clock 41 for this purpose to cause the additions to run synchronously with the supply of the weighing numbers and convolution numbers to the multiplier 51. Furthermore, the pulses from clock 41 are applied to a two-divider 55, which delivers a pulse to a part of the memory 17c each time after two clock pulses. By adding a pulse to the memory portion 17c, the contents of the adder 53 become in. the memory part 17c is stored. The address at which the content oe of counter 53 is stored is formed by the outputs of the counter 45 and of the ring counter 49, the content of the address (n, p) being supplied with a first pulse of the two-divider 55. stored, and at a second pulse from the two-divider 55 at the address (n, - ^). Of the address (n, ^) 30, the first term (n) is determined by the counter reading of the ring counter 49 and the second term by the counter reading of the counter 45. It will now be clear for what purpose the different dividers (43, 47, 55) are used.

Verder zal het duidelijk zijn dat eerst alle elementen 35 van een ring 40^ (zie fig. 2) worden doorlopen, waarna wordt overgestapt naar een aan de binnenzijde grenzende ring 4C>n . Van de ringteller 49 worden derhalve niet de normale Q-uitgangen als tellerstand genomen, maar deTTo'm- 7905282 *. * 1 1 PHN 9535 plemexitaire uitgangen Q. Zodra Q de waarde 1 bereikt, wordt de^-teller k5 gedurende enige tijd (4.M/2 pulsen) geblokkeerd, zodat een bijdrage in de absorptie in bet centraal element 40o slecbts eenmaal wordt uitgerekend.It will further be apparent that first all elements of a ring 40 ^ (see fig. 2) are run through, after which it is switched to a ring 4C> n adjacent to the inside. It is therefore not the normal Q outputs of the ring counter 49 that are taken as the counter position, but the TT0m-7905282 *. * 1 1 PHN 9535 extension outputs Q. As soon as Q reaches the value 1, the ^ counter k5 is blocked for some time (4.M / 2 pulses), so that a contribution in the absorption in the central element 40o is calculated only once .

5 Ook is in een inkrementteller 54 voorzien, die bij het bereiken van de tellerstand Q = 1 met een stap wordt achteruitgezet, zodat de matrix van elementen die in het geheugendeel 17 c is opgeslagen over een hoek wordt verdraaid (ten opzichte van de in het geheugendeel 17 4 op-10 geslagen matrix met weeggetallen en meetwegnummers) zodat de volgende groep convolutiegetallen, die bij een volgende bronpositie zijn bepaald, kunnen worden verwerkt. De in het geheugendeel 1? c opgeslagen bijdragen moeten nog in de verschillende "ringen 40" over verscheidene elementen 15 worden gesommeerd alvorens een absorptieverdeling met een uniforme informatiedichtheid kan worden weergegeven.An increment counter 54 is also provided, which is set back by one step when the counter position Q = 1 is reached, so that the matrix of elements stored in the memory part 17c is rotated by an angle (with respect to the memory part 17 4-to-10 matrix with weighing numbers and measuring path numbers), so that the next group of convolution numbers determined at a next source position can be processed. The one in the memory section 1? c stored contributions must still be summed in the different "rings 40" over several elements 15 before an absorption distribution with a uniform information density can be displayed.

Het in fig. 4 weergegeven deel van een bewerkings-inrichting 16 (uit fig. 1) is ten dele gelijk aan het in fig. 3 weergegeven blokschema. Zoals in fig. 4 te zien 20 wordt een bijdrage in de absorptie voor een bepaald element niet rechtstreeks aan het geheugendeel 17 c toegevoerd. De door de opteller 53 gemaakte sommaties worden aan een poortschakeling 57 toegevoerd, waarmee de bijdragen in de absorptiewaarde in de elementen worden geschei-25 den in bijdragen in de elementen (n,JJ>) en (n,-^) die elk respektievelijk aan sommatieschakelingen 59+ en 59 worden toegevoerd. De sommatieschakelingen 59+> sommeren de berekende bijdragen van een aantal naburige elementen in een ring (40^, fig· 2), zodat een uniforme informatievel dichtheid kan worden verkregen voor het aantal sommaties in de somraeerschakelingen 59+» afhankelijk te maken van een ringgetal n dat de straal van de ring indiceert. De uitgang van de ring-teller 49 is derhalve verbonden met een doodgeheugen 61 (ROM), waarin een tabel is opgeslagen, 35 die aangeeft hoeveel elementen € (in feite de voor de elementen £ berekende bijdragen) in een bepaalde ring moeten worden gesommeerd. Verder is een vier-deler-puls 7905282 PHN 9535 . ... . .... .......................The part of a processing device 16 (from Fig. 1) shown in Fig. 4 is partly identical to the block diagram shown in Fig. 3. As seen in Fig. 4, a contribution to the absorption for a particular element is not directly applied to the memory portion 17c. The summations made by the adder 53 are applied to a gate circuit 57, which separates the contributions in the absorbance value in the elements into contributions in the elements (n, JJ>) and (n, - ^) each of which are respectively summation circuits 59+ and 59 are applied. The summation circuits 59+> sum the calculated contributions of a number of neighboring elements in a ring (40 ^, fig. 2), so that a uniform information sheet density can be obtained for the number of summations in the summation circuits 59+ »to make dependent on a ring number n that indicates the radius of the ring. The output of the ring counter 49 is therefore connected to a death memory 61 (ROM), in which a table is stored, which indicates how many elements € (in fact the contributions calculated for the elements £) are to be summed in a given ring. Furthermore, a four divider pulse is 7905282 PHN 9535. .... .... .......................

12 teller 63 op de klok 41 aangesloten, die telt hoeveel bijdragen na een begin berekend zijn. Zodra de tellerstand van de vier-deler-puls-teller 63 overeenstemt met een door de teller-stand van ring-teller k9 opgeroepen aantal ele-5 menten, zal de comparator 65 een puls afgeven, die op de ingangen c van de sommeerschakelingen 59+ en 59 worden aangeboden. De op dat moment aanwezige sommen van alle tot dan toe berekende bijdragen, worden afzonderlijk in een deel 17 c van het geheugen 17 opgeslagen. De adressen, 10 waarop de sommen worden opgeslagen, worden bepaald door enerzijds de afzonderlijke ingangen (voor +ψ en -jP) in het geheugendeel 17 c en anderzijds door de tellerstand van een adresteller 67, die het aantal door de comparator 65 afgegeven pulsen telt.12 counter 63 is connected to clock 41, which counts how many contributions have been calculated after a start. As soon as the counter position of the four-divider pulse counter 63 corresponds to a number of elements called up by the counter position of ring counter k9, comparator 65 will emit a pulse which is applied to the inputs c of the summing circuits 59. + and 59 are offered. The sums present at that time of all contributions calculated up to that point are stored separately in a part 17c of the memory 17. The addresses at which the sums are stored are determined on the one hand by the separate inputs (for + ψ and -jP) in the memory part 17c and on the other hand by the counter reading of an address counter 67, which counts the number of pulses delivered by the comparator 65 .

^ De uitgang van de comparator 65 is verder met een terugstelingang van de vier-deler-puls-teller 63 verbonden, zodat na bepalen van de som van een aantal bij elkaar horende bijdragen de teller 63 weer op tellerstand "O” wordt gezet.The output of comparator 65 is further connected to a reset input of the four-divider pulse counter 63, so that after determining the sum of a number of associated contributions, the counter 63 is reset to counter position "0".

20 Van de ring-teller k9 wordt evenals in fig. 3 de complementaire uitgangen Q benut. Zodra de uitgangen Q de waarde 1 voeren, wordt via de verbinding 70 de^-teller geblokkeerd, zodat de bijdrage in het centraal element k00 (zie fig. 2) slechts eenmaal wordt berekend. De ring-teller 25 h9 wordt door de toch door de klok 41 gegenereerde pulsen in een uitgangstellerstand gebracht (buitenste ring van de matrix), hetgeen door een comparator 71 wordt gedetekteerd en die via een of-poort 73 een extra puls aan de P -teller 45 toevoert. Dit heeft tot gevolg, dat de adressen die 30 door de^-teller h5 in het geheugendeel 17 b worden aangestuurd, allen met een inkrement zijn verhoogd. Het doel hiervan is om de polaire matrix met weeggetallen en meet-wegnummers te verdraaien ten opzichte van de matrix in het geheugendeel 17 c, zodat een volgende groep van convolutie-35 getallen, die bij een volgende meting na verdraaiing van de röntgenbron zijn bepaald, kan worden teruggeprojecteerd.As in Fig. 3, the complementary outputs Q of the ring counter k9 are used. As soon as the outputs Q enter the value 1, the counter is blocked via the connection 70, so that the contribution in the central element k00 (see Fig. 2) is calculated only once. The ring counter 25 h9 is brought into an output counter position (outer ring of the matrix) by the pulses generated by the clock 41, which is detected by a comparator 71 and which supplies an extra pulse to the P via an or gate 73. counter 45. As a result, the addresses that are controlled by the counter h5 in the memory portion 17b are all incremented. The purpose of this is to rotate the polar matrix with weighing numbers and measuring path numbers relative to the matrix in the memory part 17c, so that a subsequent group of convolution-35 numbers, which are determined in a subsequent measurement after rotation of the X-ray source, can be projected back.

In het geval van terugprojecteren zoals aan de hand ' 790X212 ~ 13 PHN 9535In case of back projection as per '790X212 ~ 13 PHN 9535

# V# V

van fig·. U is beschreven, wordt de matrix met weeggetallen en meetwegnummers als bet ware met de röntgenbron meegedraaid. In het geval dat is beschreven aan de hand van fig. 3 wordt de matrix van elementen, waarin de groepen 5 van convolutiegetallen worden teruggeprojecteerd, in tegengestelde zin ten opzichte van de stralingsbron gedraaid, waarbij de matrix van weeggetallen en meetwegnummers en de meetwegen zelf (in feite de röntgenbron) als stilstaand worden beschouwd.of fig ·. You have been described, the matrix with weighing numbers and measuring path numbers is rotated, as it were, with the X-ray source. In the case described with reference to FIG. 3, the array of elements in which groups 5 of convolution numbers are projected back is rotated in opposite sense to the radiation source, the array of weighing numbers and measurement path numbers and the measurement paths themselves ( in fact the x-ray source) are considered to be stationary.

10 De in het geheugendeel 17 c opgeslagen absorptie- waarden (na terugprojectie van alle groepen van convolutiegetallen) kunnen via een monitor, die een circulaire of spiraalvormige schrijfstraalafbuiging heeft, worden weergegeven. Uiteraard is een vertaling (slechts één keer nodig) 15 van het polair georganiseerd geheugendeel 17 c naar een cartesisch georganiseerd geheugen mogelijk, waarna weergave van een absorptieverdeling via een conventionele monitor mogelijk is.The absorbance values stored in the memory portion 17c (after backprojection of all groups of convolution numbers) can be displayed via a monitor having a circular or spiral write beam deflection. Naturally, translation (needed only once) of the polar organized memory portion 17c to a Cartesian organized memory is possible, after which display of an absorption distribution via a conventional monitor is possible.

20 25 30 35 ygg-j-yg 220 25 30 35 ygg-j-yg 2

Claims

A method for determining a radiation absorption distribution in a plane of a body, wherein the plane is irradiated with an in-plane, planar, fan-shaped, body-penetrating and comprising beam of radiation 5 in a plurality of equiangular directions and body-passed radiation is detected to generate groups of measurement signals, representing the attenuation of radiation in the body, along groups of measurement paths, each group of measurement paths 10 extending fan-shaped from a node and dividing one fan-shaped beam of radiation, then determining the radiation absorption distribution of each group of measurement signals consecutively with a series of numbers is convolved to obtain a group of convolution numbers 15 and the group of convolution numbers are projected back onto a matrix of elements along roads defined by an associated group of measurement paths, the matrix of elements a vas has a relation to the plane of the body, and a contribution to an absorption value in an element is determined by the product of a convolution number corresponding to the measuring path passing through the element and of a weight that is inversely proportional to the square of the distance between the element and the junction of measuring paths, where an absorption distribution is formed by the absorption values to be built up in the elements by summing the contributions, characterized in that the matrix of the elements is a polar matrix, in which a element is addressed with a distance coordinate and an angular coordinate, the weight numbers being stored in a polar matrix and addressed with the same distance coordinate and the same angular coordinate, which angular coordinate after conversion or backprojection of a group of measurement signals, respectively, of a group of convolution numbers 790 5 2 82. ί · <? 35 9535 (5 is increased by one increment, which is equal to the angular increment between the directions corresponding to the groups to be processed successively.

2. Computer tomography device for determining a radiation absorption distribution in a plane of a body, which device comprises: a radiation source for generating a flat, in-plane fan-shaped beam penetrating the body, a detector device for detecting the body passed radiation 10 and supplying groups of measurement signals, a support frame for the radiation source and the detector device, driving means for moving at least the radiation source in the support frame for irradiating the body from a plurality of directions, the detector device 15 for each direction provides a group of measuring signals representing a radiation attenuation along a group of measuring paths radiating fan-shaped from the radiation source, a processing device for convolving a group of measuring signals with a series of numbers and for projecting back the thus convolved measuring signals for obtaining absorption values, a memory device for storing the absorption values and a display device for displaying the absorption values in elements of a matrix, in which the radiation absorption distribution 25 is displayed, characterized in that the memory device contains a first memory with polar addressing in which the absorbance values can be stored by addressing with a distance coordinate and an angular coordinate and a second memory with a polar addressing, which stores weight numbers for use in the backprojection of the convolution numbers. 35 7905282