Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  895 1778 (N60W11 L=287)   394     55        1 ( 872/1744- 914/1803)  50.2 W08-W14 N56-N62
      2 1399 1614 (N34W42 L=319)   992    111        1 (1371/1559-1446/1672)  36.9 W38-W47 N30-N38
      3 1505 1582 (N30W49 L=325)   788    124        1 (1454/1560-1573/1607) -40.6 W44-W53 N26-N32
      4  586 1457 (N41E25 L=252)  1105    271        0 ( 456/1434- 706/1480) -12.3 E37-E14 N37-N43
      5  238 1455 (N45E71 L=205)   882     83        2 ( 205/1430- 285/1471)   2.9 E82-E61 N43-N46
      6 1735 1316 (N09W55 L=331)   400     84        0 (1708/1276-1755/1359)  23.2 W51-W58 N06-N10
      7 1727  662 (S31W45 L=321)  2871    385        4 (1604/ 559-1844/ 829)  23.5 W32-W57 S36-S22
      8  211  629 (S07E72 L=204)  1878    115        1 ( 181/ 569- 239/ 681)  38.7 E80-E67 S11-S03
      9 1431  635 (S25W19 L=295)   448    117        0 (1365/ 622-1486/ 651) -12.5 W14-W23 S26-S23
     10  625  447 (S23E39 L=237)   469     70        0 ( 592/ 417- 665/ 477) -53.1 E40-E36 S25-S20
     11 1033  267 (S44E18 L=259)  3466    469        4 ( 842/ 185-1253/ 378)   9.5 E40-W02 S48-S38
  

    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