Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  443 1601 (N19E58 L=111)  3321    295       10 ( 339/1528- 538/1684)  49.8 E68-E50 N13-N26
      2 1642 1649 (N58W44 L=213)   678     57        1 (1623/1611-1664/1675)  32.7 W38-W54 N53-N62
      3  381 1512 (N12E58 L=111)   501     96        0 ( 353/1471- 407/1548)  78.7 E60-E56 N09-N14
      4  575 1491 (N15E41 L=128)   450     64        1 ( 543/1453- 607/1513)  90.0 E43-E38 N12-N16
      5 1708 1311 (N32W42 L=211)   778    168        1 (1672/1239-1756/1375) -41.2 W38-W46 N28-N35
      6  173 1237 (S08E70 L= 99)   767     75        0 ( 135/1211- 206/1258)   0.0 E78-E63 S09-S06
      7  354 1094 (S13E45 L=124)   282     71        0 ( 335/1054- 375/1121)  63.4 E46-E43 S15-S11
      8  148  933 (S27E70 L= 99)   747     83        1 ( 128/ 897- 173/ 974)  45.0 E74-E64 S29-S24
      9  862  725 (S24E03 L=166)   319    104        1 ( 838/ 679- 888/ 766) -78.7 E03-E01 S26-S21
     10  776  671 (S29E07 L=162)   208     76        0 ( 737/ 650- 812/ 694)  -9.5 E10-E04 S29-S28
     11 1302  623 (S19W27 L=196)   245     54        0 (1286/ 594-1316/ 654) -56.3 W26-W28 S20-S17
     12 1432  628 (S16W36 L=205)   263     57        0 (1402/ 604-1460/ 651)  63.4 W34-W36 S17-S13
     13  317  587 (S47E51 L=118)  1240    144        1 ( 264/ 562- 384/ 630)   0.0 E60-E41 S49-S44
     14  754  539 (S39E06 L=163)   302     53        0 ( 733/ 511- 777/ 562)  68.2 E07-E04 S41-S36
     15 1548  371 (S28W61 L=230)  1245    154        0 (1497/ 330-1606/ 426)  83.7 W57-W64 S31-S22
     16 1152  397 (S37W28 L=197)   413     75        1 (1110/ 377-1188/ 415)   9.5 W24-W30 S38-S35
     17  773  377 (S50W00 L=169)   657    107        1 ( 751/ 333- 812/ 412)  76.0 E02-W03 S54-S46
     18  792  280 (S59W07 L=176)   859    124        1 ( 746/ 243- 849/ 315)  41.2 W01-W10 S63-S53
  

    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