Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  609 1744 (N30E52 L=355)  1268    262        2 ( 555/1680- 647/1789) -18.4 E59-E44 N25-N33
      2  495 1639 (N20E54 L=353)   758    173        2 ( 431/1612- 546/1663)  10.3 E60-E49 N18-N21
      3  369 1543 (N12E60 L=347)  1095    207        1 ( 337/1478- 415/1617) -77.5 E63-E56 N07-N16
      4  444 1490 (N10E51 L=357)  1873    427        6 ( 296/1392- 568/1616)  53.6 E61-E44 N01-N20
      5 1083 1512 (N26E10 L= 38)   991    268        2 ( 987/1422-1197/1607) -18.4 E18-W00 N23-N30
      6 1480 1495 (N37W18 L= 65)   508     98        1 (1441/1458-1526/1531) -14.0 W14-W22 N35-N39
      7  661 1431 (N11E32 L= 15)   292     77        0 ( 644/1392- 677/1466) -36.9 E34-E29 N09-N12
      8 1634 1350 (N32W32 L= 80)   298     81        1 (1594/1326-1665/1368)   0.0 W28-W35 N31-N33
      9  569 1243 (S02E32 L= 15)   513    139        0 ( 515/1224- 631/1270)   0.0 E36-E27 S03-N00
     10  716  567 (S39E08 L= 40)   405     91        0 ( 671/ 538- 765/ 592)   0.0 E11-E03 S39-S37
  

    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