Method and assembly for recording a radiation image of an elongate body

[0002] The application claims the benefit of the U.S. Provisional application No. 60/071,630 filed Jan. 16, 1998.

[0003] The invention relates to a method and an assembly for recording a radiation image of an elongate body such as a full leg or a full spine.

[0004] In flat projection radiography examination types exist that aim at imaging a larger portion of the body than can be fit on a single, even the largest imaging cassette.

[0005] Such imaging needs arise in so-called full-leg or full-spine examinations, where clinical indication requires that e.g. the full spine or the full leg is examined at once so as to enable or to quantify diagnosis.

[0006] Analogous to conventional film-screen based radiography, in storage phosphor image acquisition wherein a radiation image is temporarily stored in a photostimulable phosphor screen, such examinations are achieved by filling a larger cassette with a plurality of storage phosphor screens, generally overlapping with one another so as to completely cover the longer, elongate cassette (an embodiment referred to as an overlapping sheet embodiment).

[0007] Such elongate cassettes are typically 35×105 cm for full-leg examinations and 30×90 cm for full-spine examinations. In the overlapping sheet embodiment, the largest single film cassette measures 35×43 cm and hence is unsuitable for imaging the spinal column or a leg. Typically, 4 screens are used, e.g. 4 (24×30 cm) screens are arranged in a 30×90 cm cassette, clearly resulting in an overlap configuration. Therefore, the part of the image corresponding to the overlap zone, will be less exposed.

[0008] Alternatively, the situation of non-overlapping images may also occur, e.g. in a configuration of 3 (35×35 cm) imaging screens in a 35×105 cm cassette, resulting in 3 pair-wise touching but otherwise non-overlapping images.

[0009] During exposure, a fixed grid is simultaneously present in the path of the X-ray beam, resulting in an image of a raster of horizontal and vertical parallel thin lines superposed to the radiation image of the elongate body. These lines aid the radiologist or operator in reconstructing the original geometry of the body, since lines on image must necessarily extend continuously and seamlessly into lines of the previous image and the next image.

[0010] In U.S. Ser. No. 09/035,528, a method is disclosed to acquire a radiation image of an elongate body by using a sequence of recording members in partially overlapping disposition.

[0011] This patent application further relates to a method to obtain a single ‘stitched’ image of an elongate body using the images read out of each of these recording members.

[0012] The image acquisition process assumes the use of an elongate cassette in which a sequence of stimulable phosphor screens are arranged in an overlapping manner so as to completely cover the cassette.

[0013] After exposure, the elongate cassette is opened, the individual screens are removed from the cassette and are each put into a smaller cassette of a size corresponding to the dimensions of an individual screen. Then, the cassettes are sequentially fed into a dedicated read-out apparatus where the cassette is again opened, the screen is taken out of the cassette, and scanned by means of stimulating radiation. The image-wise modulated light that is emitted upon stimulation is detected and converted into a digital image representation of part of the elongate body. Next, the digital image representations read out of each of the cassettes are recombined to form an image of the entire elongate body.

[0014] The handling process has several drawbacks.

[0015] First, a dark room is needed to transfer the individual screens from the elongate cassette to their corresponding individual regular size cassettes, since ambient light impinging on a screen may stimulate it, and hence may (at least partially) erase its contents.

[0016] Second, the operator handling the screens may accidentally damage them, e.g. by scratching them.

[0017] Third, the overall time-to-diagnosis is lengthened due to the additional operations and potentially increased risk of re-take.

[0018] U.S. patent application Ser. No. 09/035,528 further discloses a method of reconstructing the image of the elongate body from the partial images read out from each of the overlapping recording members, each having stored therein a part of the elongate scene. The distortions coped with in aforementioned patent application were limited to shift, overlap and rotation of each sub-image with respect to the previous and next sub-image. The method disclosed in above patent application depends on two basic algorithms: (1) detection and modelling of the periodic grid and (2) cross-correlation of overlapping image parts. The periodic grid is important to be able to align and stitch the images with respect to each other. The grid is located physically on the elongate cassette and covers its full surface in a continuous way. Hence, the grid should also be present in the resulting re-constructed image composition in a continuous way.

[0019] It is an object of the present invention to provide a method of recording a radiation image of an elongate body that does not show the drawbacks of the prior art.

[0020] It is another object of this invention to provide an assembly for recording a radiation image of an elongate body on a plurality of photostimulable phosphor screens each having a length that is smaller than the length of the elongate body, which assembly overcomes the shortcomings of the prior art.

[0021] Further objects of this invention will become apparent from the description hereafter.

[0022] The above-mentioned objects are realised by a method as disclosed in claim 1.

[0023] The prior art shortcomings are overcome by using an assembly of regular sized cassettes instead of an assembly of screens in a large size cassette.

[0024] In one embodiment the radiation sensitive recording member is photostimulable phosphor screen. Other embodiments such as a film-intensifying screen combination may be envisaged.

[0025] Several arrangements of cassettes are possible and will be described below.

[0026] In order to be able to exactly reconstruct the entire image of the elongate body, a pattern of reference marks is likewise recorded onto the radiation sensitive (radiation attenuating) members in each of the individual cassettes.

[0027] Several kinds of reference marks are possible. However, a grid made of radiation sensitive material which covers the entire length of the staggered cassette arrangement is preferred.

[0028] During exposure, the grid is present in the path of the X-ray beam, resulting in an image of a raster of horizontal and vertical parallel thin lines superposed to the radiation image of the elongate body. Lines on image must necessarily extend continuously and seamlessly into lines of the previous image and the next image. The image of the grid thus guides the radiologist in reconstructing the image and as will be explained further on.

[0029] Problems that are associated with the fact that the cassettes in the assembly of the present invention have a thickness which cannot be neglected, are also solved by means of the image of such a grid.

[0030] The present invention overcomes the prior art disadvantages since the individual screens need no longer to be handled in dark room conditions and the risk of accidental damaging of the screens or dust contamination is eliminated since the screens are no longer taken out of the cassettes.

[0031] Further advantages and embodiments of the present invention will become apparent from the following detailed description illustrated by means of the following drawings wherein

[0032] FIG. 1 depicts a mechanical set up for the acquisition of an image of an elongate body using a plurality of regular size cassettes in an alternating overlapping disposition,

[0033] FIG. 2 depicts a mechanical set-up for the acquisition of an image of an elongate body using a stack of regular size cassettes in a staircase overlapping disposition,

[0034] FIG. 3 depicts a mechanical set-up for the acquisition of an image of an elongate body using a stack of regular size cassettes in an oblique overlapping disposition,

[0035] FIG. 4 depicts a grid imprint (solid lines) in one of the acquired sub-images obtained by the method of oblique overlapping disposition. The grid may further be rotated with respect to the image borders either due to rotation of the screen in the regular size cassette, or due to rotation of the screen with respect to the scanning axis in the read-out apparatus or due to a combination of said misalignments. The dashed lines ( - - - ) denote subdivision into rectangular vertical and horizontal blocks for profile calculation, the dotted lines ( . . . ) denote the centre lines of the corresponding blocks. The hole marks (o) denote the positions of the vertical grid lines in the horizontal profiles, the cross marks (x) denote the positions of the horizontal grid lines in the vertical profiles,

[0036] FIG. 5 is a block diagram illustrating the different steps of the algorithm.

[0037] While the present invention will hereinafter be described in connection with preferred embodiments thereof, it will be understood that it is not intended to limit the invention to those embodiments.

[0038] The present invention is illustrated by means of an embodiment wherein an X-ray image of e.g. a full spine of a patient is recorded on a number of photostimulable phosphor screens. On each of these screens another part of the image of the spine is recorded.

[0039] Each of the screens are conveyed in a regular size cassette having dimensions which correspond to the dimensions of the individual screen.

[0040] The screens are supported in a staggered disposition, in this example a partially overlapping staggered disposition, so that the elongate dimension of the assembly of the cassettes is at least as long as the elongate dimension (the length) of the elongate body.

[0041] In this way an X-ray image of the elongate body can be recorded in a single exposure onto the recording members in the assembly of cassettes.

[0042] FIGS. 1, 2, 3 depict different configurations of an assembly of a plurality of Z cassettes, each cassette holding a recording member. These figures also show a mechanical set up for holding such an assembly.

[0043] An assembly comprises a rigid support 3. Several cassettes 4 can be attached to or inserted in support 3, such that each cassette partially overlaps with its predecessor and/or successor cassette.

[0044] An assembly may comprise a front panel 6 at the tube side of the support. Such a front panel may hold a density attenuating (rectangular) grid and/or an anti-scatter grid, both indicated schematically by vertical line 7.

[0045] The front side (side directed to the X-ray tube which is schematically indicated by numeral 5) of the front panel 6 is made of X-ray transparent material so as not to absorb any radiation.

[0046] A standard wall bucky X-ray recording apparatus 1 used for recording a radiation image of a patient is shown in these figures.

[0047] The rigid support 3 is mounted in the vertical position to the bucky for vertical patient postures. The support can be hooked into the wall bucky X-ray recording apparatus by hook 8.

[0048] Alternatively, the support holding the cassettes can be positioned horizontally (on an exposure table or on the floor) for horizontal postures.

[0049] The exact disposition of cassettes on the support is disclosed in detail in the following.

[0050] Also the consequences of using cassettes instead of recording members in staggered disposition are explained.

[0051] In contrast to the situation disclosed in U.S. Ser. No. 09/035,528, where a single elongate cassette is used for image recording, the cassette comprising a number of partially overlapping recording members such as photostimulable phosphor screens or screen-film combinations and where the thickness of the individual recording members is small and negligible, the thickness of the individual cassettes in this invention is not negligible with regard to deformations of an image of the attenuating (rectangular) grid.

[0052] This image of the attenuating rectangular grid is used as a guidance for reconstructing the entire image of the elongate body from the partial images read out of the individual recording members that each carry a part of the image of the elongate body.

[0053] Because the recording members are not co-planar, the projection of the rectangles of the attenuating grid is different for each screen and is governed by three-dimensional projective geometry equations.

[0054] For the purpose of clarification, a world co-ordinate system (x,y,z) is defined such that the x-axis and y-axis point in the vertical and horizontal direction and the z-axis points in the direction perpendicular to the image plane (parallel to the path of X-ray incidence).

[0055] In the context of this invention three workable dispositions are considered; other combinations may be thought of without departing from the aim and scope of the current invention.

[0056] In the first disposition, depicted in FIG. 1 the recording members are perfectly parallel to the attenuating grid plane and their distance to the attenuation grid plane alternates between successive recording member. The larger the distance between a recording members and the attenuating grid, the larger the magnification of the projected grid rectangles will be. However, there is no projective distortion (the projective grid lines remain parallel).

[0057] In order to cover the entire elongate body, a plurality of cassettes is used in an overlapping disposition.

[0058] The situation is shown in FIGS. 1b, 2b, and 3b. Hidden borders of sub-images (sub-images or partial images are images of part of the elongate body recorded on a recording member in an individual cassette) are represented by a dashed line, visible sub-image borders by a continuous line.

[0059] Because of an alternating zoom factor between odd and even numbered recording members, this situation is referred to as the alternating overlapping disposition.

[0060] The second disposition, shown in FIG. 2, shows a staircase of overlapping recording members. The situation differs from the foregoing in that, except for the last sub-image, each sub-image will have its bottom part attenuated by the next recording member. Therefore, the overlap area can be treated in a uniform way throughout the image composition process. However, the zoom factor caused by the stepwise increasing distance between recording member plane and attenuating grid plane, is different for each recording member; the sub-images acquired by the farthest recording members may therefore need relatively significant down-zooming, or conversely, the sub-images acquired by the nearest recording members may need relatively significant up-zooming. This case is referred to as the staircase overlapping disposition.

[0061] The third case, depicted in FIG. 3 is a compromise between the drawbacks of the former two cases. It shows a sequence of slightly tilted recording members, the tilt angle being defined as a function of the ratio of recording member thickness to height, so that the resulting disposition fits in a minimal and perfectly vertical housing. This disposition features both uniform overlap configuration for each recording member and an identical projective distortion geometry. The projective distortion geometry is different from a simple zooming (magnification or reduction), however. Running from top to bottom, the vertical inter-line distance undergoes a linear increasing or decreasing zooming, depending on whether the tilt angle is positive or negative respectively. However, this projective distortion is equal for all imaging screens. For the disposition case depicted in FIG. 3, the inter-line distance decreases from top to bottom as schematised in FIG. 4. This acquisition modality is referred to as the oblique overlapping disposition.

[0062] For all cases presented, the recording members are mounted in a detachable mechanical set-up to which the attenuation grid and anti-scatter grid can be fixed, enabling to put the ensemble on an exposure table or on the ground for horizontal postures. The corresponding horizontal dispositions are depicted in FIG. 1c, FIG. 2c, FIG. 3c.

[0063] All dispositions shown in FIGS. 1, 2, 3 have overlapping cassettes. Dispositions with touching or separate non-overlapping cassettes may be thought of likewise. They are sub-optimal solutions because even in case of a touching disposition, the thickness of the housing of the recording member prevents the diagnostic area from being covered completely by the area of the recording member.

[0064] In addition to the three-dimensional projective geometry distortions caused by the three-dimensional physical disposition of the screens with respect to the attenuating grid, the distortions coped with in U.S. Ser. No. 09/035,528 remain equally well present.

[0065] For sake of clarity, these classic distortions are summarised below:

[0066] shift (y-axis dislocation) of a sub-image with respect to the next or previous sub-image,

[0067] overlap (in the x-axis) of a sub-image with respect to the next or previous sub-image,

[0068] rotation of the grid of one sub-image with respect to the grid of the next or previous sub-image,

[0069] the angle of said rotation being the combined result of (1) the non-perfect alignment of the borders of the recording member with respect to the perfectly horizontal and vertical axis, in the assumption that the attenuating grid lines are placed perfectly horizontally and vertically; (2) the non-perfectly straight feed-through in the digitising apparatus of the recording member.

[0070] In order to use the aforementioned recording process depicted in FIGS. 1, 2, 3, a suitable image processing method is therefore required to compensate not only for said projective geometry distortions, but also for the classical distortions dealt with in U.S. Ser. No. 09/035,528.

[0071] FIG. 4 shows a schematic representation of the grid recorded on one of the sub-images, with definition of co-ordinate axes and key quantities. The image has been acquired by the method of oblique overlapping disposition of recording members and hence the inter-line distance of vertically running grid lines narrows from top to bottom. The inter-line distance of horizontally running grid lines remains constant. However, the algorithm disclosed hereafter is able to cope with deformations of grid lines from the perfectly vertical and the perfectly horizontal simultaneously. The co-ordinate system's origin is in the upper left corner, with co-ordinate m running inside a row from left to right along different columns (in the horizontal direction), co-ordinate n running from top to bottom along different rows (in the vertical direction). The solid lines represent the imprint of the attenuating grid onto the image. In general, the raster formed by the grid lines can be tilted with an angle α^hwith respect to the horizontal axis. The rotation is due either to a rotated position of the recording member (e.g. an image screen) with respect to the cassette physical borders, or to a non-perfectly straight screen feed-through in the digitising apparatus, or a combination of both. Obviously, the resulting rotation angle α^his small.

[0072] Basic to the operations is the detect ion and characterisation of the contrasting grid which produces a rectangular pattern of less exposed pixels in each of the sub-images S_z,z=1 . . . Z. is the total number of images to be stitched and is typically 3 or 4. The method however generalises to any actual number of sub-images. The pattern formed by the grid in sub-image S_zis fully characterised by the knowledge of (1) the period P_z^hand P_z^vbetween the gridlines in horizontal resp. vertical direction, (2) line equations of the individual horizontal and vertical grid lines as disclosed by the current invention and (3) the offset O_z^hand O_z^vof horizontal and vertical gridlines with respect to the horizontal image border and the vertical image border respectively. Knowledge of the sets of these parameters is a necessary and sufficient condition for the reconstruction of the grid in the sub-images and for forming a composed total image of an elongate body part by image stitching.

[0073] In the sequel, (i,j), (m,n) or (x,y) denote image co-ordinates, (M,N) denote the image dimensions. It is assumed without loss of generality that all sub-images have equal dimensions.

[0074] FIG. 5 depicts the general flow of operations, the steps of which are further commented.

[0075] Grid Line Matrix Extraction

[0076] This step covers all low-level operations such as disclosed in U.S. Ser. No. 09/035,528 the whole contents of which is incorporated by reference into this application to extract the positions of grid lines in blocks of integrated line profiles. The output of this step is a sequence of m respectively n co-ordinates denoting the positions of the circles respectively crosses in FIG. 4.

[0077] Grid Line Fitting

[0078] As a result of prior art stated in U.S. Ser. No. 09/035,528, two grid matrices are determined, g^hfor the horizontal resp. g^vfor the vertical direction, storing the indexes of all grid lines detected in such a manner that the resulting arrangement of non-zero entries in a column of g^hresp. a row of g^vall correspond to the row resp. column locations of the same physical grid line. By denoting with Q_hand Q_vthe number of blocks in the horizontal resp. vertical direction and denoting the number of detected grid lines in the sequences G^hand G^vby N_gresp. M_g, the dimensions of g^hresp. g^vare Q^h×N_gresp. Q^v×M_g, such that there are N_ghorizontally oriented grid lines detected and M_gvertically oriented grid lines detected. The sequences of grid line co-ordinates are input to a line fitting step so as to model each grid line with the equation of a straight line. To this purpose, the row resp. column co-ordinates of the grid matrices g^hand g^vare subjected to a linear regression, commonly known in the prior art such as T. Pavlidis, Algorithms for Graphics and Image Processing, Computer Science Press, 1982, pp. 281-292, Polygonal Approximations. Each linear regression results in the offset and slope of the fit line corresponding to a grid line. For the N_ghorizontally oriented grid lines there are Q_hdata points to a fit a line to, for the M_gvertically oriented grid lines there are Q_vdata points available for fitting. The equations of the resulting horizontally resp. vertically oriented fit lines are $\begin{array}{ccc}{L}_{l^h}^{h,z}\left(m\right)={\alpha }_{l^h}^{h,z}m+{\beta }_{l^h}^{h,z}& l^h=0\dots N_{g^h}-1& z=1\dots Z\\ L_{l^v}^{v,z}\left(n\right)={\alpha }_{l^v}^{v,z}n+{\beta }_{l^v}^{v,z}& l^v=0\dots M_{g^v}-1& z=1\dots Z\end{array}$

[0079] in which the coefficients α and β are the slope and the axis offset of the fit lines. Obviously, any other parameterizable curve could model the grid lines, depending on the nature of the projective geometry distortion. However, because of the pinhole X-ray illumination and the planar geometric disposition of recording member with respect to grid plane, straight lines remain straight lines under projection, so that straight line fitting effectively models the grid line data points.

[0080] Grid Calibration and Deformation Measurement

[0081] The line equations of the grid lines in each sub-image serve as input to determine the final position of the grid line in the stitched image. Each horizontally running grid line will be de-skewed to the perfectly horizontal, the position of which is determined in this step. Analogously, each vertically running grid line will be de-skewed to the perfectly vertical, the position of which is determined in this step.

[0082] Furthermore, the interline distance is made equal to the periodicity distance, for all grid lines running horizontally and vertically in the sub-images. The periodicity distance is derived in the “Grid line matrix extraction” step.

[0083] The result of the grid calibration step is that all previously enumerated deformations of the ideal rectangular grid are captured in an associated set of skew line equations, given in the “Grid line fitting” step. The deformations corrected for include shift of one sub-image with respect to the other sub-images, rotation of one sub-image with respect to the other sub-images, varying interline distance due to the zoom effect caused by different image screen-attenuating grid distances, and projective geometry distortions caused by non-parallel screen-grid configurations.

[0084] More specifically, the quantity P^vdenotes the minimum of all periodicity vectors between vertically running grid lines of the ensemble of Z input sub-images. $\begin{array}{cc}{P}^v=\min \left(P_z^v\right)& z=1\dots Z\end{array}$

[0085] and will be used to normalise the distance between vertically running grid lines. For a square physical grid, P^vwill be used to normalise the distance between horizontally running grid lines as well, so that a perfectly square grid will be reconstructed on completion of the warping process. For a rectangular grid, a similar quantity P^hdenotes the minimum of all periodicity vectors between horizontally running grid lines of the ensemble of Z input sub-images. $\begin{array}{cc}{P}^h=\min \left(P_z^h\right)& z=1\dots Z\end{array}$

[0086] Taking the minimum of the set of periodicity vectors ensures that each warped sub-image will fit the original sub-image size so that no diagnostic information is lost.

[0087] The position {circumflex over (β)}₀^v,zof the first vertical warped grid line and the positions of the remaining set of vertically warped grid lines are determined as follows $\begin{array}{cc}{\hat{\beta }}_0^{v,z}=\frac{1}{Z}\sum _{z=1}^Z{L}_0^{v,z}\left(N/2\right)& \\ {\hat{\beta }}_{l^v}^{v,z}={\hat{\beta }}_0^{v,z}+l^v{P}^v& l^v=1\dots M_{g^v}-1,\end{array}$

[0088] which positions are equal for all sub-images S_z,z=1 . . . Z.

[0089] The position {circumflex over (β)}₀^h,zof the first horizontally warped grid line and the positions of the remaining set of horizontally warped grid lines is determined as follows $\begin{array}{cc}{\hat{\beta }}_0^{h,z}=\max \left(L_0^{h,z}\left(0\right),L_0^{h,z}\left(M-1\right)\right)& \\ {\hat{\beta }}_{l^h}^{h,z}={\hat{\beta }}_0^{h,z}+l^h{P}^h& l^h=1\dots N_{g^h}^z-1,\end{array}$

[0090] which positions are different for all sub-images S_z,z=1 . . . Z.

[0091] The skew values for the horizontally resp. vertically oriented grid lines are given by $\begin{array}{cc}{\Delta }^{h,z}\left(m,{\hat{\beta }}_{l^h}^{h,z}\right)=L_{l^h}^{h,z}\left(m\right)-{\hat{\beta }}_{l^h}^{h,z}& l^h=0\dots N_{g^h}-1\\ {\Delta }^{v,z}\left({\hat{\beta }}_{l^v}^{v,z},n\right)=L_{l^v}^{v,z}\left(n\right)-{\hat{\beta }}_{l^v}^{v,z}& l^v=0\dots M_{g^v}-1\end{array}$

[0092] for all sub-images S_z,z=1 . . . Z.

[0093] Two 2-dimensional warping matrices are defined for the horizontal resp. vertical de-skewing operation by linearly interpolating the skew values for all positions in between two successive grid lines $\begin{array}{cc}{\Delta }^{h,z}\left(m,n\right)=\frac{\left(n-{\hat{\beta }}_{l^h}^{h,z}\right){\Delta }^{h,z}\left(m,{\hat{\beta }}_{l^h+1}^{h,z}\right)+\left({\hat{\beta }}_{l^h+1}^{h,z}-n\right){\Delta }^{h,z}\left(m,{\hat{\beta }}_{l^h}^{h,z}\right)}{P^h}& {\hat{\beta }}_{l^h}^{h,z}\le n\le {\hat{\beta }}_{l^h+1}^{h,z}\\ & l^h=0\dots N_{g^h}-2\\ {\Delta }^{v,z}\left(m,n\right)=\frac{\left(m-{\hat{\beta }}_{l^v}^{v,z}\right){\Delta }^{v,z}\left({\hat{\beta }}_{l^v+1}^{v,z},n\right)+\left({\hat{\beta }}_{l^v+1}^{v,z}-m\right){\Delta }^{v,z}\left({\hat{\beta }}_{l^v}^{v,z},n\right)}{P^v}& {\hat{\beta }}_{l^v}^{v,z}\le m\le {\hat{\beta }}_{l^v+1}^{v,z}\\ & l^v=0\dots M_{g^v}-2\end{array}$

[0094] for all sub-images S_z,z=1 . . . Z.

[0095] Sub-Image Warping by Piece-Wise Linear Separable De-Skewing

[0096] This module will transform the sub-image such that (1) each vertical fitted grid line is de-skewed to the perfectly vertical target line and the ensemble of all vertical target lines is equidistant; (2) each horizontal fitted grid line is de-skewed to the perfectly horizontal target line and the ensemble of all horizontal target lines is equidistant; (3) the skew applied to each pixel in between the vertical target lines (respectively horizontal target lines) is obtained by linear interpolation of the skew associated with the left (upper) fit line and the right (lower) fit line. Therefore, the total operation is termed image warping since the local displacements are place dependent and performed by piece-wise linear de-skewing, the skew values varying piece-wise linearly along the skew direction and their control points being defined by the skew values associated with the grid fit lines.

[0097] Since the skew values are floating point values, the problem of image re-sampling on a discrete grid needs to be addressed. Interpolation is the process of determining the values of a function at positions lying between its samples. Since interpolation reconstructs the signal lost in the sampling process by smoothing the data samples with the interpolation function, an important issue is how image sharpness can be retained at the maximum. The performance of the interpolation in the passband and stopband, the numerical accuracy and the computational cost of interpolating algorithms are therefore dependent on the interpolation kernel. In addition to the choice of the interpolation algorithm, a computational scheme has been devised to integrate the interpolation process in the aforementioned de-skew operation. Sampling theory establishes that the sinc function is the ideal interpolation kernel. Although this interpolation filter is exact, it is not practical since it is an IIR filter defined by a slowly converging infinite sum. The next to optimal interpolating kernel with local extent over 4 points is offered by the interpolating cubic B-spline. The interpolated data value is obtained as followed for the horizontal de-skew operation ${\hat{S}}_z^h\left(m,n\right)=\sum _{k=-2}^2{T}_z^h\left(m+k,n\right)h\left({\Delta }^{h,z}\left(m,n\right)+k\right)$

[0098] with the cubic B-spline interpolating kernel h given by the following piece-wise cubic polynomials, sometimes referred to as the Parzen window $h\left(x\right)=\frac{1}{6}\{\begin{array}{cc}3{\left|x\right|}^3-6{\left|x\right|}^2+4& 0\le \left|x\right|<1\\ -{\left|x\right|}^3+6{\left|x\right|}^2-12\left|x\right|+8& 1\le \left|x\right|<2\\ 0& 2\le \left|x\right|\end{array}$

[0099] with Δ^h(m,n) the horizontal 2-dimensional skew matrix defined above. The matrix T_z^h(m,n) is obtained by solving the following set of N equations, one set for each row n defined by a tri-diagonal matrix. $\left[\begin{array}{c}{S}_z\left(0,n\right)\\ S_z\left(1,n\right)\\ S_z\left(2,n\right)\\ \dots \\ S_z\left(m-2,n\right)\\ S_z\left(m-1,n\right)\end{array}\right]=\left[\begin{array}{cccccccc}4& 1& & & & & & \\ 1& 4& 1& & & & & \\ & 1& 4& 1& & & & \\ & & & \dots & & & & \\ & & & & & 1& 4& 1\\ & & & & & & 1& 4\end{array}\right]\left[\begin{array}{c}{T}_z^h\left(0,n\right)\\ T_z^h\left(1,n\right)\\ T_z^h\left(2,n\right)\\ \dots \\ T_z^h\left(m-2,n\right)\\ T_z^h\left(m-1,n\right)\end{array}\right],n=0\dots N-1$

[0100] The de-skew operation is executed independently in the horizontal and vertical direction. The result of the output of the horizontal de-skew Ŝ_z^h(m,n) is input to the vertical de-skew operation using the same de-skew algorithm to produce the final rotated image Ŝ_z(m,n)=Ŝ_z^hv(m,n). The result image will show the deformed grid geometrically corrected to its original form, that is, as two sets of mutually orthogonal equally spaced perfectly straight and parallel lines. Furthermore, since the offset (distance to the image border) of the first vertical grid line is equal for all sub-images, the operation will ensure the perfect continuity (co-linearity) of the vertically oriented de-skewed grid lines in the image composition. Analogously, the perfect continuity of periodicity of horizontal corrected grid lines is guaranteed by restrictions imposed upon the overlap computation and determination as disclosed in U.S. Ser. No. 09/035,528.

[0101] Other image interpolation schemes are devisable based on prior art such as presented in e.g. G. Wolberg, Digital Image Warping, IEEE Computer Society Press Monograph, 1990.

[0102] Overlap Image Correlation

[0103] This step proceeds as disclosed in U.S. Ser. No. 09/035,528, using the warped sub-images and the grid model parameters as input to the process.

[0104] Image Stitching

[0105] This step proceeds as disclosed in U.S. Ser. No. 09/035,528, using the warped sub-images and the overlap values as input to the process.