METHOD, DEVICE AND MACHINE FOR CALCULATING AN INDEX REPRESENTATIVE OF THE PROPERTIES OF A MATERIAL OF THE BONE TYPE

Abstract

Method, device and machine for calculating an index representative of properties of a material of the bone type of an individual to be subjected to tests, particularly wherein the method includes a first acquisition step for acquiring at least one image having a plurality of elementary units of a sample of the material, wherein a generation step is provided for generating a grid of elementary geometric elements, or cells, which is associated with, in particular superimposed on, the image, an image processing step in which it is provided to calculate at least the apparent elastic modulus and a density coefficient of the material, both as a function of characteristic values of each cell, and a calculation step for calculating the index representative of the properties of a material, wherein the index is a function of the value of the apparent elastic modulus net of the contribution of the density coefficient.

Claims

1. Method for calculating an index representative of the properties of a monophase or multiphase material, wherein the material is a material of the bone type, of an individual to be subjected to tests, the method comprising: an acquisition step for acquiring at least one image of a sample of said material, said image allowing to distinguish from each other the one or more phases object of the study inside a closed region of interest (ROI), in which the subsequent processings will be carried out, a generation step for generating a grid of elementary geometric elements, or cells, which is associated with, in particular superimposed on, said image, a step in which, for each of the cells covering the region of interest (ROI), it is provided to acquire a parameter (I.sub.CELL) indicative of the presence or absence in the cell of the phase considered (for example the bone phase); in the event in which the phase is present, said indicative parameter (I.sub.CELL) is proportional to the quantity of the phase present in the cell, a processing step for processing the image, in which it is provided to calculate a density coefficient (C) of the material as a function of characteristic values of each cell, calculated on the basis of the value(s) of said parameter (I.sub.CELL), of the phases falling within the cell; said density coefficient (C) being, for each region of interest (ROI), a univocal value calculated as the sum of said parameters (I.sub.CELL) involved; and wherein it is provided to calculate the value of the apparent elastic modulus (E*) by means of simulations of the behaviour under load resulting from the application of modes known in the art, in which the mechanical behaviour of the internal structure of the region of interest is simulated, to which each cell contributes with its own value of apparent elastic modulus (E*) and of said parameter (I.sub.CELL), said method comprising a calculation step for calculating said index (BSI) representative of the properties of a material, wherein said index (BSI) is proportional to the value of the apparent elastic modulus (E*) net of the contribution of said density coefficient (C), said method being characterised in that it comprises a rotation step for rotating said image acquired after said image acquisition step and before said acquisition step for acquiring said parameter (I.sub.CELL), said rotation step being a rotation step by such a quantity that an orientation of the structure of the bone tissue is substantially parallel to Cartesian axes of said elementary geometric elements.

2. Method according to claim 1, characterised in that said index (BSI) is calculated according to the following formula:
BSI=f[f.sub.1(E*)f.sub.2(C)] in which f, f.sub.1 and f.sub.2 are functions that are different from each other, and in particular: f is a function of the exposure parameters used during the first image acquisition step; f.sub.1 and f.sub.2 are functions whose result is a constantly positive value such that the result of f.sub.2 is subtracted from f.sub.1 in such a way as to obtain a subtraction of a contribution due to the density of the material to the contribution due to the elastic modulus.

3. Method according to claim 1, characterised in that said index (BSI) is calculated according to the following formula:
BSI=a.sub.1(b.sub.1E*b.sub.2C) where: E* is said apparent elastic modulus (apparent Young's modulus), and C is said density coefficient. a.sub.1 is a value function of the exposure parameters used during said image acquisition step; b.sub.1 and b.sub.2 are positive constants, in such a way as to obtain a subtraction of a contribution due to the density of the material to the contribution due to the elastic modulus.

4. Method according to claim 1, characterised in that in said first acquisition step it is provided to acquire three images of samples to be analysed of the individual to be subjected to tests; it being provided to repeat said processing step for each image acquired thus obtaining at least one value of apparent elastic modulus (E*) and one value of density coefficient (C) for each image; a calculation step of the respective mean values being then provided.

5. Method according to claim 1, characterised in that said phase object of the study is formed by a plurality of elementary units defining a three-dimensional matrix of uprights and crosspieces connected to one another, and in that after said acquisition step an image pre-treatment step is provided in which it is provided to filter said image and to apply said rotation step of said image by such a quantity that a direction of extension of said elementary units is substantially parallel to respective horizontal (X) and vertical (Y) directions that are arranged in such a way as to define the Cartesian axes of the plane (X-Y).

6. Method according to claim 5, characterised in that after said image pre-treatment step a definition step is provided for defining the region of interest (ROI), defined as the largest square or rectangular area that is fully inscribed in the sample of material being analysed.

7. Method according to claim 6, characterised in that said processing step, said image pre-treatment step, said definition step for defining the region of interest (ROI), said calculation step for calculating the respective mean values of the apparent elastic modulus (E*) and of the density coefficient (C), and said calculation step for calculating the index (BSI) are automatically performed by means of image processing algorithms, and/or by means of machine learning techniques of the type with neural networks; the machine learning techniques can comprise both neural networks with reinforcement learning and convolutional neural networks, optionally in conjunction with fuzzy logic techniques.

8. Method according to claim 1, characterised in that said image acquisition step is performed by means of a detection system comprising: X-ray emitting means with which screening means are associated, which are configured to screen a user or an operator from the radiation produced by the X-ray emitting means, sensor means, intended to be associated, in use, with the sample of material to be analysed, for example with the chosen anatomical site of interest of the individual to be subjected to tests, in order to detect the X-rays emitted by the X-rays emitting means.

9. Method according to claim 8, characterised in that said X-ray emitting means are an X-ray generation device with no radioactive source but generating X-rays following an activation command.

10. Method according to claim 8, characterised in that said image acquisition step is performed by means of a detection system consisting of a portable device provided with a handle for an operator, which is connected to a body within which said X-ray emitting means are arranged and with which said screening means are associated, the detection system further comprising a support base configured to support in a stable way said portable device, at least during said image acquisition step, said sensor means comprising an X-ray digital sensor built in a flexible band wrapped around one finger of the individual to be subjected to tests, in order to detect the X-rays emitted by the X-ray emitting means.

11. Method according to claim 8, characterised in that said image acquisition step is performed by means of a detection system consisting of a measuring machine provided with a frame comprising screening walls wherein the measuring machine comprises said X-ray emitting means, the measuring machine comprising a support for resting one hand of the individual to be subjected to tests, the measuring machine comprising said sensor means intended to cooperate with the X-ray emitting means to acquire said at least one image of a sample of said material, the sensor means being made in the form of a sensor plate for detecting the X-rays emitted by the X-ray emitting means.

12. Detection system for detecting images of a sample of material to be analysed suitable to be used by an operator when performing one or more steps of the method according to claim 1, characterised in that said detection system is made in the form of a portable device provided with a handle for the operator, which is connected to a body within which X-ray emitting means are arranged and with which screening means, configured to screen the operator from the radiation produced by the X-ray means, are associated, a support base configured to support in a stable way said portable device, at least during said acquisition step for acquiring the images of the sample of material to be analysed, and sensor means, of the digital type, intended to be associated with the sample of material to be analysed, for example with the chosen anatomical site of interest in case of bone material, in order to detect the X-rays emitted by the X-ray emitting means.

13. System according to claim 12, characterised in that said support base comprises at least one seat or a recess, shaped and positioned in such a way as to be able to receive at least one of either said handle or said body resting on it, so as to set the position of the portable device in a stable way.

14. Detection system for detecting images of a sample of material to be analysed suitable to be used by an individual to be subjected to tests if necessary with the help of an operator when performing one or more steps of the method as in claim 1, characterised in that said detection system is made in the form of a measuring machine provided with a frame comprising screening walls, wherein the measuring machine comprises said X-ray emitting means, the measuring machine comprising a support for resting one hand of the individual to be subjected to tests, the measuring machine comprising said sensor means intended to cooperate with the X-ray emitting means to acquire said at least one image of a sample of said material, the sensor means being made in the form of a sensor plate for detecting the X-rays emitted by the X-ray emitting means.

15. Detection system according to claim 14, characterised in that the support comprises hand positioning means selected from: positioning lines drawn on the support wherein the positioning lines are flat; light positioning lines lighting up the support by means of backlighting or projection; positioning lines arranged on the support wherein the positioning lines are raised; resting surfaces having shapes corresponding to the hand's shape; a band for fixing the hand in the optimal position on the support.

16. Detection system according to claim 13, characterised in that the support is steady and is positioned inside the measuring machine and said sensor plate for detecting the X-rays is integrated on the surface of said support.

17. Detection system according to claim 13, characterised in that the support is fixed to the measuring machine in a sliding way by means of sliding guides in such a way that the support realizes a carriage that enables a movement of the support between a first position that is pulled-out with respect to an insertion slit and a second position of insertion in which the support is almost fully or fully inserted in the measuring machine.

18. Detection system according to claim 17, characterised in that said sensor plate for detecting the X-rays is integrated on the surface of said support.

19. Detection system according to claim 17, characterised in that the insertion slit is protected by a screening cover.

20. Detection system according to claim 14, characterised in that the measuring machine includes a processing unit for calculating an index representative of the properties of the material.

21. Detection system according to claim 14, characterised in that the measuring machine includes connection means for the remote connection with an operator.

22. Detection system according to claim 14, characterised in that the measuring machine includes interface systems for interfacing with the individual to be subjected to tests.

23. Detection system according to claim 22, characterised in that the measuring machine includes one or more interface systems selected from: an interactive touch-screen display; control interface in the form of a keyboard or push-button panel; a push-button for starting measuring; a document reader for reading documents in the form of identification cards of the individual to be subjected to tests; a payment card reader by means of insertion or wireless detection, to perform service payment operations; a scanner for acquiring a medical prescription or a code authorising the required service.

24. Detection system according to claim 23, characterised in that the measuring machine includes a camera inside the machine that frames the measuring zone for displaying on a display to check the correct positioning of the hand on the support.

25. Detection system according to claim 24, characterised in that said display for displaying the measuring zone framed by said camera is said interactive touch-screen display.

26. Detection system according to claim 14, characterised in that the measuring machine includes one or more visual inspection windows protected by a lead-sealed transparent screen through which it is possible to check the correct positioning of the hand on the support.

27. Detection system according to claim 14, characterised in that the measuring machine includes a printer for issuing a medical report.

28. Use of a method according to claim 9, for calculating the index (BSI) representative of the properties of a bone material to predict, and prevent, possible alterations of the bone material of an individual to be subjected to tests; wherein said image acquisition step provides to acquire three digital radiographic images in anterior-posterior projection, by means of the detection system, of the first proximal epiphysis respectively of the three central fingers (index, middle, ring finger) of the non-dominant hand of an individual to be subjected to tests, in correspondence of which the trabecular structure comprises trabeculae oriented according to planes parallel to each other and to the directrix (X-Y) and staggered in such a way as to be arranged in succession after one another and connected by further trabeculae mainly arranged along a third direction (Z).

Description

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS

[0066] In the following an embodiment is described with reference to the drawings enclosed, which are to be considered as a non-limiting example of the present invention, in which:

[0067] FIG. 1 is a block diagram of an embodiment of a method according to the present invention;

[0068] FIG. 2 is a detail view concerning a processing step included in the block diagram of FIG. 1;

[0069] FIG. 3 and FIG. 5 are axonometric views of two examples of structures that can be analysed by means of the method according to the present invention;

[0070] FIG. 4 and FIG. 6 are front views of an image acquired in an acquisition step by means of a detection system according to the present invention of the structures as in FIG. 3 and FIG. 5;

[0071] FIG. 7, FIG. 8, FIG. 9 are schematic front views of a patient's hand arranged to perform a first image acquisition step of the method according to the present invention, wherein one can see sensor means associated with three different fingers, respectively;

[0072] FIG. 10, FIG. 11, FIG. 12 are an example of radiologic images of a portion of the fingers of FIG. 7, FIG. 8, FIG. 9, respectively, acquired by means of the detection system according to the present invention;

[0073] FIG. 13, FIG. 14, FIG. 15 are a further example of radiologic images of a portion of the fingers of FIG. 7, FIG. 8, FIG. 9, respectively, acquired by means of the detection system according to the present invention;

[0074] FIG. 16 and FIG. 17 are schematic perspective views of two embodiment variants of a detection system according to the present invention suitable to implement a method according to FIG. 1 by using a measuring system made in the form of a portable device.

[0075] FIG. 18, FIG. 19, FIG. 20, FIG. 21 are schematic perspective views of two embodiment variants of a detection system according to the present invention suitable to implement a method according to FIG. 1 by using a measuring system made in the form of a machine of the fixed type.

DETAILED DESCRIPTION OF THE INVENTION

[0076] In the figures (FIG. 1, FIG. 2) one can see a block diagram illustrating the various steps of the method according to the present invention. The method comprises an acquisition step (10) for acquiring at least one image i, preferably of a series of images i wherein i indicates an integer corresponding to an index of each of the images acquired. Each of the images acquired is a representative image of the structure to be analysed.

[0077] In a preferred embodiment the images can be acquired by means of a detection system that can be (FIG. 16, FIG. 17) a portable measuring device (20) or (FIG. 18, FIG. 19, FIG. 20, FIG. 21) a measuring machine (1), wherein the portable measuring device (20) and the measuring machine (1) operate in accordance with what is disclosed in the present invention, which will be described in further detail hereinafter with particular reference to the figures.

[0078] In an embodiment (FIG. 1) the method comprises a validation step (11) in which it is provided to validate the image or images acquired in the acquisition step (10). In a preferred embodiment the validation step (11) of the images provides to check the matching of the images with respect to reference standards, for example at least in terms of format and size.

[0079] Afterwards, the method comprises (FIG. 1, FIG. 2) a processing step (12) for processing the images acquired in the acquisition step (10), after the optional validation step (11).

[0080] According to an embodiment the processing step (12) comprises (FIG. 2) a pre-treatment sub-step (12A) for pre-treating the image acquired, a definition sub-step for defining the ROI (12B), which is the acronym for Region Of Interest, in which the subsequent processings will be carried out, and a calculation sub-step (12C), which will be described in further detail hereinafter.

[0081] According to an embodiment in the pre-treatment sub-step (12A) for pre-treating the image the image is rotated after being suitably filtered by using, for example, a sub-threshold erosion non-linear filter according to techniques known in the art.

[0082] According to an embodiment in the definition sub-step for defining the ROI (12B) the region of interest is defined as the largest square area that can be fully inscribed in the acquired image of the structure to be analysed. According to an alternative embodiment the region of interest has predetermined and constant dimensions. In an embodiment the region of interest can be rectangular and substantially coincide with the most external perimeter of the acquired image of the structure to be analysed.

[0083] Since the region of interest is generally square or rectangular, in the pre-treatment sub-step (12A) for pre-treating the image acquired the image is rotated in such a way that the two main axes of the structure to be analysed, that is to say, the Cartesian axes X and Y, which are orthogonal with respect to each other, are parallel to respective adjacent sides of the region of interest, arranged perpendicularly to each other.

[0084] According to an embodiment in the calculation sub-step (12C) the apparent elastic modulus E.sub.x*, E.sub.y* is calculated in the directions X and Y, respectively, according to what is disclosed in U.S. Pat. No. 7,386,154 B2. Thus, it is provided to add up the apparent elastic moduli E.sub.x*, E.sub.y* in the two directions X and Y.

[0085] According to an embodiment in the calculation sub-step (12C) the density coefficient (C) is also calculated, which, in the case of bone material, is indicative of the local level of mineralisation of the bone phase in the region of interest, and is directly proportional to the sum of the levels of grey of the image acquired.

[0086] In a preferred embodiment the density coefficient (C), too, is calculated according to what is disclosed in the above-mentioned U.S. Pat. No. 7,386,154 B2 as the summation of the values of the indices (I.sub.CELL) assigned to each cell of the model, for all the cells.

[0087] According to embodiments said pre-treatment sub-step (12A), ROI definition sub-step (12B), calculation sub-step (12C) can be alternatively performed manually or automatically.

[0088] After performing the processing step (12), the method (FIG. 1) provides a calculation step for calculating mean values (13), averaged over the number of images acquired. In fact, it is evident that, if several images are acquired, for each of them it is provided to repeat the validation step (11) and the processing step (12), then calculating the average of the parameters calculated.

[0089] According to embodiments it is provided to calculate the mean value respectively of the sum of the apparent elastic moduli E.sub.x*, E.sub.y* and of the density coefficient (C), which have been previously calculated.

[0090] It is evident that, in case only one image has been acquired, the calculation step for calculating the mean values (13) can be omitted.

[0091] After performing the calculation step for calculating the mean values (13), the method then provides a calculation step (14) for calculating the index representative of the properties of the materials.

[0092] According to an embodiment the representative index, herein called BSI, which is the acronym for Bone Structure Index, is calculated according to the following formula:

BSI=f[f.sub.1(E*)f.sub.2(C)]

in which f, f.sub.1 and f.sub.2 are functions that are different from each other, not necessarily linear, and in particular: [0093] f is a function of the exposure parameters used during the first image acquisition step; [0094] f.sub.1 and f.sub.2 are functions whose result (for each of them) is a constantly positive value so that the result of f.sub.2 is actually subtracted from f.sub.1.

[0095] It should be noted that, since f.sub.1 and f.sub.2 are constantly positive functions, the formula for calculating the index BSI allows to subtract the contribution due to the density of the material to the contribution due to the elastic modulus.

[0096] According to a simplified embodiment the representative index, herein called BSI, which is the acronym for Bone Structure Index, is calculated according to the following formula:

BSI=a.sub.1(b.sub.1E*b.sub.2C)

where:

[0097] a.sub.1 is a value function of the exposure parameters used during the first image acquisition step, while b.sub.1 and b.sub.2 are positive constants,

[0098] E* is the apparent elastic modulus (apparent Young's modulus), and

[0099] C is the density coefficient.

[0100] It should be noted that, since b.sub.1 and b.sub.2 are positive constants, the formula for calculating the index BSI allows to subtract the contribution due to the density of the material to the contribution due to the elastic modulus. In an embodiment, after the calculation step (14) for calculating the representative index, a comparison step is optionally provided, in which the calculated index BSI is compared with a reference value. For example, in the case of bone tissue, the index BSI is compared with the mean structural index of the young population and/or with the value of the mean structural index of the population in the same age group as the patient.

[0101] According to embodiments the above-described steps can be implemented by a processing unit (16), of the type known in the art, optionally arranged remotely with respect to the place in which the patient's images are acquired or (FIG. 18) integrated in a measuring machine (1).

[0102] According to embodiments one or more of said pre-treatment sub-step (12A), ROI definition sub-step (12B), calculation sub-step (12C) and the calculation steps for calculating the mean values (13) and for calculating the representative index (14) can be performed automatically by image processing algorithms, and/or by visual learning techniques, for example by means of machine learning techniques of the type with neural networks.

[0103] In an embodiment the machine learning techniques can comprise both neural networks with reinforcement learning and convolutional neural networks, optionally in conjunction with fuzzy logic techniques.

[0104] According to an embodiment said pre-treatment sub-step (12A), ROI definition sub-step (12B), calculation sub-step (12C) and the calculation steps for calculating the mean values (13) and for calculating the representative index (14) can be designed according to a self-learning logic with neural networks at different levels of abstraction, organised through convolutional layers followed by a pooling layer, thanks to which the system can automatically learn functional relations between input data and output data, so as to be able to operate without having to resort to specially designed characteristics.

[0105] As previously explained, the invention also relates (FIG. 16, FIG. 17, FIG. 18) to a detection system made according to the present invention, which can be made in the form of a (FIG. 16, FIG. 17) portable measuring device (20) or of a (FIG. 18) measuring machine (1). Both the portable measuring device (20) and the measuring machine (1) will be suitable to acquire the images as explained with reference to the image acquisition step (10). While in the case of the measuring machine (1) the processing unit (16) can be integrated in the machine itself, in the case of the portable measuring device (20) the processing unit (16) can be an external device to which the portable measuring device (20) is to be connected or with which the portable measuring device (20) communicates by means of a respective wireless communication channel.

[0106] The detection system made in the form of a portable measuring device (20) comprises (FIG. 16, FIG. 17) a support base (21), on which the portable measuring device (20) rests in a stable way. The portable measuring device (20) comprises a handle (22) for the operator and a body (23), within which the X-ray emitting means are arranged. It should be noted that, in the embodiments illustrated (FIG. 16, FIG. 17), the support base (21) is shaped in such a way as to receive both the handle (22) and the body (23) resting on it. To this purpose, the support base (21) comprises a seat (24) for the handle (22).

[0107] The portable measuring device (20) also comprises (FIG. 16, FIG. 17) a screening element (25), such as a screening plate, which for example is circular. The screening element is shaped and oriented in such a way as to be able to screen the operator from the X-rays emitted by the X-ray emitting means during the image acquisition step (10).

[0108] According to an embodiment the central axis of symmetry of the screening element (25) is substantially aligned with a longitudinal axis of development of the body (23) of the portable measuring device (20).

[0109] In a different embodiment (FIG. 17) of the support base (21), the latter comprises a recess (26), shaped and positioned in such a way as to be able to receive an end portion of the body (23) resting on it. The detection system made in the form of a portable measuring device (20) also comprises X-ray detecting means, for example configured as a sensor element (27) to be applied on one finger of the individual to be subjected to tests, which can be built in a flexible band that can be wrapped around one finger and can be fastened in a stable way by releasable fastening means, such as adhesive means or tear-off fastening means.

[0110] Likewise (FIG. 18, FIG. 19, FIG. 20, FIG. 21), the detection system made in the form of a measuring machine (1) comprises X-ray emitting means as well. In the case of the measuring machine (1) the screening from radiation can be obtained by means of screening walls (3) fixed to a frame (2) of the machine. The screening walls (3) screen the outside of the machine in all directions, also providing flexible or movable and adaptable screening means for screening a slit (5) adapted for the insertion of the patient's hand to be examined. The screening of the measuring machine (1) is sized and oriented in such a way as to enable the screening of the individual being subjected to tests and of the environment surrounding the measuring machine (1) itself with respect to the X-rays emitted by the X-ray emitting means during the image acquisition step (10). The measuring machine (1) also comprises X-ray detecting means, for example configured as a sensor plate (19) that can be built in a support (28) of the measuring machine (1). In an embodiment the support (28) is fixed and is positioned inside the measuring machine (1). In a preferred embodiment of the present invention (FIG. 18, FIG. 19, FIG. 20, FIG. 21) the support (28) is slidingly fixed to the measuring machine (1) by means of sliding guides (29) in such a way that the support (28) realizes a carriage (17) that enables a movement of the support (28) between a first position (FIG. 19, FIG. 20) that is pulled-out with respect to an insertion slit (5) and a second position (FIG. 21) of insertion in which the support (28) is almost fully or fully inserted in the measuring machine (1). The insertion slit (5) is protected by a screening cover (18). By this solution it is possible to facilitate the positioning of the hand (37) of the individual to be subjected to tests, because the support (28) is fully visible in the first pulled-out position (FIG. 19, FIG. 20) and the hand's positioning can be facilitated by means of positioning lines (35) drawn on the support (28), it also being possible to resort to solutions of positioning lines (35) that are flat drawn lines or raised lines or resting surfaces whose shapes are adapted to the hand's shape. One can also provide solutions in which the positioning lines (35) are light lines lighting up the support (28) by means of backlighting or projection, for example by means of a laser projector. In this case, too, one can provide a band for fixing the hand in the optimal position on the support (28).

[0111] As far as both the detection system made in the form of a portable measuring device (20) and the detection system made in the form of a measuring machine (1) are concerned, the X-ray emitting means emit radiation that crosses the sample to be analysed, and of which a radiographic image has to be taken, and that is detected by the sensor element (27) or by the sensor plate (19), respectively.

[0112] With particular reference to the detection system made in the form of a portable measuring device (20), it needs an operator to be present and to handle the portable measuring device (20), and the operator must be a trained medical operator qualified for the use of the portable measuring device (20), such as a doctor or a radiology technician. Furthermore, the portable measuring device (20) must be necessarily used inside a closed room, although no specific room screening measure is required due to the tool's low emission.

[0113] However, there is the need in the market for large-scale availability at low costs of a detection system that is able to implement a method, which will be described in further detail hereinafter, for assessing the bone's trabecular component for the purpose of preventing and monitoring the risk of fragility fracture. In order to meet such requirement, it is necessary to provide an automatic tool able to eliminate the uncertainties and the variabilities connected to the operator. The detection system made in the form of a measuring machine (1) meets such requirement, because the measuring machine (1) can be easily positioned in different venues with respect to a medical or radiology centre. For example, the measuring machine (1) can be easily installed in easily accessible places spread over the territory, such as chemist's shops, OTC pharmacies or other places. Therefore, the measuring machine (1) is an integrated detection system able to automatically perform the operations of correct acquisition of the images and of calculation of the quality index of the bone structure. Such an apparatus can be directly activated by an operator locally with the operator present at the apparatus, or remotely with the operator remotely controlling the apparatus, or in a fully automatic way directly by the individual to be subjected to tests, optionally with remote control by an operator who checks the correctness of the operations performed by the individual to be subjected to tests. The method for calculating the quality index of the bone structure, BSI, will be explained in detail in the following of the present description. In particular, the value of the index BSI can be obtained, like in example 2, from the acquisition of the three images of the proximal epiphyses of the three central fingers of the non-dominant hand, or, in a simplified version, from the image of the proximal epiphysis of only one finger of the non-dominant hand.

[0114] In this embodiment the image acquisition system is enclosed in the measuring machine (1).

[0115] For example, the measuring machine (1) may look like an ATM machine. In order to obtain an image having a resolution suitable for reading, parameters compatible with the described measuring protocol and with the reference protocols for the sector will be provided.

[0116] The measuring machine (1) has a frame (2) provided with the previously described screening means in the form of screening walls (3). For example one can provide screening walls (3) able to contain the dispersion of the radiation produced on the inside, such as walls of polycarbonate charged with lead or metal with lead-sealed screening, depending on the position of the apparatus and on the presence or non-presence of an operator.

[0117] The measuring machine (1) is further provided with interface systems (4, 6, 8, 9, 15) for interfacing with the user.

[0118] For example, the measuring machine (1) can comprise an interactive touch-screen display (4). The measuring machine (1) can comprise a keyboard integrated in the touch-screen display (4) or an external control interface (6) can be provided, such as an external keyboard or push-button panel, for inputting the required data and for reading the information from the system to the individual to be subjected to tests who is being examined. For example, one can provide the input of anthropometric data or the confirmation of a non-pregnancy status or of the presence of risk factors.

[0119] A specific push-button (39) can be provided for starting measuring or measuring can be started by means of a control on the touch-screen display (4).

[0120] Furthermore, it will be possible to acquire images from a camera (38) inside the body of the machine framing the measuring zone for the display on the touch-screen display (4) to check the correct positioning of the hand (37) on the support (28). For example, one can provide a recognition step for recognising the hand's position in order to identify a non-correct positioning with respect to the sensor plate (19) and to communicate, by displaying them on the display (4), instructions for obtaining the optimal positioning for image acquisition.

[0121] Alternatively or in combination, one can also provide one or more visual inspection windows (7) protected by a lead-sealed transparent screen through which it is possible to check the hand's correct positioning.

[0122] For example the measuring machine (1) can comprise a document reader (8) such as a reader of a health insurance card or of a tax code card or of an identification card, able to univocally identify the individual to be subjected to tests being examined and to identify other individuals, if necessary, such as any operators or maintenance technicians.

[0123] For example, the measuring machine (1) can comprise a card reader (9) by means of insertion or wireless detection, such as by RFID or NFC technology, to perform service payment operations.

[0124] For example, the measuring machine (1) can comprise a scanner (15) for acquiring a medical prescription or a code authorising the required service.

[0125] As previously explained, the measuring machine (1) comprises a support (28) optionally realized on a carriage (17) sliding on guides (29). The carriage can be moved automatically by means of an electric motor (40).

[0126] The support can comprise positioning lines (35) or can be suitably shaped to facilitate the hand's correct positioning. One can also provide solutions for lighting up or backlighting the support (28) in order to facilitate the hand's (37) correct positioning. The support comprises a digital sensor plate (19) for acquiring images, of such dimensions as to ensure the acquisition of the image of the whole zone to be analysed (one or three proximal epiphyses).

[0127] Inside the structure of the measuring machine (1) one can also find: [0128] An X-ray generation device, with settings predetermined as per protocol. The device does not contain on its inside any radioactive source, but the required field is generated only for the moment strictly required for acquisition (in the current protocol 0.17 seconds); [0129] A processing unit (16), provided with a network interface card for remote connection for diagnostics or for the connection between the individual to be subjected to tests and an operator who connects remotely for example for an interview or for controlling remote operations or for transmitting the measuring results; [0130] An optional camera (38) or, alternatively or in combination, inspection windows (7) optionally in combination with a mirror; the purpose of such devices is to check the hand's correct positioning on the support; [0131] An optional printer (42) for issuing a medical report.

[0132] The patient accesses the structure and identifies himself/herself through his/her health insurance card or tax code card. The patient inputs the required data and then places his/her hand on the support with the help of the dotted lines in order to position the hand on the sensor correctly. In an embodiment the authorised technician is present and, after checking the hand's correct positioning through the display/camera, by means of a push-button activates the X-ray apparatus by generating the emission of the minimum radiation required for taking the X-ray scan. The acquired image is displayed on the display and, after being approved by the technician, is processed for generating the medical report.

[0133] In a further embodiment the authorised technician performs the described operations remotely.

[0134] In a further embodiment all the operations described are performed in a fully automatic way.

[0135] It will be clear that it is possible to modify or add parts or steps to the above-described method and system, without departing from the scope of the present invention.

[0136] It will also be clear that, although the present invention has been described with reference to a specific example, a person skilled in the art may realize many other embodiments for calculating an index representative of the mechanical properties of a material, all included in the scope of the present invention.

[0137] Two examples of application of the method according to the present invention will now be described and will have to be considered as illustrative both as far as the detection system made in the form of a portable measuring device (20) is concerned and as far as the detection system made in the form of a measuring machine (1) is concerned.

Example 1

[0138] In the figures (FIG. 3, FIG. 5) one can see two different three-dimensional structures (30, 30) to be analysed by means of the method according to the present invention. The illustrative three-dimensional structures comprise a first three-dimensional structure (30) and a second three-dimensional structure (30).

[0139] It should be noted that the two different three-dimensional structures (30, 30) occupy the same overall volume and consist of a very similar number of elementary (geometric) units, and for this reason they are characterised by approximately equal density values but, as will be explained, their mechanical characteristics are much different from each other. This is due to the different arrangement in space of the structural elements, or elementary geometric units, included in the structure.

[0140] Each three-dimensional structure (30, 30) is made up of a plurality of elementary units (31, 32, 33) defining a three-dimensional matrix of uprights and crosspieces connected to one another. For example, a first elementary unit (31), a second elementary unit (32) and a third elementary unit (33) can be provided. According to an embodiment each three-dimensional structure (30, 30) comprises a plurality of first elementary units in the form of uprights (31), which extend parallel to a vertical direction Y, and a plurality of second elementary units (32) in the form of crosspieces and third elementary units (33) in the form of crosspieces. The crosspieces comprise: [0141] a first group of second elementary units (32) in the form of crosspieces, which extend substantially parallel to a horizontal direction X (width); [0142] a second group of third elementary units (33) in the form of crosspieces, which extend substantially parallel to a third direction Z (depth).

[0143] The directions X, Y and Z form a set of three Cartesian axes. Between the uprights (31) and the crosspieces (32, 33) in the three-dimensional structure (30, 30) a plurality of free spaces (34) is defined, that is to say, empty spaces in which the material is absent.

[0144] The above-described method according to the present invention allows to reconstruct a two-dimensional image (FIG. 4, FIG. 6) of the three-dimensional structures (30, 30) in which the different portions of the three-dimensional structure are depicted with different levels of grey, depending on their density. It should be noted that, as indicated by the dotted line in the figures (FIG. 4, FIG. 6), in this example the ROI substantially coincides with the external perimeter of the three-dimensional structures (30, 30), according to the projection on the plane X-Y.

[0145] As it is known, a radiographic image allows to reproduce on a two-dimensional image also the contribution, in terms of density, of the phase considered arranged on the subsequent planes, which are parallel and arranged after the first one. In other words, the image of the figures (FIG. 4, FIG. 6) is inscribed in the plane X-Y and the different levels of grey are a function of the development of the structural elements in the plane Z, that is to say, in the depth direction.

[0146] According to what is disclosed in U.S. Pat. No. 7,386,154 B2 a grid of cells, for example triangular, is superimposed to the image as in the figures (FIG. 4, FIG. 6), as it is schematically indicated in some portions of the figure (FIG. 6). A characteristic parameter (I.sub.CELL) is assigned to each cell according to the phase present in the space occupied by the cell itself. For example, if a cell is arranged over one of said elementary units (31, 32, 33) or said free space (34), the characteristic parameter (I.sub.CELL) of the cell will be equal to the assigned value (i.e., for example, 0, 0.2, 0.4, or 1). If a cell bridges two or more of said elementary units (31, 32, 33) or said free space (34), the characteristic parameter (I.sub.CELL) will be calculated as the weighted average of the phases falling within that cell, with weights proportional to the surface of the cell occupied by the respective phase.

[0147] More in detail, in the example described herein, light grey indicates that, in correspondence of that portion, there is one of said elementary units (31, 32, 33), in particular a structural element (33) extending in an orthogonal direction with respect to the figure (parallel to the direction Z) and, therefore, it is a portion completely full with the phase considered. According to the method, a value of the characteristic parameter (I.sub.CELL) amounting to 1 is associated with such colour. Smoke grey indicates the superimposition, one behind the other on planes parallel to the directrix X-Y and staggered along Z, of two of said elementary units (31, 32, 33) or structural elements, be they uprights (31) or crosspieces (32). According to the method, a value of the characteristic parameter (I.sub.CELL) amounting to 0.4 is associated with such colour. Black indicates the presence of only one of said elementary units (31, 32, 33) or structural elements, be it an upright (31) or a crosspiece (32). According to the method, a value of the characteristic parameter (I.sub.CELL) amounting to 0.2 is associated with such colour. Finally, according to the method, a value of the characteristic parameter (I.sub.CELL) amounting to 0 is associated with each free space (34), in correspondence of which the phase is completely absent, such free spaces not contributing to the mechanical characteristics of the structure.

[0148] The results of the third processing step (12) of the method according to the present invention provide the following.

[0149] As far as the first three-dimensional structure (30) is concerned (FIG. 3), one obtains an apparent elastic modulus (E*) amounting to 202 MPa, and an apparent elastic modulus (E*) amounting to 228 MPa calculated on the basis of the model constructed by the above-described cell method starting (FIG. 4) from the two-dimensional image. With such values the index BSI, calculated according to what is disclosed in the present invention, amounts to 174.

[0150] As far as the second three-dimensional structure (30) is concerned (FIG. 5), one obtains an apparent elastic modulus (E*) amounting to 251 MPa, with a 20% variation with respect to the previous structure, and an apparent elastic modulus (E*) amounting to 289 MPa calculated on the basis of the model constructed by the above-described cell method starting (FIG. 6) from the two-dimensional image, with a 21% variation with respect to the previous structure.

[0151] With such values the index BSI, calculated according to what is disclosed in the present invention, amounts to 225, with a 23% variation with respect to the previous structure, greater than the variation detected with reference to the elastic modulus. Therefore, the method according to the present invention allows to emphasize the differences between structures having analogous density, by means of a numerical index that is readily and easily interpreted and that facilitates the classification thereof.

Example 2

[0152] In this example the method according to the present invention is used to assess the structural resistance of the bone tissue of an individual to be subjected to tests.

[0153] In the first acquisition step (10) three images of respective anatomical sites of an individual to be subjected to tests are acquired. In this example the chosen anatomical sites are the first phalanges of the index (FIG. 7), middle (FIG. 8) and ring finger (FIG. 9), respectively. In this case (FIG. 7, FIG. 8, FIG. 9) the left hand of an individual to be subjected to tests is shown, wherein, in correspondence of the first phalanx the digital sensor (27), built in a flexible band wrapped around the finger, is arranged.

[0154] Three radiographic images of the first proximal epiphysis of the three central fingers (index, middle, ring finger) of the non-dominant hand of an individual to be subjected to tests are then acquired, in accordance with what has been previously disclosed in the present invention.

[0155] In an embodiment the acquired images are radiograms, for example acquired by means of a detection system according to the present invention, which can be a detection system made in the form (FIG. 16, FIG. 17) of a portable measuring device (20) or made in the form (FIG. 18, FIG. 19, FIG. 20, FIG. 21) of a measuring machine (1). As a non-exhaustive example, the portable measuring device (20) can be the device called NOMAD Pro 2 marketed by Aribex, or the device called EzRay Air marketed by Vatech.

[0156] According to embodiments during the first acquisition step (10) the exposure parameters are kept constant and are related to the type of device used.

[0157] In an embodiment the digital sensor (27) can be a sensor of the known type, such as the sensor called GXS700 marketed by Gendex or EZ Sensor Classic Slim marketed by Vatech.

[0158] Two sets of three radiographic images acquired in the first acquisition step (10) are shown in FIG. 10, FIG. 11, FIG. 12 and FIG. 13, FIG. 14, FIG. 15, respectively.

[0159] In each of these Figures the square box indicates the region of interest in which the planned processing is carried out. It should be noted that this is the largest square area that is fully inscribed in the trabecular zone of the bone phase being analysed.

[0160] We would like to point out that, by suitably enlarging the ROIs, it would be possible to highlight the bone structures in which the trabeculae are substantially oriented like the elementary units (31, 32, 33) or structural elements and like the free spaces (34) of the previously discussed example (FIG. 3, FIG. 4, FIG. 5, FIG. 6). In other words, in these anatomical sites, the trabeculae are oriented on parallel planes (lying on the plane directrix X-Y) reciprocally connected by trabeculae mainly arranged in a direction that is orthogonal to these planes, the trabeculae extending parallel to the third direction Z (depth). The method advantageously provides a rotation step for rotating the image(s) acquired after the image acquisition step and before the parameter (I.sub.CELL) acquisition step. The rotation step consists of a rotation by such a quantity that an orientation of the structure of the bone tissue is substantially parallel to Cartesian axes forming said elementary geometric elements. In this way one of the plane directrices X-Y, according to which the trabeculae are oriented, coincides with the axes of the elementary units (31, 32, 33) or structural elements, in particular with the axes of uprights (31) and crosspieces (32, 33), in the structure (30, 30).

[0161] In this case, too, like in the above-described previous example, the structural elements forming the structure comprise portions depicted with different levels of grey, which, as already said, take into account the quantity of the considered phase arranged in the depth direction Z. As stated above, a different value of the characteristic parameter (I.sub.CELL) is associated with each level of grey, by re-processing which according to the method disclosed in the above-mentioned U.S. Pat. No. 7,386,154 B2 shortly mentioned again above, it is possible to obtain the results below.

[0162] As far as the bone structure as in FIG. 10, FIG. 11, FIG. 12 is concerned, one obtains an apparent elastic modulus (E*) amounting to 588 MPa, which is the result of the mean of the six values E.sub.x*, E.sub.y*, calculated for the three images acquired. With such values the index BSI, calculated according to what is disclosed in the present invention, amounts to 154.

[0163] As far as the bone structure as in FIG. 13, FIG. 14, FIG. 15 is concerned, one obtains an apparent elastic modulus (E*) amounting to 604 MPa, which is the result of the mean of the six values E.sub.x*, E.sub.y*, calculated for the three images acquired. With such values the index BSI, calculated according to what is disclosed in the present invention, amounts to 194.

[0164] The first set of radiographic images was acquired from the non-dominant hand of a female individual to be subjected to tests, aged eighty-three, whose bone structure is essentially weak due to old age.

[0165] The second set of radiographic images was acquired from the non-dominant hand of a female individual to be subjected to tests, aged twenty-five, whose bone structure is resistant as expected from the individual's young age.

[0166] The method according to the present invention allows to emphasize the quality of the bone structure by means of a numerical index that is readily and easily interpreted.

[0167] In fact, it should be noted that, against a difference by about 2% between the apparent elastic moduli (E*) calculated on the basis of the first and of the second set of radiographic images, the corresponding index BSI shows a difference by about 20% in both cases.

[0168] The present invention has been described with reference to the figures enclosed in a preferred embodiment thereof, but it is evident that many possible changes, modifications and variants will be readily understood by a person skilled in the art in the light of the previous description. Thus, it should be noted that the present invention is not limited to the present description, but it includes any changes, modifications and variants in compliance with the appended claims.

NOMENCLATURE USED

[0169] With reference to the identification numbers in the enclosed figures, the following nomenclature has been used: [0170] 1. Measuring machine [0171] 2. Frame [0172] 3. Screening wall [0173] 4. Display [0174] 5. Slit [0175] 6. Control interface [0176] 7. Inspection window [0177] 8. Document reader [0178] 9. Card reader [0179] 10. Image acquisition step [0180] 11. Image validation step [0181] 12. Image processing step [0182] 12A. Image pre-treatment sub-step [0183] 12B. ROI definition sub-step [0184] 12C. Calculation sub-step for calculating the elastic module and characteristic coefficient [0185] 13. Calculation step for calculating the mean value of BSI [0186] 14. Calculation step for calculating the index BSI [0187] 15. Scanner [0188] 16. Processing unit [0189] 17. Carriage [0190] 18. Cover [0191] 19. Sensor [0192] 20. Portable measuring device [0193] 21. Support base [0194] 22. Handle [0195] 23. Body [0196] 24. Seat [0197] 25. Screening means [0198] 26. Recess [0199] 27. Sensor [0200] 28. Support [0201] 29. Guide [0202] 30. First three-dimensional structure [0203] 30. Second three-dimensional structure [0204] 31. First elementary unit [0205] 32. Second elementary unit [0206] 33. Third elementary unit [0207] 34. Free space [0208] 35. Positioning lines [0209] 36. Movement direction [0210] 37. Hand [0211] 38. Camera [0212] 39. Push-button [0213] 40. Motor [0214] 41. X-ray generation device [0215] 42. Printer