Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  877 1758 (N46E29 L=145)   386     57        2 ( 850/1736- 907/1778)  78.7 E32-E26 N43-N48
      2 1020 1679 (N42E14 L=161)   422     63        0 ( 986/1657-1051/1703)   0.0 E16-E11 N40-N43
      3  483 1549 (N23E49 L=126)   336     83        1 ( 462/1510- 510/1593) -80.5 E50-E46 N20-N26
      4  281 1455 (N14E64 L=110)  5569    593       14 ( 180/1307- 379/1550)  56.3 E75-E55 N03-N21
      5  408 1351 (N10E47 L=127)   711    165        1 ( 380/1314- 461/1393)  78.7 E49-E42 N07-N12
      6  899 1291 (N14E12 L=163)   221     62        1 ( 877/1267- 924/1315) -18.4 E13-E10 N12-N15
      7  484 1239 (N04E38 L=136)   287     87        0 ( 472/1193- 500/1288) -63.4 E39-E36 N01-N07
      8  405 1015 (S10E40 L=134)  1408    272        3 ( 375/ 910- 460/1136) -59.7 E43-E34 S16-S03
      9  857 1104 (N02E11 L=163)   391     80        0 ( 819/1082- 898/1132) -11.3 E13-E08 N00-N02
     10 1059 1018 (S00W02 L=177)   380    105        1 (1008/ 996-1119/1036)  31.0 W00-W05 S02-N01
     11 1811  738 (S05W63 L=238)   949    156        2 (1789/ 666-1823/ 802) -63.4 W59-W67 S09-S01
     12 1290  699 (S16W22 L=197)   217     90        1 (1267/ 658-1331/ 725)  59.0 W21-W24 S18-S13
     13 1409  407 (S33W41 L=215)  1098    288        0 (1298/ 378-1516/ 452)   3.0 W30-W49 S34-S29
  

    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