Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  630 1689 (N40E45 L=123)   594     91        2 ( 587/1656- 664/1732) -12.1 E53-E39 N38-N43
      2 1243 1603 (N42W10 L=178)   928    232        3 (1160/1551-1334/1682) -21.8 W02-W17 N38-N47
      3  418 1614 (N32E62 L=106)   576     75        0 ( 399/1575- 442/1647) -24.4 E67-E55 N29-N34
      4  636 1452 (N24E32 L=135)   616    195        4 ( 590/1378- 694/1531)  73.3 E33-E29 N18-N29
      5 1479 1219 (N18W28 L=196)   262    131        0 (1455/1173-1514/1280) -56.3 W26-W30 N15-N21
      6  271  872 (S16E54 L=114)  3817    532        7 ( 180/ 666- 377/1046)  38.5 E78-E43 S31-S04
      7 1669  676 (S14W50 L=218)   810    203        1 (1636/ 578-1691/ 768) -76.0 W47-W54 S20-S08
      8 1201  752 (S14W14 L=182)   249     69        1 (1169/ 738-1238/ 765)   0.0 W11-W16 S14-S12
      9 1633  510 (S25W54 L=222)  1015    120        2 (1583/ 483-1689/ 544)  24.0 W49-W59 S27-S23
     10 1243  440 (S34W24 L=191)   587    110        1 (1196/ 421-1297/ 460)  23.2 W20-W27 S35-S32
  

    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