Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  721 1696 (N26E38 L=119)   738    203        2 ( 659/1654- 793/1743) -12.1 E45-E31 N23-N28
      2  376 1634 (N16E63 L= 94)   857    196        1 ( 342/1589- 400/1672) -56.3 E67-E60 N12-N18
      3  480 1592 (N15E51 L=107)   345     68        1 ( 456/1557- 508/1613)  76.0 E52-E48 N12-N16
      4  452 1481 (N07E48 L=109)  1086    353        1 ( 336/1419- 585/1534)  27.9 E55-E37 N01-N10
      5  934 1447 (N15E16 L=142)   260     91        0 ( 888/1422- 979/1475)  -9.5 E18-E12 N13-N15
      6  834  786 (S25E06 L=152)   217     50        0 ( 820/ 753- 851/ 821) -56.3 E07-E04 S26-S23
      7  196  739 (S41E57 L=100)   473     96        1 ( 172/ 723- 226/ 758)   9.5 E59-E53 S42-S39
      8 1751  623 (S09W60 L=218)   624     96        1 (1711/ 594-1780/ 653)  45.0 W57-W62 S11-S06
      9 1599  355 (S27W63 L=220)   557     76        0 (1561/ 340-1629/ 374)  45.0 W59-W65 S29-S25
     10  860  210 (S65W24 L=181)   611     69        2 ( 832/ 188- 892/ 235)  18.4 W19-W27 S67-S62
  

    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