Gradovi (Srbija) SQL Ready

[webarto]

Izvučeno iz MaxMind baze podataka, i pripremljeno za SQL.

http://static.webarto.com/download/maxmind/serbia.zip

[lurker]

ehm

VALUES('rs','Bacsszollos','46.133333','19.866667') ;
INSERT INTO gradovi(drzava,grad,lat,lon) VALUES('rs','Bacsszolos','46.133333','19.866667');
INSERT INTO gradovi(drzava,grad,lat,lon) VALUES('rs','Bactovaros','45.353611','19.325');
.
.
.
pa Baja i jos mnogo gradova u Madjarskoj
mnogo smo velika drzava ili se sprema neki novi rat?

wtf?

[webarto]

"Backi Vinogradi is also known as Bacsszollos"... ne kontam? To je izvučeno iz MaxMind baze podataka a za te ratove ne znam

[ivanhoe]

jel se krsi neki copyright ovim? samo pitam, inace hvala na trudu u svakom slucaju

[webarto]

Redistribution and use with or without modification, are permitted provided
that the following conditions are met:
1. Redistributions must retain the above copyright notice, this list of
conditions and the following disclaimer in the documentation and/or other
materials provided with the distribution.
2. All advertising materials and documentation mentioning features or use of
this database must display the following acknowledgement:
"This product includes data created by MaxMind, available from
http://www.maxmind.com/"
3. "MaxMind" may not be used to endorse or promote products derived from this
database without specific prior written permission.

Napisao sam gore da je izvučeno iz MaxMind, valjda je dovoljno

Nije me još ViaMichelin zvao
http://www.google.com/search?q=viamichelin+php

[webarto]

Možda nekom bude korisno, kako odabrati sva naselja u krugu od nekog grada npr.

PHP kôd:

function radius($grad, $radius = "10"){

    $radius = round($radius / 111, 7);

        $sql = mysql_query("SELECT lat,lon FROM gradovi WHERE LOWER(grad) = '$grad' LIMIT 1");
        $grad = mysql_fetch_assoc($sql);
        $lat = $grad["lat"]; $lon = $grad["lon"];
    
        $sql = mysql_query("SELECT grad FROM gradovi WHERE lat BETWEEN ($lat - $radius) AND ($lat + $radius) AND lon BETWEEN ($lon - $radius) AND ($lon + $radius)");
        while($red = mysql_fetch_assoc($sql)){
            $gradovi[] = $red["grad"];
        }

    return $gradovi;

} 


PHP kôd:

$gradovi = radius("banja luka");
$gradovi = join(", ", $gradovi);
echo $gradovi; 


Anusici, Banialuca, Banja Luka, Banja Luka-Vrbanja, Baralici, Barica, Bastahi, Bukvalek, Bukvaluk, Cesma, Cokori, Curlici, Cvisici, Debeljaci, Delibasino Selo, Dikevci, Diljevici, Donji Ducipolje, Donji Ponir, Dragicevici, Dragocaj, Drakulic, Gavranici, Glamocani, Gornje Ducipolje, Gornje Presnace, Gornji Ponir, Gornji Seher, Grabljani, Grijecani, Hiseti, Jagare, Jajcevici, Jakobasici, Joldici, Kajkuti, Kola Donja, Kolibiste, Kolonija, Koprene, Krcma, Krcmarice, Krecari, Krndija, Kuljani, Kumala, Laus, Ljevari, Madir, Magljani, Makivici, Mali Prnjavor, Matosevci, Meljani, Moconji, Motike, Nazaret, Novi Kovici, Novoselija, Orlovac, Ostrike, Petricevac, Peulje, Pistelici, Ponir, Posalici, Potkucnice, Prijecani, Prnjavor Mali, Pustahije, Rebrovac, Saracica, Sargovac, Sedici, Spahinjci, Srpske Toplice, Starcevica, Stranjani, Sumonje, Trapisti, Trn, Udovicici, Visekruna, Vrbanja, Vuklisevici, Zaluzani, Zisci

Nažalost nisam importovao bazu Srbije, ali je svejedno npr, samo se unese "beograd", i izlistaće se sva naselja, iz baze naravno u krugu od onoliko kilometara koliko navedete.

[jablan]

^ Ako se ne varam, nije u pitanju krug oko nekog mesta, već kvadrat.

[jablan]

Pod pretpostavkom da se podelom sa 111 dovoljno dobro aproksimira konverzija kilometara u stepene:

select m1.* from mesta m1
INNER JOIN mesta m2 ON (m1.lat - m2.lat) ^ 2 + (m1.lon - m2.lon) ^ 2 < (10.0/111.0) ^ 2
WHERE m2.mesto = 'Krusevac';

Ovo je za Postgres, možda treba malo promeniti da bi radilo u MySQL.

Takođe, lat i lon valja importovati kao NUMERIC, ne VARCHAR.

[zira]

Za precizniji racun treba malo vise matematike, ima detaljno objasnjeno za slucaj MySQL-a na http://www.arubin.org/files/geo_search.pdf

Na primjer, za nalazenje 10 najbilizih "tacaka" zadatoj koordinati:

set @orig_lat=121.9763; set @orig_lon=37.40445;
set @dist=10;

SELECT *, 3956 * 2 * ASIN(SQRT(
POWER(SIN((@orig_lat - abs(dest.lat)) * pi()/180 / 2),
2) + COS(@orig_lat * pi()/180 ) * COS(abs(dest.lat) *
pi()/180) * POWER(SIN((@orig_lon - dest.lon) *
pi()/180 / 2), 2) )) as distance
FROM hotels dest
having distance < @dist
ORDER BY distance limit 10


(obrati paznju da je dist u miljama, za kilometre konvertovati 3956 u kilometre)

[webarto]

^ Haversine (half sinus versus) funkcija, tj razdaljina između dvije tačke u koordinatom sistemu sa zakrivljenošću.

PHP kôd:

function haversine($phi1, $lambda1, $phi2, $lambda2){
    $radius = 6371;
    
    $dPhi = deg2rad($phi2 - $phi1);
    $dLambda = deg2rad($lambda2 - $lambda1);
    
    $a = sin($dPhi/2) * sin($dPhi/2) + cos(deg2rad($phi1)) * cos(deg2rad($phi2)) * sin($dLambda/2) * sin($dLambda/2);
    $c = 2 * asin(sqrt($a));
    $d = $radius * $c;
    
    return $d;
} 


@jablan, kontam šta hoćeš reći mada mislim da nije (probao sam iscrtavajući radius na mapi)...



Zelene su iste, kvadrat se može posmatrati kao 2 jednakostranična trougla i vidi da hipotenuza nije ista kao dužine stranica (normalno).
Jesi li na to mislio?