Modulo:Sandbox/M.casanova/DatiMappa

require('Module:No globals')
local getArgs = require('Module:Arguments').getArgs

local p = {}

-- Calcola la scala adeguata per mappa OSM
-- esempio chiamata di funzione {{#invoke:Mappa OSM|scala|w|h|x1|y1|x2|y2|...}}
function p.scala(frame)
	local t = getArgs(frame)
	local amn, amx, bmn, bmx = 400, -400, 400, -400
	local dimx = tonumber(t[1]) or 1 --larghezza
	local dimy = tonumber(t[2]) or 1 --altezza
	local vlr
	for i,v in ipairs(t) do 
		if i>2 then
			vlr = tonumber(t[i])
			if (i % 2 == 0) then --calcolo latitudine minima e massima
				if (amn > vlr) then amn = vlr end
				if (amx < vlr) then amx = vlr end
			else --calcolo longitudine minima e massima
				if (bmn > vlr) then bmn = vlr end
				if (bmx < vlr) then bmx = vlr end
			end
		end
	end
	-- per evitare punti troppo vicini al margine
	local dx = 1.1*(bmx-bmn)
	local dy = 1.1*(amx-amn)
	if (dimx>dimy) then
		if (dx<180) then dx = dx*dimy/dimx end
	else
		if (dy<90) then dy = dy*dimx/dimy end
	end
	-- calcolo scale per coordinate
	local scalax, scalay
	if (dx == 0) then scalax = 18 else scalax = math.floor(math.log(360/dx)/math.log(2)) end
	if (dy == 0) then scalay = 18 else scalay = math.floor(math.log(180/dy)/math.log(2)) end
	local scala
    if ((dx == 0) and (dy == 0)) then
    	scala = 10 --valore default per singolo punto
   	else
   		scala = math.max(0,math.min(18,scalax, scalay)) --calcolo scala minima compresa tra 0 e 18
   	end
	return scala
end

-- Calcola la scala adeguata per mappa OSM
-- esempio chiamata di funzione {{#invoke:Mappa OSM|scala|w|h|x1|y1|x2|y2|...}}
function p.scala2(frame)
	local t = getArgs(frame)
	local lat1, lat2, lon1, lon2 = 400, -400, 400, -400
	local dimx = tonumber(t[1]) or 1 --larghezza
	local dimy = tonumber(t[2]) or 1 --altezza
	local vlr
	for i,v in ipairs(t) do 
		if i>2 then
			vlr = tonumber(t[i])
			if (i % 2 == 0) then --calcolo latitudine minima e massima
				if (lat1 > vlr) then lat1 = vlr end
				if (lat2 < vlr) then lat2 = vlr end
			else --calcolo longitudine minima e massima
				if (lon1 > vlr) then lon1 = vlr end
				if (lon2 < vlr) then lon2 = vlr end
			end
		end
	end

	-- problemi con latitudine oltre 85°
	if ((lat2>85) or (lat1<-85)) then return 0 end

	-- calcolo punti proiezione
--	local p_lon1 = (1+lon1/180)/2
--	local p_lat1 = (1-math.log(math.tan( math.pi*(1 + lat1/90)/4))/math.pi)/2
--	local p_lon2 = (1+lon2/180)/2
--	local p_lat2 = (1-math.log(math.tan( math.pi*(1 + lat2/90)/4))/math.pi)/2 ]]

	-- per evitare punti troppo vicini al margine
	local dx = 1.05*(lon2-lon1)/360
	local dy = 1.05*(math.log(math.tan( math.pi*(1 + lat2/90)/4)) - math.log(math.tan( math.pi*(1 + lat1/90)/4)))/(2*math.pi)
	-- calcolo scale per coordinate
	local scalax, scalay
	if (dx == 0) then scalax = 18 else scalax = math.floor(-math.log(dx)/math.log(2)) end
	if (dy == 0) then scalay = 18 else scalay = math.floor(-math.log(dy)/math.log(2)) end
	local scala
    if ((dx == 0) and (dy == 0)) then
    	scala = 10 --valore default per singolo punto
   	else
   		scala = math.max(0,math.min(18,scalax, scalay)) --calcolo scala minima compresa tra 0 e 18
   	end
	return scala
end

-- Calcola centro per singola coordinata
-- esempio chiamata di funzione {{#invoke:Mappa OSM|centro|x1|x2|...}}
function p.centro(frame)
	local t = getArgs(frame)
	local mn, mx = 400, -400
	local xval
	for i,v in ipairs(t) do
		xval = tonumber(t[i])
		if (mn > xval) then mn = xval end
		if (mx < xval) then mx = xval end
	end
    return (mn+mx)/2
end

return p