Search in sources :

Example 1 with SVG

use of com.github.lindenb.jvarkit.util.svg.SVG in project jvarkit by lindenb.

the class VcfToSvg method doWork.

@Override
public int doWork(final List<String> args) {
    if (this.outputFile != null && !outputFile.getName().contains(SEGMENT)) {
        LOG.error("output file must contain the word " + SEGMENT + " :" + this.outputFile);
        return -1;
    }
    TabixKnownGeneFileReader tabix = null;
    VcfIterator r = null;
    OutputStream outputStream = null;
    XMLStreamWriter w = null;
    PrintWriter manifestW = null;
    try {
        LOG.info("opening knownGene ");
        tabix = new TabixKnownGeneFileReader(knownGeneUri);
        if (manifestFile != null && this.outputFile != null) {
            manifestW = new PrintWriter(manifestFile);
        } else {
            manifestW = new PrintWriter(new NullOuputStream());
        }
        final Set<String> chromosomes = tabix.getChromosomes();
        final XMLOutputFactory xof = XMLOutputFactory.newInstance();
        r = super.openVcfIterator(super.oneFileOrNull(args));
        final VCFHeader header = r.getHeader();
        while (r.hasNext()) {
            final VariantContext ctx = r.next();
            String tabixContig = ctx.getContig();
            if (!chromosomes.contains(tabixContig)) {
                if (tabixContig.startsWith("chr")) {
                    tabixContig = tabixContig.substring(3);
                } else if (!tabixContig.startsWith("chr")) {
                    tabixContig = "chr" + tabixContig;
                }
                if (!chromosomes.contains(tabixContig)) {
                    while (r.hasNext()) {
                        final VariantContext ctx2 = r.peek();
                        if (!ctx2.getContig().equals(ctx.getContig()))
                            break;
                        r.next();
                    }
                    LOG.error("No chromosome " + ctx.getContig() + " in " + knownGeneUri + ". Check the chromosome nomenclature.");
                    continue;
                }
            }
            final List<VariantContext> variants = new ArrayList<>();
            final List<KnownGene> genes = new ArrayList<>();
            variants.add(ctx);
            int chromStart = ctx.getStart() - 1;
            int chromEnd = ctx.getEnd();
            /* walk over know gene, loop until there is no overapping transcript 
			 * over that region */
            for (; ; ) {
                genes.clear();
                /* the max chromEnd, let's see if we can get a bigger */
                int newStart = chromStart;
                int newEnd = chromEnd;
                final Iterator<KnownGene> kgr = tabix.iterator(tabixContig, chromStart, chromEnd);
                while (kgr.hasNext()) {
                    final KnownGene g = kgr.next();
                    if (this.removeNonCoding && g.isNonCoding())
                        continue;
                    genes.add(g);
                    newStart = Math.min(g.getTxStart(), newStart);
                    newEnd = Math.max(g.getTxEnd(), newEnd);
                }
                if (newStart >= chromStart && newEnd <= chromEnd) {
                    break;
                }
                chromStart = newStart;
                chromEnd = newEnd;
            }
            // intergenic, no gene over that variant
            if (genes.isEmpty())
                continue;
            // fill the variant for that region
            while (r.hasNext()) {
                final VariantContext ctx2 = r.peek();
                if (!ctx2.getContig().equals(ctx.getContig()))
                    break;
                if (ctx2.getStart() > chromEnd)
                    break;
                variants.add(r.next());
            }
            if (this.variantsInExonOnly) {
                variants.removeIf(V -> {
                    for (final KnownGene gene : genes) {
                        for (final KnownGene.Exon exon : gene.getExons()) {
                            if (V.getEnd() < exon.getStart() || V.getStart() >= exon.getEnd()) {
                            // rien
                            } else {
                                return false;
                            }
                        }
                    }
                    return true;
                });
            }
            if (this.variantFILTEREDOpacity <= 0) {
                variants.removeIf(V -> V.isFiltered());
            }
            if (this.variantIndelOpacity <= 0) {
                variants.removeIf(V -> V.isIndel());
            }
            if (variants.isEmpty())
                continue;
            LOG.info("Variants (" + variants.size() + ") Transcripts (" + genes.size() + ") " + tabixContig + ":" + chromStart + "-" + chromEnd);
            if (outputFile != null) {
                File fname = new File(outputFile.getParentFile(), outputFile.getName().replaceAll("__SEGMENT__", ctx.getContig() + "_" + chromStart + "_" + chromEnd));
                LOG.info("saving as " + fname);
                outputStream = IOUtils.openFileForWriting(fname);
                w = xof.createXMLStreamWriter(outputStream);
                manifestW.println(ctx.getContig() + "\t" + chromStart + "\t" + chromEnd + "\t" + genes.stream().map(G -> G.getName()).collect(Collectors.joining(",")) + "\t" + genes.size() + "\t" + variants.size() + "\t" + fname);
            } else {
                w = xof.createXMLStreamWriter(stdout());
            }
            double featureHeight = 10;
            double TRANSCRIPT_HEIGHT = featureHeight;
            final int all_genotypes_width = variants.size() * this.genotype_width;
            if (trimToVariants) {
                chromStart = variants.stream().map(V -> V.getStart() - 1).min((A, B) -> A.compareTo(B)).get();
                chromEnd = variants.stream().map(V -> V.getEnd() + 1).max((A, B) -> A.compareTo(B)).get();
            }
            final int drawinAreaWidth = Math.max(all_genotypes_width, 1000);
            final Interval interval = new Interval(ctx.getContig(), chromStart, chromEnd);
            final int interline_weight = 6;
            final int margin_top = 10;
            final int margin_bottom = 10;
            final int margin_right = 100;
            final int margin_left = 100;
            w.writeStartDocument("UTF-8", "1.0");
            w.writeStartElement("svg");
            w.writeDefaultNamespace(SVG.NS);
            w.writeNamespace("xlink", XLINK.NS);
            w.writeAttribute("version", "1.1");
            w.writeAttribute("width", String.valueOf(margin_right + margin_right + drawinAreaWidth));
            w.writeAttribute("height", String.valueOf(margin_top + margin_bottom + genes.size() * TRANSCRIPT_HEIGHT + interline_weight * featureHeight + header.getSampleNamesInOrder().size() * this.genotype_width));
            title(w, ctx.getContig() + ":" + chromStart + "-" + chromEnd);
            w.writeStartElement("desc");
            w.writeCharacters("generated with " + getProgramName() + "\n" + "Author: Pierre Lindenbaum PhD. @yokofakun .");
            w.writeEndElement();
            // defs
            w.writeStartElement("defs");
            // genotypes
            w.writeStartElement("g");
            // 
            w.writeAttribute("id", "g_" + GenotypeType.HOM_REF);
            w.writeEmptyElement("rect");
            w.writeAttribute("style", "fill:lime;stroke;none;");
            w.writeAttribute("x", "0");
            w.writeAttribute("y", "0");
            w.writeAttribute("width", String.valueOf(this.genotype_width));
            w.writeAttribute("height", String.valueOf(this.genotype_width));
            w.writeEndElement();
            w.writeStartElement("g");
            // 
            w.writeAttribute("id", "g_" + GenotypeType.NO_CALL);
            w.writeEmptyElement("rect");
            w.writeAttribute("style", "fill:silver;stroke;gray;");
            w.writeAttribute("x", "0");
            w.writeAttribute("y", "0");
            w.writeAttribute("width", String.valueOf(this.genotype_width));
            w.writeAttribute("height", String.valueOf(this.genotype_width));
            w.writeEndElement();
            w.writeStartElement("g");
            // 
            w.writeAttribute("id", "g_" + GenotypeType.HOM_VAR);
            w.writeEmptyElement("rect");
            w.writeAttribute("style", "fill:crimson;stroke;none;");
            w.writeAttribute("x", "0");
            w.writeAttribute("y", "0");
            w.writeAttribute("width", String.valueOf(this.genotype_width));
            w.writeAttribute("height", String.valueOf(this.genotype_width));
            w.writeEndElement();
            w.writeStartElement("g");
            // 
            w.writeAttribute("id", "g_" + GenotypeType.MIXED);
            w.writeEmptyElement("rect");
            w.writeAttribute("style", "fill:pink;stroke;none;");
            w.writeAttribute("x", "0");
            w.writeAttribute("y", "0");
            w.writeAttribute("width", String.valueOf(this.genotype_width));
            w.writeAttribute("height", String.valueOf(this.genotype_width));
            w.writeEndElement();
            w.writeStartElement("g");
            // 
            w.writeAttribute("id", "g_" + GenotypeType.UNAVAILABLE);
            w.writeEmptyElement("rect");
            w.writeAttribute("style", "fill:gray;stroke;none;");
            w.writeAttribute("x", "0");
            w.writeAttribute("y", "0");
            w.writeAttribute("width", String.valueOf(this.genotype_width));
            w.writeAttribute("height", String.valueOf(this.genotype_width));
            w.writeEndElement();
            w.writeStartElement("g");
            // 
            w.writeAttribute("id", "g_" + GenotypeType.HET);
            w.writeEmptyElement("rect");
            w.writeAttribute("style", "fill:lime;stroke;black;");
            w.writeAttribute("x", "0");
            w.writeAttribute("y", "0");
            w.writeAttribute("width", String.valueOf(genotype_width));
            w.writeAttribute("height", String.valueOf(genotype_width));
            w.writeEmptyElement("polygon");
            w.writeAttribute("style", "fill:crimson;stroke;black;");
            w.writeAttribute("points", "0,0 " + genotype_width + ",0 0," + genotype_width + " 0,0");
            w.writeEndElement();
            // strand
            w.writeEmptyElement("polyline");
            w.writeAttribute("id", "strandF");
            w.writeAttribute("points", "-5,-5 0,0 -5,5");
            w.writeEmptyElement("polyline");
            w.writeAttribute("id", "strandR");
            w.writeAttribute("points", "5,-5 0,0 5,5");
            // gradients
            w.writeStartElement("linearGradient");
            w.writeAttribute("id", "grad01");
            w.writeAttribute("x1", "50%");
            w.writeAttribute("x2", "50%");
            w.writeAttribute("y1", "0%");
            w.writeAttribute("y2", "100%");
            w.writeEmptyElement("stop");
            w.writeAttribute("offset", "0%");
            w.writeAttribute("style", "stop-color:black;stop-opacity:1;");
            w.writeEmptyElement("stop");
            w.writeAttribute("offset", "50%");
            w.writeAttribute("style", "stop-color:white;stop-opacity:1;");
            w.writeEmptyElement("stop");
            w.writeAttribute("offset", "100%");
            w.writeAttribute("style", "stop-color:black;stop-opacity:1;");
            w.writeEndElement();
            // defs
            w.writeEndElement();
            w.writeStartElement("style");
            w.writeCharacters("svg {fill:none; stroke:black;}\n" + "text {fill:black;stroke:none;font-size:" + (featureHeight / 1.5) + "px;}\n" + ".ruler-label { stroke:red;}\n" + ".frame { stroke:black;fill:none;}\n" + ".kgexon {fill:url(#grad01);stroke:black;}\n" + ".gcpercent {fill:url(#grad02);stroke:black;}" + ".coverage {fill:url(#grad03);stroke:black;}" + ".kgcds {fill:yellow;stroke:black;opacity:0.7;}\n" + ".variant{stroke:none;fill:red;opacity:0.2;}\n" + ".xaxis{stroke:gray;fill:none;opacity:0.2;}\n" + ".postick{font-size:9px;stroke:black;stroke-width:1;}");
            // style
            w.writeEndElement();
            final Function<Integer, Integer> trim = new Function<Integer, Integer>() {

                @Override
                public Integer apply(final Integer t) {
                    return Math.max(interval.getStart(), Math.min(interval.getEnd(), t));
                }
            };
            final Function<Integer, Double> baseToPixel = new Function<Integer, Double>() {

                @Override
                public Double apply(final Integer t) {
                    return margin_left + drawinAreaWidth * (t - (double) interval.getStart()) / ((double) interval.length());
                }
            };
            final Function<Integer, Double> variantIndexToPixel = new Function<Integer, Double>() {

                @Override
                public Double apply(final Integer idx) {
                    final double variant_width = drawinAreaWidth / (double) variants.size();
                    final double midx = variant_width * idx + variant_width / 2.0;
                    return margin_left + midx - genotype_width / 2.0;
                }
            };
            final Function<VariantContext, String> variantTitle = V -> (V.getContig().startsWith("chr") ? V.getContig().substring(3) : V.getContig()) + ":" + V.getStart() + " " + V.getReference().getDisplayString();
            /**
             * title
             */
            double y = 0;
            w.writeStartElement("text");
            w.writeAttribute("x", "0");
            w.writeAttribute("y", String.valueOf(featureHeight));
            w.writeCharacters(interval.toString());
            w.writeEndElement();
            y += featureHeight;
            for (final KnownGene g : genes) {
                int cdsHeigh = 5;
                double exonHeight = TRANSCRIPT_HEIGHT - 5;
                double midY = TRANSCRIPT_HEIGHT / 2;
                w.writeStartElement("g");
                w.writeAttribute("transform", "translate(0," + y + ")");
                title(w, g.getName());
                w.writeStartElement("text");
                w.writeAttribute("x", String.valueOf(margin_left - 10));
                w.writeAttribute("y", String.valueOf(featureHeight));
                w.writeAttribute("style", "text-anchor:end;");
                w.writeCharacters(g.getName());
                w.writeEndElement();
                /* transcript line */
                w.writeEmptyElement("line");
                w.writeAttribute("class", "kgtr");
                w.writeAttribute("x1", String.valueOf(baseToPixel.apply(trim.apply(g.getTxStart()))));
                w.writeAttribute("y1", String.valueOf(midY));
                w.writeAttribute("x2", String.valueOf(baseToPixel.apply(trim.apply(g.getTxEnd()))));
                w.writeAttribute("y2", String.valueOf(midY));
                /* strand symbols */
                for (double pixX = 0; pixX < drawinAreaWidth; pixX += 30) {
                    double pos0 = interval.getStart() + (pixX / (double) drawinAreaWidth) * interval.length();
                    if (pos0 + 1 < g.getTxStart())
                        continue;
                    if (pos0 > g.getTxEnd())
                        break;
                    w.writeEmptyElement("use");
                    w.writeAttribute("class", "kgstrand");
                    w.writeAttribute("xlink", XLINK.NS, "href", "#strand" + (g.isPositiveStrand() ? "F" : "R"));
                    w.writeAttribute("x", String.valueOf(margin_left + pixX));
                    w.writeAttribute("y", String.valueOf(midY));
                }
                /* exons */
                for (KnownGene.Exon exon : g.getExons()) {
                    if (exon.getStart() + 1 >= interval.getEnd())
                        continue;
                    if (exon.getEnd() <= interval.getStart())
                        continue;
                    w.writeStartElement("rect");
                    w.writeAttribute("class", "kgexon");
                    w.writeAttribute("x", String.valueOf(baseToPixel.apply(trim.apply(exon.getStart()))));
                    w.writeAttribute("y", String.valueOf(midY - exonHeight / 2));
                    w.writeAttribute("width", String.valueOf(baseToPixel.apply(trim.apply(exon.getEnd())) - baseToPixel.apply((trim.apply(exon.getStart())))));
                    w.writeAttribute("height", String.valueOf(exonHeight));
                    title(w, exon.getName());
                    w.writeEndElement();
                }
                /* coding line */
                if (!g.isNonCoding()) {
                    w.writeEmptyElement("rect");
                    w.writeAttribute("class", "kgcds");
                    w.writeAttribute("x", String.valueOf(baseToPixel.apply(trim.apply(g.getCdsStart()))));
                    w.writeAttribute("y", String.valueOf(midY - cdsHeigh / 4.0));
                    w.writeAttribute("width", String.valueOf(baseToPixel.apply(trim.apply(g.getCdsEnd())) - baseToPixel.apply((trim.apply((g.getCdsStart()))))));
                    w.writeAttribute("height", String.valueOf(cdsHeigh / 2.0));
                }
                // String label=String.format("%15s", g.getName());
                // w.writeEmptyElement("path");
                // double fontHeight=Math.min(10,0.8*TRANSCRIPT_HEIGHT);
                // w.writeAttribute("d",this.hershey.svgPath(label,-insets.left,midY-fontHeight/2,insets.left*0.9,fontHeight));
                w.writeEndElement();
                w.writeCharacters("\n");
                y += featureHeight;
            }
            /* draw lines to variants */
            for (int vidx = 0; vidx < variants.size(); ++vidx) {
                final VariantContext vc = variants.get(vidx);
                double x1 = baseToPixel.apply(vc.getStart());
                double x2 = baseToPixel.apply(vc.getEnd());
                final double y2 = y + featureHeight * interline_weight;
                w.writeStartElement("polygon");
                w.writeAttribute("style", "fill:" + (vidx % 2 == 0 ? "ghostwhite" : "lavender") + ";stroke:black;opacity:0.6;stroke-width:0.5;");
                w.writeAttribute("points", "" + x1 + "," + (y - featureHeight / 2.0) + " " + x2 + "," + (y - featureHeight / 2.0) + " " + variantIndexToPixel.apply(vidx) + "," + y2 + " " + (variantIndexToPixel.apply(vidx) + this.genotype_width) + "," + y2);
                title(w, variantTitle.apply(vc));
                w.writeEndElement();
            }
            for (int vidx = 0; vidx < variants.size(); ++vidx) {
                final VariantContext vc = variants.get(vidx);
                final double y2 = y + featureHeight * interline_weight;
                w.writeStartElement("text");
                w.writeAttribute("transform", "translate(" + (String.valueOf(variantIndexToPixel.apply(vidx) + genotype_width / 2.0)) + "," + String.valueOf(y2 - 5) + ") " + "rotate(-45)");
                w.writeAttribute("x", "0");
                w.writeAttribute("y", "0");
                w.writeAttribute("class", "postick");
                w.writeCharacters(variantTitle.apply(vc));
                w.writeEndElement();
                w.writeCharacters("\n");
            }
            y += featureHeight * interline_weight;
            w.writeStartElement("g");
            for (final String sample : header.getSampleNamesInOrder()) {
                for (int vidx = 0; vidx < variants.size(); ++vidx) {
                    final VariantContext vc = variants.get(vidx);
                    final Genotype g = vc.getGenotype(sample);
                    double opacity = 1.0;
                    if (vc.isIndel())
                        opacity *= this.variantIndelOpacity;
                    if (vc.isFiltered())
                        opacity *= this.variantFILTEREDOpacity;
                    if (opacity > 1)
                        opacity = 1;
                    if (opacity <= 0)
                        continue;
                    if (opacity < 1) {
                        w.writeStartElement("g");
                        w.writeAttribute("style", "opacity:" + opacity + ";");
                    }
                    w.writeEmptyElement("use");
                    w.writeAttribute("x", "" + variantIndexToPixel.apply(vidx));
                    w.writeAttribute("y", String.valueOf(y));
                    w.writeAttribute("xlink", XLINK.NS, "href", "#g_" + g.getType());
                    if (opacity < 1) {
                        w.writeEndElement();
                    }
                }
                w.writeCharacters("\n");
                w.writeStartElement("text");
                w.writeAttribute("x", String.valueOf(margin_left - 10));
                w.writeAttribute("y", String.valueOf(y + this.genotype_width / 2.0));
                w.writeAttribute("style", "text-anchor:end;");
                w.writeCharacters(sample);
                w.writeEndElement();
                y += this.genotype_width;
            }
            w.writeCharacters("\n");
            w.writeEndDocument();
            w.writeCharacters("\n");
            w.flush();
            w.close();
            if (outputFile != null) {
                outputStream.flush();
                outputStream.close();
                outputStream = null;
            }
            if (stop_at_first) {
                LOG.info("Stop after first SVG document");
                break;
            }
        }
        r.close();
        manifestW.flush();
        manifestW.close();
        manifestW = null;
        return 0;
    } catch (Throwable err) {
        LOG.error(err);
        return -1;
    } finally {
        CloserUtil.close(r);
        CloserUtil.close(tabix);
        CloserUtil.close(outputStream);
        CloserUtil.close(manifestW);
    }
}
Also used : Genotype(htsjdk.variant.variantcontext.Genotype) XLINK(com.github.lindenb.jvarkit.util.ns.XLINK) Program(com.github.lindenb.jvarkit.util.jcommander.Program) Parameter(com.beust.jcommander.Parameter) NullOuputStream(com.github.lindenb.jvarkit.io.NullOuputStream) VCFHeader(htsjdk.variant.vcf.VCFHeader) Function(java.util.function.Function) SVG(com.github.lindenb.jvarkit.util.svg.SVG) ArrayList(java.util.ArrayList) Interval(htsjdk.samtools.util.Interval) XMLStreamException(javax.xml.stream.XMLStreamException) IOUtils(com.github.lindenb.jvarkit.io.IOUtils) Launcher(com.github.lindenb.jvarkit.util.jcommander.Launcher) XMLStreamWriter(javax.xml.stream.XMLStreamWriter) CloserUtil(htsjdk.samtools.util.CloserUtil) OutputStream(java.io.OutputStream) PrintWriter(java.io.PrintWriter) Iterator(java.util.Iterator) Logger(com.github.lindenb.jvarkit.util.log.Logger) VcfIterator(com.github.lindenb.jvarkit.util.vcf.VcfIterator) GenotypeType(htsjdk.variant.variantcontext.GenotypeType) Set(java.util.Set) Collectors(java.util.stream.Collectors) KnownGene(com.github.lindenb.jvarkit.util.ucsc.KnownGene) TabixKnownGeneFileReader(com.github.lindenb.jvarkit.util.ucsc.TabixKnownGeneFileReader) File(java.io.File) XMLOutputFactory(javax.xml.stream.XMLOutputFactory) List(java.util.List) VariantContext(htsjdk.variant.variantcontext.VariantContext) XMLOutputFactory(javax.xml.stream.XMLOutputFactory) OutputStream(java.io.OutputStream) ArrayList(java.util.ArrayList) VariantContext(htsjdk.variant.variantcontext.VariantContext) Function(java.util.function.Function) VcfIterator(com.github.lindenb.jvarkit.util.vcf.VcfIterator) XMLStreamWriter(javax.xml.stream.XMLStreamWriter) NullOuputStream(com.github.lindenb.jvarkit.io.NullOuputStream) VCFHeader(htsjdk.variant.vcf.VCFHeader) PrintWriter(java.io.PrintWriter) Genotype(htsjdk.variant.variantcontext.Genotype) TabixKnownGeneFileReader(com.github.lindenb.jvarkit.util.ucsc.TabixKnownGeneFileReader) KnownGene(com.github.lindenb.jvarkit.util.ucsc.KnownGene) File(java.io.File) Interval(htsjdk.samtools.util.Interval)

Example 2 with SVG

use of com.github.lindenb.jvarkit.util.svg.SVG in project jvarkit by lindenb.

the class CytobandToSvg method doWork.

@Override
public int doWork(final List<String> args) {
    XMLStreamWriter w = null;
    FileOutputStream fout = null;
    try {
        final List<Cytoband> cytobands;
        if (StringUtil.isBlank(this.cytobandsUri) || this.cytobandsUri.equals("-")) {
            cytobands = parseCytobands(IOUtils.openStdinForBufferedReader());
        } else {
            cytobands = parseCytobands(IOUtils.openURIForBufferedReading(this.cytobandsUri));
        }
        cytobands.stream().map(C -> C.getContig()).collect(Collectors.toCollection(LinkedHashSet::new)).stream().map(C -> new Contig(cytobands.stream().filter(T -> T.getContig().equals(C)).sorted((A, B) -> Integer.compare(A.getStart(), B.getStart())).collect(Collectors.toList()))).forEach(C -> this.name2contig.put(C.getContig(), C));
        ;
        this.max_contig_length = this.name2contig.values().stream().mapToLong(C -> C.getSequenceLength()).max().orElse(0L);
        if (this.name2contig.isEmpty() || this.max_contig_length <= 0) {
            LOG.error("no contig found");
            return -1;
        }
        final Rectangle drawingArea = new Rectangle(0, 0, this.width, this.height);
        final XMLOutputFactory xof = XMLOutputFactory.newFactory();
        if (this.outputFile == null) {
            w = xof.createXMLStreamWriter(stdout(), "UTF-8");
        } else {
            fout = new FileOutputStream(this.outputFile);
            w = xof.createXMLStreamWriter(fout, "UTF-8");
        }
        w.writeStartDocument("UTF-8", "1.0");
        w.writeStartElement("svg");
        w.writeAttribute("width", format(drawingArea.width));
        w.writeAttribute("height", format(drawingArea.height));
        w.writeAttribute("version", "1.1");
        w.writeDefaultNamespace(SVG.NS);
        w.writeNamespace("xlink", XLINK.NS);
        w.writeStartElement("title");
        w.writeCharacters(StringUtil.isBlank(this.title) ? CytobandToSvg.class.getName() : this.title);
        w.writeEndElement();
        w.writeStartElement("description");
        w.writeCharacters("Cmd:" + getProgramCommandLine() + "\n");
        w.writeCharacters("Version:" + getVersion() + "\n");
        w.writeCharacters("Author: Pierre Lindenbaum\n");
        w.writeEndElement();
        w.writeStartElement("defs");
        w.writeStartElement("linearGradient");
        w.writeAttribute("id", "grad01");
        w.writeAttribute("x1", "0%");
        w.writeAttribute("x2", "100%");
        w.writeAttribute("y1", "50%");
        w.writeAttribute("y2", "50%");
        w.writeEmptyElement("stop");
        w.writeAttribute("offset", "0%");
        w.writeAttribute("style", "stop-color:white;stop-opacity:1;");
        w.writeEmptyElement("stop");
        w.writeAttribute("offset", "50%");
        w.writeAttribute("style", "stop-color:white;stop-opacity:0;");
        w.writeEmptyElement("stop");
        w.writeAttribute("offset", "75%");
        w.writeAttribute("style", "stop-color:black;stop-opacity:0;");
        w.writeEmptyElement("stop");
        w.writeAttribute("offset", "100%");
        w.writeAttribute("style", "stop-color:black;stop-opacity:0.9;");
        w.writeEndElement();
        w.writeStartElement("filter");
        w.writeAttribute("id", "filter01");
        w.writeEmptyElement("feGaussianBlur");
        w.writeAttribute("in", "SourceGraphic");
        w.writeAttribute("stdDeviation", "10");
        // filter
        w.writeEndElement();
        double x = drawingArea.getX();
        double contig_width = drawingArea.getWidth() / (double) this.name2contig.size();
        for (final Contig contig : this.name2contig.values()) {
            contig.bounds.x = x;
            contig.bounds.y = drawingArea.getY();
            contig.bounds.width = contig_width;
            contig.bounds.height = drawingArea.getHeight();
            contig.def(w);
            x += contig_width;
        }
        w.writeEndElement();
        w.writeStartElement("style");
        w.writeCharacters(".ctgborderp {fill:url(#grad01);stroke:green;}" + ".ctgborderq {fill:url(#grad01);stroke:green;}" + ".ctglabel {text-anchor:middle;stroke:none;fill:darkgrey;font: bold 10px Verdana, Helvetica, Arial, sans-serif;}" + ".cytoband {fill:silver;stroke:none;}" + ".bedlabel {stroke:red;fill:none;text-anchor:start;font: normal 7px Verdana, Helvetica, Arial, sans-serif;}" + ".maintitle {stroke:none;fill:darkgrey;text-anchor:middle;font: normal 12px Verdana, Helvetica, Arial, sans-serif;}" + ".ctgback {fill:gainsboro;stroke:none;filter:url(#filter01);}");
        // style
        w.writeEndElement();
        w.writeStartElement("g");
        w.writeAttribute("id", "body");
        /* title */
        if (!StringUtil.isBlank(this.title)) {
            w.writeStartElement("text");
            w.writeAttribute("class", "maintitle");
            w.writeAttribute("x", format(this.width / 2.0));
            w.writeAttribute("y", format(12));
            w.writeCharacters(this.title);
            w.writeEndElement();
        }
        w.writeStartElement("g");
        w.writeAttribute("id", "karyotype");
        for (final Contig contig : this.name2contig.values()) {
            contig.paint(w);
        }
        // g@id=karyotype
        w.writeEndElement();
        final ContigNameConverter ctgNameConverter = ContigNameConverter.fromContigSet(this.name2contig.keySet());
        ctgNameConverter.setOnNotFound(ContigNameConverter.OnNotFound.SKIP);
        for (final String filename : args) {
            final Pattern tab = Pattern.compile("[\t]");
            w.writeStartElement("g");
            w.writeAttribute("id", "input" + filename);
            BufferedReader br = IOUtils.openURIForBufferedReading(filename);
            String line;
            while ((line = br.readLine()) != null) {
                if (line == null || StringUtil.isBlank(line) || BedLine.isBedHeader(line))
                    continue;
                final String[] tokens = tab.split(line, 4);
                if (tokens.length < 3)
                    continue;
                final String ctgName = ctgNameConverter.apply(tokens[0]);
                if (StringUtil.isBlank(ctgName))
                    continue;
                final Contig ctg = this.name2contig.get(ctgName);
                if (ctg == null)
                    continue;
                final Rectangle2D.Double r = ctg.getContigBounds();
                final int chromStart = Integer.parseInt(tokens[1]);
                final int chromEnd = Integer.parseInt(tokens[2]);
                final Map<String, String> attributes = new HashMap<>();
                if (tokens.length > 3) {
                    final StreamTokenizer st = new StreamTokenizer(new StringReader(tokens[3]));
                    st.wordChars('_', '_');
                    String key = null;
                    while (st.nextToken() != StreamTokenizer.TT_EOF) {
                        String s = null;
                        switch(st.ttype) {
                            case StreamTokenizer.TT_NUMBER:
                                s = String.valueOf(st.nval);
                                break;
                            case '"':
                            case '\'':
                            case StreamTokenizer.TT_WORD:
                                s = st.sval;
                                break;
                            case ';':
                                break;
                            default:
                                break;
                        }
                        if (s == null)
                            continue;
                        if (key == null) {
                            key = s;
                        } else {
                            attributes.put(key, s);
                            key = null;
                        }
                    }
                }
                w.writeStartElement("g");
                if (attributes.containsKey("href")) {
                    w.writeStartElement("a");
                    w.writeAttribute("href", attributes.get("href"));
                }
                w.writeStartElement("rect");
                w.writeAttribute("style", "stroke:red;");
                w.writeAttribute("x", format(r.getMaxX() + 2));
                w.writeAttribute("width", format(2));
                w.writeAttribute("y", format(ctg.pos2pixel(chromStart, r)));
                w.writeAttribute("height", format(ctg.pos2pixel(chromEnd, r) - ctg.pos2pixel(chromStart, r)));
                if (attributes.containsKey("title")) {
                    w.writeStartElement("title");
                    w.writeCharacters(attributes.get("title"));
                    w.writeEndElement();
                }
                w.writeEndElement();
                if (attributes.containsKey("label")) {
                    w.writeStartElement("text");
                    w.writeAttribute("class", "bedlabel");
                    w.writeAttribute("x", format(r.getMaxX() + 4));
                    w.writeAttribute("y", format(3 + ctg.pos2pixel(chromStart, r)));
                    w.writeCharacters(attributes.get("label"));
                    w.writeEndElement();
                }
                if (attributes.containsKey("href")) {
                    w.writeEndElement();
                }
                // g
                w.writeEndElement();
            }
            br.close();
            w.writeEndElement();
        }
        // g@id=body
        w.writeEndElement();
        // svg
        w.writeEndElement();
        w.writeEndDocument();
        w.flush();
        w.close();
        if (fout != null) {
            fout.flush();
            fout.close();
            fout = null;
        }
        return 0;
    } catch (final Exception e) {
        LOG.error(e);
        return -1;
    } finally {
        CloserUtil.close(fout);
    }
}
Also used : LinkedHashSet(java.util.LinkedHashSet) Rectangle(java.awt.Rectangle) ContigNameConverter(com.github.lindenb.jvarkit.util.bio.fasta.ContigNameConverter) XLINK(com.github.lindenb.jvarkit.util.ns.XLINK) Program(com.github.lindenb.jvarkit.util.jcommander.Program) Parameter(com.beust.jcommander.Parameter) Rectangle2D(java.awt.geom.Rectangle2D) HashMap(java.util.HashMap) SVG(com.github.lindenb.jvarkit.util.svg.SVG) LinkedHashMap(java.util.LinkedHashMap) StringUtil(htsjdk.samtools.util.StringUtil) Map(java.util.Map) XMLStreamException(javax.xml.stream.XMLStreamException) IOUtils(com.github.lindenb.jvarkit.io.IOUtils) Launcher(com.github.lindenb.jvarkit.util.jcommander.Launcher) StreamTokenizer(java.io.StreamTokenizer) XMLStreamWriter(javax.xml.stream.XMLStreamWriter) LinkedHashSet(java.util.LinkedHashSet) CloserUtil(htsjdk.samtools.util.CloserUtil) Locatable(htsjdk.samtools.util.Locatable) Logger(com.github.lindenb.jvarkit.util.log.Logger) DecimalFormat(java.text.DecimalFormat) FileOutputStream(java.io.FileOutputStream) IOException(java.io.IOException) Collectors(java.util.stream.Collectors) File(java.io.File) XMLOutputFactory(javax.xml.stream.XMLOutputFactory) List(java.util.List) StringReader(java.io.StringReader) BufferedReader(java.io.BufferedReader) Pattern(java.util.regex.Pattern) BedLine(com.github.lindenb.jvarkit.util.bio.bed.BedLine) Pattern(java.util.regex.Pattern) XMLOutputFactory(javax.xml.stream.XMLOutputFactory) HashMap(java.util.HashMap) LinkedHashMap(java.util.LinkedHashMap) Rectangle(java.awt.Rectangle) Rectangle2D(java.awt.geom.Rectangle2D) XMLStreamException(javax.xml.stream.XMLStreamException) IOException(java.io.IOException) XMLStreamWriter(javax.xml.stream.XMLStreamWriter) FileOutputStream(java.io.FileOutputStream) BufferedReader(java.io.BufferedReader) StringReader(java.io.StringReader) ContigNameConverter(com.github.lindenb.jvarkit.util.bio.fasta.ContigNameConverter) StreamTokenizer(java.io.StreamTokenizer)

Aggregations

Parameter (com.beust.jcommander.Parameter)2 IOUtils (com.github.lindenb.jvarkit.io.IOUtils)2 Launcher (com.github.lindenb.jvarkit.util.jcommander.Launcher)2 Program (com.github.lindenb.jvarkit.util.jcommander.Program)2 Logger (com.github.lindenb.jvarkit.util.log.Logger)2 XLINK (com.github.lindenb.jvarkit.util.ns.XLINK)2 SVG (com.github.lindenb.jvarkit.util.svg.SVG)2 CloserUtil (htsjdk.samtools.util.CloserUtil)2 File (java.io.File)2 List (java.util.List)2 Collectors (java.util.stream.Collectors)2 XMLOutputFactory (javax.xml.stream.XMLOutputFactory)2 XMLStreamException (javax.xml.stream.XMLStreamException)2 XMLStreamWriter (javax.xml.stream.XMLStreamWriter)2 NullOuputStream (com.github.lindenb.jvarkit.io.NullOuputStream)1 BedLine (com.github.lindenb.jvarkit.util.bio.bed.BedLine)1 ContigNameConverter (com.github.lindenb.jvarkit.util.bio.fasta.ContigNameConverter)1 KnownGene (com.github.lindenb.jvarkit.util.ucsc.KnownGene)1 TabixKnownGeneFileReader (com.github.lindenb.jvarkit.util.ucsc.TabixKnownGeneFileReader)1 VcfIterator (com.github.lindenb.jvarkit.util.vcf.VcfIterator)1