Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  968 1887 (N61E39 L=117)  1150    107        2 ( 934/1852-1009/1933) -11.7 E60-E30 N56-N65
      2  642 1580 (N27E39 L=118)   453    154        0 ( 599/1558- 703/1624)  -5.7 E44-E34 N25-N29
      3  317 1552 (N20E70 L= 86)   776     99        2 ( 277/1529- 356/1570)   4.4 E78-E63 N17-N20
      4  539 1470 (N18E42 L=115)   486     78        1 ( 518/1440- 565/1505) -33.7 E44-E38 N16-N20
      5 1845 1313 (N32W65 L=221)  2005    220        2 (1766/1245-1902/1402) -23.6 W54-W77 N28-N37
      6 1549 1256 (N23W31 L=187)   312     91        0 (1532/1208-1563/1297) -80.5 W29-W32 N19-N25
      7 1809  924 (N07W58 L=214)   457     98        0 (1787/ 873-1829/ 966)  63.4 W56-W59 N03-N09
      8  227  853 (S25E58 L= 98)   877    148        0 ( 206/ 782- 262/ 916)  55.0 E61-E53 S29-S19
      9  362  807 (S25E43 L=113)   343     97        1 ( 325/ 787- 407/ 836)  -7.1 E46-E38 S26-S23
     10 1032  620 (S25W08 L=165)  2559    580        5 ( 811/ 513-1246/ 753)  43.7 E05-W18 S35-S12
     11 1533  626 (S16W42 L=198)   578    152        0 (1485/ 579-1598/ 661)  45.0 W39-W45 S19-S12
  

    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