Update area_on_face.m for special of no intersection area.

scripts/area_on_face.m

@@ -24,7 +24,7 @@ for i = 1:NVertices
     allOneSide(i)       = all(northWest(:,i)) || all(~northWest(:,i));
 end
 
-% [nw(:,N),on(:,N)]   = check_points(pol,a,b,isv); % hva med p linje??
+% [nw(:,N),on(:,N)]   = check_points(pol,a,b,isv); % hva med på linje??
 % alloneside(N)       = all(nw(:,N)) || all(~nw(:,N));
 
 
@@ -65,5 +65,14 @@ end
 
 inside          = inpolygon(polygon(:,1),polygon(:,2),Polygon(:,1),Polygon(:,2));
 insidePoints    = polygon(inside,:);
-[~, a]          = convhull([intersections(1:j,:);insidePoints]);
+try
+    [~, a]          = convhull([intersections(1:j,:);insidePoints]);
+catch
+    % The exception should only occur if the intersection is empty. Displays the polygons
+    % for the user to inspect.
+    disp('Returned zero intersection area for coarse and fine polygons')
+    disp(Polygon)
+    disp(polygon)
+    a = 0;
+end
 end
\ No newline at end of file