Changeset 467 for trunk

Timestamp:

12/15/11 15:07:54 (12 years ago)

Author:

lelod

Message:

validation faite

Location:

trunk/src

Files:

• extract_amsua.pro (modified) (7 diffs)
• interpolswath.pro (modified) (4 diffs)

trunk/src/extract_amsua.pro

@@ -93 +93 @@
 ; ==========
 ;
+; lelod 20111215
+; tests de fonctionnement effectues: la correction nadir est
+;correctement implementee
+; l'interpolation de la fauchee fonctionne
+; mais attention test sur la longitude code en dur:
+; pbs avec les fauchees qui traversent le meridien 180deg!!!!
+; l'interpolation ajoute des points parasites dans la zone...
 ; $Id$
 ;
@@ -192 +199 @@
 ; TODO
 ; ====
+;
+; traiter les longitudes autour de 180 degres (passage de 180 a -180)
+; interpolswath ne gere pas cette situation
 ;
 ; comprendre les points bizarres autour de -40,30 visibles avec l'exemple
@@ -265 +275 @@
 openr, lun_a, list_filea, /get_lun
 filea=''
-nlist=1
+nlist=0
 ilist=0
-   ; while (not eof(lun_a)) do begin
 loadct,39
 numwindow=1
 ;while ilist le nlist do begin
- while (not eof(lun_a)) do begin
+while (not eof(lun_a)) do begin
    ilist=ilist+1
    onefile = ''
@@ -370 +379 @@
 ;   colscale[j]=(tb[j]-min(tb))*255/(max(tb)-min(tb))
 ;endfor
-;plots,lolo,lala,psym=5,symsize=2,color=colscale
+;plots,lolo,lala,psym=5,symsize=1,color=colscale
 ;numwindow=numwindow+1
 
@@ -376 +385 @@
 ; reguliere dans la fauchee
 ; et selection de la zone conservee
-;tek_color
-;resol=1
+; attention test sur la longitude code en dur:
+; pbs avec les fauchees qui traversent le meridien 180deg!!!!
+; l'interpolation ajoute des points parasites dans la zone...
    nbgrid=0
+cont=0L
    for i=0L,n_scan-1L do begin
       tb=ch_nadir[*,i]
@@ -384 +395 @@
       lat=amalat[*,i]
       mask=landseamask[*,i]
-      nbgrid1=nbgrid
-      interpolswath,tb,lat,lon,mask,resol,nbgrid,tbgrid,latgrid,longrid,maskgrid
-      if i eq 0L then begin
-         fovgrid=indgen(nbgrid)+1
-      endif
-      ;ind=where(tbgrid gt tbmin and tbgrid lt tbmax,nbon)
-     ; if (nbon ne 0) then begin
-         longit=longrid
-         latit=latgrid
-      ;   oplot,longit,latit,psym=3,color=0
-     ; endif else begin
-     ;    print,'pas de point valable dans la fauchee'
-     ; endelse
+      ind=where(tb gt 0 and lon gt -60 and lon lt 100,nbon)
+      if nbon eq nbpix then begin
+         cont=cont+1L
+         nbgrid1=nbgrid
+         interpolswath,tb,lat,lon,mask,resol,nbgrid,tbgrid,latgrid,longrid,maskgrid
+         if cont eq 1L then begin
+            fovgrid=indgen(nbgrid)+1
+         endif
+         ;longit=longrid
+         ;latit=latgrid
+         ;oplot,longit,latit,psym=3,color=0
 
 ; selection des taches au sol situees dans la zone d'interet
-      zone=where((longrid ge lon_min) and (longrid le lon_max) $
-                 and (latgrid ge lat_min) and (latgrid le lat_max) and (tbgrid gt tbmin), npt)
-      if npt ne 0 then begin
-         if (n_elements(chint) EQ 0) then begin
-            timeint=replicate(tt[i],npt)
-            chint=tbgrid[zone]
-            latint=latgrid[zone]
-            lonint=longrid[zone]
-            fovint=fovgrid[zone]
-            maskint=maskgrid[zone]
+         zone=where((longrid ge lon_min) and (longrid le lon_max) $
+                    and (latgrid ge lat_min) and (latgrid le lat_max) and (tbgrid gt tbmin), npt)
+         if npt ne 0 then begin
+            if (n_elements(chint) EQ 0) then begin
+               timeint=replicate(tt[i],npt)
+               chint=tbgrid[zone]
+               latint=latgrid[zone]
+               lonint=longrid[zone]
+               fovint=fovgrid[zone]
+               maskint=maskgrid[zone]
+            endif else begin
+               chint=[chint,tbgrid[zone]]
+               latint=[latint,latgrid[zone]]
+               lonint=[lonint,longrid[zone]]
+               fovint=[fovint,fovgrid[zone]]
+               maskint=[maskint,maskgrid[zone]]
+               timeint=[timeint,replicate(tt[i],npt)]
+            endelse
          endif else begin
-            chint=[chint,tbgrid[zone]]
-            latint=[latint,latgrid[zone]]
-            lonint=[lonint,longrid[zone]]
-            fovint=[fovint,fovgrid[zone]]
-            maskint=[maskint,maskgrid[zone]]
-            timeint=[timeint,replicate(tt[i],npt)]
+            print, 'www : aucun point dans la zone'
          endelse
       endif else begin
-         print, 'www : aucun point dans la zone'
+         print,'fauchee incomplete',nbon
       endelse
 
@@ -485 +497 @@
 plot,lonint,latint,xrange=[lon_min,lon_max],yrange=[lat_min,lat_max],xstyle=1,ystyle=1,/nodata
 ind=where(maskint eq 1)
-oplot,lonint[ind],latint[ind],psym=1
+;oplot,lonint[ind],latint[ind],psym=1
 jnd=where(chint gt tbmin and chint lt tbmax)
 tb=chint[jnd]

trunk/src/interpolswath.pro

@@ -40 +40 @@
 ;   * minimal rest header
 ;
+; TO DO
+; =====
+; regler le pb des longitudes autour de 180 degres
+; l'interpolation ne les gere pas et cree des points dans la zone AMMA...
 ;-
 pro interpolswath, tb, latlu, lonlu, masklu, resol, nbgrid, tbint, latgrid, longrid, mask
@@ -70 +74 @@
 ;print,pixatot
 ;window,1
-;ploterr,ifov,fova,pixatot
+ ;ploterr,ifov,fova,pixatot
 
 distmax=-min(fova)+max(fova)+max(pixa)
@@ -79 +83 @@
 if (dmax/delta-fix(dmax/delta) le 0.5) then ndemigrid=ndemigrid-1
 nbgrid=ndemigrid*2
-print,'nb pts',nbgrid, ' dmax',dmax
+;print,'nb pts',nbgrid, ' dmax',dmax
 
 grid2=indgen(ndemigrid)*delta+fova[nnadir]
@@ -116 +120 @@
 longrid=interpol(lonlu,fova,grid)
 
+ind=where(longrid gt 35 and longrid lt 45 and latgrid gt 10 and latgrid lt 25,nbind)
+if nbind ne 0 then begin
+  ; print,'longitude intiniale',lonlu
+  ; print,'long. interpolee',longrid
+endif
 ;print,tbint
 end