Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  907 1783 (N60W13 L=287)   369     29        1 ( 886/1760- 925/1805)  26.6 W10-W16 N57-N62
      2 1486 1607 (N31W49 L=323)  1843    257        2 (1400/1572-1609/1648) -28.1 W40-W59 N27-N36
      3  600 1464 (N41E23 L=251)  1014    254        1 ( 490/1442- 720/1486) -13.4 E34-E13 N38-N44
      4  225 1448 (N44E73 L=201)   654     36        1 ( 208/1422- 245/1470) -14.0 E82-E67 N42-N45
      5 1751  671 (S31W47 L=321)  2765    366        5 (1636/ 568-1862/ 835)  35.8 W35-W59 S36-S22
      6  227  627 (S07E70 L=204)  1198    104        3 ( 197/ 573- 252/ 666)  80.5 E74-E64 S10-S04
      7 1487  648 (S25W23 L=297)   346     87        0 (1441/ 633-1526/ 666)   0.0 W19-W26 S26-S23
      8  643  459 (S22E37 L=237)   325     67        0 ( 618/ 431- 682/ 486) -68.2 E38-E34 S24-S19
      9 1055  273 (S44E16 L=258)  3389    482        6 ( 860/ 182-1279/ 386)  10.5 E39-W04 S48-S38
  

    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