| Re: IsoSurface [message #39798 is a reply to message #39725] |
Thu, 10 June 2004 10:45   |
siliconcube
Messages: 11
Registered: June 2004
|
Junior Member |
|
|
David Fanning <davidf@dfanning.com> wrote in message news:<MPG.1b3136406ab74403989792@news.frii.com>...
> Aleks writes:
>
>> ok here is the problem I have: I'm trying to run a test on different
>> images to figure out how things work. Previously I was helped to find
>> the size of my tiff image. I modified the code accordingly.
>> Unfortunatelly, when I try to run this code on my images the first
>> Window function shows me all my images but when the execution gets to
>> the second Window Function I lose IDL ie it crashes. I have 81 nearly
>> identical image. (if you would like to see the image you can check it
>> out here http://www.geocities.com/siliconcube/tree_01.tif
>> images were modified with MatLab).
>
> Oh, oh. There's your problem. IDL is not going to
> work with anything created in MatLab!!
>
> (No, I'm kidding. It's just a joke.)
>
>> Can anyone help me identify the
>> problem, is my image not suitable for Iso Surface or am I doing
>> somethign wrong =/.
>
> What makes you think you are doing something wrong?
> Did you get an error when you ran this program? Or
> are you just not seeing anything on your display?
>
> Do you have any idea of the range of values in your
> data? (You could print the min and max of your volume.)
> Is 81 a smart choice for the threshold? Where did 81 come
> from? It seems an odd choice. Maybe you could plot a
> histogram of your volume data and find a value that
> makes more sense for the specific data your have.
>
> Cheers,
>
> David
This is the compilation log that I have when I run the program:
IDL Version 6.0, Microsoft Windows (Win32 x86 m32). (c) 2003, Research Systems, Inc.
IDL> .compile test
% Compiled module: $MAIN$.
IDL> .go
% Compiled module: BSORT.
% Loaded DLM: TIFF.
% Compiled module: LOADCT.
% Compiled module: FILEPATH.
% Compiled module: PATH_SEP.
% LOADCT: Loading table B-W LINEAR
% Compiled module: SCALE3.
% Compiled module: T3D.
% POLYSHADE: Expression must be an array in this context: POLYGONS.
% Execution halted at: $MAIN$ 19 C:\RSI\IDL60\test.pro
And this is the code
files=findfile('*.tif')
index=bsort(files, sortedfiles)
volume=bytarr(761,671,81)
for j=0,80 do begin
image=read_tiff(sortedfiles[j])
volume[0,0,j] = image
endfor
;Window, Title='Original Image Slices', XSize=800, YSize=500, 1, XPos=0, YPos=0
LoadCT, 0
Device, Decomposed=0
FOR j=0,80 DO TV, volume[*,*,j], j
Window, 6, XSize=800, YSize=800, Title='Test IsoSurface'
Scale3, XRange=[0,760], YRange=[0,670], ZRange=[0,80], AZ=-150
;Shade_Volume, volume, 81, vertices, polygons, /Low
theHead = PolyShade(vertices, polygons, /T3D)
TV, theHead
end
Thank you
Aleks
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
