Queria que meu mapa aparecesse em uma janela de largura 600px e altura 430px, com o zoom que fizesse uma certa feature aproveitar esse espaço da melhor maneira, ficando totalmente visível nessa área.
Encontrei a seguinte fórmula:
Map resolution = 156543.04 meters/pixel * cos(latitude) / (2 ^ zoomlevel)
Fui jogando as coisas de um lado pro outro:
(2 ^ zoomlevel) = 156543.04 meters/pixel * cos(latitude)/Map resolution
E enfim, considerando que se 2^x = a, então x = (ln a/ln 2):
zoomlevel = ln(156543.04 meters/pixel * cos(latitude)/Map resolution)/ln 2
Como latitude, usei a função st_ymin (no meu caso, seria praticamente igual a usar o st_ymax) do PostGIS, passando a feature como parâmetro. Para calcular a resolução, tive que descobrir se a feição era maior no eixo x ou no eixo y usando a seguinte query:
select st_distance_sphere(st_makepoint(st_xmax(the_geom), st_ymin(the_geom)), st_makepoint(st_xmin(the_geom), st_ymin(the_geom))) as dist_x,
st_distance_sphere(st_makepoint(st_xmin(the_geom), st_ymax(the_geom)), st_makepoint(st_xmin(the_geom), st_ymin(the_geom))) as dist_y
from <tabela> where <cláusula where>
Como já sei o tamanho da janela onde ficará o mapa (600 x 430), então se a feature é maior no eixo x, a resolução é dist_x/600. Se a feature é maior no eixo y, a resolução é dist_y/430.
Encontrei a fórmula para calcular o zoom em:
http://msdn.microsoft.com/en-us/library/aa940990.aspx
Obrigado ao Cristiano Pereira da Silva (ardagh) pela ajuda com a matemática.
Nenhum comentário:
Postar um comentário