Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1100 1603 (N30E05 L=284)   560    140        0 (1059/1540-1141/1655)  77.5 E06-E03 N25-N34
      2  467 1553 (N18E47 L=242)  1177    149        1 ( 398/1525- 545/1573)   0.0 E53-E40 N15-N19
      3 1051 1452 (N19E06 L=284)   251     53        0 (1035/1420-1069/1490) -63.4 E06-E03 N17-N21
      4  229 1382 (N05E63 L=226)   420     51        0 ( 214/1349- 249/1408)  71.6 E65-E61 N03-N06
      5 1682 1225 (N17W39 L=328)   594    179        1 (1588/1193-1747/1249)   8.7 W31-W44 N16-N18
      6  808 1201 (N00E15 L=274)   402    168        0 ( 713/1169- 873/1237)  -5.7 E21-E11 N00-N01
      7 1912 1039 (N12W64 L=354)  3728    237       10 (1869/ 934-1961/1152)  85.9 W58-W73 N05-N19
      8 1835 1086 (N13W54 L=344)   704     94        3 (1811/1038-1867/1128)  63.4 W51-W57 N09-N15
      9 1664 1071 (N08W39 L=329)   252     57        2 (1646/1036-1681/1102) -76.0 W37-W40 N05-N09
     10 1723 1036 (N07W44 L=334)   390     54        0 (1702/ 999-1746/1062)  63.4 W42-W46 N04-N08
     11  841  804 (S22E08 L=282)   227     63        0 ( 815/ 772- 868/ 838)  68.2 E08-E06 S24-S19
     12 1013  696 (S26W05 L=294)  1392    408        1 ( 842/ 649-1218/ 743)  16.9 E05-W18 S31-S21
     13  899  550 (S37W00 L=290)   446    179        0 ( 833/ 522- 974/ 609) -10.3 E05-W06 S38-S34
     14 1608  489 (S26W54 L=343)  3619    419        4 (1490/ 412-1748/ 600)  -1.5 W40-W78 S29-S20
     15  846  441 (S46E02 L=288)  1397    335        2 ( 695/ 400-1000/ 503)   0.0 E14-W11 S47-S42
     16 1226  354 (S43W32 L=321)  1160    304        0 (1101/ 298-1373/ 413)   0.0 W19-W48 S45-S41
  

    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