Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1316 1317 (N13W25 L=340)   279    109        0 (1287/1277-1372/1347)   8.1 W22-W29 N11-N14
      2  181 1102 (N30E58 L=258)  2228    407        0 ( 113/ 987- 252/1209) -53.1 E64-E52 N21-N37
      3 1869 1029 (S19W57 L= 12)   733    138        1 (1849/ 968-1894/1088)  50.2 W54-W60 S23-S16
      4  833  963 (N07E12 L=303)   861    253        1 ( 717/ 922- 959/ 990) -23.2 E18-E04 N03-N10
      5  457  913 (N13E37 L=278)  1261    489        1 ( 269/ 839- 629/ 968)  -5.9 E53-E24 N11-N15
      6 1300  879 (S10W11 L=326)   306    105        1 (1246/ 850-1334/ 922)   8.1 W07-W14 S10-S07
      7 1507  836 (S18W23 L=338)   260     54        1 (1490/ 810-1523/ 862) -71.6 W21-W24 S19-S16
      8  193  783 (N11E65 L=250)   515     62        0 ( 175/ 756- 217/ 822) -76.0 E66-E62 N09-N13
      9 1164  766 (S12W01 L=316)   192     97        0 (1148/ 718-1182/ 821)  63.4 W00-W03 S15-S09
     10 1004  508 (S23E16 L=299)  1418    307        6 ( 973/ 375-1056/ 642)  47.0 E22-E08 S31-S16
     11  353  475 (S10E63 L=252)  5038    672       10 ( 285/ 368- 432/ 560)  49.8 E74-E55 S17-S04
     12  485  441 (S15E53 L=263)   473    120        2 ( 461/ 407- 511/ 464)  23.2 E56-E49 S16-S13
  

    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