Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1240 1637 (N35E03 L=306)   752    140        3 (1170/1598-1310/1674)  12.5 E07-W02 N32-N36
      2  873 1466 (N14E19 L=289)   248     58        0 ( 835/1448- 905/1480)   0.0 E22-E17 N13-N15
      3 1517 1331 (N23W22 L=331)   219     57        0 (1506/1299-1532/1365) -45.0 W20-W23 N20-N24
      4 1785 1345 (N33W45 L=354)   414     38        0 (1754/1330-1810/1363)  23.2 W41-W48 N31-N34
      5 1382 1322 (N18W13 L=322)   481    146        0 (1313/1288-1442/1360)   0.0 W08-W17 N17-N19
      6 1778 1219 (N24W45 L=354)   458     66        1 (1761/1179-1798/1261) -53.1 W42-W47 N22-N27
      7  116 1081 (S21E72 L=236)   770     96        4 (  98/1043- 133/1129)  68.2 E76-E69 S22-S17
      8  788 1093 (S09E15 L=294)   211     69        2 ( 767/1067- 817/1127) -26.6 E16-E12 S10-S07
      9  633 1055 (S14E23 L=285)   442    177        0 ( 612/ 972- 667/1116) -82.9 E24-E20 S18-S10
     10 1788 1044 (N14W48 L=357)   944    143        3 (1747/ 976-1825/1098) -50.2 W44-W53 N11-N17
     11 1611  996 (N06W35 L=344)   613    182        2 (1523/ 957-1681/1041)  45.0 W29-W40 N01-N10
     12  662  824 (S27E17 L=292)   375     84        0 ( 620/ 796- 701/ 849)   0.0 E20-E14 S27-S25
     13 1876  694 (S01W72 L= 20)   754     62        1 (1849/ 679-1916/ 712)   7.6 W66-W81 S02-N01
     14  386  620 (S46E40 L=269)  1805    385        3 ( 234/ 554- 563/ 667)   0.0 E61-E21 S47-S42
     15 1787  532 (S12W72 L= 20)  1259    141        6 (1751/ 467-1822/ 590) -45.0 W67-W78 S15-S08
     16 1542  423 (S26W55 L=  4)   930    182        2 (1442/ 388-1606/ 485)   2.9 W43-W63 S27-S23
     17 1463  255 (S37W66 L= 15)   971     95        1 (1427/ 233-1505/ 287)  68.2 W60-W72 S39-S34
  

    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