Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  624 1567 (N35E36 L=133)   461    111        1 ( 578/1519- 655/1598) -26.6 E41-E31 N31-N37
      2 1611 1409 (N30W43 L=212)  3168    650        6 (1392/1240-1836/1550) -27.2 W25-W64 N18-N39
      3  295 1470 (N25E65 L=104)   495     64        1 ( 271/1446- 326/1492)  45.0 E68-E62 N23-N27
      4  978 1301 (N19E04 L=165)   212     96        0 ( 945/1269-1017/1323) -31.0 E06-E01 N17-N20
      5  225  967 (S06E59 L=110)   507    123        0 ( 196/ 912- 244/1009)  50.2 E63-E56 S09-S03
      6 1753  858 (S05W53 L=222)   515     60        0 (1723/ 836-1778/ 878)  26.6 W49-W55 S06-S03
      7  218  749 (S20E64 L=105)   853    148        0 ( 163/ 727- 285/ 777)  -3.2 E73-E55 S21-S19
      8  247  651 (S27E64 L=105)  1136    149        2 ( 191/ 596- 283/ 706) -37.9 E71-E60 S30-S23
      9 1677  638 (S19W51 L=220)   465     95        0 (1651/ 603-1702/ 694)  80.5 W48-W52 S21-S15
     10 1480  646 (S19W33 L=203)   281    115        1 (1436/ 603-1526/ 680)  51.3 W31-W36 S22-S16
     11 1679  555 (S25W55 L=224)   385     53        1 (1652/ 541-1706/ 576)  11.3 W52-W58 S25-S23
     12  605  523 (S33E29 L=140)   455    100        0 ( 556/ 503- 651/ 545)  -8.1 E33-E25 S34-S31
  

    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