Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  839 1609 (N26E28 L=223)  1906    425        3 ( 728/1501- 967/1690)  71.6 E31-E22 N16-N34
      2 1038 1599 (N30E15 L=236)  1272    194        4 ( 979/1522-1086/1667) -26.6 E22-E10 N25-N33
      3 1492 1333 (N26W21 L=272)  2600    519        4 (1358/1098-1562/1502) -51.0 W09-W28 N13-N35
      4  153 1141 (S16E68 L=183)  3044    302        4 (  81/1027- 206/1269)  34.3 E86-E63 S23-S08
      5  407  981 (S21E38 L=213)   779    199        1 ( 339/ 918- 500/1025)  36.9 E43-E31 S25-S16
      6 1770  957 (N12W50 L=301)   418     93        1 (1756/ 911-1782/1006) -80.5 W48-W51 N09-N15
      7  843  853 (S18E06 L=245)   381    165        0 ( 812/ 798- 897/ 898)  63.4 E07-E03 S20-S14
      8  421  858 (S28E35 L=216)   450    103        0 ( 377/ 845- 479/ 876)   6.3 E39-E30 S28-S26
      9 1301  791 (S10W22 L=273)   394    139        0 (1243/ 758-1363/ 813)   7.1 W18-W26 S11-S08
     10  736  396 (S51E01 L=250)  1704    270        4 ( 631/ 316- 831/ 495) -17.7 E12-W10 S54-S45
     11  941  246 (S55W26 L=277)   375     24        0 ( 919/ 227- 966/ 262)  76.0 W23-W28 S57-S53
  

    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