-
Notifications
You must be signed in to change notification settings - Fork 125
/
sql_mapper.go
46 lines (38 loc) · 1.47 KB
/
sql_mapper.go
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
package geo
import (
"database/sql"
"fmt"
)
// A Mapper that uses Standard SQL Syntax to perform mapping functions and queries
type SQLMapper struct {
conf *SQLConf
sqlConn *sql.DB
}
// NewSQLMapper: Creates and returns a pointer to a new geo.SQLMapper.
func NewSQLMapper(filename string, conn *sql.DB) (*SQLMapper, error) {
conf, confErr := GetSQLConfFromFile(filename)
if confErr != nil {
return nil, confErr
}
return &SQLMapper{conf: conf, sqlConn: conn}, nil
}
// SqlDbConn returns a pointer to the SQLMapper's SQL Database Connection.
func (s *SQLMapper) SqlDbConn() *sql.DB {
return s.sqlConn
}
// PointsWithinRadius: Uses SQL to retrieve all points within the radius (in meters)
// passed in from the origin point passed in.
// Original implemenation from : http://www.movable-type.co.uk/scripts/latlong-db.html
// Returns a pointer to a sql.Rows as a result, or an error if one occurs during the query.
func (s *SQLMapper) PointsWithinRadius(p *Point, radius float64) (*sql.Rows, error) {
select_str := fmt.Sprintf("SELECT * FROM %v a", s.conf.table)
lat1 := fmt.Sprintf("sin(radians(%f)) * sin(radians(a.lat))", p.lat)
lng1 := fmt.Sprintf("cos(radians(%f)) * cos(radians(a.lat)) * cos(radians(a.lng) - radians(%f))", p.lat, p.lng)
where_str := fmt.Sprintf("WHERE acos(%s + %s) * %f <= %f", lat1, lng1, float64(EARTH_RADIUS), radius)
query := fmt.Sprintf("%s %s", select_str, where_str)
res, err := s.sqlConn.Query(query)
if err != nil {
panic(err)
}
return res, err
}