Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1022 1886 (N53E40 L=155)  1607    153        2 ( 959/1844-1078/1921)   9.5 E50-E30 N48-N56
      2 1184 1760 (N47E15 L=180)   397     65        0 (1152/1742-1222/1776)  38.7 E17-E11 N44-N49
      3  371 1509 (N08E59 L=136)  2163    272        3 ( 311/1385- 422/1620) -39.5 E71-E50 N01-N15
      4 1600 1478 (N40W28 L=222)   744    137        2 (1534/1457-1662/1497)  17.1 W21-W34 N37-N43
      5  177 1100 (S20E65 L=130)  1976    176        2 ( 136/1036- 213/1176)  45.0 E71-E60 S23-S14
      6  331  901 (S28E44 L=150)  1498    318        2 ( 258/ 761- 376/1013)  67.6 E53-E41 S38-S21
      7  785  757 (S26E07 L=187)   445    163        0 ( 721/ 727- 851/ 802)  -5.7 E12-E02 S26-S24
      8 1542  558 (S16W47 L=242)   272     59        1 (1515/ 536-1564/ 581)  76.0 W45-W48 S17-S13
      9 1338  509 (S25W36 L=230)   538     84        0 (1298/ 476-1381/ 542)  71.6 W33-W37 S27-S21
     10  486  498 (S52E28 L=167)   353     43        1 ( 468/ 471- 507/ 532) -78.7 E29-E25 S54-S49
  

    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