Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  837 1597 (N35E11 L=324)   284     69        0 ( 818/1568- 851/1633) -78.7 E11-E09 N32-N37
      2 1687 1464 (N22W50 L= 25)   944    292        4 (1593/1360-1809/1558) -40.9 W43-W59 N15-N28
      3 1643 1299 (N12W42 L= 17)   249     34        0 (1625/1276-1660/1323) -45.0 W41-W43 N10-N13
      4 1262 1151 (N04W15 L=350)  1674    464        6 (1098/ 964-1336/1273) -56.3 W05-W18 S07-N12
      5  835 1111 (N04E11 L=324)   602    229        2 ( 796/1053- 909/1163) -56.3 E13-E06 N00-N06
      6  341  915 (S05E46 L=289)   386     49        1 ( 320/ 891- 364/ 941) -71.6 E47-E43 S06-S03
      7 1781  814 (S17W52 L= 27)   606    124        1 (1746/ 752-1818/ 863)  54.5 W49-W55 S21-S14
      8  261  746 (S15E56 L=279)  1768    314        3 ( 222/ 644- 326/ 820)  32.7 E62-E48 S20-S10
      9  949  622 (S27E07 L=328)   301     88        0 ( 922/ 583- 969/ 668)  59.0 E09-E05 S28-S23
     10 1503  434 (S42W36 L= 11)  1979    414        7 (1327/ 378-1674/ 522)  14.0 W21-W49 S46-S36
  

    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