Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  308 1470 (N14E57 L=352)  4779    409       12 ( 215/1325- 442/1598)  71.6 E65-E47 N04-N23
      2  464 1484 (N16E44 L=  5)  1458    300        8 ( 395/1412- 540/1566)  78.7 E48-E37 N11-N21
      3 1839 1080 (N10W55 L=104)   835    133        1 (1818/1013-1857/1150)  77.5 W53-W57 N05-N14
      4  224 1010 (S14E56 L=353)  1203    200        6 ( 184/ 949- 265/1120)  55.0 E60-E52 S17-S07
      5  149 1079 (S10E65 L=344)   663     79        3 ( 128/1035- 166/1118) -68.2 E68-E62 S12-S07
      6  529  984 (S13E30 L= 19)   314    115        1 ( 507/ 956- 572/1032) -36.9 E31-E27 S15-S10
      7  966  853 (S17E01 L= 48)   625    214        0 ( 918/ 758- 992/ 945) -65.6 E05-W00 S22-S11
      8 1569  854 (S08W36 L= 85)   507    152        0 (1532/ 782-1594/ 929)  90.0 W35-W37 S12-S03
      9  282  854 (S23E51 L=358)  2930    537        4 (  94/ 774- 544/ 901)   8.1 E84-E28 S27-S18
     10 1203  858 (S13W13 L= 62)   291    113        0 (1147/ 845-1262/ 876)   8.1 W09-W16 S14-S12
     11  569  711 (S30E27 L= 22)  2302    469        7 ( 357/ 628- 720/ 856)  28.7 E47-E16 S36-S19
     12 1328  816 (S14W21 L= 70)   193     72        0 (1307/ 782-1361/ 846)  68.2 W20-W22 S16-S11
     13 1466  687 (S19W33 L= 82)   646    226        0 (1432/ 624-1530/ 780) -48.8 W29-W39 S22-S13
     14 1113  409 (S43W17 L= 66)  2547    450        8 ( 897/ 366-1329/ 459)  10.0 W00-W34 S45-S38
     15  797  389 (S49E08 L= 41)   268     41        0 ( 773/ 369- 819/ 411)  45.0 E10-E06 S51-S47
     16 1509  387 (S37W51 L=100)   305     13        1 (1493/ 367-1528/ 409) -21.8 W48-W53 S37-S35
     17 1589  373 (S36W61 L=110)   733     86        1 (1553/ 348-1628/ 398)  23.2 W54-W64 S37-S33
  

    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