Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1106 1517 (N28W11 L= 72)   381    171        0 (1017/1488-1182/1535) -21.8 W05-W16 N25-N30
      2  801 1419 (N25E10 L= 50)   778    184        2 ( 696/1379- 872/1445) -22.6 E17-E05 N22-N28
      3  237 1256 (N22E56 L=  5)  1826    288        2 ( 162/1160- 347/1342) -26.6 E69-E43 N14-N28
      4  642 1046 (N04E23 L= 37)   425    148        0 ( 577/1010- 700/1100) -40.6 E26-E19 N02-N08
      5 1704  822 (S20W44 L=105)  1661    426        3 (1647/ 639-1782/1017)  73.1 W39-W51 S31-S08
      6  388  976 (N03E41 L= 19)   312     53        0 ( 374/ 936- 405/1003)  63.4 E42-E39 N00-N04
      7 1208  904 (S10W10 L= 70)   251     97        1 (1171/ 872-1255/ 929)  31.0 W07-W12 S11-S08
      8 1766  671 (S30W54 L=114)   956    119        1 (1716/ 637-1806/ 719)  18.4 W48-W57 S32-S27
      9  371  489 (S26E58 L=  2)  1855    214        3 ( 327/ 374- 410/ 585)  31.3 E75-E50 S33-S19
  

    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