Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1415 1467 (N21W32 L=197)   378     85        2 (1387/1422-1439/1511) -90.0 W30-W33 N18-N24
      2  846 1480 (N30E06 L=159)   232     85        0 ( 811/1445- 883/1504) -51.3 E07-E03 N27-N32
      3  229 1214 (N21E56 L=108)   990    125        1 ( 187/1170- 290/1251)  15.3 E60-E49 N19-N23
      4 1720  972 (S12W46 L=210)   472    155        4 (1669/ 949-1781/1010)  -5.7 W41-W51 S13-S09
      5 1393  844 (S15W20 L=185)  1121    382        1 (1342/ 666-1433/ 992)  81.5 W17-W22 S26-S06
      6 1791  903 (S17W52 L=217)   743     98        1 (1760/ 866-1841/ 939) -12.5 W49-W58 S19-S14
      7 1765  832 (S21W50 L=214)   870    195        0 (1702/ 767-1830/ 874) -25.0 W43-W58 S25-S17
      8 1285  748 (S20W13 L=177)   210    115        0 (1250/ 705-1302/ 803)  56.3 W10-W14 S22-S16
      9  260  538 (S19E70 L= 94)  1206     87        3 ( 232/ 502- 280/ 591)  45.0 E76-E64 S21-S16
     10 1671  421 (S49W54 L=219)   932     92        2 (1633/ 403-1719/ 451)  -9.0 W46-W65 S51-S46
  

    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