Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  425 1332 (N14E45 L=244)   412     77        1 ( 402/1296- 451/1372) -33.7 E48-E41 N12-N16
      2  715 1258 (N13E22 L=267)   357    151        1 ( 673/1224- 769/1280) -18.4 E24-E18 N11-N14
      3  832 1031 (N00E12 L=277)   302     47        1 ( 809/1010- 858/1052)   0.0 E13-E10 S01-N01
      4  862  863 (S10E08 L=280)   353     83        1 ( 835/ 837- 897/ 887) -14.0 E10-E06 S11-S08
      5  896  749 (S17E05 L=283)   549    191        1 ( 814/ 729- 963/ 797)   0.0 E10-E01 S17-S14
      6 1559  716 (S13W39 L=328)   860    205        0 (1526/ 640-1619/ 795)  73.3 W37-W42 S17-S07
      7 1322  692 (S16W23 L=312)   468    172        1 (1262/ 656-1372/ 762)  54.5 W19-W25 S19-S11
      8  220  738 (S25E63 L=225)   591     69        0 ( 181/ 726- 249/ 760)  13.0 E71-E58 S26-S23
      9 1650  609 (S19W51 L=340)   985    283        1 (1546/ 514-1729/ 681)  50.2 W43-W57 S26-S13
     10  737  543 (S32E16 L=273)   282    107        1 ( 678/ 528- 780/ 571)  -7.1 E20-E12 S33-S31
     11 1496  361 (S38W51 L=339)  1030    115        2 (1454/ 326-1541/ 412)  74.1 W45-W54 S41-S33
  

    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