Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  768 1547 (N17E31 L=346)   373     85        1 ( 711/1516- 807/1573)  51.3 E32-E27 N13-N18
      2 1692 1451 (N41W34 L= 51)   710    129        3 (1652/1424-1738/1481)  12.5 W29-W38 N39-N43
      3  763 1417 (N09E26 L=351)   511     84        2 ( 725/1378- 787/1455)  90.0 E28-E23 N06-N11
      4  329 1303 (S07E52 L=324)   339     69        0 ( 314/1268- 355/1329) -76.0 E53-E50 S09-S05
      5 1031 1166 (N02E03 L= 13)   182     76        1 ( 992/1143-1074/1188)   9.5 E06-W00 N01-N02
      6  232  994 (S27E53 L=324)   520    148        0 ( 195/ 954- 279/1041) -16.7 E57-E47 S28-S24
      7 1769  920 (N11W50 L= 66)   319     66        1 (1753/ 896-1789/ 949)  90.0 W48-W52 N09-N13
      8  936  885 (S16E01 L= 16)   553    105        2 ( 892/ 847- 987/ 924)  -9.5 E04-W02 S17-S13
      9 1478  318 (S30W57 L= 74)   628     98        1 (1446/ 290-1507/ 345)   8.1 W53-W61 S30-S28
  

    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