Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  888 1804 (N43E28 L= 91)  1547    178        4 ( 821/1731- 944/1859) -29.1 E38-E20 N37-N47
      2 1349 1667 (N41W12 L=132)   289     40        0 (1328/1644-1377/1687)  51.3 W10-W14 N38-N43
      3 1186 1628 (N34W00 L=120)   362     75        1 (1164/1590-1214/1667) -59.0 E01-W02 N31-N37
      4  634 1628 (N26E38 L= 81)   680    187        0 ( 551/1566- 696/1666) -14.9 E45-E30 N22-N28
      5 1671 1590 (N44W44 L=164)   621     71        0 (1646/1555-1693/1626)  38.7 W39-W49 N39-N47
      6 1768 1323 (N26W48 L=168)  1771    365        2 (1658/1173-1852/1447) -37.4 W38-W57 N17-N32
      7  159 1212 (S04E68 L= 51)  2043    230        6 ( 106/1138- 207/1296)  36.0 E76-E63 S08-N00
      8  451 1120 (S08E38 L= 82)   551    175        0 ( 430/1041- 483/1199) -60.9 E39-E34 S12-S03
      9  168 1016 (S16E64 L= 55)  5932    462       14 ( 104/ 833- 233/1165) -90.0 E75-E56 S27-S07
     10  274  873 (S24E52 L= 67)  2527    507        6 ( 245/ 619- 301/1086)  80.2 E58-E49 S40-S11
     11 1720  937 (N01W47 L=166)   515    161        1 (1658/ 908-1780/ 970)  16.7 W41-W51 N00-N03
     12 1713  828 (S05W49 L=168)   680    200        0 (1630/ 767-1771/ 859)  -4.1 W41-W55 S07-S03
     13 1662  619 (S18W52 L=171)   705    104        1 (1619/ 576-1706/ 660)  80.5 W49-W54 S20-S14
     14 1670  533 (S23W58 L=177)   781    107        1 (1623/ 510-1715/ 565)  18.4 W52-W63 S23-S20
     15 1067  481 (S39W14 L=133)   265     59        0 (1051/ 448-1083/ 510) -53.1 W12-W15 S40-S36
     16 1521  422 (S33W52 L=171)   772    138        3 (1452/ 395-1576/ 447)  24.0 W46-W57 S35-S30
     17  743  422 (S49E11 L=108)   334     41        1 ( 714/ 400- 768/ 446)  21.8 E13-E08 S51-S47
  

    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