Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  398 1616 (N29E55 L=137)  1394    150        5 ( 334/1573- 451/1658)  32.0 E60-E49 N25-N32
      2 1174 1515 (N26W06 L=198)   801    248        0 (1071/1440-1255/1569) -22.6 E01-W11 N22-N29
      3  447 1360 (N12E40 L=152)  2415    444        9 ( 332/1181- 577/1514) -45.0 E51-E29 N01-N21
      4  510 1423 (N16E37 L=155)   239     82        0 ( 493/1390- 526/1468) -80.5 E37-E35 N13-N19
      5  722 1397 (N15E22 L=171)   410     88        1 ( 685/1361- 752/1440)  78.7 E23-E19 N13-N18
      6  645 1405 (N15E27 L=165)   225     34        1 ( 629/1381- 669/1428) -56.3 E28-E25 N13-N16
      7 1695 1240 (N14W43 L=235)   370     41        1 (1660/1219-1725/1259) -11.3 W39-W45 N12-N15
      8 1780 1106 (N07W50 L=242)   412     69        1 (1765/1069-1796/1155) -68.2 W48-W51 N04-N10
      9  410 1107 (S03E40 L=152)   489    103        0 ( 354/1086- 460/1129)   7.1 E44-E36 S04-S02
     10  192 1039 (S08E59 L=133)  1389    199        5 ( 153/ 977- 225/1118) -69.4 E64-E55 S11-S03
     11  306  878 (S17E49 L=144)  2504    524        9 ( 234/ 642- 382/1052)  64.4 E54-E41 S31-S06
     12  411  888 (S16E39 L=153)   534    112        0 ( 374/ 844- 437/ 939)  50.2 E42-E37 S19-S13
     13  638  882 (S16E23 L=169)  1511    521        3 ( 477/ 793- 859/ 933)  17.7 E34-E09 S21-S11
     14 1475  478 (S34W40 L=232)   497    127        0 (1424/ 448-1528/ 517) -14.0 W33-W45 S35-S31
     15  985  459 (S41W02 L=195)   360     42        0 ( 957/ 432-1009/ 487) -36.9 W00-W04 S42-S38
     16 1445  212 (S54W64 L=256)   867    102        0 (1398/ 189-1497/ 239)  -1.8 W52-W83 S54-S52
  

    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