[et_pb_section fb_built=”1″ _builder_version=”3.22.7″][et_pb_row _builder_version=”3.22.7″][et_pb_column type=”4_4″ _builder_version=”3.22.7″][et_pb_text _builder_version=”3.22.7″]<\/p>\n
In my last my post I discussed this topic in the context of linear algebra and differential geometry. We also briefly discussed root finding method. Although the method works properly but you are probably looking for something faster! <\/p>\n
Here I would like to use a method from Revit API to improve the performance. <\/span><\/p>\n [\/et_pb_text][et_pb_text _builder_version=”3.22.7″ background_color=”#f2f2f2″ border_width_all=”1px” border_color_all=”#d3d3d3″ custom_margin=”-2px|||||” custom_padding=”11px|15px|0px|15px|false|true” saved_tabs=”all”]<\/p>\n [\/et_pb_text][et_pb_text _builder_version=”3.22.7″]<\/p>\n Above method compute the intersection of curve and face. If we can generate a planar face in the given plane in which intersect the curve then we can use the method to compute the Plane – Curve intersection. The question is in which position in the plane and in what size we should create the face? Let us assume that curve c<\/em><\/strong> intersect the plane Q<\/strong><\/em> at least in one point P<\/em><\/strong>. It is obvious that P<\/em><\/strong> lies in the curve c’<\/em><\/strong> the projection of Curve c<\/strong> <\/em>on the plane Q<\/em><\/strong>. Consequently if curve c<\/em><\/strong> intersect plane Q<\/em><\/strong> at a any point P<\/em><\/strong> then it also intersect curve c’<\/em><\/strong>, meaning the normal projection of curve c<\/em><\/strong> on plane Q<\/strong> <\/em>holds all possible intersection points. Now let us take any point on the curve c<\/strong><\/em> , let’s say point r<\/strong><\/em> and find it’s projection on plane Q<\/strong><\/em> which is r’<\/em><\/strong> a point on curve c’<\/strong><\/em>. From point r’<\/strong><\/em> on the plane Q<\/em><\/strong> we draw the circle s<\/strong><\/em> with the radius of l<\/em> <\/strong>the length of the curve c<\/strong><\/em>. We proof that circle s<\/strong><\/em> always contains curve c’<\/em><\/strong>. <\/p>\n [\/et_pb_text][et_pb_image src=”https:\/\/www.parametriczoo.com\/wp-content\/uploads\/2020\/04\/projection-of-curve-in-circle-s.jpg” _builder_version=”3.22.7″][\/et_pb_image][et_pb_text _builder_version=”3.22.7″]<\/p>\n We proof by contradiction. Let’s assume the Curve c’<\/strong><\/em> intersect the circle s<\/strong><\/em> in point x<\/strong><\/em>. This means the length of the line r’-x <\/b><\/em>equals to l<\/strong> <\/em>the length of the curve<\/b> c. <\/b><\/em>Hence we can write:<\/p>\n [\/et_pb_text][et_pb_image src=”https:\/\/www.parametriczoo.com\/wp-content\/uploads\/2020\/04\/image001.png” align=”center” _builder_version=”3.22.7″ width=”35%”][\/et_pb_image][et_pb_text _builder_version=”3.22.7″]<\/p>\n In the other hand the length of the line segment r-x’ <\/b><\/em>spanning between two points on the curve c’<\/strong> <\/em><\/b><\/em>is shorter than the curve :<\/i><\/b><\/p>\n [\/et_pb_text][et_pb_image src=”https:\/\/www.parametriczoo.com\/wp-content\/uploads\/2020\/04\/image003.png” align=”center” _builder_version=”3.22.7″ width=”37%”][\/et_pb_image][et_pb_text _builder_version=”3.22.7″]<\/p>\n From (1) and (2) we conclude :<\/p>\n [\/et_pb_text][et_pb_image src=”https:\/\/www.parametriczoo.com\/wp-content\/uploads\/2020\/04\/image005-1.png” align=”center” _builder_version=”3.22.7″ width=”37%”][\/et_pb_image][et_pb_text _builder_version=”3.22.7″]<\/p>\n Above statement is incorrect , as you can see in the previous post<\/a> the length of the normal projection of the curve c<\/strong><\/em> is less or equal to the length of the curve c<\/strong><\/em>. <\/p>\n Now that we know a circle with the radius of L(c)<\/strong><\/em> around any point of the c’<\/strong><\/em> includes the entire curve c’,<\/strong><\/em> we shall use this to construct a face of this 2l X 2l<\/strong><\/em> centered on a point (for example first point) of curve c’<\/strong><\/em> ti insure no intersection will falls within the face. This is nothing but the bounding rectangle of the circle s <\/strong><\/em>as shown in below<\/p>\n [\/et_pb_text][et_pb_image align=”center” src=”https:\/\/www.parametriczoo.com\/wp-content\/uploads\/2020\/04\/bounding-rectangle-1.jpg” _builder_version=”3.22.7″][\/et_pb_image][et_pb_text _builder_version=”3.22.7″]<\/p>\n In order to build a face out of the bounding rectangle we need to use a method from GeometryCreationUtilities<\/a> class. [\/et_pb_text][et_pb_text _builder_version=”3.22.7″ background_color=”#f2f2f2″ border_width_all=”1px” border_color_all=”#d3d3d3″ custom_margin=”-2px|||||” custom_padding=”11px|15px|0px|15px|false|true” saved_tabs=”all”]<\/p>\n [\/et_pb_text][et_pb_text _builder_version=”3.22.7″]<\/p>\n Using above method we can create an open shell with a single face . To do so , we require a list of CurveLoops. <\/a>To build a loft surface which looks like a rectangle we consider two lines either on the top and bottom or left and right side of the rectangle, CreateLoftGeometry()<\/a> then generate a surface spanning between the two lines. Below is the complete implementation. The method returns IntersectionResultArray<\/a> holding the information about the intersection events.<\/p>\n [\/et_pb_text][et_pb_text _builder_version=”3.22.7″ background_color=”#f2f2f2″ border_width_all=”1px” border_color_all=”#d3d3d3″ custom_margin=”-2px|||||” custom_padding=”11px|15px|0px|15px|false|true” saved_tabs=”all”]<\/p>\npublic<\/span> FaceIntersectionFaceResult<\/a> Intersect<\/span>(\n\tFace<\/a> face<\/span>,\n\tout<\/span> Curve<\/a> result<\/span>\n)<\/code><\/pre>\n
\nAll the methods in this class create solid geometry which makes it difficult to simply create a bounded surface (Face<\/a>) , Except one method! <\/p>\npublic<\/span> static<\/span> Solid<\/a> CreateLoftGeometry<\/span>(\n\tIList<\/a><CurveLoop<\/a>> profileLoops<\/span>,\n\tSolidOptions<\/a> solidOptions<\/span>\n)<\/code><\/pre>\n\/\/\/<\/span> <\/span><<\/span>summary<\/span>><\/span>\n\/\/\/<\/span> Compute the intersections of the plane with given curve<\/span>\n\/\/\/<\/span> <\/span><\/<\/span>summary<\/span>><\/span>\n\/\/\/<\/span> <\/span><<\/span>param<\/span> name<\/span>=<\/span>\"<\/span>p<\/span>\"<\/span>><\/span>Plane of the intersection<\/span><\/<\/span>param<\/span>><\/span>\n\/\/\/<\/span> <\/span><<\/span>param<\/span> name<\/span>=<\/span>\"<\/span>c<\/span>\"<\/span>><\/span>Curve to intersect with<\/span><\/<\/span>param<\/span>><\/span>\n\/\/\/<\/span> <\/span><<\/span>returns<\/span>><\/<\/span>returns<\/span>><\/span>\npublic<\/span> static<\/span> IntersectionResultArray<\/span> PlaneCurveIntersection<\/span>(<\/span>this<\/span> Plane<\/span> p<\/span>,<\/span>Curve<\/span> c<\/span>)<\/span>\n{<\/span>\n \/\/ Find the projection point<\/span>\n p<\/span>.<\/span>DistanceTo<\/span>(<\/span>c<\/span>.<\/span>GetEndPoint<\/span>(<\/span>0<\/span>),<\/span> out<\/span> XYZ<\/span> projection<\/span>);<\/span>\n \/\/ construct two parallel lines l1 and l2 <\/span>\n Line<\/span> l1<\/span> =<\/span> Line<\/span>.<\/span>CreateBound<\/span>(<\/span>projection<\/span> +<\/span> c<\/span>.<\/span>Length<\/span> *<\/span> (<\/span>p<\/span>.<\/span>XVec<\/span> +<\/span> p<\/span>.<\/span>YVec<\/span>),<\/span>\n projection<\/span> +<\/span> c<\/span>.<\/span>Length<\/span> *<\/span> (<\/span>p<\/span>.<\/span>XVec<\/span> -<\/span> p<\/span>.<\/span>YVec<\/span>));<\/span>\n Line<\/span> l2<\/span> =<\/span> Line<\/span>.<\/span>CreateBound<\/span>(<\/span>projection<\/span> +<\/span> c<\/span>.<\/span>Length<\/span> *<\/span> (<\/span>-<\/span>p<\/span>.<\/span>XVec<\/span> +<\/span> p<\/span>.<\/span>YVec<\/span>),<\/span>\n projection<\/span>