Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1002 1732 (N35E23 L= 26)   275     76        1 ( 983/1703-1024/1760)  90.0 E24-E21 N33-N38
      2  272 1215 (S12E55 L=354)   564     62        0 ( 247/1178- 298/1244) -76.0 E57-E52 S13-S09
      3  304 1164 (S14E50 L=359)   894    183        1 ( 245/1095- 357/1242)  59.0 E54-E47 S18-S08
      4  650 1035 (S15E22 L= 27)   476    182        0 ( 597/ 954- 695/1121) -35.0 E27-E17 S18-S11
      5 1717 1061 (N13W42 L= 91)   285     71        2 (1692/1028-1743/1089) -76.0 W39-W43 N11-N15
      6 1331  976 (S02W18 L= 67)   967    228        1 (1229/ 938-1435/1052)   0.0 W10-W24 S03-N00
      7 1850  999 (N15W55 L=104)   369    108        1 (1824/ 980-1875/1025)   9.5 W52-W58 N13-N16
      8 1894  936 (N13W62 L=111)   513     84        2 (1873/ 904-1912/ 965) -36.9 W60-W65 N11-N14
      9 1761  839 (N03W52 L=101)   965    162        1 (1731/ 759-1792/ 911) -54.5 W48-W56 N00-N07
     10  644  849 (S26E19 L= 30)   233     44        1 ( 619/ 826- 670/ 877) -33.7 E20-E16 S27-S24
     11  388  839 (S32E38 L= 11)  1059    331        1 ( 260/ 811- 500/ 865)   8.1 E49-E28 S35-S29
     12 1160  724 (S20W16 L= 65)   282     91        0 (1109/ 706-1210/ 739)  15.9 W12-W19 S22-S19
     13 1400  676 (S16W32 L= 81)   674    239        0 (1310/ 651-1510/ 719)  26.6 W25-W37 S18-S11
     14 1683  596 (S12W56 L=104)   955    199        2 (1616/ 556-1725/ 650)  63.4 W50-W59 S16-S08
     15 1122  561 (S31W19 L= 68)   517    132        0 (1052/ 540-1192/ 583)  40.6 W15-W22 S33-S27
     16 1701  517 (S15W62 L=111)  1062    168        2 (1645/ 477-1747/ 560)  58.0 W58-W66 S19-S11
     17 1632  524 (S17W55 L=104)   692    137        0 (1586/ 498-1683/ 559)  45.0 W51-W57 S19-S13
  

    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