Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  791 1864 (N44E56 L=321)   863     75        1 ( 759/1839- 819/1893) -10.6 E65-E49 N41-N46
      2  892 1595 (N27E25 L=352)   517    155        0 ( 864/1538- 938/1657) -73.3 E26-E22 N22-N32
      3  398 1605 (N16E64 L=313)  1251    185        3 ( 341/1528- 436/1649) -13.5 E80-E55 N11-N17
      4 1278 1545 (N35W02 L= 18)   238     71        1 (1254/1506-1303/1573) -36.9 W00-W04 N32-N36
      5  632  943 (S18E21 L=355)   267     90        1 ( 602/ 920- 674/ 979)  -9.5 E24-E18 S19-S16
      6  561  929 (S20E26 L=351)   225     78        1 ( 544/ 906- 590/ 966) -56.3 E27-E23 S21-S18
      7  821  731 (S26E05 L= 12)   538    173        1 ( 740/ 706- 906/ 754)  33.7 E09-W00 S28-S22
      8 1289  353 (S36W43 L= 59)  3429    538       12 (1114/ 279-1458/ 466)  60.3 W30-W50 S45-S23
  

    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