Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  460 1663 (N37E59 L=302)  2066    157        4 ( 423/1596- 497/1742) -15.7 E81-E49 N32-N42
      2  661 1610 (N35E35 L=327)  1675    245        5 ( 564/1572- 778/1644)   6.3 E43-E24 N31-N38
      3  563 1472 (N25E37 L=324)   408    116        0 ( 549/1420- 589/1527) -81.9 E39-E35 N21-N28
      4  627 1402 (N21E30 L=331)   599    219        0 ( 550/1355- 671/1485)  56.3 E35-E28 N17-N26
      5 1242 1390 (N26W11 L= 13)   351    141        0 (1210/1342-1303/1420)  39.8 W09-W15 N22-N28
      6 1724 1297 (N24W50 L= 51)   366     69        1 (1707/1258-1739/1333) -78.7 W48-W50 N21-N26
      7 1772 1098 (N12W52 L= 54)   411     73        0 (1755/1057-1788/1135) -68.2 W50-W54 N09-N14
      8  190  975 (S09E62 L=300)   478     80        2 ( 170/ 951- 222/1016) -31.0 E64-E57 S11-S06
      9 1454  992 (N03W27 L= 29)   466    155        0 (1394/ 958-1525/1015)  18.4 W23-W32 N00-N04
     10  168  898 (S14E66 L=296)   914    125        2 ( 140/ 850- 195/ 949) -71.6 E70-E61 S17-S11
     11  250  869 (S15E55 L=306)   416     59        1 ( 232/ 841- 274/ 894)  38.7 E57-E52 S17-S13
     12 1760  698 (S13W57 L= 59)  3138    390        3 (1629/ 581-1850/ 862)  55.6 W48-W64 S21-S02
     13  436  588 (S32E41 L=321)  1246    226        2 ( 369/ 508- 525/ 643)  28.8 E51-E31 S38-S27
  

    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