Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  519 1742 (N23E63 L= 67)   757     67        2 ( 494/1710- 549/1765) -16.7 E68-E58 N20-N24
      2  839 1609 (N22E29 L=101)  1157    252        3 ( 806/1497- 884/1706) -55.3 E33-E23 N15-N28
      3 1457 1622 (N44W10 L=140)   886    140        1 (1393/1588-1524/1645)  14.0 W04-W16 N41-N46
      4 1290 1529 (N30W01 L=130)   469     74        2 (1258/1502-1333/1555)   0.0 E01-W03 N28-N32
      5 1301 1319 (N18W07 L=137)   596    258        2 (1239/1265-1404/1396)   0.0 W02-W14 N15-N22
      6 1822 1313 (N37W50 L=180)   894    102        2 (1784/1269-1849/1360) -29.7 W45-W54 N34-N39
      7 1455 1302 (N22W17 L=147)   290     97        0 (1420/1276-1481/1344)  78.7 W15-W18 N19-N24
      8 1778 1160 (N25W45 L=175)   524     91        0 (1739/1130-1815/1183)  40.6 W41-W49 N21-N27
      9 1698 1029 (N14W41 L=171)  1407    327        1 (1578/ 986-1808/1105)  10.8 W29-W50 N11-N19
     10  214 1021 (S27E59 L= 71)   684    155        0 ( 186/ 953- 237/1086) -58.0 E63-E55 S30-S22
     11  801  797 (S25E07 L=123)  1403    219        7 ( 710/ 735- 894/ 858)  50.7 E11-E01 S31-S19
     12 1702  632 (S06W56 L=186)   492    102        0 (1665/ 604-1737/ 648)  33.7 W53-W59 S08-S04
     13  763  612 (S37E04 L=126)   642    185        1 ( 701/ 584- 848/ 641)  24.0 E08-W01 S39-S33
     14  583  555 (S47E16 L=113)  1412    208        6 ( 536/ 502- 652/ 632) -26.6 E20-E09 S51-S41
     15 1487  527 (S19W45 L=175)   644    112        1 (1450/ 490-1528/ 579) -45.0 W41-W49 S21-S16
     16 1676  483 (S14W65 L=194)   953    140        4 (1632/ 445-1735/ 518)  49.4 W59-W70 S17-S10
     17 1715  443 (S15W74 L=204)  2834    179        5 (1660/ 384-1779/ 504)  73.3 W68-W81 S19-S09
  

    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