Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  990 1749 (N45E06 L= 99)   353     39        0 ( 965/1728-1016/1765)  45.0 E08-E03 N43-N46
      2 1143 1717 (N43W07 L=111)   374     69        0 (1104/1704-1184/1731)  -8.1 W03-W10 N41-N44
      3  581 1525 (N27E33 L= 71)   268     76        0 ( 546/1491- 607/1549) -33.7 E36-E30 N24-N28
      4 1149 1508 (N27W07 L=111)   212     68        0 (1122/1477-1175/1537) -53.1 W04-W08 N24-N28
      5  479 1477 (N23E40 L= 64)   456     91        0 ( 461/1420- 502/1519) -81.9 E41-E38 N19-N26
      6 1839 1362 (N22W64 L=169)  1735    157        2 (1810/1307-1898/1409)  48.8 W59-W74 N17-N25
      7 1590 1249 (N12W36 L=141)   284     53        1 (1558/1222-1614/1266)  21.8 W33-W38 N10-N13
      8 1668 1193 (N09W42 L=147)   294     71        1 (1637/1170-1702/1214)  26.6 W39-W45 N07-N10
      9 1704 1137 (N06W45 L=150)   443     80        2 (1677/1107-1743/1173) -21.8 W42-W48 N04-N08
     10  221 1137 (N02E58 L= 47)   521     90        1 ( 176/1122- 262/1154)  -5.7 E63-E53 N01-N03
     11 1667 1022 (S01W42 L=147)   748    125        0 (1643/ 965-1690/1085)  90.0 W40-W44 S04-N02
     12  178  961 (S08E63 L= 41)  2098    166        7 ( 133/ 879- 236/1015)  35.0 E70-E55 S13-S04
     13  875  974 (S07E09 L= 96)   197     52        0 ( 844/ 954- 904/ 991) -33.7 E10-E07 S08-S06
     14  322  816 (S17E49 L= 55)   506     73        0 ( 307/ 770- 343/ 851)  68.2 E51-E47 S20-S15
     15  768  788 (S19E16 L= 89)   435     92        1 ( 748/ 744- 794/ 846)  71.6 E16-E13 S21-S15
     16  391  748 (S22E44 L= 60)   280     70        0 ( 368/ 715- 414/ 776) -45.0 E45-E42 S23-S19
  

    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