A projection device configured to project a first calibration pattern onto a projection screen, wherein the first calibration pattern comprises a plurality of first feature points (u, v), said projection device comprising a processor circuit configured to carry out the steps of: acquiring a second calibration pattern by capturing, from a viewer position, an image of the first calibration pattern projected onto the projection screen, wherein the second calibration pattern comprises a plurality of second feature points (x, y), generating a similarity matrix based on the plurality of first feature points (u, v) and on the plurality of second feature points (x, y), wherein the step of generating the similarity matrix comprises: acquiring a largest circumscribed quadrilateral according to an arrangement of the plurality of first feature points (u, v), acquiring a largest inscribed quadrilateral according to an arrangement of the plurality of second feature points (x, y), acquiring a scaling relationship (s w, s h) and a translation relationship (t w, t h) between the largest circumscribed quadrilateral and the largest inscribed quadrilateral, based on four vertex coordinates of the largest circumscribed quadrilateral and four vertex coordinates of the largest inscribed quadrilateral, and generating the similarity matrix based on the scaling relationship (s w, s h) and the translation relationship (t w, t h); converting by the processor circuit, the plurality of second feature points (x, y) to a plurality of warped feature points (s, t) based on the similarity matrix; converting a first mapping table to a second mapping table according to the w plurality of warped feature points (s, t), wherein the first mapping table and second mapping table are two-dimensional matrices, wherein the first mapping table has the same resolution as a resolution of the projection device, and wherein each value of the first mapping table represents a position where each pixel of an input image is to be mapped; and generating a reverse-warped image from the input image based on the second mapping table; wherein said projection device is further configured to project the reverse-warped image onto the projection screen. acquiring a second calibration pattern by capturing, from a viewer position, an image of the first calibration pattern projected onto the projection screen, wherein the second calibration pattern comprises a plurality of second feature points (x, y), generating a similarity matrix based on the plurality of first feature points (u, v) and on the plurality of second feature points (x, y), wherein the step of generating the similarity matrix comprises: acquiring a largest circumscribed quadrilateral according to an arrangement of the plurality of first feature points (u, v), acquiring a largest inscribed quadrilateral according to an arrangement of the plurality of second feature points (x, y), acquiring a scaling relationship (s w, s h) and a translation relationship (t w, t h) between the largest circumscribed quadrilateral and the largest inscribed quadrilateral, based on four vertex coordinates of the largest circumscribed quadrilateral and four vertex coordinates of the largest inscribed quadrilateral, and generating the similarity matrix based on the scaling relationship (s w, s h) and the translation relationship (t w, t h); converting by the processor circuit, the plurality of second feature points (x, y) to a plurality of warped feature points (s, t) based on the similarity matrix; converting a first mapping table to a second mapping table according to the w plurality of warped feature points (s, t), wherein the first mapping table and second mapping table are two-dimensional matrices, wherein the first mapping table has the same resolution as a resolution of the projection device, and wherein each value of the first mapping table represents a position where each pixel of an input image is to be mapped; and generating a reverse-warped image from the input image based on the second mapping table; wherein said projection device is further configured to project the reverse-warped image onto the projection screen. acquiring a largest circumscribed quadrilateral according to an arrangement of the plurality of first feature points (u, v), acquiring a largest inscribed quadrilateral according to an arrangement of the plurality of second feature points (x, y), acquiring a scaling relationship (s w, s h) and a translation relationship (t w, t h) between the largest circumscribed quadrilateral and the largest inscribed quadrilateral, based on four vertex coordinates of the largest circumscribed quadrilateral and four vertex coordinates of the largest inscribed quadrilateral, and generating the similarity matrix based on the scaling relationship (s w, s h) and the translation relationship (t w, t h); converting by the processor circuit, the plurality of second feature points (x, y) to a plurality of warped feature points (s, t) based on the similarity matrix; converting a first mapping table to a second mapping table according to the w plurality of warped feature points (s, t), wherein the first mapping table and second mapping table are two-dimensional matrices, wherein the first mapping table has the same resolution as a resolution of the projection device, and wherein each value of the first mapping table represents a position where each pixel of an input image is to be mapped; and generating a reverse-warped image from the input image based on the second mapping table; wherein said projection device is further configured to project the reverse-warped image onto the projection screen.