Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  597 1661 (N20E48 L=133)  1322    235        4 ( 553/1608- 660/1698)  74.1 E52-E42 N15-N22
      2  248 1467 (N01E71 L=110)   563     78        0 ( 231/1440- 264/1492) -15.9 E74-E66 N00-N02
      3 1740 1179 (N25W41 L=221)   531     75        0 (1702/1144-1771/1217) -15.9 W37-W44 N22-N26
      4  896  895 (S16E04 L=177)   606    178        0 ( 812/ 844- 984/ 947)  55.0 E07-W00 S21-S11
      5 1836  921 (N13W56 L=237)   452     70        0 (1814/ 894-1857/ 946)  63.4 W54-W58 N11-N15
      6  236  874 (S35E51 L=129)   460     63        1 ( 216/ 846- 253/ 908)  63.4 E53-E49 S36-S32
      7 1613  730 (S04W44 L=225)  1578    337        3 (1572/ 582-1642/ 905) -59.5 W37-W48 S12-N05
      8 1175  659 (S22W20 L=200)   263     88        1 (1127/ 638-1220/ 681)   0.0 W16-W23 S22-S20
      9  707  561 (S42E06 L=174)  1025    284        2 ( 575/ 537- 849/ 592)  15.9 E17-W04 S45-S38
     10 1470  366 (S28W53 L=234)   443     98        0 (1421/ 347-1513/ 383)  38.7 W50-W55 S29-S25
  

    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