You are adding lines to the input dataset at the end i.e. adding new cells to the dataset. However the cell-data which has to match the number of cells was never updated to have valid values for the newly added cells. Are you interested in the cell-data at all? Otherwise simply call output->GetCellData()->Initialize() to get rid of all cell-data arrays. You may want to refer to the VTK book or the VTK Users guide if you are unclear about how the cell-data, point-data etc. work.
Utkarsh 2009/9/23 young_jh123 <[email protected]>: > Thank you for your advice.I find it was the input problem,but I got a new > one.When I run my filter ,I got an > error:Cell array cell types with 1 components, has only 18 tuples but there > are 180 cells > . > > The below is how my RequestData() write.I am sorry to type long boring code > . I have tried to settle it by myself,but it seems a long time job. I write > the filter almost follow the source code vtkAppendFilter .I make my own code > bold-faced.I just want to add some lines to the input dataset. I know I > should add the cell data when I insert the line cells,but I don't know how > to add. > > code: > //get the output info object > vtkInformation *outInfo= outputVector->GetInformationObject(0); > //get the output > vtkUnstructuredGrid* output = > vtkUnstructuredGrid::SafeDownCast(outInfo->Get(vtkDataObject::DATA_OBJECT())); > vtkIdType numPts,numCells,ptOffset; > int tenth,count,abort = 0; > float decimal; > vtkPoints *newPts; > vtkPointData *pd; > vtkCellData *cd; > vtkIdList *ptIds,*newPtIds; > int idx; > vtkDataSet *ds; > vtkIdType ptId,cellId,newCellId; > vtkPointData *outputPD=output->GetPointData(); > vtkCellData *outPutCD = output ->GetCellData(); > vtkDebugMacro(<<"Appending data together"); > count = 0 ; > decimal = 0 ; > numPts = 0; > numCells = 0; > int numInputs = inputVector[0]->GetNumberOfInformationObjects(); > vtkDataSetAttributes::FieldList ptList(numInputs); > vtkDataSetAttributes::FieldList cellList(numInputs); > int firstPD = 1; > int firstCD = 1; > vtkInformation *inInfo = 0; > for(idx = 0;idx<numInputs;++idx) > { > inInfo = inputVector[0]->GetInformationObject(idx); > ds = 0; > if(inInfo) > { > ds = vtkDataSet::SafeDownCast(inInfo->Get(vtkDataObject::DATA_OBJECT())); > } > if(ds!=NULL) > { > { > if(ds->GetNumberOfPoints()<=0&&ds->GetNumberOfCells()<=0) > continue; > } > numPts +=ds->GetNumberOfPoints(); > numCells+=ds->GetNumberOfCells(); > pd = ds->GetPointData(); > if(firstPD) > { > ptList.InitializeFieldList(pd); > firstPD = 0; > } > else ptList.IntersectFieldList(pd); > cd = ds->GetCellData(); > if(firstCD) > { > cellList.InitializeFieldList(cd); > firstCD = 0; > } > else cellList.IntersectFieldList(cd); > }//if non-empty dataset > }//for all inputs > if(numPts<1) > { > vtkDebugMacro(<<"No data to append!"); > return 1; > } > // > > int numRays; > vector<WaveRay*> rayVector; > numRays = IcosahedronModel::GenerateRays(&rayVector,TF); > > //now can allocate > output->Allocate(numCells+numRays); > outputPD->CopyGlobalIdsOn(); > outputPD->CopyAllocate(ptList,numPts); > outPutCD->CopyGlobalIdsOn(); > outPutCD->CopyAllocate(cellList,numCells); > newPts = vtkPoints::New(); > newPts->SetNumberOfPoints(numPts+2*numRays); > ptIds = vtkIdList::New();ptIds->Allocate(VTK_CELL_SIZE); > newPtIds = vtkIdList::New();newPtIds->Allocate(VTK_CELL_SIZE); > //append each input dataset together > tenth = (numPts+numCells)/10+1; > ptOffset = 0; > int inputCount = 0; > for(idx = 0;idx <numInputs&&!abort;++idx) > { > inInfo = inputVector[0]->GetInformationObject(idx); > ds = 0; > if(inInfo) > { > ds = vtkDataSet::SafeDownCast(inInfo->Get(vtkDataObject::DATA_OBJECT())); > } > if(ds!=NULL&&(ds->GetNumberOfCells()>0||ds->GetNumberOfPoints()>0)) > { > numPts = ds->GetNumberOfPoints(); > numCells = ds->GetNumberOfCells(); > pd = ds->GetPointData(); > for(ptId = 0;ptId<numPts&&!abort;ptId++) > { > newPts->SetPoint(ptId+ptOffset,ds->GetPoint(ptId)); > outputPD->CopyData(ptList,pd,inputCount,ptId,ptId+ptOffset); > //Update progress > count++; > if(!(count%tenth)) > { > decimal+=0.1; > this->UpdateProgress(decimal); > abort= this->GetAbortExecute(); > } > } > cd = ds->GetCellData(); > for(cellId = 0;cellId<numCells&&!abort;cellId++) > { > ds->GetCellPoints(cellId,ptIds); > newPtIds->Reset(); > for(int i =0;i<ptIds->GetNumberOfIds();i++) > { > newPtIds->InsertId(i,ptIds->GetId(i)+ptOffset); > } > newCellId = output->InsertNextCell(ds->GetCellType(cellId),newPtIds); > outPutCD->CopyData(cellList,cd,inputCount,cellId,newCellId); > //Update progress > count++; > if(!(count%tenth)) > { > decimal+=0.1; > this->UpdateProgress(decimal); > abort = this->GetAbortExecute(); > } > } > ptOffset+=numPts; > ++inputCount; > } > } > //Add lines to the input dataset > vtkLine* line; > for(idx = 0;idx<numRays;++idx) > { > //newPts->InsertNextPoint((); > //newPts->InsertNextPoint(); > newPts->SetPoint(ptOffset++,((WaveRay*)rayVector[idx])->GetPosition()); > newPts->SetPoint(ptOffset++,((WaveRay*)rayVector[idx])->GetLinePosition(2)); > line = vtkLine::New(); > line->GetPointIds()->SetId(0,ptOffset-2); > line->GetPointIds()->SetId(1,ptOffset-1); > output->InsertNextCell(line->GetCellType(),line->GetPointIds()); > //I know I should add cell data attributes,but how? > > line->Delete(); > } > > //Update ourselves and release memory > > output->SetPoints(newPts); > newPts->Delete(); > ptIds->Delete(); > newPtIds->Delete(); > > > > > > > 2009-09-24 > ________________________________ > young_jh123 > ________________________________ > 发件人： Utkarsh Ayachit > 发送时间： 2009-09-22 22:25:12 > 收件人： young_jh123 > 抄送： Moreland, Kenneth; [email protected] > 主题： Re: [Paraview] Problems with user defined filter > So you have written a vtk filter, then I am assuming you have some main() in > which you are creating a VTK-pipeline to use this filter? Are you setting up > the inputs correctly? Try calling Update() directly on the filter, does > that result in calling of your RequestData() method? > Utkarsh > > On Tue, Sep 22, 2009 at 9:04 AM, young_jh123 <[email protected]> wrote: >> >> I want to write VTK filter first.Because I think it will be easier coding >> in windows than in linux,and if a VTK filter can work I know how to change a >> it to a paraview plugin filter. >> I'll try to write it as a paraview plugin directly.Thank you ! >> >> 2009-09-22 >> ________________________________ >> young_jh123 >> ________________________________ >> 发件人： Moreland, Kenneth >> 发送时间： 2009-09-22 06:59:07 >> 收件人： young_jh123; [email protected] >> 抄送： >> 主题： Re: [Paraview] Problems with user defined filter >> How are you trying to use your filter? Are you trying to define it in a >> ParaView plugin? If so, you should be able to leave out any of the _EXPORT >> declarations and it should work OK. >> >> -Ken >> >> >> On 9/21/09 9:02 AM, "young_jh123" <[email protected]> wrote: >> >> I write a filter class just like the vtkAppendFilter ,and change a little >> on the mothod RequestData ..My class also inherited from >> vtkUnstructuredGridAlgorithm .I did the job all on windows vs2008.The >> vtkAppendFilter.h defines the class with VTK_GRAPHICS_EXPORT which means >> __declspec( dllimport ) on my computer,but if I used that macro to define my >> own filter class it would met some link errors .If I use the macro >> VTK_EXPORT or did not use macro to define my filter class ,it can complie >> well,but it did not work.I tried the debug.I found it would not invoke my >> RequestData function at all,and the it would not invoke >> vtkUnstructuredGridAlgorithm'RequestData methord either.The program invoked >> the method RequestData form vtkAlgorithm class directly.I don't know why? >> what should I do if want to define my own filter? Thanks a lot! >> >> 2009-09-21 >> ________________________________ >> young_jh123 >> >> >> >> **** Kenneth Moreland >> *** Sandia National Laboratories >> *********** >> *** *** *** email: [email protected] >> ** *** ** phone: (505) 844-8919 >> *** web: http://www.cs.unm.edu/~kmorel >> >> >> _______________________________________________ >> Powered by www.kitware.com >> >> Visit other Kitware open-source projects at >> http://www.kitware.com/opensource/opensource.html >> >> Please keep messages on-topic and check the ParaView Wiki at: >> http://paraview.org/Wiki/ParaView >> >> Follow this link to subscribe/unsubscribe: >> http://www.paraview.org/mailman/listinfo/paraview >> > > _______________________________________________ Powered by www.kitware.com Visit other Kitware open-source projects at http://www.kitware.com/opensource/opensource.html Please keep messages on-topic and check the ParaView Wiki at: http://paraview.org/Wiki/ParaView Follow this link to subscribe/unsubscribe: http://www.paraview.org/mailman/listinfo/paraview
