Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1639 1537 (N16W55 L=104)  2145    337        3 (1506/1428-1776/1635) -41.4 W49-W68 N08-N25
      2 1339 1548 (N25W34 L= 82)  2511    578        4 (1107/1434-1579/1602) -32.2 W18-W46 N12-N32
      3 1217 1424 (N21W22 L= 70)   259     68        0 (1177/1413-1255/1437) -11.3 W18-W23 N19-N22
      4  284 1082 (N23E48 L=  0)  1191    225        2 ( 188/1002- 357/1134) -36.3 E58-E42 N16-N27
      5 1592 1082 (S07W35 L= 83)   344     74        1 (1560/1056-1627/1107) -26.6 W32-W36 S08-S05
      6 1589 1005 (S11W33 L= 81)   424    130        1 (1523/ 981-1638/1045) -14.0 W28-W36 S12-S09
      7 1222  920 (S07W08 L= 57)   252     56        1 (1186/ 895-1243/ 939) -45.0 W06-W10 S08-S05
      8 1119  848 (S08W01 L= 49)   883    270        4 (1002/ 792-1180/ 916) -41.2 E05-W05 S12-S05
      9  706  779 (S03E24 L= 25)  2323    601        6 ( 523/ 606- 908/ 910) -48.2 E33-E15 S16-N05
     10  229  849 (N10E58 L=351)   749    141        3 ( 178/ 823- 302/ 884) -18.4 E62-E50 N08-N12
     11 1011  783 (S10E06 L= 42)   324    100        2 ( 971/ 760-1056/ 819) -21.8 E08-E03 S11-S07
     12 1082  740 (S14E03 L= 45)   344    138        2 (1021/ 715-1141/ 777)  -9.5 E06-W00 S14-S11
     13  442  520 (S12E51 L=358)  5086    602       15 ( 365/ 327- 523/ 713)  46.2 E71-E41 S24-N00
     14 1262  599 (S27W04 L= 52)  1573    314        3 (1157/ 560-1397/ 675) -39.8 E02-W12 S29-S19
     15  989  370 (S35E21 L= 28)  4077    678        7 ( 843/ 150-1100/ 660)  43.3 E48-E10 S48-S14
     16  304  562 (S07E62 L=347)   423     92        1 ( 281/ 539- 339/ 600)  12.5 E65-E56 S08-S04
     17  572  541 (S14E40 L=  8)   565    128        2 ( 536/ 476- 608/ 579)  32.0 E44-E36 S17-S11
  

    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