Lagi

Merencanakan peta greenland utara di R


Pengetahuan saya di GIS dan R diuji lagi ketika saya ingin memetakan sebagian kecil Greenland utara (disebut Washington Land) di R (dan data analisis plot di atas nanti).

Saya sudah berhasil merencanakannya dengan

library(maps) library(mapdata) library(rworldmap) map("worldHires", "Greenland", xlim=c(-67,-63), ylim=c(78,82)) newmap <- getMap(resolusi = " tinggi") plot(peta baru, xlim=c(-67, -63), ylim=c(78,82), asp=1)

Tetapi resolusi dan proyeksi sangat buruk untuk area ini dan kebutuhan saya.

Saya memiliki peta yang bagus di MapInfo, tetapi bagaimana saya harus menampilkan peta, sehingga dapat memuat di R dan R dapat memplot data secara akurat di atasnya?


Anda dapat menggunakan data shapefile GADM Level 2 yang tersedia gratis untuk Greenland, yang hadir dalam EPSG:4326 (proyeksi Lintang/Bujur) secara default, dan mengimpornya ke R menggunakan fungsi GDAL yang ditawarkan olehrgdalkemasan. Kemudian, tentukan luas geografis untuk Washington Land dan potong shapefile yang sebelumnya diimpor. Pastikan Anda telah menginstal semua dependensi yang diperlukan (mis.rgeos)

# Paket yang diperlukan library(rgdal) library(raster) # Impor data shapefile Greenland shp.greenland <- readOGR(dsn = "/path/to/shapefile/folder", layer = "GRL_adm2") # Tentukan luas tanam untuk Washington Land ext .washington.land <- extent(-67, -63, 78, 82) # Subset Greenland shapefile shp.washington.land <- crop(shp.greenland, ext.washington.land)

Memperbarui:
Sekarang, jika Anda ingin memasukkan area glasiasi, Anda dapat mengunduh mis. shapefile 10m yang sesuai dari Natural Earth, impor dan potong sesuai dengan tingkat yang telah ditentukan, dan tambahkan ke plot Anda.

# Impor area glaciated shp.glaciers <- readOGR(dsn = "/path/to/shapefile/folder", layer = "ne_10m_glaciated_areas") # Subset (dan plot) area glaciated menurut luas Washington Land shp.glaciers <- crop( shp.glaciers, ext.washington.land) plot(shp.washington.land, col = "darkolivegreen3") plot(shp.glaciers, col = "azure", add = TRUE)


Tonton videonya: Peta Konsep Benua Eropa Utara (September 2021).