Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1368 1642 (N45W05 L=169)   384    104        0 (1310/1620-1410/1659)  12.5 W00-W09 N42-N46
      2  485 1608 (N17E55 L=109)   578     78        3 ( 451/1579- 521/1637)  -6.3 E59-E50 N15-N19
      3  466 1536 (N12E52 L=112)   497     78        3 ( 426/1504- 496/1566)  31.0 E54-E48 N10-N14
      4  746 1504 (N17E30 L=134)   278     50        1 ( 719/1485- 773/1526)  18.4 E31-E28 N15-N18
      5 1661 1298 (N30W35 L=199)  1850    364        3 (1524/1198-1793/1403)  -8.4 W22-W49 N27-N33
      6  187 1073 (S20E63 L=101)   729    123        1 ( 159/1017- 224/1127)  45.0 E67-E59 S23-S16
      7 1105 1059 (N00W04 L=168)   451    146        0 (1066/ 999-1141/1122) -29.7 W00-W07 S02-N02
      8 1138  948 (S05W08 L=172)   423     97        1 (1110/ 908-1168/ 996)  78.7 W07-W09 S07-S02
      9  191  827 (S35E63 L=101)   523    113        0 ( 178/ 785- 207/ 883)  49.4 E66-E60 S38-S31
     10 1045  826 (S14W06 L=170)   348     81        0 (1004/ 805-1087/ 850)  45.0 W04-W08 S16-S11
     11 1455  548 (S19W41 L=206)   389     69        0 (1414/ 528-1490/ 563)   9.5 W38-W44 S20-S17
  

    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