Радиус сферы: proj4, WKT

Вопросы по статистическому пакету R. Не обязательно гео.
Ответить
nickleb
Гуру
Сообщения: 847
Зарегистрирован: 22 май 2010, 20:20
Репутация: 93

Радиус сферы: proj4, WKT

Сообщение nickleb » 01 апр 2017, 08:01

Здравствуйте.
Радиус сферы в строке proj4 6371116, а после showWKT 6370997.
Делал так:

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

library(sp)
library(rgdal)

OurProj <- "+proj=laea +lat_0=90 +lon_0=0 +ellps=sphere +a=6371116 +b=6371116 +units=m +no_defs"

print(showWKT(OurProj))

В консоли печатается:
PROJCS[\"Lambert_Azimuthal_Equal_Area\",GEOGCS[\"GCS_Normal Sphere (r=6370997)\",DATUM[\"D_unknown\",SPHEROID[\"sphere\",6370997,0]],PRIMEM[\"Greenwich\",0],UNIT[\"Degree\",0.017453292519943295]],PROJECTION[\"Lambert_Azimuthal_Equal_Area\"],PARAMETER[\"latitude_of_origin\",90],PARAMETER[\"central_meridian\",0],PARAMETER[\"false_easting\",0],PARAMETER[\"false_northing\",0],UNIT[\"Meter\",1]]"
Последний раз редактировалось nickleb 01 апр 2017, 15:39, всего редактировалось 1 раз.

Ariki
Гуру
Сообщения: 730
Зарегистрирован: 12 янв 2011, 22:40
Проекты: 1
Репутация: 303
Ваше звание:

Re: Радиус сферы: proj4, WKT

Сообщение Ariki » 01 апр 2017, 13:29

Нужно убрать +ellps=sphere из строки проекции. Имя "sphere" задаёт конкретную сферу:

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

sphere a=6370997.0      b=6370997.0      Normal Sphere (r=6370997)

nickleb
Гуру
Сообщения: 847
Зарегистрирован: 22 май 2010, 20:20
Репутация: 93

Re: Радиус сферы: proj4, WKT

Сообщение nickleb » 01 апр 2017, 14:24

Ariki писал(а):Нужно убрать +ellps=sphere из строки проекции. Имя "sphere" задаёт конкретную сферу:

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

sphere a=6370997.0      b=6370997.0      Normal Sphere (r=6370997)
Ariki, спасибо Вам большое!

nickleb
Гуру
Сообщения: 847
Зарегистрирован: 22 май 2010, 20:20
Репутация: 93

Re: Радиус сферы: proj4, WKT

Сообщение nickleb » 01 апр 2017, 18:51


Ответить

Вернуться в «R»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 2 гостя