Hi Jim,<br><br>I get the solution, the problem is that i am dont select position correctly with mouse. I stop use integer and use now float for get position the mouse cilck and i now i get the correct values.<br><br>One observation, the positions the ITKsnap and vtk are be differents. The Y axi is inverted.<br>
<br>Thank you for everything.<br><br><div class="gmail_quote">2010/9/13 Jim Peterson <span dir="ltr"><<a href="mailto:jimcp@cox.net">jimcp@cox.net</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;">
Italo,<br>
I am not sure what I am supposed to "look very close" at, nor what actual problem we are addressing. What you have raised is a technical need to access the scalar data for an image in vtk, could easily be a classroom assignment, and not a problem I have any great interest in.<br>
Have you implemented a point picker in vtk to interact with the 3d model in vtk to access the scalar values for selected points? I believe you will find each structured grid "point" to be a 3d space location derived from the row and column location in a dicom image file, the combination of row, column and file number times the sample spacing results in a location is space for the point. the radio intensity value results in a scalar. at a sufficient scaling a point scalar is a voxel scalar is a pixel scalar, depending on whether we are talking the structured grid in 3d space, a surface element in the space or a row and column in a single dicom image file.<br>
<br>
The image of the dicom data loaded is generated in itk and vtk by similar "surface generating" algorithms that construct polygon defined surfaces connecting points with similar scalar values. So if ITKSnap and vtk both generate similar images, they are both finding similar scalar data. <br>
I am afraid I am at the end if my interest in this problem,<br>
Jim<br>
<br>
Italo Ribeiro wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div class="im">
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>
2010/9/12 Bill Lorensen <<a href="mailto:bill.lorensen@gmail.com" target="_blank">bill.lorensen@gmail.com</a> <mailto:<a href="mailto:bill.lorensen@gmail.com" target="_blank">bill.lorensen@gmail.com</a>>><div class="im">
<br>
<br>
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>
<br>
On Sun, Sep 12, 2010 at 11:37 AM, Jim Peterson <<a href="mailto:jimcp@cox.net" target="_blank">jimcp@cox.net</a><br></div><div><div></div><div class="h5">
<mailto:<a href="mailto:jimcp@cox.net" target="_blank">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<br>
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<br>
a dental<br>
> scan. I believe the problem is, in spite of expectations and<br>
intentions, the<br>
> location tested are not from the same input data, maybe the<br>
coordinate<br>
> schemes in vtk and itk do not match. Can you find corresponding<br>
values in<br>
> the ITKSnap? perhaps by transposing a coordinate location? can<br>
we find<br>
> locations to "sample" in the vtk model that produce the dense<br>
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<br>
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" target="_blank">jimcp@cox.net</a> <mailto:<a href="mailto:jimcp@cox.net" target="_blank">jimcp@cox.net</a>><br></div></div><div><div></div><div class="h5">
<mailto:<a href="mailto:jimcp@cox.net" target="_blank">jimcp@cox.net</a> <mailto:<a href="mailto:jimcp@cox.net" target="_blank">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<br>
display the<br>
>> point coordinates that have that value in ITKSnap? Then see<br>
if vtk<br>
>> has the same value at that point? I have trouble believing<br>
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<br>
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<br>
i get<br>
>> dimension the image, vtk show 512 x 512, same dimesion<br>
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" target="_blank">jimcp@cox.net</a><br>
<mailto:<a href="mailto:jimcp@cox.net" target="_blank">jimcp@cox.net</a>> <mailto:<a href="mailto:jimcp@cox.net" target="_blank">jimcp@cox.net</a> <mailto:<a href="mailto:jimcp@cox.net" target="_blank">jimcp@cox.net</a>>><br>
>> <mailto:<a href="mailto:jimcp@cox.net" target="_blank">jimcp@cox.net</a> <mailto:<a href="mailto:jimcp@cox.net" target="_blank">jimcp@cox.net</a>><br>
<mailto:<a href="mailto:jimcp@cox.net" target="_blank">jimcp@cox.net</a> <mailto:<a href="mailto:jimcp@cox.net" target="_blank">jimcp@cox.net</a>>>>><br>
>><br>
>><br>
>> Personally, I am sorry, but there is not enough<br>
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<br>
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,<br>
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" target="_blank">italo.ribeiro@gmail.com</a><br>
<mailto:<a href="mailto:italo.ribeiro@gmail.com" target="_blank">italo.ribeiro@gmail.com</a>><br>
>> <mailto:<a href="mailto:italo.ribeiro@gmail.com" target="_blank">italo.ribeiro@gmail.com</a><br>
<mailto:<a href="mailto:italo.ribeiro@gmail.com" target="_blank">italo.ribeiro@gmail.com</a>>><br>
>> <mailto:<a href="mailto:italo.ribeiro@gmail.com" target="_blank">italo.ribeiro@gmail.com</a><br>
<mailto:<a href="mailto:italo.ribeiro@gmail.com" target="_blank">italo.ribeiro@gmail.com</a>><br>
>> <mailto:<a href="mailto:italo.ribeiro@gmail.com" target="_blank">italo.ribeiro@gmail.com</a><br>
<mailto:<a href="mailto:italo.ribeiro@gmail.com" target="_blank">italo.ribeiro@gmail.com</a>>>><br>
>> <mailto:<a href="mailto:italo.ribeiro@gmail.com" target="_blank">italo.ribeiro@gmail.com</a><br>
<mailto:<a href="mailto:italo.ribeiro@gmail.com" target="_blank">italo.ribeiro@gmail.com</a>><br>
>> <mailto:<a href="mailto:italo.ribeiro@gmail.com" target="_blank">italo.ribeiro@gmail.com</a><br>
<mailto:<a href="mailto:italo.ribeiro@gmail.com" target="_blank">italo.ribeiro@gmail.com</a>>><br>
>> <mailto:<a href="mailto:italo.ribeiro@gmail.com" target="_blank">italo.ribeiro@gmail.com</a><br>
<mailto:<a href="mailto:italo.ribeiro@gmail.com" target="_blank">italo.ribeiro@gmail.com</a>><br>
>> <mailto:<a href="mailto:italo.ribeiro@gmail.com" target="_blank">italo.ribeiro@gmail.com</a><br>
<mailto:<a href="mailto:italo.ribeiro@gmail.com" target="_blank">italo.ribeiro@gmail.com</a>>>>>><br>
>><br>
>><br>
>> I am expecting MIP of image point. In my<br>
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" target="_blank">drescherjm@gmail.com</a><br>
<mailto:<a href="mailto:drescherjm@gmail.com" target="_blank">drescherjm@gmail.com</a>><br>
>> <mailto:<a href="mailto:drescherjm@gmail.com" target="_blank">drescherjm@gmail.com</a> <mailto:<a href="mailto:drescherjm@gmail.com" target="_blank">drescherjm@gmail.com</a>>><br>
>> <mailto:<a href="mailto:drescherjm@gmail.com" target="_blank">drescherjm@gmail.com</a><br>
<mailto:<a href="mailto:drescherjm@gmail.com" target="_blank">drescherjm@gmail.com</a>> <mailto:<a href="mailto:drescherjm@gmail.com" target="_blank">drescherjm@gmail.com</a><br>
<mailto:<a href="mailto:drescherjm@gmail.com" target="_blank">drescherjm@gmail.com</a>>>><br>
>> <mailto:<a href="mailto:drescherjm@gmail.com" target="_blank">drescherjm@gmail.com</a><br>
<mailto:<a href="mailto:drescherjm@gmail.com" target="_blank">drescherjm@gmail.com</a>><br>
>> <mailto:<a href="mailto:drescherjm@gmail.com" target="_blank">drescherjm@gmail.com</a><br>
<mailto:<a href="mailto:drescherjm@gmail.com" target="_blank">drescherjm@gmail.com</a>>> <mailto:<a href="mailto:drescherjm@gmail.com" target="_blank">drescherjm@gmail.com</a><br>
<mailto:<a href="mailto:drescherjm@gmail.com" target="_blank">drescherjm@gmail.com</a>><br>
>> <mailto:<a href="mailto:drescherjm@gmail.com" target="_blank">drescherjm@gmail.com</a><br>
<mailto:<a href="mailto:drescherjm@gmail.com" target="_blank">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<br>
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>
>> <<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>><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>
>> _______________________________________________<br>
>> Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
<<a href="http://www.kitware.com" target="_blank">http://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<br>
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>><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>
>><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>
> 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>
><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>
<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>
</div></div></blockquote>
<br>
</blockquote></div><br><br clear="all"><br>-- <br>Ítalo Mendes<br><a href="http://ogre.irados.org">ogre.irados.org</a><br>