Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  397 1549 (N30E58 L=144)  2069    200        4 ( 337/1486- 474/1621) -15.4 E75-E46 N26-N34
      2  247 1346 (N15E64 L=137)   507     70        0 ( 224/1316- 266/1377)  76.0 E66-E61 N13-N17
      3  256  899 (S12E56 L=146)   559    116        0 ( 228/ 852- 286/ 955)  45.0 E59-E51 S15-S08
      4  190  862 (S15E65 L=137)   941    150        3 ( 167/ 797- 218/ 946) -78.7 E67-E61 S19-S09
      5  731  776 (S16E17 L=185)   246     44        3 ( 708/ 756- 756/ 796)  33.7 E18-E15 S16-S14
      6  437  755 (S20E39 L=162)   280     85        2 ( 396/ 733- 459/ 782)  26.6 E42-E36 S21-S17
      7 1155  744 (S15W10 L=212)   594    267        1 (1033/ 706-1259/ 775)  -3.8 W02-W17 S16-S12
      8 1647  624 (S19W49 L=251)  1812    318        8 (1565/ 548-1723/ 683)  42.0 W42-W57 S24-S15
      9  753  383 (S44E17 L=185)   284     56        0 ( 729/ 363- 779/ 410)  38.7 E19-E14 S45-S41
     10  517  352 (S50E46 L=156)   534     68        1 ( 485/ 330- 550/ 378) -26.6 E49-E41 S51-S47
  

    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