Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  880 1556 (N28E13 L=224)   227     60        0 ( 864/1523- 894/1587)  90.0 E13-E11 N26-N30
      2  755 1483 (N23E20 L=216)   550    148        1 ( 676/1452- 818/1503) -11.3 E25-E15 N21-N24
      3  668 1406 (N18E25 L=212)   340    110        0 ( 621/1372- 714/1430) -23.2 E28-E21 N15-N19
      4  558 1360 (N14E32 L=205)   289     78        0 ( 548/1317- 572/1404) -80.5 E32-E30 N11-N17
      5 1865 1231 (N15W61 L=297)  2797    384        2 (1769/1128-1954/1349) -24.4 W51-W74 N08-N21
      6 1164 1256 (N10W07 L=244)   425     90        0 (1125/1222-1214/1305) -33.7 W04-W10 N08-N13
      7 1556 1210 (N10W33 L=269)   463    133        0 (1492/1192-1619/1229)  12.5 W28-W37 N08-N11
      8  418 1042 (S05E39 L=197)   691    125        1 ( 397/ 970- 446/1108) -82.9 E41-E37 S09-S01
      9  810 1059 (S03E13 L=223)   238     39        0 ( 779/1042- 837/1077)   0.0 E14-E11 S04-S02
     10 1839  982 (S01W58 L=295)   452    106        2 (1816/ 951-1871/1009)  26.6 W55-W62 S02-N01
     11  139  899 (S13E68 L=169)  1323    187        4 ( 116/ 839- 179/ 963) -74.1 E71-E62 S16-S09
     12  794  889 (S14E14 L=223)   228     48        1 ( 772/ 866- 823/ 909)  33.7 E14-E11 S14-S12
     13 1098  757 (S20W06 L=243)  1075    254        0 ( 998/ 657-1223/ 837)  42.7 W00-W13 S26-S14
     14 1770  546 (S28W62 L=298)  1548    288        3 (1682/ 513-1830/ 600)  30.3 W52-W67 S30-S23
     15  923  553 (S34E05 L=232)   212     35        1 ( 905/ 530- 946/ 576) -33.7 E05-E02 S35-S32
     16  840  496 (S39E11 L=226)   258     61        0 ( 807/ 476- 875/ 520)  31.0 E13-E08 S40-S36
  

    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