Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  634 1708 (N28E48 L=139)  1444    249        4 ( 553/1654- 706/1741)  45.0 E52-E42 N22-N30
      2  367 1581 (N14E63 L=123)  3123    442       15 ( 309/1485- 440/1668)  81.9 E69-E57 N07-N21
      3  487 1575 (N16E51 L=135)   353     83        1 ( 466/1548- 510/1599)  90.0 E53-E49 N14-N17
      4 1034  823 (S15W06 L=192)   273     77        1 (1003/ 795-1070/ 847)  53.1 W04-W07 S16-S12
      5 1627  736 (S04W45 L=232)   683    112        2 (1602/ 676-1648/ 789) -63.4 W42-W47 S07-S01
      6 1273  715 (S15W23 L=209)   440    143        0 (1207/ 646-1331/ 772)  74.7 W20-W24 S20-S09
      7  736  558 (S39E06 L=180)  1140    265        3 ( 595/ 532- 905/ 584)  17.7 E17-W05 S42-S34
      8  320  528 (S54E49 L=138)   757    228        1 ( 268/ 509- 374/ 545)  14.9 E55-E40 S55-S50
      9 1224  415 (S34W32 L=219)   386    131        1 (1181/ 392-1294/ 444)  20.6 W28-W36 S35-S30
     10 1568  392 (S25W58 L=245)  1949    449        1 (1457/ 334-1652/ 426)   7.8 W48-W71 S27-S22
  

    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