Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  776 1801 (N45E33 L=260)  1035    104        1 ( 725/1774- 835/1823)  11.3 E38-E26 N42-N47
      2 1126 1700 (N40E00 L=293)   379    111        0 (1081/1677-1181/1722)   8.1 E03-W04 N37-N42
      3  443 1474 (N18E45 L=248)  1311    217        0 ( 384/1398- 527/1554) -29.1 E54-E36 N13-N23
      4 1689 1292 (N19W43 L=336)   815    189        1 (1625/1242-1775/1348) -18.4 W37-W50 N16-N21
      5 1383 1189 (N08W20 L=314)   214     29        0 (1357/1176-1412/1203)  14.0 W18-W22 N07-N09
      6 1882  886 (S02W66 L=359)  1290    144        2 (1821/ 841-1940/ 919)  20.6 W59-W75 S06-N00
      7 1116  837 (S15W07 L=301)   247     57        0 (1081/ 820-1153/ 857)   0.0 W05-W09 S16-S14
      8  626  694 (S29E24 L=269)  1044    357        0 ( 557/ 557- 733/ 836)  51.7 E31-E16 S38-S19
      9 1765  603 (S21W63 L=357)   997    158        0 (1719/ 577-1821/ 623)  20.0 W58-W71 S23-S18
  

    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