Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  964 1718 (N38E25 L=131)  1516    215        7 ( 889/1626-1030/1780) -31.0 E34-E18 N31-N42
      2  552 1712 (N27E56 L=100)   393     55        0 ( 531/1679- 573/1736) -18.4 E60-E51 N25-N28
      3 1043 1622 (N33E15 L=141)   659    119        1 (1015/1565-1071/1680) -49.4 E17-E11 N30-N37
      4 1613 1596 (N50W33 L=189)  1189    150        3 (1546/1554-1693/1627)   3.2 W25-W43 N47-N52
      5  438 1317 (N01E44 L=113)  1394    212        4 ( 400/1231- 492/1406) -48.8 E47-E38 S03-N05
      6 1313 1317 (N21W10 L=166)   454    168        0 (1248/1252-1361/1371) -29.1 W05-W14 N17-N22
      7  566  955 (S17E26 L=130)   222     45        1 ( 547/ 933- 592/ 977)   0.0 E27-E24 S18-S16
      8 1721  835 (N03W49 L=205)   290     92        0 (1702/ 789-1742/ 876) -33.7 W46-W52 N01-N05
      9 1060  820 (S14W07 L=163)   941    317        0 ( 923/ 750-1195/ 867)  42.7 W00-W13 S20-S08
     10 1512  788 (S05W35 L=191)   303    120        0 (1480/ 730-1531/ 829) -38.7 W31-W37 S07-S02
     11 1187  593 (S24W21 L=177)  1118    257        1 (1081/ 572-1319/ 630)  18.4 W13-W29 S26-S20
     12  470  561 (S45E31 L=126)   419     74        2 ( 434/ 531- 500/ 587)  45.0 E34-E27 S47-S42
     13 1545  360 (S28W61 L=217)   402     65        0 (1532/ 335-1568/ 391) -14.0 W57-W65 S29-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