Форумы GIS-Lab.info

Vavik писал(а): Нашёл вот такого вида запросик, может посмотрев на него поймёте чего я хочу) правда он не работает(

SELECT ST_Distance(a.location, b.location)
FROM points_tbl a, points_tbl b
WHERE a.id = 1000 AND b.id = 3029;

SELECT ST_DISTANCE(a.the_geom,b.the_geom) FROM city a, city b WHERE a.city_id=333 AND b.city_id=3444;

ещё, какой запрос выведет мне координаты точки, если я укажу её название или id ?

SELECT ST_X(the_geom),ST_Y(the_geom) FROM city WHERE city_id=333;

SELECT AsText(the_geom) FROM city WHERE city_id=333;

Vavik писал(а):Спасибо, всё работает:) Правда у меня расстояния неизвестно в чём указываются. Получил вот такое число 0,01887057232, а реальное расстояние около 1.75-1.95 км. В каких велечинах я получил результат ?) SRID 26986, ранее был 4326, но ничего не изменилось.

Ещё нужен запрос, который выдавал все точки, находящиеся в радиусе другой. Делаю так

Код: Выделить всё

 SELECT  b.name,st_dwithin(a.the_geom,b.the_geom, 1000)  from table_p a, table_p b 
WHERE a.name~'point' AND b.name~'mod point' 

слово "mod" присутствует в нескольких точках, то есть мне надо, чтобы был результат с такими точками попавшими в радиус, а он мне выводит, все которые только есть.

Нашёл ещё такой запрос,который должен это делать

Код: Выделить всё

SELECT * FROM geotable 
	WHERE ST_DWithin(geocolumn, 'POINT(1000 1000)', 100.0);

Но не пойму что такое geocolumn ?

KolesovDmitry писал(а): Не видя данные, и не зная, откуда они взялись, как вы их положили в PostGIS, сложно что-то сказать... Но очень подозрительно то, что раньше SRID был 4326, а теперь 26986, но результат не изменился (поскольку у первой системы единицы измерения градусы, а у второй - метры). Вы как данные трансформировали из одной системы в другую?

KolesovDmitry писал(а): Честно сказать, я не понял, что нужно сделать...

SELECT  b.name,st_dwithin(a.the_geom,b.the_geom, 1000)  from table_p a, table_p b 
WHERE a.name~'point' AND b.name~'tip1 point'

Vavik писал(а): Вот так я импортировал данные в Postgis

shp2pgsql -i -D -s 4326 PolyanaKrasnaya_A.shp PolyanaKrasnaya_A > PolyanaKrasnaya_A.sql
psql -U postgres -f PolyanaKrasnaya_A.sql -d PolyanaKrasnaya(map)

Тут мне не понравились результаты запросов с непонятными цифрами, полазил на одном форуме, там советывали перевести SRID на 26986. Ну я по новой загнал Данные только поменяв SRID, а толку 0.

shp2pgsql -i -D -s 26986 PolyanaKrasnaya_A.shp PolyanaKrasnaya_A > PolyanaKrasnaya_A.sql
psql -U postgres -f PolyanaKrasnaya_A.sql -d PolyanaKrasnaya(map)

Вообщем есть такая функция
ST_DWithin(geometry, geometry, float)
Которая возвращает истину, если геометрии находятся в пределах указанного расстояния одна от другой.

У меня же имеются точки - точки tip1 точки tip2 и tip3. Я хочу в результате получить, например точки tip1 попадающие в радиус отдельно взятой точки. Только, мне нужно, чтобы выводились название этих точек, а не лишь логическое значение, то что они попадают в радиус.
Я сделал такой запрос

Код: Выделить всё

SELECT  b.name,st_dwithin(a.the_geom,b.the_geom, 1000)  from table_p a, table_p b 
WHERE a.name~'point' AND b.name~'tip1 point'

Но он мне выводит все точки tip1, которые имеются, в не зависимости, попадают они в радиус или нет.

SELECT b.name,a.name FROM table_p a, table_p b WHERE ST_DWithin(a.geometry, b.geometry, 1000) AND b.name='чему он там равняется'

Vavik писал(а):Вроде SRID 4326, а как его поменять на 26986 ?

• 1) преобразовать ваш shp перед экспортом в SQL. Это можно сделать, например, при помощи утилиты ogr2ogr или любым другим способом.
  2) экспрортировать shp в PostGIS со SRID 4326, а потом уже в PostGIS (его средствами) преобразовать в требуемую систему. Для этого воспользуйтесь функцией ST_TRANSFORM (функция очень простая, подробности об использовании смотрите в документации).

ProFF писал(а):Помогите пишу запрос в Postgis - е

shp2pgsql -i -D -s 4326 wayline.shp wayline > wayline.sql

Выдает такую ошибку:
ERROR: ошибка синтаксиса в или рядом "shp2pgsql"
СТРОКА 1:shp2pgsql -i -D -s 4326 wayline.shp wayline > wayline.sql
^

********** Ошибка **********

ERROR: ошибка синтаксиса в или рядом "shp2pgsql"
SQL state: 42601
Характеристика:1

Что делать???

KolesovDmitry писал(а): Тогда есть два пути:
1) преобразовать ваш shp перед экспортом в SQL. Это можно сделать, например, при помощи утилиты ogr2ogr или любым другим способом.
...

KolesovDmitry писал(а): Тогда есть два пути:
1) преобразовать ваш shp перед экспортом в SQL. Это можно сделать, например, при помощи утилиты ogr2ogr или любым другим способом.
...

ProFF писал(а):Подскажите как изменить кодировку у базы с win1251 на UTF8

Vavik писал(а):

ProFF писал(а):Подскажите как изменить кодировку у базы с win1251 на UTF8

Кажись никак. Можно сделать копию базы, потом создать новую в нужной кодировке и восстановить туда сохранённую бд) А вообще, чтто меняет кодировка ? а то я не знаю для чего она))