Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  750 1651 (N54E01 L=356)   466     64        2 ( 716/1622- 780/1680)  -9.5 E04-W02 N50-N55
      2 1579 1594 (N23W58 L= 55)  1274    164        5 (1504/1554-1648/1629) -32.5 W50-W64 N18-N25
      3 1519 1540 (N22W49 L= 46)   313     43        1 (1490/1519-1544/1561) -45.0 W47-W50 N20-N23
      4 1623 1515 (N17W56 L= 53)   387     52        1 (1600/1494-1654/1538) -53.1 W52-W58 N15-N19
      5 1428 1446 (N20W37 L= 34)  1027    325        2 (1290/1358-1551/1503) -49.4 W30-W42 N12-N26
      6 1584 1458 (N16W49 L= 46)   282     54        0 (1571/1427-1597/1486)  21.8 W46-W51 N14-N17
      7 1280 1312 (N17W23 L= 19)   693    204        3 (1223/1243-1334/1357)  12.5 W17-W26 N14-N19
      8 1553 1127 (S01W34 L= 31)  1036    183        3 (1504/1049-1596/1190)  32.0 W29-W38 S04-N02
      9 1695  889 (S19W39 L= 36)   337     56        1 (1657/ 868-1723/ 903) -31.0 W36-W41 S20-S16
     10 1521  677 (S27W22 L= 19)   671    182        1 (1458/ 636-1586/ 744) -51.3 W18-W27 S30-S20
     11 1620  699 (S28W30 L= 27)   280     66        0 (1588/ 676-1641/ 724)  11.3 W27-W32 S29-S27
     12 1801  682 (S37W50 L= 47)   986     98        1 (1757/ 653-1848/ 718)  -4.4 W44-W57 S38-S34
     13 1722  588 (S41W41 L= 38)   509     56        0 (1688/ 570-1755/ 602) -14.0 W36-W45 S42-S38
     14 1559  484 (S42W24 L= 20)   475     79        2 (1524/ 468-1600/ 515) -23.2 W20-W27 S44-S38
     15  967  336 (S34E24 L=332)   644    120        0 ( 927/ 312-1023/ 379) -45.0 E26-E21 S36-S30
     16  571  307 (S27E57 L=300)  1626    196        3 ( 485/ 264- 638/ 363) -35.0 E64-E51 S31-S21
  

    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