This is a solution in python version (3.0+) using a package called pyBigWig
to extract a given genomic region from a whole genome bigwig file.
Prepare your input bigwig file:
import pyBigWig
# First open bigwig file
bwIn = pyBigWig.open('input.bw')
# check bigwig file header
print(bwIn.header())
Prepare output, since your output doesn’t have a header, you need to add the header using the chosen chromosome size, here I’m using a region from chr18 as an example.
bwOutput = pyBigWig.open('output.bw','w')
bwOutput.addHeader([('chr18',80373285)]) # chromosome size
for x in bwIn.intervals('chr18',62926563,63516911):
bwOutput.addEntries(['chr18'],[x[0]],ends=[x[1]],values=[x[2]])
bwOutput.close()