Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  690 1683 (N37E29 L= 82)  3890    708        9 ( 489/1593- 913/1779) -17.2 E53-E10 N30-N45
      2  725 1421 (N19E21 L= 90)   237    117        0 ( 689/1371- 759/1459) -50.2 E23-E18 N15-N21
      3  579 1343 (N13E30 L= 81)   371     86        0 ( 563/1302- 596/1389) -80.5 E31-E28 N10-N16
      4 1330 1286 (N13W18 L=129)   505    173        1 (1239/1232-1402/1346) -28.6 W11-W22 N09-N16
      5 1863 1042 (N02W60 L=171)   467     76        1 (1847/ 999-1883/1076)  53.1 W58-W62 N00-N04
      6  943 1027 (S04E05 L=106)   256     67        0 ( 910/ 995- 967/1060)  63.4 E06-E03 S06-S02
      7  765  693 (S25E16 L= 96)  1005    274        0 ( 659/ 601- 904/ 796)  37.4 E23-E06 S31-S18
      8 1449  678 (S23W30 L=141)   250     50        0 (1421/ 663-1481/ 691)  14.0 W27-W32 S23-S21
      9 1090  550 (S33W07 L=118)   439     93        1 (1034/ 534-1141/ 569)   0.0 W03-W10 S34-S32
     10  704  321 (S53E27 L= 84)   295     43        0 ( 690/ 295- 718/ 349)  63.4 E28-E25 S54-S50
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y