Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1195 1693 (N38W34 L=286)  1217    246        2 (1089/1662-1324/1712) -22.4 W25-W43 N34-N41
      2 1516 1682 (N28W57 L=310)   386     71        0 (1475/1666-1548/1699)  -6.3 W52-W61 N26-N29
      3 1688 1640 (N20W72 L=324)  1911    120        1 (1627/1602-1750/1675) -30.3 W65-W80 N16-N23
      4 1449 1484 (N19W39 L=291)  1614    393        0 (1347/1356-1596/1615) -79.2 W35-W44 N07-N28
      5 1562 1540 (N19W50 L=303)  1768    347        1 (1419/1472-1683/1610) -49.1 W43-W57 N11-N26
      6  267 1304 (N38E52 L=201)  2089    305        2 ( 125/1285- 418/1325)  -6.8 E84-E34 N34-N41
      7  598 1299 (N30E20 L=232)   607    237        0 ( 499/1267- 696/1324) -29.7 E27-E13 N25-N33
      8  908 1103 (N11E05 L=248)   417    124        1 ( 862/1051- 941/1149) -74.1 E06-E02 N07-N14
      9 1689 1104 (S08W43 L=295)   311     67        0 (1660/1071-1721/1127)  18.4 W39-W45 S09-S07
     10 1910  893 (S28W65 L=317)  1120     89        0 (1860/ 874-1957/ 918)  -5.7 W57-W77 S29-S25
     11 1731  782 (S29W42 L=294)   866    246        0 (1668/ 677-1782/ 879)  36.9 W35-W47 S33-S24
     12 1609  653 (S33W30 L=282)   596    144        1 (1565/ 617-1656/ 718) -55.0 W26-W33 S37-S27
     13  571  629 (S08E37 L=215)   571    144        0 ( 536/ 556- 596/ 696)  83.7 E39-E35 S12-S03
     14  832  511 (S21E25 L=227)   622    137        1 ( 800/ 454- 856/ 583)  36.9 E29-E21 S24-S16
     15 1295  396 (S42W00 L=253)   897    211        1 (1205/ 327-1359/ 465)  19.7 E08-W06 S44-S38
     16  774  200 (S42E51 L=201)   458     64        1 ( 744/ 178- 803/ 229)   5.7 E56-E46 S42-S39
     17  960  172 (S50E40 L=212)  1542    201        3 ( 869/ 125-1057/ 221)   6.1 E56-E26 S52-S47
  

    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