Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1305 1765 (N57W20 L= 20)   347     81        1 (1253/1751-1343/1781)  14.0 W13-W25 N54-N58
      2  467 1661 (N37E58 L=302)  2007    163        6 ( 427/1590- 504/1739) -17.9 E80-E48 N32-N42
      3  676 1604 (N35E34 L=327)  1595    255        3 ( 578/1563- 791/1639)  10.0 E41-E23 N31-N38
      4  663 1411 (N22E28 L=332)   562    201        0 ( 617/1349- 691/1512)  90.0 E30-E26 N17-N28
      5  577 1445 (N23E35 L=325)   247     77        0 ( 563/1412- 587/1493) -78.7 E36-E34 N21-N26
      6 1723 1349 (N27W51 L= 51)   596    160        0 (1680/1284-1748/1426) -71.6 W48-W52 N23-N32
      7 1267 1393 (N27W13 L= 13)   438    105        0 (1228/1355-1326/1421)  29.7 W10-W17 N23-N28
      8 1301 1277 (N19W16 L= 16)   204     83        1 (1257/1259-1343/1300)  -9.5 W12-W18 N18-N20
      9 1485 1190 (N15W28 L= 28)   256     93        0 (1437/1179-1530/1210)   0.0 W25-W31 N14-N16
     10 1798 1062 (N10W55 L= 55)   463    116        0 (1777/1003-1818/1115) -63.4 W52-W57 N06-N12
     11 1587 1009 (N05W37 L= 37)   259     58        0 (1572/ 974-1600/1043)  90.0 W35-W37 N02-N07
     12 1187 1031 (N03W10 L= 10)   197     59        0 (1154/1019-1228/1042)   0.0 W07-W12 N02-N04
     13 1463  984 (N03W28 L= 28)   446    147        0 (1412/ 935-1528/1011)  20.6 W24-W32 N00-N04
     14  178  926 (S12E64 L=296)  1611    186        4 ( 141/ 846- 209/ 995)  78.7 E69-E59 S17-S07
     15  258  854 (S16E54 L=306)   483     63        1 ( 238/ 815- 282/ 887)  68.2 E57-E51 S18-S13
     16 1777  695 (S13W59 L= 59)  3547    385        5 (1651/ 582-1872/ 854)  48.4 W50-W67 S20-S02
     17  730  755 (S18E16 L=344)   235     69        0 ( 697/ 742- 775/ 769)  21.8 E18-E13 S18-S16
     18  455  584 (S32E39 L=321)  1180    200        1 ( 390/ 507- 542/ 636)  31.4 E48-E30 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