Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1535 1714 (N28W65 L=259)   855    103        1 (1499/1680-1567/1760)  56.3 W60-W68 N25-N31
      2 1216 1693 (N38W38 L=233)   277     33        1 (1191/1676-1241/1712) -56.3 W36-W39 N36-N39
      3 1399 1614 (N27W45 L=239)   825    188        0 (1307/1597-1491/1639) -36.0 W39-W51 N23-N31
      4 1012 1611 (N40W19 L=214)   346     87        1 ( 966/1588-1051/1627) -15.9 W15-W22 N38-N41
      5 1451 1471 (N18W40 L=234)   758    366        0 (1271/1451-1601/1495) -26.6 W28-W50 N12-N23
      6 1116 1361 (N22W15 L=210)   723    210        2 (1030/1333-1232/1398) -13.0 W09-W22 N20-N24
      7  656 1185 (N24E18 L=177)   227     89        1 ( 620/1152- 690/1202) -59.0 E20-E15 N21-N26
      8  376 1137 (N28E40 L=154)   254     76        0 ( 354/1109- 399/1173) -53.1 E41-E38 N26-N30
      9  485 1018 (N19E33 L=161)   672    219        2 ( 408/ 985- 605/1047) -19.7 E38-E24 N16-N22
     10  371  962 (N18E43 L=151)   487     55        1 ( 335/ 932- 400/ 993) -53.1 E46-E40 N15-N20
     11 1787  933 (S22W47 L=242)   355    102        0 (1748/ 918-1823/ 951) -14.0 W43-W51 S22-S19
     12 1518  825 (S19W24 L=218)   216     84        0 (1494/ 786-1538/ 857)  36.9 W21-W25 S20-S17
     13  375  759 (N06E48 L=146)   668    105        1 ( 324/ 735- 426/ 796) -20.6 E52-E44 N04-N08
     14 1162  657 (S18E02 L=192)   257    100        1 (1148/ 611-1185/ 700)  59.0 E04-W00 S21-S15
     15 1669  496 (S48W34 L=229)  1389    325        1 (1549/ 410-1785/ 586)   5.9 W20-W49 S50-S45
     16  495  324 (S23E63 L=131)  1094    137        2 ( 458/ 264- 529/ 390)  51.3 E68-E58 S27-S17
     17 1358  304 (S51E01 L=194)   759    118        1 (1311/ 279-1411/ 341)   6.3 E05-W04 S53-S48
  

    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