Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1529 1607 (N30W46 L= 60)   556    100        1 (1477/1570-1575/1635) -33.7 W41-W48 N27-N32
      2 1059 1470 (N26W07 L= 21)   297     89        1 (1015/1450-1109/1487)   0.0 W04-W10 N25-N26
      3  770 1340 (N21E13 L=  1)   231     69        0 ( 749/1314- 799/1378) -68.2 E13-E10 N18-N23
      4 1708  645 (S31W48 L= 62)   322     36        0 (1688/ 620-1730/ 667)  45.0 W46-W50 S32-S29
      5  339  608 (S18E55 L=319)   505     65        1 ( 308/ 587- 378/ 634) -18.4 E58-E51 S19-S16
      6  391  521 (S24E54 L=320)   481    116        1 ( 360/ 467- 415/ 575)  31.0 E61-E50 S27-S20
      7 1695  413 (S49W62 L= 76)   633     66        0 (1668/ 384-1718/ 442) -16.7 W53-W73 S51-S45
      8 1631  390 (S50W53 L= 67)   455     53        0 (1601/ 371-1659/ 405) -12.1 W46-W60 S51-S48
      9 1504  354 (S51W37 L= 51)   624     74        1 (1481/ 315-1530/ 397) -36.0 W32-W43 S55-S47
     10 1553  343 (S53W45 L= 59)   516     59        0 (1535/ 319-1577/ 371) -26.6 W40-W51 S55-S50
  

    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