Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1405 1573 (N29W42 L=157)   371     48        1 (1384/1550-1432/1594) -56.3 W39-W44 N27-N30
      2  286 1068 (N24E52 L= 63)   430     71        1 ( 267/1028- 310/1102)  76.0 E53-E48 N21-N26
      3  345 1064 (N23E46 L= 69)   359     72        0 ( 325/1032- 363/1097)  45.0 E47-E43 N21-N24
      4 1577  920 (S12W31 L=145)   445    103        0 (1525/ 906-1628/ 934) -14.0 W27-W35 S14-S10
      5 1762  893 (S20W46 L=161)   388     69        0 (1744/ 853-1781/ 928)  36.9 W44-W48 S21-S17
      6 1042  559 (S21E10 L=104)  2007    482        2 ( 876/ 413-1163/ 726)  26.6 E25-W01 S27-S13
      7 1290  565 (S27W06 L=120)   588    216        1 (1238/ 494-1317/ 661)  51.3 W00-W09 S30-S20
      8  766  404 (S25E33 L= 81)  1885    416        2 ( 701/ 330- 884/ 545) -86.6 E38-E25 S32-S15
      9  679  278 (S33E49 L= 65)   647    108        1 ( 654/ 243- 704/ 321)  50.2 E54-E45 S35-S29
  

    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