Il giorno lunedì 16 ottobre 2017 11:20:59 UTC+2, Markus Schmassmann ha scritto:
> On 10/15/2017 06:08 PM, termybel@gmail.com wrote:
>> Il giorno giovedì 12 ottobre 2017 16:27:18 UTC+2, Markus Schmassmann ha scritto:
>>> On 10/12/2017 02:10 PM, termybel@gmail.com wrote:
>>>> Thanks to answer but I have a probelm.
>>>>
>>>> You puth this: if n_elements(li) ne 1 then message, 'not exactly 1 contour'
>>>>
>>>> and my program say: not exactly 1 contour
>>>>
>>>> How I can solve it?
>>>>
>>>> If I ask in idl " help,/str,li" and the answer is:
>>>>
>>>> Structure CONTOUR_DBL_PATH_STRUCTURE, 6 tags, length=24, data length=20:
>>>> TYPE BYTE 1
>>>> HIGH_LOW BYTE 1
>>>> LEVEL INT 0
>>>> N LONG 7
>>>> OFFSET LONG 0
>>>> VALUE DOUBLE 9.0000000
>>>>
>>>> and for " help,/str,lines"
>>>> LINES DOUBLE = Array[2, 5447]
>>>>
>>>> My contour isn't a perfect circle. Maybe is this the problem?
>>>>
>>>> This is my program where the image is a flat field panel minus dark:
>>>>
>>>> file_ff1='ff_100s_3.fits' ;immagine
>>>> immagine_ff1=readfits (file_ff1, header1) ; leggo l'immagine del flat field panel
>>>> file_dark='ff_100s_dark_3.fits' ;immagine
>>>> immagine_dark=readfits (file_dark, header1) ; leggo l'immagine della dark
>>>>
>>>> n=100. ; secondi di esposizione dell'immagine
>>>> ffp=(immagine_ff1-immagine_dark)/n
>>>>
>>>> ; contour per selezionare soglie di equivalore
>>>>
>>>> speriamo=fltarr(2048,2048) ; creo una matrice 2048x2048
>>>> dimensioni=size(speriamo,/dimensions)
>>>> cubo=where(ffp lt 9.5 or ffp gt 10)
>>>>
>>>>
>>>>
>>>> wrong_matrix=array_indices(dimensioni, cubo,/dimensions)
>>>> indxw=reform(wrong_matrix(0,*))
>>>> indyw=reform(wrong_matrix(1,*))
>>>> ;ffp(indxw,indyw)=0.
>>>> ;c = CONTOUR(ffp, dimensions=[512,512], Title='prova cubo')
>>>>
>>>> ; DATA being the data to be contoured
>>>> level=9
>>>> contour, smooth(ffp(250:1700, 250:1700),3), path_info=li,closed=1 , path_xy=lines, /path_data_coord, levels=[level], /path_double
>>>> ; lix=lines(0,*)
>>>> ; liy=lines(1,*)
>>>> ; liyd=deriv(lix,liy)
>>>> ;ind=where(abs(liyd) le 0.0001)
>>>>
>>>> if n_elements(li) ne 1 then message, 'not exactly 1 contour'
>>>>
>>>>
>>>> cont_obj =obj_new('IDLanROI',lines)
>>>> void= cont_obj.ComputeGeometry(centroid=center)
>>>> ;fit_ellipse(
>>>> print, center[0:1]
>>>>
>>>> end
>>>
>>> contour, smooth(ffp(250:1700, 250:1700),3), path_info=li,closed=1 , $
>>> path_xy=lines, /path_data_coord, levels=[level], /path_double
>>> contour, smooth(ffp(250:1700, 250:1700),3), levels=[level]
>>> ; shows you there are more than 1 contour,
>>> ; you need to identify the correct one
>>> ; often the best is the longest
>>> void=max(li.n,j)
>>> line=[*,li[j].offset+lindgen(li[j].n)]
>>> plot, line[0,*],line[1,*]
>>> cont_obj =obj_new('IDLanROI',line)
>>
>> I don't know why but IDL says this message (it doesn't read the *):
>>
>> line=[*,li[j].offset+lindgen(li[j].n)]
>> ^
>> % Syntax error.
>>
> ; sorry, it should be
> line=lines[*,li[j].offset+lindgen(li[j].n)]
Thank you so much! :D :D :D
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
