WO2008059752A1 - Procédé d'encodage et procédé de décodage de signal d'image, procédé d'encodage et procédé de décodage de source d'information, dispositifs pour ceux-ci, leurs programmes, et support de mémoire avec programme enregistré - Google Patents

Procédé d'encodage et procédé de décodage de signal d'image, procédé d'encodage et procédé de décodage de source d'information, dispositifs pour ceux-ci, leurs programmes, et support de mémoire avec programme enregistré Download PDF

Info

Publication number
WO2008059752A1
WO2008059752A1 PCT/JP2007/071719 JP2007071719W WO2008059752A1 WO 2008059752 A1 WO2008059752 A1 WO 2008059752A1 JP 2007071719 W JP2007071719 W JP 2007071719W WO 2008059752 A1 WO2008059752 A1 WO 2008059752A1
Authority
WO
WIPO (PCT)
Prior art keywords
information source
integer value
integer
dimensional
value
Prior art date
Application number
PCT/JP2007/071719
Other languages
English (en)
French (fr)
Inventor
Seishi Takamura
Original Assignee
Nippon Telegraph And Telephone Corporation
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nippon Telegraph And Telephone Corporation filed Critical Nippon Telegraph And Telephone Corporation
Priority to US12/446,818 priority Critical patent/US8300966B2/en
Priority to BRPI0718239A priority patent/BRPI0718239B1/pt
Priority to EP20070831450 priority patent/EP2083560B1/en
Priority to CA2667727A priority patent/CA2667727C/en
Priority to ES07831450.7T priority patent/ES2450265T3/es
Priority to CN2007800412601A priority patent/CN101536484B/zh
Priority to JP2008544115A priority patent/JP4769305B2/ja
Publication of WO2008059752A1 publication Critical patent/WO2008059752A1/ja

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N1/00Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
    • H04N1/41Bandwidth or redundancy reduction
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M7/00Conversion of a code where information is represented by a given sequence or number of digits to a code where the same, similar or subset of information is represented by a different sequence or number of digits
    • H03M7/30Compression; Expansion; Suppression of unnecessary data, e.g. redundancy reduction
    • H03M7/40Conversion to or from variable length codes, e.g. Shannon-Fano code, Huffman code, Morse code
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T9/00Image coding
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M7/00Conversion of a code where information is represented by a given sequence or number of digits to a code where the same, similar or subset of information is represented by a different sequence or number of digits
    • H03M7/30Compression; Expansion; Suppression of unnecessary data, e.g. redundancy reduction
    • H03M7/40Conversion to or from variable length codes, e.g. Shannon-Fano code, Huffman code, Morse code
    • H03M7/42Conversion to or from variable length codes, e.g. Shannon-Fano code, Huffman code, Morse code using table look-up for the coding or decoding process, e.g. using read-only memory
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N19/00Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
    • H04N19/10Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using adaptive coding
    • H04N19/102Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using adaptive coding characterised by the element, parameter or selection affected or controlled by the adaptive coding
    • H04N19/13Adaptive entropy coding, e.g. adaptive variable length coding [AVLC] or context adaptive binary arithmetic coding [CABAC]
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N19/00Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
    • H04N19/90Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using coding techniques not provided for in groups H04N19/10-H04N19/85, e.g. fractals
    • H04N19/91Entropy coding, e.g. variable length coding [VLC] or arithmetic coding
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N19/00Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
    • H04N19/90Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using coding techniques not provided for in groups H04N19/10-H04N19/85, e.g. fractals

Definitions

  • Image signal encoding method and decoding method information source encoding and decoding method, apparatus thereof, program thereof, and storage medium storing the program
  • the present invention relates to an image signal encoding method for simply and efficiently encoding an image signal indicating a Gaussian integer signal, and an image signal decoding method for decoding encoded data generated by the image signal encoding method
  • An information source encoding method and apparatus for encoding a Gaussian integer signal simply and efficiently, an information source decoding method and apparatus for decoding encoded data generated by the information source encoding method, and Information source encoding program used for realizing the information source encoding method, a computer-readable recording medium storing the program, an information source decoding program used for realizing the information source decoding method, and the program And a computer-readable recording medium on which is recorded.
  • a Gaussian signal is a signal whose occurrence probability follows a normal distribution (also called a Gaussian distribution), and this normal distribution is an extremely important distribution that appears in various mathematical / engineering situations.
  • the Gaussian signal taken up as an encoding target here is assumed to have an integer signal value. Assuming that generality is not lost, the average signal value is 0 and iid (independent and identically distributed).
  • the Golomb code is a signal that follows an exponential distribution (also called a Laplace distribution or a geometric distribution; unless otherwise noted, it refers to a bilateral exponential distribution). It is a code that can be encoded efficiently, does not require a table for encoding and decoding, is extremely simple to process, and is widely used. Furthermore, there is a feature that an arbitrarily large integer input value can be encoded.
  • the Golomb code is a Golomb code parameter (g) that takes an integer value of 1 or more. ) Changes the code representation.
  • the Golomb code has a relationship in which, for each code parameter g, the code length power increases as the value power 1 ⁇ 2 of the encoded integer z increases.
  • the Golomb code parameter g increases, the increase in the code length becomes gradual.For this reason, it is recommended to increase the Golomb code parameter g to encode a gentle distribution. Conversely, as the Golomb code parameter g decreases, the code length increases. Since the increase is steep, the Golomb code parameter g is suitable for encoding a steep distribution that concentrates on 0.
  • the Huffman code is a code (compact code) that can be instantaneously decoded and has the shortest average code amount among all the variable length codes using a code table
  • the arithmetic code is a variable length code. It is a code that can be compressed to the theoretical limit (sometimes exceeding the compact code) by a different method.
  • Patent Document 1 shown below describes a power S that describes an invention in which a feature quantity of an image is vector-quantized according to a binary search tree by a neural network.
  • This invention uses dynamic Huffman coding. This technique encodes an image.
  • FIG. 21 shows the frequency distribution of normal distribution and exponential distribution with mean 0 and variance 16.
  • the vertical axis is the frequency probability, and the logarithm is displayed so that both tails of the distribution are easy to understand.
  • the normal distribution is a parabola and the exponential distribution is a triangle.
  • each value must be converted to an integer of 0 or more.
  • a method of separately expressing positive / negative information and absolute value information, and the following relatively simple conversion are often used.
  • the encoding efficiency is significantly lower than that of the Huffman code or the arithmetic code. This is an essential problem because the signal occurrence probability assumed by the Golomb code follows an exponential distribution rather than a normal distribution.
  • both encoder and decoder require a frequency table that is not required for Golomb codes.
  • Non-Patent Document 2 shown below proposes a hybrid golomb code capable of instantaneous decoding of a generalized Gaussian signal source. This is capable of encoding / decoding arbitrarily large integer input values, but the structure is complicated, and it targets distributions that are steeper than exponential distributions! /, And! / was there. [0019]
  • Golomb encoding is performed at the expense of encoding efficiency. The method to do was taken. However, this method has a problem that it is necessary to optimize Golomb code parameters.
  • Patent Document 1 Japanese Patent Laid-Open No. 2001-5967
  • Non-patent 1 P Boldi, 3 ⁇ 4 Vigna: ompressed perfect embedded skip lists for quick ⁇ nvertea-index lookups, Proceedings of String Processing and Information Retrieval, 12th International Conference, pp.1-15, 2005
  • Non-Patent Document 2 S Xue, Y Xu, B Oelmann: "Hybrid Golomb codes for a group of quant ised GG sources, IEE Proceedings, Vision Image and signal Processing, Vol.150, No.4, pp.256-260 , 2003
  • the Golomb code is a code that can efficiently encode a signal according to an exponential distribution, does not require a table for encoding and decoding, and can be instantaneously decoded with extremely simple processing.
  • the feature is that an arbitrarily large integer input value can be encoded.
  • the Golomb code when used for coding a Gaussian signal, assumes the signal generation probability of exponential distribution, and the coding efficiency is higher than that of the Huffinan code or arithmetic code. It will drop greatly.
  • Non-Patent Document 2 proposes a hybrid golomb code.
  • this method has a problem that the configuration is complicated, and also has a problem that it targets a steeper distribution than a normal distribution.
  • Non-Patent Document 1 the Golomb code parameter is optimized and the Golomb code is used by giving priority to the simplicity of the encoding process without considering the encoding efficiency.
  • this method has the problem that the coding efficiency is lower than that of Huffman codes and arithmetic codes, and the Golomb code parameters need to be optimized.
  • the present invention has been made in view of power and circumstances, and is a new information source encoding / decoding technique that enables simple and efficient encoding / decoding of a Gaussian integer signal.
  • the purpose is to provide.
  • the information source encoding device of the present invention provides (i) input means for inputting a signal value sequence of a Gaussian integer signal to be encoded in order to realize simple and efficient encoding of a Gaussian integer signal. (Port) Integer value pair conversion means for converting the signal values included in the signal value series input by the input means into two integer value pairs in the input order; and (c) Conversion of the integer value pair conversion means.
  • a mapping means that obtains an integer value of 0 or more by applying a 2D ⁇ 1D mapping that maps each integer value pair as a grid point on 2D coordinates and maps the grid point closer to the origin to a smaller value
  • the encoding means may encode the integer value obtained by the mapping means using a Golomb code suitable for encoding the information source according to the exponential distribution.
  • a Golomb code suitable for encoding the information source according to the exponential distribution.
  • variance calculation means for calculating the variance of the signal value input by the input means
  • a value proportional to the variance value calculated by the variance calculation means Code parameter determining means for determining a code parameter of a Golomb code having.
  • the information source encoding apparatus of the present invention can also encode an image signal indicating a Gaussian integer signal. In this case, the information source encoding apparatus of the present invention It will function as an image signal encoding device.
  • the information source encoding method of the present invention realized by the operation of each processing means described above can also be realized by a computer program.
  • This computer program is recorded on an appropriate computer-readable recording medium.
  • the present invention is realized by operating on a control means such as a CPU, which is provided when the present invention is implemented and installed on the network.
  • the information source decoding apparatus provides (i) an integer value generated by the information source encoding apparatus of the present invention in order to realize decoding of the encoded data generated by the information source encoding apparatus of the present invention.
  • a one-dimensional ⁇ two-dimensional mapping that is a reverse mapping of the two-dimensional mapping
  • a restoring means that restores the integer value pair that is the mapping source of the integer value
  • a restored integer value pair is constructed.
  • an output means for outputting the integer values in that order.
  • the decoding means When adopting this configuration, when the information source encoding device of the present invention generates encoded data of integer values using Golomb codes, the decoding means performs Golo mb decoding of the encoded data of integer values. By decoding the integer value.
  • the information source encoding apparatus of the present invention determines the code parameter of the Golomb code having a value proportional to the variance of the signal value to be encoded, it is determined as such as the Golomb code parameter used for decoding.
  • a Golomb code parameter input means for inputting the Golomb code parameters thus obtained is provided.
  • the information source decoding device of the present invention when the information source encoding device of the present invention generates encoded data using an image signal indicating a Gaussian integer signal as an encoding target, the information source decoding device of the present invention performs image signal decoding. It will function as a device.
  • the information source decoding method of the present invention realized by the operation of each processing means described above can also be realized by a computer program.
  • the present invention is realized by being provided on a computer-readable recording medium or provided via a network, installed when the present invention is carried out, and operating on a control means such as a CPU. become.
  • the encoding efficiency of the variable-length code is improved by expanding the information source.
  • the integer-to-one integer value conversion process used in the present invention is nothing but secondary expansion of the information source. From now on, according to the present invention, it is possible to improve the encoding efficiency! RU
  • FIG. 1 is a diagram showing an example of a two-dimensional ⁇ one-dimensional mapping executed in the present invention.
  • FIG. 2 is an explanatory diagram of a table that stores the correspondence between integer value pairs and integer values.
  • FIG. 3 shows an apparatus configuration of an information source encoding device and an information source decoding device according to the present invention.
  • FIG. 4 is a flowchart executed by the information source encoding apparatus of the above embodiment.
  • FIG. 5 is an explanatory diagram of an algorithm for realizing a two-dimensional to one-dimensional mapping executed in the present invention.
  • FIG. 6 is an explanatory diagram of a result of an experiment performed to verify the effectiveness of the present invention.
  • FIG. 7] is also an explanatory diagram of the results of experiments conducted to verify the effectiveness of the present invention.
  • FIG. 8 is also an explanatory diagram of the results of experiments conducted to verify the effectiveness of the present invention.
  • FIG. 9 Similarly, it is an explanatory diagram of the results of experiments conducted to verify the effectiveness of the present invention.
  • FIG. 10 Similarly, it is an explanatory diagram of the results of experiments conducted to verify the effectiveness of the present invention.
  • FIG. 11 is a flowchart executed by the information source decoding apparatus of the above embodiment.
  • FIG. 12 is a detailed flowchart executed by the information source encoding apparatus of the embodiment.
  • FIG. 13 is a detailed flowchart executed by the information source encoding apparatus of the embodiment.
  • FIG. 14 is a detailed flowchart executed by the information source encoding apparatus of the embodiment.
  • FIG. 15 is a detailed flowchart executed by the information source encoding apparatus of the embodiment.
  • FIG. 16 is a detailed flowchart executed by the information source encoding apparatus according to the above-described embodiment.
  • FIG. 16 is a detailed flowchart executed by the information source decoding apparatus according to the above-described embodiment.
  • 18 This is a detailed flowchart executed by the information source decoding apparatus of the above embodiment.
  • 19 This is a detailed flowchart executed by the information source decoding apparatus of the above embodiment.
  • FIG. 20 is an explanatory diagram of Golomb codes.
  • Gaussian signals are input in the order of a, a, a, a,.
  • a two-element integer in the order of its input such as (a, a), (a, a), ...
  • an integer value z greater than or equal to 0 is obtained by applying a two-dimensional to one-dimensional mapping that maps the integer value pair (X, y) to a smaller value at a lattice point closer to the origin.
  • This two-dimensional to one-dimensional mapping is such that a grid point closer to the origin maps to a smaller value.
  • the integer value pair (X, y) is mapped to the integer value z in the form shown in Figure 1.
  • mapping process for example, as shown in FIG. 2, a table for storing the correspondence between integer value pairs and integer values is prepared in advance, and the integer value pair (X, y) is used as a key. By referring to the table, it is possible to obtain the integer value z that is the mapping result for the integer value pair (X, y).
  • the table prepared at this time is a general one that does not depend on the information source to be encoded. Can be prepared.
  • the minimum value of the distance between the unprocessed (not listed) grid point and the origin is obtained from the grid point of the origin, and the grid point of the minimum distance is defined. It is possible to obtain an integer value z that is the mapping result for an integer value pair (X, y) by repeating enumeration and assigning integer values according to the order of .
  • the z value of a lattice point that is separated from the origin by a distance L is approximately equal to the number of lattice points in the circumference of radius L, which is substantially equal to the area of the circumference. From now on, it becomes "z ⁇ ". Note that this is a global approximation (the actual number of grid points is quite large).
  • the probability f of this lattice point can be expressed as “f ex p ( —a L 2 ) (where a is a constant)” by a normal distribution. Therefore, from the relation of “z TT L 2 ”, “f cc exp (— az / ⁇ ) "Exponential distribution for z.
  • a two-dimensional to one-dimensional mapping is performed by using a Gaussian signal source as an exponential distribution signal source. Realizes a reversible mapping that converts to.
  • the Golomb code is a code that can efficiently encode a signal according to an exponential distribution, does not require a table for encoding and decoding, and can be instantaneously decoded with extremely simple processing.
  • the feature is that an arbitrarily large integer input value can be encoded.
  • the integer value ⁇ encoded by the information encoding method of the present invention is decoded in the input order of the encoded data.
  • the information source encoding method of the present invention when the information source encoding method of the present invention generates encoded data of an integer value using a Golomb code, the integer value is decoded by Golomb decoding the encoded data.
  • the decoded integer value z is subjected to a 1D ⁇ 2D mapping, which is a reverse mapping of the 2D ⁇ 1D mapping used by the information source coding method of the present invention,
  • the integer value pair (X, y) from which the integer value is mapped is restored.
  • This one-dimensional to two-dimensional mapping is performed by executing the mapping process of the two-dimensional to one-dimensional mapping used by the information source encoding method of the present invention, and when the decoded integer value z appears at that time. Then, the integer value pair (X, y) associated with it is specified.
  • this mapping process is similar to the mapping process performed by the information source encoding method of the present invention.
  • a table that stores beforehand the correspondence between integer value pairs and integer values as shown in FIG.
  • you can obtain an integer value pair (X, y) that is the mapping result for the integer value z. is there.
  • this mapping process starts from the grid point of the origin, obtains the minimum value of the distance between the grid point not listed and the original point, and lists the grid points of the minimum distance according to a prescribed order. Then, by assigning integer values repeatedly, it is possible to obtain an integer value pair (X, y) that is a mapping result for the integer value z. [0061] Next, the information source decoding method of the present invention outputs the integer values X and y constituting the restored integer value pair (X, y) in that order.
  • a Gaussian integer signal is obtained by decoding encoded data encoded with a code used for encoding an information source according to an exponential distribution. Decoding can be performed easily and efficiently.
  • the Gaussian signal source is converted into an exponential signal source through quadratic expansion and one-dimensional serialization of the Gaussian signal source, and this is encoded by a Golomb code having a property suitable for the exponential signal source. According to the embodiment, the present invention will be described in detail.
  • FIG. 3 shows an embodiment of the device configuration of the information source encoding device 1 and the information source decoding device 2 provided with the present invention.
  • the information source encoding device 1 of the present invention has a signal input unit 10 for inputting a signal value sequence of a Gaussian integer signal to be encoded, and an input of the signal input unit 10.
  • the integer value pair conversion unit 11 that converts signal values into two integer value pairs in the order of input, and the integer value pair converted by the integer value pair conversion unit 11, the smaller the lattice point closer to the origin. 2D ⁇ 1D mapping to obtain an integer value greater than or equal to 0 2D ⁇ 1D mapping unit 12 and Golomb code, 2D ⁇ 1D mapping unit 12 obtained integer value And a Golomb encoding unit 13.
  • the information source decoding device 2 of the present invention includes an encoded data input unit 20 for inputting encoded data of an integer value generated by the information source encoding device 1 of the present invention, and an encoded data input unit 20
  • the Golomb decoding unit 21 that decodes the integer value by performing Golomb decoding on the encoded data of the integer value input in step 2 and the mapping of the two-dimensional to one-dimensional mapping unit 12 included in the information source encoding device 1 of the present invention.
  • the Golomb decoding unit 21 When the integer value decoded by the Golomb decoding unit 21 appears at that time, by specifying the integer value pair associated therewith, the Golomb decoding unit 21 1 The integer values that make up the restored integer value pairs in the 1D ⁇ 2D inverse mapping unit 22 and the 1D ⁇ 2D inverse mapping unit 22 are output in that order. A sequential output unit 23 is provided.
  • FIG. 4 shows an example of a flow chart executed by the information source encoding apparatus 1 of the present invention configured as shown in FIG. Next, processing executed by the information source encoding device 1 of the present invention will be described in detail according to this flowchart.
  • step S101 the information source encoding apparatus 1 of the present invention inputs a signal value of a Gaussian integer signal two elements at a time from the top.
  • step S102 the Gaussian properties are inputted in the order of a, a ⁇ a, a ⁇ ,.
  • a numerical signal is converted into an integer value pair of two elements in the order of input, such as (a, a), (a, a), ...
  • step S 103 a two-dimensional to one-dimensional mapping is performed on this integer value pair (X, y), which is mapped to a smaller value at a lattice point closer to the origin. get z.
  • step S 104 the integer value z is Golomb-coded using a separately provided Golomb code parameter g, and in the subsequent step S 105, the encoded data is output.
  • step S 106 it is determined whether or not the input of the Gaussian integer signal has been exhausted. If it has not been exhausted, the process returns to step S 101 again, and if it has been exhausted, the process ends.
  • step S 103 Next, the 2D ⁇ 1D mapping executed in step S 103 will be described.
  • l (x, y) shown in Fig. 5 gives the distance between the origin and the point (x, y) and is usually the Euclidean distance.
  • the distance function l (x, y) is not necessarily isotropic, and the input Gaussian source is not strictly iid (ind. Endent and identically distributed). If there is, use an appropriate constant ⁇ ,
  • L l norm or more generally the L y norm ( ⁇ of each signal element) shown in the following equation.
  • FIG. 6 shows a frequency distribution of a signal source that is very close to a Gaussian signal obtained by processing an actual image.
  • the vertical axis is the logarithm of probability. It can be seen that the frequency distribution of this signal source is parabolic, that is, follows a normal distribution.
  • Fig. 7 shows the signals of this signal source sequentially taken out as two numerical pairs (X, y) and their frequency distributions displayed in three dimensions. Similarly, the vertical axis is the logarithm of probability. The graph is a paraboloid that rotates the parabola.
  • FIG. 8 shows the contour lines of the paraboloid of revolution. Although a slight positive correlation is observed, it is almost concentric.
  • the exponential distribution shows a triangular force signal as shown in Fig. 21.
  • the one-sided distribution by this conversion becomes a straight line that descends to the right.
  • FIG. 9 shows the frequency distribution of the two-dimensional to one-dimensional mapping applied to (X, y) according to the present invention and converted to the one-dimensional value z as a vertical logarithmic graph similar to FIG.
  • FIG. 10 shows a comparison of the obtained encoding efficiencies.
  • datal corresponds to the above data.
  • the coding efficiency of the present invention has achieved at least about 90%.
  • the present invention has an average of 94.1% and the conventional method has a difference of about 7 points, 87.4%.
  • FIG. 11 shows an example of a flow chart executed by the information source decoding apparatus 2 of the present invention configured as shown in FIG.
  • the information source decoding device 2 of the present invention firstly, at step S201, the encoded data (Golomb code of integer value z) generated by the information source encoding device 1 of the present invention from the head. Enter one element at a time according to the order.
  • step S202 the input Golomb code is decoded using the Golomb code parameter g to obtain an integer value z of 0 or greater.
  • step S203 the integer value z is converted to an integer value pair (X
  • this integer value pair (X, y) is output in the order of X, y.
  • the one-dimensional ⁇ two-dimensional mapping executed in step S203 executes the mapping process of the two-dimensional ⁇ one-dimensional mapping used by the information source coding apparatus 1 of the present invention, and is then decoded.
  • an integer value z appears, specify the integer value pair (X, y) associated with it.
  • step S205 it is determined whether or not the input of the Golomb code has been exhausted. If V, Na! /, Return to Step S 201, and if it is exhausted! /, End.
  • FIG. 12 shows a detailed flowchart of the process executed in step S102 of the flowchart of FIG.
  • step S102 when the signal values of the two Gaussian integer signals input in step S101 are received, as shown in the flowchart of Fig. 12, first, in step S301, of them. The first signal value X is input and stored in the memory X. Subsequently, in step S302, the other signal value y is input and stored in the memory y.
  • step S303 the signal X is read from the memory X, the signal value y is read from the memory y, and these X and y are output as an integer value pair (X, y).
  • step S102 by executing the flowchart of FIG. 12, the process of bundling and outputting two signal values in the signal sequence of the Gaussian integer signal is performed. .
  • FIGS. 13 to 16 show detailed flowcharts of the processing executed in step S103 of the flowchart of FIG.
  • step S103 the two-dimensional to one-dimensional mapping realized by the algorithm shown in Fig. 5 is performed to map the integer value pair (X, y) to the integer value z. Performs processing to output the numerical value z.
  • step S401 when the integer values xO and yO to be mapped are input, the lattice point storage memory Z is emptied in the subsequent step S402. At the same time, the variable z is initialized to 0.
  • step S403 procedure A defined by the processing of the flowchart shown in FIG. Run and continue
  • step S404 procedure B defined by the processing of the flowchart shown in FIG. 15 is executed.
  • procedure A is set by procedure B "dmin ⁇ x ⁇ dmin,- For dmin ⁇ y ⁇ dmin ", the minimum value dmin of the distance to the origin is determined for the unenumerated grid point (x, y) to which no integer value is assigned. At this time, there are usually a plurality of grid points that realize this minimum value dmin.
  • procedure B the grid points are selected one by one in the prescribed order, and an integer value one larger than the previous one is assigned to them! /, And so on.
  • X and y are both dmin and less than dmin (dmin is an integer, and is incremented by 1 by procedure B).
  • the minimum value of the distance from the origin is obtained as a processing target, and it is set as dmin.
  • step S501 the integer value dmin is initialized to zero. Subsequently, in step S502, one set is generated for each loop (S502 to S508) in accordance with the specified order without any interval from X to dmin and y from dmin to dmin.
  • step S503 the generated (X, y) is already stored in the lattice point storage memory Z. If it is already stored, the process proceeds to step S508. If it is not stored, the process proceeds to step S504, and the variable d is processed by step S505. Substitute the distance l (x, y) between the required grid point (X, y) and the origin.
  • the distance Kx, y) used here is, for example,
  • step S506 the magnitudes of d and dmin are compared. If d is greater than or equal to dmin, the process immediately proceeds to step S508. If d is less than dmin, the process proceeds to step S507. Proceed to set the value of dmin to d, and then proceed to step S508.
  • step S508 it is determined whether all combinations of X and y are within the range from dmin to dmin. If not, the process of step S502 is performed. Return to, and if it is exhausted, the processing as the procedure A ends.
  • step S601 the flag found is set to 0. Subsequently, in step S602, X is generated from dmin to dmin and y from dmin to dmin, one set for each one crepe (S602 to S609) in the prescribed order.
  • step S603 it is determined whether or not the generated (x, y) has already been stored in the lattice point storage memory Z. Proceeding to S609, if not stored, proceeding to step S604, and for the variable d, the distance l (x, y) between the lattice point (X, y) obtained by the processing of step S605 and the origin. Substitute y).
  • the distance Kx, y) used here is, for example,
  • step S606 d and dmin are compared. If d and dmin do not match! /, Proceed immediately to step S609, and if d and dmin match, proceed to step S607 and execute procedure X.
  • step S608 a new integer value pair (x, y) is stored in the lattice point storage memory Z, the value of the separately defined integer z is incremented by 1, and the flag Set the found value to 1 before proceeding to step S609.
  • step S609 it is determined whether all combinations of X and y are within the range from dmin to dmin. If not, the process of step S602 is performed. If it returns and runs out, the process proceeds to step S610 to determine whether the value of the flag found is equal to 0! /.
  • step S610 When it is determined that the value of flag found is not equal to 0 according to the determination process of step S610, the process returns to step S601, and when it is determined that the value of flag found is equal to 0, step Proceeding to S611, the value of dmin is incremented by 1, and the processing as procedure B ends.
  • step S103 of the flowchart of Fig. 4 the two-dimensional to one-dimensional mapping realized by the algorithm shown in Fig. 5 is performed. By doing so, the process of mapping the integer value pair (X, y) to the integer value z is performed, and the process of outputting the integer value z is performed.
  • step S203 the integer value z decoded in the process of step S202 is converted into an integer value pair (X, y) by applying a one-dimensional to two-dimensional mapping realized by the algorithm shown in FIG. The mapping process is performed, and the integer value pair (X, y) is output.
  • step S801 when the integer value ⁇ to be mapped is input, in the subsequent step S802, the lattice point storage memory Z is emptied, Initialize variable z to 0.
  • step S803 procedure A defined by the process of the flowchart shown in FIG. 14 is executed, and in step S804, procedure B defined by the process of the flowchart shown in FIG. 15 is executed ( However, “Procedure X” is changed to “Procedure X '” described below)
  • step S803 if the condition is not satisfied, nothing is output and the process returns to step A by returning to step S803.
  • the numerical value (X, y) is output as the mapping result, the process proceeds to step S903, and the one-dimensional ⁇ two-dimensional mapping process is terminated.
  • the procedure X ′ is terminated (that is, the procedure returns to step S608 of procedure B).
  • step S203 of the flowchart of Fig. 11 the one-dimensional to two-dimensional copy realized by the algorithm shown in Fig. 5 is performed.
  • the process of mapping the integer value z to the integer value pair (x, y) is performed, and the process of outputting the integer value pair (X, y) is performed.
  • FIG. 19 shows a detailed flowchart of the process executed in step S204 of the flowchart of FIG.
  • step S204 when the integer value pair (X, y) obtained in the process of step S203 is received, as shown in the flowchart of FIG. 19, first, in step S 1001, the top of them is displayed. Is stored in the memory X, and the other signal value y is stored in the memory y.
  • step S1002 the first signal value X is read out from the memory X and output, and in the subsequent step S1003, the other signal value y is read out from the memory y and output.
  • step S204 by executing the flowchart of FIG. 19, the integer value X, y constituting the integer value pair (X, y) of the mapping result is output according to the order. It does.
  • the present invention encodes and decodes a Gaussian integer signal, and uses a configuration that realizes a reversible mapping that converts a Gaussian signal source to a signal source having an exponential distribution. Although it appears in various scenes, it is possible to simply and efficiently encode and decode Gaussian integer signals that could not be encoded efficiently by conventional Golomb codes.

Landscapes

  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Compression Or Coding Systems Of Tv Signals (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)
  • Compression Of Band Width Or Redundancy In Fax (AREA)
  • Image Processing (AREA)

Description

明 細 書
画像信号符号化方法及び復号方法、情報源符号化及び復号方法、それ らの装置、及びそれらのプログラム並びにプログラムを記録した記憶媒体
技術分野
[0001] 本発明は、ガウス性整数信号を示す画像信号を簡易かつ効率よく符号化する画像 信号符号化方法と、その画像信号符号化方法により生成された符号化データを復号 する画像信号復号方法と、ガウス性整数信号を簡易かつ効率よく符号化する情報源 符号化方法およびその装置と、その情報源符号化方法により生成された符号化デー タを復号する情報源復号方法およびその装置と、その情報源符号化方法の実現に 用いられる情報源符号化プログラムおよびそのプログラムを記録したコンピュータ読 み取り可能な記録媒体と、その情報源復号方法の実現に用いられる情報源復号プロ グラムおよびそのプログラムを記録したコンピュータ読み取り可能な記録媒体とに関 する。
本願 (ま、 2006年 11月 14曰 ίこ出願された特願 2006— 307512号 ίこ基づき優先権 を主張し、その内容をここに援用する。
背景技術
[0002] ガウス性信号は、生起確率が正規分布(ガウス分布とも呼ばれる)に従う信号であり 、この正規分布は数学 ·工学の様々な場面で現れる極めて重要な分布である。
[0003] ここで符号化対象として取り上げるガウス性信号は、信号値が整数のものであると する。また、一般性を失わない仮定として、信号平均値は 0で、 iid(independent and id entically distributed:独立かつ同一分布)であるとする。
[0004] 整数信号の符号化方式は多数知られているが、それらのうち Golomb符号は、指数 分布 (ラプラス分布や幾何分布とも呼ばれる。特に断らない限り両側指数分布のこと を指す)に従う信号を効率よく符号化でき、かつ符号化'復号のための表が不要で、 処理が極めて簡易な瞬時復号可能な符号であり、広く用いられている。さらに、任意 に大きな整数入力値を符号化することが可能という特長がある。
[0005] Golomb符号は、 1以上の整数値をとる Golomb符号パラメータ(gと呼ぶことにする )により、その符号表現が変化する。
[0006] 図 20に示す表は、整数 z = 0,...,10に対する Golomb符号(Golomb符号パラメータ
§ = 1,...,6)でぁる。 Golomb符号には、各符号パラメータ gにおいて、被符号化整数 z の値力 ½の値だけ増すと符号長力 増す、という関係がある。
Golomb符号パラメータ gが増すにつれ符号長の増加は緩やかになるため、緩やか な分布の符号化には Golomb符号パラメータ gを大きくすることが向き、逆に、 Golomb 符号パラメータ gが減少するにつれて符号長の増加は急になるため、 0に集中する 急峻な分布の符号化には Golomb符号パラメータ gを小さくすることが向く。
[0007] ここで、指数分布に従う信号としては、画像信号中、時空間的に隣接する画素間の 輝度差や、画素輝度値の直交変換係数などがあり、この符号化に Golomb符号が用 いられること力 Sfcる。
[0008] 一方、 Huffman符号は、符号表を用い、瞬時復号可能、かつ平均符号量が全ての 可変長符号中最短となる符号 (コンパクト符号)であり、また、算術符号は、可変長符 号と異なる方式により定常的信号源を (コンパクト符号をときとして上回る)理論限界ま で圧縮可能な符号である。
[0009] 下記に示す特許文献 1には、ニューラルネットワークによって画像の特徴量を 2進 探索木に従ってベクトル量子化するという発明が記載されている力 S、この発明は動的 ハフマン符号化を利用して画像を符号化するという技術である。
[0010] 図 21に、平均 0、分散 16の正規分布および指数分布の頻度分布を示す。縦軸は 頻度確率であり、分布の両裾がわ力、りやすいようにと、対数で表示している。このよう に、縦対数グラフでは正規分布は放物線、指数分布は三角形となる。
[001 1] このような左右対称な整数分布は、符号化するために、各値を 0以上の整数に変換 する必要がある。これには、正負情報と絶対値情報とに分離して表現する方法や、次 のような比較的単純な変換などがよく用いられる。
[0012] この変換は、変換前の値を a、変換後を bとすると
b = 2a - l a〉0の場合
b = - 2a a≤0の場合
とするものであり、この変換により、 aとして例えば 3 , - 2 , - 1 , 0 , 1 , 2 , 3は、それ ぞれ bとして 6, 4, 2, 0, 1 , 3, 5に変換される。
[0013] Golomb符号をガウス性信号の符号化に用いた場合、符号化効率が Huffman符号 や算術符号のそれより大きく低下してしまう。これは、 Golomb符号が仮定している信 号生起確率が、正規分布ではなく指数分布に従うことに起因する、本質的な問題で ある。
[0014] Golomb符号のように符号表を必要としない符号には Fibonacci符号、 Elias符号、 E xp- Golomb符号など様々な種類が存在するが、信号の生起確率が正規分布に従 うことを仮定した符号は存在しなレ、。
[0015] 下記に示す非特許文献 1力 その 8頁で、 "Contrarily to what happens for geometri c and double-sided geometric distributions, there is no simple, instantaneous code f or the normal distribution.",すなわち、「指数分布や両側指数分布の場合と異なり、 正規分布のための簡易かつ瞬時復号可能な符号はない」と述べているように、信号 の生起確率が正規分布に従うことを仮定した符号は存在しない。
[0016] そのため、ガウス性信号の符号化において、「符号化効率」を優先する場合には、 従来、 Huffman符号あるいは算術符号などが用いられて!/、た。
[0017] しかしながら、これには、
•算術符号では、符号化器 '復号器ともに、 Golomb符号では必要としない頻度表 が必要となる
•Huffman符号では、符号化器'復号器ともに、 Golomb符号では必要としない符 号表あるいは頻度表が必要となる
•両者ともに、例外処理なしには任意入力値の符号化ができない、すなわち入力 値の範囲を予め知りうる必要がある
•両者ともに、処理量が Golomb符号より多ぐ特に算術符号は多い
などの問題があった。
[0018] また、下記に示す非特許文献 2では、一般化ガウス性信号源の瞬時復号可能な Hy brid Golomb符号を提案している。これは任意に大きい整数入力値の符号化復号が 可能であるが、構成が複雑であり、正規分布より急峻な指数分布よりさらに急峻な分 布を対象として!/、ると!/、う問題があった。 [0019] また、ガウス性信号の符号化において、「符号化処理の簡易さ」を優先する場合に は、従来、非特許文献 1に記載されるように、符号化効率を犠牲にし Golomb符号化 する方法が採られていた。ただし、この方法では、 Golomb符号パラメータを最適化 する必要があるという問題があった。
特許文献 1 :特開 2001— 5967号公報
非特許文 1: P Boldi, ¾ Vigna: し ompressed perfect embedded skip lists for quick ι nvertea-index lookups, Proceedings of String Processing and Information Retrieval, 12th International Conference, pp.1-15, 2005
非特許文献 2 : S Xue, Y Xu, B Oelmann: "Hybrid Golomb codes for a group of quant ised GG sources, IEE Proceedings, Vision Image and signal Processing, Vol.150, N o.4, pp.256- 260, 2003
発明の開示
発明が解決しょうとする課題
[0020] 正規分布に従う一般的な画像や音声などの種々の信号源を効率的に符号化する ユニバーサル符号が広く求められているものの、非特許文献 1に記載されているよう に、信号の生起確率が正規分布に従うことを仮定した符号は存在しない。
[0021] そこで、 Golomb符号を使って、ガウス性信号を符号化するという方法を用いること が考えられる。
[0022] 前述したように、 Golomb符号は、指数分布に従う信号を効率よく符号化でき、かつ 符号化'復号のための表が不要で、処理が極めて簡易な瞬時復号可能な符号であり 、さらに、任意に大きな整数入力値を符号化することが可能という特長がある。
[0023] しかしながら、 Golomb符号は、指数分布の信号生起確率を仮定してレ、ること力 、 ガウス性信号の符号化に用いた場合には、符号化効率が Huffinan符号や算術符号 のそれより大きく低下してしまう。
[0024] これから、従来では、ガウス性信号を符号化する場合には、 Huffman符号や算術符 号を用いるようにしている。
[0025] しかしながら、 Huffman符号や算術符号を用いる場合には、 Golomb符号では必要 としない頻度表が必要になるという問題や、入力値の範囲を予め知りうる必要がある という問題や、処理量が Golomb符号よりも多いという問題がある。
[0026] このような問題を考慮して、非特許文献 2では、 Hybrid Golomb符号を提案してレ、 る。し力、しながら、この方法は構成が複雑であるという問題があるとともに、正規分布よ り急峻な指数分布よりさらに急峻な分布を対象としているという問題がある。
[0027] また、非特許文献 1では、符号化効率のことを考えずに、符号化処理の簡易さを優 先して、 Golomb符号パラメータを最適化して Golomb符号を用いるようにしている。 しかしながら、この方法は符号化効率が Huffman符号や算術符号よりも低下するとい う問題があるとともに、 Golomb符号パラメータを最適化する必要があるという問題が ある。
[0028] 本発明は力、かる事情に鑑みてなされたものであって、ガウス性整数信号を簡易か つ効率よく符号化 ·復号することができるようにする新たな情報源符号化 ·復号技術 の提供を目的とする。
課題を解決するための手段
[0029] 〔 1〕本発明の情報源符号化装置
本発明の情報源符号化装置は、ガウス性整数信号を簡易かつ効率よく符号化する ことを実現するために、(ィ)符号化対象のガウス性整数信号の信号値系列を入力す る入力手段と、(口)入力手段の入力した信号値系列に含まれる信号値を、その入力 順に二個ずつの整数値対とする整数値対変換手段と、(ハ)整数値対変換手段の変 換した整数値対の各々を 2次元座標上の格子点と見て、原点に近い格子点ほど小さ な値に写像する 2次元→1次元写像を施すことで 0以上の整数値を得る写像手段と、 (二)指数分布に従う情報源の符号化に用いられる符号を使って、写像手段の得た 整数値を符号化する符号化手段とを備えるように構成する。
[0030] この構成を採るときに、符号化手段は、指数分布に従う情報源の符号化に好適な G olomb符号を使って写像手段の得た整数値を符号化することがあり、このときには、 Golomb符号パラメータの自動設定を可能とするために、(i)入力手段の入力した信 号値の分散を算出する分散算出手段と、(ii)分散算出手段の算出した分散値に比 例する値を持つ Golomb符号の符号パラメータを決定する符号パラメータ決定手段と を備えることがある。 [0031] ここで、本発明の情報源符号化装置は、ガウス性整数信号を示す画像信号を符号 化対象とすることも可能であり、この場合には、本発明の情報源符号化装置は画像 信号符号化装置として機能することになる。
[0032] 以上の各処理手段が動作することで実現される本発明の情報源符号化方法はコン ピュータプログラムでも実現できるものであり、このコンピュータプログラムは、適当な コンピュータ読み取り可能な記録媒体に記録して提供されたり、ネットワークを介して 提供され、本発明を実施する際にインストールされて CPUなどの制御手段上で動作 することにより本発明を実現することになる。
[0033] 〔2〕本発明の情報源復号装置
本発明の情報源復号装置は、本発明の情報源符号化装置の生成した符号化デー タを復号することを実現するために、(ィ)本発明の情報源符号化装置の生成した整 数値の符号化データを復号することで、その整数値を復号する復号手段と、(口)復 号手段の復号した整数値に対して、本発明の情報源符号化装置が用いた 2次元→1 次元写像の逆写像である 1次元→2次元写像を施すことで、その整数値の写像元と なった整数値対を復元する復元手段と、(ハ)復元手段の復元した整数値対を構成 する整数値を、その順番に出力する出力手段とを備えるように構成する。
[0034] この構成を採るときに、本発明の情報源符号化装置が Golomb符号を使って整数 値の符号化データを生成するときには、復号手段は、整数値の符号化データを Golo mb復号することで整数値を復号する。
このとき、本発明の情報源符号化装置が、符号化対象の信号値の分散に比例する 値を持つ Golomb符号の符号パラメータを決定するときには、復号に用いる Golomb 符号パラメータとして、そのようにして決定された Golomb符号パラメータを入力する Golomb符号パラメータ入力手段を備えることになる。
[0035] ここで、本発明の情報源符号化装置がガウス性整数信号を示す画像信号を符号化 対象として符号化データを生成する場合には、本発明の情報源復号装置は画像信 号復号装置として機能することになる。
[0036] 以上の各処理手段が動作することで実現される本発明の情報源復号方法はコンビ ユータプログラムでも実現できるものであり、このコンピュータプログラムは、適当なコ ンピュータ読み取り可能な記録媒体に記録して提供されたり、ネットワークを介して提 供され、本発明を実施する際にインストールされて CPUなどの制御手段上で動作す ることにより本発明を実現することになる。
発明の効果
[0037] 以上に説明したように、本発明によれば、数学 ·工学の様々な場面で現れるにもか かわらず、従来の Golomb符号などでは効率的に符号化ができなかったガウス性整 数信号を、簡易かつ効率よく符号化し復号することができるようになる。
[0038] また、一般に情報源の拡大により可変長符号の符号化効率は改善する。本発明で 用いる整数値対一整数値の変換処理は情報源の二次拡大に他ならず、これから、 本発明によれば、符号化効率を改善できると!/、う利点を持って!/、る。
図面の簡単な説明
[0039] 図 1]本発明で実行する 2次元→1次元写像の一例を示す図である。
図 2]整数値対と整数値との対応関係について記憶するテーブルの説明図である。 図 3]本発明の情報源符号化装置および情報源復号装置の装置構成の一実施形態 ί列である。
図 4]上記実施形態例の情報源符号化装置の実行するフローチャートである。
図 5]本発明で実行する 2次元→1次元写像を実現するアルゴリズムの説明図である 図 6]本発明の有効性を検証するために行った実験結果の説明図である。
図 7]同様に、本発明の有効性を検証するために行った実験結果の説明図である。 図 8]同様に、本発明の有効性を検証するために行った実験結果の説明図である。 図 9]同様に、本発明の有効性を検証するために行った実験結果の説明図である。 図 10]同様に、本発明の有効性を検証するために行った実験結果の説明図である。 図 11]上記実施形態例の情報源復号装置の実行するフローチャートである。
図 12]上記実施形態例の情報源符号化装置の実行する詳細なフローチャートである
[図 13]上記実施形態例の情報源符号化装置の実行する詳細なフローチャートである [図 14]上記実施形態例の情報源符号化装置の実行する詳細なフローチャートである
[図 15]上記実施形態例の情報源符号化装置の実行する詳細なフローチャートである
[図 16]上記実施形態例の情報源符号化装置の実行する詳細なフローチャートである 園 17]上記実施形態例の情報源復号装置の実行する詳細なフローチャートである。 園 18]上記実施形態例の情報源復号装置の実行する詳細なフローチャートである。 園 19]上記実施形態例の情報源復号装置の実行する詳細なフローチャートである。
[図 20]Golomb符号の説明図である。
園 21]正規分布および指数分布の説明図である。
符号の説明
1 情報源符号化装置
2 情報源復号装置
10 信号入力部
11 整数値対変換部
12 2次元→1次元写像部
13 Golomb符号化部
20 符号化データ入力部
21 Golomb復号部
22 1次元→2次元逆写像部
23 順次出力部
発明を実施するための最良の形態
[0041] 本発明の情報源符号化方法では、 a ,a ,a ,a,····というような順に入力されるガウス
1 2 3 4
性整数信号を、(a , a ), (a , a ), ·····というように、その入力順に 2要素ずつの整数
1 2 3 4
値対に変換する。この整数値対を (X, y)とする。
[0042] 続いて、この整数値対 (X, y)に対して、原点に近い格子点ほど小さな値に写像す る 2次元→1次元写像を施すことで 0以上の整数値 zを得る。 [0043] この 2次元→1次元写像は、原点に近い格子点ほど小さな値に写像するものであり
、例えば、図 1に示すような形で、整数値対 (X, y)を整数値 zに写像する。
[0044] この写像処理は、例えば、図 2に示すように、あらかじめ整数値対と整数値との対応 関係について記憶するテーブルを用意しておいて、整数値対 (X, y)をキーにしてそ のテーブルを参照することで、整数値対 (X, y)に対しての写像結果となる整数値 zを 得ることで fiうこと力 S可倉である。
[0045] ここで、このとき用意するテーブルは、 Huffman符号や算術符号を用いる場合に必 要とされる頻度表と異なって、符号化対象の情報源に依存せずに一般的なものとし て用意すること力できる。
[0046] また、この写像処理は、原点の格子点を起点にして、未処理の(列挙されていない) 格子点と原点との距離の最小値を求めて、その最小距離の格子点を規定の順番に 従って列挙して整数値を割り当てることを繰り返していくことで、整数値対 (X, y)につ いての写像結果となる整数値 zを得る、という方法で行うことが可能である。
[0047] ここで、この方法を用いると、信号値の上下限を仮定することなぐ整数値対 (X, y) と整数値 zとを対応付けることが可能である。
[0048] 本発明の情報源符号化方法で、原点に近い格子点ほど小さな値に写像するという
2次元→1次元写像を用いる理由は次の通りである。
[0049] すなわち、原点から距離 Lだけ離れた格子点の z値は、半径 Lの円周内の格子点数 にほぼ等しぐそれは円周の面積にほぼ等しい。これから、 "z 兀 "となる。なお、 これは大局的に見ての近似であり(実際の格子点の数はかなり大きい)、例えば z = 数百レベルにおいてほぼ等しくなる。
[0050] この格子点の確率 fは元々正規分布により" f exp (— aL2 ) (但し aは定数)"のよう に表せ、従って "z TT L2 "という関係から、 "f cc exp (— az/兀) "というように zについ ての指数分布となる。
[0051] このように、本発明の情報源符号化方法で用いる、原点に近い格子点ほど小さな 値に写像するとレ、う 2次元→1次元写像は、ガウス性信号源を指数分布の信号源に 変換する可逆写像を実現する。
[0052] これから、本発明の情報源符号化方法では、整数値対 (X, y)に対して、この 2次元 →1次元写像を施すことで整数値 zを得ると、 Golomb符号などのような指数分布に従 う情報源の符号化に用いられる符号を使って、その整数値 zを符号化する。
[0053] 前述したように、 Golomb符号は、指数分布に従う信号を効率よく符号化でき、かつ 符号化'復号のための表が不要で、処理が極めて簡易な瞬時復号可能な符号であり 、さらに、任意に大きな整数入力値を符号化することが可能という特長がある。
[0054] これから、本発明の情報源符号化方法によれば、ガウス性整数信号を簡易かつ効 串よく符号ィ匕することカでさるようになる。
[0055] また、本発明の情報源復号方法では、 z ,ζ ,ζ ,ζ,····のように入力される整数値の符
1 2 3 4
号化データを復号することで、本発明の情報符号化方法により符号化した整数値 Ζを 符号化データの入力順に復号する。
[0056] このとき、本発明の情報源符号化方法が Golomb符号を使って整数値の符号化デ ータを生成するときには、その符号化データを Golomb復号することで整数値を復号 する。
[0057] 続いて、その復号した整数値 zに対して、本発明の情報源符号化方法が用いた 2次 元→1次元写像の逆写像である 1次元→2次元写像を施すことで、その整数値の写 像元となった整数値対 (X, y)を復元する。
[0058] この 1次元→2次元写像は、本発明の情報源符号化方法が用いた 2次元→1次元 写像の写像処理を実行して、そのときに、復号した整数値 zが出現する場合に、それ に対応付けられる整数値対 (X, y)を特定することで行う。
[0059] これから、この写像処理は、本発明の情報源符号化方法の行う写像処理と同様に、 例えば、図 2に示すように、あらかじめ整数値対と整数値との対応関係について記憶 するテーブルを用意しておいて、整数値 zをキーにしてそのテーブルを参照すること で、整数値 zに対しての写像結果となる整数値対 (X, y)を得ることで行うことが可能で ある。
[0060] また、この写像処理は、原点の格子点を起点にして、列挙されていない格子点と原 点との距離の最小値を求めて、その最小距離の格子点を規定の順番に従って列挙 して整数値を割り当てることを繰り返していくことで、整数値 zに対しての写像結果とな る整数値対 (X, y)を得る、という方法で行うことが可能である。 [0061] 続!/、て、本発明の情報源復号方法は、復元した整数値対 (X, y)を構成する整数値 X, yを、その順番に出力する。
[0062] このように、本発明の情報源復号方法によれば、指数分布に従う情報源の符号化 に用いられる符号で符号化された符号化データを復号することで、ガウス性整数信 号を簡易かつ効率よく復号することができるようになる。
[0063] 以下、ガウス性信号源の二次拡大と一次元系列化を経て、ガウス性信号源を指数 性信号源に変換し、これを指数性信号源に適した性質のある Golomb符号により符 号化するとレ、う実施の形態に従って、本発明につ!/、て詳細に説明する。
[0064] 図 3に、本発明を具備する情報源符号化装置 1および情報源復号装置 2の装置構 成の一実施形態例を図示する。
[0065] この図に示すように、本発明の情報源符号化装置 1は、符号化対象となるガウス性 整数信号の信号値系列を入力する信号入力部 10と、信号入力部 10の入力した信 号値を、その入力順に二個ずつの整数値対に変換する整数値対変換部 11と、整数 値対変換部 11の変換した整数値対に対して、原点に近い格子点ほど小さな値に写 像する 2次元→1次元写像を施すことで 0以上の整数値を得る 2次元→1次元写像部 12と、 Golomb符号を使って、 2次元→1次元写像部 12の得た整数値を符号化する Golomb符号化部 13とを備える。
[0066] 一方、本発明の情報源復号装置 2は、本発明の情報源符号化装置 1の生成した整 数値の符号化データを入力する符号化データ入力部 20と、符号化データ入力部 20 の入力した整数値の符号化データを Golomb復号することで、その整数値を復号す る Golomb復号部 21と、本発明の情報源符号化装置 1の備える 2次元→1次元写像 部 12の写像処理を実行して、そのときに、 Golomb復号部 21の復号した整数値が出 現する場合に、それに対応付けられる整数値対を特定することで、 Golomb復号部 2 1の復号した整数値の写像元となった整数値対を復元する 1次元→2次元逆写像部 22と、 1次元→2次元逆写像部 22の復元した整数値対を構成する整数値を、その順 番に出力する順次出力部 23とを備えるように構成する。
[0067] 図 4に、図 3のように構成される本発明の情報源符号化装置 1の実行するフローチ ヤートの一例を図示する。 [0068] 次に、このフローチャートに従って、本発明の情報源符号化装置 1の実行する処理 について詳細に説明する。
[0069] 本発明の情報源符号化装置 1は、先ず最初に、ステップ S 101で、ガウス性整数信 号の信号値を、先頭からの順番に従って 2要素ずつ入力する。
[0070] 続いて、ステップ S 102で、 a ,a→a ,a→,····というような順に入力されるガウス性整
1 2 3 4
数信号を、(a , a ), (a , a ), ·····というように、その入力順に 2要素ずつの整数値対に
1 2 3 4
変換する。この整数値対を (X, y)とする。
[0071] 続いて、ステップ S 103で、この整数値対(X, y)に対して、原点に近い格子点ほど 小さな値に写像する 2次元→1次元写像を施すことで 0以上の整数値 zを得る。
[0072] 続いて、ステップ S 104で、別途与えられる Golomb符号パラメータ gを用いて、整数 値 zを Golomb符号化し、続くステップ S 105で、その符号化データを出力する。
[0073] 続いて、ステップ S 106で、ガウス性整数信号の入力が尽きたかどうかを判定し、尽 きていない場合は再びステップ S 101へ戻り、尽きている場合は終了する。
[0074] ここで、 Golomb符号パラメータ gについては別途与えるものとした力 符号化前に 入力信号の統計的性質を一旦調査し、その結果から適切に定めることも可能である
[0075] 入力信号の分散が σ 2であれば、例えば
g = 2 , log 2 · π · σ
e
により Golomb符号パラメータ gを定めることができる。但し、適宜整数丸めを施すもの とする。
[0076] 図 20で説明したように、 Golomb符号で符号化する場合、緩やかな分布の符号化 には Golomb符号パラメータ gを大きくすることが向き、 0に集中する急峻な分布の符 号化には Golomb符号パラメータ gを小さくすることが向くので、このような法則で Golo mb符号パラメータ gを定めると、適切なパラメータの値を自動的に設定できるという利 点、かある。
[0077] 次に、ステップ S 103で実行する 2次元→1次元写像について説明する。
[0078] この 2次元→1次元写像では、図 1にその一例を示すように、原点に近い格子点ほ ど小さい整数値 zに写像する。したがって、最初の 9個の zに対応する 1次元→2次元 逆写像は図 2に示すようになる。
[0079] 原点からの距離が等しレ、格子点は一般に複数存在するが、それらの順序付けは任 意である。図 1に示すように、反時計回りを保つように選ぶ、などが考えられる力 特 にそうでなくとも、復号側でも逆変換ができるように、順序付けが一意に決まってさえ いればよい。
[0080] 入力値 X, yの上下限値を仮定しておく必要はあるが、図 2に示すような適切な表を 用意しておけば、二次拡大'一次元系列化、およびその逆処理を表参照のみで簡易 に fiうことができる。
[0081] 実用上、入力の上下限値を仮定できる場合は多ぐ例えば、 8bit画像差分信号が 対象であれば、 - 255≤x, y≤ 255としてよい。この場合、
(255 - ( - 255) + 1) 2 = 5112 = 261 , 121
とレ、う要素の表を持つこととなる。
[0082] 一方、入力の上下限を仮定できない場合も、例えば、図 5に示すようなアルゴリズム により発生する擬似コードにより、格子点を原点に近い順に限りなく列挙していくこと で、(X, y) zの無限の対応づけが可能である。このアルゴリズムの詳細については 後述する。
[0083] ここで、図 5中に示す l (x, y)は、原点と点(x, y)の距離を与えるもので、距離として 通常はユークリッド距離
Kx, y) = (x2 + y ) 1/2
を用いる。また、より処理が簡易な二乗ユークリッド距離
Hx, y) =χ + y
を用いても、順序付けの結果は変わらない。
[0084] また、距離関数 l (x, y)は、必ずしも等方的でなくてもよぐ入力ガウス情報源が厳 密に iid(ind印 endent and identically distributed:独立かつ同一分布)でなく相関があ る場合、適当な定数 αを用い、
1 (x, y) = x + y + xy
としてもよい。ここで、 x, yに正の相関があれば αく 0、負の相関があれば α〉0とな [0085] また、入力信号対の二次元分布の形状に応じ、
Figure imgf000016_0001
という L lノルムであったり、より一般的に、下式に示す L yノルム(各信号要素の γ乗 )を用いてもよい。
[0086] [数 1]
!(xyy) - W +- [>·'
[0087] 次に、本発明の有効性を検証するために行った実験の結果について説明する。
[0088] 図 6に、実際の画像を処理して得られたガウス性信号に非常に近い信号源の頻度 分布を示す。縦軸は確率の対数である。この信号源の頻度分布が放物線状に分布 している、すなわち、正規分布に従っていることがわかる。
[0089] 図 7に、この信号源の信号を (X, y)の 2数値対として順次取り出し、その頻度分布 を 3次元表示したものを示す。同様に縦軸は確率の対数である。グラフは放物線を回 転した回転放物線面となる。
[0090] 図 8に、その回転放物線面の等高線を示す。わずかな正相関が観察されるがほぼ 同心円状となっている。
[0091] 縦対数グラフにおいて、指数分布は図 21に示したように三角形となる力 信号を、 刖述の
b = 2a - l a〉0の場合
b = - 2a a≤0の場合
という変換により片側分布としたものは右下がりの直線となる。
[0092] 図 9に、本発明に従って (X, y)に 2次元→1次元写像を施し 1次元値 zと変換したも のの頻度分布を、図 6と同様の縦対数グラフで示す。
グラフはほぼ右下がりの直線となっている、つまり指数分布がよくあてはまることがわ かる。従って、 Golomb符号による効率よい符号化が期待できる。
[0093] 実際の符号化結果を示す。用いた情報源はサンプル数 13, 509,440個、エントロピ 一から求まる情報量は 64,035,731 [bit]であった。これは理論値であり、実際の符号 化において符号量はこの値を必ず上回る。
[0094] この情報源に本発明の 2次元→1次元写像を施し、 Golomb符号パラメータ g = 7の Golomb符号で表現すると、 64,503, 706[bit]となり、 0.73%の増加(符号化効率 0.992
7)で符号化が可能であった。
[0095] 比較のための従来法として、同じ情報源に特段の写像を施さずに Golomb符号化 すると、 g= lの Golomb符号が最小値を与え 68,898,893[bit]必要であった。符号化 効率は 0.9294であった。
[0096] その他 4種類のデータについて同様の実験を行った。
[0097] 図 10に、得られた符号化効率の比較を示す。この図で data;!〜 data5のうち、 datal が上記のデータに対応している。本発明の符号化効率は最低でも約 90%を達成して おり、平均では、本発明は平均 94.1 %、従来方法は 87.4%と約 7ポイントの差があつ た。
[0098] このようにして、図 4のフローチャートに従ってガウス性整数信号を符号化する構成 を採る本発明の情報源符号化装置 1の符号化処理の有効性を検証できた。
[0099] 図 11に、図 3のように構成される本発明の情報源復号装置 2の実行するフローチヤ ートの一例を図示する。
[0100] 次に、このフローチャートに従って、本発明の情報源復号装置 2の実行する処理に ついて詳細に説明する。
[0101] 本発明の情報源復号装置 2は、先ず最初に、ステップ S201で、本発明の情報源 符号化装置 1により生成された符号化データ(整数値 zの Golomb符号)を、先頭から の順番に従って 1要素ずつ入力する。
[0102] 続いて、ステップ S202で、入力した Golomb符号を Golomb符号パラメータ gを用い て復号し、 0以上の整数値 zを得る。
[0103] 続いて、ステップ S203で、 1次元→2次元写像を施すことで整数値 zを整数値対(X
, y)に写像し、続くステップ S204で、この整数値対 (X, y)を X, yの順に出力する。
[0104] ここで、ステップ S203で実行する 1次元→2次元写像は、本発明の情報源符号化 装置 1が用いた 2次元→1次元写像の写像処理を実行して、そのときに、復号した整 数値 zが出現する場合に、それに対応付けられる整数値対 (X, y)を特定することで 行う。
[0105] 続いて、ステップ S205で、 Golomb符号の入力が尽きたかどうかを判定し、尽きて V、な!/、場合はステップ S 201に戻り、尽きて!/、る場合は終了する。
実施例
[0106] 次に、図 4のフローチャートのステップ S102で実行する処理の実施例と、図 4のフ ローチャートのステップ S 103で実行する処理の実施例と、図 11のフローチャートの ステップ S203で実行する処理の実施例と、図 11のフローチャートのステップ S204 で実行する処理の実施例について説明する。
[0107] 〔1〕図 4のフローチャートのステップ S 102で実行する処理の実施例
図 12に、図 4のフローチャートのステップ S102で実行する処理の詳細なフローチヤ ートを図示する。
[0108] このステップ S102では、ステップ S101の処理で入力された 2つのガウス性整数信 号の信号値を受け取ると、図 12のフローチャートに示すように、先ず最初に、ステツ プ S301で、そのうちの先頭の信号値 Xを入力してメモリ Xに格納し、続いて、ステップ S302で、もう一方の信号値 yを入力してメモリ yに格納する。
[0109] 続いて、ステップ S303で、メモリ Xから信号ィ直 Xを読み出すとともに、メモリ yから信号 値 yを読み出して、それらの X, yを整数値対 (X, y)として出力する。
[0110] このようにして、ステップ S 102では、図 12のフローチャートを実行することで、ガウ ス性整数信号の信号系列中の信号値を 2つずつ束ねて出力するという処理を行うの である。
[0111] 〔2〕図 4のフローチャートのステップ S 103で実行する処理の実施例
図 13〜図 16に、図 4のフローチャートのステップ S103で実行する処理の詳細なフ ローチャートを図示する。
[0112] このステップ S103では、図 5に示すアルゴリズムにより実現される 2次元→1次元写 像を施すことで整数値対 (X, y)を整数値 zに写像する処理を行って、その整数値 zを 出力する処理を行う。
[0113] すなわち、図 13のフローチャートに示すように、先ず最初に、ステップ S401で、写 像対象となる整数値 xO, yOを入力すると、続くステップ S402で、格子点記憶メモリ Z を空にするとともに、変数 zを 0に初期化する。
[0114] 続いて、ステップ S403で、図 14に示すフローチャートの処理で定義される手続き A を実行し、続く
ステップ S404で、図 15に示すフローチャートの処理で定義される手続き Bを実行す
[0115] この手続き Bの中で、図 16に示すフローチャートの処理で定義される手続き Xを実 行することになる力 この手続き Xでは、ある条件を満たすと、整数値対 (xO, yO)の 写像結果となる整数値 zが求められたことを判断して、その整数値 zを出力して、ステ ップ S405に進み、処理を終了する。
一方、その条件が満たされない場合には、何も出力せずに、ステップ S403に戻る ことで手続き Aに戻る処理を行う。
[0116] 以下に、これらの手続き A, B, Xについて詳細に説明する力 S、ここで簡単に説明す るならば、手続き Aでは、手続き Bで設定される " dmin≤x≤dmin, - dmin≤y≤ dmin"について、整数値を割り当てていない未列挙の格子点(x, y)を処理対象とし て、原点との間の距離の最小値 dminを求める。このとき、この最小値 dminを実現す る格子点は通常の場合複数となる。
これから、手続き Bでは、それらの格子点を規定の順番に従って 1つずつ選び出し て、それに対して前回よりも 1つ大きな整数値を割り当てて!/、く。
そして、手続き Xで、その割り当てに際して、入力した (xO, yO)が格子点 (X, y)とし て現れたのか否かを判断して、現れた場合には、それに対して割り当てた整数値 zを 写像結果として決定し、現れない場合には、手続き Bで dminの値を 1増やしてしてい くことで、入力した (xO, yO)が現れるまで、この処理を繰り返していくことになる。
[0117] 次に、図 14のフローチャートに従って、手続き Aの処理について説明する。
[0118] この手続き Aでは、 X, yがともに dmin以上 dmin以下(dminは整数で、手続き B により 1増やされていく)の 2次元範囲内の格子点(X, y)のうち、未列挙のものを処理 対象として原点からの距離の最小値を求めてそれを dminとする、という処理を行う。
[0119] 先ず最初に、ステップ S 501で、整数値 dminを 0に初期化する。続いて、ステップ S 502で、 Xを dminから dmin、 yを dminから dminまで 1間隔でもれなく規定の順 番に従って、 1ループ(S 502〜S508)に対して 1組ずつ発生させる。
[0120] 続いて、ステップ S503で、発生した (X, y)が格子点記憶メモリ Zに既に記憶されて いるのか否かを判断して、既に記憶されている場合には、ステップ S 508に進み、記 憶されていない場合には、ステップ S504に進んで、変数 dに対して、ステップ S505 の処理により求められる格子点 (X, y)と原点との間の距離 l (x, y)を代入する。
[0121] このとき用いる距離 Kx, y)としては、例えば、
l (x, y) =χ +y
という関数を用いる。
[0122] 続いて、ステップ S506で、 dと dminとの大小を比較して、 dが dmin以上である場合 には、直ちにステップ S508に進み、 dが dmin未満である場合には、ステップ S507 に進んで、 dminの値を dに設定してから、ステップ S508に進む。
[0123] そして、ステップ S508で、 Xと yがそれぞれ dminから dminまでの範囲で全ての 組み合わせが出尽くしたのか否かを判断して、まだ出尽くしていない場合には、ステ ップ S 502の処理に戻り、出尽くしていれば、手続き Aとしての処理を終了する。
[0124] 次に、図 15のフローチャートに従って、手続き Bの処理について説明する。
[0125] この手続き Bでは、 X, yがともに dmin以上 dmin以下(dminは 1増やしていく整 数)の 2次元範囲内の格子点(X, y)のうち、未列挙でかつ l (x, y) = dmin (dmin :手 続き Aで求められたもの)となるものについて、図 16のフローチャートに示す「手続き X」を実行する、という処理を行う。
[0126] 先ず最初に、ステップ S601で、フラグ foundを 0に設定する。続いて、ステップ S60 2で、 Xを dminから dmin、 yを dminから dminまで 1間隔でもれなく規定の順番 に従って、 1ノレープ(S602〜S609)に対して 1組ずつ発生させる。
[0127] 続!/、て、ステップ S603で、発生した (x, y)が格子点記憶メモリ Zに既に記憶されて いるのか否かを判断して、既に記憶されている場合には、ステップ S609に進み、記 憶されていない場合には、ステップ S604に進んで、変数 dに対して、ステップ S605 の処理により求められる格子点 (X, y)と原点との間の距離 l (x, y)を代入する。
[0128] このとき用いる距離 Kx, y)としては、例えば、
l (x, y) =χ +y
という関数を用いる。
[0129] 続!/、て、ステップ S606で、 dと dminとを比較し、 dと dminとが一致しな!/、場合には、 直ちにステップ S609に進み、 dと dminとが一致する場合には、ステップ S607に進ん で、手続き Xを実行する。
[0130] 続!/、て、ステップ S608で、格子点記憶メモリ Zに対して新たに整数値対 (x, y)を記 憶させ、別途定義される整数 zの値を 1増加させ、フラグ foundの値を 1に設定してか ら、ステップ S609に進む。
[0131] そして、ステップ S609で、 Xと yがそれぞれ dminから dminまでの範囲で全ての 組み合わせが出尽くしたのか否かを判断して、まだ出尽くしていない場合には、ステ ップ S602の処理に戻り、出尽くしていれば、ステップ S610に進んで、フラグ found の値が 0と等し!/、のか否かを判断する。
[0132] このステップ S610の判断処理に従って、フラグ foundの値が 0と等しくないことを判 断するときには、ステップ S601の処理に戻り、フラグ foundの値が 0と等しいことを判 断するときには、ステップ S611に進んで、 dminの値を 1増やして、手続き Bとしての 処理を終了する。
[0133] 次に、図 16のフローチャートに従って、手続き Xの処理について説明する。
[0134] この手続き Xでは、写像対象として入力された χθおよび yOの値力 S、現在ループで 回している Xおよび yの値とそれぞれ等しくなつた場合に、それに対応付けられる 1次 元写像値 zを出力する、という処理を行う。
[0135] 先ず最初に、ステップ S701で、 x = xO力、つ y=yOであるの力、を半 IJ断して、 x = xO力、 つ y=yOであることを判断する場合は、ステップ S702に進んで、(x, y)に割り当てら れた整数値 zを写像結果として出力して、ステップ S703に進んで、 2次元→1次元写 像処理を終了する。一方、 x = xOかつ y=yOでないことを判断する場合には、手続き Xの処理を終了する(即ち、手続き Bのステップ S608に戻る)。
[0136] このようにして、図 13〜図 16のフローチャートを実行することで、図 4のフローチヤ ートのステップ S 103では、図 5に示すアルゴリズムにより実現される 2次元→1次元写 像を施すことで整数値対 (X, y)を整数値 zに写像する処理を行って、その整数値 zを 出力する処理を行うのである。
[0137] 〔3〕図 11のフローチャートのステップ S203で実行する処理の実施例
図 17〜図 18に、図 11のフローチャートのステップ S203で実行する処理の詳細な フローチャートを図示する。
[0138] このステップ S203では、図 5に示すアルゴリズムにより実現される 1次元→2次元写 像を施すことで、ステップ S202の処理で復号された整数値 zを整数値対 (X, y)に写 像する処理を行って、その整数値対 (X, y)を出力する処理を行う。
[0139] すなわち、図 17のフローチャートに示すように、先ず最初に、ステップ S801で、写 像対象となる整数値 ζθを入力すると、続くステップ S802で、格子点記憶メモリ Zを空 にするとともに、変数 zを 0に初期化する。
[0140] 続いて、ステップ S803で、図 14に示すフローチャートの処理で定義される手続き A を実行し、続くステップ S804で、図 15に示すフローチャートの処理で定義される手 続き Bを実行する(ただし、「手続き X」力 以下に説明する「手続き X'」に変更される)
[0141] 即ち、この手続き Bの中で、図 18に示すフローチャートの処理で定義される手続き X'を実行することになる力 この手続き X'では、ある条件を満たすと、整数値 ζθの写 像結果となる整数値対 (X, y)が求められたことを判断して、その整数値対 (X, y)を出 力して、ステップ S805に進み、処理を終了する。
一方、その条件が満たされない場合には、何も出力せずに、ステップ S803に戻る ことで手続き Aに戻る処理を行う。
[0142] 次に、図 18のフローチャートに従って、手続き X'の処理について説明する。
[0143] この手続き X'では、写像対象として入力された ζθの値力 S、現在ループで回している zの値と等しくなつた場合に、それに対応付けられる 2次元写像値 (X, y)を出力する、 という処理を行う。
[0144] 先ず最初に、ステップ S901で、 z = zOであるのかを判断して、 z = zOであることを判 断する場合は、ステップ S902に進んで、整数値 zに対応付けられた整数値 (X, y)を 写像結果として出力して、ステップ S903に進んで、 1次元→2次元写像処理を終了 する。一方、 z = zOでないことを判断する場合には、手続き X'の処理を終了する(即 ち、手続き Bのステップ S608に戻る)。
[0145] このようにして、図 17〜図 18のフローチャートを実行することで、図 11のフローチヤ ートのステップ S203では、図 5に示すアルゴリズムにより実現される 1次元→2次元写 像を施すことで整数値 zを整数値対 (x, y)に写像する処理を行って、その整数値対( X, y)を出力する処理を行うのである。
[0146] 〔4〕図 11のフローチャートのステップ S204で実行する処理の実施例
図 19に、図 11のフローチャートのステップ S204で実行する処理の詳細なフローチ ヤートを図示する。
[0147] このステップ S204では、ステップ S203の処理で得られた整数値対(X, y)を受け取 ると、図 19のフローチャートに示すように、先ず最初に、ステップ S 1001で、そのうち の先頭の信号値 Xをメモリ Xに格納するとともに、もう一方の信号値 yをメモリ yに格納 する。
[0148] 続いて、ステップ S1002で、メモリ Xから先頭の信号値 Xを読み出して出力し、続くス テツプ S 1003で、メモリ yからもう一方の信号値 yを読み出して出力する。
[0149] このようにして、ステップ S204では、図 19のフローチャートを実行することで、写像 結果の整数値対 (X, y)を構成する整数値 X, yをその順番に従って出力するといぅ処 理を行うのである。
産業上の利用可能性
[0150] 本発明は、ガウス性整数信号を符号化'復号対象とするものであり、ガウス性信号 源を指数分布の信号源に変換する可逆写像を実現する構成を用いることで、数学や 工学の様々な場面で現れるにもかかわらず、従来の Golomb符号などでは効率的に 符号化ができなかったガウス性整数信号を簡易かつ効率よく符号化し復号すること ができる。

Claims

請求の範囲
[1] ガウス性整数信号を示す画像信号を符号化する画像信号符号化方法であって、 符号化対象の画像信号の信号値系列を入力するステップと、
前記入力した信号値系列に含まれる信号値を、その入力順に二個ずつの整数値 対とするステップと、
前記整数値対の各々を 2次元座標上の格子点と見て、原点に近!/、格子点ほど小さ な値に写像する 2次元→1次元写像を施すことで 0以上の整数値を得るステップと、
Golomb符号を使って、前記整数値を符号化するステップと
を有する画像信号符号化方法。
[2] 符号化対象のガウス性整数信号を示す画像信号の信号値系列をその入力順に二 個ずつの整数値対として、その整数値対の各々を 2次元座標上の格子点と見て、原 点に近い格子点ほど小さな値に写像する 2次元→1次元写像を施すことで 0以上の 整数値を得て、それを Golomb符号を使って符号化することにより生成された符号化 データを復号する画像信号復号方法であって、
前記整数値の符号化データを Golomb復号することで、前記整数値を復号するス 前記復号した整数値に対して、前記 2次元→1次元写像の逆写像である 1次元→2 次元写像を施すことで、前記整数値対を復元するステップと、
前記復元した整数値対を構成する整数値を、その順番に出力するステップと を有する画像信号復号方法。
ガウス性整数信号を符号化する情報源符号化方法であって、
符号化対象のガウス性整数信号の信号値系列を入力するステップと、
前記入力した信号値系列に含まれる信号値を、その入力順に二個ずつの整数値 対とするステップと、
前記整数値対の各々を 2次元座標上の格子点と見て、原点に近!/、格子点ほど小さ な値に写像する 2次元→1次元写像を施すことで 0以上の整数値を得るステップと、 指数分布に従う情報源の符号化に用いられる符号を使って、前記整数値を符号化 を有する情報源符号化方法。
[4] 請求項 3に記載の情報源符号化方法にお!/、て、
前記符号化するステップでは、 Golomb符号を使って前記整数値を符号化する情 報源符号化方法。
[5] 請求項 4に記載の情報源符号化方法にお!/、て、
前記入力した信号値の分散を算出するステップと、
前記算出した分散値に比例する値を持つ Golomb符号の符号パラメータを決定す を有する情報源符号化方法。
[6] 請求項 3に記載の情報源符号化方法にお!/、て、
前記整数値を得るステップでは、あらかじめ作成された前記整数値対と前記整数値 との対応関係を記憶するテーブルを参照することで、前記整数値対にっレ、ての写像 結果となる前記整数値を得る情報源符号化方法。
[7] 請求項 3に記載の情報源符号化方法にお!/、て、
前記整数値を得るステップでは、原点の格子点を起点にして、列挙されていない格 子点と原点との距離の最小値を求めて、その最小距離の格子点を規定の順番に従 つて列挙して整数値を割り当てることを繰り返して!/、くことで、前記整数値対につ!/、て の写像結果となる前記整数値を得る情報源符号化方法。
[8] 符号化対象のガウス性整数信号の信号値系列をその入力順に二個ずつの整数値 対として、その整数値対の各々を 2次元座標上の格子点と見て、原点に近い格子点 ほど小さな値に写像する 2次元→1次元写像を施すことで 0以上の整数値を得て、そ れを指数分布に従う情報源の符号化に用いられる符号を使って符号化することによ り生成された符号化データを復号する情報源復号方法であって、
前記整数値の符号化データを復号することで、前記整数値を復号するステップと、 前記復号した整数値に対して、前記 2次元→1次元写像の逆写像である 1次元→2 次元写像を施すことで、前記整数値対を復元するステップと、
前記復元した整数値対を構成する整数値を、その順番に出力するステップと を有する情報源復号方法。 請求項 8に記載の情報源復号方法において、
前記指数分布に従う情報源の符号化に用いられる符号として、 Golomb符号が用 いられる情報源復号方法。
請求項 9に記載の情報源復号方法において、
前記整数値の符号化データの復号に、前記 Golomb符号の Golomb符号パラメ一 タを用い、
前記信号値系列に含まれる信号値の分散に比例する値を持つものとして設定され た Golomb符号パラメータを入力するステップを有する情報源復号方法。
請求項 8に記載の情報源復号方法において、
前記整数値対を復元するステップでは、あらかじめ作成された前記整数値対と前記 整数値との対応関係を記憶するテーブルを参照することで、前記整数値にっレ、ての 逆写像結果となる前記整数値対を復元する情報源復号方法。
請求項 8に記載の情報源復号方法において、
前記整数値対を復元するステップでは、前記 2次元→1次元写像として、原点の格 子点を起点にして、列挙されていない格子点と原点との距離の最小値を求めて、そ の最小距離の格子点を規定の順番に従って列挙して整数値を割り当てることを繰り 返して!/、くことで、前記整数値対につ!/、ての写像結果となる前記整数値を得るとレ、う 写像が用いられる場合には、その写像結果に基づいて、前記整数値についての逆 写像結果となる前記整数値対を復元する情報源復号方法。
ガウス性整数信号を符号化する情報源符号化装置であって、
符号化対象のガウス性整数信号の信号値系列を入力する手段と、
前記入力した信号値系列に含まれる信号値を、その入力順に二個ずつの整数値 対とする手段と、
前記整数値対の各々を 2次元座標上の格子点と見て、原点に近!/、格子点ほど小さ な値に写像する 2次元→1次元写像を施すことで 0以上の整数値を得る手段と、 指数分布に従う情報源の符号化に用いられる符号を使って、前記整数値を符号化 する手段と
を有する情報源符号化装置。 [14] 請求項 13に記載の情報源符号化装置において、
前記符号化する手段は、 Golomb符号を使って前記整数値を符号化する情報源符 号化装置。
[15] 符号化対象のガウス性整数信号の信号値系列をその入力順に二個ずつの整数値 対として、その整数値対の各々を 2次元座標上の格子点と見て、原点に近い格子点 ほど小さな値に写像する 2次元→1次元写像を施すことで 0以上の整数値を得て、そ れを指数分布に従う情報源の符号化に用いられる符号を使って符号化することによ り生成された符号化データを復号する情報源復号装置であって、
前記整数値の符号化データを復号することで、前記整数値を復号する手段と、 前記復号した整数値に対して、前記 2次元→1次元写像の逆写像である 1次元→2 次元写像を施すことで、前記整数値対を復元する手段と、
前記復元した整数値対を構成する整数値を、その順番に出力する手段と を有する情報源復号装置。
[16] 請求項 15に記載の情報源復号装置において、
前記指数分布に従う情報源の符号化に用いられる符号として、 Golomb符号が用 いられる情報源復号装置。
[17] 請求項 3に記載の情報源符号化方法の実現に用いられる処理をコンピュータに実 行させるための情報源符号化プログラム。
[18] 請求項 3に記載の情報源符号化方法の実現に用いられる処理をコンピュータに実 行させるための情報源符号化プログラムを記録したコンピュータ読み取り可能な記録 媒体。
[19] 請求項 8に記載の情報源復号方法の実現に用いられる処理をコンピュータに実行 させるための情報源復号プログラム。
[20] 請求項 8に記載の情報源復号方法の実現に用いられる処理をコンピュータに実行 させるための情報源復号プログラムを記録したコンピュータ読み取り可能な記録媒体
PCT/JP2007/071719 2006-11-14 2007-11-08 Procédé d'encodage et procédé de décodage de signal d'image, procédé d'encodage et procédé de décodage de source d'information, dispositifs pour ceux-ci, leurs programmes, et support de mémoire avec programme enregistré WO2008059752A1 (fr)

Priority Applications (7)

Application Number Priority Date Filing Date Title
US12/446,818 US8300966B2 (en) 2006-11-14 2007-11-08 Image signal encoding method and decoding method, information source encoding method and decoding method, apparatuses therefor, programs therefor, and storage media which store the programs to process gaussian integer signal
BRPI0718239A BRPI0718239B1 (pt) 2006-11-14 2007-11-08 processo de codificação de fonte de informação, processo de decodificação de fonte de informação, aparelho e codificação de fonte de informação e aparelho de decodificação de fonte de informação
EP20070831450 EP2083560B1 (en) 2006-11-14 2007-11-08 Image signal coding method and decoding mehod, information source coding method and decoding mehod, devices for them, their programs, and memory medium with recorded program
CA2667727A CA2667727C (en) 2006-11-14 2007-11-08 Information source encoding method and decoding method, apparatuses therefor, programs therefor, and storage media which store the programs
ES07831450.7T ES2450265T3 (es) 2006-11-14 2007-11-08 Método de codificación y método de decodificación de señal de imagen, método de codificación y método de decodificación de fuente de información, dispositivos para ellos, sus programas, y medio de memoria con programa registrado
CN2007800412601A CN101536484B (zh) 2006-11-14 2007-11-08 图像信号编码方法及解码方法、信息源编码及解码方法、它们的装置
JP2008544115A JP4769305B2 (ja) 2006-11-14 2007-11-08 画像信号符号化方法及び復号方法、情報源符号化及び復号方法、それらの装置、及びそれらのプログラム並びにプログラムを記録した記憶媒体

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2006-307512 2006-11-14
JP2006307512 2006-11-14

Publications (1)

Publication Number Publication Date
WO2008059752A1 true WO2008059752A1 (fr) 2008-05-22

Family

ID=39401559

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2007/071719 WO2008059752A1 (fr) 2006-11-14 2007-11-08 Procédé d'encodage et procédé de décodage de signal d'image, procédé d'encodage et procédé de décodage de source d'information, dispositifs pour ceux-ci, leurs programmes, et support de mémoire avec programme enregistré

Country Status (11)

Country Link
US (1) US8300966B2 (ja)
EP (1) EP2083560B1 (ja)
JP (1) JP4769305B2 (ja)
KR (1) KR101010938B1 (ja)
CN (1) CN101536484B (ja)
BR (1) BRPI0718239B1 (ja)
CA (1) CA2667727C (ja)
ES (1) ES2450265T3 (ja)
RU (1) RU2406222C1 (ja)
TW (1) TWI335147B (ja)
WO (1) WO2008059752A1 (ja)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010005359A1 (en) * 2008-07-10 2010-01-14 Ali Majeed Data compression involving spiral transformation
WO2018235418A1 (ja) * 2017-06-22 2018-12-27 日本電信電話株式会社 符号化装置、復号装置、符号化方法、復号方法、およびプログラム
WO2019021552A1 (ja) * 2017-07-25 2019-01-31 日本電信電話株式会社 符号化装置、復号装置、符号列のデータ構造、符号化方法、復号方法、符号化プログラム、復号プログラム
WO2019198383A1 (ja) * 2018-04-13 2019-10-17 日本電信電話株式会社 符号化装置、復号装置、符号化方法、復号方法、プログラム、および記録媒体
WO2020095706A1 (ja) * 2018-11-07 2020-05-14 日本電信電話株式会社 符号化装置、復号装置、符号列のデータ構造、符号化方法、復号方法、符号化プログラム、復号プログラム

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8817882B2 (en) 2010-07-30 2014-08-26 Qualcomm Incorporated Coding blocks of data using a generalized form of golomb codes
EP3621206B1 (en) * 2011-10-24 2024-06-05 Velos Media International Limited Significant map decoding using partition selection
US8891630B2 (en) 2011-10-24 2014-11-18 Blackberry Limited Significance map encoding and decoding using partition set based context assignment
CN104038234B (zh) 2013-03-07 2017-09-29 华为技术有限公司 极性码的译码方法和译码器
GB201321851D0 (en) 2013-12-10 2014-01-22 Canon Kk Run length parameters coding for palette mode
CN105009461B (zh) * 2013-12-24 2017-12-22 华为技术有限公司 极性码的译码方法和译码装置
EP3761313B1 (en) * 2018-03-02 2023-01-18 Nippon Telegraph And Telephone Corporation Encoding device, encoding method, program, and recording medium
CN114124413A (zh) * 2020-08-10 2022-03-01 中国电信股份有限公司 媒体文件加密方法、解密方法、传输方法和装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10243395A (ja) * 1997-02-28 1998-09-11 Graphics Commun Lab:Kk 画像の符号化および復号方法と装置
JP2001005967A (ja) 1999-06-21 2001-01-12 Matsushita Electric Ind Co Ltd 画像送信装置及びニューラルネットワーク

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5303265A (en) 1990-06-29 1994-04-12 Digital Equipment Corporation Frequency independent encoding technique and apparatus for digital communications
US5471207A (en) * 1994-02-23 1995-11-28 Ricoh Company Ltd. Compression of palettized images and binarization for bitwise coding of M-ary alphabets therefor
CA2156889C (en) 1994-09-30 1999-11-02 Edward L. Schwartz Method and apparatus for encoding and decoding data
US6137589A (en) * 1995-09-20 2000-10-24 Hewlett-Packard Company Efficient storage of dithered raster image data in a printer or the like
US5764374A (en) 1996-02-05 1998-06-09 Hewlett-Packard Company System and method for lossless image compression having improved sequential determination of golomb parameter
JP2000115782A (ja) 1998-10-06 2000-04-21 Canon Inc 符号化装置及び方法及び記憶媒体
KR100317279B1 (ko) 1998-11-04 2002-01-15 구자홍 무손실부호화방법및장치
US6668092B1 (en) * 1999-07-30 2003-12-23 Sun Microsystems, Inc. Memory efficient variable-length encoding/decoding system
KR100399932B1 (ko) 2001-05-07 2003-09-29 주식회사 하이닉스반도체 메모리의 양을 감소시키기 위한 비디오 프레임의압축/역압축 하드웨어 시스템
US7342965B2 (en) 2002-04-26 2008-03-11 Nokia Corporation Adaptive method and system for mapping parameter values to codeword indexes
US6987468B1 (en) * 2004-10-29 2006-01-17 Microsoft Corporation Lossless adaptive encoding and decoding of integer data
US7015837B1 (en) 2004-10-29 2006-03-21 Microsoft Corporation Lossless adaptive encoding and decoding of integer data
US7580585B2 (en) 2004-10-29 2009-08-25 Microsoft Corporation Lossless adaptive Golomb/Rice encoding and decoding of integer data using backward-adaptive rules
JP4468858B2 (ja) 2005-06-01 2010-05-26 オリンパスイメージング株式会社 データ符号化装置、データ符号化方法、プログラム

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10243395A (ja) * 1997-02-28 1998-09-11 Graphics Commun Lab:Kk 画像の符号化および復号方法と装置
JP2001005967A (ja) 1999-06-21 2001-01-12 Matsushita Electric Ind Co Ltd 画像送信装置及びニューラルネットワーク

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
P BOLDI; S VIGNA: "Compressed perfect embedded skip lists for quick inverted-index lookups", PROCEEDINGS OF STRING PROCESSING AND INFORMATION RETRIEVAL, 12TH INTERNATIONAL CONFERENCE, 2005, pages 1 - 15
S XUE; Y XU; B OELMANN: "Hybrid Golomb codes for a group of quantised GG sources", IEE PROCEEDINGS, VISION IMAGE AND SIGNAL PROCESSING, vol. 150, no. 4, 2003, pages 256 - 260, XP006021077, DOI: doi:10.1049/ip-vis:20030510
See also references of EP2083560A4 *

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010005359A1 (en) * 2008-07-10 2010-01-14 Ali Majeed Data compression involving spiral transformation
WO2018235418A1 (ja) * 2017-06-22 2018-12-27 日本電信電話株式会社 符号化装置、復号装置、符号化方法、復号方法、およびプログラム
CN110771045A (zh) * 2017-06-22 2020-02-07 日本电信电话株式会社 编码装置、解码装置、编码方法、解码方法、以及程序
US12028093B2 (en) 2017-06-22 2024-07-02 Nippon Telegraph And Telephone Corporation Encoder, decoder, encoding method, decoding method and program
CN110771045B (zh) * 2017-06-22 2024-03-29 日本电信电话株式会社 编码装置、解码装置、编码方法、解码方法、以及记录介质
JPWO2018235418A1 (ja) * 2017-06-22 2020-05-21 日本電信電話株式会社 符号化装置、復号装置、符号化方法、復号方法、およびプログラム
US10840944B2 (en) 2017-07-25 2020-11-17 Nippon Telegraph And Telephone Corporation Encoding apparatus, decoding apparatus, data structure of code string, encoding method, decoding method, encoding program and decoding program
WO2019021552A1 (ja) * 2017-07-25 2019-01-31 日本電信電話株式会社 符号化装置、復号装置、符号列のデータ構造、符号化方法、復号方法、符号化プログラム、復号プログラム
JPWO2019021552A1 (ja) * 2017-07-25 2020-04-09 日本電信電話株式会社 符号化装置、復号装置、符号列のデータ構造、符号化方法、復号方法、符号化プログラム、復号プログラム
WO2019198383A1 (ja) * 2018-04-13 2019-10-17 日本電信電話株式会社 符号化装置、復号装置、符号化方法、復号方法、プログラム、および記録媒体
JPWO2019198383A1 (ja) * 2018-04-13 2021-04-08 日本電信電話株式会社 符号化装置、復号装置、符号化方法、復号方法、プログラム、および記録媒体
JP7173134B2 (ja) 2018-04-13 2022-11-16 日本電信電話株式会社 符号化装置、復号装置、符号化方法、復号方法、プログラム、および記録媒体
JP7139897B2 (ja) 2018-11-07 2022-09-21 日本電信電話株式会社 符号化装置、復号装置、符号化方法、復号方法、符号化プログラム、復号プログラム
JP2020077969A (ja) * 2018-11-07 2020-05-21 日本電信電話株式会社 符号化装置、復号装置、符号列のデータ構造、符号化方法、復号方法、符号化プログラム、復号プログラム
WO2020095706A1 (ja) * 2018-11-07 2020-05-14 日本電信電話株式会社 符号化装置、復号装置、符号列のデータ構造、符号化方法、復号方法、符号化プログラム、復号プログラム

Also Published As

Publication number Publication date
KR101010938B1 (ko) 2011-01-25
EP2083560B1 (en) 2013-12-04
KR20090084841A (ko) 2009-08-05
US8300966B2 (en) 2012-10-30
ES2450265T3 (es) 2014-03-24
TWI335147B (en) 2010-12-21
JPWO2008059752A1 (ja) 2010-03-04
EP2083560A1 (en) 2009-07-29
CN101536484B (zh) 2011-10-05
JP4769305B2 (ja) 2011-09-07
RU2406222C1 (ru) 2010-12-10
CA2667727C (en) 2014-07-08
CN101536484A (zh) 2009-09-16
BRPI0718239A2 (pt) 2013-11-19
CA2667727A1 (en) 2008-05-22
EP2083560A4 (en) 2011-12-28
BRPI0718239B1 (pt) 2019-12-24
US20100027903A1 (en) 2010-02-04
TW200830732A (en) 2008-07-16

Similar Documents

Publication Publication Date Title
WO2008059752A1 (fr) Procédé d'encodage et procédé de décodage de signal d'image, procédé d'encodage et procédé de décodage de source d'information, dispositifs pour ceux-ci, leurs programmes, et support de mémoire avec programme enregistré
US11632560B2 (en) Methods and apparatuses for encoding and decoding a bytestream
TWI239489B (en) System and method for progressively transforming and coding digital data
US20190228050A1 (en) Bezier volume representation of point cloud attributes
Kavitha A survey on lossless and lossy data compression methods
JP6242074B2 (ja) 信号データの圧縮および圧縮解除のための方法および装置(信号データの圧縮および圧縮解除)
US9543980B2 (en) Systems and methods for model-free compression and model-based decompression
US20130159263A1 (en) Lossy compression of data points using point-wise error constraints
JPH06224777A (ja) 符号化方法、符号化装置、復号方法、復号器、データ圧縮装置、ビットストリーム生成方法及び遷移マシン生成方法
JP6246233B2 (ja) 頂点誤差訂正のための方法および装置
US10366698B2 (en) Variable length coding of indices and bit scheduling in a pyramid vector quantizer
CN116156072B (zh) 隐写图像生成方法、隐写信息提取方法及相关装置
CN116128985A (zh) 基于双分支神经网络的大规模点云几何压缩方法
Balado et al. Asymptotically optimum perfect universal steganography of finite memoryless sources
JP3603070B2 (ja) 圧縮コード生成方法、および圧縮コード展開方法
Maleika et al. Lossless compression method for digital terrain model of seabed shape
Al-Azawi et al. Compression of Audio Using Transform Coding.
Negi et al. A review on fractal image compression
JP4694511B2 (ja) 符号装置、復号装置、符号化復号化システム、符号化復号化方法、符号化プログラム、復号化プログラム、符号化復号化プログラム及び記録媒体
WO2023132329A1 (ja) 点群復号装置、点群復号方法及びプログラム
Nandi et al. Fractal image compression using Fast Context Independent HV partitioning scheme and MRBH, MRBHM coding
WO2002058038A1 (fr) Appareil, procede et programme de cryptage a compression arithmetique
Bakulina Efficient image compression by coding of low-entropy sources
KR20240009371A (ko) 메모리 장치, 데이터 압축 및 압축 해제 방법
Gupta et al. Adaptive Byte Compression and Decompression: A New Approach for Fractal Image Compression

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 200780041260.1

Country of ref document: CN

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 07831450

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2008544115

Country of ref document: JP

WWE Wipo information: entry into national phase

Ref document number: 2007831450

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 2667727

Country of ref document: CA

WWE Wipo information: entry into national phase

Ref document number: 12446818

Country of ref document: US

Ref document number: 2253/CHENP/2009

Country of ref document: IN

WWE Wipo information: entry into national phase

Ref document number: 1020097008940

Country of ref document: KR

ENP Entry into the national phase

Ref document number: 2009116946

Country of ref document: RU

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: PI0718239

Country of ref document: BR

Kind code of ref document: A2

Effective date: 20090505