Hi Jim and Bill,<br><br>Look very closed, in ITKsnap the slice 1 is slice 107 in my code, but when i get the value, i am getting the same slice. <br><br>I are begining think that i must get values intensity of voxel and not point. Look for cell data in vtk, i dont find one function that "GetScalarCell()", so how i get cell data in vtk? I must use "GetCell(100)->Derivatives()"? Cell = Voxel?<br>
<br>thanks<br><br><div class="gmail_quote">2010/9/12 Bill Lorensen <span dir="ltr"><<a href="mailto:bill.lorensen@gmail.com">bill.lorensen@gmail.com</a>></span><br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
vtk images are flipped in the y -direction. This is because vtk takes<br>
a "computer graphics" view of images while itk takes a "medical" view.<br>
<br>
In itk, the first pixel in memory is the first pixel on disk. In vtk,<br>
the first pixel in memory is the first pixel of the last y "line".<br>
<div><div></div><div class="h5"><br>
On Sun, Sep 12, 2010 at 11:37 AM, Jim Peterson <<a href="mailto:jimcp@cox.net">jimcp@cox.net</a>> wrote:<br>
> ok, I believe we have a registration problem between the vtk model<br>
> coordinates and the ITKSnap model coordinates. if I understand the data you<br>
> have given, we have some point data like:<br>
><br>
> tested point<br>
> ITKSnap value<br>
> ITK Snap Material<br>
> VTK Value<br>
> Vtk Material<br>
> pt1<br>
> 2318<br>
> dense bone<br>
> -1024<br>
> air (-1000 is air at STP calibration)<br>
> pt2<br>
> 2210<br>
> dense bone<br>
> 1023<br>
> soft bone<br>
> pt3<br>
> 2410<br>
> dense bone<br>
> 41<br>
> watery liquid (0 is water at STP calibration)<br>
> pt4<br>
> 2130<br>
> dense bone<br>
> 993<br>
> soft bone<br>
><br>
><br>
> We agree all of these values are reasonable for some location in a dental<br>
> scan. I believe the problem is, in spite of expectations and intentions, the<br>
> location tested are not from the same input data, maybe the coordinate<br>
> schemes in vtk and itk do not match. Can you find corresponding values in<br>
> the ITKSnap? perhaps by transposing a coordinate location? can we find<br>
> locations to "sample" in the vtk model that produce the dense bone values?<br>
><br>
> does that make sense?<br>
> Hope that helps,<br>
> Jim<br>
><br>
> Italo Ribeiro wrote:<br>
>><br>
>> Hi Jim,<br>
>><br>
>> The value not is same the entire tooth, the diferents values are the same<br>
>> point. In point (226, 100, 1) the values are:<br>
>><br>
>> ITKsnap : 2130<br>
>> Mycode: 993<br>
>><br>
>> 2010/9/12 Jim Peterson <<a href="mailto:jimcp@cox.net">jimcp@cox.net</a> <mailto:<a href="mailto:jimcp@cox.net">jimcp@cox.net</a>>><br>
>><br>
>> Hi Italo,<br>
>><br>
>> I am not an ITKSnap user myself, you are suggesting ITKSnap<br>
>> reports the same value for the entire tooth? I think you would<br>
>> have to iterate over the points in the tooth to find the maximum<br>
>> value in vtk to duplicate that if it is true. Can you display the<br>
>> point coordinates that have that value in ITKSnap? Then see if vtk<br>
>> has the same value at that point? I have trouble believing vtk has<br>
>> incorrect pixel data by point for imaging content.<br>
>><br>
>> Jim<br>
>><br>
>> Italo Ribeiro wrote:<br>
>><br>
>> Hi Jim,<br>
>><br>
>> I have dicom image with teeth, i am looking intensity values<br>
>> of 2 front teeth.<br>
>><br>
>> But some happen very strange, the tooth the left have greater<br>
>> intensity in ITKsnap 2518 and my code -1024. The tooth the<br>
>> right side have in ITKsnap 2402 and in my code 41. I do not<br>
>> know why the values are different.<br>
>><br>
>> In my code i get 255 points of heigh and width, but when i get<br>
>> dimension the image, vtk show 512 x 512, same dimesion show in<br>
>> ITKsnap. ITKsnap get dimension of voxels, so i must get<br>
>> information the intensity value voxel of volume and no pixel<br>
>> the image?<br>
>><br>
>> thanks.<br>
>><br>
>> 2010/9/11 Jim Peterson <<a href="mailto:jimcp@cox.net">jimcp@cox.net</a> <mailto:<a href="mailto:jimcp@cox.net">jimcp@cox.net</a>><br>
>> <mailto:<a href="mailto:jimcp@cox.net">jimcp@cox.net</a> <mailto:<a href="mailto:jimcp@cox.net">jimcp@cox.net</a>>>><br>
>><br>
>><br>
>> Personally, I am sorry, but there is not enough information in<br>
>> your description of the problem to be sure you have a<br>
>> legitimate<br>
>> complaint. does your dicom data consist of only one point?<br>
>> how are<br>
>> you certain you have selected the exact same location?<br>
>><br>
>> in my mainframe programming life we have a saying:<br>
>><br>
>> "when expected results do not match observed results, check the<br>
>> library concatenation"<br>
>><br>
>> .... in my opinion, chances are you are looking at<br>
>> different points.<br>
>><br>
>> sorry, Hoping to be helpful not critical.<br>
>> Jim<br>
>><br>
>> Italo Ribeiro wrote:<br>
>><br>
>> No more ideas?<br>
>><br>
>> 2010/9/10 Italo Ribeiro <<a href="mailto:italo.ribeiro@gmail.com">italo.ribeiro@gmail.com</a><br>
>> <mailto:<a href="mailto:italo.ribeiro@gmail.com">italo.ribeiro@gmail.com</a>><br>
>> <mailto:<a href="mailto:italo.ribeiro@gmail.com">italo.ribeiro@gmail.com</a><br>
>> <mailto:<a href="mailto:italo.ribeiro@gmail.com">italo.ribeiro@gmail.com</a>>><br>
>> <mailto:<a href="mailto:italo.ribeiro@gmail.com">italo.ribeiro@gmail.com</a><br>
>> <mailto:<a href="mailto:italo.ribeiro@gmail.com">italo.ribeiro@gmail.com</a>><br>
>> <mailto:<a href="mailto:italo.ribeiro@gmail.com">italo.ribeiro@gmail.com</a><br>
>> <mailto:<a href="mailto:italo.ribeiro@gmail.com">italo.ribeiro@gmail.com</a>>>>><br>
>><br>
>><br>
>> I am expecting MIP of image point. In my program the<br>
>> value<br>
>> for one<br>
>> point is 1023 and ITKsnap show 2210.<br>
>><br>
>> Why the value is different?<br>
>><br>
>> 2010/9/10 John Drescher <<a href="mailto:drescherjm@gmail.com">drescherjm@gmail.com</a><br>
>> <mailto:<a href="mailto:drescherjm@gmail.com">drescherjm@gmail.com</a>><br>
>> <mailto:<a href="mailto:drescherjm@gmail.com">drescherjm@gmail.com</a> <mailto:<a href="mailto:drescherjm@gmail.com">drescherjm@gmail.com</a>>><br>
>> <mailto:<a href="mailto:drescherjm@gmail.com">drescherjm@gmail.com</a><br>
>> <mailto:<a href="mailto:drescherjm@gmail.com">drescherjm@gmail.com</a>> <mailto:<a href="mailto:drescherjm@gmail.com">drescherjm@gmail.com</a><br>
>> <mailto:<a href="mailto:drescherjm@gmail.com">drescherjm@gmail.com</a>>>>><br>
>><br>
>><br>
>><br>
>> > Using GetScalarType the type returned is<br>
>> short, i was<br>
>> using<br>
>> int. But the<br>
>> > value maximum showed is 1024 and i see the<br>
>> maximum short<br>
>> value is 32767 and<br>
>> > the intensity value is not corret.<br>
>> ><br>
>><br>
>> Depending on how you read the dicom file the value<br>
>> should be in<br>
>> hounsfield units. I am confused on whether you<br>
>> mean you<br>
>> were<br>
>> expecting<br>
>> 32767 or hounsfield.<br>
>><br>
>> John<br>
>><br>
>><br>
>><br>
>><br>
>> -- Ítalo Mendes<br>
>> <a href="http://ogre.irados.org" target="_blank">ogre.irados.org</a> <<a href="http://ogre.irados.org" target="_blank">http://ogre.irados.org</a>><br>
>> <<a href="http://ogre.irados.org" target="_blank">http://ogre.irados.org</a>><br>
>> <<a href="http://ogre.irados.org" target="_blank">http://ogre.irados.org</a>><br>
>><br>
>><br>
>><br>
>><br>
>> -- Ítalo Mendes<br>
>><br>
>> <a href="http://ogre.irados.org" target="_blank">ogre.irados.org</a> <<a href="http://ogre.irados.org" target="_blank">http://ogre.irados.org</a>><br>
>> <<a href="http://ogre.irados.org" target="_blank">http://ogre.irados.org</a>> <<a href="http://ogre.irados.org" target="_blank">http://ogre.irados.org</a>><br>
>><br>
>> ------------------------------------------------------------------------<br>
>><br>
>> _______________________________________________<br>
>> Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a> <<a href="http://www.kitware.com" target="_blank">http://www.kitware.com</a>><br>
>> <<a href="http://www.kitware.com" target="_blank">http://www.kitware.com</a>><br>
>><br>
>><br>
>> Visit other Kitware open-source projects at<br>
>> <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
>><br>
>> Please keep messages on-topic and check the VTK FAQ at:<br>
>> <a href="http://www.vtk.org/Wiki/VTK_FAQ" target="_blank">http://www.vtk.org/Wiki/VTK_FAQ</a><br>
>><br>
>> Follow this link to subscribe/unsubscribe:<br>
>> <a href="http://www.vtk.org/mailman/listinfo/vtkusers" target="_blank">http://www.vtk.org/mailman/listinfo/vtkusers</a><br>
>><br>
>><br>
>><br>
>><br>
>> -- Ítalo Mendes<br>
>> <a href="http://ogre.irados.org" target="_blank">ogre.irados.org</a> <<a href="http://ogre.irados.org" target="_blank">http://ogre.irados.org</a>> <<a href="http://ogre.irados.org" target="_blank">http://ogre.irados.org</a>><br>
>><br>
>><br>
>><br>
>><br>
>><br>
>> --<br>
>> Ítalo Mendes<br>
>> <a href="http://ogre.irados.org" target="_blank">ogre.irados.org</a> <<a href="http://ogre.irados.org" target="_blank">http://ogre.irados.org</a>><br>
><br>
> _______________________________________________<br>
> Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
><br>
> Visit other Kitware open-source projects at<br>
> <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
><br>
> Please keep messages on-topic and check the VTK FAQ at:<br>
> <a href="http://www.vtk.org/Wiki/VTK_FAQ" target="_blank">http://www.vtk.org/Wiki/VTK_FAQ</a><br>
><br>
> Follow this link to subscribe/unsubscribe:<br>
> <a href="http://www.vtk.org/mailman/listinfo/vtkusers" target="_blank">http://www.vtk.org/mailman/listinfo/vtkusers</a><br>
><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>Ítalo Mendes<br><a href="http://ogre.irados.org">ogre.irados.org</a><br>