# R switch between point of segment and points of square in Hilberts Curve

Happy new year! Below I have my code in R which creates the Hilber curve of the 1st chromosome from a Bam file.

``````library(chromstaR)
library(GenomicRanges)
library(HilbertCurve)
library(circlize)
bampath <-"C:/Users/User/Documents/R/win-library/4.0/GenomicRanges/extdata/wgEncodeUwRepliSeqBg02esG1bAlnRep1.bam"
bam <- bampath

hc = GenomicHilbertCurve(chr = "chr1", level = 5)
hc_points(hc, reads, gp = gpar(col = "blue"))
``````

Now what I want to create is the following. I want from a point of a segment to find the coordinates in the Hilbert curve and the opposite i.e from a coordinate from the Hilbert Curve to find the point in the line segment.

I found from the web the following solution but is in js.

``````function hilbertMap(quadits) {
if ( quadits.length === 0 ) {
return new Point(0, 0);
} else {
return (function() {
var pt, t, x, y;
x = pt.x;
y = pt.y;
switch(t) {
case 0:
return new Point(y * 0.5 + 0, x * 0.5 + 0);
case 1:
return new Point(x * 0.5 + 0, y * 0.5 + 0.5);
case 2:
return new Point(x * 0.5 + 0.5, y * 0.5 + 0.5);
case 3:
return new Point(y * -0.5 + 1.0, x * -0.5 + 0.5);
}
})();
}
}

function hilbertMapInverse(x, y, depth) {
var t;
if ( depth > 1 ) {
return 0;
}
else {
if ( x < 0.5 ) {
if ( y < 0.5 ) {
return (hilbertMapInverse(y * 2, x * 2, depth * 4) + 0) / 4;
}
else {
return (hilbertMapInverse(x * 2, y * 2 - 1, depth * 4) + 1) / 4;
}
}
else {
if ( y >= 0.5 ) {
return (hilbertMapInverse(x * 2 - 1, y * 2 - 1, depth * 4) + 2) / 4;
}
else {
return (hilbertMapInverse(1 - y * 2, 2 - x * 2, depth * 4) + 3) / 4;
}
}
}
}

})();
``````

Please in need help to transform in R the code in js and to attach on my Hilbert curve. Thank you

This topic was automatically closed 21 days after the last reply. New replies are no longer allowed.

If you have a query related to it or one of the replies, start a new topic and refer back with a link.