LDBlockShow icon indicating copy to clipboard operation
LDBlockShow copied to clipboard

Blank png files

Open CNLove-260 opened this issue 1 year ago • 5 comments

I am attempting to plot regions along a genome. I know I have dense SNP data, but some regions work, and others don't? An example of the command for a region that does not work is:

LDBlockShow -InVCF concat.recode.vcf.gz -OutPut chr25_out_LD -Region chr25:46949317-47474693 -OutPng -NoShowLDist 1000000 -SeleVar 2
#Compress the SVG file
ShowLDSVG -InPreFix chr25_out_LD -OutPut chr25_out_LD_compressed.svg

Output File:

#warning skip Indel site, there are total skip Indel sites number is : 1309 #Warning skip non bi-allelic(Singleton/ThreeMulti allelic) site, and total skip allelic sites number is :3169 #Warning skip high missing site, and total skip allelic sites number is :22 #Warning skip low Minor Allele Frequency site, and total skip allelic sites number is :1164 ##Start Region Cal... :chr25 46949317 47474693; In This Region TotalSNP Number is 3217 In Big SNP Number :3217 ,Para -NumGradien suggest be maxValue : 10 ,auto be it convert SVG ---> PNG ... ALL done

Error File:

Start draw... SVG info: SNPNumber :3217 , SVG (width,height) = (5919.28,4375.12) Start draw... SVG info: SNPNumber :3217 , SVG (width,height) = (5919.28,4375.12) Start draw... SVG info: SNPNumber :4813 , SVG (width,height) = (6641.94,4909.26) Start draw... SVG info: SNPNumber :4813 , SVG (width,height) = (6641.94,4909.26)

Thank you for the help!

CNLove-260 avatar Aug 28 '24 03:08 CNLove-260

You can run the command convert xxx.svg xxx.png yourself convert xxx.svg xxx.png

hewm2008 avatar Aug 28 '24 03:08 hewm2008

Thank you for the reply! Sadly this did not working either. I still end up with a blank file. The file size increases a bit? Originally it was 7.3K but now it is 13K.

I have sence tried running a number of other regions as well, none of which are plotting.

Some throw an error:

Can't Found the Site file : chr25_out_LD.site.gz Too many SNP site 4813 , I suggest that you add the Para [-OutPng] to see the png File In Big SNP Number :4813 ,Para -NumGradien suggest be maxValue : 10 ,auto be it ALL done

And again, the command I used that resulted in this error was: LDBlockShow -InVCF dp8_concat.recode.vcf -OutPut chr25_whole_out_LD -Region chr25:1-700000 -OutPng -NoShowLDist 108720060 -SeleVar 2

CNLove-260 avatar Aug 28 '24 04:08 CNLove-260

svg is only 17k, which is not too big. There may be something wrong with the convert command in your linux system. I suggest you download svg to your local computer, open it with a browser, and then use other screenshot software to convert it to png format.

hewm2008 avatar Aug 28 '24 04:08 hewm2008

The svg file is actually 308M. I have tried bringing it locally, but does not like to open in the browser. The convert function seems to work for some regions, but not most.

CNLove-260 avatar Aug 28 '24 05:08 CNLove-260

You can retry running this command manually on a computer node with large memory.convert xxx.svg xxx.png

hewm2008 avatar Aug 28 '24 05:08 hewm2008