Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  840 1702 (N51W10 L=279)   402     69        2 ( 804/1683- 880/1717) -26.6 W07-W13 N49-N53
      2 1136 1412 (N23W17 L=286)   558    132        1 (1102/1362-1202/1458)  20.6 W13-W22 N20-N24
      3  942 1371 (N26W04 L=273)   741    265        2 ( 790/1340-1034/1401) -18.4 E05-W10 N23-N29
      4  229 1331 (N42E58 L=211)   854    112        1 ( 167/1307- 283/1347)  -8.1 E70-E49 N39-N43
      5 1721  888 (S22W41 L=310)   705    145        0 (1687/ 825-1743/ 961)  69.4 W37-W43 S25-S17
      6 1500  878 (S16W24 L=293)   454    141        0 (1457/ 845-1551/ 933) -54.5 W22-W27 S18-S11
      7  319  779 (N06E52 L=217)   667     90        2 ( 284/ 742- 364/ 811)  -7.1 E55-E47 N04-N08
      8  332  655 (S01E55 L=214)   787    123        2 ( 307/ 597- 359/ 704)  51.3 E57-E51 S04-N01
      9  978  526 (S23E16 L=252)   851    288        1 ( 853/ 486-1078/ 587) -15.9 E22-E08 S25-S18
     10 1296  473 (S36W02 L=270)   311     35        1 (1274/ 449-1318/ 497)  45.0 W00-W03 S37-S33
     11  728  382 (S26E38 L=230)   489    121        0 ( 713/ 320- 746/ 443)  52.1 E41-E34 S30-S21
     12 1222  383 (S40E06 L=262)   251     30        0 (1204/ 361-1239/ 406)  56.3 E07-E04 S41-S37
  

    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