11,483,586 members (68,912 online)

# Detect image skew angle and deskew image

, 13 Feb 2013 CPOL 95.1K 10.2K 123
 Rate this:
Deskew an image by converting it to grayscale + edges and rotating by the detected skew angle.

## So! What is image skew?

Using the above drawn figure as reference, I can say - the theta pointed out by that bluish arrow sign is the image skew .

To put it in proper words, the angle by which the image seems to be deviated from its perceived steady position is the image skew.

## Some pre-requisites for detecting the skew

Before detecting the skew, the first step is to differentiate between a text image and a regular image. By text image, I mean images which are scanned documents or screenshot of a text document, or in other words, images that contains letters and texts.

And a regular image in my terms is a picture, or photo of a scenery, or something drawn by Mr. Van Gogh and Co.

In one of my articles, I have covered the topic of a basic algorithm to differentiate between regular image and text images.

However, although the skew detection of a text image is relatively easier to implement, for a regular image, it's quite difficult and next to impossible in some cases. Suppose I have a picture of only the head part of my favourite celebrity Ms. Julia Roberts, ... there is no way to know from the image whether she really has her head slanted towards one direction, or the camera man who took the beautiful picture had a problem with the camera angle.

Therefore, for regular images, we have to make some sort of assumptions. The case we are going to discuss here will assume, the pictures to be deskewed will be framed. Which means, the pictures that we are dealing with are pictures taken from another picture having frames/border around them. You can refer to the application screenshot at the very top as an example. The screenshot you see is of a camera captured image of another picture having a thick white-ish border around it.

## The basic concept of detecting the skew

The main idea to detect skew for a regular image and a text image is the same. We have to first turn the pictures to Gray Scale.

Then, for regular (non text) images, we have to find the edges or feature lines of the image. There are many algorithms to detect the edges of an image, and in our case, we used the canny edge detection algorithm. We can wish to skip this part for text images.

The next step and the most important step is to cast rays from one side of the picture to the other, and by using the intersection info of the rays with the various parts of the images, we come up with a good skew angle.

As you can see from the above image, we have first converted the image (the one you see in the application screenshot) to its very basic feature lines/edges and have casted horizontal rays from one side to the other. The rays which we have casted (I mean the horizontal lines) intersect with the border/frame of the image. If I put the intersection points of the rays and the left edge of the picture in an X, Y plot (where X is horizontal and Y is vertical), I will end up with some sort of straight line curve that has a slope angle close to 2 degrees with respect to the Y axis. So, there you have it, your skew angle for the image in the application screenshot is close to 2 degrees or 0.035 radians.

This is more or less the basic concept used for detecting the skew of a regular framed image. You have to keep in mind this idea is applicable only to regular images and not text images.

By the way, the uploaded source code's skew detection algorithm does a little more than what I have explained above. What it does is - cast rays from left to right to intersect with the left edge, then top to bottom, then right to left, and then bottom to top. The algorithm considers the one with the most intersection points, and creating consistently occurring slopes as the edge to consider as the reference edge. For example, if the bottom to top rays intersects most with the bottom edge of the image and more frequently occurring slopes are observed, then that edge has the priority and the points of the bottom to top rays would be considered to be put in the X,Y plot. All these extra effort is just for the fact that, in some images, the border edges might be broken and not as beautifully framed fully around the image, like you are seeing in the application screenshot image. So, we'll have a set of rays with different slopes for such images, because the rays will go beyond the border edge wherever there is a broken edge line, and will create a different angle with another intersected point that lies somewhere inside the image bounded by the frame. So, the most occurring slopes (if two slopes have less than 0.0349 radians / 2.0 degrees difference between them - we'd consider the slopes as same valued) will be taken into consideration to detect the skew.

In order to detect the skew of text images, we don't have to worry too much about the picture having border edges or not. As mentioned before, we cast rays from left to right and find how many blackish pixels the rays intersect with. Then, we rotate the angle of the rays by a small amount and cast the rays again, and do the same process again and again.

The angle at which most white space is encountered (in other words, most black pixels are shot down, i.e., intersected), we consider that to be the skew angle. Most likely, we would start from 0 degree and go down to 90 degrees, and then again, from 0 degree go up to 90 degrees with the angle of the rays. Whenever the white ratio reaches the maximum peak and then starts to come down again, we stop immediately, and we know we have found our skew angle. This is the basic idea of detecting skew for text images.

The deskewing of an image is less of a trouble. All you need to do is rotate the image by the same amount of the skew, only to the reverse direction.

The code for the rotation is available in the ImageFunctions.h header file of the uploaded source.

However, if you want a real nice quality rotated image, you can use the wonderful `aarot` class by Mark Gordon (http://codeguru.earthweb.com/cpp/cpp/algorithms/math/article.php/c10677/Anti-Aliased-Image-Rotation-Aarot.htm). Since the rotation code inside this class performs antialiasing, the performance is a little slow, but the output image quality is very good.

## Using the code

In the uploaded source code, mainly two classes are used for detecting the skew. One class, the main class, is named `SkewManager`, which has code to detect the skew, and the other class is `ImageData`, which is just a wrapper for the bitmap bits and also keeps the width and the height of the image to be used by the `SkewManager` to apply its algorithms to the bitmap bits.

The `SkewManager` expects the bitmap that will be passed to be gray scaled. The blurring of the image to its edges (if required) is done inside the `ImageData` class.

Following is the example of the usage of the two classes (the `m_` prefixed variables as well as the `skewMan` variable are member variables of the `CScannedDocTestDoc` class).

```BOOL CScannedDocTestDoc::OnOpenDocument(LPCTSTR lpszPathName)
{
m_Image.Destroy();

if(SUCCEEDED(hr))
{
int nPitch = m_Image.GetPitch();
int nWidth = m_Image.GetWidth();
int nHeight = m_Image.GetHeight();

int nBytesPerPixel = m_Image.GetBPP() / 8;
if(nBytesPerPixel)
{
byte* pGrayScaleBits = NULL;
// Get a gray scaled bitmap and bitmap bits from the original image
m_hBmpGrayScale = ImageFunctions::GetGrayScaleImage((HBITMAP)m_Image,
nWidth, nHeight, nBytesPerPixel, nPitch, &pGrayScaleBits);
if(pGrayScaleBits != NULL)
{
// Find out if the image is a regular image or text image
m_bTextImage = ImageFunctions::IsTextImage(pGrayScaleBits, nWidth, nHeight);

// get the skew angle
double angle = 0.0;
skewMan.GetSkewAngle(m_bTextImage, m_imgData, angle);
m_dSkewAngle = angle;
}
}
return TRUE;
}
return FALSE;
}```

## Points of interest

In order to achieve skew detection and deskewing of images, I had to go through a lot of trouble to perform basic image operations such as creating a gray scaled image, rotating an image etc., and I have put all those functionalities inside one header, ImageFunctions.h. Most of it was improvisation on already available materials in the Internet and MSDN. The canny edge detection algorithm can be found in the the header canny_edge.h.

## CSharpClient

This is a cut down version of the deskewimage code built to compile into a DLL.  A small C# demo application is included in the solution showing how to call the code contained in within the DLL,  pass in a file name to be analysed and the skew angle in radians is returned and displayed on the screen.

## Acknowledgements

• Kathey Marsden, Professor Richard J. Fateman
• for the OCRchie project - the foundation upon which, the text image skew detection algorithms of this article's codes are based upon.

• Also, thanks to Archie Russell, James Hopkin, and Cynthia Tian, who contributed significantly to the original design.
• Heath, M., Sarkar, S., Sanocki, T., and Bowyer, K. Comparison of edge detectors: a methodology and initial study, Computer Vision and Image Understanding 69 (1), 38-54, January 1998.
• Heath, M., Sarkar, S., Sanocki, T. and Bowyer, K.W. A Robust Visual, Method for Assessing the Relative Performance of Edge Detection Algorithms, IEEE Transactions on Pattern Analysis and Machine Intelligence 19 (12), 1338-1359, December 1997.
• for the excellent canny edge detection algorithm.

• for providing with reading materials related to deskewing text and framed images.

## History

• Article uploaded: 24 August, 2010.
• Article Updated 29 May 2012 CsharpClient Added
• Article Updated 14 February 2013 CSharpClient.zip 1.1 Added
Updated to compile in Visual Studio 2008 (Professional Edition)
Combined into a single VS Solution containing two projects (for easier debugging)
Minor Bug Fixes
Angle of 999 is now returned if the document could not be analysed (not a 8bpp image)
Added a Text box to show the name of the file analysed

## Share

Chief Technology Officer Rational Technologies
Your Lord does not like people who waste by extravagance, and only the righteous C s, and does so in increments, and they are of those people who are mindful of what memory and time has in offer and they waste none of it.

Australia
No Biography provided

 First PrevNext
 Machine learning may be more appropriate Maruf Maniruzzaman8-Mar-15 10:58 Maruf Maniruzzaman 8-Mar-15 10:58
 Memory Leak David_Pollard4-Apr-13 15:11 David_Pollard 4-Apr-13 15:11
 Re: Memory Leak Mukit, Ataul5-Apr-13 8:57 Mukit, Ataul 5-Apr-13 8:57
 Re: Memory Leak David_Pollard5-Apr-13 13:57 David_Pollard 5-Apr-13 13:57
 Re: Memory Leak Mukit, Ataul5-Apr-13 19:40 Mukit, Ataul 5-Apr-13 19:40
 Re: Memory Leak David_Pollard5-Apr-13 19:54 David_Pollard 5-Apr-13 19:54
 I'm only testing now with the C# demo application that is attached to this article.... OK I found the problem. There may be others but this is certainly the bigest leak to start with. The file "ImageFunctions.h" contains the method named GetGrayScaleImage(). The line within the if statement shown below byte* p8Bit = new byte [grayScaleSize]; Aloocates one byte of memory for every pixel in the image being tested. I can't just dealocate this a the end of the method because obviously the calling method wants to use it. I have done some reading today about constructors and destructors and they don't appear to be used. It seems that the method that calls GetGrayScaleImage() should create the p8bit thing and then destroy it when it is finished. I'm easily confused by the C++ syntax How can I destroy p8Bit thing? Even if I do it just before it is created again would ge good enough. Any ideas? Thanks David ```static HBITMAP GetGrayScaleImage(HBITMAP hBmp, int width, int height, unsigned int bytesPerPixel, int pitch, byte** ppGrayScaleBits = NULL) { if(bytesPerPixel == 0 || bytesPerPixel == 2) return NULL; int width_in_bytes = width * bytesPerPixel; int unusedBytes = absolute(pitch) - width_in_bytes; if(!(unusedBytes % 3)) unusedBytes = 0; HBITMAP hBmp8 = Create8bppBitmap(NULL, width + unusedBytes, height, NULL);   HDC hdcScr = ::GetDC(NULL);   HDC hdc = ::CreateCompatibleDC(hdcScr); HBITMAP hBmpOld = (HBITMAP)::SelectObject(hdc, hBmp);   HDC hdc8 = ::CreateCompatibleDC(hdcScr); HBITMAP hBmpOld8 = (HBITMAP)::SelectObject(hdc8, hBmp8);   ::BitBlt(hdc8, 0, 0, width, height, hdc, 0, 0, SRCCOPY); if(ppGrayScaleBits) { int grayScaleSize = width * height;   byte* p8Bit = new byte [grayScaleSize]; ::GetBitmapBits(hBmp8, grayScaleSize, p8Bit);   *ppGrayScaleBits = p8Bit;   } ::SelectObject(hdc, hBmpOld); ::SelectObject(hdc8, hBmpOld8);   ::DeleteDC(hdc); ::DeleteDC(hdc8);   ::ReleaseDC(NULL, hdcScr); return hBmp8; }```
 Re: Memory Leak Mukit, Ataul6-Apr-13 9:22 Mukit, Ataul 6-Apr-13 9:22
 Re: Memory Leak David_Pollard10-Apr-13 20:32 David_Pollard 10-Apr-13 20:32
 Compile Errors Member 97749616-Feb-13 13:12 Member 9774961 6-Feb-13 13:12
 Re: Compile Errors Mukit, Ataul7-Feb-13 6:28 Mukit, Ataul 7-Feb-13 6:28
 Re: Compile Errors Member 97749617-Feb-13 15:56 Member 9774961 7-Feb-13 15:56
 Re: Compile Errors Member 977496111-Feb-13 17:29 Member 9774961 11-Feb-13 17:29
 Re: Compile Errors Mukit, Ataul11-Feb-13 18:07 Mukit, Ataul 11-Feb-13 18:07
 Re: Compile Errors Member 977496111-Feb-13 19:20 Member 9774961 11-Feb-13 19:20
 Re: Compile Errors Mukit, Ataul11-Feb-13 21:08 Mukit, Ataul 11-Feb-13 21:08
 Re: Compile Errors David_Pollard11-Feb-13 22:04 David_Pollard 11-Feb-13 22:04
 Re: Compile Errors Mukit, Ataul12-Feb-13 7:21 Mukit, Ataul 12-Feb-13 7:21
 Re: Compile Errors Mukit, Ataul13-Feb-13 10:25 Mukit, Ataul 13-Feb-13 10:25
 Re: Compile Errors David_Pollard13-Feb-13 12:55 David_Pollard 13-Feb-13 12:55
 Re: Compile Errors Mukit, Ataul13-Feb-13 18:10 Mukit, Ataul 13-Feb-13 18:10
 Re: Compile Errors David_Pollard13-Feb-13 18:16 David_Pollard 13-Feb-13 18:16
 Re: Compile Errors Mukit, Ataul13-Feb-13 18:36 Mukit, Ataul 13-Feb-13 18:36
 Re: Compile Errors David_Pollard13-Feb-13 18:43 David_Pollard 13-Feb-13 18:43
 Re: Compile Errors David_Pollard13-Feb-13 20:32 David_Pollard 13-Feb-13 20:32
 thanks for the code but.. Christian Lin7-Jan-13 16:43 Christian Lin 7-Jan-13 16:43
 Last Visit: 31-Dec-99 19:00     Last Update: 25-May-15 0:32 Refresh 123 Next »