plim.dream...@gmail.com wrote:
>> Yes, so what you need to do is to compute the ellipse
>> in a way that does not depend on the graphic device used and
>> then plot it. Can I just quickly ask what is your starting data
>> again (that is, how exactly you define your points in the first
>> place)?
>>
>> Paolo> B
>
>
> Hi Paolo and David,
> So, I went ahead an managed to create ellipses for the blobs!! I
> changed set_plot,'x' so now !D.X_size gives me 640 instead which the
> computer can crunch and produce nice figures on the screen. I would
> still like to convert the values for x,y, semimajor, semiminor, and
> center of the ellipse back to data coordinates such that I can place
> them on the postscript file again but convert_coord( ... /to_data)
> doesn't do the job.
>
> Some code to make this more interesting:
>
> contour_array = [s2c2(0,*),s2c2(1,*)] ; this is the x,y positions of
> the vertices of the polygon which defines the contour
May ask again what are your staring points, before you do the contour?
How are these points defind? I have this feeling that there may be a
much
simpler way...
Paolo
>
> blob_x = s2c2(0,*)
> blob_y = s2c2(1,*)
>
> set_plot,'x'
>
> polyfill,blob_x,blob_y,color=100,/data
>
> D = convert_coord(blob_x, blob_y, /to_device)
>
> result = polyfillv(d(0,*), d(1,*),640,640)
>
> Ellipse_out = fit_ellipse(result,semiaxes=semiaxes,center=center)
>
> plots,Ellipse_out,/device
>
> set_plot,'ps
> device,filename='Surface_Plots.ps'
> New_coords = convert_coord(Ellipse_out, /to_device)
> plots,s2c2(0,*),s2c2(1,*)
> plots,New_coords,/device
> device,/close
>
> Much thanks!
> Bruno
>
> p.s. I haven't banged my head on the IDLanROI more than being capable
> of making the object:
> blob_roi = OBJ_NEW('IDLanROI', blob_x,blob_y)
> blob_mask = blob_roi->ComputeMask(initialize=0)
> because I still don't see how to fit an ellipse and obtain its
> parameters using the object functions.
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
