Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  445 1579 (N26E56 L=241)   675     78        1 ( 423/1538- 467/1619) -63.4 E58-E53 N22-N28
      2 1418 1586 (N43W22 L=319)   380     62        2 (1388/1553-1441/1608) -76.0 W18-W24 N40-N45
      3  678 1485 (N23E31 L=266)   290    107        0 ( 639/1445- 707/1504) -18.4 E34-E28 N20-N24
      4 1701 1302 (N27W46 L=343)   315     67        1 (1679/1262-1721/1333) -63.4 W43-W47 N25-N29
      5 1427 1189 (N16W23 L=320)   416    169        0 (1369/1139-1477/1253) -29.1 W18-W27 N13-N19
      6  136 1065 (S10E73 L=224)  2996    168        7 ( 101/ 988- 180/1131) -56.3 E85-E65 S15-S06
      7  310  972 (S13E49 L=248)   420     85        1 ( 290/ 940- 329/1009)  90.0 E50-E46 S15-S11
      8 1561  773 (S07W39 L=336)   248     89        0 (1520/ 751-1601/ 797)  -8.1 W35-W42 S08-S06
      9  635  749 (S23E21 L=276)   236     66        0 ( 617/ 727- 667/ 779) -45.0 E22-E19 S23-S20
     10  255  747 (S29E57 L=240)   464     56        1 ( 240/ 721- 281/ 770)  23.2 E60-E53 S30-S27
     11 1375  741 (S12W27 L=324)   270    117        0 (1317/ 726-1425/ 766)  29.7 W22-W29 S13-S09
     12 1655  654 (S13W51 L=348)   771    176        2 (1619/ 581-1688/ 729) -71.6 W48-W53 S17-S08
     13 1559  647 (S15W43 L=340)   999    239        4 (1498/ 571-1599/ 699) -41.2 W36-W46 S20-S11
     14 1718  603 (S15W60 L=357)   664    130        0 (1664/ 567-1768/ 649)  36.9 W55-W63 S18-S11
     15  933  551 (S31W02 L=299)  1445    445        2 ( 746/ 501-1132/ 612)  27.5 E11-W14 S37-S24
     16 1240  301 (S44W35 L=332)  1414    260        2 (1140/ 273-1360/ 331)  26.6 W26-W43 S48-S39
  

    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