Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1020 1547 (N38W04 L=301)   516    129        0 ( 956/1510-1084/1578)  18.4 E01-W08 N36-N40
      2 1773 1465 (N25W67 L=  5)  1422    134        3 (1730/1414-1812/1527)  45.0 W61-W76 N21-N28
      3 1839 1424 (N21W76 L= 13)  1162     58        0 (1820/1396-1862/1457)  76.0 W71-W81 N19-N23
      4 1175 1163 (N12W10 L=308)   386    149        0 (1108/1129-1238/1183) -15.9 W06-W13 N09-N13
      5  190 1109 (N12E64 L=234)   613     80        0 ( 165/1071- 216/1151) -29.7 E67-E59 N10-N14
      6  860  890 (S02E11 L=287)  1978    484        1 ( 786/ 682- 986/1102) -85.4 E14-E03 S15-N10
      7  529 1051 (N09E32 L=266)   326    102        0 ( 511/1005- 560/1090)  63.4 E33-E29 N05-N11
      8 1673  968 (S04W43 L=340)   785    205        0 (1593/ 924-1738/1045)  18.4 W36-W49 S05-N00
      9 1692  832 (S12W45 L=342)   479    105        0 (1669/ 780-1729/ 872)  51.3 W42-W47 S15-S10
     10  398  668 (S15E46 L=251)   617    120        1 ( 362/ 624- 431/ 721) -56.3 E48-E44 S17-S11
     11 1464  589 (S26W28 L=325)   519     98        0 (1424/ 554-1505/ 627)  38.7 W24-W30 S28-S23
     12  397  596 (S20E49 L=249)   404     81        0 ( 369/ 574- 421/ 622) -26.6 E51-E46 S21-S17
     13  680  477 (S29E29 L=269)   408    105        0 ( 663/ 424- 697/ 532)  58.0 E31-E25 S32-S24
     14 1246  346 (S43W13 L=311)  1767    116        5 (1188/ 285-1299/ 409) -74.7 W08-W17 S48-S37
     15  898  289 (S45E17 L=280)   428     56        0 ( 875/ 262- 918/ 319)  45.0 E20-E14 S47-S42
     16  703  252 (S48E39 L=259)  2718    290        4 ( 572/ 212- 837/ 295)   7.1 E63-E23 S52-S43
  

    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