Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1548 1694 (N58W29 L=212)   709     85        1 (1515/1658-1572/1719)  53.1 W24-W34 N53-N61
      2  408 1653 (N21E68 L=116)  1088    116        1 ( 372/1605- 440/1695) -74.1 E74-E64 N17-N24
      3  348 1584 (N16E68 L=115)  2411    196        6 ( 308/1511- 383/1660) -65.6 E81-E59 N10-N21
      4  702 1611 (N25E37 L=146)   453    109        0 ( 661/1563- 731/1655) -82.9 E38-E36 N20-N28
      5 1416 1569 (N41W13 L=196)   469    109        1 (1376/1522-1460/1605) -20.6 W09-W17 N39-N43
      6 1579 1329 (N29W29 L=213)   342     76        1 (1557/1307-1612/1363) -21.8 W27-W32 N27-N31
      7  228 1015 (S21E57 L=126)   525    107        1 ( 193/ 991- 257/1051) -26.6 E61-E53 S21-S18
      8  599  743 (S29E21 L=163)   237     90        0 ( 558/ 729- 646/ 756)  18.4 E23-E17 S30-S28
      9 1110  699 (S19W13 L=196)   270    106        0 (1096/ 651-1130/ 744) -51.3 W11-W15 S21-S16
     10  637  460 (S48E14 L=169)   506     93        2 ( 601/ 441- 677/ 489)  -7.1 E17-E09 S50-S44
     11  643  397 (S53E12 L=171)   455     51        0 ( 627/ 365- 662/ 426)  80.5 E14-E10 S56-S50
     12  687  313 (S60E06 L=177)   512     64        0 ( 661/ 282- 720/ 346)  63.4 E08-E03 S63-S55
  

    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