Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  842 1864 (N55E47 L=204)   676     67        1 ( 803/1843- 878/1877)   0.0 E53-E40 N53-N56
      2 1189 1727 (N50E02 L=249)   317     44        0 (1169/1701-1214/1746)  56.3 E03-W00 N47-N52
      3  494 1674 (N32E59 L=192)   752     86        3 ( 466/1640- 524/1711) -14.9 E67-E52 N30-N34
      4 1352 1625 (N45W15 L=265)   394     33        1 (1329/1603-1380/1646) -21.8 W12-W17 N43-N46
      5  801 1551 (N29E25 L=225)   289    109        0 ( 756/1531- 854/1569)  26.6 E27-E21 N27-N30
      6 1560 1289 (N25W32 L=283)  1053    453        1 (1383/1169-1692/1426) -21.8 W18-W43 N19-N31
      7  324  926 (S17E47 L=204)   408     91        1 ( 291/ 903- 365/ 943)   0.0 E50-E43 S18-S15
      8  182  816 (S26E66 L=184)  1001    103        2 ( 165/ 763- 206/ 867)  49.4 E71-E61 S30-S23
      9  751  467 (S41E10 L=240)   255     87        0 ( 727/ 450- 785/ 502) -21.8 E12-E07 S42-S38
     10 1077  182 (S60W34 L=284)   558     56        1 (1064/ 150-1096/ 223) -29.7 W26-W40 S63-S55
  

    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