| Re: is there a 3D equivalent of CONTOUR? [message #68240] |
Thu, 08 October 2009 06:50  |
David Fanning
Messages: 11724
Registered: August 2001
|
Senior Member |
|
|
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.")
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
