R/BTC.R

In hexbin: Hexagonal Binning Routines

BTC <- function(n, beg = 1, end = 256)
{
    if(beg < 1 || end < 1 || beg > 256 || end > 256)
	stop("`beg' and `end' must be numbers in the interval [1,256]")

    M <- rbind(c(0,0,0),
	       c(0,0,40),
	       c(0,4,56),
	       c(0,9,61),
	       c(0,12,64),
	       c(0,14,66),
	       c(0,17,69),
	       c(0,20,73),
	       c(0,22,74),
	       c(0,25,78),
	       c(0,27,79),
	       c(0,30,83),
	       c(0,31,85),
	       c(0,33,86),
	       c(0,36,90),
	       c(0,38,91),
	       c(0,39,93),
	       c(0,41,95),
	       c(0,43,96),
	       c(0,46,100),
	       c(0,47,102),
	       c(0,49,103),
	       c(0,51,105),
	       c(0,52,107),
	       c(0,54,108),
	       c(0,55,110),
	       c(0,57,112),
	       c(0,57,112),
	       c(0,58,113),
	       c(0,60,115),
	       c(0,62,117),
	       c(0,63,119),
	       c(0,65,120),
	       c(0,66,122),
	       c(0,68,124),
	       c(0,70,125),
	       c(0,71,127),
	       c(0,73,129),
	       c(0,73,129),
	       c(0,74,130),
	       c(0,76,132),
	       c(0,78,134),
	       c(0,79,136),
	       c(0,81,137),
	       c(0,82,139),
	       c(0,84,141),
	       c(0,86,142),
	       c(0,87,144),
	       c(0,89,146),
	       c(0,90,147),
	       c(0,92,149),
	       c(0,94,151),
	       c(0,94,151),
	       c(0,95,153),
	       c(0,97,154),
	       c(0,98,156),
	       c(0,100,158),
	       c(0,102,159),
	       c(0,103,161),
	       c(0,105,163),
	       c(0,106,164),
	       c(0,108,166),
	       c(0,109,168),
	       c(0,111,170),
	       c(0,113,171),
	       c(0,114,173),
	       c(0,116,175),
	       c(0,117,176),
	       c(0,119,178),
	       c(0,121,180),
	       c(0,121,180),
	       c(0,122,181),
	       c(0,124,183),
	       c(0,125,185),
	       c(0,127,187),
	       c(0,129,188),
	       c(0,130,190),
	       c(0,132,192),
	       c(0,133,193),
	       c(0,135,195),
	       c(0,137,197),
	       c(0,138,198),
	       c(0,140,200),
	       c(0,141,202),
	       c(0,143,204),
	       c(0,143,204),
	       c(0,145,205),
	       c(0,146,207),
	       c(0,148,209),
	       c(0,149,210),
	       c(0,151,212),
	       c(0,153,214),
	       c(0,154,215),
	       c(0,156,217),
	       c(0,157,219),
	       c(0,159,221),
	       c(0,160,222),
	       c(0,160,222),
	       c(0,162,224),
	       c(0,164,226),
	       c(0,165,227),
	       c(0,167,229),
	       c(0,168,231),
	       c(0,170,232),
	       c(0,172,234),
	       c(0,173,236),
	       c(0,175,238),
	       c(0,175,238),
	       c(0,176,239),
	       c(0,178,241),
	       c(0,180,243),
	       c(0,181,244),
	       c(0,183,246),
	       c(2,184,248),
	       c(4,186,249),
	       c(4,186,249),
	       c(4,186,249),
	       c(6,188,251),
	       c(6,188,251),
	       c(9,189,253),
	       c(9,189,253),
	       c( 11,191,255),
	       c( 11,191,255),
	       c( 13,192,255),
	       c( 13,192,255),
	       c( 13,192,255),
	       c( 16,194,255),
	       c( 18,196,255),
	       c( 20,197,255),
	       c( 20,197,255),
	       c( 23,199,255),
	       c( 25,200,255),
	       c( 27,202,255),
	       c( 30,204,255),
	       c( 32,205,255),
	       c( 34,207,255),
	       c( 37,208,255),
	       c( 37,208,255),
	       c( 39,210,255),
	       c( 41,211,255),
	       c( 44,213,255),
	       c( 46,215,255),
	       c( 48,216,255),
	       c( 51,218,255),
	       c( 53,219,255),
	       c( 53,219,255),
	       c( 55,221,255),
	       c( 57,223,255),
	       c( 60,224,255),
	       c( 62,226,255),
	       c( 64,227,255),
	       c( 67,229,255),
	       c( 67,229,255),
	       c( 69,231,255),
	       c( 71,232,255),
	       c( 74,234,255),
	       c( 76,235,255),
	       c( 78,237,255),
	       c( 81,239,255),
	       c( 81,239,255),
	       c( 83,240,255),
	       c( 85,242,255),
	       c( 88,243,255),
	       c( 90,245,255),
	       c( 92,247,255),
	       c( 95,248,255),
	       c( 95,248,255),
	       c( 97,250,255),
	       c( 99,251,255),
	       c(102,253,255),
	       c(104,255,255),
	       c(106,255,255),
	       c(106,255,255),
	       c(108,255,255),
	       c(111,255,255),
	       c(113,255,255),
	       c(115,255,255),
	       c(115,255,255),
	       c(118,255,255),
	       c(120,255,255),
	       c(122,255,255),
	       c(122,255,255),
	       c(125,255,255),
	       c(127,255,255),
	       c(129,255,255),
	       c(129,255,255),
	       c(132,255,255),
	       c(134,255,255),
	       c(136,255,255),
	       c(136,255,255),
	       c(139,255,255),
	       c(141,255,255),
	       c(143,255,255),
	       c(143,255,255),
	       c(146,255,255),
	       c(148,255,255),
	       c(150,255,255),
	       c(150,255,255),
	       c(153,255,255),
	       c(155,255,255),
	       c(155,255,255),
	       c(157,255,255),
	       c(159,255,255),
	       c(159,255,255),
	       c(162,255,255),
	       c(164,255,255),
	       c(164,255,255),
	       c(166,255,255),
	       c(169,255,255),
	       c(171,255,255),
	       c(171,255,255),
	       c(173,255,255),
	       c(176,255,255),
	       c(176,255,255),
	       c(178,255,255),
	       c(180,255,255),
	       c(180,255,255),
	       c(183,255,255),
	       c(185,255,255),
	       c(185,255,255),
	       c(187,255,255),
	       c(190,255,255),
	       c(190,255,255),
	       c(192,255,255),
	       c(194,255,255),
	       c(197,255,255),
	       c(197,255,255),
	       c(199,255,255),
	       c(201,255,255),
	       c(204,255,255),
	       c(204,255,255),
	       c(206,255,255),
	       c(208,255,255),
	       c(210,255,255),
	       c(210,255,255),
	       c(213,255,255),
	       c(215,255,255),
	       c(217,255,255),
	       c(217,255,255),
	       c(220,255,255),
	       c(222,255,255),
	       c(224,255,255),
	       c(227,255,255),
	       c(229,255,255),
	       c(229,255,255),
	       c(231,255,255),
	       c(234,255,255),
	       c(236,255,255),
	       c(238,255,255),
	       c(241,255,255),
	       c(243,255,255),
	       c(243,255,255),
	       c(245,255,255),
	       c(248,255,255),
	       c(250,255,255),
	       c(255,255,255))[ round(seq(beg,end, length = n)), ]

    rgb(M[,1]/255,
	M[,2]/255,
	M[,3]/255)
}