Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  574 1694 (N57E38 L= 75)   682     92        0 ( 526/1671- 618/1712) -10.6 E46-E30 N54-N58
      2  397 1310 (N30E43 L= 69)   267    102        0 ( 367/1269- 418/1338) -45.0 E46-E40 N27-N31
      3 1127 1245 (N18W09 L=122)   403    114        0 (1099/1185-1145/1298) -90.0 W08-W10 N14-N21
      4  311 1215 (N24E50 L= 62)   379    106        1 ( 269/1184- 355/1253) -21.8 E55-E45 N22-N26
      5 1584 1138 (N05W37 L=150)   666    184        0 (1549/1051-1609/1209)  83.7 W35-W39 N00-N09
      6 1800  765 (S22W54 L=167)  1705    183        3 (1752/ 687-1838/ 864)  79.7 W49-W58 S27-S16
      7 1782  585 (S35W59 L=171)  1568    132        0 (1752/ 523-1811/ 643) -28.6 W51-W73 S41-S29
      8  408  598 (S15E49 L= 64)   525     85        0 ( 367/ 576- 450/ 619)  12.5 E53-E44 S15-S13
      9 1377  545 (S29W18 L=130)  1023    181        1 (1302/ 505-1455/ 595)  20.0 W12-W23 S31-S26
     10 1146  463 (S31W01 L=113)   656    119        1 (1097/ 422-1201/ 493)  23.2 E03-W04 S33-S29
     11  943  325 (S40E17 L= 95)   374     97        0 ( 895/ 309- 992/ 344)  12.5 E21-E12 S40-S38
     12  729  266 (S42E39 L= 74)   621     70        1 ( 684/ 251- 772/ 284)   0.0 E44-E34 S43-S40
  

    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