Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  711 1790 (N39E51 L=138)  1409    174        2 ( 651/1726- 779/1845) -12.5 E66-E39 N35-N43
      2  760 1629 (N29E34 L=154)   894    215        2 ( 647/1591- 834/1670)  45.0 E40-E30 N23-N33
      3 1593 1513 (N43W32 L=221)  1477    211        1 (1520/1422-1637/1590) -45.0 W24-W36 N36-N47
      4 1533 1321 (N27W27 L=215)  2339    732        2 (1350/1029-1665/1534) -38.9 W10-W41 N12-N37
      5  901 1508 (N24E19 L=169)   236     87        0 ( 868/1475- 937/1529) -21.8 E21-E16 N22-N25
      6  434 1481 (N12E51 L=137)   709    118        1 ( 400/1427- 471/1528)  90.0 E53-E49 N08-N16
      7 1138 1407 (N23E02 L=187)   483    167        1 (1072/1360-1194/1469) -24.0 E07-W02 N20-N25
      8 1836 1173 (N26W58 L=246)  1194    116        1 (1810/1111-1857/1231) -81.9 W54-W60 N22-N29
      9  316 1060 (S15E49 L=140)  1194    226        4 ( 254/ 983- 372/1131)  48.0 E54-E44 S20-S09
     10  385  937 (S21E41 L=148)   374    100        0 ( 366/ 898- 416/ 983) -51.3 E42-E37 S22-S17
     11  542  895 (S20E28 L=161)   284    125        0 ( 496/ 857- 572/ 944)  60.3 E30-E25 S23-S16
     12 1622  925 (N05W40 L=228)   246     47        0 (1594/ 907-1649/ 943)  45.0 W37-W41 N03-N06
     13 1300  830 (S07W21 L=209)   204     74        0 (1261/ 813-1328/ 854)  45.0 W18-W22 S08-S05
     14 1514  785 (S05W36 L=224)   294    107        1 (1479/ 758-1571/ 806)  33.7 W33-W39 S07-S02
     15 1471  714 (S10W35 L=223)   272     43        1 (1447/ 688-1493/ 737) -14.0 W32-W36 S11-S08
     16 1392  622 (S17W33 L=221)   691    252        2 (1285/ 602-1487/ 668)  33.7 W26-W39 S20-S12
     17  376  608 (S43E43 L=146)  1151    148        1 ( 299/ 589- 445/ 631)  11.3 E54-E34 S45-S40
  

    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