Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  370 1578 (N11E60 L=206)   573     68        2 ( 344/1548- 399/1599)  45.0 E64-E57 N09-N12
      2 1694 1416 (N35W35 L=302)  1508    275        1 (1576/1344-1798/1486)  -2.6 W24-W46 N32-N36
      3  330 1332 (S04E52 L=215)   497    122        0 ( 295/1286- 378/1376)  60.3 E54-E49 S07-N00
      4  335 1238 (S10E49 L=218)   376     66        1 ( 318/1208- 356/1275) -36.9 E50-E46 S11-S07
      5  959 1181 (S00E07 L=260)   199     55        0 ( 945/1152- 972/1217) -45.0 E08-E05 S01-N01
      6 1869 1054 (N18W57 L=324)   435     37        1 (1848/1030-1894/1078)  56.3 W54-W59 N16-N19
      7 1221  791 (S15W17 L=284)   236    108        0 (1194/ 733-1255/ 824)  80.5 W16-W18 S18-S12
      8 1225  653 (S22W22 L=289)   742    174        1 (1166/ 601-1290/ 728) -15.3 W16-W27 S24-S19
      9 1181  597 (S27W21 L=288)   211     59        0 (1158/ 569-1199/ 627) -26.6 W18-W23 S28-S25
     10  730  488 (S46E05 L=262)  1833    422        3 ( 566/ 420- 925/ 565)   3.6 E20-W12 S48-S43
  

    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