Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  957 1777 (N38E32 L= 99)  2726    464        7 ( 829/1657-1036/1868) -22.9 E47-E21 N30-N43
      2  458 1674 (N18E60 L= 70)  3784    565        9 ( 325/1604- 556/1755)  32.9 E71-E51 N11-N24
      3  535 1542 (N11E46 L= 84)  1503    331        5 ( 452/1465- 604/1602)   4.1 E54-E40 N06-N14
      4  462 1454 (N04E48 L= 83)   289     71        0 ( 443/1421- 479/1483) -31.0 E49-E44 N02-N05
      5 1323 1336 (N20W08 L=138)   749    136        1 (1272/1289-1386/1382) -12.5 W04-W13 N16-N22
      6 1601  573 (S13W51 L=181)   299     56        1 (1582/ 549-1624/ 599) -56.3 W48-W52 S14-S11
      7 1157  299 (S43W38 L=168)   965    192        2 (1083/ 277-1242/ 344)  17.1 W31-W44 S45-S40
  

    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