{"id":292,"date":"2018-12-04T13:29:42","date_gmt":"2018-12-04T17:29:42","guid":{"rendered":"http:\/\/dominionsw.com\/wordpress\/?p=292"},"modified":"2019-11-25T07:47:17","modified_gmt":"2019-11-25T11:47:17","slug":"classifying-tree-seedlings-with-machine-learning-part-1","status":"publish","type":"post","link":"https:\/\/www.dominionsw.com\/?p=292","title":{"rendered":"Classifying tree seedlings with Machine Learning, Part 1."},"content":{"rendered":"
\"\"<\/figure>\r\n\r\n\r\n\r\n

I’ve been working with a non-profit, GreenStand<\/a>, helping them with image processing tasks and related machine learning features for the technology part of their mission.<\/p>\r\n\r\n\r\n\r\n

GreenStand is (really) making the world a better place – walkin’ the walk.\u00a0 In short, they’re investing in people and local communities by encouraging them, financially and otherwise, to plant trees, because many of these communities have been deforested, for fuel and other reasons. You can read more about their mission here<\/a>.\u00a0<\/p>\r\n\r\n\r\n\r\n

When I approached them about their needs, one of the first things that came up was how to identify “bad” photos. First, one has to define what a “bad” photo is. It turns out that the planters upload a photo of the newly planted seedling to demonstrate to the group that it has been planted. Since the cell phones used are often old and with low resolution cameras, and submitted over a typically spotty cell phone connection, the photos can be out of focus or distorted. Sometimes people upload the wrong photo, so the photo might be of a car, or some random scene. In addition, the photos are not taken in such a way as to make identification easier – by having the same pose, focal length, or removing other stray vegetation around the seedling.\u00a0<\/p>\r\n

Duplicate photos<\/h2>\r\n

Another problem to be solved – and this is where the machine learning \/ deep learning comes in –\u00a0 is caused by people uploading in batches and the same photos are uploaded over and over again. Knowing the basic type of tree (if not the exact species) and its geolocation in the meta data of the photo, we can flag it as a potential duplicate, and keep from overloading the server with data.<\/p>\r\n\r\n\r\n\r\n

Image processing for any task is completely dependent on the data – what approaches you take, the algorithms you choose, the results that you can attain – depend entirely upon the data you are given. Fortunately, the folks at GreenStand have been up and running long enough to give me over 50000 photos to start with. I’m going to be doing prototyping for this project with Matlab. Ultimately, to run on the server, likely to be a cloud based solution, I will have to port the code over to an open source solution – probably Python and OpenCV. But for rapid prototyping of image processing and machine learning, you really cannot beat Matlab, in my opinion. It’s pricey but it saves a lot of time.<\/p>\r\n\r\n\r\n\r\n

The latest version of Matlab that I’m using, 2018b has an integrated image browser “app”, that let’s the user browse a folder of images and load them into various other apps or just import them into the user’s matlab workspace. I loaded up the images and browsed through them, getting an idea of what the “good” or “bad” image characteristics might be.<\/p>\r\n\r\n\r\n\r\n

\"\"\"\"<\/a>\r\n
Image Browser with out of focus images<\/figcaption>\r\n<\/figure>\r\n\r\n\r\n\r\n

The image browser let’s you load images directly into other “apps”, such as the color segmentation app. This will come in very handy later. But for now I’m just looking at the out of focus images.<\/p>\r\n\r\n\r\n\r\n

 <\/p>\r\n\r\n\r\n\r\n

“Focus” algorithms<\/h2>\r\n\r\n\r\n\r\n

Determining if an image is out of focus is really determining the sharpness of the image – how “crisp” the image is. This means measuring contrast. There are quite a few algorithms that have been developed for this purpose. After reading up on the literature, and testing a few, I decided upon using the “Brenner Focus Algorithm”.\u00a0

https:\/\/cs.uwaterloo.ca\/~vanbeek\/Publications\/spie2014.pdf<\/a><\/p>\r\n\r\n\r\n\r\n

I implemented a variation on this algorithm by taking the average of the maximum of the square of the vertical and horizontal gradients:<\/p>\r\n

 <\/p>\r\n Rr = \\{1...rows\\} <\/span>\r\n Rc = \\{1...cols\\} <\/span>\r\n \\forall r \\in Rr ,\\forall c \\in Rc <\/span>\r\n \\bigtriangledown H=I(r,c) - I(r,c-2) <\/span>\r\n \\bigtriangledown V=I(r,c) - I(r-2,c) <\/span>\r\n M = max(\\bigtriangledown H,\\bigtriangledown V) <\/span>\r\n F = mean(M) <\/span>\r\n

Matlab Code<\/p>\r\n

function [Measure] = brenners(Image)
\u00a0 \u00a0 [M, N] = size(Image);
\u00a0 \u00a0 DH = zeros(M, N,’single’);
\u00a0 \u00a0 DV = zeros(M, N,’single’);
\u00a0 \u00a0 DHG = gpuArray(DH);
\u00a0 \u00a0 DVG = gpuArray(DV);
\u00a0 \u00a0 IG = gpuArray(Image);
\u00a0 \u00a0 DVG(1:M-2,:) = IG(3:end,:)-IG(1:end-2,:);
\u00a0 \u00a0 DHG(:,1:N-2) = IG(:,3:end)-IG(:,1:end-2);
\u00a0 \u00a0 Mx = max(abs(DHG), abs(DVG));
\u00a0 \u00a0 M2 = Mx.^2;
\u00a0 \u00a0 M2CPU = gather(M2);

\u00a0 \u00a0 Measure = mean2(M2CPU);
end<\/p>\r\n\r\n\r\n\r\n

Java Code<\/p>\r\n

\/\/ Compute average of sum of squares of the gradient in H and V directions.\r\n\tprivate static double brennersFocusMetric(int[][] input,int rows, int cols ) {\r\n\r\n\t\tint[][] V = new int[rows][cols];\r\n\t\tint[][] H = new int[rows][cols];\r\n\t\tfor(int row = 0; row < rows; row++)\r\n\t\t{\r\n\t\t\tfor (int col = 0; col < cols-2; col++) {\r\n\t\t\t\tint grad = input[row][col+2] - input[row][col];\r\n\t\t\t\tH[row][col] = grad;\r\n\t\t\t}\r\n\t\t}\r\n\r\n\t\tfor(int row = 0; row < rows-2; row++)\r\n\t\t{\r\n\t\t\tfor (int col = 0; col < cols; col++) {\r\n\t\t\t\tint grad = input[row+2][col] - input[row][col];\r\n\t\t\t\tV[row][col] = grad;\r\n\t\t\t}\r\n\t\t}\r\n\r\n\t\tdouble sum = 0;\r\n\t\tfor(int row = 0; row < rows; row++)\r\n\t\t{\r\n\t\t\tfor (int col = 0; col < cols; col++) {\r\n\t\t\t\tdouble HRC = H[row][col];\r\n\t\t\t\tdouble VRC = V[row][col];\r\n\t\t\t\tsum  += Math.abs(HRC) > Math.abs(VRC) ? HRC * HRC : VRC * VRC;\r\n\t\t\t}\r\n\t\t}\r\n\t\treturn sum \/ (double)(rows * cols);\r\n\t}<\/pre>\r\n\r\n\r\n\r\n
The above code finds the “bad” photos well – based on my test, it eliminated all the out of focus photos, with a less than 5 percent false positive rate. Some of these false positives were the result of a dual focus – the seedling was in focus, but the background was not. The algorithm returns a number that describes a measure, somewhat subjective, of the “goodness” of focus. One will need to experiment to find the range values best suited for their own application. Of course, more parameters of what is “bad” can be used to refine the meaning of the term. The goal is to warn the user that their image of the seedling is not good, and they should retake the photo.<\/p>\r\n
Next….<\/h2>\r\n
In the next of this series of 4 posts, I’ll go into segmentation of the images and preparation of them for testing a Support Vector Machine that classifies the images into four basic types…<\/p>\r\n
 <\/p>\r\n
 <\/p>\r\n
 <\/p>\r\n","protected":false},"excerpt":{"rendered":"
I’ve been working with a non-profit, GreenStand, helping them with image processing tasks and related machine learning features for the technology part of their mission. GreenStand is (really) making the world a better place – walkin’ the walk.\u00a0 In short, they’re investing in people and local communities by encouraging them, financially and otherwise, to plant […]<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"closed","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[1],"tags":[],"class_list":["post-292","post","type-post","status-publish","format-standard","hentry","category-uncategorized"],"_links":{"self":[{"href":"https:\/\/www.dominionsw.com\/index.php?rest_route=\/wp\/v2\/posts\/292","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.dominionsw.com\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.dominionsw.com\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.dominionsw.com\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/www.dominionsw.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=292"}],"version-history":[{"count":102,"href":"https:\/\/www.dominionsw.com\/index.php?rest_route=\/wp\/v2\/posts\/292\/revisions"}],"predecessor-version":[{"id":607,"href":"https:\/\/www.dominionsw.com\/index.php?rest_route=\/wp\/v2\/posts\/292\/revisions\/607"}],"wp:attachment":[{"href":"https:\/\/www.dominionsw.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=292"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.dominionsw.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=292"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.dominionsw.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=292"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}