;SE METTRE EN IDL ;IDL ;LANCER IDL SESSION ;@IDLSESSION ;CATALOGUE DE DARKS CD,'DATAVS5:[LASCO.IMAGES.MV_NRL.NRLOCT]' DARK_CATA ;TRAITEMENT IMAGES ACQUISITIONS EN COMPATIBLE IDL ;CAS DES LAMPES A PRO_CAL,'C2CSTI300.IMG',DARK,HDARK,STI300,HSTI300 PRO_CAL,'C2CSTI301.IMG',DARK,HDARK,STI301,HSTI301 CD,'DATAVS5:[LASCO.IMAGES.MV_NRL.NRLNOV]' DARK_CATA PRO_CAL,'C2CSTI400.IMG',DARK,HDARK,STI400,HSTI400 PRO_CAL,'C2CSTI401.IMG',DARK,HDARK,STI401,HSTI401 CD,'DATAVS5:[LASCO.IMAGES.MV_NRL.NRLOCT]' WSET,0 TVSCL,BYTSCL(STI300,MIN=MIN(STI300),MAX=MAX(STI300)) WSET,1 TVSCL,BYTSCL(STI400,MIN=MIN(STI400),MAX=MAX(STI400)) ;VISUALISER LE DECALAGE: DIF300=STI400-STI300 ---> LISERET EN HAUTEUR #1.5 PX TVSCL,BYTSCL(DIF300,MIN=MIN(DIF300),MAX=MAX(DIF300)) ;DETERMINER LE DECALAGE : ;PAR EXAMEN DE LA POSITION DU BORD SUPERIEUR ET DU BORD GAUCHE ;DE L'OMBRE DES 2 IMAGES AVEC: PROFILES,STI400 ET PROFILES,STI300 (DECROCHEMENT A MI HAUTEUR) DECALER L'IMAGE STI400 PAR STIR400=ROT_INT(STI400,0.,1.,255.5,254.0) ;RECONTROLER LA DIFFERENCE DIFR300=STIR400-STI300 ;VISUALISER: TVSCL,BYTSCL(DIFR300,MIN=MIN(DIFR300),MAX=MAX(DIFR300)) ;FAIRE UN STAT_CURSOR ET AJUSTER LES CUTS TVSCL,BYTSCL(DIFR300,MIN=-100.,MAX=300.) :REITERER SI NECESSAIRE ;VERIFIER PAR RDPIX LA DISPARITION DU LISERET. OPERER DE MEME POUR LAMPES B WSET,0 TVSCL,BYTSCL(STI301,MIN=MIN(STI301),MAX=MAX(STI301)) WSET,1 TVSCL,BYTSCL(STI401,MIN=MIN(STI401),MAX=MAX(STI401)) DIF301=STI401-STI301 TVSCL,BYTSCL(DIF301,MIN=MIN(DIF301),MAX=MAX(DIF301)) ;REITERER STIR401=ROT_INT(STI401,0.,1.,255.5,254.0) DIFR301=STIR401-STI301 TVSCL,BYTSCL(DIFR301,MIN=MIN(DIFR301),MAX=MAX(DIFR301)) STAT_CURSOR,DIFR301 RDPIX,DIFR301 ;LE DECALAGE EST: STIR401=ROT_INT(STI401,0.,1.,255.22,254.15) TVSCL,BYTSCL(DIFR301,MIN=-100.,MAX=300.) LOADCT,4 ADJCT STAT_IMA,DIFR301,0,0,511,511 ; ;DECALAGE TROUVE : EN X 255.5-255.22=0.28 ; EN Y 255.5-254-15=1.45 ;; CONTROLE PAR LA METHODE DES CENTRES ;UTILISER SIMCIRCLE,NOM DE L'IMAGE POUR TROUVER LE CENTRE D'UN CERCLE ;OCCULTEUR EXTERNE AVANT ET APRES LE TEST ; TRACE UN CERCLE A MI HAUTEUR ;EXEMPLE : ;WSET,0 ;SIMCIRCLE,STI300 X_CENTRE = 201.077 Y_CENTRE = 326.861 RAYON= 172.336 ;WSET,1 ;SIMCIRCLE,STI400 ;X_CENTRE = 200.711 Y_CENTRE = 325.226 RAYON= 172.624 PRINT,201.077-200.711 0.365997 PRINT , PRINT,326.861-325.226 1.63498 DECALAGE TROUVE : EN X=0.36 EN Y=1.63 *********************** DEPLACEMENT DES MARQUES *********************** WSET,0 STAT_CURSOR,STI300,64,64 For area [ 431, 415: 494, 478] WSET,1 STAT_CURSOR,STI400,64,64 For area [ 430, 416: 493, 479] IL YA # -1 PX EN X ET +1 PX EN Y FENETRE GLOBALE 430:494 EN X ET 415:479 EN Y 1-CROSSCORR,STI300,STI400,430,494,415,479,IMA_CROS,RESMAX >>> begin study of crosscorrelation's function <<< >>> research of the off-centering (precision : 1 pixel) <<< deltax : 0.500000 deltay : 0.500000 >>> research of the off-centering (precision : fractionnal pixel) <<< >>> start of surface fitting <<< >>> procedure of cubic interpolation <<< deltax : 0.500000 deltay : 0.500000 2-CROSSCORRBIS,STI300,STI400,430,494,415,479,IMA_CROS2,RESMAX2 >>> begin study of crosscorrelation's function <<< >>> research of the off-centering (precision : 1 pixel) <<< deltax : 0.500000 deltay : 0.500000 >>> research of the off-centering (precision : fractionnal pixel) <<< >>> start of surface fitting <<< >>> procedure of sinc interpolation <<< ... FFT given (source) image ... ... insert source into destination ... ... inverse fft ... deltax : 0.400000 deltay : 0.600000 IL EST IMPOSSIBLE DE CONCLURE :LIMITE DE LA PRECISION PB:UN ESSAI FAIT AVEC LA MEME IMAGE ---> 0.5 ET 0.5 AVEC CROSSCORR ET CROSSCORRBIS !!!:oui car il faut travailler avec une fenetre en x et y pairs MARQUES RELATIVES A STI301 STI401: EN HAUT A GAUCHE STAT_CURSOR,STI301,64,64 For area [ 46, 417: 109, 480] MEME ZONE. CROSSCORR,STI301,STI401,46,109,417,480,IMA_CROS1,RESMAX1 deltax : 0.000000 deltay : 0.000000 procedure of cubic interpolation <<< deltax : 0.000000 deltay : 0.000000 EN BAS A DROITE STAT_CURSOR,STI401,64,64 For area [ 427, 26: 490, 89] MEME ZONE CROSSCORR,STI301,STI401,427,490,26,89,IMA_CROS2,RESMAX2 deltax : 0.000000 deltay : 0.000000 procedure of cubic interpolation <<< deltax : 0.000000 deltay : 0.000000