Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1086 1533 (N27E11 L= 33)   515     81        2 (1057/1498-1123/1575) -33.7 E13-E07 N25-N29
      2  662 1392 (N08E31 L= 13)   520    167        0 ( 595/1352- 738/1429)  48.8 E34-E27 N03-N11
      3 1143 1397 (N20E03 L= 41)   266     61        1 (1119/1370-1171/1422)  71.6 E04-E01 N18-N21
      4  988 1298 (N10E09 L= 35)   551    121        1 ( 944/1242-1038/1347)  82.9 E10-E07 N06-N14
      5 1831  988 (N16W55 L= 99)  2704    359        5 (1765/ 837-1863/1111) -42.9 W45-W61 N07-N22
      6  954 1093 (S02E06 L= 38)   209     55        0 ( 924/1076- 993/1104)  26.6 E07-E03 S03-N00
      7 1608 1009 (N11W35 L= 79)   488    192        0 (1549/ 941-1653/1057) -18.4 W31-W40 N07-N12
      8  705  984 (S14E17 L= 26)   293    135        1 ( 637/ 934- 760/1014)  49.4 E20-E14 S18-S11
      9  310  942 (S26E46 L=358)   704    155        2 ( 258/ 922- 365/ 977)  33.7 E50-E40 S28-S22
     10  716  810 (S24E13 L= 31)   307     73        0 ( 679/ 788- 758/ 835)   0.0 E15-E09 S24-S23
     11 1298  731 (S13W24 L= 67)   584    162        1 (1261/ 646-1351/ 798) -84.3 W22-W25 S18-S08
     12  401  706 (S39E35 L=  8)  1427    235        2 ( 318/ 618- 482/ 792) -23.6 E43-E27 S43-S35
     13  612  639 (S38E17 L= 27)   644    146        2 ( 547/ 617- 681/ 666)  10.3 E22-E11 S39-S34
     14 1090  638 (S24W15 L= 58)   339    143        0 (1018/ 623-1163/ 648)  18.4 W10-W19 S26-S22
     15 1492  390 (S26W53 L= 96)  1637    334        3 (1343/ 323-1633/ 453)  48.8 W44-W61 S35-S18
     16 1594  356 (S25W64 L=107)   561     75        0 (1558/ 323-1626/ 379)  68.2 W62-W66 S27-S22
  

    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