Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  932 1874 (N72W06 L=152)   694     68        1 ( 911/1844- 956/1899)  45.0 W01-W11 N67-N76
      2 1594 1565 (N32W54 L=201)   471     81        0 (1562/1545-1636/1581) -14.0 W50-W58 N30-N33
      3 1130 1530 (N37W14 L=160)   371     95        0 (1083/1505-1177/1553)  14.0 W09-W17 N35-N37
      4 1210 1521 (N35W19 L=166)   277     69        1 (1190/1484-1229/1550)  63.4 W17-W21 N32-N36
      5 1558 1492 (N28W46 L=193)   607    134        0 (1527/1427-1586/1540)  76.0 W44-W48 N23-N31
      6 1711 1476 (N25W62 L=208)   659     77        1 (1669/1441-1737/1508)  13.0 W54-W67 N23-N26
      7  815 1350 (N27E11 L=135)   536    173        0 ( 784/1285- 865/1417) -66.0 E13-E08 N22-N32
      8  402 1280 (N25E44 L=103)  2916    460        4 ( 289/1155- 576/1347)  11.3 E54-E28 N17-N29
      9 1292 1249 (N17W19 L=166)   216     69        1 (1249/1231-1324/1268) -11.3 W16-W21 N15-N18
     10 1347 1220 (N14W23 L=169)   204     75        0 (1303/1207-1374/1240)   0.0 W19-W24 N13-N15
     11  413 1071 (N12E41 L=105)   718    131        2 ( 359/1031- 461/1110) -35.5 E45-E37 N09-N14
     12 1751  857 (S13W50 L=196)   567     71        0 (1725/ 826-1774/ 882) -53.1 W47-W52 S15-S11
     13  689  797 (S06E23 L=123)   280     97        0 ( 643/ 770- 729/ 814) -18.4 E26-E20 S07-S04
     14  255  745 (S08E60 L= 86)   643     51        1 ( 232/ 715- 285/ 775)  63.4 E63-E56 S09-S05
     15 1061  477 (S31E03 L=143)   265    103        1 (1012/ 466-1114/ 489) -18.4 E06-W00 S31-S29
     16  448  405 (S32E57 L= 90)   672    113        0 ( 425/ 354- 481/ 462)  66.0 E59-E53 S36-S27
     17  868  403 (S35E18 L=128)   449    130        1 ( 808/ 388- 937/ 417)  -5.7 E22-E12 S35-S32
     18 1253  250 (S54W12 L=159)   433     35        0 (1232/ 235-1279/ 272) -31.0 W10-W15 S56-S51
  

    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