Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  485 1562 (N53E37 L= 73)   606     83        0 ( 462/1520- 509/1594) -45.0 E41-E33 N49-N55
      2  721 1563 (N47E10 L= 99)   767    133        1 ( 658/1534- 774/1586) -31.0 E16-E05 N44-N50
      3 1792 1451 (N12W71 L=180)   765     74        0 (1776/1411-1806/1487)  24.0 W67-W77 N09-N13
      4  680 1278 (N29E17 L= 92)   241     95        0 ( 658/1237- 715/1307)  36.9 E19-E14 N26-N30
      5  273 1131 (N26E54 L= 56)  1914    368        0 ( 180/1051- 365/1198) -26.6 E65-E44 N19-N30
      6  139 1056 (N22E73 L= 37)   775    118        0 ( 119/1000- 162/1113)  45.0 E77-E68 N18-N25
      7  799  818 (S01E17 L= 92)   635    195        2 ( 744/ 731- 851/ 889) -77.5 E19-E15 S06-N03
      8 1072  745 (S10E03 L=107)   252     57        1 (1045/ 714-1101/ 773)  14.0 E04-W00 S11-S08
      9 1716  720 (S28W42 L=152)   494     44        1 (1694/ 690-1736/ 755)  63.4 W40-W44 S29-S25
     10  464  385 (S24E59 L= 50)  1822    218        2 ( 380/ 312- 542/ 439) -55.0 E65-E53 S30-S18
  

    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