Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1528 1758 (N43W60 L=123)   966    149        0 (1498/1719-1559/1795)  71.6 W58-W61 N39-N45
      2 1855 1364 (N11W72 L=136)  1127    170        1 (1817/1310-1892/1407) -45.0 W67-W78 N07-N14
      3  539 1355 (N32E30 L= 33)   345    112        0 ( 491/1339- 586/1375) -20.6 E34-E26 N31-N34
      4  261 1278 (N29E55 L=  8)  2138    319        4 ( 172/1220- 346/1344)  20.6 E64-E46 N25-N33
      5 1584 1294 (N13W41 L=105)   383    166        1 (1513/1263-1624/1335)  11.3 W35-W45 N12-N15
      6  666 1109 (N15E22 L= 42)   704    241        2 ( 569/1038- 764/1192) -42.5 E27-E15 N10-N21
      7  214  707 (S06E67 L=356)  2302    343        8 ( 176/ 592- 252/ 824)  75.1 E73-E63 S14-N01
      8  352  709 (S06E51 L= 12)   775    167        2 ( 306/ 676- 414/ 744) -18.4 E56-E46 S08-S04
      9 1773  633 (S32W52 L=115)  1219    277        2 (1691/ 588-1852/ 706)  10.0 W43-W60 S34-S28
     10  326  619 (S12E57 L=  7)   604    110        1 ( 309/ 571- 343/ 672)  66.8 E58-E54 S15-S08
     11 1220  606 (S23W07 L= 70)   452    110        1 (1162/ 573-1279/ 637)   7.1 W02-W10 S24-S20
  

    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