Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  592 1745 (N31E51 L= 64)  1515    208        1 ( 524/1682- 665/1798) -15.6 E65-E40 N26-N34
      2  957 1682 (N33E20 L= 95)   788     96        1 ( 917/1625- 990/1732) -48.8 E23-E15 N28-N36
      3 1059 1503 (N22E08 L=107)   526    143        3 (1018/1454-1112/1573) -49.4 E11-E03 N19-N26
      4  893 1300 (N07E13 L=102)   901    312        2 ( 868/1174- 920/1436) -60.3 E16-E08 N00-N14
      5 1713 1395 (N32W41 L=156)   344     95        0 (1676/1350-1735/1424) -31.0 W37-W43 N28-N33
      6 1626 1356 (N26W33 L=147)   615    109        3 (1598/1305-1664/1413) -50.2 W30-W36 N23-N29
      7 1221 1373 (N17W05 L=119)   363     82        2 (1193/1335-1257/1402)  53.1 W02-W06 N15-N19
      8 1270 1345 (N17W08 L=123)   208     47        0 (1241/1323-1290/1364) -18.4 W06-W09 N15-N17
      9 1854 1148 (N20W56 L=171)   590     88        1 (1814/1115-1882/1179) -12.5 W51-W60 N18-N21
     10 1408  919 (S05W25 L=139)   290     54        1 (1387/ 891-1438/ 940)  56.3 W23-W26 S06-S03
     11 1230  852 (S13W16 L=130)   415    166        0 (1158/ 820-1310/ 885)  32.0 W11-W19 S15-S09
     12 1049  803 (S19W06 L=121)   282     54        1 (1031/ 771-1076/ 832) -71.6 W04-W07 S20-S17
     13  761  790 (S25E12 L=103)   328    114        0 ( 706/ 759- 809/ 822)  -8.1 E15-E08 S26-S24
     14  700  652 (S35E14 L=101)   332     76        1 ( 674/ 622- 743/ 683)  -9.5 E16-E10 S37-S33
     15 1712  646 (S12W56 L=170)   534     71        2 (1686/ 610-1739/ 672)  90.0 W53-W58 S14-S10
     16 1457  467 (S29W45 L=159)  1123    301        3 (1328/ 401-1561/ 501)   4.6 W33-W58 S32-S26
     17  807  461 (S46E02 L=112)   443     69        2 ( 763/ 436- 837/ 490)   0.0 E06-W00 S47-S43
     18  792  221 (S68W05 L=120)   858    158        1 ( 715/ 188- 857/ 246)   0.0 E07-W17 S69-S65
  

    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