Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  955 1821 (N62E12 L= 68)   405     61        3 ( 936/1790- 979/1842)  71.6 E15-E09 N58-N64
      2  810 1697 (N49E22 L= 57)   360     44        1 ( 788/1676- 846/1718) -15.9 E25-E18 N46-N50
      3  368 1569 (N36E63 L= 16)   710    105        1 ( 357/1513- 383/1618) -28.3 E70-E57 N32-N39
      4 1525 1564 (N39W42 L=122)   528     95        1 (1486/1526-1562/1596) -51.3 W39-W44 N36-N41
      5 1750 1477 (N32W65 L=145)   701     90        2 (1710/1444-1782/1500) -26.6 W59-W69 N29-N33
      6  282 1431 (N26E64 L= 15)   516     63        1 ( 260/1401- 304/1460) -16.7 E70-E60 N24-N27
      7 1369 1366 (N25W23 L=103)   232     94        0 (1330/1336-1401/1394)  26.6 W20-W26 N23-N27
      8 1607 1274 (N19W41 L=121)   570    161        0 (1548/1218-1662/1324) -41.2 W36-W44 N15-N22
      9 1746  856 (S07W52 L=132)   601     62        1 (1718/ 822-1773/ 883)  90.0 W49-W54 S09-S05
     10 1240  835 (S08W14 L= 94)   278    123        1 (1194/ 812-1282/ 869)   9.5 W10-W16 S09-S05
     11 1092  734 (S15W05 L= 85)   685    223        0 (1064/ 644-1120/ 841) -81.3 W02-W06 S20-S07
     12  792  719 (S16E14 L= 65)   818    299        4 ( 698/ 649- 860/ 800) -42.0 E20-E10 S20-S11
     13  263  710 (S20E59 L= 20)  1021    175        0 ( 214/ 656- 335/ 744)  17.5 E69-E50 S23-S16
     14 1379  674 (S18W24 L=104)   228     68        0 (1351/ 649-1400/ 699) -45.0 W22-W25 S20-S16
     15 1606  576 (S25W45 L=125)  1360    172        5 (1566/ 500-1646/ 654)  90.0 W41-W47 S30-S19
     16  880  578 (S26E09 L= 71)   911    392        2 ( 689/ 534-1036/ 613)  14.0 E23-W01 S29-S22
     17  437  538 (S30E46 L= 34)   436    136        1 ( 386/ 484- 475/ 575)  28.3 E54-E40 S35-S27
  

    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