Gradovi (Srbija) SQL Ready

[webarto]

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

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

http://pokit.etf.ba/get/d653d36ef96bec9a0f9f06dfa2affc9c.png

[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.

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;

}

$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.


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)...

http://pokit.etf.ba/get/4326b65f86d515365d0636d0cf3920d2.png

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?

[webarto]

//Beograd
$grad = mysql_query("SELECT * FROM srbija WHERE grad = 'Beograd' LIMIT 1");
$grad = mysql_fetch_assoc($grad);
$lat = $grad["lat"]; $lon = $grad["lon"];

$radius = 10; // 10 kilometara
$radius = round($radius / 111, 7); // 1 stepen = 60 nautičkih milja ~ 111km, zemlja nije krug ali je razlika 0.3%

$sql = mysql_query("SELECT * FROM srbija WHERE lat BETWEEN ($lat - $radius) AND ($lat + $radius) AND lon BETWEEN ($lon - $radius) AND ($lon + $radius)");
while($red = mysql_fetch_assoc($sql)){
echo $grad["grad"]." - ".$red["grad"]." = ".haversine($lat, $lon, $red["lat"], $red["lon"])."<br />";
}

Beograd - Alt-Borscha = 5.79316506111
Beograd - Bárányos = 9.35219525075
Beograd - Belgrad = 0
Beograd - Belgrade = 0
Beograd - Belgrado = 0
Beograd - Beograd = 0
Beograd - Bezanija = 7.20118432749
Beograd - Borca = 5.79316506111
Beograd - Borcsa = 5.79316506111
Beograd - Bubanj Potok = 11.085803226
Beograd - Cukarica = 5.50759465291
Beograd - Dorcol = 0.242342119582
Beograd - Filmski Grad = 7.29055150674
Beograd - Gisellenhain = 9.35219525075
Beograd - Jajince = 8.8637130281
Beograd - Jajinci = 8.8637130281
Beograd - Kaluderica = 10.0123750758
Beograd - Kaludirica = 10.0123750758
Beograd - Kalugjerica = 10.0123750758
Beograd - Karaburma = 2.66566666921
Beograd - Kumodraz = 9.39940017819
Beograd - Mala Ciganti = 4.38736003825
Beograd - Mali Mokri Lug = 7.03790686951
Beograd - Mokri Lug = 8.23222134719
Beograd - Novi Beograd = 3.75416324228
Beograd - Ovca = 9.35219525075
Beograd - Ovcsa = 9.35219525075
Beograd - Reva = 5.46912292882
Beograd - Rospi Cuprija = 4.62338905561
Beograd - Sajmiste = 1.93995055969
Beograd - Savski Venac = 4.52617380096
Beograd - Semlin = 5.93785726157
Beograd - Singidunum = 0
Beograd - Stara Borca = 5.79316506111
Beograd - Stari Grad = 1.68389610582
Beograd - Taurunum = 5.93785726157
Beograd - Tosin Bunar = 5.03993677578
Beograd - Veliki Mokri Lug = 8.23222134719
Beograd - Vojna Basta = 9.27635855118
Beograd - Vozdivac = 4.52055022048
Beograd - Vozdovac = 4.52055022048
Beograd - Vracar = 3.96300110998
Beograd - Zarkovo = 7.83593874715
Beograd - Zeleznicka Kolonja = 6.9968568556
Beograd - Zemlén = 5.93785726157
Beograd - Zemun = 5.93785726157
Beograd - Zimony = 5.93785726157
Beograd - Zvezdara = 5.76970836174

Novi Sad

Novi Sad - Kamanc = 3.18673675777
Novi Sad - Kamancz = 3.18673675777
Novi Sad - Kamenicz = 3.18673675777
Novi Sad - Kamenitz = 3.18673675777
Novi Sad - Kamonc = 3.18673675777
Novi Sad - Neoplanta = 0
Novi Sad - Neusatz = 0
Novi Sad - Novi Ledinci = 5.23735837334
Novi Sad - Novi Rakovac = 7.48661578121
Novi Sad - Novi Sad = 0
Novi Sad - Nový Sad = 0
Novi Sad - Ó-péterváradja = 0
Novi Sad - Pasuljiste = 4.70016657699
Novi Sad - Pétervárad = 3.37320414099
Novi Sad - Peterwardein = 3.37320414099
Novi Sad - Petrovaradin = 3.37320414099
Novi Sad - Sremska Kamenica = 3.18673675777
Novi Sad - Újvidék = 0
Novi Sad - Varadinum Petri = 0
Novi Sad - Vásáros-Várad = 0
Novi Sad - Veternik = 6.06948254759

Izvinjavam se na dužem postu ;)

[uros]

preciscena verzija... za slucaj da jos neko ne zeli Rospi Cupriju, Tosin Bunar etc...

gradove sam pokupio sa belih strana...



CREATE TABLE `gradovi` (
`id` int(11) not null auto_increment,
`grad` varchar(100) not null,
`lat` varchar(10) not null,
`lon` varchar(10) not null,
PRIMARY KEY (`id`),
KEY `grad` (`grad`)
) ENGINE=MyISAM DEFAULT CHARSET=utf8 AUTO_INCREMENT=34;

INSERT INTO `gradovi` (`id`, `grad`, `lat`, `lon`) VALUES
('1', 'Beograd', '44.818611', '20.468056'),
('2', 'Bor', '44.078333', '22.095278'),
('3', 'Cacak', '43.891389', '20.349722'),
('4', 'Gnjilane', '42.468889', '21.463333'),
('5', 'Jagodina', '43.981389', '21.262222'),
('6', 'Kikinda', '45.641111', '20.414167'),
('7', 'Kosovska Mitrovica', '42.883333', '20.866667'),
('8', 'Kragujevac', '44.016667', '20.916667'),
('9', 'Kraljevo', '43.725833', '20.689444'),
('10', 'Krusevac', '43.58', '21.333889'),
('11', 'Leskovac', '42.998056', '21.946111'),
('12', 'Nis', '43.324722', '21.903333'),
('13', 'Novi Pazar', '43.136667', '20.512222'),
('14', 'Novi Sad', '45.251667', '19.836944'),
('15', 'Pancevo', '44.870833', '20.640278'),
('16', 'Pec', '42.66', '20.292222'),
('17', 'Pirot', '43.153056', '22.586111'),
('18', 'Pozarevac', '44.615278', '21.1825'),
('19', 'Prijepolje', '43.371944', '19.640556'),
('20', 'Pristina', '42.666667', '21.166667'),
('21', 'Prizren', '42.213889', '20.739722'),
('22', 'Prokuplje', '43.234167', '21.588056'),
('23', 'Sabac', '44.746667', '19.69'),
('24', 'Smederevo', '44.662778', '20.93'),
('25', 'Sombor', '45.774167', '19.112222'),
('26', 'Sremska Mitrovica', '44.976389', '19.612222'),
('27', 'Subotica', '46.1', '19.666667'),
('28', 'Urosevac', '42.370556', '21.155278'),
('29', 'Uzice', '43.855833', '19.841111'),
('30', 'Valjevo', '44.270833', '19.884167'),
('31', 'Vranje', '42.551389', '21.900278'),
('32', 'Zajecar', '43.904167', '22.284722'),
('33', 'Zrenjanin', '45.383611', '20.381944');

[webarto]

Sa Wikipedia ali bez LL :)

INSERT INTO gradovi(grad) VALUES('Beograd');
INSERT INTO gradovi(grad) VALUES('Niš');
INSERT INTO gradovi(grad) VALUES('Kragujevac');
INSERT INTO gradovi(grad) VALUES('Čačak');
INSERT INTO gradovi(grad) VALUES('Kraljevo');
INSERT INTO gradovi(grad) VALUES('Šabac');
INSERT INTO gradovi(grad) VALUES('Smederevo');
INSERT INTO gradovi(grad) VALUES('Valjevo');
INSERT INTO gradovi(grad) VALUES('Kruševac');
INSERT INTO gradovi(grad) VALUES('Zaječar');
INSERT INTO gradovi(grad) VALUES('Užice');
INSERT INTO gradovi(grad) VALUES('Vranje');
INSERT INTO gradovi(grad) VALUES('Novi Pazar');
INSERT INTO gradovi(grad) VALUES('Požarevac');
INSERT INTO gradovi(grad) VALUES('Pirot');
INSERT INTO gradovi(grad) VALUES('Bor');
INSERT INTO gradovi(grad) VALUES('Jagodina');
INSERT INTO gradovi(grad) VALUES('Prokuplje');
INSERT INTO gradovi(grad) VALUES('Paraćin');
INSERT INTO gradovi(grad) VALUES('Smederevska Palanka');
INSERT INTO gradovi(grad) VALUES('Aranđelovac');
INSERT INTO gradovi(grad) VALUES('Gornji Milanovac');
INSERT INTO gradovi(grad) VALUES('Lazarevac');
INSERT INTO gradovi(grad) VALUES('Obrenovac');
INSERT INTO gradovi(grad) VALUES('Mladenovac');
INSERT INTO gradovi(grad) VALUES('Loznica');
INSERT INTO gradovi(grad) VALUES('Ćuprija');
INSERT INTO gradovi(grad) VALUES('Priboj');
INSERT INTO gradovi(grad) VALUES('Novi Sad');
INSERT INTO gradovi(grad) VALUES('Subotica');
INSERT INTO gradovi(grad) VALUES('Zrenjanin');
INSERT INTO gradovi(grad) VALUES('Pančevo');
INSERT INTO gradovi(grad) VALUES('Sombor');
INSERT INTO gradovi(grad) VALUES('Kikinda');
INSERT INTO gradovi(grad) VALUES('Sremska Mitrovica');
INSERT INTO gradovi(grad) VALUES('Vršac');
INSERT INTO gradovi(grad) VALUES('Ruma');
INSERT INTO gradovi(grad) VALUES('Bačka Palanka');
INSERT INTO gradovi(grad) VALUES('Inđija');
INSERT INTO gradovi(grad) VALUES('Vrbas');
INSERT INTO gradovi(grad) VALUES('Bečej');
INSERT INTO gradovi(grad) VALUES('Senta');
INSERT INTO gradovi(grad) VALUES('Kula');
INSERT INTO gradovi(grad) VALUES('Apatin');
INSERT INTO gradovi(grad) VALUES('Temerin');
INSERT INTO gradovi(grad) VALUES('Priština');
INSERT INTO gradovi(grad) VALUES('Prizren');
INSERT INTO gradovi(grad) VALUES('Peć');
INSERT INTO gradovi(grad) VALUES('Đakovica');
INSERT INTO gradovi(grad) VALUES('Kosovska Mitrovica');
INSERT INTO gradovi(grad) VALUES('Gnjilane');
INSERT INTO gradovi(grad) VALUES('Podujevo');
INSERT INTO gradovi(grad) VALUES('Uroševac');
INSERT INTO gradovi(grad) VALUES('Kosovo Polje');
INSERT INTO gradovi(grad) VALUES('Orahovac');

[jablan]

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?

Hmm, ne. Tvoj query gleda praktično pripadnost kvadratu opisanom oko kružnice radijusa R. To praktično znači da ako tražiš radijus od 100 km oko Beograda, u njega će upasti i tačka koja je 140km od beograda (u pravcu dijagonale). Nema veze.

[vidak]

evo linka (http://websoft.me/images/stories/exYU.tar) sa oko 1600 mesta u Republici Srbiji.
Osim Srbije u bazi su države exYU i njihovi gradovi i mesta.

[ivanhoe]

za ovakve stvari postoji podrska za spatial indekse u mysqlu, to je mnogo brze od ove DIY trigonometrije
http://dev.mysql.com/doc/refman/5.0/en/creating-a-spatially-enabled-mysql-database.html

[webarto]

http://pokit.etf.ba/get/47707ce31359897bf9a58dc43edfc506.png

U pravu si jablane, ide od -1 do 1 :)

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

$sql = mysql_query("
SELECT g1 . *
FROM gradovi AS g1
INNER JOIN gradovi AS g2 ON POW( (
g1.lat - g2.lat
), 2 ) + POW( (
g1.lon - g2.lon
), 2 ) < POW( ( $radius / 111.12 ) , 2 )
WHERE g2.grad = '$grad'");

return mysql_fetch_assoc($sql);

}

[webarto]

class Radius{

public $table = "routes";
public $column_lat = "lat";
public $column_lon = "lon";
public $order = "distance";
public $order_direction = "ASC";
public $kilometers = true;

public function build_query($lat, $lon, $distance){

if($kilometers){
$multiplier = 112.12;
}else{
$multiplier = 69.0467669;
}

$query =
"SELECT *, (SQRT(POW(($this->column_lat - $lat), 2) + POW(($this->column_lon - $lon), 2)) * $multiplier) AS distance
FROM $this->table
WHERE POW(($this->column_lat - $lat), 2) + POW(($this->column_lon - $lon), 2) < POW(($distance / $multiplier), 2)
ORDER BY $this->order $this->order_direction";

return $query;

}

}

$lat = 44.818611;
$lon = 20.468056;

$radius = new Radius;
$query = $radius->build_query($lat, $lon, 10);
/**
SELECT *, (SQRT(POW((lat - 44.818611), 2) + POW((lon - 20.468056), 2)) * 112.12) AS distance
FROM routes
WHERE POW((lat - 44.818611), 2) + POW((lon - 20.468056), 2) < POW((10 / 112.12), 2)
ORDER BY distance ASC */

$result = mysql_query($query);
while($row = mysql_fetch_assoc($result)){

}

http://i.imgur.com/MjCYU.png
(Picture unrelated)

I tako to :)