Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  777 1751 (N60W11 L= 74)   462     94        0 ( 735/1713- 808/1782)  14.0 W04-W16 N58-N62
      2  635 1688 (N60E08 L= 55)   610     60        2 ( 596/1666- 667/1709) -29.1 E13-E04 N57-N63
      3  504 1645 (N62E27 L= 36)   390     21        0 ( 481/1626- 525/1663) -26.6 E31-E23 N59-N63
      4 1165 1560 (N32W27 L= 90)   810    332        2 (1022/1493-1308/1603) -49.1 W19-W32 N24-N39
      5  348 1536 (N57E50 L= 14)   432     37        0 ( 332/1512- 366/1562) -26.6 E55-E46 N55-N59
      6  643 1449 (N42E13 L= 50)   558    178        0 ( 608/1380- 685/1523) -80.5 E15-E11 N35-N47
      7  985 1384 (N27W08 L= 72)   600    202        2 ( 882/1362-1087/1412) -32.5 W02-W14 N23-N31
      8  757 1163 (N21E12 L= 52)  2048    510        7 ( 600/1004- 898/1269) -57.0 E20-E04 N08-N29
      9  531 1033 (N19E30 L= 34)   572    192        0 ( 451/ 981- 619/1067) -36.0 E35-E24 N14-N22
     10 1752  900 (S22W44 L=107)   849    212        0 (1691/ 819-1808/ 993)  23.2 W36-W50 S25-S18
     11 1539  833 (S19W25 L= 89)   273     81        1 (1512/ 795-1574/ 867) -71.6 W23-W27 S22-S16
     12 1633  800 (S24W32 L= 95)   504    146        1 (1592/ 777-1684/ 856) -56.3 W28-W35 S26-S20
     13 1477  761 (S21W19 L= 83)   289    118        0 (1439/ 718-1512/ 810)  26.6 W16-W22 S22-S19
     14  966  512 (S21E18 L= 45)   663    161        2 ( 917/ 446-1024/ 560)  32.0 E23-E14 S24-S19
     15  430  463 (S12E57 L=  6)   530     81        1 ( 405/ 419- 453/ 500)  63.4 E59-E55 S15-S09
     16  837  351 (S29E33 L= 30)   559    223        1 ( 770/ 296- 868/ 449)  34.7 E41-E27 S31-S22
  

    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