Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  380 1542 (N19E54 L=343)  1171    140        4 ( 332/1485- 430/1590)  66.8 E59-E50 N15-N22
      2  459 1476 (N15E44 L=353)  3555    412       12 ( 302/1386- 612/1548)   9.1 E57-E32 N09-N20
      3 1796 1465 (N35W56 L= 93)   538     53        2 (1767/1443-1825/1491)  11.3 W51-W61 N32-N37
      4  224 1363 (N07E63 L=334)  3394    263        9 ( 169/1287- 291/1440) -39.3 E72-E55 N02-N11
      5 1078 1072 (S03W02 L= 40)   371    139        1 (1059/1004-1094/1144) -81.9 W01-W03 S06-N01
      6  373 1051 (S11E42 L=355)   749    200        1 ( 310/ 974- 421/1136)  59.0 E47-E38 S15-S05
      7  258 1066 (S11E53 L=344)   619     77        1 ( 232/1025- 277/1110)  59.0 E55-E51 S13-S08
      8 1940 1039 (N10W70 L=107)  1570    146        1 (1911/ 962-1963/1104)  66.0 W65-W74 N05-N14
      9  988  905 (S13E01 L= 36)   257    107        1 ( 941/ 884-1043/ 927)  38.7 E03-W02 S15-S11
     10  292  826 (S25E50 L=347)  2520    391        0 ( 135/ 763- 466/ 872)   4.6 E71-E34 S29-S21
     11  678  850 (S20E19 L= 18)   299    116        0 ( 615/ 835- 727/ 866)   8.1 E23-E16 S21-S18
     12 1172  740 (S21W13 L= 50)   598    166        0 (1154/ 656-1189/ 817) -66.0 W10-W14 S25-S16
     13  721  671 (S31E15 L= 22)  1959    463        2 ( 508/ 624- 884/ 783)  29.1 E32-E05 S36-S21
     14 1607  628 (S20W46 L= 83)   594    159        2 (1585/ 572-1660/ 700) -45.0 W41-W51 S23-S16
     15 1398  648 (S23W30 L= 67)   379     90        0 (1351/ 621-1433/ 683)  45.0 W27-W31 S24-S19
     16 1213  377 (S43W27 L= 64)  2487    407        6 (1034/ 330-1401/ 433)   7.1 W11-W43 S45-S39
  

    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