Vylepšení funkce pro nadmořskou výšku – MS SQL – Fórum – Programujte.com
 x   TIP: Přetáhni ikonu na hlavní panel pro připnutí webu
Reklama
Reklama

Vylepšení funkce pro nadmořskou výšku – MS SQL – Fórum – Programujte.comVylepšení funkce pro nadmořskou výšku – MS SQL – Fórum – Programujte.com

 

Hledá se programátor! Plat 1 800 € + bonusy (firma Boxmol.com)
soadfan0
Duch
4. 5. 2015   #1
-
0
-

Zdravím Vás
dostal jsem úkol - musím "vylepšit" funkci, kterou přikládám níže:

cd ~/test/freemap/elevation
scp -p -q -P 21122 92.240.244.41:/home/michal/elevation.sql.bz2 .

psql:
create table public.fm_elevation (
x int, y int,
p00 real, p01 real, p02 real, p03 real, 
p10 real, p11 real, p12 real, p13 real,
p20 real, p21 real, p22 real, p23 real,
p30 real, p31 real, p32 real, p33 real);


bzcat elevation.sql.bz2 | grep -v '/*!' |grep -v 'LOCK' | sed 's/`//g' |sql mapnik

-- x = floor($Lat * 1200);  => lat = x/1200

create table public.elevation as 
select geography(st_makepoint(y/1200.0, x/1200.0)) as p, p11 as ele
from public.fm_elevation;
create index on elevation using gist(p);
analyze elevation;

select count(*) from elevation  where st_covers((select st_buffer(geography(waym),1000) from oblasti where typ='sr'), p);

create table t_elevation as 
select geography(st_makepoint(y/1200.0, x/1200.0)) as p, p11 as ele
from public.fm_elevation where st_covers((select st_buffer(geography(waym),2000) from oblasti where typ='sr'), geography(st_makepoint(y/1200.0, x/1200.0)));
create index on t_elevation using gist(p);
analyze t_elevation;


create or replace function getz(geography) returns real 
LANGUAGE sql IMMUTABLE  AS $$
--select ele from elevation where p && st_buffer($1, 150) order by st_distance(p, $1) asc limit 1;
select case when sum(100/(0.01+st_distance(p, $1))) > 0.5 then (sum(ele*100/(0.01+st_distance(p, $1)))/sum(100/(0.01+st_distance(p, $1))))::real else NULL end from t_elevation where p && st_buffer($1, 100);
$$;

create or replace function public.st_asgeojson3d(geography) returns text 
LANGUAGE sql IMMUTABLE  AS $$
select st_asgeojson(st_makeline(st_makepoint(st_x(geom), st_y(geom), round(getz(geom))))) from (select (st_dumppoints(geometry($1))).geom) as t;
$$;


tahle funkce vrací hodnotu nadm. výšky v požadovaném bodě pomocí spruměrování 8 okolních pixlu, což je v připadě vrcholových bodu nepřesné.
dá se to nějak vylepšit?
děkuji moc

Nahlásit jako SPAM
IP: 85.216.217.–
Reklama
Reklama
Zjistit počet nových příspěvků

Přidej příspěvek

Toto téma je starší jak čtvrt roku – přidej svůj příspěvek jen tehdy, máš-li k tématu opravdu co říct!

Ano, opravdu chci reagovat → zobrazí formulář pro přidání příspěvku

×Vložení zdrojáku

×Vložení obrázku

Vložit URL obrázku Vybrat obrázek na disku
Vlož URL adresu obrázku:
Klikni a vyber obrázek z počítače:

×Vložení videa

Aktuálně jsou podporována videa ze serverů YouTube, Vimeo a Dailymotion.
×
 
Podporujeme Gravatara.
Zadej URL adresu Avatara (40 x 40 px) nebo emailovou adresu pro použití Gravatara.
Email nikam neukládáme, po získání Gravatara je zahozen.
-
Pravidla pro psaní příspěvků, používej diakritiku. ENTER pro nový odstavec, SHIFT + ENTER pro nový řádek.
Sledovat nové příspěvky (pouze pro přihlášené)
Sleduj vlákno a v případě přidání nového příspěvku o tom budeš vědět mezi prvními.
Reaguješ na příspěvek:

Uživatelé prohlížející si toto vlákno

Uživatelé on-line: 0 registrovaných, 5 hostů

Podobná vlákna

Funkce pro web — založil stoker303

Swing + funkce pro kresleni — založil David

Funkce pro přepínání kódu — založil Thermopolis

 

Hostujeme u Českého hostingu       ISSN 1801-1586       ⇡ Nahoru Webtea.cz logo © 20032016 Programujte.com
Zasadilo a pěstuje Webtea.cz, šéfredaktor Lukáš Churý