Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  499 1658 (N21E52 L= 25)   868     94        3 ( 453/1628- 535/1695)  59.0 E55-E48 N18-N23
      2  391 1633 (N18E61 L= 16)   676     93        1 ( 362/1607- 429/1672) -11.3 E65-E55 N16-N20
      3  651 1584 (N18E37 L= 41)   263     72        0 ( 625/1550- 683/1615)  78.7 E37-E35 N15-N20
      4 1640 1382 (N29W32 L=110)   283     64        0 (1613/1358-1672/1406)  51.3 W30-W35 N26-N31
      5  270 1387 (N01E60 L= 18)   839    139        0 ( 202/1366- 341/1404)   3.8 E67-E52 N00-N02
      6  591 1308 (N01E31 L= 46)   226    127        0 ( 576/1244- 621/1348) -78.7 E32-E29 S02-N03
      7  282 1168 (S12E52 L= 26)   829    175        1 ( 254/1086- 304/1242) -83.7 E55-E49 S16-S07
      8 1089 1146 (N01W01 L= 79)   187     96        0 (1046/1113-1125/1178) -14.0 E01-W03 N00-N02
      9  243 1065 (S18E55 L= 23)   431     49        0 ( 226/1033- 258/1101) -63.4 E56-E52 S20-S16
     10 1518 1069 (N07W28 L=105)   228     76        0 (1477/1053-1557/1086)   0.0 W24-W30 N05-N07
     11  287  952 (S25E49 L= 29)  1157    279        1 ( 205/ 854- 394/1031)  34.4 E59-E40 S31-S18
     12  741 1010 (S14E16 L= 62)   205     41        1 ( 712/ 993- 767/1031)   0.0 E17-E14 S14-S12
     13 1793  993 (N10W50 L=128)   345     74        1 (1776/ 951-1815/1022) -53.1 W48-W52 N07-N11
     14 1235  798 (S15W17 L= 95)  1164    642        1 ( 987/ 725-1497/ 883)  25.8 W02-W31 S21-S04
     15 1592  663 (S13W44 L=122)   265     46        1 (1564/ 646-1613/ 683)  14.0 W42-W46 S14-S12
     16 1737  564 (S14W62 L=140)   392     72        0 (1707/ 533-1755/ 585) -21.8 W59-W65 S15-S13
     17 1298  517 (S29W31 L=109)   255    103        0 (1255/ 504-1354/ 541)  38.7 W28-W33 S31-S26
  

    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