R mapproj lambert-азимутална проекция с еднаква площ

Как мога да направя ламбертова азимутална проекция с еднаква площ с пакета mapproj? http://cran.r-project.org/web/packages/mapproj/mapproj.pdf

Има

azequalarea() равна площ

и

lambert(lat0,lat1) конформен, истински мащаб на lat0 и lat1

Коя е най-близка до азимуталната проекция с еднаква площ на Ламберт?


person flohump    schedule 04.10.2012    source източник


Отговори (1)


В пакет mapproj и функцията mapproject() изглежда, че азимуталната проекция с еднаква площ на Ламберт ще бъде azequalarea(), тъй като lambert(lat0,lat1) е равномерната конична проекция на Ламберт (тъй като е посочена в ръководството в Polar conic projections symmetric about the Prime Meridian).

В пакет GEOmap проекцията се извършва с помощта на функция lamaz.eqarea().

В пакет rgdal това се прави с извикване на функция spTransform(), съдържаща низ proj4, започващ с +proj=laea.

person plannapus    schedule 04.10.2012
comment
Хубав отговор. На всичко друго освен много малък проект, силно бих препоръчал да се придържате към proj4strings за последователност и гъвкавост. - person Ari B. Friedman; 04.10.2012
comment
Благодаря за вашият отговор! Тъй като искам да начертая карти с ggplot2, разчитам на mapproject(). - person flohump; 06.10.2012