Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1147 1877 (N59W04 L=302)   571     58        2 (1112/1860-1180/1892)   0.0 W00-W07 N56-N61
      2  785 1856 (N53E34 L=264)   658     61        2 ( 760/1826- 810/1882) -71.6 E37-E29 N49-N56
      3  989 1866 (N56E13 L=285)   422     74        2 ( 952/1849-1026/1880)  23.2 E16-E09 N53-N57
      4 1121 1645 (N36W02 L=300)   278     62        0 (1097/1608-1136/1675) -59.0 W00-W03 N33-N38
      5  376 1628 (N31E60 L=239)   576     57        2 ( 350/1603- 397/1649)  71.6 E64-E56 N29-N33
      6  726 1568 (N28E25 L=274)   713    189        1 ( 665/1513- 809/1591)  24.0 E28-E18 N23-N29
      7  188  981 (S10E61 L=237)   366     71        1 ( 176/ 940- 206/1020)  68.2 E62-E58 S12-S07
      8  902  846 (S16E06 L=292)   233     54        0 ( 878/ 812- 927/ 877)  63.4 E07-E04 S18-S14
      9  619  743 (S24E25 L=273)   586    157        0 ( 559/ 704- 699/ 769)  21.8 E29-E19 S26-S21
     10  995  480 (S39W02 L=300)   680    118        3 ( 953/ 442-1042/ 534) -23.2 E01-W06 S41-S35
     11 1533  375 (S41W52 L=351)   525    128        1 (1487/ 348-1604/ 396)  14.0 W46-W59 S43-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