Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1597 1532 (N21W53 L=133)  1938    346        2 (1454/1492-1719/1589) -35.8 W44-W63 N15-N28
      2  300 1437 (N47E51 L= 29)  1246    319        2 ( 236/1386- 348/1481) -25.0 E59-E44 N43-N50
      3 1284 1283 (N16W22 L=101)   376    124        3 (1227/1245-1317/1336)  23.2 W17-W24 N14-N18
      4  309 1127 (N27E47 L= 32)   629    100        0 ( 274/1093- 347/1165) -39.8 E50-E44 N24-N29
      5  888  613 (S14E17 L= 62)   956    229        1 ( 838/ 517- 937/ 713)  47.7 E22-E12 S19-S08
      6 1687  602 (S37W38 L=117)   429     90        1 (1662/ 577-1720/ 627) -39.8 W34-W41 S39-S33
      7  305  553 (S09E64 L= 16)  2927    311        6 ( 243/ 505- 395/ 625) -18.4 E72-E55 S12-S03
      8 1693  531 (S43W39 L=119)  1160    132        2 (1660/ 495-1727/ 578) -48.4 W34-W43 S46-S37
      9  633  359 (S26E44 L= 35)   848    317        2 ( 595/ 278- 658/ 463)  54.5 E51-E39 S32-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