JPH0827835B2 - Parallel projection information generation method - Google Patents

Parallel projection information generation method

Info

Publication number
JPH0827835B2
JPH0827835B2 JP61171211A JP17121186A JPH0827835B2 JP H0827835 B2 JPH0827835 B2 JP H0827835B2 JP 61171211 A JP61171211 A JP 61171211A JP 17121186 A JP17121186 A JP 17121186A JP H0827835 B2 JPH0827835 B2 JP H0827835B2
Authority
JP
Japan
Prior art keywords
projection
value
projection matrix
residual
point
Prior art date
Legal status (The legal status 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 status listed.)
Expired - Lifetime
Application number
JP61171211A
Other languages
Japanese (ja)
Other versions
JPS6326781A (en
Inventor
哲 高山
千香 小野寺
成夫 島崎
桂 川上
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Panasonic Holdings Corp
Original Assignee
Matsushita Electric Industrial Co Ltd
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 Matsushita Electric Industrial Co Ltd filed Critical Matsushita Electric Industrial Co Ltd
Priority to JP61171211A priority Critical patent/JPH0827835B2/en
Publication of JPS6326781A publication Critical patent/JPS6326781A/en
Publication of JPH0827835B2 publication Critical patent/JPH0827835B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Landscapes

  • Image Processing (AREA)

Description

【発明の詳細な説明】 産業上の利用分野 本発明は計算機による画像の拡大縮小の処理を制御す
る情報を、高速に算出する投影情報並列生成方法に関す
る。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a projection information parallel generation method for calculating information for controlling image enlargement / reduction processing by a computer at high speed.

従来の技術 画像の拡大縮小は、原画A(a1、a2……、am)に対し
て結果B(b1、b2……bn)を得る処理、即ち、原画素の
インデックスであるiに対応すべき、結果画素のインデ
ックスjを算出し、aiの値をbjに代入する作業である。
(ただし、ai、bjは、それぞれ原画及び結果の各画素の
値である。)この拡大縮小の処理において、原画素aiを
直線に写影することにより、結果画素bjを求める方法が
ある。以下、これを説明する。
Scaling of the prior art image original A (a 1, a 2 ...... , a m) for the results B (b 1, b 2 ...... b n) processing to obtain, i.e., the index of the original pixel This is an operation of calculating an index j of a result pixel that should correspond to a certain i and substituting the value of ai for bj.
(However, ai and bj are the values of each pixel of the original image and the result.) In this scaling process, there is a method of obtaining the result pixel bj by mapping the original pixel ai into a straight line. This will be described below.

第11図に、原画素mビットを、結果画素nビットに縮
小(m≧n)する場合を示す。x軸方向に原画素A
(a1、a2……、am)の並びを、y軸方向に結果画素B
(b1、b2……、bn)の並びを、結果画素間の境界に、水
平格子線hを想定する。hのインデックスjは、すぐ上
に存在する結果画素と同一の値を持つものとする。図中
の直線Lは、縮小率n/mを傾きとし、量子化の誤差を表
わす値(以下、残差と呼ぶ)の初期値R0をy切片に持つ
直線である。原画Aから結果Bへの対応は、この直線L
による写影から求めることができる。
FIG. 11 shows a case where the original pixel m bits are reduced to the result pixel n bits (m ≧ n). Original pixel A in the x-axis direction
(A 1 , a 2, ..., Am ) are arranged in the y-axis direction as a result pixel B.
A horizontal grid line h is assumed at the boundary between the result pixels with the arrangement of (b 1 , b 2, ..., B n ). It is assumed that the index j of h has the same value as the result pixel immediately above. A straight line L in the drawing is a straight line having a reduction ratio n / m as an inclination and an initial value R 0 of a value (hereinafter, referred to as a residual) representing a quantization error as a y intercept. The correspondence from the original picture A to the result B is this straight line L
It can be obtained from the map.

原画素aiの中心を通る垂線と、直線Lとの交点をPiと
する。点Piを通る水平線が貫く結果画素をbjとすると
き、原画素のインデックスiが、結果画素のインデック
スjに対応するものとし、aiの値をbjに転送する。aiが
bjに対応することが予めわかっているとき、aiの隣接画
素ai+1に対応する結果画素を求める方法は簡略化でき
る。即ち、ai+1は、bjまたはbj+1のどちらかに対応
するが、そのどちらであるかは、点Pi+1との間に水平
格子線hj+1があるか否かを確めることによって決定で
きる。
The intersection of the perpendicular line passing through the center of the original pixel ai and the straight line L is Pi. When the result pixel penetrated by the horizontal line passing through the point Pi is bj, the index i of the original pixel corresponds to the index j of the result pixel, and the value of ai is transferred to bj. ai
When it is known in advance that it corresponds to bj, the method of obtaining the result pixel corresponding to the adjacent pixel ai + 1 of ai can be simplified. That is, ai + 1 corresponds to either bj or bj + 1, which can be determined by ascertaining whether or not there is a horizontal grid line hj + 1 between the point and Pi + 1.

このように、直線Lに写影する方法で、原画素と結果
画素の対応関係を求めることが可能であるが、その対応
関係の情報を、原画または結果画のうち、画素数の多い
方と同じ要素を持つベクトルT(t1、t2、……、tm
(以下投影ベクトルと呼ぶ)によって表わすことができ
る。投影ベクトルTの要素tiは、原画素aiがどの結果画
素に対応するかの情報をもっており、以下の手順により
定まるものとする。
In this way, it is possible to obtain the correspondence relationship between the original pixel and the result pixel by the method of mapping on the straight line L. Vector T (t 1 , t 2 , ……, t m ) with the same elements
(Hereinafter referred to as a projection vector). The element ti of the projection vector T has information about which result pixel the original pixel ai corresponds to, and is determined by the following procedure.

一担、aiとbjが対応していることがわかっている場
合、ai+1はbjか、あるいはbj+1のどちらかに対応す
る。たとえば、水平格子線hj+1が点Piと点Pi+1の間
に介在する場合、ai+1はbj+1に対応するものとす
る。また水平格子線hj+1が点Piと点Pi+1の間に存在
しない場合、ai+1はbjに対応するものとする。前者の
場合にti+1の値を1とし、後者の場合にti+1の値を
0とする。すなわち、ti+1の値が1となるときは、点
pi+1が水平格子線hj+1の上位に、もしくは、水平格
子線hj+1上に存在する場合であり、ti+1の値が0と
なるときは、その他の場合である。t1の値が1となると
きは、点p1が水平格子線h1の上位に、もしくは、水平格
子線h1上に存在する場合であり、t1の値が0となるとき
はその他の場合である。例えば、原画素7ビットを結果
画素3ビットに縮小する場合を、第12図をもとに説明す
る。前述したように、x軸、y軸それぞれに原画素A
(a1、a2、……、a7)結果画素B(b1、b2、b3)の並び
を想定する。又、結果画素間に、水平格子線h1、h2、h3
を描く。水平線hのインデックスjは、第11図に示した
ように、すぐ上に存在する結果画素と同一の値をもつ。
図中の直線Lの傾きは3/7である。この傾き3/7は、原画
素数を分母に、結果画素数を分子とした値である。又、
直線Lのy切片であるR0は、前述した量子化の誤差を表
わす値の初期値、つまり残差初期値である。
If it is known that ai and bj correspond to each other, ai + 1 corresponds to either bj or bj + 1. For example, if the horizontal grid line hj + 1 is interposed between the points Pi and Pi + 1, then ai + 1 corresponds to bj + 1. When the horizontal grid line hj + 1 does not exist between the points Pi and Pi + 1, ai + 1 corresponds to bj. In the former case, the value of ti + 1 is 1, and in the latter case, the value of ti + 1 is 0. That is, when the value of ti + 1 becomes 1, the point
This is a case where pi + 1 exists above the horizontal grid line hj + 1 or on the horizontal grid line hj + 1, and when the value of ti + 1 becomes 0, it is the other case. When the value of t 1 is 1, it means that the point p 1 is above the horizontal grid line h 1 or on the horizontal grid line h 1 , and when the value of t 1 is 0 Is the case. For example, a case where 7 bits of original pixels are reduced to 3 bits of result pixels will be described with reference to FIG. As described above, the original pixel A is set on each of the x-axis and the y-axis.
(A 1 , a 2 , ..., A 7 ) Assume the arrangement of the result pixels B (b 1 , b 2 , b 3 ). In addition, horizontal grid lines h 1 , h 2 , h 3 between the result pixels
Draw. The index j of the horizontal line h has the same value as the result pixel immediately above, as shown in FIG.
The slope of the straight line L in the figure is 3/7. This slope 3/7 is a value with the original pixel number as the denominator and the resulting pixel number as the numerator. or,
R 0 , which is the y-intercept of the straight line L, is the initial value of the value representing the quantization error described above, that is, the residual initial value.

原画素の中心を通る垂線と直線Lとの交点をそれぞれ
左からp1、p2、……、p7とする。t1の値は、点p1が、水
平格子線h1の上位に位置するか否かで判断する。第12図
において点p1は、水平格子線h1の上位に位置している。
故に t1=1 である。
The intersections of the perpendicular line passing through the center of the original pixel and the straight line L are p 1 , p 2 , ..., P 7 from the left. The value of t 1 is determined by whether or not the point p 1 is located above the horizontal grid line h 1 . In FIG. 12, the point p 1 is located above the horizontal grid line h 1 .
Therefore t 1 = 1.

次にt2の値を定める。投影ベクトルTの第2番目以降
の値を求める場合は、前述したようにその1つ前の原画
素が、どの結果画素に対応しているかという情報を必要
とする。第12図において点p1は、水平格子線h1の上位に
位置していることより、原画素a4は、結果画素b1に対応
していることがわかっている。このことより、t2の値
は、点p2が水平格子線h2の上位に位置しているか、また
は水平格子線h2上にある場合に1、その他の場合に0と
なる。故に t2=0 である。
Next, determine the value of t 2 . When obtaining the second and subsequent values of the projection vector T, as described above, it is necessary to have information as to which result pixel the original pixel immediately before the projection vector T corresponds to. In FIG. 12, the point p 1 is located above the horizontal grid line h 1 , and it is known that the original pixel a 4 corresponds to the result pixel b 1 . From this, the value of t 2 is 1, the other becomes zero in the case of when the point p 2 is either located on top of the horizontal grid lines h 2, or on the horizontal grid lines h 2. Therefore t 2 = 0.

次にt3の値を考える。原画素a2は、結果画素b1に対応
しているので、点p3が、水平格子線h2の上位に位置して
いるか否かで判断すればよい。第12図に於いて点p3は、
水平格子線h2の上位に位置している。故に t3=1 である。
Next, consider the value of t 3 . Since the original pixel a 2 corresponds to the result pixel b 1 , it may be determined whether or not the point p 3 is located above the horizontal grid line h 2 . In Fig. 12, the point p 3 is
It is located above the horizontal grid line h 2 . Therefore t 3 = 1.

以上のように順にt7まで求めてゆくと、投影ベクトル
Tは T=(1010010) ……(1) となる。
When t 7 is sequentially obtained as described above, the projection vector T becomes T = (1010010) (1).

点Piのy座標の値をyiとすると、この点に於ける残差
Riは、 Ri=yi−i ……(2) である。(但し、iは、点Piの直下にある水平線hiの
高さである。また、第12図のR5は残差を示す。) 以上のように投影ベクトル、及び残差を算出する。拡
大の場合は、原画素をB(b1、b2、……、bn)、結果画
素をA(a1、a2、……、am)として対応させる場合であ
り、投影ベクトルは、同一である。
Let yi be the y-coordinate value of point Pi, and the residual at this point
Ri is Ri = yi-i (2). (However, i is the height of the horizontal line hi immediately below the point Pi. In addition, R 5 in FIG. 12 represents the residual.) The projection vector and the residual are calculated as described above. In the case of enlargement, the original pixel is made to correspond to B (b 1 , b 2 , ..., B n ) and the resulting pixel is made to correspond to A (a 1 , a 2 , ..., Am ), and the projection vector is , Are the same.

投影ベクトルの計算は、原画素aiの中心からひいた垂
線と直線Lとの交点piを、格子点に対応させ量子化する
計算であり、これは従来DDA(デジタル微分解析)アル
ゴリズムにより1個の演算器を使用して、逐次処理によ
り求めていた。実際の計算過程を、以下に説明する。
The projection vector is calculated by quantizing the intersection point pi of the perpendicular line drawn from the center of the original pixel ai and the straight line L in correspondence with the grid point, which is calculated by the conventional DDA (digital differential analysis) algorithm. It was obtained by sequential processing using an arithmetic unit. The actual calculation process will be described below.

縮小率n/mの場合には格子線間の長さをmとし、残差R
iにnを加算した値と、mとの大小比較を行う。上記演
算結果であるRi+nがmより大きいか、又は等しい場
合、投影ベクトルTの第i+1番目の要素であるti+1
の値を1とし、次の残差Ri+1の値をRi+n−mとす
る。小さい場合は、ti+1の値を0とし、次の残差Ri+
1の値をRi+nとする。(なお、このアルゴリズムの詳
細は、ケー カワカミ アンド エス シマザキ(K.Ka
wakami and S.Simazaki)、“エー スペシャル パー
パス エルエスアイ プロセッサー ユージイング ザ
デーデーエー アルゴリズム フォー イメージ ト
ランスフォーメーション(A Special Purpose LSI Proc
essor Using The DDA Algrithm For Image Transformat
ion)”アイイーイーイー(IEEE)、コム アーク(Com
Arc)、11th、1984に詳しく述べられている。) 次に、第13図に、このフローを示す。80〜82は要素カ
ウンタ、及び残差レジスタの初期値設定を行う処理であ
る。83〜86は投影ベクトルの要素1個と残差を算出する
処理である。87は原画素カウンタ、投影ベクトル要素ポ
インタを1だけ増加する処理である。88は投影ベクトル
の要素をm個算出したかどうかの判断を行う処理であ
る。否であれば、82にもどり、前記の処理を繰り返す。
When the reduction ratio is n / m, the length between the grid lines is m and the residual R
The value obtained by adding n to i is compared with m. When the calculation result Ri + n is larger than or equal to m, the i + 1-th element of the projection vector T, ti + 1
Is 1, and the value of the next residual Ri + 1 is Ri + n−m. If it is smaller, the value of ti + 1 is set to 0 and the next residual Ri +
The value of 1 is Ri + n. (For details of this algorithm, see K. Kawakami & S. Shimazaki (K.Ka
wakami and S. Simazaki), "A Special Purpose LSI Processor Equipped the Dayda Algorithm for Image Transformation (A Special Purpose LSI Proc
essor Using The DDA Algrithm For Image Transformat
ion) ”IEEE, Com Ark
Arc), 11th, 1984. ) Next, FIG. 13 shows this flow. 80 to 82 are processes for setting initial values of the element counter and the residual register. 83 to 86 are processes for calculating one element of the projection vector and the residual. 87 is a process of incrementing the original pixel counter and the projection vector element pointer by 1. Reference numeral 88 is a process of determining whether or not m elements of the projection vector have been calculated. If not, the process returns to 82 and the above process is repeated.

以上のように、従来は1個毎に投影ベクトルの要素と
残差をもとめるという処理により行なうことができる。
As described above, conventionally, the process can be performed by obtaining the elements and residuals of the projection vector for each one.

発明が解決しようとする問題点 しかし従来の上述したような方法では、前述の投影ベ
クトルの生成をプログラムにより逐次的に行なっていた
ために、多数の演算回路を必要とし、時間がかかり、拡
大縮小の処理の高速化を図るのが著しく困難であった。
Problems to be Solved by the Invention However, in the conventional method as described above, since the above-described projection vector generation is sequentially performed by a program, a large number of arithmetic circuits are required, and it takes time and enlargement / reduction is required. It was extremely difficult to increase the processing speed.

本発明は上記従来技術に鑑み、拡大縮小の処理の高速
化を図るものである。
In view of the above-mentioned conventional technique, the present invention aims to speed up the scaling process.

問題点を解決するための手段 本発明は、J個の画素からなる第1の画素列と、K個
の画素からなる第2の画素列との直線への写影を介した
対応関係を表すLビット(LはJ,Kのうち大きな数)か
らなる投影ベクトルを、一度にPビット(L>P)づつ
並列に求める処理を、第1の画素列と第2の画素列との
縮小率n/m<1を決定するn,m、及び与えられた量子化残
差から、P行P列の投影行列を求める投影行列生成処理
と、前記投影行列生成処理に続く、次の投影行列生成処
理に用いる量子化残差のP+1個の候補をn,m、及び前
記投影行列生成処理で用いられた量子化残差から求める
残差候補生成処理と、前記投影行列生成処理によって得
られた投影行列からP個の投影ベクトルの要素を決定す
る投影ベクトル決定処理と、前記残差候補生成処理によ
って生成されたP+1個の残差の候補と、前記投影行列
生成処理により得られたP行P列の投影行列の第P行の
要素から次の投影行列生成処理に用いる量子化残差を決
定する残差決定処理とを、Q回(Q≧1)行うことによ
り実現することを特徴とするものである。
Means for Solving the Problems The present invention represents a correspondence relationship between a first pixel column composed of J pixels and a second pixel column composed of K pixels via a mapping to a straight line. A process of obtaining a projection vector composed of L bits (L is a large number of J and K) P bits (L> P) at a time in parallel is performed by a reduction ratio between the first pixel row and the second pixel row. A projection matrix generation process for obtaining a projection matrix of P rows and P columns from n, m for determining n / m <1 and a given quantization residual, and the next projection matrix generation following the projection matrix generation process. Residual candidate generation processing for obtaining P + 1 candidates of quantized residuals used in processing from n, m and the quantized residual used in the projection matrix generation processing, and projection obtained by the projection matrix generation processing. The projection vector determination process for determining the elements of the P projection vectors from the matrix and the residual candidate generation process are performed. Therefore, the P + 1 residual candidates generated and the quantized residual used in the next projection matrix generation process are determined from the elements in the Pth row of the projection matrix of P rows and P columns obtained by the projection matrix generation process. The residual determination process is performed Q times (Q ≧ 1).

作用 本発明は、前記構成により、残差Rを与えると、P個
の投影ベクトルの要素を同時に算出し、高速化を図るも
のである。
Operation The present invention is to increase the speed by simultaneously calculating the elements of the P projection vectors when the residual R is given by the above configuration.

実施例 以下、本発明の一実施例について説明する。Example One example of the present invention will be described below.

なお、説明を簡明にするために、縮小に於ける投影情
報を4ビット並列にもとめる際を、説明する。
In order to simplify the explanation, a description will be given of the case where the projection information in the reduction is set to 4 bits in parallel.

第7図において、x軸方向に原画素の並びA(a1
a2、a3、a4)を、y軸方向に結果画素の並びB(b1
b2、b3、b4)を、更に、各原画素の中心を貫く垂直格子
線gと、結果画素間の境界を通る水平格子線hを想定す
る。垂直格子線gのインデックスiは、giの貫く原画素
の値と同じであり、水平格子線hのインデックスjは、
hjのすぐ上の結果画素のインデックスと同一である。格
子線の間隔は、水平及び垂直の両方向において等しく、
その長さをlとする。傾きn/m(m≧n)、初期残差R0
をy切片とする直線Lと垂直交子線giとの交点をPiとす
る。点p1のy座標y1は、残差の初期値R0にn/m lを加算
した値である。またyiとyi+1の差はn/mであるから、 である。
In FIG. 7, an array of original pixels A (a 1 ,
a 2 , a 3 , a 4 ), and the resulting pixel array B (b 1 ,
b 2 , b 3 , b 4 ), and a vertical grid line g passing through the center of each original pixel and a horizontal grid line h passing through the boundary between the result pixels. The index i of the vertical grid line g is the same as the value of the original pixel penetrated by gi, and the index j of the horizontal grid line h is
It is the same as the index of the result pixel just above hj. The grid lines are equally spaced in both horizontal and vertical directions,
Let l be its length. Slope n / m (m ≧ n), initial residual R 0
Let Pi be the point of intersection of the straight line L with y as the y-intercept and the vertical cross line gi. Y-coordinate y 1 of the point p 1 is a value obtained by adding n / ml to the initial value R 0 of the residual. Also, since the difference between yi and yi + 1 is n / m, Is.

次に、aiの中心を貫く垂直格子線giと水平格子線hjと
の交点をQijとすると、その高さQijは、 Qij=j×l=hj −(4) となる。ただしhjは、hjの高さを表わす。
Next, when the intersection of the vertical grid line gi penetrating the center of ai and the horizontal grid line hj is Qij, the height Qij is Qij = j × l = hj- (4). However, hj represents the height of hj.

次に、以下のような要素Cijを持つ、行列C(以下、
投影行列と呼ぶ)を考える。
Next, the matrix C (hereinafter,
Call the projection matrix).

これは、前記第(3)、(4)の式より、 と書ける。即ち、要素Cijは、点piが格子点Qijの上にあ
るか、下にあるかを表わしている。つまり、要素Cijの
値が1ならば、点piは、水平格子線hj上に、もしくは上
位に位置しており、要素Cijの値が0ならば、点piは水
平格子線hjの下位に位置している。
From the above equations (3) and (4), Can be written. That is, the element Cij represents whether the point pi is above or below the grid point Qij. That is, if the value of the element Cij is 1, the point pi is located on or above the horizontal grid line hj, and if the value of the element Cij is 0, the point pi is located below the horizontal grid line hj. are doing.

次に、縮小の場合は、n/m1であるから、投影行列
Cは次の条件を満たす。
Next, in the case of reduction, since n / m1, the projection matrix C satisfies the following condition.

Cij=0;i<j 従って、投影行列Cは、 となる。Cij = 0; i <j Therefore, the projection matrix C is Becomes

第6図の投影行列C6である。この投影行列C6の第1行より点P1は、格子点Q
11の上位にあり、Q12の下位にあることがわかる。第2
行により点P2は、格子点Q21の上位にあり、かつQ22の下
位にあることがわかる。第3行により点P3は、格子点Q
32の上位にあり、かつQ33の下位にあることがわかる。
第4行により点P4は、格子点Q42の上位にあり、かつQ43
の下位にあることがわかる。
The projection matrix C 6 in FIG. 6 is Is. From the first row of this projection matrix C 6 , the point P 1 is the grid point Q
It can be seen that it is above 11 and below Q 12 . Second
The row shows that point P 2 is above grid point Q 21 and below grid point Q 22 . By the third line, the point P 3 is the grid point Q.
You can see that it is above 32 and below Q 33 .
By the fourth line, point P 4 is above grid point Q 42 , and Q 43
You can see that it is in the lower rank of.

このように、投影行列Cの第i行は、原画素aiの直線
L上への投影点piがどの格子点間に存在するかという情
報をもっている。この投影行列Cから、投影ベクトルT
を求めることが可能である。この方法を、以下に説明す
る。
As described above, the i-th row of the projection matrix C has information about which grid point the projection point pi of the original pixel ai on the straight line L exists. From this projection matrix C, the projection vector T
Is possible. This method will be described below.

投影ベクトルの第1番目の要素t1は、C11の値に等し
い。なぜなら、C11の値が1ならば、点P1は、水平格子
線h1の上位に位置している。よって、投影ベクトルの要
素t1の値は1である。C11の値が0ならば、点P1は、水
平格子線h1の下位に位置している。よって、投影ベクト
ルの要素t1の値は0である。
The first element t 1 of the projection vector is equal to the value of C 11 . Because, if the value of C 11 is 1, the point P 1 is located above the horizontal grid line h 1 . Therefore, the value of the element t 1 of the projection vector is 1. If the value of C 11 is 0, the point P 1 is located below the horizontal grid line h 1 . Therefore, the value of the element t 1 of the projection vector is 0.

次に、ti+1の求め方は、tiの値が0であるか1であ
るかによって異なる。
Next, the method of obtaining ti + 1 differs depending on whether the value of ti is 0 or 1.

(i) ti=1の場合(第9図に示す) 投影行列第i行中の1の値をもつ要素の中で、インデ
ックスjの最大のものをCi、jmaxとする。即ち、点pi
は、格子点Qi、jmaxの上位にあり、格子点Qi、jmax+1
の下位にある。よって、投影ベクトルの第i+1番目の
要素であるti+1の値は、Ci+1、jmax+1の値に等し
い。即ち、Ci+1、jmax+1は、点pi+1が格子点Qi+
1、jmax+1の上位にあるか否かを判断する要素であ
り、Ci+1、jmax+1の値が1の場合は、点pi+1と点
piの間に、水平格子線h jmax+1が存在していることを
意味しており、投影ベクトルの要素ti+1の値は1とな
る。Ci+1、jmax+1の値が0の場合、点pi+1は水平
格子線h jmax+1の下位にある。つまり、点pi+1と点
piの間には、水平格子線が存在しないということを、意
味しているから、投影ベクトルの要素ti+1の値は0と
なる。
(I) When ti = 1 (shown in FIG. 9) Of the elements having a value of 1 in the i-th row of the projection matrix, the maximum index j is defined as Ci and jmax. That is, the point pi
Is above the grid points Qi, jmax, and the grid points Qi, jmax + 1
Is in the lower part of. Therefore, the value of ti + 1 which is the i + 1th element of the projection vector is equal to the value of Ci + 1 and jmax + 1. That is, for Ci + 1 and jmax + 1, the point pi + 1 is the grid point Qi +
It is an element to judge whether it is higher than 1, jmax + 1. If the value of Ci + 1, jmax + 1 is 1, point pi + 1 and point
This means that the horizontal grid line h jmax + 1 exists between pi, and the value of the element ti + 1 of the projection vector is 1. When the values of Ci + 1 and jmax + 1 are 0, the point pi + 1 is below the horizontal grid line h jmax + 1. That is, point pi + 1 and point
Since it means that there is no horizontal grid line between pi, the value of the element ti + 1 of the projection vector is zero.

(ii) ti=0の場合(第10図に示す) 投影行列第i行中の0の値をもつ要素の中で、インデ
ックスjの最小値をjminとする。即ち、点piは、格子点
Qi、jminの下位にあり、格子点Qi、jmin−1の上位にあ
る。よって、投影ベクトルの第i+1番目の要素ti+1
は、Ci+1、jminの値に等しい。即ち、Ci+1、jmin
は、点pi+1が、格子点Qi+1、jminの上位にあるか否
かを判断する要素であり、Ci+1、jminの値が1の場合
は、点pi+1と点piの間に、水平格子線h jminが存在し
ていることを意味しており、投影ベクトルの要素ti+1
の値は1である。Ci+1、jminの値が0の場合、点pi+
1は、水平格子線h jminの下位に位置している。つま
り、点pi+1と点piの間には、水平格子線は存在しない
ということを意味しているから、投影ベクトルの要素ti
+1の値は、0である。
(Ii) When ti = 0 (shown in FIG. 10) Among the elements having a value of 0 in the i-th row of the projection matrix, the minimum value of the index j is jmin. That is, the point pi is the grid point
It is lower than Qi and jmin and higher than the lattice point Qi and jmin-1. Therefore, the i + 1th element ti + 1 of the projection vector
Is equal to the value of Ci + 1, jmin. That is, Ci + 1, jmin
Is an element for determining whether or not the point pi + 1 is above the lattice points Qi + 1, jmin. When the value of Ci + 1, jmin is 1, the horizontal lattice line h jmin is between the points pi + 1 and pi. Means that there exists, and the projection vector element ti + 1
Has a value of 1. If Ci + 1 and jmin are 0, point pi +
1 is located below the horizontal grid line h jmin. That is, it means that there is no horizontal grid line between the point pi + 1 and the point pi.
The value of +1 is 0.

以上のようにして、投影行列Cより、投影ベクトルT
を、一義的に定めることができる。
As described above, from the projection matrix C, the projection vector T
Can be uniquely determined.

さて、投影ベクトルTから、投影行列Cを求めること
も可能である。
Now, it is also possible to obtain the projection matrix C from the projection vector T.

例えば、4ビットの投影ベクトルT7(1010)となる場
合の投影行列C7の求め方を、第7図を用いて以下に説明
する。
For example, a method of obtaining the projection matrix C 7 in the case of a 4-bit projection vector T 7 (1010) will be described below with reference to FIG. 7.

まず第1の処理として、投影行列第1行の要素をもと
める。
First, as the first processing, the element of the first row of the projection matrix is obtained.

点P1は、格子点Q11より上位にあり、格子点Q12
Q13、Q14の下位にある。よって (C11、C12、C13、C14)=(1、0、0、0)−(8) 次に第2の処理として、投影行列第2行の要素をもと
める。
The point P 1 is higher than the lattice point Q 11 , and the lattice point Q 12 ,
Lower in the Q 13, Q 14. Therefore, (C 11 , C 12 , C 13 , C 14 ) = (1, 0, 0, 0) − (8) Next, as the second processing, the element of the second row of the projection matrix is obtained.

点P2は、格子点Q21より上位にあり、格子点Q22
Q23、Q24の下位にある。よって (C21、C22、C23、C24)=(1、0、0、0)−(9) そして第3の処理として、投影行列第3行の要素をも
とめる。
Point P 2 is higher than grid point Q 21 , and grid point Q 22 is
Lower in the Q 23, Q 24. Therefore, (C 21 , C 22 , C 23 , C 24 ) = (1, 0, 0, 0) − (9) Then, as the third processing, the element of the third row of the projection matrix is obtained.

点P3は、格子点Q31、Q32より上位にあり、格子点
Q33、Q34の下位にある。よって (C31、C32、C33、C34)=(1、1、0、0)−(10) そして第4の処理として、投影行列の第4行の要素を
もとめる。
Point P 3 is higher than grid points Q 31 and Q 32 , and
It is lower than Q 33 and Q 34 . Therefore, (C 31 , C 32 , C 33 , C 34 ) = (1, 1, 0, 0) − (10) Then, as the fourth processing, the element in the fourth row of the projection matrix is obtained.

点P4は、格子点Q41、Q42の上位にあり、格子点Q43、Q
44の下位にある。よって (C41、C42、C43、C44)=(1、1、0、0)−(11) 以上(8)、(9)、(10)、(11)より、推定行列
C6は、 となる。
Point P 4 is above grid points Q 41 and Q 42 , and grid point Q 43 and Q 42
It is under 44 . Therefore, (C 41 , C 42 , C 43 , C 44 ) = (1,1,0,0) − (11) From the above (8), (9), (10), (11), the estimation matrix
C 6 is Becomes

以上のように、投影ベクトルTは、投影行列Cから求
めることが可能である。また、逆に、投影行列Cは、投
影ベクトルTから求めることも可能である。即ち、投影
ベクトルTと、投影行列Cは、一対一の対応関係を持
つ。
As described above, the projection vector T can be obtained from the projection matrix C. On the contrary, the projection matrix C can be obtained from the projection vector T. That is, the projection vector T and the projection matrix C have a one-to-one correspondence.

次に、第7図の投影行列C6より、投影ベクトルT
6(t1、t2、t3、t4)を求める場合を説明する。
Next, from the projection matrix C 6 in FIG. 7, the projection vector T
A case of obtaining 6 (t 1 , t 2 , t 3 , t 4 ) will be described.

前述したように、t1の値は要素C11に等しい。よって t1=C11=1 −(13) である。As mentioned above, the value of t 1 is equal to element C 11 . Therefore, t 1 = C 11 = 1- (13).

t2の値は、t1の値が1より、前記(i)の場合を用い
る。jmaxは1であるから、t2の値は要素C2の値に等
しい。よって t2=C22=0 −(14) である。
As the value of t 2 , the case of (i) above is used since the value of t 1 is 1. Since jmax is 1, the value of t 2 is equal to the value of the element C 2, 2. Therefore, t 2 = C 22 = 0 − (14).

t3の値は、t2の値が0より、前記(ii)の場合を用い
る。jminは2であるから、t3の値は、要素C32の値に等
しい。よって t3=C32=1 −(15) である。
As the value of t 3 , the case of (ii) above is used since the value of t 2 is 0. Since jmin is 2, the value of t 3 is equal to the value of element C 32 . Therefore, t 3 = C 32 = 1- (15).

t4の値は、t3の値が1より、前記(i)の場合を用い
る。jmaxは2であるから、t4の値は、要素C43の値に等
しい。よって t4=C43=0 −(16) である。
As the value of t 4 , the case of (i) above is used since the value of t 3 is 1. Since jmax is 2, the value of t 4 is equal to the value of element C 43 . Therefore, t 4 = C 43 = 0 − (16).

以上、第(13)、(14)、(15)、(16)式より、投
影ベクトルT6は、 T6(t1、t2、t3、t4)=(1010) となる。
As described above, from the expressions (13), (14), (15), and (16), the projection vector T 6 is T 6 (t 1 , t 2 , t 3 , t 4 ) = (1010).

ところで、投影ベクトルT6を求める際に用いた要素
は、C11、C22、C32、C43である。その他の要素C12
C13、C14、C23、C34は常に0であり、算出する必要はな
い。更に、投影ベクトルT6を求める際に用いた投影行列
の要素がC11、C22、C32、C43であるように投影ベクトル
Tの要素を求める際に必要な投影行列の要素は、前記
(i)、(ii)よりjmaxまたはjminのインデックスをも
つもののみで、その他は関与しない。よって他の要素は
算出する必要はない。即ち、投影ベクトルTを求める
際、投影行列Cのいくつかの要素は不要となる。このこ
とは、回路の簡略化に有効である。
By the way, the elements used when obtaining the projection vector T 6 are C 11 , C 22 , C 32 , and C 43 . Other elements C 12 ,
C 13 , C 14 , C 23 , and C 34 are always 0, and need not be calculated. Furthermore, the elements of the projection matrix used when obtaining the projection vector T 6 are C 11 , C 22 , C 32 , and C 43 so that the elements of the projection matrix required when obtaining the elements of the projection vector are From (i) and (ii), only those having an index of jmax or jmin are not involved. Therefore, it is not necessary to calculate other elements. That is, when obtaining the projection vector T, some elements of the projection matrix C are unnecessary. This is effective in simplifying the circuit.

以上の説明に於ては、投影ベクトルの次元は4であっ
たが、上記の方法を一般的なlの次元を持った場合に拡
張すると次のようになる。
In the above description, the dimension of the projection vector is 4, but the above method is expanded to have a general dimension of l as follows.

投影ベクトルT(t1、t2、……、tl)の第1番目の要
素t1は、投影行列の要素C11に等しいことは、すべてに
述べた。ここで、要素C11の値が、1の場合、第1列中
の要素C21、C31、……、Cl1の値は、必ず1となる。な
ぜなら、必ず点P2、P3、……、Plは、点P1の上位に存在
するからである。Cijが1ならば、第j列中の要素Ci+
1、j、Ci+2、j、……Cl、jの値は、上記と同様に
必ず1となる。
It has been mentioned above that the first element t 1 of the projection vector T (t 1 , t 2 , ..., T l ) is equal to the projection matrix element C 11 . Here, when the value of the element C 11 is 1, the value of the elements C 21 , C 31 , ..., C 11 in the first column is always 1. This is because the points P 2 , P 3 , ..., P 1 always exist above the point P 1 . If Cij is 1, element Ci + in the j-th column
The values of 1, j, Ci + 2, j, ... C l , j are always 1 as in the above case.

また、投影する直線の傾きn/mは、1に等しいか、1
より小さい場合であるから、投影行列の行番号が1増え
る度に、その行中の1の値をもつ要素の数は、たかだが
1増加するだけである。つまり、投影ベクトルの第i番
目の要素tiを決定した、投影行列の要素をCijとする
と、以下のことが、示される。
Also, the slope n / m of the projected straight line is equal to 1 or 1
Since it is smaller, each time the row number of the projection matrix increases by 1, the number of elements having a value of 1 in the row increases by 1 at most. That is, when the element of the projection matrix that determines the i-th element ti of the projection vector is Cij, the following is shown.

ti=Cij=1のとき、ti+1=Ci+1、j+1 −(17) ti=Ci、j=0のとき、ti+1=Ci+1、j −(18) これより、第8図の投影ベクトルT(t1、t2、……、
tl)を求める。t1は、C11の値に等しい。よってt1は、
1である。t2は、第(17)式より、C22の値に等しい。
よってt2は、0である。t3は、第(18)式より、C32
値に等しい。よってt3は1である。t4は第(17)式より
C43の値に等しい。よって、t4は、1である。以上の処
理を繰り返す。
When ti = Cij = 1, ti + 1 = Ci + 1, j + 1- (17) ti = Ci, j = 0, ti + 1 = Ci + 1, j- (18) From this, the projection vector T (t 1 , t 1 in FIG. t 2 , ……,
t l ). t 1 is equal to the value of C 11 . Therefore t 1 is
It is 1. From the equation (17), t 2 is equal to the value of C 22 .
Therefore, t 2 is 0. From Equation (18), t 3 is equal to the value of C 32 . Therefore t 3 is 1. t 4 is calculated from equation (17)
Equal to the value of C 43 . Therefore, t 4 is 1. The above process is repeated.

tl-1は、tl-2の値を決定した投影行列の要素がC
l−2,k−2で、その値が0の場合、第(18)式よりC
l−1,k−2の値に等しく、ll-2の値を決定した投影行
列の要素がCl−2,k−2で、その値が1の場合、第(1
7)式よりCl−1,k−1の値に等しい。
The element of the projection matrix that determines the value of t l-2 is t l-1.
If it is l-2, k-2 and its value is 0, then from the equation (18), C
If the element of the projection matrix that determines the value of l 1-2 is equal to the value of 1−1 , k−2 and is C 1−2, k−2 , and the value is 1, then the first (1
From the equation (7) , it is equal to the value of C l−1, k−1 .

同様に、tの値は、tl-1=Cl−1,k−1=0の場合
は、Cl,k−1の値に等しく、t1-1=Cl−1,k−1=1
の場合は、Cl,kの値に等しい。
Similarly, the value of t is equal to the value of C l, k−1 when t l−1 = C l−1, k−1 = 0, and t 1-1 = C l−1, k− 1 = 1
, Then equal to the value of C l, k .

次にm個の要素をもつ投影ベクトルT(t1、t2、…
…、tm)を求める場合には、 (イ) 以上の方式をm次元に拡大する (ロ) l個づつ逐次的に求める の2つの方法がある。(イ)は、l=mとするだけであ
るから、以下(ロ)の場合について述べる。説明を簡単
にするため l=4 4(K−1)<m≦4K の場合を説明する。
Next, a projection vector T (t 1 , t 2 , ...
, T m ), there are two methods of (a) expanding the above method into m dimensions and (b) sequentially obtaining l pieces each. In (a), since only 1 = m is set, the case of (b) will be described below. In order to simplify the description, the case of l = 44 (K-1) <m ≦ 4K will be described.

第7図に於いて、投影行列の要素C11、C21、C22
C31、C32、C33、C41、C42、C44を算出し、投影ベクトル
T(t1、t2、t3、t4)を求める。次の投影ベクトルT
(t5、t6、t7、t8)を求める際に必要な、投影ベクトル
Cの要素を算出する場合、初期残差R0と同等の意味をも
つ、残差R4が必要となる。
In FIG. 7, elements C 11 , C 21 , C 22 of the projection matrix,
C 31 , C 32 , C 33 , C 41 , C 42 , and C 44 are calculated, and the projection vector T (t 1 , t 2 , t 3 , t 4 ) is obtained. Next projection vector T
When calculating the elements of the projection vector C necessary for obtaining (t 5 , t 6 , t 7 , t 8 ), a residual R 4 having the same meaning as the initial residual R 0 is required. .

残差R4の値は、以下の第(19)式により求まる5つの
値、 のいづれかである。つまり、次の4個を求める際に必要
な、残差の初期値は、前段の最終段差である。この次の
残差を決めるのが、C41、C42C43、C44である。これは、
次のような対応関係(第(20)式)がある。
The value of the residual R 4 is five values obtained by the following equation (19), It is either That is, the initial value of the residual required to obtain the next four is the final step of the preceding stage. It is C 41 , C 42 C 43 , and C 44 that determine the next residual. this is,
There is the following correspondence (formula (20)).

C41=0のときR4=R0+4n C41=1かつC42=0のときR4=R0+4n−m C42=1かつC43=0のときR4=R0+4n−2m C43=1かつC44=0のときR4=R0+4n−3m C44=1のときR4=R0+4n−4m 以上のように、残差をもとめ、次のひとまとまりの投
影行列Cの要素を算出し、4ビットの投影ベクトルTを
求める。
When C 41 = 0, R 4 = R 0 + 4n When C 41 = 1 and C 42 = 0, R 4 = R 0 + 4n−m When C 42 = 1 and C 43 = 0, R 4 = R 0 + 4n−2m When C 43 = 1 and C 44 = 0, R 4 = R 0 + 4n−3m When C 44 = 1 R 4 = R 0 + 4n−4m As shown above, the residual is calculated and the next set of projection matrix is obtained. The element of C is calculated to obtain the 4-bit projection vector T.

mビットの投影ベクトルTをもとめる場合、第6図の
ように、mビットを4ビット毎のブロックに分割し、上
記の処理をm個の要素をもつ投影ベクトルが決定される
まで、すなわちk回繰り返すことになる。
When obtaining the m-bit projection vector T, as shown in FIG. 6, the m-bit is divided into blocks each having 4 bits, and the above process is repeated k times until a projection vector having m elements is determined. I will repeat.

以上の方法は、lビットを同時にもとめる場合に対し
ても容易に拡張可能である。
The above method can be easily extended to the case of obtaining 1 bit at the same time.

また、拡大の場合は、第7図において、x軸方向に結
果画素の並びを、y軸方向に原画素の並びを想定すれば
よく、投影ベクトルの算出方法と、投影ベクトルは変わ
らない。つまり、拡大の場合は、結果画素mビットを原
画素nビット(m≧n)に対応させる場合と考えればよ
い。
Further, in the case of enlargement, in FIG. 7, the arrangement of the result pixels in the x-axis direction and the arrangement of the original pixels in the y-axis direction may be assumed, and the projection vector calculation method and the projection vector do not change. That is, in the case of enlargement, it can be considered that the result pixel m bits correspond to the original pixel n bits (m ≧ n).

次に、初期残差R0、原画素をmビット,結果画素をn
ビットとし、原画素を結果画素に縮小する場合のmビッ
トからなる投影ベクトルの要素を、l個同時に算出する
場合のフローを第5図に示す。
Next, the initial residual R 0 , the original pixel is m bits, and the result pixel is n.
FIG. 5 shows a flow for simultaneously calculating 1 element of a projection vector consisting of m bits in the case where the original pixel is reduced to the result pixel in bits.

41、42は原画素カウンタおよび残差レジスタに初期値
として、0、R0を設定する。43は原画素数m,初期残差
R0、結果画素数nから第(5)式に従い、l行l列の投
影行列を求める投影行列生成処理と、第(19)式に従い
l+1種の残差を求める残差候補生成処理を行う。
41 and 42 set 0 and R 0 as initial values in the original pixel counter and the residual register. 43 is the number of original pixels m, initial residual
A projection matrix generation process for obtaining a projection matrix of l rows and 1 columns from R 0 and the result pixel number n according to the formula (5) and a residual candidate generation process for obtaining l + 1 types of residuals according to the formula (19) .

44は処理43で算出したl行l列の投影行列から第(1
7)式、第(18)式に基づきl個の投影ベクトルの要素
を決定する投影ベクトル決定処理である。
44 is the first (1
This is a projection vector determination process for determining the elements of l projection vectors based on the equations (7) and (18).

45は次のl個投影ベクトルの要素を計算するのに必要
な残差を、43の残差候補生成処理で求めた残差候補か
ら、同じく43の投影行列生成処理で求めた投影行列の要
素を用いて第(20)式に従い決定する残差決定処理であ
る。46は原画素カウンタと投影情報ポインタをlだけ増
化する処理である。47は、投影ベクトルのすべての要素
の算出が終了したかどうかの判断を行う処理である。以
上のようにしてmビットの投影ベクトルをもとめる。
45 denotes the residuals required to calculate the elements of the next l projection vectors, the elements of the projection matrix also obtained by the projection matrix generation processing of 43 from the residual candidates obtained by the residual candidate generation processing of 43. Is a residual determination process that is performed according to the equation (20) using. 46 is a process for increasing the original pixel counter and the projection information pointer by 1. 47 is a process of determining whether or not the calculation of all the elements of the projection vector has been completed. As described above, the m-bit projection vector is obtained.

次に、以上の原理を用いた実際の装置を説明する。第
1図は4次元の投影ベクトルを生成する演算装置の一実
施例を示すモジュール結線図である。001、002、003、0
04はそれぞれ2n、4n、2m、4mを出力する桁上演算器(以
下、シフターと呼ぶ)である。300、302はそれぞれ3n、
3mを計算する加算器である。501、502、503、504はそれ
ぞれn、m、R、及び第6の処理で得られた次の残差用
の信号線である。701は第6の処理で得られた次の残差
用のレジスタである。
Next, an actual device using the above principle will be described. FIG. 1 is a module connection diagram showing an embodiment of an arithmetic unit for generating a four-dimensional projection vector. 001, 002, 003, 0
04 is a carry calculator (hereinafter referred to as a shifter) that outputs 2n, 4n, 2m, and 4m, respectively. 300 and 302 are 3n,
It is an adder that calculates 3m. 501, 502, 503, and 504 are signal lines for n, m, R, and the next residual obtained in the sixth processing, respectively. Reference numeral 701 is a register for the next residual obtained in the sixth processing.

上記構成において、前述したように、拡大でも縮小で
も投影ベクトルは変わらないので、説明を簡単にするた
め、以下、倍率n/mの縮小(m≧n)の場合について、
この動作を説明する。まず残差Rと倍率n/mが与えられ
ると、501、502、503の信号線を通って、第1から第4
演算ユニットI〜IVにそれぞれn、m、Rが入力され
る。
In the above configuration, as described above, since the projection vector does not change between enlargement and reduction, in order to simplify the description, the following will be described in the case of reduction of magnification n / m (m ≧ n).
This operation will be described. First, when the residual R and the magnification n / m are given, the first to fourth
N, m, and R are input to the arithmetic units I to IV, respectively.

まず第1演算ユニットIでは、加算器101でR+nが
加算され、その出力とmとの大小比較を比較器201が行
い、投影行列の要素C11を出力する。
First, in the first arithmetic unit I, R + n is added by the adder 101, the output thereof is compared with m by the comparator 201, and the element C 11 of the projection matrix is output.

第2演算ユニットIIでは、シフター001、003より、2
n、2mがそれぞれ算出され、加算器201で、R+2nが計算
される。比較器211、212によって、mと2n+R、及び2m
と2n+Rの大小が比較され、投影行列の要素C21、C22
算出される。
In the second operation unit II, 2 from shifter 001, 003
n and 2m are calculated, and the adder 201 calculates R + 2n. M, 2n + R, and 2m by the comparators 211 and 212
And the magnitude of 2n + R are compared, and the elements C 21 and C 22 of the projection matrix are calculated.

第3演算ユニットIIIでは、加算器300、302で3n、及
び3mが、出力され、加算器301でR+3nが計算され、比
較器311、312、313で、mとR+3n、2mとR+3n、及び3
mとR+3nの大小比較が行なわれる。その結果投影行列
の要素C31、C32、C33が、算出される。
In the third arithmetic unit III, 3n and 3m are output by the adders 300 and 302, R + 3n is calculated by the adder 301, and m and R + 3n, 2m and R + 3n and 3 are calculated by the comparators 311, 312 and 313.
The magnitude comparison of m and R + 3n is performed. As a result, the elements C 31 , C 32 , and C 33 of the projection matrix are calculated.

第4演算ユニットIVでは、シフター002、004より、4
n、4mが、それぞれ算出され、加算器401でR+4nが算出
される。減算器411、412、413、414で、R+4n−m、R
+4n−2m、R+4n−3m及び、R+4n−4mが算出される。
これらは、次のフェイズに於ける残差となる可能性のあ
る値であり、この計算を行なうと、投影行列の要素
C41、C42、C43、C44も、同時に算出される。
In the 4th arithmetic unit IV, from shifter 002,004, 4
n and 4m are calculated respectively, and the adder 401 calculates R + 4n. Subtractors 411, 412, 413, 414 use R + 4n-m, R
+ 4n-2m, R + 4n-3m, and R + 4n-4m are calculated.
These are the values that may be the residuals in the next phase, and when this calculation is performed, the elements of the projection matrix are
C 41 , C 42 , C 43 and C 44 are also calculated at the same time.

投影ベクトル選択ユニット601の内容を、第2図に示
す。第1から第4演算処理でもとめた投影行列の要素を
入力すると、それに対応した4ビットの投影ベクトルが
出力されるマトリックス回路となっている。この回路
は、前述した推定行列Cの性質から簡略化してある。即
ち、投影ベクトルの要素を求める際に必要な投影行列の
要素は、jmax、またはjminのインデックスを持つものだ
けであるという性質と、ある行中の1の値を持つ要素数
は、行番号が1つ増える度に、たかだか1増加するだけ
という性質である。このことにより、投影ベクトルを4
ビット求めるために必要な投影行列の要素は、たかだか
4個でよいということになる。
The contents of the projection vector selection unit 601 are shown in FIG. When the elements of the projection matrix determined in the first to fourth arithmetic processes are input, the matrix circuit outputs a 4-bit projection vector corresponding thereto. This circuit is simplified from the property of the estimation matrix C described above. That is, the property that the elements of the projection matrix required to obtain the elements of the projection vector are only those having an index of jmax or jmin, and the number of elements having a value of 1 in a certain row is It is a property that only increments by 1 for each increment of 1. This makes the projection vector 4
This means that the number of elements of the projection matrix required to obtain bits is at most four.

次に残差選択ユニット603の内容を第3図に示す。第
4演算ユニットIVでもとめた投影行列の要素C41、C42
C43、C44の内容によりR+4n、R+4n−m、R+4n−2
m、R+4n−3m、R+4n−4mのうちの一つが、次に使用
される残差として出力されるようになっている。この回
路も、前記の性質を利用し、簡略化してある。
Next, the content of the residual selection unit 603 is shown in FIG. Elements C 41 , C 42 of the projection matrix determined by the fourth arithmetic unit IV,
C 43, the content of C 44 R + 4n, R + 4n-m, R + 4n-2
One of m, R + 4n-3m, and R + 4n-4m is output as the residual to be used next. This circuit is also simplified by utilizing the above property.

第1図に於いて、入力n、m、Rから4ビットの投影
ベクトルが出力されるまでの経路は、桁上演算器001を
通り、加算器300、301を通り、比較器311を通り、第6
の処理を行う601、602を通過するものが最長であるか
ら、高速に投影ベクトルをもとめることができる。
In FIG. 1, the path from the inputs n, m, and R to the output of the 4-bit projection vector passes through the carry arithmetic unit 001, the adders 300 and 301, and the comparator 311. Sixth
Since the longest one passes through 601 and 602, the projection vector can be obtained at high speed.

第4図に、第1図に対応したl次元の投影ベクトルを
もとめる場合のモジュール図を示す。このように、l個
同時に求める場合に、容易に拡張可能である。
FIG. 4 shows a module diagram for obtaining an l-dimensional projection vector corresponding to FIG. In this way, it is possible to easily expand when 1 is obtained at the same time.

発明の効果 本発明は、以上のように同時に複数個の要素を持つ、
画像の拡大縮小を制御する投影ベクトルを求めることに
より、従来要していた時間を大幅に減少することができ
る。
As described above, the present invention has a plurality of elements at the same time,
By obtaining the projection vector for controlling the scaling of the image, the time required conventionally can be greatly reduced.

【図面の簡単な説明】[Brief description of drawings]

第1図は本発明の一実施例における投影情報並列生成方
法に対するその演算装置のモジュール結線図、第2図、
第3図は第1図の要部ブロック結線図、第4図は第1図
構成をlビットに拡張した際のモジュール結線図、第5
図は本発明の一実施例における投影情報並列生成方法の
流れを説明するフロー図、第6図は残差の流れを表わす
概念図、第7図は本発明の4次元の投影ベクトルを生成
するための考え方を示した概念図、第8図は第7図をl
ビットに拡張した際の概念図、第9図、第10図は、それ
ぞれjmax、jminを用いた投影ベクトルの要素の求め方を
示した概念図、第11図は従来の投影ベクトルを生成する
ための基本的な考え方を示した概念図、第12図は従来の
縮小率3/7の場合の例を示した概念図、第13図は従来の
方法の流れを示したフロー図である。 001、002、003……桁上演算器、101、201、301、401…
…加算器、201、211、212、311、312、313……比較器、
411、412、413、414……減算器、601……投影ベクトル
選択ユニット、603……残差選択ユニット、701……残差
レジスタ。
FIG. 1 is a module connection diagram of an arithmetic unit for a projection information parallel generation method according to an embodiment of the present invention, FIG.
FIG. 3 is a block diagram of a main part of FIG. 1, FIG. 4 is a diagram of a module diagram when the configuration of FIG. 1 is expanded to 1 bit, and FIG.
FIG. 6 is a flow chart for explaining the flow of the projection information parallel generation method in one embodiment of the present invention, FIG. 6 is a conceptual diagram showing the flow of residuals, and FIG. 7 is for generating a four-dimensional projection vector of the present invention. 8 is a conceptual diagram showing the concept for
Conceptual diagrams when expanded to bits, FIGS. 9 and 10 are conceptual diagrams showing how to obtain the elements of the projection vector using jmax and jmin, respectively, and FIG. 11 is for generating conventional projection vectors. FIG. 12 is a conceptual diagram showing a basic concept of FIG. 12, FIG. 12 is a conceptual diagram showing an example of a conventional reduction ratio of 3/7, and FIG. 13 is a flow chart showing a flow of a conventional method. 001, 002, 003 ... Carry calculator, 101, 201, 301, 401 ...
… Adder, 201, 211, 212, 311, 312, 313 …… Comparator,
411, 412, 413, 414 ... subtractor, 601 ... projection vector selection unit, 603 ... residual selection unit, 701 ... residual register.

───────────────────────────────────────────────────── フロントページの続き (72)発明者 川上 桂 神奈川県川崎市多摩区東三田3丁目10番1 号 松下技研株式会社内 (56)参考文献 特開 昭60−189579(JP,A) ─────────────────────────────────────────────────── ─── Continuation of the front page (72) Inventor Katsura Kawakami 3-10-10 Higashisanda, Tama-ku, Kawasaki City, Kanagawa Prefecture Matsushita Giken Co., Ltd. (56) Reference JP-A-60-189579 (JP, A)

Claims (1)

【特許請求の範囲】[Claims] 【請求項1】J個の画素からなる第1の画素列と、K個
の画素からなる第2の画素列との直線への写影を介した
対応関係を表すLビット(LはJ,Kのうち大きな数)か
らなる投影ベクトルを、 第1の画素列と第2の画素列との縮小率n/m<1を決定
するn,m、及び与えらえた量子化残差から、P行P列
(L>P)の投影行列の要素を並列に求める投影行列生
成処理と、 前記投影行列生成処理に続く、次の投影行列生成処理に
用いる量子化残差のP+1個の候補をn,m、及び前記投
影行列生成処理で用いられた量子化残差から並列に求め
る残差候補生成処理と、 前記投影行列生成処理によって得られた投影行列からP
個の投影ベクトルの要素を並列に決定する投影ベクトル
決定処理と、 前記残差候補生成処理によって生成されたP+1個の残
差の候補と、前記投影行列生成処理により得られたP行
P列の投影行列の第P行の要素から次の投影行列生成処
理に用いる量子化残差を決定する残差決定処理とからな
る、一度にPビット(L>P)づつ並列に求める処理
を、Q回(Q≧1)行うことにより求めることを特徴と
する投影情報並列生成方法。
1. A L-bit (L is J, which represents a correspondence relationship between a first pixel column consisting of J pixels and a second pixel column consisting of K pixels through a mapping to a straight line. A projection vector consisting of a large number of K) is set to P from the n, m that determines the reduction ratio n / m <1 between the first pixel row and the second pixel row, and the given quantization residual. A projection matrix generation process for obtaining the elements of the projection matrix of row P columns (L> P) in parallel, and P + 1 candidates of quantized residuals used in the next projection matrix generation process subsequent to the projection matrix generation process are n , m, and a residual candidate generation process that is obtained in parallel from the quantized residuals used in the projection matrix generation process, and P from the projection matrix obtained by the projection matrix generation process.
Projection vector determination processing for determining the elements of the projection vectors in parallel, P + 1 residual candidates generated by the residual candidate generation processing, and P rows and P columns obtained by the projection matrix generation processing. A process for determining in parallel P bits (L> P) at a time, which consists of a residual determination process for determining a quantized residual used in the next projection matrix generation process from the elements in the Pth row of the projection matrix, is performed Q times. (Q ≧ 1) A projection information parallel generation method characterized by being obtained by performing.
JP61171211A 1986-07-21 1986-07-21 Parallel projection information generation method Expired - Lifetime JPH0827835B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP61171211A JPH0827835B2 (en) 1986-07-21 1986-07-21 Parallel projection information generation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP61171211A JPH0827835B2 (en) 1986-07-21 1986-07-21 Parallel projection information generation method

Publications (2)

Publication Number Publication Date
JPS6326781A JPS6326781A (en) 1988-02-04
JPH0827835B2 true JPH0827835B2 (en) 1996-03-21

Family

ID=15919094

Family Applications (1)

Application Number Title Priority Date Filing Date
JP61171211A Expired - Lifetime JPH0827835B2 (en) 1986-07-21 1986-07-21 Parallel projection information generation method

Country Status (1)

Country Link
JP (1) JPH0827835B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7848581B2 (en) * 2007-03-07 2010-12-07 Taiwan Semiconductor Manufacturing Company, Ltd. System and method for decoding and viewing of image files

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS60189579A (en) * 1984-03-09 1985-09-27 Matsushita Electric Ind Co Ltd Processor of picture data

Also Published As

Publication number Publication date
JPS6326781A (en) 1988-02-04

Similar Documents

Publication Publication Date Title
Monteiro et al. Interior path following primal-dual algorithms. Part I: Linear programming
Dyer et al. Parallel image processing by memory-augmented cellular automata
JP2022003501A (en) Quantum circuit simulation method, device, facility, storage medium, and program
JP3637923B2 (en) Method for operating a processing device
JPH03127282A (en) Radon transforming method of digital image
Hekstra et al. Floating point CORDIC
EP0744054A4 (en) High speed function generating apparatus and method
US5402533A (en) Method and apparatus for approximating a signed value between two endpoint values in a three-dimensional image rendering device
JPH0827835B2 (en) Parallel projection information generation method
JPH06266747A (en) Dct/inverse dct arithmetic unit
Li et al. Batched trajectory compression algorithm based on hierarchical grid coordinates
RU2698413C1 (en) Device for comparing numbers in a system of residual classes
US4484298A (en) Method and device for generation of quadratic curve signal
JP3231811B2 (en) Matrix operation circuit
Schrack et al. The spatial U-order and some of its mathematical characteristics
JP2713400B2 (en) Projection information generation device
JPH0314128A (en) Arithmetic unit utilizing logarithmic expression numerical value
CN113222154A (en) Method and device for determining amplitude of quantum state
JP2605792B2 (en) Arithmetic processing unit
US3324288A (en) Data processing apparatus including means for correcting codes arranged in a packed format
SU697994A1 (en) Device for computing elementary functions
JP2889244B2 (en) Image processing device
Nielsen et al. Algorithms On Continued And Multicontinued Fractions
CN112799845A (en) Graph algorithm parallel acceleration method and device based on GRAPE framework
SU1305671A1 (en) Device for calculating values of function z=arccos y/x