[Paraview] how to show the co-ordinates of the a point in paraview
hi, I have an unstructured grid of 758257 points. How do I show get paraview to show the co-ordinates of the points ..? ___ 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
Re: [Paraview] VTK_Edge for LIC
Zhanping, Will the options to modify the density, thickness, and colour of the streaklines be available to the user? It is definitely a good tool, but sometimes we find it better to visualise using the standard streamtracer as we can modify their appearance for clarity. -Bastil: If you use extract surface, you'll need to make sure you have a shear stress vector available on the surface (or create one before-hand), assuming you are doing aero flows. However, if you use extract block on the surfaces within the ensight volume data, you should be able to use the velocity vector as input for LIC too (I think it takes the data just off the surface for visualising, although I'm not sure). It isn't as nice as using shear stress, but works well enough. Regards, Adriano === Adriano Gagliardi MEng PhD Project Scientist Computational Aerodynamics Aircraft Research Association Ltd. Manton Lane Bedford Tel: 01234 32 4644 E-mail: agaglia...@ara.co.uk Url: www.ara.co.uk _ From: paraview-boun...@paraview.org [mailto:paraview-boun...@paraview.org] On Behalf Of Zhanping Liu Sent: 27 January 2010 22:24 To: bastil2...@yahoo.de Cc: ParaView Subject: Re: [Paraview] VTK_Edge for LIC In fact we are now integrating LIC to VTK. One of the features is that flow streaks are clearer than before in the resulting image. -Zhanping On Wed, Jan 27, 2010 at 5:02 PM, bastil2...@yahoo.de bastil2...@yahoo.de wrote: Dear group, I have compiled VTK_Edge against paraview 3.6.2. Everything worked fine except that make install did not install libVTKEdge.so to the install directory - I had to copy this manually. Is this a bug or a feature? Afterwards I can load libVTKEdge.so to paraview without problems. However, where are the new features? I am especially interested in LIC and I only found lic on 2D structured grids and LIC on image data in the filters menu. Both are greyed off with my Ensight-Reader multiblock unstructured grids. Does it only run on structured grid? Am I missing something? Regards BastiL ___ Powered by www.kitware.com http://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 -- Zhanping Liu, PhD Kitware, Inc. 28 Corporate Drive Clifton Park, NY 12065-8662 Phone: 518-371-3971 x 138 http://www.zhanpingliu.org --- This email contains information that is private and confidential and is intended only for the addressee. If you are not the intended recipient please delete it and notify us immediately by e-mailing the sender. Note: All email sent to or from this address may be accessed by someone other than the recipient, for system management and security reasons. Aircraft Research Association Ltd. Registered in England, Registration No 503668 Registered Office: Manton Lane, Bedford MK41 7PF England VAT No GB 196351245 ___ 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
Re: [Paraview] block merging using Paraview
Andy, Thanks for your help. I have attached two simple files I am trying to merge. But looking at the problem, I am afraid you will tell me that there is no solution. As you will see there is no duplicate points, thus Paraview cannot guess it. What is the best way to make them paraview-compliant? Add duplicate cells? Thanks. Julien 2010/1/27 Andy Bauer andy.ba...@kitware.com To be clear, I meant first use the group data set filter or merge block filter to put all of the cells in the same data set and then use the clean to grid filter to get rid of duplicate points. Points won't get merged if they are not within a certain distance from each other. So again the real problem isn't that all of the cells aren't in the same data set but the fact that the grid connectivity is not correct. If you are doing this in parallel this can also have an affect on appearance due to not having the correct connectivity between cells on different partitions/processors. If you can give me a more detailed description of the problem and send the grid I can take a closer look at it. Andy On Wed, Jan 27, 2010 at 10:04 AM, Julien Bodart julien.bod...@gmail.comwrote: Well, actually I think I tried all the filters, included the clean to grid one. In my case there is actually no duplicate point between files: for example: file-1, z-coord: 0.0,0.1,0.2 file-2, z-coord: 0.3,0.4,0.5 And the grid is structured. Maybe in this case there is no way to make a single dataset? 2010/1/27 Andy Bauer andy.ba...@kitware.com Try using the clean to grid filter. Although all of the cells are in the same data set (probably an unstructured grid), the filter does not realize that there are duplicate points and because of this the grid connectivity is probably not what you're expecting. Andy On Wed, Jan 27, 2010 at 8:03 AM, Julien Bodart julien.bod...@gmail.comwrote: Hi, I am new to Paraview, thanks to the new Netcdf Reader of 3.6 release. (Thank you very much for that) Therefore my problem is probably really trivial but I can't get it resolved... I am reading 2 or more netcdf file describing several part of the grid. The complete domain is a cube(rectilinear grid), sliced in a given number of files. There is no ghost cells so there is no obvious match between files. When I try to merge the different block (using the group dataset filter or merge block), I end up with a split domain, leading to a visual gap between each sub-domain, whether I am plotting contour , vertical plane or whatever. Is there a way to really merge those subdomain to end up with a single domain. Thanks in advance! Julien ___ 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 test.1.nc Description: Cdf file test.2.nc Description: Cdf file ___ 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
[Paraview] Problems writing an .vtu file
Hello everybody, I've been trying to write a c++ function to convert my data into .vtu files. Unfortunately, paraview isn't able to read my .vtu files. I get the following error-message: ERROR: In /home/zam/domini/paraview/ParaView3src/VTK/Filtering/vtkExecutive.cxx, line 757 vtkCompositeDataPipeline (0x8cf34b0): Algorithm vtkXMLUnstructuredGridReader(0x8d915d0) returned failure for request: vtkInformation (0x90c34e8) Debug: Off Modified Time: 383066 Reference Count: 1 Registered Events: (none) Request: REQUEST_INFORMATION FORWARD_DIRECTION: 0 ALGORITHM_AFTER_FORWARD: 1 I'm obviously doing something wrong... can anyone help? I really can't find the error. Below, you find a sample .vtu file. Thanks in advance!! Sonja ?xml version=1.0? VTKFile type=UnstructuredGrid version=0.1 byte_order=LittleEndian UnstructuredGrid Piece NumberOfPoints=96 NumberOfCells=12 PointData Scalars=process_id DataArray type=Float32 name=process_id format=ascii 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 4 4 4 4 4 4 4 4 5 5 5 5 5 5 5 5 6 6 6 6 6 6 6 6 0 0 0 0 0 0 0 0 7 7 7 7 7 7 7 7 8 8 8 8 8 8 8 8 10 10 10 10 10 10 10 10 9 9 9 9 9 9 9 9 11 11 11 11 11 11 11 11 /DataArray /PointData CellData /CellData Points DataArray type=Float32 NumberOfComponents=3 format=ascii 0 -37.5 -225 37.5 -37.5 -225 0 0 -225 37.5 0 -225 0 -37.5 -75 37.5 -37.5 -75 0 0 -75 37.5 0 -75 -37.5 0 -225 0 0 -225 -37.5 37.5 -225 0 37.5 -225 -37.5 0 -75 0 0 -75 -37.5 37.5 -75 0 37.5 -75 0 0 -225 37.5 0 -225 0 37.5 -225 37.5 37.5 -225 0 0 -75 37.5 0 -75 0 37.5 -75 37.5 37.5 -75 -37.5 -37.5 -75 0 -37.5 -75 -37.5 0 -75 0 0 -75 -37.5 -37.5 75 0 -37.5 75 -37.5 0 75 0 0 75 0 -37.5 -75 37.5 -37.5 -75 0 0 -75 37.5 0 -75 0 -37.5 75 37.5 -37.5 75 0 0 75 37.5 0 75 -37.5 0 -75 0 0 -75 -37.5 37.5 -75 0 37.5 -75 -37.5 0 75 0 0 75 -37.5 37.5 75 0 37.5 75 -37.5 -37.5 -225 0 -37.5 -225 -37.5 0 -225 0 0 -225 -37.5 -37.5 -75 0 -37.5 -75 -37.5 0 -75 0 0 -75 0 0 -75 37.5 0 -75 0 37.5 -75 37.5 37.5 -75 0 0 75 37.5 0 75 0 37.5 75 37.5 37.5 75 -37.5 -37.5 75 0 -37.5 75 -37.5 0 75 0 0 75 -37.5 -37.5 225 0 -37.5 225 -37.5 0 225 0 0 225 -37.5 0 75 0 0 75 -37.5 37.5 75 0 37.5 75 -37.5 0 225 0 0 225 -37.5 37.5 225 0 37.5 225 0 -37.5 75 37.5 -37.5 75 0 0 75 37.5 0 75 0 -37.5 225 37.5 -37.5 225 0 0 225 37.5 0 225 0 0 75 37.5 0 75 0 37.5 75 37.5 37.5 75 0 0 225 37.5 0 225 0 37.5 225 37.5 37.5 225 /DataArray /Points Cells DataArray type=Int32 Name=connectivity format=ascii 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 /DataArray DataArray type=Int32 Name=offsets format=ascii 8 16 24 32 40 48 56 64 72 80 88 96 /DataArray DataArray type=UInt8 Name=types format=ascii 11 11 11 11 11 11 11 11 11 11 11 11 /DataArray /Cells /Piece /UnstructuredGrid /VTKFile Forschungszentrum Juelich GmbH 52425 Juelich Sitz der Gesellschaft: Juelich Eingetragen im Handelsregister des Amtsgerichts Dueren Nr. HR B 3498 Vorsitzende des Aufsichtsrats: MinDir'in Baerbel Brumme-Bothe Geschaeftsfuehrung: Prof. Dr. Achim Bachem (Vorsitzender), Dr. Ulrich Krafft (stellv. Vorsitzender), Prof. Dr.-Ing. Harald Bolt, Prof. Dr. Sebastian M. Schmidt ___ 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
Re: [Paraview] X Error: BadLength on OpenSUSE11.2
All right. I solved the issue. My graphic card drivers where not properly installed. Re-installation of these have solved the problem and now I have direct rendering again and can use volume. Thanks for your hint. Robert ___ 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
Re: [Paraview] Problems writing an .vtu file
s.habbinga wrote: Hello everybody, I've been trying to write a c++ function to convert my data into .vtu files. Unfortunately, paraview isn't able to read my .vtu files. I get the following error-message: there is a typographic error in your file. Replace the line DataArray type=Float32 name=process_id format=ascii by DataArray type=Float32 Name=process_id format=ascii ___ 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
[Paraview] compare two datasets in paraview
hi me again, how do I compare the co-ordinates of two datasets in paraview.. I have two dataset..of the same set of points but their cordinates differ slightly. I need the number of these co-ordinates ___ 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
Re: [Paraview] pvserver : cannot connect to X server
Well it does not solve my problem ... I never had this behaviour before / connecting with X was not necessary. This problem started this week / last week cvs was fine. Any idea why it's causing that ? Thanks PO On 2010-01-27, at 2:43 PM, Pierre-Olivier Dallaire wrote: I never used the the -X option when connecting with ssh on the server in order to initialize pvserver. I think now the solution is to use the -X but start pvserver with option -dc PO On 2010-01-27, at 2:35 PM, Utkarsh Ayachit wrote: Are you sure the display is accessible on the server machine? What happens when you open something like gvim? Utkarsh On Wed, Jan 27, 2010 at 10:04 AM, Pierre-Olivier Dallaire pierre-olivier.dalla...@videotron.ca wrote: Hi, when trying to start pvserver on the my ParaView server linux machine, I'm getting this new error with the latest cvs version, error that I never had before : pvserver : cannot connect to X server Thanks! PO ___ 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 ___ 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
Re: [Paraview] how to show the co-ordinates of the a point in paraview
Hey Kalpana, You can get ParaView to display scalar or vector quantities associated with points while they are selected. If you select some points, then under the View menu open the Selection Inspector, at the bottom of this panel there is a section called Display Style. Click on the Point Label tab and check Visible. You will see a combo-box called Label Mode, and in this list are all of the scalars and vectors associated with your point data. The point coordinates are not listed, but you can can create a new vector that contains a copy of the coordinates with the Calculator filter. In the Calculator Properties panel there is a combo box/button labeled Vectors, and one of the entries should be coords (or you can just type coords in the equation line edit. You can rename the result to something more meaningful like Coordinates. Now when you select points from this new data, the new Coordinates vector should show up under the Label Mode of the Point Label section of the Selection Inspector. Let us know if you have any problems, -Eric -- Eric E Monson Duke Visualization Technology Group On Jan 28, 2010, at 4:31 AM, Kalpana Kanthasamy wrote: hi, I have an unstructured grid of 758257 points. How do I show get paraview to show the co-ordinates of the points ..? ___ 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
[Paraview] how to load whole trajectory data into paraview
Dear friends, I have a trajectory data of granular particles from a simulation. I can visualize easily single configuration of particles using .vtk file. But when I load the whole trajectory data having 100 frames (or particle's configurations) then Paraview shows some error Please guide me how can I do that. Thanks. Ashish ___ 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
Re: [Paraview] compare two datasets in paraview
Hey again, There may be other ways, too, but you can do this with the Python Programmable Filter (PPF). To have two inputs to the PPF you need to select both data sets in the pipeline (on mac it's command-click, others may be control-click...) and then apply the filter. Here is a sample script which subtracts the coordinates of the two data sets from each other and makes a new vector data attribute out of that difference: # === from paraview.vtk import vtkFloatArray in0 = self.GetInputDataObject(0,0) in1 = self.GetInputDataObject(0,1) pti0 = in0.GetPoints() pti1 = in1.GetPoints() numPts = in0.GetNumberOfPoints() vec = vtkFloatArray() vec.SetName('coordDiff') vec.SetNumberOfComponents(3) vec.SetNumberOfTuples(numPts) for ii in range(numPts): # compute velocity and put in vec tmpDiff = [0.0, 0.0, 0.0] pos1 = pti1.GetPoint(ii) pos0 = pti0.GetPoint(ii) for kk in range(len(tmpDiff)): tmpDiff[kk] = pos1[kk] - pos0[kk] vec.SetTuple3(ii,tmpDiff[0],tmpDiff[1],tmpDiff[2]) out1 = self.GetOutputDataObject(0) out1.ShallowCopy(in0) out1.GetPointData().AddArray(vec) # === This doesn't do any checking for matching IDs of the points or whether both sets have the same number of points (I have an extended script which searches for points with matching ID if you need help doing that), but if they really match it will work. (Others can let you know if there are better ways of doing this or if I've done anything non-standard in my script.) Good luck, -Eric -- Eric E Monson Duke Visualization Technology Group On Jan 28, 2010, at 7:10 AM, Kalpana Kanthasamy wrote: hi me again, how do I compare the co-ordinates of two datasets in paraview.. I have two dataset..of the same set of points but their cordinates differ slightly. I need the number of these co-ordinates ___ 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
Re: [Paraview] how to show the co-ordinates of the a point in paraview
If you just want to see a list of point coordinates, you could create a spreadsheet view and view your grid in that. -Ken On 1/28/10 7:10 AM, Eric E. Monson emon...@cs.duke.edu wrote: Hey Kalpana, You can get ParaView to display scalar or vector quantities associated with points while they are selected. If you select some points, then under the View menu open the Selection Inspector, at the bottom of this panel there is a section called Display Style. Click on the Point Label tab and check Visible. You will see a combo-box called Label Mode, and in this list are all of the scalars and vectors associated with your point data. The point coordinates are not listed, but you can can create a new vector that contains a copy of the coordinates with the Calculator filter. In the Calculator Properties panel there is a combo box/button labeled Vectors, and one of the entries should be coords (or you can just type coords in the equation line edit. You can rename the result to something more meaningful like Coordinates. Now when you select points from this new data, the new Coordinates vector should show up under the Label Mode of the Point Label section of the Selection Inspector. Let us know if you have any problems, -Eric -- Eric E Monson Duke Visualization Technology Group On Jan 28, 2010, at 4:31 AM, Kalpana Kanthasamy wrote: hi, I have an unstructured grid of 758257 points. How do I show get paraview to show the co-ordinates of the points ..? ___ 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 Kenneth Moreland *** Sandia National Laboratories *** *** *** *** email: kmo...@sandia.gov ** *** ** 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
Re: [Paraview] pvserver : cannot connect to X server
I've just run into the same problem (CVS updated at earlier in the week) and I cannot run pvserver without a X display. It means I can't run it on my cluster (not even when using OSMesa and the --use-offscreen-rendering option) Any ideas? Thanks, Paul On 28 January 2010 12:57, Pierre-Olivier Dallaire pierre-olivier.dalla...@videotron.ca wrote: Well it does not solve my problem ... I never had this behaviour before / connecting with X was not necessary. This problem started this week / last week cvs was fine. Any idea why it's causing that ? Thanks PO On 2010-01-27, at 2:43 PM, Pierre-Olivier Dallaire wrote: I never used the the -X option when connecting with ssh on the server in order to initialize pvserver. I think now the solution is to use the -X but start pvserver with option -dc PO On 2010-01-27, at 2:35 PM, Utkarsh Ayachit wrote: Are you sure the display is accessible on the server machine? What happens when you open something like gvim? Utkarsh On Wed, Jan 27, 2010 at 10:04 AM, Pierre-Olivier Dallaire pierre-olivier.dalla...@videotron.ca wrote: Hi, when trying to start pvserver on the my ParaView server linux machine, I'm getting this new error with the latest cvs version, error that I never had before : pvserver : cannot connect to X server Thanks! PO ___ 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 ___ 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
Re: [Paraview] block merging using Paraview
Please keep the discussions on the mailing list so that others can contribute and benefit from the information. I think the best solution is to create a custom filter that stitches the two rectilinear grids together in the way you want. Most other solutions will convert your rectilinear grids to unstructured grids and thus use a lot more memory. Usually when you split a grid into multiple grids you'll get visual artifacts. Andy On Thu, Jan 28, 2010 at 12:12 PM, Julien Bodart julien.bod...@gmail.comwrote: I am more working on the first option. Actually those two files share a single uniform grid, so the space between the two blocks is exactly the cell size. (so in this simple case it is quite huge). The reverse question works as well: what happens if I split one single data file (with a rectilinear grid) in two files (with no connections), just two subsets and two subgrids, how to load them and see a single grid. In my case loading the two files gives you all the grid nodes and then all the needed information, but it looks like you have two different blocks. In a finite volume point of view, maybe it means adding one cell between the two blocks, but if the information is on the nodes like here, nothing more is needed. Your first solution works for me. I mean what you get using these simple files is exactly what I want. The problem there is the memory consumption (for example on my pc I can open two regular files using just 22% of the memory, the MergeBlock filters goes to 66% and then every action make paraview crash). It looks the brute force for something I was expecting to be just tricky. (no extra nodes needed) I hope it is more descriptive. 2010/1/28 Andy Bauer andy.ba...@kitware.com I'm still not completely understanding what you want done. Since the two data sets have significant space between them do you want to fill in that space with extra cells or do you want to translate one of the data sets to get rid of the space between? If you want to add in cells then it might be easiest to merge the data sets (GroupDataSet filter and then MergeBlocks filter like below) and then use the vtkDelauney3D filter to create a new grid that keeps the point data. If you want to translate one of the grids to get rid of the space between them you'll need to use something like the CleanToGrid filter to get an unstructured grid from your rectilinear grid. Then you can use the Transform filter to move it to the proper position (in this case -0.100148999 in the z direction for test.2.nc since the vtkCleantoGridFilter is quite finicky). Then use GroupDataSets filter to merge the data sets into a multiblock (right click on the GroupDataSets filter to select one of the test.*.nc files and the the Transform filter). Finally, you'll need to use the vtkCleantoGrid filter again to get the proper connectivity. You'll want to check what it does with the point data at the merged points though. If neither of these is what you want then you'll need to be very descriptive of what you hope to accomplish in order for me to help you any more. Andy On Thu, Jan 28, 2010 at 5:51 AM, Julien Bodart julien.bod...@gmail.comwrote: Andy, Thanks for your help. I have attached two simple files I am trying to merge. But looking at the problem, I am afraid you will tell me that there is no solution. As you will see there is no duplicate points, thus Paraview cannot guess it. What is the best way to make them paraview-compliant? Add duplicate cells? Thanks. Julien 2010/1/27 Andy Bauer andy.ba...@kitware.com To be clear, I meant first use the group data set filter or merge block filter to put all of the cells in the same data set and then use the clean to grid filter to get rid of duplicate points. Points won't get merged if they are not within a certain distance from each other. So again the real problem isn't that all of the cells aren't in the same data set but the fact that the grid connectivity is not correct. If you are doing this in parallel this can also have an affect on appearance due to not having the correct connectivity between cells on different partitions/processors. If you can give me a more detailed description of the problem and send the grid I can take a closer look at it. Andy On Wed, Jan 27, 2010 at 10:04 AM, Julien Bodart julien.bod...@gmail.com wrote: Well, actually I think I tried all the filters, included the clean to grid one. In my case there is actually no duplicate point between files: for example: file-1, z-coord: 0.0,0.1,0.2 file-2, z-coord: 0.3,0.4,0.5 And the grid is structured. Maybe in this case there is no way to make a single dataset? 2010/1/27 Andy Bauer andy.ba...@kitware.com Try using the clean to grid filter. Although all of the cells are in the same data set (probably an unstructured grid), the filter does not realize that there are duplicate points and because of this the
Re: [Paraview] block merging using Paraview
Thanks for your help and sorry for the mailing list. I will investigate the custom filters to get rid of the visual artefact. I still have your solution. I was just afraid to miss an obvious way to do it. Thanks again for your time. Julien 2010/1/28, Andy Bauer andy.ba...@kitware.com: Please keep the discussions on the mailing list so that others can contribute and benefit from the information. I think the best solution is to create a custom filter that stitches the two rectilinear grids together in the way you want. Most other solutions will convert your rectilinear grids to unstructured grids and thus use a lot more memory. Usually when you split a grid into multiple grids you'll get visual artifacts. Andy On Thu, Jan 28, 2010 at 12:12 PM, Julien Bodart julien.bod...@gmail.comwrote: I am more working on the first option. Actually those two files share a single uniform grid, so the space between the two blocks is exactly the cell size. (so in this simple case it is quite huge). The reverse question works as well: what happens if I split one single data file (with a rectilinear grid) in two files (with no connections), just two subsets and two subgrids, how to load them and see a single grid. In my case loading the two files gives you all the grid nodes and then all the needed information, but it looks like you have two different blocks. In a finite volume point of view, maybe it means adding one cell between the two blocks, but if the information is on the nodes like here, nothing more is needed. Your first solution works for me. I mean what you get using these simple files is exactly what I want. The problem there is the memory consumption (for example on my pc I can open two regular files using just 22% of the memory, the MergeBlock filters goes to 66% and then every action make paraview crash). It looks the brute force for something I was expecting to be just tricky. (no extra nodes needed) I hope it is more descriptive. 2010/1/28 Andy Bauer andy.ba...@kitware.com I'm still not completely understanding what you want done. Since the two data sets have significant space between them do you want to fill in that space with extra cells or do you want to translate one of the data sets to get rid of the space between? If you want to add in cells then it might be easiest to merge the data sets (GroupDataSet filter and then MergeBlocks filter like below) and then use the vtkDelauney3D filter to create a new grid that keeps the point data. If you want to translate one of the grids to get rid of the space between them you'll need to use something like the CleanToGrid filter to get an unstructured grid from your rectilinear grid. Then you can use the Transform filter to move it to the proper position (in this case -0.100148999 in the z direction for test.2.nc since the vtkCleantoGridFilter is quite finicky). Then use GroupDataSets filter to merge the data sets into a multiblock (right click on the GroupDataSets filter to select one of the test.*.nc files and the the Transform filter). Finally, you'll need to use the vtkCleantoGrid filter again to get the proper connectivity. You'll want to check what it does with the point data at the merged points though. If neither of these is what you want then you'll need to be very descriptive of what you hope to accomplish in order for me to help you any more. Andy On Thu, Jan 28, 2010 at 5:51 AM, Julien Bodart julien.bod...@gmail.comwrote: Andy, Thanks for your help. I have attached two simple files I am trying to merge. But looking at the problem, I am afraid you will tell me that there is no solution. As you will see there is no duplicate points, thus Paraview cannot guess it. What is the best way to make them paraview-compliant? Add duplicate cells? Thanks. Julien 2010/1/27 Andy Bauer andy.ba...@kitware.com To be clear, I meant first use the group data set filter or merge block filter to put all of the cells in the same data set and then use the clean to grid filter to get rid of duplicate points. Points won't get merged if they are not within a certain distance from each other. So again the real problem isn't that all of the cells aren't in the same data set but the fact that the grid connectivity is not correct. If you are doing this in parallel this can also have an affect on appearance due to not having the correct connectivity between cells on different partitions/processors. If you can give me a more detailed description of the problem and send the grid I can take a closer look at it. Andy On Wed, Jan 27, 2010 at 10:04 AM, Julien Bodart julien.bod...@gmail.com wrote: Well, actually I think I tried all the filters, included the clean to grid one. In my case there is actually no duplicate point between files: for example: file-1, z-coord: 0.0,0.1,0.2 file-2, z-coord: 0.3,0.4,0.5 And the grid is structured. Maybe
[Paraview] Kitware blog
Kitware launched http://www.kitware.com/news/home/browse/267 its first developer blog http://www.kitware.com/blog/ today with contributions from Company technical and business leaders. Following the blog is a great way to obtain early access to evolving technology and provide commentary feedback. Blog postings may be followed by topic area, or on an individual basis. You may wish to subscribe to blog notifications via emailhttp://www.kitware.com/blog/home/subscribe or RSS http://www.kitware.com/blog/home/rss feed or you may subscribe to posts written by specific members of our RD team. This can be done through the blog directly or through a specific developer's Team Pagehttp://kitware.com/company/ourteam.html. As always feedback and comments are welcome. We now return you to your regularly scheduled program. David E DeMarle Kitware, Inc. RD Engineer 28 Corporate Drive Clifton Park, NY 12065-8662 Phone: 518-371-3971 x109 ___ 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