Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1342 1744 (N56W19 L=166)   537     56        0 (1322/1713-1365/1770)  45.0 W16-W23 N53-N59
      2 1452 1744 (N59W35 L=182)   797     91        1 (1410/1724-1501/1764)  13.2 W27-W44 N55-N60
      3  748 1676 (N40E34 L=113)   334     52        1 ( 719/1662- 773/1695)   0.0 E36-E31 N38-N41
      4 1278 1656 (N47W11 L=158)   302     47        0 (1259/1634-1298/1683) -56.3 W09-W13 N44-N48
      5  459 1447 (N20E46 L=100)  1805    371        3 ( 358/1345- 533/1553)  72.6 E51-E41 N11-N27
      6 1027 1486 (N30E06 L=140)   338     80        1 (1009/1442-1040/1527) -71.6 E07-E04 N26-N32
      7  346 1457 (N19E58 L= 88)  1803    133        6 ( 281/1420- 399/1506)  26.6 E65-E51 N15-N22
      8  559 1456 (N21E39 L=108)   575     87        3 ( 525/1424- 600/1488)  45.0 E40-E35 N19-N23
      9 1616 1431 (N34W40 L=187)   804    185        1 (1540/1381-1683/1463) -17.1 W33-W46 N31-N35
     10  549 1389 (N17E37 L=109)   373    101        1 ( 502/1360- 575/1417)  53.1 E39-E35 N14-N19
     11  336 1376 (N14E55 L= 91)   540     96        1 ( 308/1331- 366/1401)  78.7 E58-E52 N10-N15
     12  297  852 (S19E51 L= 96)   989    153        1 ( 258/ 790- 351/ 890)  37.9 E55-E44 S23-S16
     13  785  712 (S21E12 L=135)   274     96        0 ( 741/ 692- 809/ 751)  51.3 E14-E10 S23-S18
     14 1703  645 (S14W55 L=201)   503     98        1 (1676/ 608-1725/ 675)  63.4 W52-W57 S16-S12
     15 1700  521 (S23W62 L=209)   471     40        0 (1672/ 506-1725/ 537)  -6.3 W57-W66 S23-S21
     16 1241  434 (S34W25 L=172)   369    151        0 (1188/ 395-1308/ 476)  53.1 W22-W28 S38-S30
  

    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