Pesatnya perkembangan suatu wilayah, maka bertambah juga penduduk yang ada, sehingga berdampak terhadap meningkatnya kebutuhan akan tempat tinggal. Perubahan lahan merupakan hal yang umum terjadi, baik dalam bidang perencanaan wilayah dan kota, maupun pengelolaan sumber daya alam. Perubahan lahan dapat diartikan sebagai praktik konversi lahan menjadi lahan perkebunan atau peternakan, perluasan fungsi lahan pertanian, penggundulan hutan, penanaman kembali fungsi lahan hutan, dan ekspansi lahan perkotaan [urban sprawl]. Salah satu dampak dari perubahan lahan yang tidak teratur adalah berkurangnya fungsi produktivitas biologis dan keberagaman ekosistem. Berdasarkan hal tersebut, isu terkait perubahan lahan masih relevan untuk dianalisis, dalam konteks identifikasi perubahan tutupan lahan [land cover] dari periode waktu tertentu.
KLASIFIKASI TUTUPAN LAHAN
Definisi tutupan lahan mengacu dari SNI 7645:2010 tentang Klasifikasi Penutup Lahan adalah tutupan biofisik pada permukaan bumi yang dapat diamati merupakan suatu hasil pengaturan, aktivitas, dan perlakuan manusia yang dilakukan pada jenis penutup lahan tertentu untuk melakukan kegiatan produksi, perubahan, ataupun perawatan pada penutup lahan tersebut. Kelas penutup lahan dibagi menjadi dua bagian besar, yaitu daerah bervegetasi dan daerah tak bervegetasi. Daerah bervegetasi merupakan daerah dengan liputan vegetasi minimal 4% setidaknya selama 2 bulan, atau dengan liputan Lichens/ Mosses lebih dari 25% jika tidak terdapat vegetasi lain. Terdapat kelas daerah pertanian dan daerah bukan pertanian. Daerah tak bervegetasi merupakan daerah dengan total liputan vegetasi kurang dari 4% selama lebih dari 10 bulan, atau daerah dengan liputan Lichens/ Mosses kurang dari 25% jika tidak terdapat vegetasi kayu atau herba. Terdapat kelas lahan terbuka, permukiman dan lahan bukan pertanian yang berkaitan, dan perairan. Penjelasan lebih lengkap dapat diakses melalui laman berikut SNI 7645:2010 tentang Klasifikasi Penutup Lahan.
Tutorial ini berusaha untuk memberikan demonstrasi secara sederhana terkait pemanfaatan Sistem Informasi Geografis [SIG] dengan menggunakan data penginderaan jauh berupa citra Landsat. Dalam rangka mengidentifikasi perubahan lahan di Kawasan Metropolitan Malang Raya [Kota Malang, Kota Batu, dan Kabupaten Malang] selama dua periode yaitu tahun 2013 dan tahun 2019. Luas wilayah Kawasan Metropolitan Malang Raya yang mencakup 2 kota dan 1 kabupaten adalah sebesar 3.769 km2.
Daerah perkotaan dipilih berdasarkan perubahan yang drastis akibat pengaruh urbanisasi dan pertumbuhan penduduk yang memberikan pengaruh signifikan terhadap kondisi lingkungan, ekonomi, dan kehidupan sosial perkotaan. Kawasan Metropolitan Malang Raya merupakan salah satu daerah perkotaan yang terdapat di Indonesia yang mengalami perkembangan perkotaan yang cukup pesat. Tingkat urbanisasi yang tinggi dan adanya proses menjadi kota, khususnya di daerah lingkar perkotaan telah memberikan dampak yang besar terhadap perubahan di Kawasan Metropolitan Malang Raya. Informasi terbaru dan akurat tentang kondisi dan kecenderungan perubahan wilayah berupa identifikasi perubahan lahan, dapat dimanfaatkan untuk mengembangkan beberapa strategi pembangunan berkelanjutan dan peningkatan kehidupan wilayah perkotaan.
PENGGUNAAN TEKNOLOGI SIG [SISTEM INFORMASI GEOGRAFIS]
Secara sederhana, sangat mudah untuk melakukan identifikasi tutupan lahan. Metode yang bisa dilakukan adalah digitasi manual on screen berdasarkan kelas tutupan lahan yang terlihat degan citra satelit, kemudian mengkoreksi dengan kondisi eksisting di lapangan. Tetapi, hal ini bisa menjadi kompleks dan memakan banyak waktu apabila daerah yang menjadi area of interest [AOI] memiliki luas yang besar. Untuk itu diperlukan metode yang lebih efektif dan efisien, dengan menggunakan salah satu analisis [toolbox] pada perangkat lunak ArcMap. ArcMap merupakan salah satu produk perangkat lunak SIG berbasis desktop yang diproduksi oleh ESRI. ArcMap memiliki kemampuan utama untuk visualisasi, membangun database spasial yang baru, memilih data [query], editing, menciptakan desain-desain peta, analisis dan pembuatan tampilan akhir dalam laporan-laporan kegiatan.
Tutorial ini menggunakan ArcGIS Desktop yaitu ArcMap 10.8, perangkat lunak tersebut dapat diperoleh melalui laman berikut ArcGIS Desktop - ArcMap. Metode pengolahan data citra landsat dalam ArcGIS menggunakan metode klasifikasi tersupervisi [Supervised Classification]. Metode ini mengkelaskan citra berdasarkan pengenalan spektral [nilai reflektan] yang didapatkan dari sample piksel [poligon yang merepresentasikan sampel area untuk setiap jenis tutupan lahan yang berbeda]. Sampel ini dikoleksi secara manual dan analisa citra digunakan untuk menghitung klasifikasi citra. Metode ini membutuhkan ekstensi Spatial Analyst dan tools [ArcToolbox] antara lain Composite, Clip, Pan-sharpened, Maximum likelihood classification.
Data wilayah yang digunakan adalah shapefile batas wilayah dari Kawasan Metropolitan Malang Raya yang terdiri dari Kota Malang, Kota Batu, dan Kabupaten Malang. Sumber data batas wilayah administrasi dari portal resmi Badan Informasi Geospasial [BIG] yang dapat diakses oleh umum melalui laman berikut Indonesia Geospatial Portal.
Data citra satelit yang digunakan adalah citra Landsat, Landsat merupakan program penangkapan citra bumi dengan satelit Landsat. Satelit Landsat menghasilkan citra berkualitas tinggi, untuk seluruh dunia, setiap 16 hari. Citra ini disediakan oleh United States Geological Surveys [USGS] yang dapat diakses oleh umum melalui laman berikut EarthExplorer. Citra Landsat 8 memiliki resolusi piksel 28,5m, dengan satu band yang resolusi lebih tinggi dengan ukuran piksel 15m. Sekali melintas, satelit Landsat ini menangkap jalur citra selebar 185km, diukur di permukaan bumi. Jalur citra ini dipotong untuk memudahkan distribusi dan pengelolaan data. Setiap potongan jalur [scene] diberikan nomor jalur [path] dan nomor barisan [row], Kawasan Metropolitan Malang Raya memiliki nomor jalur WRS_PATH = 118 dan nomor barusan WRS_ROW = 66. Data citra satelit Landsat 8 Kawasan Metropolitan Malang Raya diunduh pada bulan Agustus 2013 dan Juli 2019, dengan mempertimbangkan tutupan awan [cloud cover] paling sedikit.
Citra Landsat Kawasan Metropolitan Malang Raya
Satelit landsat 8 memiliki sensor onboard Operational Land Imager [OLI] dan Thermal Infrared Sensor [TIRS] dengan jumlah saluran sebanyak 11 buah. Diantara saluran tersebut, 9 saluran [band 1-9] berada pada OLI dan 2 lainnya [band 10 dan 11] pada TIRS. Data citra Landsat 8 bisa dikombinasikan dengan ArcGIS untuk memiliki ketajaman gambar sesuai dengan kebutuhan analisa. Berdasarkan sumber dari web ESRI terdapat 10 kombinasi band Landsat yang bisa dibuat berdasarkan karakteristik dari masing-masing band untuk keperluan analisa. Kombinasi band Landsat dapat diakses melalui laman berikut Band Combination for Landsat 8.
A. Input Data
1. Pada tampilan awal ArcMap, perlu ditentukan sistem koordinat yang akan digunakan. Penentuan sistem koordinat ini penting untuk dilakukan, agar data shapefile yang diinput ke dalam ArcMap sesuai dengan posisi sebenarnya dan nantinya dapat dilakukan perhitungan geometri. Pada menu Table of Contents, klik kanan pada Layers, kemudian klik Properties. Akan muncul kotak dialog Data Frame Properties, pilih tab Coordinate System, dan klik opsi WGS 1984 UTM Zone 49S [zona untuk Kawasan Metropolitan Malang Raya].
4. Data citra yang diinput memiliki koordinat WGS 1984 UTM Zone 49N. Seharusnya koordinat yang benar adalah WGS 1984 UTM Zone 49S, karena posisi Kawasan Metropolitan Malang Raya berada di bawah garis ekuator. Oleh sebab itu, perlu dilakukan proyeksi ulang sistem koordinat untuk mengoreksi sistem koordinat data citra.
B. Proyeksi Sistem Koordinat
1. Sebelum dikoreksi sistem koordinatnya, dilakukan penggabungan citra dari band 1 sampai band 7, menggunakan tool Composite Bands. Band 8 tidak digabungkan, karena merupakan citra pankromatik dengan resolusi spasial 15m x 15m untuk penajaman citra. Fungsi Composite Bands terletak di menu ArcToolBox > Data Management Tools > Raster > Raster Processing > Composite Bands. Setelah muncul kotak dialog Composite Bands, kemudian input data citra band 1 sampai band 7, kemudian tentukan folder penyimpanan, klik Save dan klik OK.
C. Area of Interest Pemetaan [Clip]
1. Citra landsat yang telah dikomposit dan dikoreksi sistem koordinatnya, akan diambil sesuai wilayah yang dianalisis, yakni Kawasan Metropolitan Malang Raya. Klik icon Add Data kemudian input shapefile batas administrasi Kota Malang, Kabupaten Malang, dan Kota Batu yang sudah di Union ke dalam layer kerja. Kemudian langkah clipping data citra sesuai bentuk wilayah administrasi, dilakukan dengan tools Clip. Klik icon ArcToolBox > Data Management Tools > Raster > Raster Processing > Clip. Pada kotak dialog Clip, input data citra hasil komposit yang sudah dikoreksi sistem koordinatnya ke dalam kolom Input Raster, dan masukkan data shapefile wilayah administrasi ke dalam kolom Output Extent. Centang Use Input Features for Clipping Geometry, kemudian klik OK.
Data Citra Landsat RGB 4-3-2 Tahun 2013 [Kiri] dan Tahun 2019 [Kanan]
D. Penajaman Citra Landsat [Pan Sharpening]
1. Langkah selanjutnya adalah penajaman data citra Landsat, dengan menambahkan band 8 sebagai data pankromatik dengan resolusi 15m x 15m. Penambahan data pankromatik ini dilakukan dengan menggunakan ArcToolBox > Data Management Tools > Raster > Raster Processing > Create Pan-sharpened Raster Dataset. Masukkan data citra landsat Kawasan Metropolitan Malang Raya di kolom Input Raster, masukkan data citra band 8 ke dalam kolom Panchromatic Image, kemudian klik OK.
E. Metode Klasifikasi Terbimbing [Supervised Classification]
1. Identifikasi tutupan lahan dilakukan dengan menggunakan metode Klasifikasi Terbimbing [Supervised Classification]. Metode ini merupakan cara mengklasifikasi berdasarkan sampel yang dibuat untuk setiap jenis tutupan lahan. Sebelum melakukan klasifikasi, langkah yang harus dilakukan terlebih dahulu adalah mengaktifkan ekstensi Spatial Analyst melalui tab menu Customize > Extensions > centang kolom Spatial Analyst > kemudian klik Close.
4. Poin penting yang harus diperhatikan dalam pembuatan training sample adalah harus representatif. Banyaknya training sample bukan jaminan menghasilkan model yang baik jika training sample yang dibuat tidak akurat. Jumlah training sample yang dibutuhkan disesuaikan dengan keberagaman objek. Idealnya semakin banyak training sample akan memberikan hasil yang semakin baik.
6. Langkah selanjutnya adalah membuat sampel dalam format .gsg sebagai bahan dalam klasifikasi terbimbing. Klik menu Create a signature file pada toolbar Training Sample Manager > tentukan folder penyimpanan > tulis nama sampel > save as signature files > klik Save. Karena menggunakan dua periode waktu, maka lakukan digitasi training sample yang dibuat untuk masing-masing pada tahun 2013 dan tahun 2019.
Tutupan Lahan Tahun 2013 [Kiri] dan Tutupan Lahan Tahun 2019 [Kanan]
F. Identifikasi Perubahan Lahan
1. Perubahan lahan dapat diidentifikasi dengan membandingkan data tutupan lahan tahun 2013 dengan data tutupan lahan tahun 2019. Salah satu cara yang bisa dilakukan adalah dengan menghitung luasan setiap jenis tutupan lahan. Langkah yang dilakukan adalah merubah data raster menjadi data vektor berupa format shapefile. Klik ikon ArcToolbox > Conversion Tools > From Raster > Raster to Polygon. Setelah muncul kotak dialog Raster to Polygon, isi kolom Input raster dengan data raster hasil klasfikasi terbimbing, kemudian pilih folder penyimpanan, lalu klik OK.
3. Klik kanan pada data shapefile hasil klasifikasi > Open Attribute Table > Kotak dialog Table > Table Options > Add Field > Isi nama untuk keterangan jenis tutupan lahan > Type Text > Isi kolom Length [jumlah huruf] sesuai kebutuhan > Klik OK. Lakukan hal yang sama untuk field luas, hanya saja yang berbeda setelah isi nama untuk keterangan luasan dalam hektar > Type Double > Isi kolom Precision dan Scale sesuai kebutuhan > Klik OK.
5. Langkah selanjutnya adalah menghitung luasan tiap jenis tutupan lahan menggunakan field dengan nama Luas yang telah dibuat. Caranya adalah dengan klik kanan pada field luas > Calculate Geometry > Property Area > Coordinate System > Units hectares [ha] > klik OK.
TAHUN 2013 |
TAHUN 2019 |
|||||||
No |
Tutupan Lahan |
Luas [ha] |
Persentase |
No |
Tutupan Lahan |
Luas [ha] |
Persentase |
|
1 |
Permukiman |
31683.65 |
8.40% |
1 |
Permukiman |
33016.65 |
8.76% |
|
2 |
Industri |
3186.95 |
0.85% |
2 |
Industri |
3350.66 |
0.89% |
|
3 |
Hutan |
176036.99 |
46.70% |
3 |
Hutan |
161598.50 |
42.87% |
|
4 |
Pertanian |
145987.24 |
38.73% |
4 |
Pertanian |
158016.78 |
41.92% |
|
5 |
Perairan |
3919.90 |
1.04% |
5 |
Perairan |
4743.48 |
1.26% |
|
6 |
Lahar |
1774.00 |
0.47% |
6 |
Lahar |
2576.54 |
0.68% |
|
7 |
Awan |
8055.16 |
2.14% |
7 |
Awan |
8878.71 |
2.36% |
|
8 |
Bayangan |
6330.84 |
1.68% |
8 |
Bayangan |
4793.41 |
1.27% |
|
Total |
376974.73 |
100% |
Total |
376974.73 |
100% |
8. Dapat diketahui bahwa terjadi pengurangan tutupan lahan hutan sebesar 4,1 %, dan terjadi peningkatan tutupan lahan permukiman sebesar 0,3%, industri sebesar 0,04%, dan pertanian sebesar 3,2%. Secara umum untuk luasan hijau masih memenuhi persentase minimum sebesar 30%.
KESIMPULAN
Dapat dilihat bahwa perkembangan wilayah, dalam hal ini lahan terbangun di Kawasan Malang Raya terpusat di daerah dekat perkotaan [Kota Malang dan Kota Batu]. Seiring dengan bertambahnya jumlah penduduk, terlihat bahwa kawasan terbangun di Kawasan Metropolitan Malang Raya terjadi peningkatan yang cukup signifikan. Perkembangan daerah perkotaan awalnya hanya berada di pusat inti kota dan sebagian terdapat di kawasan perkotaan. Namun, seiring bertambahnya waktu, perkembangannya meningkat terutama di kawasan lingkar kota.
Pada tahun 2013, ekspansi penduduk terlihat mulai berkembang ke arah utara dan barat. Bagian utara merupakan akses menuju Ibu Kota Provinsi Jawa Timur yaitu Kota Surabya, sehingga perkembangan mengikuti jalan nasional. Akses transportasi jalan telah membentuk perkembangan pola permukiman sepanjang jalan. Adanya jalan tol Malang - Pandaan yang menghubungkan wilayah metropolitan Malang Raya dan Kota Surabaya juga telah mempengaruhi perkembangan perkotaan. Bagian barat berkembang karena pengaruh Kota Batu sebagai Kota Wisata yang menjangkau wilayah regional, sehingga banyak terjadi alih fungsi lahan pertanian menjadi pemukiman, begitu juga dengan lahan hutan yang berubah menjadi lahan pertanian, karena kebutuhan masyarakat yang semakin.
Hasil identifikasi perubahan lahan ini dapat dimanfaatkan sebagai monitoring wilayah untuk kebijakan pemerintah ke depan, serta dengan analisis lebih lanjut dapat digunakan untuk memprediksi perkembangan perkotaan di tahun berikutnya. Pemantauan ekspansi daerah perkotaan menggunakan data citra satelit dan SIG merupakan suatu hal yang penting dalam kaitannya dengan perencanaan tata ruang dan lingkungan. Penggunaan citra satelit yang memiliki resolusi yang tinggi memungkinkan untuk memantau perubahan guna lahan secara rinci.
Beberapa hal yang perlu dilakukan dalam rangka peningkatan hasil studi dan kemungkinan studi lebih lanjut, diantaranya yaitu perlunya melakukan proses validasi dan koreksi data, baik melalui data observasi lapangan maupun data dari sumber lain; pendeteksian dengan menggunakan data citra satelit yang memiliki resolusi yang tinggi dan metode analisis klasifikasi yang lebih baik; identifikasi faktor yang mempengaruhi perubahan lahan, misalnya pertumbuhan jumlah penduduk, harga lahan, kondisi ekonomi, kondisi infrastruktur, dan perubahan kondisi sosial.