P.S.Pudyastuti's Lovely Rooms…

Learning is a long-life process, so it's never late to learn…

Bermain dengan Maptools

Leave a comment

Kali ini saya ingin mengajak bermain R dengan menggunakan package Maptools, untuk menggambar peta. Package Maptools bisa diunduh dari laman CRAN. Untuk menggunakan package Maptools, file input data yang bisa digunakan antara lain adalah file shp (dari ArcMap) dan file csv.

Berikut ini contoh script sederhana menggambar peta lokasi stasiun pengukur curah hujan di suatu sub DAS (Daerah Aliran Sungai) dengan menggunakan package Maptools.

setwd(“J:/Prog/”)
library(maptools)

#membaca shp file peta sub DAS Hilir
lower <- readShapeSpatial(“hilir”)

#membaca shp file peta sungai
rivers <- readShapeLines(“sungai1”)

#membaca csv file yang berisi lokasi stasiun pengukur curah hujan
sta1 <- read.csv(“rainsta-1.csv”,header=TRUE)

#menggambar peta sub DAS Hilir
plot (lower,col=”wheat1″,border=”brown”,lwd=1.5,axes=TRUE,
xlim=c(500000,700000),ylim=c(9150000,9250000),cex.axis=0.8)

#menambahkan gambar sungai ke dalam peta
lines(rivers, col=”dodgerblue3″, lwd=0.8)

#menambahkan lokasi stasiun pengukur curah hujan ke dalam peta
coordinates(sta1) <- c(“X”,”Y”,”Elev”)
points(coordinates(sta1),pch=20,col=”red”)

#menambahkan judul peta
title(main=”Sub DAS Hilir”,cex.main=0.9)
title(xlab=”X”,ylab=’Y’,cex.lab=1.0,line=2.25)

#menambahkan tanda arah Utara
SpatialPolygonsRescale(layout.north.arrow(1), offset= c(500000,9270000), scale = 15000,
col=”grey”,plot.grid=F)

#menambahkan Legend pada peta
legend(“bottomright”, title=”Legend”,legend=”Stasiun Pengukur Hujan”,cex=0.8,pch=20,
col=”red”, bg=”white”)

Jika script di atas di-run, maka hasil outputnya adalah seperti gambar di bawah ini.

Screenshot running R.

Screenshot running R.

Jika ingin menyimpan hasil running dalam file tersendiri, letakkan kursor pada peta hasil running, dan klik tab File pada jendela RGui, lalu pilih format file yang diinginkan. Pilihan format file antara lain Metalife, Postcript, PDF, Png, Bmg, TIFF, dan Jpeg.

Gambar di bawah ini adalah output hasil running script di atas yang saya simpan dalam format Jpeg.

Peta hasil running script di atas.

Peta hasil running script di atas.

Selamat mencoba.

Salam,

PSP

Author: P.S.Pudyastuti

I am an academic staff at Civil Engineering Department, Faculty of Engineering, Universitas Muhammadiyah Surakarta - Indonesia. I got my Bachelor Degree (ST) from Gadjah Mada University, Indonesia, and Master Degree (M.Sc) from Loughborough University, England. I am married and have got 3 children.

Leave a Reply

Please log in using one of these methods to post your comment:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s