Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1181 1508 (N33W21 L= 60)   643    156        3 (1122/1441-1238/1546)  20.0 W15-W26 N29-N34
      2 1363 1057 (N02W21 L= 60)   258     62        1 (1332/1030-1388/1093)  14.0 W18-W22 N00-N04
      3 1872  953 (S16W60 L= 99)  2479    423       10 (1825/ 855-1930/1061)  69.0 W53-W68 S22-S09
      4  190  946 (N13E63 L=336)  1079    211        2 ( 153/ 847- 216/1020) -78.7 E66-E60 N07-N17
      5  259  752 (N01E60 L=339)  1122    268        6 ( 198/ 723- 330/ 782)  -7.6 E67-E52 S01-N02
      6 1712  626 (S34W43 L= 81)   489    105        2 (1685/ 598-1749/ 657)  21.8 W39-W46 S35-S31
      7  317  555 (S12E62 L=337)  1114    227        3 ( 290/ 494- 362/ 595)  41.2 E67-E56 S16-S09
      8 1590  526 (S37W32 L= 70)   403     76        1 (1550/ 503-1623/ 553)   9.5 W28-W34 S38-S35
      9 1413  388 (S44W16 L= 55)   401     83        1 (1381/ 373-1449/ 405) -26.6 W13-W19 S45-S41
  

    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