Module:Sandbox/AlexNB/nmColor: Difference between revisions

Content deleted Content added
AlexNB (talk | contribs)
another absorption algorithm: color wheel
AlexNB (talk | contribs)
No edit summary
 
(66 intermediate revisions by the same user not shown)
Line 28:
greenvalue = (645 - wavelength) / 65
bluevalue = 0
elseif (wavelength >= 645) and (wavelength <= 780) then
redvalue = 1
greenvalue = 0
Line 60:
-- code based on a color wheel drawn at: http://www2.chemistry.msu.edu/faculty/reusch/virttxtjml/spectrpy/uv-vis/spectrum.htm
local wavelength_vis
local minVioletIndigo = 380 -- originally, it was 400
local maxVioletIndigo = 430
local widthVioletIndigo = maxVioletIndigo - minVioletIndigo
Line 73:
local widthYellow = maxYellow - minYellow
local minOrange = 580
local maxOrange = 650 -- originally, it was 620
local widthOrange = maxOrange - minOrange
local minRed = 650 -- originally, it was 620
local maxRed = 780 -- originally, it was 800
local widthRed = maxRed - minRed
if (wavelength_abs>=minVioletIndigo) and (wavelength_abs<=maxVioletIndigo) then
wavelength_vis = (wavelength_abs - minVioletIndigo)*widthYellow/widthVioletIndigo + minYellow
elseif (wavelength_abs>minBlue) and (wavelength_abs<=maxBlue) then
Line 89:
wavelength_vis = (wavelength_abs - minOrange)*widthBlue/widthOrange + minBlue
elseif (wavelength_abs>minRed) and (wavelength_abs<=maxRed) then
wavelength_vis = (wavelength_abs - minRed)*widthGreen/(widthRed + 20) + minGreen
-- one of the "corrected" versions: wavelength_vis = (wavelength_abs - minRed - 30)*widthGreen/widthRed + minGreen
end
return wavelength_vis
end
 
local minAntiMagenta = 550
local maxAntiMagenta = 570
 
local function process_magenta(wavelength_abs)
local redsum = 0redvalue_RED
local bluevalue_RED
local greenvalue_RED
local redvalue_BLUE
local bluevalue_BLUE
local greenvalue_BLUE
local widthAntiMagenta = maxAntiMagenta - minAntiMagenta
local bias_BLUE = (wavelength_abs - minAntiMagenta)/widthAntiMagenta
local bias_RED = (maxAntiMagenta - wavelength_abs)/widthAntiMagenta
-- nm2RGB(wlcycle700, false)
redvalue_RED = redvalue*bias_RED
bluevalue_RED = bluevalue*bias_RED
greenvalue_RED = greenvalue*bias_RED
-- nm2RGB(wavelength400, false)
redvalue_BLUE = redvalue*bias_BLUE
bluevalue_BLUE = bluevalue*bias_BLUE
greenvalue_BLUE = greenvalue*bias_BLUE
redvalue = 0.6+ 0.4*(redvalue_RED + redvalue_BLUE)
bluevalue = 0.7 + 0.3*(bluevalue_RED + bluevalue_BLUE)
greenvalue = greenvalue_RED + greenvalue_BLUE
end
 
Line 97 ⟶ 124:
-- this function returns the string "#RRGGBB" with approximate RGB value for a complementary color (reflection) to a color, defined as a wavelength passed as a first argument
local wavelength = tonumber(frame.args[1])
-- adding magentas to smoothen red->blue transition. Trying 540..570 nm
nm2RGB(absorption_to_visible(wavelength), false)
if (wavelength>=minAntiMagenta) and (wavelength<=maxAntiMagenta) then
-- if wlcycle < process_magenta(wavelength then)
else nm2RGB(absorption_to_visible(wavelength), false)
-- end
local result='#' .. string.format("%.2X%.2X%.2X", 255*redvalue, 255*greenvalue, 255*bluevalue)
return result
end
 
function p.absorption_old(frame)
-- this function returns the string "#RRGGBB" with approximate RGB value for a complementary color (reflection) to a color, defined as a wavelength passed as a first argument
-- The following formula is used for the complimentary color: {0xFF - RR, 0xFF - (RR + BB)/2, 0xFF - BB}
-- local wavelength = tonumber(frame.args[1])
-- nm2RGB(wavelength)
-- local result='#' .. string.format("%.2X%.2X%.2X", 255*(1-redvalue), 255*(1-greenvalue), 255*(1-bluevalue))
local wavelength = tonumber(frame.args[1])
local redsum = 0
local greensum = 0
local bluesum = 0
local intensityfactor = 0
local numofitems = 400
-- for wlcycle = 380,780,1 do
-- nm2RGB(wlcycle, false)
-- if wlcycle < wavelength then
-- intensityfactor = (wlcycle - 380) / (wavelength - 380)
-- else
-- intensityfactor = (wlcycle - wavelength) / (380 - wavelength) + 1
-- end
-- if intensityfactor < 0 then
-- intensityfactor = 0
-- numofitems = numofitems - 1
-- end
-- redsum = redsum + redvalue*intensityfactor*1.8
-- greensum = greensum + greenvalue*intensityfactor*1.5
-- bluesum = bluesum + bluevalue*intensityfactor*1
-- end
-- local result='#' .. string.format("%.2X%.2X%.2X", 255*(redsum/numofitems), 255*(greensum/numofitems), 255*(bluesum/numofitems))
nm2RGB(1160-wavelength, true)
local result='#' .. string.format("%.2X%.2X%.2X", 255*redvalue, 255*bluevalue, 255*greenvalue)
return result
end
 
return p