Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1596 1285 (N07W42 L=152)   273    101        2 (1575/1258-1613/1324)  56.3 W39-W43 N05-N08
      2  218 1207 (N34E55 L= 55)   427     71        1 ( 196/1183- 243/1236) -90.0 E57-E52 N32-N36
      3  239  976 (N20E56 L= 54)   335    116        1 ( 208/ 953- 280/ 999) -23.2 E58-E51 N18-N21
      4 1790  741 (S31W46 L=156)  2109    409        7 (1682/ 621-1896/ 879)  20.6 W35-W59 S36-S26
      5 1358  706 (S19W12 L=122)   289    109        0 (1298/ 679-1399/ 726) -39.8 W08-W14 S22-S17
      6  937  634 (S13E15 L= 95)   444    134        2 ( 890/ 591- 983/ 688)  15.9 E17-E10 S15-S10
      7 1589  638 (S31W27 L=137)   461     71        1 (1560/ 600-1611/ 679)  45.0 W24-W29 S33-S28
      8  990  573 (S18E13 L= 97)   224     98        0 ( 950/ 549-1012/ 610)  11.3 E16-E11 S19-S15
      9 1116  560 (S22E06 L=104)   231    101        1 (1083/ 526-1158/ 604) -74.1 E07-E04 S25-S18
     10  301  573 (S05E64 L= 46)   441     68        1 ( 282/ 545- 319/ 598)  38.7 E66-E61 S07-S03
  

    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