Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  886 1860 (N42E45 L=149)  1061    130        2 ( 830/1822- 948/1893)   0.0 E52-E36 N39-N44
      2  544 1680 (N19E54 L=140)  1439    145        4 ( 481/1653- 610/1713)   8.7 E61-E47 N16-N22
      3 1150 1624 (N32E11 L=183)  1162    171        1 (1096/1533-1194/1701) -33.7 E17-E05 N26-N36
      4  801 1557 (N17E29 L=165)   701    182        1 ( 756/1507- 865/1613) -21.8 E33-E23 N15-N21
      5  729 1496 (N12E31 L=163)   484    162        0 ( 682/1445- 785/1562) -20.0 E36-E25 N10-N14
      6  430 1391 (S01E47 L=147)   284     38        0 ( 412/1361- 444/1416) -33.7 E49-E45 S02-N00
      7 1570 1289 (N24W25 L=220)   879    242        0 (1526/1200-1595/1393) -76.0 W20-W27 N17-N29
      8 1736 1287 (N31W40 L=234)   758    159        0 (1685/1213-1761/1335) -32.0 W34-W42 N26-N32
      9  408 1277 (S08E45 L=149)   566    114        0 ( 382/1218- 445/1324)  81.9 E46-E43 S11-S04
     10 1018 1272 (N07E07 L=187)   606    169        0 ( 969/1217-1097/1310)  54.5 E08-E03 N02-N09
     11  505 1219 (S09E36 L=158)   552    165        0 ( 438/1165- 542/1282)  63.4 E39-E34 S13-S04
     12 1593 1187 (N19W29 L=224)   243     22        0 (1577/1164-1610/1206)  71.6 W27-W30 N17-N20
     13 1053  913 (S12W05 L=199)   194     32        0 (1037/ 890-1070/ 939) -33.7 W03-W06 S13-S10
     14  189  864 (S36E63 L=132)   863    152        0 ( 140/ 828- 248/ 904)  16.7 E74-E53 S39-S33
     15 1245  818 (S12W18 L=213)   367    129        3 (1198/ 783-1288/ 848)  38.7 W15-W20 S13-S09
     16 1003  668 (S27W10 L=204)   376    100        2 ( 959/ 644-1055/ 690)  18.4 W06-W12 S28-S25
     17 1452  347 (S30W57 L=251)  1113    183        0 (1355/ 328-1529/ 374)  17.4 W48-W65 S33-S27
  

    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