FR2768888A1 - Data processing method for - Google Patents

Abstract

The data is divided into subsets with size varying inversely with complexity. For each basic subset, two or more working subsets are used and a group of parameters is selected from a finite set of such groups, to test and optimize selection criteria values by use of a multilinear function which associates each working subset with at least one parameter.

Description

 The present invention relates, in general, to a method and a device for processing data, and in particular for compressing and decompressing image data.

 More particularly, the present invention relates to the transmission (that is to say the remote representation) and / or the storage (that is to say the representation in a memory) of data, with or without compression of those -this. These data can advantageously consist of digital or analog data of an image or sequence of images, and / or digital or analog data relating to sound (music, speech, etc.), and / or digital data or analog signals relating to any mono or multidimensional signal.

Before setting out the objectives and means of the invention, it is proposed below to give the following definitions:
"Data set": any data set representing physical quantities (most often voltage levels) which may themselves represent other physically perceptible or measurable quantities. In privileged applications, the data sets concerned are images, sound or mono or multidimensional signals, etc. In the application concerning image processing, we sometimes aim at the "image" instead of the "data set" relating to the image. What is explained below with regard to "images" or "sub-images" is respectively applicable to "data sets" or "data subsets" and vice versa.

 "Data representation": any "processing" of a set of data of a given type, resulting in the perfect or imperfect transformation of said set of data, into another type. Example: the data may consist of the 256 gray levels of the pixels of a monochrome image, their representation consisting of a series of binary words capable of being stored or transmitted; conversely, the representation of the data constituted by said binary words can consist of their transformation to find data of the same type as the initial data.

 "Primary representation" by convention, any processing resulting in the transformation of data of a first type to data of a second type. This is called "primary processing" of the data, such as image compression.

 "Secondary representation": by convention, the transformation of data of the second type resulting from primary processing. We speak in this case of "secondary processing", such as for example an image decompression.

 "Tertiary data": this is data resulting from secondary processing, such as for example an image decompression. We observe that the tertiary data represent, perfectly or not, the primary data from which they come.

 Various methods of image representation, with compression, using the so-called "fractal" technique are known in the prior art.

 Document US-A-5,065,447 discloses a method and a device allowing image processing or storage with compression of the data relating to the initial image. In this method, the data relating to the initial image are divided, that is to say that the image is divided into a plurality of elementary sub-images (called "domain blocks"). This method then consists in generating an ordered dictionary of a set of reference sub-images (called "range blocks") constituted from image portions of predetermined size having undergone a certain number of predetermined operations such as rotation and reversal. along various axes of symmetry. Then, for each elementary sub-image, a comparison is made with the set of reference sub-images of the dictionary and the most similar reference sub-image of the elementary sub-image considered is selected from the dictionary. Finally, the method consists in processing or storing parametric information relating to the addresses of the reference sub-images selected in the dictionary, to represent each of the original elementary sub-images.

 It is by this set of operations that the method described in this document makes it possible to obtain a first representation of the image with data compression.

 To carry out, from said parametric information, the restitution of the initial image, the method and the device described in this document, perform the following operations: starting from an arbitrary starting image, steps similar to those are carried out above, that is to say that we partition the arbitrary starting image and that we constitute a dictionary from the elementary sub-images thereof. However, the dictionary is only partially constituted by performing for each sub-image only the predetermined operation corresponding to the relevant address of the dictionary for the position of the sub-image considered.

 The data thus obtained is used to repeat this process until the difference between two consecutively obtained images is less than a predetermined threshold. The image obtained last is then an approximation of the initial image.

 This method, which is interesting from a theoretical point of view, has the major drawback that it is sometimes difficult to practice on an industrial level with the means known to date. Indeed, each image must give rise to a long analysis with the constitution of a particularly important dictionary, the elementary sub-images of the image to be stored or to be transmitted being all compared with each of the reference sub-images present in the dictionary. The inventors, who carried out simulations, thus found that for certain images of medium size and resolution (512 x 512 pixels at 3 x 256 color levels) the processing time for compression was of the order of 1000 2000 seconds on a workstation, or around half an hour. Such processing time is obviously prohibitive for almost any industrial application.

 Document WO-A-93/17519 discloses another method and device for processing images with compression of the data relating to the initial image. In this process, the data relating to the initial image are also partitioned, that is to say that the image is divided into a plurality of elementary sub-images (called "domains"). In this process, we then determine a set of fractal transforms linking the sub-images to the set of initial data, that is to say that we search in the initial image for a portion of image similar to the image. elementary considered, so as to minimize the error between the initial data and the approximate data thus obtained. The coefficients of the set of fractal transformations thus produced constitute the first representation of the initial image.

 To restore the image, the process and the device described in this document search for the fixed points of the elementary fractal transformations whose coefficients have been received without however explaining how this search is carried out. The initial image is thus restored by bringing together the restored sub-images.

 This method is also interesting from a theoretical point of view.

However, the inventors have found that processing times can be very long, especially during data compression.

 The inventors have sought to develop a data processing method of the same kind as the prior art methods analyzed above, but having better performance both in terms of the quality of the data restored and in terms of processing speed.

 The methods and devices according to the invention solve this problem.

In general, the present invention relates, according to a first aspect, to a method of primary processing of a set of so-called "primary" data representative of physical quantities to determine a set of so-called "secondary" data representative of the primary data, process characterized in that it comprises:
a step of partitioning the primary data set into so-called "elementary"subsets;
- For each elementary subset, a step of determining at least two so-called "working" subsets of the primary data set, each associated, in a predetermined manner, with said elementary subset;
- for each elementary subset, an operation of selecting a group of parameters, from a predetermined set of groups of parameters, during which:
one determines for each group of parameters of said
predetermined set, the value of a criterion of
selection of a parameterized combination of sub
work packages associated with said subset
elementary, each subset of work being
associated, by said combination, with at least one said
parameter of said group of parameters, and
o we select the group of parameters which optimizes the
value of said selection criterion, the parameter values which optimize the value of the selection criterion constituting secondary data representative of the primary data of the elementary subset considered.

Jointly, according to a second aspect, the present invention relates to a device for primary processing of a set of so-called "primary" data representative of physical quantities to determine a set of so-called "secondary" data representative of the primary data, method characterized in that 'it comprises:
- primary data entry,
a means of partitioning the primary data set into so-called "elementary"subsets;
- suitable calculation means:
for each elementary subset, to be determined
at least two so-called "working" subsets of
the primary data set, associated, each, of
predetermined way, said subset
elementary,
# for each elementary subset, to be selected
a group of parameters, among a set
predetermined groups of parameters, by determining,
for each group of parameters in said set
predetermined, the value of a selection criterion of a
parametric combination of working sub-assemblies
associated with said elementary subset, each sub
work set being associated, by said
combination, to at least one said parameter of said group
parameters, and selecting the group of
parameters which optimizes the value of said criterion of
selection, the parameter values which optimize the value of the selection criterion constituting secondary data representative of the primary data of the elementary subset considered.

 Thanks to these provisions, because the working sub-assemblies are associated in a predetermined manner, the method according to the invention does not have to test all the possible working sub-assemblies, as in the prior art. This process can therefore work very quickly.

In addition, since:
- at least two working sub-assemblies are associated with each elementary sub-assembly and
- at least one parameter of the combination concerns each work subset, the variety of combinations of the work subsets associated in a predetermined manner with each elementary subset is very important and makes it possible to fairly represent said elementary subset .

 As, moreover, by choosing a limited number of groups of parameters from which a test and optimization procedure of the value of a selection criterion is carried out, the speed of the primary processing is further increased.

 Finally, the parameters being independent of the dimension of the elementary subset with which they are associated, a reconstruction can be carried out with the same simplicity and with the same efficiency whatever the dimensions of the reconstructed image.

 According to particular characteristics, the method as succinctly explained above comprises an operation of determining at least one value representative of a difference between said combination and the data of the elementary subset, said representative value of difference constituting a representative secondary value primary data of the considered subset.

 Thanks to these provisions, the secondary data represent the primary data more precisely.

According to a third aspect, the present invention relates to a method of secondary processing of a set of so-called "secondary" data representative of physical quantities to determine a set of so-called "tertiary" data, characterized in that it comprises:
a step of defining a so-called "initial" tertiary data set having the same number of data as the tertiary data set to be determined,
a step of partitioning the initial data set into so-called "elementary" subsets and allocating secondary data to each elementary subset,
- For each elementary subset, a step of determining at least two so-called "working" subsets of the initial data set, each associated, in a predetermined manner, with said elementary subset; and
- iteratively and, during each iteration, successively for each elementary subset:
a sub image determination operation
work packages associated with said subset
elementary, by a parameterized multilinear application
by secondary data assigned to said sub
elementary set, to constitute a new sub
elementary set, the number of parameters of said
multilinear application associated with said subset
elementary, being at least equal to the number of sub
work packages associated with said subset
elementary.

Jointly, according to a fourth aspect, the present invention relates to a device for secondary processing of a set of so-called "secondary" data representative of physical quantities to determine a set of so-called "tertiary" data, characterized in that it comprises:
- a secondary data entry,
a means of defining a so-called "initial" tertiary data set having the same number of data as the tertiary data set to be determined,
a means of partitioning the initial data set into so-called "elementary" subsets and allocating secondary data to each elementary subset,
- suitable calculation means:
r for each elementary subset, to be determined
at least two so-called "working" subsets of
the initial, associated tertiary data set,
each, in a predetermined manner to said subset
elementary; and
to be done, iteratively and, during each
iteration, successively for each subset
elementary:
a sub image determination operation
work packages associated with said subset
elementary, by a multilinear application
parameterized by assigned secondary data
audit elementary subset, to constitute a
new elementary subset, the number of
associated multilinear application parameters
said elementary subset, being at least equal
the number of associated work subsets
audit elementary subset, and
- an output of tertiary data consisting of the data of the elementary subsets resulting from said iterations.

 Thanks to these provisions, it is a method of successive approximations which provides the tertiary data.

 It is observed here that, when all the elementary subsets have the same dimensions, a compromise must be made. Indeed, when these dimensions are too small, the compression rate is limited since it is directly a function of the number of pixels which are represented by the five parameters associated with each elementary subset. On the other hand, when the dimensions of the elementary subsets are large, the quality of the reconstructed image after a secondary processing of these parameters, may be too limited for a comfortable exploitation of the reconstructed image, the smallest details may have been poorly processed. during primary treatment.

 In an image, it is common for certain areas to have very little detail and to be practically uniform and for other areas of the image to be very complex (complexity which can be represented by the variance of the values of each component of the pixels of the area considered).

 The invention proposes, in addition to the characteristics set out above, a method making it possible to use elementary sub-assemblies of large dimensions in the most uniform areas of the image and elementary sub-assemblies of smaller dimension in the less uniform image areas. To this end, for each elementary subset being processed, a complexity criterion is determined and an estimation of the partition opportunity relating to the complexity of representation of the primary data in the elementary subset is carried out. When this estimate does not meet the complexity criterion, that is to say when the representation of the elementary subset considered is supposed to present a certain complexity, it is partition born into elementary subsets of smaller dimension which are processed in turn, unless they have a minimum dimension.

The characteristics, advantages and other objects of the invention will moreover emerge during the description of the appended drawings in which
- Figure 1 schematically illustrates the arrangement of an elementary subset Kj and four working subsets L1, L2, L3 and L4, each of which is associated with it, each dimension of which is twice that of the sub -associated elementary assembly, and which are offset in both directions, with respect to said elementary subset, in a predetermined manner,
FIG. 2 is a simplified block diagram of a particular embodiment of the device for primary processing of a set of data,
FIG. 3 is a simplified block diagram of a particular embodiment of the device for secondary processing of a set of data,
FIG. 4 is a simplified flowchart for compressing a color image of 512 × 512 pixels used in the device of FIG. 2,
FIGS. 5 and 7 are respectively flow diagrams illustrating certain operations of the flow diagram of FIG. 4 and a first of its variants,
FIG. 6 is a flow diagram illustrating certain operations of the flow diagram of FIG. 5,
FIG. 8 is a simplified flowchart for decompressing a color image implemented in the device of FIG. 3,
FIG. 9 is a flow diagram illustrating an operation of the flow diagram of FIG. 8,
FIG. 10 is a flow diagram illustrating an operation of the flow diagram of FIG. 8,
FIG. liA is a representation of a succession of partitions of elementary subsets,
FIG. 11B is a representation of the structure of partition information and of parameters representing this succession of partitions,
FIG. liC is a representation of the organization of the parameters and of the partition information resulting from the primary processing of a single image,
FIG. il D is a representation of a second variant of the organization of the parameters and of the partition information resulting from the primary processing of a single image,
FIG. 11 E is a representation of the parameters and of the partition information resulting from the primary processing and actually used during a secondary processing, according to the second variant illustrated in FIG. 11D, and
- Figure 12 is a representation of an organization of the parameters and partition information resulting from primary processing of several images.

 FIG. 1 illustrates how, in this first embodiment, each working sub-assembly 4 is linked to at least one given elementary sub-assembly K1. The elementary subset Kj, of size T, having the coordinates (xj, y,) is illustrated in bold lines, while the corresponding working subsets are illustrated in dotted lines.

 The values of the coordinates of the lower left corners of the elementary subsets constitute the set of pairs of multiples of 32 less than 511. The dimensions of the elementary subsets are 32 pixels in each direction of the image. The set of elementary subsets therefore constitutes a partition of the image into 16 x 16 elementary subsets.

 It is observed that, in this first embodiment, to each of the subsets Kj are linked four working subsets Li (numbered L1 ... L4) whose dimensions are each twice that of the subset elementary considered and whose centers are offset by 16 pixels in both directions, compared to that of the elementary subset considered. The expressions below show the coordinates of a subset Kj and the working subsets L, ... L4.

Kj: [xj, xj +31] x [yi, yi + 31]
L1: [xj - 32, xj + 31] x [y1, y + 63]
L2: [x1 xi + 63] x [yi, yi + 63]
L3: [xi - 32, xi +31] x [yj - 32, yj +31]
L4: [xi, xi + 63] x [y1-32, yj +31]
It is noted that, in this embodiment, the division of the image into working sub-assemblies constitutes a covering by slates (in English "tiling").

 The points of the work subsets defined above which do not belong to the image, that is to say those whose at least one of the coordinates is negative, or whose at least one of the coordinates is greater to 512, are fictitiously associated with an average value, for example 127, for each component.

 We will now describe, with the support of FIG. 2, a first embodiment of a device for processing a set of data representative of physical quantities, and intended to allow a primary representation of said set by the implementation of the process of the invention.

 In this embodiment, which relates to a device for processing color images, the image data in 512 x 512 pixels are accepted as input (referenced 591, in FIG. 2), while at output (referenced 581 , FIG. 2), the device delivers the parameters aj, i = 1, ..., 4 and the value b and partition information, relating to each of the subsets constituting a partition of the image.

 It is noted here that the present invention is certainly not limited to a particular image format but that it allows the processing of images of all sizes. However, in an embodiment described and shown, the processed and restored image is an image of 512 x 512 pixels at 256 levels on each of the three components Y, luminance, Cr, red chrominance, and Cb, blue chrominance.

 In the block diagram of FIG. 2, the coding device has the general reference 500.

It comprises, connected to a data and address bus, designated under the reference 501:
- a controller (ALU + CTR) 550 incorporating an arithmetic logic unit. This controller is constituted, in this embodiment, by a microprocessor of the type of that marketed by the company INTEL under the reference i486,
a random access memory RAM, bearing the reference 510. This memory includes a certain number of registers which are described below and a rapidly addressable part (cache) 514,
a program memory 530, of ROM type,
- a first direct access memory register 560 (REG 1),
- a second direct access memory register 570 (REG 2),
- an input buffer (IBFR) 590, connected to an input connector 591 and
- an output buffer (OBFR) 580, connected to an output connector 581.

 In the ROM 530 is recorded a program allowing the processing of the input data 591, with compression of these in accordance with the method of the invention. The compression program, which relates to FIG. 2, the reference 531, is described below in support of FIGS. 4, 5, 6, 7, 11A, 11B, liC, nest and 12.

 All the means 510, 530, 531, 550, 560 and 570 here constitute an embodiment of "first logical transformation means" within the meaning of the invention.

In this embodiment, the RAM 510 includes, among other things, the following registers intended for the implementation of the program 531:
- a register 511 (IM), intended to store the image data,
- a counter 513 (CNT) used as explained below (Figure 4).

 - part 514 (CACHE) of the RAM. This part is a part of very quickly addressable memory, in which are, as is known, stored the most commonly used data at a given instant.

The (CACHE) part 514 of the memory 510 comprises the following registers:
- registers f (L1) f (L4> referenced 516 to 519, in which the values [f (L1)] (x, y) corresponding to each pixel of coordinates (x, y) of the elementary subset are recorded, these values being obtained as detailed later, from values of several pixels of the working subset Li,
a register 520 intended to store the selected parameters at, ..., a4 and the calculated value b, relating to an elementary subset considered Kj, and
- a register 522 (T) intended to keep the size of the current elementary subset Kj,
the register 521, intended to keep an intermediate variable i identifying the references of elementary subsets forming an initial partition of an image to be compressed,
- registers 523 (Dmin) and 524 (D) which store distortion measurements of the representation of elementary subsets,
- registers 525 and 526 intended to keep the coordinates xj and yj, respectively, of the lower left corner of the current elementary subset,
registers 540 to 543 respectively containing variables "height ~ strip", "ideal height", "height ~ total" and "nbligneSlues" used as described below,
- register 560 (REG 1), intended to cumulate, for each elementary subset, the parameters a1 ..., a4 selected and the calculated value b, as well as any partition information, relating to each working subset ,
- Register 570 (REG 2), intended to store the list of addresses (list ~ under ~ images) of all the elementary subsets Kj, as well as their size (T).

Connected to the input means 591, an uncompressed data source 592 can include many means within the reach of those skilled in the art. For example, it may consist of a hard disk or floppy drive or compact disc (known as CD-ROM). The image data to be compressed can therefore advantageously be saved on the corresponding media: hard disk, floppy disk or CD
ROM (compact disc). This recording can be done in any manner known to those skilled in the art.

 The source of uncompressed data can also be constituted by a video interface, capable of delivering the data of the image to be compressed at the input 591. This source can also be constituted by means of reception of the data connected to a network of transmission, the source being capable of transforming the data received into an adequate format so that they can be presented in different definition or resolution. These data at the output of the secondary processing device are called "tertiary data", whether they correspond to an image of the same resolution as the image having undergone the primary processing or to an image having another resolution.

 In the block diagram of FIG. 3, the iterative rendering or calculation device bears the general reference 600.

It comprises, connected to a data and address bus, designated under the reference 601:
- a controller (ALU + CTR) 650 incorporating an arithmetic logic unit. This controller is constituted, in this embodiment, by a microprocessor of the type of that marketed by the Company INTEL under the reference i486;
a random access memory RAM, bearing the reference 610, which includes a certain number of registers which are described below and a rapidly addressable part (called "cache memory") 614;
- a program memory 630, of ROM type. In this memory, there is recorded in particular a program 631 for decompressing the data of the image presented at input 681 for the purpose of reconstructing an approximation of the initial image of identical or different resolution.

- a first register 560 with direct access (REG 1), and
a second direct access register 570 (REG 2) identical to those described in support of FIG. 2,
- an output buffer (OBFR) 690, connected to an output connector 691, and
- an input buffer (IBFR) 680, connected to an input connector 681.

In this embodiment, the RAM 610 has, among other things, the following registers intended for implementing the program 631:
a register 511 (IM), described in support of FIG. 2,
a register 621 (RE), intended to store the data of an image,
a register (IT) 612 intended to store the current value of the number of iterations carried out to perform the secondary representation,
a register (p) 623 intended to store a number of iterations as explained below (FIG. 8),
- a scale factor z in a register "z" 628.

a counter 513 (CNT) described in support of FIG. 2. This counter is used here as explained below (FIG. 8),
- the part 614 (CACHE) of the RAM 610. This part is a part of very quickly addressable memory, in which are, as is known, stored the most commonly used data at a given instant, the part ( CACHE) 614 of memory 610 includes registers 516 to 520, 540 to 543 and 625 to 627 respectively identical to registers 516 to 520, 540 to 543 and 525, 526 and 522 described in support of FIG. 2,
the register 560 (REG 1), intended to store, for each elementary subset, the selected parameters a1 ..., a4 and the calculated value b, and the possible partition information, associated with this subset, and
- the register 570 (REG 2), intended to memorize the list of addresses of all the elementary subsets Kj.

 Connected to the input means 681, a source of compressed data (not shown) may consist of means for storing the compressed image (hard disk, floppy disk, CD-ROM, etc.) or a simple interface intended for delivering these. This source can also consist of a data reception system connected to a telecommunications network. Means for using decompressed data (not shown), connected to output 691 may consist, for example, of means for storing decompressed data or in a simple interface delivering on its output the decompressed data. They can also consist of a device for transmitting decompressed data over a telecommunications network.

Description of a specific compression flowchart
a 512 x 512 pixel color image
This flow diagram is illustrated in FIG. 4, FIG. 6 being a flow diagram explaining certain operations of the compression flow diagram of FIG. 4.

 In this case, to each pair of values (x, y) (column, row coordinates) taken in the definition domain (in this case in the plane image 512 x 512 pixels) correspond three values f (x, y ) from 256 levels relating to each of the three components Y, Cr and Cb.

We organize three sequences of numbers each corresponding to the level of a component (i.e. 512 x 512 = 262,144 f (x, y) values per component), by scanning the image data as follows:
- the lines are scanned from left to right;
- the lines are incremented from bottom to top.

 These three sequences are recorded here in the uncompressed data source. They can be read on the input port 591 at the command of the microprocessor 550.

In 200A, the registers (IM) 511, read lines and two registers (REG 1) 560 and (REG 2) 570 are initialized to the uniform zero value. During operation 200A, the so-called "ideal" height of the independent strips is also determined, this height corresponding, for example, to a dimension:
smaller than the size of the available part of memory
hide and
corresponding to a multiple of 32 lines of the image, with a
minimum of 64 lines.

 We observe, in particular, that the ideal height can thus be equal to the height of the image (here 512 lines).

The files relating to the components Y, Cr and Cb are accessible separately and are provided in a predetermined order, namely:
- firstly the data of the image considered and relating to the luminance component (Y),
- then the image data relating to the red chrominance component (Cr), and
- finally the image data relating to the blue chrominance component (Cb).

 The 2008 test consists of determining whether the variable n lines read is equal to the total height of the image or not. When the result of the test 200B is negative, the operation 200C is carried out and consists in initializing to zero the value of the counter (CNT) 513, which is used to identify a color component of the image. Then, the operation 200D determines the height of the next stripe constituting a subset of the image to be processed. During this operation, the height of the part of the image remaining to be processed is compared to twice the so-called "ideal" height. When the height of the part of the image remaining to be processed is greater than twice the so-called "ideal" height, the height of the next strip of the image is equal to the so-called "ideal" height. When this is not the case, the height of the next strip is equal to the height of the part of the image remaining to be processed. The value of the height of the next strip to be processed thus determined is loaded into the strip height register and the value of the number of lines read is incremented by the value of strip height.

 Operation 200E consists in incrementing the counter (CNT) 513.

Then, during the test 201, it is determined whether the counter (CNT) 513 is strictly greater than 3, or not.

 When the result of the test 201 is negative, the program branches to step 202, during which the microprocessor 550 will read from the input port 591 the next component of the image strip to be processed. each component, Y, Cr or Cb, is presented as a series of 262,144 values t (x, y) corresponding to the level of the component considered for each of the pixels of the image. A part of the sequence which corresponds to the component being processed and to the image band to be processed is therefore loaded into the register (IM) 511 during step 202.

 The operations 200B to 202 thus constitute a so-called assembly division step during which the primary data set is divided into primary data subsets, here in horizontal bands which are parallel and equal to each other, with the possible exception from the last.

 The program then connects to step 204, during which the initial partition of the band of the image of band height x 512 pixels is carried out into elementary subsets Ki of 32 x 32 pixels, as illustrated in figure 1. This partition is done by determining the coordinates (x, y) of the lower left corners of each of the m subsets K1 considered. Each of the addresses in the lower left corners of the 32 × 32 pixel blocks corresponds to a particular address in the list recorded in the register (IM) 511 of the RAM 510.

 The list of these addresses, each accompanied by the initial size 32, is loaded into the register (REG 2).

 During test 205, it is checked whether the list recorded in the register (REG 2) is empty or not. When the result of test 205 is negative, this means that there are still elementary subsets K1 to be processed and the program branches to step 206.

 In step 206, the coordinates xj, yj of the first elementary subset of the list recorded in the register (REG 2) are deleted, as well as the size T of the elementary subset considered to be copied into the registers 525. , 526 and 522, respectively. The program then connects to step 207 detailed in FIG. 5 or, according to the first variant, in FIG. 7. As illustrated in FIG. 5, during this step 207, a test 255 determines whether the elementary subset considered has the minimum dimensions, that is to say in the embodiment described and shown, 2 times 2 pixels, or not.

When the result of the test 255 is positive, the operation 250 consists in determining the four work sub-images Li associated with the work subset considered, as illustrated with reference to FIG. 1. During this operation 250, with starting from the coordinates kept in registers 525 and 526 of the memory RAM 510 and the size T, it is proceeded to the search of the four neighboring working sub-assemblies L1 ..., L4 (see figure 1 for the sub-assemblies associated with the initial partition). To this end, the rule recalled above is applied:
each working sub-assembly has dimensions twice that of the elementary sub-assembly to which it is linked, that is to say in each direction, the size T; and
each center of a working sub-assembly is offset by a half-dimension (Tl2) of the elementary sub-assembly, in each direction, relative to the center of the associated elementary sub-assembly.

Thus, the expressions below show the coordinates of a subset Ki and the working subsets L1, ..., L4, for all values of size T, the latter taking into account the complementary partitions which have been performed (see operation 256 below):
Ki: [x, xi + T- 1] x [y *, yj + T-1]
L1: [x - T, x1 + T-1] x [y *, y1 + 2T-I]
L2 [x, x1 + 2T-1] x [ys, yj + 2T-1]
L3: [x1 - T, xi + T-1] x [yj - T, yj + T-1]
L4: [xj, xi + 2T-1] x [yj - T, yj + T-1]
The working subsets having the same shape and four times more pixels than the corresponding elementary subset, it is understood that, by homothety, each pixel of Ki corresponds to four pixels of L1 (respectively L2, L3 and L4).

 For each elementary subset Ki considered, the controller 550 reads from the register 570 (REG 2) the coordinates (xj, y,) of the lower left pixel of K1 and the size T of the elementary subset Kj, from the register ( T) 522.

Then, the arithmetic unit 550:
- will search, for the component, Y, Cr or Cb, during processing, and successively for each of the pixels of the elementary subset considered, the values corresponding to each of the 4 pixels of L, (respectively L2, L3 and L4 ) which correspond to it, and
- stores the average of these four pixels in the register f (l,) 516 (respectively flL2) 517, f (L3) 518 and fiL4) 519) of the cache memory 514).

 Following operation 250 (FIG. 5), the program performs operation 251 during which, in general, the parameters a (in this case a1 ..., a4) are selected and the b value is calculated.

 A simplified description of operation 251 (and also of operations 258 and 260 illustrated below) is given in the flow chart represented in FIG. 6.

 As illustrated in FIG. 6, these operations consist in determining the parameters aj from a restricted number of different values, and preferably from a selection of power inverses of 2. During these stages, the parameters ai (in this case a1 ..., a4) are selected and the value b is calculated by measuring a representation defect, here called a distortion D, for each of the combinations of values taken from the restricted choice.

 It is observed here that, in base 2, the multiplication by inverses of powers of 2 returns, to the shift to the right of the binary values.

The use of a limited number of power reversals of 2, therefore makes it possible to carry out multiplication operations very easy to execute and to reduce the duration of the treatments carried out, because these divisions and multiplications by offset avoid the manipulation of numbers. real whose representation, in binary numbers, in memory, is complex and costly.

In the embodiment described and represented in FIG. 6, the restricted choice is limited for each of the four parameters a1 ..., a4 to +1/4 or -1/4. Consequently, sixteen 4-uples are tested, and constitute table 1
a1 a2 a3 a4
+1/4 +1/4 +1/4 +1/4
+1/4 +1/4 +1/4 -1/4
+1/4 +1/4 -1/4 +1/4
+1/4 +1/4 -1/4 -1/4
+1/4 -1/4 +1/4 +1/4
+1/4 -1/4 +1/4 -1/4
+1/4 -1/4 -1/4 +1/4
+1/4 -1/4 -1/4 -1/4
-1/4 +1/4 +1/4 +1/4
-1/4 +1/4 +1/4 -1/4
-1/4 +1/4 -1/4 +1/4
-1/4 +1/4 -1/4 -1/4
-1/4 -1/4 +1/4 +1/4
-1/4 -1/4 +1/4 -1/4
-1/4 -1/4 -1/4 +1/4
-1/4 -1/4 -1/4 -1/4
We observe here that the choice of the values +1/4 and -1/4 is linked to the number (equal to 4) of working subsets which are linked to the elementary subset. Indeed, preferably, according to the invention, for a given number e of work subset associated with a given elementary subset, the restricted choice of the values of the parameters includes at least the values + 1 / e and -1 / e. Additionally, if you want each parameter value to be chosen from two values, so that the data compression rate is maximum, the restricted choice will be limited to the values + 1 / e and -1 / e. In addition, if one wishes to use the multiplication operations which are the simplest and fastest to perform by microprocessor processing means, the value e is chosen from the powers of 2 (a multiplication by the inverse d 'a power of two corresponding, in base two, to a simple shift to the left of the binary data translating the value of a pixel, shift of a number of bits equal to the power of 2 considered).

During operation 910, the intermediate variable Dmin which represents a distortion between the primary data and the multilinear combination on the working subsets, is initialized to the value 2562 = 65.536
During the 911 test, it is determined whether the last 4-uple of all possible combinations of values from the restricted choice of values has been used or not.

 When the result of the test 911 is negative, the operation 912 consists in taking the next set of possible values of the parameters (a1 ..., a4) from the restricted choice described above. When operation 912 is performed for the first time, the first set is taken. Then operation 913 consists in determining the value of b which corresponds to the 4-up current of values of aj, by carrying out the following calculation: b is equal to the sum of the values of the pixels in the elementary subset minus the four products parameters aj by the means of the pixel values of the working subset to which aj is attached, this sum being divided by the number of points of the elementary subset.

 Operation 915 then consists in determining a selection criterion based on the distortion D corresponding to the four parameters determined during operation 913. The distortion D is here determined according to the following formula: D is equal to the sum of the pixel values in the elementary subset minus the four products of the parameters aj by the means of the pixel values of the working subset to which aj relates and minus the calculated value b, the whole being squared and divided by the number of points of the elementary subset.

 It is observed that D is a measure of the error of the representation of the points of the elementary subset by the five parameters considered and more precisely, it is the mean square error.

The test 916 consists in determining whether the value of the distortion D which was determined during the operation 915 is less than the value Dmin previously kept in RAM. When the result of the 916 test is negative, the 911 test is repeated. When the result of the 916 test is positive,
Operation 917 consists in memorizing in random access memory, on the one hand, the parameters (a1, a2, a3, a4) and the calculated value b, on the other hand, the new value of Dmin which takes the value of D determined during operation 915 and, finally, the value D associated with the memorized parameters. Then operation 912 is performed.

When the result of the 911 test is positive, the 919 test consists in reading from memory the selected group of parameters a ,, ..., a4 and the value b and in determining whether the group of parameters is equal, term by term, to group (1/4, 1/4, 1/4, 1/4, 0), or not. Indeed, in the case where an area of the image is uniform or almost uniform, that is to say in the case where each pixel is almost worth a value u, an elementary subset within this area uniform is very well represented by the group of parameters (1/4, 1/4, 1/4, 1/4, 0) if the working subsets to which they relate are also included in this uniform area (ie if all their pixels are also of value u). In this case of uniform area, we effectively observe that each pixel of the elementary subset has a value perfectly equal to:
1/4 f (Le) + 1/4 f (L2) + 1/4 f (L3) + 1/4 f (L4) + 0, i.e.
1/4 (u + u + u + u) 14 + 1/4 (u + u + u + u) 14 + 1/4 (u + u + u + u) 14 + 1/4 (utu + u + u) / 4 + 0 = u.

 The problem appears with the restitution of the data when such parameters are applied to the initial arbitrary image, usually chosen uniformly gray. We immediately understand that this gray value is restored for all the pixels of the elementary subset, and this is perpetuated with each iteration, and does not change if the uniform area is very large (case of a blue sky in a landscape, for example) unless the number of iterations is also very large.

 The inventors observed, moreover, that the parameters (1/4, 1 / 4.1 / 4.1 / 4.0.0) were selected according to the method described almost only in the case of uniform bottoms.

 By replacing the parameters (1 / 4,1 / 4,1 / 4,1 / 4, 0) by (1/4, -1/4, 1/4, -1/4, u), we obtain, on the other hand, a very rapid convergence towards the uniform value u. This is what is done during Operation 914.

 When the result of test 919 is negative, operation 251 (FIG. 5) is finished. When the result of test 919 is positive, operation 914 is carried out.

 Operation 914 consists in replacing the parameters (+1/4, +1/4, +1/4, +1/4) by (+1/4, -1/4, +1/4, -1 / 4) and to determine the values of b and D which correspond to the 4-uple of values of aj: (+1/4, -1/4, +1/4, -1/4) (see operations 913 and 915 above).

 We observe here that, The set of multilinear combinations of the work subsets in the elementary subsets is not a contraction, in the mathematical sense, because some of the multilinear combinations do not meet the convergence criterion that the sum of the multiplicative coefficients (here the parameters a1) of this multilinear combination is, in absolute value, strictly less than 1.

 However, when several iterations of the multilinear combinations are applied to any set of data called "initials", the successive images thus obtained converge on a fixed point. According to the invention, this fixed point is an approximation of a homothety of the primary data set. When the initial data set has the same dimensions as the primary data set, the fixed point is an approximation of the primary data set.

 We also observe here that operation 914 is used to avoid that for large uniform flat surfaces, such as for example blue skies, the convergence of the successive approximation is slowed down, if not erroneously.

 Following operation 914, operation 920 consists in storing in place of the old group of parameters a1 ..., a4 and at the value b considered and the associated distortion, the new group of parameters and the corresponding b value determined during operation 914 and the associated distortion. Operation 251 (Figure 5) is then completed.

 It is observed here that, insofar as the elementary subset has minimal dimensions, no partition information is necessary, since no additional partition can be performed (see operation 256 below).

 It is understood that by implementing the flow diagram illustrated in FIG. 6, each of the parameters a1 ..., a4 can be represented by a single binary datum representing its sign, since it can take only two possible predefined and opposite values .

When the result of test 255 (FIG. 5) is negative, operation 252 consists in measuring the variance of the levels of the component considered, in the elementary subset considered. Next, operation 253 consists in determining a so-called "complexity" criterion applicable to the elementary subset considered. Preferably, the definition of the complexity criterion consists of the definition of a threshold below which the complexity criterion is considered to be met. Preferably also, the threshold value takes into account the dimensions of the elementary subset considered. In the embodiment described and shown, this threshold is the function corresponding to the following table:
dimensions of the subset in threshold pixels
2x2,200
4x4 160
8x8 100
16x 16 60
32 x 32 20
We observe here that it is preferable to use an array corresponding to a decreasing function of the size T.

We note here that according to variants not shown:
the value of the threshold is constant whatever the dimensions of the elementary subsets, or
- the threshold value is an inverse function of a dimension of the elementary subset considered.

 Test 254 determines whether the complexity criterion is met or not. In the embodiment described and shown, it is respected when the measure of the variance of the levels of the component considered is less than the value of the threshold which corresponds, in the table below, to the dimensions of the elementary subset considered.

 When the result of the test 254 is positive, the operation 259 consists in determining the four working sub-images as explained above concerning the operation 250. Then, the operation 258 is identical to the operation 251, if this is only stored, in register 560 (REG 1), successively the partition information equal to "0" since no additional partition has been performed and the parameters selected (a, ..., a4) and the calculated value b. Then, test 205 (FIG. 2) is repeated.

 When the result of test 254 is negative, operation 256 consists in performing a complementary partition of the elementary subset considered. In the embodiment described and shown, this complementary partition consists in dividing each elementary subset into four equal square elementary subsets. To this end, during operation 256, the coordinates of the four points which define the lower left corners of the four elementary sub-sets resulting from the additional partition are added to the list list ~ under ~ images. With these coordinates are associated a size T which is equal to the size of the elementary subset that is partitioned, divided by 2. In addition, partition information equal to "1", since an additional partition has been performed, is stored in register 560 (REG 1).

 Following operation 256, test 205 (FIG. 2) is repeated.

 According to the first variant illustrated in FIG. 7, the complexity criterion used to determine the need to carry out an additional partition is not, as presented with regard to FIG. 5, the comparison, on the one hand, of a measurement statistics carried out on the levels of the component considered in an elementary subset with, on the other hand, a threshold, but the comparison of a measure of difference between, on the one hand, the image obtained after applying the combinations multilinear to the work subsets associated with the elementary subset, with the initial image in this elementary subset.

 This first variant is illustrated in FIG. 7, in which the tests and operations common with FIG. 5 bear the same references.

In this first variant, following operation 206,1, operation 250 is carried out, taking into account the present value of the size T. Operation 250 determines, it will be recalled, the associated working sub-assemblies to the elementary subset considered.

 Then, the operation 260 is identical to the operation 251 (FIG. 5) with the exception of the storage in the memory register 560 (REG 1) which is not carried out. Operation 260 therefore provides multiplicative (a ,, ..., a4) and additive (b) parameters or coefficients of a multilinear combination of the working subsets, in the elementary subset considered as well as a measure of the distortion D. The test 255 follows operation 260 and determines, as explained above, whether the elementary subset has minimal dimensions or not.

 When the result of the test 255 is positive, the operation 263 consists in memorizing in the register 560 (REG 1) the four parameters a1 ..., a4 selected during the operation 260 and the associated value b, without associating them partition information since no new additional partition can take place. Then test 205 (FIG. 4) is repeated.

When the result of test 255 is negative, operation 262 consists in determining a so-called "complexity" criterion applicable to the elementary subset considered. Preferably, the definition of the complexity criterion consists of the definition of a threshold below which the complexity criterion is considered to be met. Preferably also, the threshold value takes into account the dimensions of the elementary subset considered. In the embodiment described and shown, this threshold is the function corresponding to the following table:
dimensions of the subset in threshold pixels
2x2,800
4x4 640
8x8,400
16 x 16 240
32 x 32 80
We observe here that it is preferable to use an array corresponding to a decreasing function of the size T.

 Test 254 then determines whether the complexity criterion is met or not. In the embodiment described and shown, it is respected when the measurement of the mean quadratic difference between the component considered and its image by the parameterized multilinear combinations, on the elementary subset considered, is less than the value of the threshold which corresponds, in the table above to the dimensions of the elementary subset considered. This mean square deviation corresponds, in fact, to the distortion D.

 When the result of test 254 is negative, operation 256 consists in performing a complementary partition of the elementary subset considered. In the embodiment described and shown, this complementary partition consists in dividing each elementary subset into four equal square elementary subsets. To this end, during operation 256, the coordinates of the four points which define the lower left corners of the four resulting elementary subsets are added at the head of the list of lists ~ images in register 570 (REG 2). of the additional partition. With these coordinates is associated a size T which is equal to the size of the elementary subset that is partitioned, divided by 2. Finally, a partition information equal to "1" is stored in the register 560 (REG 1) since an additional partition has been performed.

 When the result of test 254 is positive, during operation 261, the four parameters a1 ..., a4 selected during operation 260 and the associated value b are stored in register 560 (REG 1) with partition information equal to "0" since no additional partition has taken place.

 In the first variant exposed in FIG. 7, it is advantageous to measure the average of the values of the mean square deviation of all the elementary subsets treated for which this deviation is measured during operation 254, to determine the usefulness of additional treatment during secondary treatment: when this average is greater than a predetermined overall critical value, it is useful to provide additional treatment (see operations 416 and 417, figure 8).

 It is observed here that the flowcharts presented with reference to FIGS. 5 and 7 combine particularly well, making the operation 262 of FIG. 7 succeed one another at the end of the operations 258 of FIG. 5.

 We also observe that, in the case illustrated in FIG. 5, the memorization of the value of the distortion D (operations 917 and 920, FIG. 6) associated with the memorized parameters, is not necessary, whereas in the case of the FIG. 7, it avoids having to carry out its determination again.

 FIGS. 11A and 11B illustrate the effect of the complementary partitions which are carried out on the elementary subsets originating from the initial partition or from previous complementary partitions.

 For explanatory purposes, and as an exception to the general rules of the description, it is assumed that, in the embodiment shown in FIGS. 11A to 12, the minimum dimensions of the elementary sub-assemblies are four times smaller than their maximum dimensions, that is to say 8 by 8 pixels.

FIG. 11A represents a part of an image, corresponding to twelve elementary subsets referenced Ko to K 1 resulting from the initial partition, in which:
the elementary subset K1 has undergone a complementary partition into four elementary subsets K1,0 to K13,
- the elementary subset K.1.3 has undergone a complementary partition into four elementary subsets K13 o to K1,3,3,
the elementary subset K8 has undergone a complementary partition into four elementary subsets K6.0 to K6.3,
the elementary subset K6 1 has undergone a complementary partition into four elementary subsets K6,1,0 to K6,13, and
- the elementary subset K11 has undergone a complementary partition into four elementary subsets K11,0 to K1t, 3.

 FIG. 11B illustrates, for the elementary subsets thus partitioned, the partition information which is associated with the elementary subsets resulting from the initial partition (subsets for which the size T is the initial size 32 × 32 pixels) or the elementary subsets resulting from complementary partitions.

 Thus, the elementary subsets K0, K2, K3, K4, K5, K7, K8, Kg and K10 are associated with zero partition information, the elementary subsets K1, K6 and K11 are associated with partition information being equal to 1, their elementary subsets K1.0, K1.1, K1.2, K6.0, K62, K63, K11.0, K11.1, K112 and K11.3 are associated with zero partition information, the elementary subsets K13 and K6.1, are associated with partition information equal to "1" and the elementary subsets resulting from their complementary partition are not associated with any partition information since the dimensions of these elementary subsets are the minimum dimensions.

The dimensions of all these elementary subsets can thus be unambiguously defined by the following partition information
01000100001010000001 0000
Figure I IC provides a representation of the parameters a1, a2, a3, a4 and the calculated value b, relating to a single image and which, together, are referenced here P and partition information associated with them, using the continued as set out above.

It is observed there that the sequence of the compressed data corresponding to the image part represented in FIG. 11A comprises:
- a header, representing in particular the dimensions of the image (worth 512 x 512 pixels here) the initial dimensions of the elementary subsets (corresponding here to 32 x 32 pixels), the ideal height and the definition of the minimum dimensions elementary subsets (corresponding here to 8 × 8 pixels) and additional processing information and the possible dimension ratio between the resolution of the luminance component and that of the chrominance components;
- the following partition information presented above: 01000100001010000001 0000 associated with any secondary data corresponding to them, during the storage carried out during steps 251 and 263 (without partition information), 258 and 261 (with partition information equal to "0") and 256 (with partition information equal to "1") knowing that the partition information equal to 1 is not associated with any secondary data.

 FIG. 1d illustrates a second variant, relating to the organization of the parameters, associated calculated values and partition information resulting from the primary processing of a single image.

We observe that the sequence of compressed data corresponding to the image part represented in FIG. 11A comprises:
- a header, representing in particular the dimensions of the image (worth 512 x 512 pixels here) the initial dimensions of the elementary subsets (corresponding here to 32 x 32 pixels), the ideal height and the definition of the minimum dimensions elementary subsets (corresponding here to 8 × 8 pixels) and additional processing information and the possible dimension ratio between the resolution of the luminance component and that of the chrominance components;
- the rest of the partition information presented above:
0 1 000 1 0000 1 01000000 1 0000, each of them being associated with the data P (parameters a1 ..., a4 and calculated value b) which correspond to it, during the storage carried out during steps 251 and 263 ( without partition information), 258 and 261 (with partition information equal to "0") and 256 (with partition information equal to "1"), knowing that, according to this second variant, even the partition information equal to 1 is associated with P data and are recorded with them in register (REG 1) 560.

 Thanks to this arrangement, as illustrated in FIG. IiE, by knowing the scale factor z (obtained during operation 420, FIG. 8, infra), it is possible to select the data P and the partition information to be used during secondary treatment.

In this case, for a scale factor equal to 1 / v, v being an integer, only the data are kept for secondary processing.
P corresponding to elementary subsets of dimensions greater than or equal to 2v x 2v pixels in the primary image, that is to say which will give subsets of dimensions greater than or equal to 2 x 2 pixels in the tertiary image.

 In the example of FIG. 11 E, when the minimum dimensions of the elementary subsets of the primary image, that is to say of the image which is processed during the primary processing, are 8 × 8 pixels and the scale factor is strictly less than 1/4, the elementary subsets of minimum dimension are not to be taken into account, because they would otherwise correspond to elementary subsets of strictly smaller dimensions at 2 x 2 pixels.

 In the example illustrated in FIG. 11 E, it is assumed that the scale factor is equal to one eighth, which amounts to dividing by sixty-four the number of points of the image resulting from the secondary processing. It is observed, in FIG. 11E, that none of the information concerning the elementary subsets of smaller dimensions, which resulted from two successive complementary partitions during the primary processing, is taken into account during the secondary processing.

 The implementation of this second variant (FIGS. 1 1 D and 11E) makes it possible to reduce the duration of the secondary processing, because it is not, here, necessary to process all the data P resulting from the primary image processing in order to reconstruct a smaller image.

Figure 12 is a representation of the data organization
P and partition information resulting from primary processing of several images. The P data and additional partition data concerning different images are organized as follows:
a file 900 has a common header 901 in which various general data common to all the images represented in this file 900 are recorded,
an index table 909, in which the addresses of the data sets relating to the images represented in said file 900 are represented. In the embodiment described and represented here, the addresses for the start of description of the images are represented in the global header in the form of translation data, that is to say data indicating how many successive addresses it is necessary to move in the file to reach the description of said image, this starting with a specific header.

- the specific headers 902, 904 and 906 represent, in particular:
the number of color components (here 3),
. the dimensions of the image (worth 512 x 512 pixels for a complete image or less for sub-images),
the ideal height of the strips,
the initial dimensions of the elementary subsets (corresponding here to 32 × 32 pixels),
the definition of the minimum dimensions of the elementary subsets (corresponding here to 8 × 8 pixels),
additional processing information, and
the possible ratio of dimensions between the resolution of the luminance component and that of the chrominance components;
It is noted here that preferably, in the case where some of this information is common to all the images in the file, then this information is contained in the common header 901 and not in the specific headers 902, 904 and 906.

 The information blocks 903, 905 and 907 comprise the sequences of partition information, each of them being associated with the data P which correspond to it, during the storage carried out during steps 251 (without partition information) and 258 (with partition information), knowing that the partition information worth 1 is not associated with any data (except, of course, if the second variant illustrated in FIGS. Il D and il E is implemented).

 Thanks to the arrangements for organizing files relating to a plurality of images illustrated in FIG. 12, the number of files used is limited as well as the operations of opening and closing of files, operations which have a particularly long duration. It is recalled here that a file is a set of data which jointly have the same address in memory. Secondary processing is thus accelerated.

 In particular, when the simultaneous display of several images of reduced size is provided to have a global view of these images, the use of a single file to represent several images considerably increases the speed of display of each reduced image.

 When the operation 207 (FIG. 4) is completed, the test 205 is repeated and all the operations 205 to 207 are repeated as long as the memory (REG 2) is not empty.

 When the result of test 205 is positive, operation 200E is repeated. When the result of test 201 is positive, test 200B is repeated.

When the result of the test 200B is positive, the operation 209A for determining additional treatment is carried out. During this operation 209A, it is carried out:
- on the one hand, when determining whether or not it is useful
additional smoothing treatment on the image resulting from a treatment
secondary P data and partition information, and
- on the other hand, when memorizing the description of this
image, namely its specific header (including in particular the
image dimensions, maximum and minimum dimensions of sub
elementary sets, the possible ratio of dimensions between the resolution
of the luminance component and that of the chrominance components
and the ideal height), its P data and partition information and
the complementary processing information representing the utility of said
additional treatment.

 To determine this utility, several criteria can be used. The first consists in comparing the minimum dimension of the elementary subsets with a dimension which serves as threshold. When the minimum dimension is greater than the threshold dimension, this means that relatively large elementary subsets (for example 8 × 8 pixels) will be used even if the complexity criterion, defined in step 253, is not not respected and that it will therefore be beneficial to reduce the edge feeling affecting such elementary sub-assemblies. The second, introduced above (operation 254, FIG. 10), consists in comparing the thresholds used in the table defining the complexity criterion with so-called "critical" values, of predetermined thresholds. When the thresholds actually used are higher than these critical values, this means that the complexity criterion is easily met on elementary large-sized subsets, a subset whose edges may be visible. It is then again useful to provide additional treatment.

 Next, operation 209B includes updating, in a common file (FIG. 12), a header common to several data sets resulting from the primary processing of several images (referenced 901, FIG. 12), updating date of an index table (referenced 909, FIG. 12) which contains addresses of data sets relating to images, in the file considered. Then the end step 210 is carried out.

 Thus, for at least two sets of primary data representative of physical quantities, the device according to the invention implements, for each set of primary data, steps 200A to 207, described with reference to FIGS. 4 to 7, then a step organization of the computer file 900 comprising secondary data representative of each set of primary data. This organizational step includes operations 209A and 209B.

 According to a variant not shown, each file resulting from the primary processing of an image is arranged in an individual specific file.

 Description of an embodiment of a decompression flowchart.

 We will describe, in support of Figures 8, 9 and 10,1'brogram of program for reconstructing the data of the original color image, or at least an approximate representation.

During operation 400, FIG. 8, it is carried out:
- the reading, in the specific header concerning the secondary data to be processed, of the total height of the image and the ideal height of the strips,
- an initialization of the registers (IM), (RE) and (IT) of the RAM 610, of the two memories (REG 1) 560 and (REG 2) 570, and of the intermediate variable n lines read, at the uniform value nothing.

 During operation 420, the scale factor z is read at the input, by means of a keyboard, by means of technical provisions not shown but well known to those skilled in the art.

During test 421, it is determined whether the variable ndlines ~ read is equal to the total height of the image, or not. When the result of test 421 is negative, operation 422 consists in initializing the variable CNT, which indicates the component number, at the value 1, in the register 513 (CNT). Then,
Operation 423 consists, in an identical manner to operation 200D (FIG. 4), of determining the strip height and of updating the variable ndlignes ~ read.

 The program then performs operation 403 during which the compressed data from the input memory for the current band are loaded into the memory (REG 1) 560. It is recalled that the compressed data are formed from a series of five or six numbers, namely the selected parameters a "..., a4 and the calculated value b and, optionally, of partition information, relating to each elementary subset Kj.

 Operation 403 consists first of all in opening the file 900, in which a plurality of images are represented (FIG. 12), then, in the common header 901 of the file 900, in reading the data relating to all the files and, in index table 909, the address of the data relating to the image which one wishes to process in order to load this data into the register 560.

 The program then performs operation 404, during which the list of references for each elementary subset is created. This is a list of coordinates (xj, yi) of the lower left pixel of each elementary subset Kj, as well as the size T of this subset, this size T being recalculated as a function partition information as exposed with the primary processing as well as the scale factor z which is associated with the entire image resulting from the secondary processing. This list is kept in the register (REG 2).

 Thus, during operation 404, when partition information representative of a complementary partition, that is to say having a predetermined value ("1" here), is assigned to an elementary subset whose the dimension is not minimal, the processor 650 performs a step of complementary partitioning of said elementary subset and allocates secondary data to each elementary subset originating from the complementary partition. The additional partition is performed as indicated above.

 It is noted here that the scale factor z is a multiplicative factor which applies to the dimensions of the image reconstructed during secondary processing from data P coming from an image processed during primary compression processing . This multiplicative factor also applies to the dimensions of each elementary subset, whether it comes from the initial partition or from complementary partitions, and to each dimension of working subset associated with one of these elementary subsets.

 When a ratio between the resolutions of the chrominance components and the luminance component has been transmitted in the header, the scale factor z is applied to the luminance and the product of the scale factor z by the said ratio is applied to the chrominance components, so that all the data sets on which a component is reconstructed have the same resolution.

 According to the invention, no scale factor z causes a pixel duplication forming the appearance, in the reconstructed image, of uniform blocks.

The program then connects to step 406 during which:
- in accordance with the embodiment of the invention described here, an initial arbitrary image having a number of lines equal to the strip height multiplied by the scale factor z, each line having a number of pixels equal to 512 multiplied by the scale factor z, is stored in the register (IM) 511 of the RAM 610. In the present case of decompression of a color image where each component has 256 levels having respectively the values 0 to 255, we initialize the register (IM) at the digital value 128 (corresponding substantially to the average value between 0 and 255) for each pixel for the first component processed. For the following components, the result obtained for the preceding component is used as initialization (Figure 9).

 Register (IT) 612 is initialized with a value p read from register 623 of RAM 610. The number of iterations is recorded in register 623 (p), taking into account the desired RMSE ("RMSE" is the notation classic of the square root of the mean square deviation, averaged over all pixels and all components, the efficiency of a compression algorithm being expressed by the compression rate obtained for a given RMSE). The inventors have indeed observed during tests with the prototype that, to obtain a given RMSE from a compression ratio, a limited number of successive iterations during the approximation calculation was necessary. Typically: 2 <p <10; however for very complex images p can have a higher value, -.

 The program then executes a number p of iterations of the loop referenced 450 (operations 407, 408, 409 and 410) in order to allow the restitution of a band of the component of the image.

 The program connects first to test 407 during which it is checked whether the content of register (IT) 612 is positive and not zero.

 If test 407 is positive, the program connects to test 408 during which it is checked whether the list recorded in (REG 2) 570 is empty or not. If the test 408 is negative, the program branches to step 410 which itself includes steps 411 to 415 (FIG. 10).

 In step 411, the coordinates (xj, yi) corresponding to the elementary subset Kj considered as well as its size are removed and they are copied into the registers 525 and 526 of the RAM RAM 510. At the same time, the selected parameters (jazz ..., a4) and the calculated value b, relating to the elementary subset considered are copied from the register (REG 1) to the register 520 of the part 614 (CACHE) of the RAM 610.

 The purpose of operation 412 is to determine the work subsets geometrically and predeterminedly linked to the elementary subset K1 considered, taking into account the fact that each work subset has double dimensions of the elementary subset considered (thus taking account of the scale factor z) and has a center offset, with respect to this elementary subset, by the half dimension of the elementary subset considered, in each direction.

 This operation 412 includes the same operations as the operation 250 described above. The values f (L), ..., f (L4) of the means, for each pixel of the elementary subset, of the four pixels which correspond to it in each of the working subsets L, ..., L4, are respectively copied into registers 516, 517, 518 and 519 of part 614 (CACHE) of the RAM 610.

à l'aide des paramètres ai,..., a4 et de la valeur b stockés dans le registre 520 et des données 1(L1), ..., gL4) stockées dans les registres 516 à 519, (étape 413) et on mémorise (étape 414) cette valeur dans le registre (RE) 621 de la
RAM 610.The program is then connected to steps 413 and 414 during which the value, for each pixel (x, y) of the elementary subset Ki, is calculated.

using the parameters ai, ..., a4 and the value b stored in the register 520 and the data 1 (L1), ..., gL4) stored in the registers 516 to 519, (step 413) and this value is stored (step 414) in the register (RE) 621 of the
RAM 610.

 If the test 408 is positive, the content of the register (RE) 621 is copied into the register (IM) 511, and the content of the register (IT) 612 is reduced by one.

 Test 407 is again carried out. As long as this test is positive, this means that the number of iterations provided for by the number recorded in the register 623 (p) is not finished and the operations 408 to 410 are repeated.

When the test 407 is negative or zero, the number of iterations planned is finished, for the component being processed, and the program leaves loop 450. This component is therefore decompressed.

 We then proceed to step 405 which consists of copying the content of the register (IM) 511 into the output buffer (OBFR) 690. Operation 424 then consists in incrementing the value of the variable CNT in the register by 1 (CNT) 513. Next, test 401 determines whether or not the content of counter 513 (CNT) is strictly greater than 3. As long as the result

 We observe that, according to general characteristics of the invention, the fixed point of the multilinear applications associated with all the elementary subsets (and with the working subsets, by the intermediary of the parameters a1) constitutes an approximation of all or part of this set, result of decompression, or secondary treatment, in so-called "tertiary" data.

 In FIG. 9, during the test 430, the counter 513 (CNT) is tested. If its content is equal to zero, this means that it is the first component, that of luminance, and the register (IM) 511 is initialized to 128 for each pixel (step 431). The register (IT) is initialized at the value p (step 432).

 If the content of the counter 513 (CNT) is greater than zero, this means that at least one component has been decompressed and, in step 433, the register (IM) 511 is loaded with the content of the register (RE) 621 , that is to say with the values corresponding to the component previously decompressed. This makes it possible to decrease the number of iterations p by two units (step 434 where we initialize (IT) to the value p-2) because the initial image is very close to the final image. The calculation time is reduced thanks to this characteristic. The counter 513 (CNT) is then incremented by one unit (step 435).

When the three components have been decompressed, the test 401 (FIG. 8) is positive, the decompression program is then terminated for the three components of the band of the image considered. The program then performs the test 421 described above. When the result of test 421 is positive, the decompression of the entire image is complete and test 416 then consists in determining whether
either the additional processing information present in the header of the data P and the partition information is representative of a utility for carrying out additional processing,
- either the scale factor z has a value which makes additional processing useful (z being, for example, greater than or equal to 4) or not.

 When the result of the test 416 is negative, the end step 402 is carried out directly. When the result of test 416 is positive, operation 417 of further processing is carried out. This operation 417 consists of an operation known to a person skilled in the art under the name of "post-processing" (which can be translated into French by "additional processing"), which refers to a processing carried out after decompression image illustrated in FIG. 8. This processing comprises, in known manner, spatial filtering, intended to reduce the possible border sensation of the elementary sub-assemblies.

 Preferably, the parameters of this filtering are a function of the dimensions of the subset processed. For example, the dimensions of the filter used are proportional to those of the subset considered, on the one hand, and the number of pixels to which this filter is applied is proportional to each dimension of the elementary subset considered, on the other hand go. According to each of these examples, the scale factor z intervenes in the configuration of the complementary processing. During operation 417, the data relating to each of the three components are, in turn, copied into the output buffer (OBFR) 690, at three different addresses each corresponding to one of the three components, replacing the stored data. during operation 405, from where they are transferred to the user means of decompressed data.

 For the same component, the successive transfer of the tertiary data resulting from the processing of the secondary data, during the successive operations 405 has the effect of assembling these bands in random access memory (610) of the device.

 According to a variant not shown, the independent bands are not formed from complete lines of the image but form a partition of the latter into rectangles, squares or any geometric shapes adapted to the method of the invention.

 It is observed here that, to obtain an image reconstructed in gray levels, without color, it suffices to carry out a primary processing only on the data coming from the primary processing of the luminance of the color image.

According to a variant not shown, the additional partitions are made by performing a non-homothetic division of the elementary subsets, such as for example a division into triangles, into two equal rectangles, successively by forming two rectangles side by side horizontally or vertically,
According to a variant not shown, the scale factor is not uniform. It can thus, for example, have a horizontal value, which serves as a multiplicative coefficient for the horizontal dimensions of the elementary subsets and a vertical value which serves as a multiplicative coefficient for the vertical dimensions of the elementary subsets. The scale factor can also have a value that varies depending on the distance from the center of the image.

 According to a preferred variant, the scale factor is different for the different components of the image processed. For example, the scale factor is double for each chrominance component than it is for the luminance component. In this case, a filtering is carried out preliminary to the primary processing (that is to say here to the compression), filtering which divides the resolution of the chrominance components by a factor of two in each direction, by giving each component of each pixel of the lower resolution image the average value of the corresponding component of four pixels of the initial image. This factor is transmitted in the header in the form of the aforementioned ratio.

 Conjointly, the number of initial data used for the reconstruction of the luminance component is, in each direction, double what it is for each of the chrominance components. Indeed, it is common that the complexity of the chrominance components is less than that of the luminance component, on the one hand and that this complexity is less perceived by the human eye.

 In combination with this variant, the primary processing device is adapted to perform, before the operations illustrated in FIG. 7, a filtering of the chrominance components in order to halve their resolution in each direction of the image.

 According to a variant not shown, each elementary subset is associated with two working subsets, two multiplicative parameters and an additive coefficient characterizing the multilinear application of the two working subsets in the elementary subset considered.

Preferably, the two working sub-assemblies cover the elementary sub-assembly.

 In this variant, and in accordance with considerations set out above concerning the restricted choice of values of the multiplicative parameters, the possible values of each multiplicative parameter are +1/2 and -1/2.

According to another variant not shown, the device for primary processing of primary data:
- performs, during operation 910 (FIG. 6), the initialization of the list of sixteen 4-uples of parameters given by table 1, supra,
- performs, in place of operation 914, an operation to remove the group of parameter values, here (1/4, 1/4, 1/4, 1/4), from the list of sixteen 4-uples possible and we repeat the successions of operations 911 to 917, with the new table thus made up of the last 15 4-uples of table 1.

* According to this variant, the second predetermined group of parameters which replaces the first group of parameters when this is equal to a first group of predetermined values is the one which, among the 15 4-uples tested, minimizes the value of the distortion . This variant has the advantage of better representing the primary data and the disadvantage of being slower than the method illustrated with regard to FIG. 6.

 According to a last variant not shown, the dimensions of the working sub-assemblies are equal to the dimensions of the elementary sub-assemblies.

 Of course, the present invention is in no way limited to the embodiments chosen and shown, but encompasses, on the contrary, all variants within the reach of the skilled person.

Claims (48)

 1. Method for primary processing of a set of so-called "primary" data representative of physical quantities to determine a set of so-called "secondary" data representative of the primary data, method characterized in that it comprises:  - a step of partitioning (204) of the primary data set into so-called "elementary" subsets (K1);  for each elementary subset, a step of determining (250, 259) at least two so-called "working" subsets (L,) of the primary data set, each associated in a predetermined manner, said elementary subset;  - for each elementary subset, an operation for selecting (915 to 917) a group of parameters (ai), from a predetermined set of groups of parameters, during which:  o one determines for each group of parameters of said  predetermined set, the value of a selection criterion (D)  of a parameterized combination of working sub-assemblies  associated with said elementary subset, each sub  work set being associated, by said combination, with  minus a said parameter of said group of parameters, and  . we select the group of parameters that optimizes the value  of said selection criterion, the parameter values which optimize the value of the selection criterion constituting secondary data representative of the primary data of the elementary subset considered.  2. A primary treatment method according to claim 1, characterized in that it comprises an operation of determining (913, 914) of at least one value (b) representative of a difference between said combination and the data of the sub- elementary set (K,), said representative value of difference constituting a secondary value representative of the primary data of the considered subset.  3. A primary treatment method according to any one of claims 1 or 2, characterized in that said combination is a multilinear application of the working subassemblies associated with an elementary subassembly therein.  4. A primary treatment method according to claim 3, characterized in that the parameters (a,) associated with the elementary subsets are coefficients of the multilinear application.  5. Primary treatment method according to any one of claims 1 to 4, characterized in that each datum of an elementary subset (K) is associated with at least one datum of each working subset (Ld associated with said subset - elementary unit, by said combination.  6. Primary processing method according to any one of claims 1 to 5, characterized in that the number of data of the working subsets (L,) associated with each elementary subset is greater than the number of data of said subset elementary (K,).  7. A primary treatment method according to any one of claims 1 to 6, characterized in that, during the selection operation (915 to 917), use is made, as a predetermined set of groups of parameters (a ,), groups of integers of a second predetermined value.  8. Primary treatment method according to claim 7, characterized in that, during the selection operation (915 to 917), groups of parameters are used as predetermined set of groups of parameters (ai) binary, that is to say likely to take only two values.  9. A primary treatment method according to claim 8, characterized in that the parameters are capable of taking as the value on the one hand the inverse value of the number of working subsets associated with each elementary subset and, on the other hand, the opposite of this reverse.  10. Primary treatment method according to any one of claims 1 to 9, characterized in that, during the selection operation (915 to 917), a selection criterion (D) is used whose value is a function of the mean square deviation of the differences between the primary data of the elementary subset (K,) and the data resulting from the combination of the working subsets (him).  11. A primary treatment method according to any one of claims 1 to 10, characterized in that, during the selection operation (915 to 917), the group of parameters which optimizes the value of the selection criterion (D ) is the one for which a distance between, on the one hand, the data values from the combination applied to the working subsets (him) and, on the other hand, the data values of the elementary subset (K1 ), is minimal.  12. A primary treatment method according to any one of claims 1 to 11, characterized in that:  - it implements a dimension of elementary subsets known as "minimum" and  - it comprises, for each elementary subset (K1) which does not have said minimum dimension, at least one iteration of a sequence of steps comprising:  o a division opportunity estimation stage (252 to 254,  262) during which an opportunity for division is estimated  of said elementary subset, and  . when the opportunity estimate does not meet criteria  predetermined so-called "of complexity", a stage of "partition  complementary "(256) of said subset into subsets  elementary.  13. A primary processing method according to claim 12, characterized in that the selected parameters associated with the elementary subsets (Ki) for which the opportunity estimate does not meet said predetermined complexity criteria constitute secondary data.  14. A primary processing method according to claim 12, characterized in that the selected parameters associated with the elementary subsets (K1) for which the opportunity estimate does not meet said predetermined complexity criteria do not constitute secondary data.  15. A primary treatment method according to any one of claims 12 to 14, characterized in that, during the step of estimating division opportunity (252 to 254, 262), said complexity criterion takes into account counts the dimension of said elementary subset (K1).  16. A primary processing method according to any one of claims 12 to 15, characterized in that at each additional partition step (256), so-called "partition" information is generated, said information constituting secondary data representative of the elementary subset of primary data (K1).  17. A primary treatment method according to any one of claims 1 to 16, characterized in that it comprises:  a step called "set division" (200D to 202) during which the primary data set is divided into primary data subsets, and  - The steps of partition (204), of determination of working subsets (250, 259), of determination of selection criterion (D), and of selection are applied successively to each of said primary subsets of data.  18. A method of primary processing of at least two sets of primary data representative of physical quantities, characterized in that:  - it implements, for each set of primary data, the steps of a primary processing method according to any one of claims 1 to 17, and  - It includes a stage of organization (209A, 209B) of a computer file (900) comprising secondary data representative of each set of primary data.  19. A method of secondary processing of a set of so-called "secondary" data representative of physical quantities to determine a set of so-called "tertiary" data, characterized in that it comprises:  a step of defining a so-called "initial" tertiary data set (406) having the same number of data as the tertiary data set to be determined,  a step of partitioning the initial data set (404) into so-called "elementary" subsets (K1) and allocating secondary data to each elementary subset,  for each elementary subset, a step of determining (412) at least two so-called "working" subsets (L) of the initial data set, each associated, in a predetermined manner, with said subset elementary set; and  - iteratively and, during each iteration, successively for each elementary subset:  an operation for determining the image of the subsets  of work associated with said elementary subset (413), by  a multilinear application parameterized by data  secondary assigned to said elementary subset, for  constitute a new elementary subset, the number  parameters of said multilinear application associated with said  elementary subset, being at least equal to the number of  working subsets associated with said subset  elementary.  20. A secondary processing method according to claim 19, characterized in that it implements secondary data capable (a,) of taking only two different values.  21. A method of secondary processing of a set of secondary data, resulting from a primary processing of a set of primary data, according to any one of claims 19 or 20, characterized in that the number of tertiary data is different the number of primary data.  22. A secondary treatment method according to any one of claims 19 to 21, characterized in that:  - it implements a minimum dimension,  - the secondary data may include partition information, and  - during the partition step (404), when partition information having a predetermined value is assigned to an elementary subset whose dimension is not minimal, a complementary partition step is carried out of said sub- elementary set and secondary data is assigned to each elementary subset from the complementary partition.  23. A method of secondary processing of secondary data sets, characterized in that a processing method according to any one of claims 19 to 22 is implemented on each set of secondary data, and in that it comprises an assembly step (405), during which the tertiary data from the secondary processing carried out on each of the secondary data sets are assembled.  24. Device for primary processing of a set of so-called "primary" data representative of physical quantities to determine a set of so-called "secondary" data representative of the primary data, device characterized in that it comprises:  - a primary data entry (591),  a means of partitioning (510, 530, 550) of the primary data set into so-called "elementary" subsets (Kj);  - suitable calculation means (510, 530, 550):  for each elementary subset, to be determined (250,  259) at least two so-called "working" subsets (Li) of  the primary data set, associated, each, so  predetermined, said elementary subset,  . for each elementary subset, to be selected (915 to  917) a group of parameters, from a set  predetermined groups of parameters, determining, for  each group of parameters of said predetermined set, the  value of a selection criteria for a parameterized combination  working subsets associated with said subset  elementary, each subset of work being associated,  by said combination, to at least one said parameter of said  group of parameters, and selecting the group of  parameters which optimize the value of said selection criterion, the parameter values which optimize the value of the selection criterion constituting secondary data representative of the primary data of the elementary subset considered.  25. Primary treatment device according to claim 24, characterized in that:  the calculation means (510, 530, 550) are adapted to determine at least one value (b) representative of a difference between said combination and the data of the elementary subset, and  - Said representative value of difference constitutes a secondary value on an output (581) of secondary data.  26. Primary treatment device according to any one of claims 24 or 25, characterized in that the calculation means (510, 530, 550) are adapted so that said combination is a multilinear application of the working sub-assemblies ( L) associated with an elementary subset (K1), in that.  27. Primary processing device according to claim 26, characterized in that the calculation means (510, 530, 550) are adapted so that the parameters (Q) associated with the elementary subsets are coefficients of the multilinear application .  28. Primary treatment device according to any one of claims 24 to 27, characterized in that the calculation means (510, 530, 550) are adapted so that each datum of an elementary subset (K,) is associated with at least one datum of each work subset (it) associated with said elementary subset, by said combination.  29. Primary processing device according to any one of claims 24 to 28, characterized in that the calculation means (510, 530, 550) are adapted so that the number of data of the working sub-assemblies (Li) associated with each elementary subset (Ki) is greater than the number of data of said elementary subset.  30. Primary processing device according to any one of claims 24 to 29, characterized in that the calculation means (510, 530, 550) are adapted to implement, as a predetermined set of groups of parameters, groups of integer multiples of a second predetermined value.  31. Primary processing device according to claim 30, characterized in that the calculation means (510, 530, 550) are adapted to implement, as a predetermined set of groups of parameters, groups of binary parameters, that is that is, likely to take only two values.  32. Primary treatment device according to claim 31, characterized in that the calculation means (510, 530, 550) are adapted to implement values of parameters which are, on the one hand, the value of the reverse the number of working subsets (Li) associated with each elementary subset (ksi) and, on the other hand, the opposite of this inverse.  33. Primary treatment device according to any one of claims 24 to 32, characterized in that the calculation means (510, 530, 550) are adapted to implement a selection criterion (D) whose value is a function of the mean square deviation of the differences between the primary data of the elementary subset (Ki) and the data resulting from the combination of the working subsets (Li).  34. Primary treatment device according to any one of claims 24 to 33, characterized in that the calculation means (510, 530, 550) are adapted so that the group of parameters which optimizes the value of the selection criterion ( D) is that for which a distance between the data values from the combination applied to the working subsets and the data values of the elementary subset is minimum.  35. Primary treatment device according to any one of claims 24 to 34, characterized in that the calculation means (510, 530, 550) are suitable:  - to implement a dimension of elementary subsets called "minimum" and  - for each elementary subset (Ki) which does not have said minimum dimension, to perform at least one iteration of a sequence of steps comprising:  o a division opportunity estimation stage (252 to 254,  262) during which they consider an opportunity to  division of said elementary subset, and  o when the opportunity estimate does not meet criteria  predetermined so-called "of complexity", a stage of "partition  complementary "(256) of said subset into subsets  elementary.  36. Primary processing device according to claim 35, characterized in that the calculation means (510, 530, 550) are adapted so that the selected parameters associated with the elementary subsets for which the opportunity estimate does not respond not to said predetermined complexity criteria constitute secondary data.  37. Primary processing device according to claim 35, characterized in that the calculation means (510, 530, 550) are adapted so that the selected parameters associated with the elementary subsets for which the opportunity estimate does not respond not to said predetermined complexity criteria do not constitute secondary data.  38. Primary processing device according to any one of claims 35 to 37, characterized in that the calculation means (510, 530, 550) are adapted so that said complexity criterion takes into account the dimension of said subset elementary.  39. Primary treatment device according to any one of claims 35 to 38, characterized in that the calculation means (510, 530, 550) are adapted so that at each complementary partition step (256), one generates so-called "partition" information, said information constituting secondary data representative of the set of primary data on the output (281) of secondary data.  40. Primary treatment device according to any one of claims 35 to 39, characterized in that the calculation means (510, 530, 550) are adapted to:  - divide the primary data set into primary data subsets, and  - determining the elementary (K1) and working (Li) subsets, and selecting a group of parameters (Q) successively for each of said primary data subsets.  41. Device for primary processing of at least two sets of primary data representative of physical quantities, characterized in that it comprises:  a primary processing device according to any one of claims 24 to 40, adapted to successively process said primary data sets and  a means of organization (510, 530, 550, 209A, 209B), of a computer file (900) comprising secondary data representative of each set of primary data.  42. Device for secondary processing of a set of so-called "secondary" data representative of physical quantities to determine a set of so-called "tertiary" data, characterized in that it comprises:  - a secondary data entry (681),  a means of defining a so-called "initial" tertiary data set (610, 630, 650, 400) having the same number of data as the tertiary data set to be determined,  a means of partitioning (610, 630, 650, 404) of the initial data set into so-called "elementary" subsets (Ki) and of allocating secondary data to each elementary subset,  - suitable calculation means (610, 630, 650):  for each elementary subset, to be determined (412) at  minus two so-called "working" subsets (Li) of the set  initial tertiary data, associated, each, in a way  predetermined to said elementary subset; and  o to perform, iteratively and, during each iteration,  successively for each elementary subset:  an operation for determining (413) the image of the sub  work packages associated with said subset  elementary, by a multilinear application parameterized by  secondary data allocated to said subset  elementary, to constitute a new subset  elementary, the number of parameters of said application  multilinear associated with said elementary subset,  being at least equal to the number of subsets of  work associated with said elementary subset, and  - an output (691) of tertiary data consisting of the data of the elementary subsets resulting from said iterations.  43: secondary processing device according to claim 42, characterized in that the calculation means (610, 630, 650) are adapted to implement secondary data likely to take only two different values.  44. Device for secondary processing of a set of secondary data, resulting from primary processing of a set of primary data, according to any one of claims 42 or 43, characterized in that the calculation means (610, 630, 650) are adapted so that the number of tertiary data on the tertiary data output (691) is different from the number of primary data.  45. A secondary treatment device according to any one of claims 42 to 44, characterized in that:  - it implements a minimum dimension,  - secondary data arriving at the secondary data entry (681) may include partition information, and  - the calculation means (610, 630, 650, 404) are suitable, when partition information having a predetermined value is assigned to an elementary subset (K1) whose dimension is not minimum, to perform a complementary partition of said elementary subset and to assign secondary data to each elementary subset from the complementary partition.  46. Device for secondary processing of secondary data sets, characterized in that it comprises:  a secondary processing device according to any one of claims 42 to 45, adapted to successively process each set of secondary data, and  - an assembly means (610, 630, 650, 405) adapted to assemble the tertiary data resulting from the secondary processing carried out on each of the secondary data sets, the tertiary data output (691) being adapted to transmit the assembled tertiary data by means of assembly.  47. Primary processing device according to any one of claims 24 to 41, characterized in that the calculation means (510, 530, 550) comprise a processor (550) and at least one memory (510, 530).  48. A secondary processing device according to any one of claims 42 to 46, characterized in that the calculation means (610, 630, 650) comprise a processor (650) and at least one memory (610, 630).