Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  457 1598 (N38E50 L=237)   494    117        3 ( 432/1573- 485/1625) -18.4 E52-E46 N36-N40
      2  678 1196 (N13E23 L=264)   818    317        1 ( 561/1142- 777/1300)  31.6 E30-E16 N09-N19
      3  164  966 (S04E67 L=220)   565    138        3 ( 143/ 928- 188/1002) -29.7 E70-E63 S06-S01
      4  615  880 (S07E26 L=261)   311    127        1 ( 571/ 851- 671/ 903)   0.0 E28-E22 S08-S05
      5 1848  826 (S09W64 L=351)  1171    223        0 (1785/ 769-1898/ 899)  31.6 W56-W70 S12-S04
      6  715  781 (S13E19 L=268)   258     68        1 ( 680/ 746- 745/ 805) -36.9 E21-E17 S14-S11
      7  948  647 (S21E04 L=283)   300     95        1 ( 929/ 603- 964/ 701)  80.5 E05-E03 S23-S17
      8 1710  588 (S25W54 L=341)   415     58        0 (1686/ 565-1728/ 620)  45.0 W51-W55 S26-S22
      9  875  572 (S26E09 L=278)   546    118        1 ( 804/ 553- 927/ 598)   0.0 E13-E05 S27-S24
     10  412  401 (S42E53 L=234)  1303    302        2 ( 363/ 374- 463/ 439)  15.9 E61-E47 S45-S38
  

    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