| Re: is there a 3D equivalent of CONTOUR? [message #68301 is a reply to message #68240] |
Tue, 13 October 2009 03:02  |
dorthe
Messages: 12
Registered: April 2007
|
Junior Member |
|
|
On 8 Oct, 15:50, David Fanning <n...@dfanning.com> wrote:
> Dorthe Wildenschild writes:
>> excellent, - but I can't get it to work.. not very good at working
>> with objects I'm afraid - I do it once every 3 years or so :-)
>
>> Here's what I do:
>
>> I use LABEL_REGION to separate all the objects in my image and then I
>> pick one in a loop structure, representing "t" below
>
>> ISOSURFACE, label_image,
>> t, Outverts, Outconn
>> ; print, outverts
>> Roi=Obj_New('IDLanROI', Outverts, type=2)
>
>> Status = Roi-> ComputeGeometry(Area=area,Perimeter=perimeter,
>> Centroid=centroid)
>
>> t_perimeter(r)=perimeter
>> t_area(r-1)=area
>> t_x_center(r)=centroid[0]
>> t_y_center(r)=centroid[1]
>> blob_siz(r)=blob_size[1]
>
>> I get lots of values in Outverts, but area, perimeter etc. are all
>> zero?
>> Any ideas for the slightly incompetent...?
>
> Humm, well, it must have something to do with what you are
> passing to isosurface. I have no problems with this code:
>
> IDL> vol = loaddata(8)
> IDL> minmax, vol
> 0 232
> IDL> isosurface, vol, 150, verts, conn
> IDL> help, verts
> VERTS FLOAT = Array[3, 178729]
> IDL> roi = obj_new('idlanroi', verts)
> IDL> status = roi -> computegeometry(perimeter=p, area=a, centroid=c)
> IDL> print, p
> 361882.75
> IDL> print, a
> 5674.0855
> IDL> print, c
> -35.847264 107.88590 42.979445
>
>> BTW, nice to see you're still very active on this site.
>
> "Active" is probably too strong a verb. Most of what is
> discussed here these days is over my head. :-)
>
> Cheers,
>
> David
>
> --
> David Fanning, Ph.D.
> Fanning Software Consulting, Inc.
> Coyote's Guide to IDL Programming:http://www.dfanning.com/
> Sepore ma de ni thui. ("Perhaps thou speakest truth.")- Hide quoted text -
>
> - Show quoted text -
Ok, so I need to step back and make sure I'm passing correct data to
isosurface.
Or bail and use Avizo Fire instead - oops perhaps I lose my right to
post here if I do that :-)
Thanks for helping, I'll give it another try.
Dorthe
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
