Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  898 1611 (N36E19 L= 77)  1357    284        3 ( 849/1501- 936/1705) -51.8 E26-E14 N28-N42
      2  566 1572 (N28E42 L= 54)  1410    448        2 ( 496/1440- 631/1696) -33.1 E55-E32 N20-N36
      3 1702 1556 (N46W52 L=147)  1375    284        2 (1663/1515-1743/1600)  32.7 W45-W59 N41-N50
      4  281 1449 (N16E64 L= 32)  3803    533       11 ( 230/1295- 336/1589) -52.1 E72-E54 N07-N25
      5  924 1357 (N19E11 L= 85)   610    221        0 ( 856/1269- 958/1424) -84.3 E14-E08 N13-N23
      6 1148 1340 (N21W04 L= 99)   264    134        0 (1116/1281-1166/1392) -50.2 W00-W05 N17-N24
      7  681 1296 (N12E25 L= 70)   460    134        1 ( 622/1253- 720/1345) -26.6 E30-E22 N09-N14
      8  495 1212 (N05E37 L= 59)   358    121        1 ( 437/1189- 543/1238)  29.7 E40-E33 N02-N06
      9 1557 1144 (N15W33 L=129)   350     68        1 (1534/1114-1584/1182) -36.9 W31-W35 N12-N16
     10 1301 1069 (N07W16 L=112)   562    210        1 (1256/ 960-1335/1162) -57.5 W12-W19 N00-N11
     11 1523 1059 (N09W31 L=127)   282    196        2 (1486/1036-1586/1089) -14.0 W28-W36 N07-N10
     12  152 1048 (S10E68 L= 28)   598    179        2 ( 116/1028- 180/1062)  10.3 E74-E63 S10-S08
     13  743  924 (S09E16 L= 80)   205     65        1 ( 715/ 904- 776/ 942) -26.6 E17-E13 S10-S08
     14 1164  802 (S11W12 L=107)   289    117        0 (1135/ 777-1219/ 846) -26.6 W09-W15 S12-S08
     15 1557  685 (S13W40 L=136)   605    122        1 (1518/ 643-1583/ 722) -29.7 W36-W43 S15-S10
     16 1237  468 (S32W24 L=120)   262     62        1 (1204/ 448-1258/ 493)  14.0 W21-W25 S33-S29
     17  776  454 (S40E10 L= 86)   547    105        2 ( 722/ 438- 824/ 477)  12.5 E14-E05 S42-S37
  

    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