| 95 | } |
| 96 | |
| 97 | void UCSCWiggleTrackWriter::process(const Wiggle& wiggle) { |
| 98 | int sp, ep; |
| 99 | |
| 100 | if (wiggle.read_depth.empty()) return; |
| 101 | |
| 102 | sp = ep = -1; |
| 103 | for (size_t i = 0; i < wiggle.length; i++) { |
| 104 | if (wiggle.read_depth[i] >= 0.0095) { |
| 105 | ep = i; |
| 106 | } |
| 107 | else { |
| 108 | if (sp < ep) { |
| 109 | ++sp; |
| 110 | fprintf(fo, "fixedStep chrom=%s start=%d step=1\n", wiggle.name.c_str(), sp + 1); |
| 111 | for (int j = sp; j <= ep; j++) fprintf(fo, "%.2f\n", wiggle.read_depth[j]); |
| 112 | } |
| 113 | sp = i; |
| 114 | } |
| 115 | } |
| 116 | if (sp < ep) { |
| 117 | ++sp; |
| 118 | fprintf(fo, "fixedStep chrom=%s start=%d step=1\n", wiggle.name.c_str(), sp + 1); |
| 119 | for (int j = sp; j <= ep; j++) fprintf(fo, "%.2f\n", wiggle.read_depth[j]); |
| 120 | } |
| 121 | } |
| 122 | |
| 123 | ReadDepthWriter::ReadDepthWriter(std::ostream& stream) |
| 124 | : stream_(stream) { |
no test coverage detected