<div dir="ltr"><span style="font-size:12.8000001907349px">Yes, I can store it in PointData(), but splineData->GetDimensions(</span><span style="font-size:12.8000001907349px">numControlPointsPerAxis) will not return the correct number of control points. It will return {0,0,0} if we hold the control points in PointData(). Though we can also hold a simple 1*3 array in vtkFieldData just like we store the length of knot vectors.</span></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jun 30, 2015 at 2:10 PM, David Thompson <span dir="ltr"><<a href="mailto:david.thompson@kitware.com" target="_blank">david.thompson@kitware.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Lin,<br>
<span class=""><br>
> I didn't separate the homogenous coordinate with point coordinates but stored them all in a vtkDataArray with 4 components per tuple just as you suggest. In this way, vtkStructuredGrid cannot hold them as vtkPoints. The current API is like<br>
><br>
> void vtkNURBSPatchAdaptor::GetPatchShape(vtkUnstructuredGrid* outputSp, vtkDataArray* ctrlPts, int* ctrlPtsNum)<br>
> the outputSp holds the surface mesh, ctrlPts and ctrlPtsNum hold information and data of nurbs control points, and the vtkSturcturedGrid only holds knot vectors and length.<br>
<br>
</span>Actually, the vtkStructuredGrid can hold the control points as PointData() (as opposed to Point  coordinates). There is no constraint on the number of components per tuple. The idea is to access the control points and knot data like this:<br>
<br>
  vtkStructuredGrid* splineData;<br>
  vtkDataArray* controlPoints = splineData->GetPointData()->GetArrayByName("control polygon");<br>
  vtkDataArray* knotVector = splineData->GetFieldData()->GetArrayByName("knot vector");<br>
  int numControlPointsPerAxis[3];<br>
  splineData->GetDimensions(numControlPointsPerAxis);<br>
<br>
That way all of the spline data is grouped into one dataset; it's just that the splineData's vtkPoints would not hold the control points. We could also store additional information in the arrays of splineData->GetFieldData(). For example, flags indicating which axes (if any) were periodic could be stored in a separate field-data array since there is no constraint on the number of components or tuples of field-data arrays.<br>
<br>
Does that make sense?<br>
<span class="HOEnZb"><font color="#888888"><br>
        David<br>
<br>
</font></span></blockquote></div><br></div>