Hi Samuel,
After tinkering around with workingSubset, I found that ellipses #4 and #18 were causing the error when integrated with others and the code works for workingSubset = [1:3, 5;17, 19]. I started further analyzing #3 and #4 together by plotting the ellipses using the following code:
figure, hold on;
axis equal;
xlim([-halfSide, halfSide]);
ylim([-halfSide, halfSide]);
theta = linspace(0, 2*pi);
for idx = workingSubset
a = s(idx).MajorAxisLength/2; % Semi-major axis
b = s(idx).MinorAxisLength/2; % Semi-minor axis
X0 = s(idx).Centroid(1)-halfSide; % Adjusted X centroid
Y0 = s(idx).Centroid(2)-halfSide; % Adjusted Y centroid
phi = deg2rad(s(idx).Orientation); % Orientation in radians
% Parametric equation of the ellipse
x = X0 + a*cos(theta)*cos(phi) - b*sin(theta)*sin(phi);
y = Y0 + a*cos(theta)*sin(phi) + b*sin(theta)*cos(phi);
plot(x, y, 'r', 'LineWidth', 2);
end
hold off;
On plotting the ellipses, I could see that the ellipses #3 and #4 are very close to each other. This could be a cause for the error of invalid geometry.
To resolve this:
- If possible, increase the distance between close geometries. If your geometries are very close to each other but not overlapping or intersecting, increasing the distance between them may help resolve the issue.
- You can also try simplyfying the geometries. This could involve reducing the number of points in your geometries, or simplifying the shapes themselves.
According to the documentation of decsg, specifically this part:
https://www.mathworks.com/help/pde/ug/decsg.html?s_tid=doc_ta#:~:text=Data%20Types%3A%20double-,Limitations,-In%20rare%20cases error can occur due to geometries being almost tangent.
Hope this help!
