Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1000 1597 (N37W11 L= 53)   303    138        0 ( 934/1579-1063/1613) -18.4 W05-W15 N34-N38
      2  591 1558 (N43E22 L= 20)   268     29        1 ( 564/1539- 613/1574)   0.0 E23-E19 N41-N44
      3 1647  595 (S37W38 L= 80)   419     55        0 (1609/ 575-1673/ 617)   0.0 W33-W40 S38-S35
      4  366  467 (S20E62 L=340)  2072    215        6 ( 311/ 388- 425/ 545)  79.7 E66-E57 S25-S14
      5  464  454 (S22E53 L=349)   797     91        2 ( 422/ 418- 507/ 492) -45.0 E55-E49 S24-S19
      6  631  340 (S33E46 L=356)   939    218        0 ( 573/ 279- 695/ 429)  80.5 E49-E43 S38-S26
      7 1212  381 (S43E01 L= 41)   521     91        1 (1177/ 346-1247/ 420) -81.9 E03-W01 S46-S39
  

    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