Keterbatasan Data Hujan Pos dan Alternatif GPM IMERG

Mengapa data satelit semakin relevan dalam analisis hidrologi Indonesia

Penulis

Irpan Chumaedi, M.T.

Diterbitkan

10 Mei 2026

1 Pendahuluan

Ketersediaan data curah hujan yang memadai merupakan fondasi utama dalam analisis hidrologi — mulai dari estimasi debit banjir, pemodelan erosi, hingga penyusunan dokumen AMDAL. Namun di Indonesia, praktisi hidrologi seringkali dihadapkan pada keterbatasan data pos hujan yang signifikan, terutama di wilayah Kalimantan, Sulawesi, Maluku, dan Papua.

Tulisan ini membahas keterbatasan tersebut secara kuantitatif dan memperkenalkan GPM IMERG sebagai alternatif data curah hujan berbasis satelit yang dapat melengkapi atau bahkan menggantikan data pos hujan dalam kondisi tertentu.

Konteks analisis

Pembahasan ini relevan untuk analisis hidrologi dalam kerangka studi AMDAL, penyusunan SID, serta pemodelan HEC-HMS pada DAS yang kekurangan data pos hujan.

2 Keterbatasan Data Pos Hujan di Indonesia

2.1 Distribusi Spasial Stasiun BMKG

Berdasarkan inventarisasi BMKG, terdapat 6.487 stasiun hujan aktif di seluruh Indonesia yang terdiri dari tiga jenis: Pos Hujan Kerjasama (PHK), Automatic AWS (AAWS), dan UPT.

Tampilkan kode R

# Data stasiun — embed langsung, tidak perlu file eksternal
df_sta <- do.call(rbind.data.frame, list(
  list('Genteng I','Pos Hujan Kerjasama','Jawa Timur','Kab. Banyuwangi', 171,113.94690, -8.34778),
  list('Huta Bagasan','Pos Hujan Kerjasama','Sumatera Utara','Kab. Asahan',NA_real_, 99.32000,  2.75000),
  list('Surat/Mojo','Pos Hujan Kerjasama','Jawa Timur','Kab. Kediri', 123,111.95030, -7.88194),
  list('Niniari','Pos Hujan Kerjasama','Maluku','Kab. Seram Bagian Barat',NA_real_,128.40000, -3.07000),
  list('Kerta Sari','Pos Hujan Kerjasama','Jawa Barat','Kab. Bekasi',NA_real_,107.28414, -6.20928),
  list('Wanadadi','Pos Hujan Kerjasama','Jawa Tengah','Kab. Banjarnegara',NA_real_,109.62580, -7.36860),
  list('Jonggol','Pos Hujan Kerjasama','Jawa Barat','Kab. Bogor', 137,107.06590, -6.46790),
  list('Jatibarang 2','Pos Hujan Kerjasama','Jawa Tengah','Kab. Brebes',  22,109.05470, -6.96297),
  list('Debut','Pos Hujan Kerjasama','Maluku','Kab. Maluku Tenggara',NA_real_,132.70000, -5.74000),
  list('Stageof Onowembo','Pos Hujan Kerjasama','Sumatera Utara','Kota Gunungsitoli',NA_real_, 97.55200,  1.29000),
  list('Bakunase','Pos Hujan Kerjasama','Nusa Tenggara Timur','Kota Kupang',   8,123.58722,-10.18754),
  list('Lemong','Pos Hujan Kerjasama','Lampung','Kab. Lampung Barat',NA_real_,105.20167, -4.60000),
  list('Bangko','Pos Hujan Kerjasama','Riau','Kab. Rokan Hilir',   4,100.80010,  2.17000),
  list('Gumbil','Pos Hujan Kerjasama','Kalimantan Selatan','Kab. Hulu Sungai Selatan',NA_real_,115.34970, -2.72127),
  list('Tinjowan Bibitan','Pos Hujan Kerjasama','Sumatera Utara','Kab. Simalungun',NA_real_, 99.48000,  3.09000),
  list('Kaliwringin','Pos Hujan Kerjasama','Jawa Tengah','Kab. Kendal', 205,110.19660, -7.10410),
  list('Tarsanabaru','Pos Hujan Kerjasama','Jawa Barat','Kab. Cirebon',   5,108.37330, -6.88330),
  list('Pakkat Hauagong','Pos Hujan Kerjasama','Sumatera Utara','Kab Humbang Hasundutan', 531, 98.46900,  2.17700),
  list('','Pos Hujan Kerjasama','Nanggroe Aceh Darussalam','Kab. Aceh Barat',NA_real_, 96.00000,  4.32000),
  list('Talaga','Pos Hujan Kerjasama','Jawa Barat','Kab. Majalengka', 636,108.36670, -6.96667),
  list('Bener Kelipah','Pos Hujan Kerjasama','Nanggroe Aceh Darussalam','Kab. Bener Meriah',NA_real_, 96.91570,  4.77421),
  list('Payung','Pos Hujan Kerjasama','Jawa Tengah','Kab. Pemalang',NA_real_,109.46210, -6.99827),
  list('Borong Rapoa','Pos Hujan Kerjasama','Sulawesi Selatan','Kab. Bulukumba', 707,120.04420, -5.39417),
  list('Batang Peranap','Pos Hujan Kerjasama','Riau','Kab. Indragiri Hulu',NA_real_,101.84170, -0.74842),
  list('Baning','Pos Hujan Kerjasama','Kalimantan Barat','Kab. Sintang',  29,111.50030, -0.06083),
  list('Sidaharja','Pos Hujan Kerjasama','Jawa Barat','Kab. Ciamis',  65,108.69536, -7.42385),
  list('Bah Jambi Afd 1','Pos Hujan Kerjasama','Sumatera Utara','Kab. Simalungun', 181, 99.16400,  2.93200),
  list('Unib','Pos Hujan Kerjasama','Bengkulu','Kota Bengkulu',  29,102.26890, -2.24111),
  list('Surantih','Pos Hujan Kerjasama','Sumatera Barat','Kab. Pesisir Selatan',  16,100.65080,  1.59889),
  list('Sumber Rejo','Pos Hujan Kerjasama','Lampung','Kota Bandar Lampung',NA_real_,105.09675, -4.67849),
  list('Rantau Laban','Pos Hujan Kerjasama','Sumatera Utara','Kota Tebing Tinggi',NA_real_, 99.16700,  3.38500),
  list('Pakembaran','Pos Hujan Kerjasama','Jawa Tengah','Kab. Tegal',NA_real_,109.13380, -6.98115),
  list('Setia Bakti','Pos Hujan Kerjasama','Lampung','Kab. Lampung Tengah',NA_real_,105.27500, -4.40000),
  list('Beduai','Pos Hujan Kerjasama','Kalimantan Barat','Kab. Sanggau',   6,110.39000,  0.70218),
  list('Karangploso','Pos Hujan Kerjasama','Jawa Timur','Kab. Malang', 575,112.59670, -7.90139),
  list('Situnggaling','Pos Hujan Kerjasama','Sumatera Utara','Kab. Karo',NA_real_, 98.52800,  2.91600),
  list('Binjai Barat','Pos Hujan Kerjasama','Sumatera Utara','Kota Binjai',  32, 98.49000,  3.62000),
  list('Bantar Kalong','Pos Hujan Kerjasama','Jawa Barat','Kab. Tasikmalaya', 540,108.08000, -7.07000),
  list('Losarang','Pos Hujan Kerjasama','Jawa Barat','Kab. Indramayu',   4,108.16670, -6.40000),
  list('Kuta Cot Glie','Pos Hujan Kerjasama','Nanggroe Aceh Darussalam','Kab. Aceh Besar',  49, 95.56000,  5.31000),
  list('Ngargoyoso','Pos Hujan Kerjasama','Jawa Tengah','Kab. Karanganyar', 830,111.11320, -7.60780),
  list('Nawangan','Pos Hujan Kerjasama','Jawa Timur','Kab. Pacitan', 662,111.19030, -7.99917),
  list('Tamansari/Kemuningsari','Pos Hujan Kerjasama','Jawa Timur','Kab. Jember',  36,113.61530, -8.28722),
  list('Pabuaran','Pos Hujan Kerjasama','Banten','Kab. Serang',NA_real_,106.07190, -6.20525),
  list('Pulerejo/Muneng','Pos Hujan Kerjasama','Jawa Timur','Kab. Madiun',  73,111.61610, -7.49889),
  list('Sekerak','Pos Hujan Kerjasama','Nanggroe Aceh Darussalam','Kab. Aceh Tamiang',NA_real_, 97.92510,  4.34230),
  list('Loji','Pos Hujan Kerjasama','Jawa Barat','Kab. Majalengka',  60,108.30000, -6.71667),
  list('Menes','Pos Hujan Kerjasama','Banten','Kab. Pandeglang',NA_real_,105.91990, -6.37667),
  list('Pesanggrahan/ Depok','Pos Hujan Kerjasama','DKI Jakarta','Kota Adm. Jakarta Selatan',  86,106.77230, -6.39718),
  list('Tenggarong Sebrang / Teluk Dalam','Pos Hujan Kerjasama','Kalimantan Timur','Kab. Kutai Kertanegara',  29,117.45560,  2.16353),
  list('Lawa','Pos Hujan Kerjasama','Kalimantan Selatan','Kab. Banjar',  39,115.00880, -3.24751),
  list('Salam','Pos Hujan Kerjasama','Jawa Tengah','Kab. Magelang',NA_real_,110.32100, -7.62500),
  list('Bonegunu','Pos Hujan Kerjasama','Sulawesi Tenggara','Kab. Buton Utara',NA_real_,122.98580, -4.75515),
  list('Spp Kalasey','Pos Hujan Kerjasama','Sulawesi Utara','Kab. Minahasa',  46,124.76830,  1.44639),
  list('Cukurguling/Lumbang','Pos Hujan Kerjasama','Jawa Timur','Kab. Pasuruan', 251,112.97860, -7.78472),
  list('Ciasem','Pos Hujan Kerjasama','Jawa Barat','Kab. Subang',  24,107.65000, -6.35000),
  list('Bonto Lebang','Pos Hujan Kerjasama','Sulawesi Selatan','Kab. Takalar',NA_real_,119.38690, -5.42444),
  list('Krueng Sabe','Pos Hujan Kerjasama','Nanggroe Aceh Darussalam','Kab. Aceh Jaya',  13, 95.65000,  4.61000),
  list('Mesuji Raya','Pos Hujan Kerjasama','Sumatera Selatan','Kab. Ogan Komering Ilir',NA_real_,105.01240, -3.20900),
  list('Rantau Pandan','Pos Hujan Kerjasama','Jambi','Kab. Bungo',NA_real_,102.31190, -1.70127),
  list('Pulau Rimau','Pos Hujan Kerjasama','Sumatera Selatan','Kab. Banyuasin',  15,104.46390, -2.60789),
  list('Pernantian','Pos Hujan Kerjasama','Sumatera Utara','Kab. Labuhanbatu Utara',NA_real_, 99.78500,  2.21400),
  list('Beji','Pos Hujan Kerjasama','DI Yogyakarta','Kab. Kulon Progo',NA_real_,110.16000, -7.87000),
  list('Wadaslintang','Pos Hujan Kerjasama','Jawa Tengah','Kab. Wonosobo', 317,109.81930, -7.54364),
  list('Pancur Bpp','Pos Hujan Kerjasama','Jawa Tengah','Kab. Rembang',  30,111.47200, -6.71669),
  list('Rakumpit','Pos Hujan Kerjasama','Kalimantan Tengah','Kota Palangkaraya',  49,113.67000, -1.75060),
  list('Tapung Hulu','Pos Hujan Kerjasama','Riau','Kab. Kampar',  88,100.93880,  0.66305),
  list('Batukliang','Pos Hujan Kerjasama','Nusa Tenggara Barat','Kab. Lombok Tengah',   7,116.30000, -8.57000),
  list('Condong','Pos Hujan Kerjasama','Jawa Timur','Kab. Probolinggo',  95,113.36970, -7.87083),
  list('Pulau Banyak','Pos Hujan Kerjasama','Nanggroe Aceh Darussalam','Kab. Aceh Singkil',  20, 97.28000,  2.19000),
  list('Cubadak','Pos Hujan Kerjasama','Sumatera Barat','Kab. Tanah Datar',NA_real_,100.54190,  0.46583),
  list('Prailiu','Pos Hujan Kerjasama','Nusa Tenggara Timur','Kab. Sumba Timur',   7,120.27970, -9.65778),
  list('Cimanuk','Pos Hujan Kerjasama','Banten','Kab. Pandeglang', 223,106.04380, -6.35035),
  list('Junjungan Sirih (Muaro Pingai)','Pos Hujan Kerjasama','Sumatera Barat','Kab. Solok',NA_real_,100.40920, -0.71222),
  list('Bunder','Pos Hujan Kerjasama','Jawa Barat','Kab. Indramayu',NA_real_,108.26065, -6.47600),
  list('Muara Beliti 1','Pos Hujan Kerjasama','Sumatera Selatan','Kab. Musi Rawas',NA_real_,103.41640, -3.41636),
  list('Wanasabakidul','Pos Hujan Kerjasama','Jawa Barat','Kab. Cirebon', 116,108.56670, -6.71670),
  list('Teritip','Pos Hujan Kerjasama','Kalimantan Timur','Kota Balikpapan',  10,116.87080, -1.26611),
  list('Tualang','Pos Hujan Kerjasama','Riau','Kab. Siak',  37,101.81670,  0.98330),
  list('Satrian/Kanigoro','Pos Hujan Kerjasama','Jawa Timur','Kab. Blitar', 168,112.22030, -8.14472),
  list('Sp. Mamplam','Pos Hujan Kerjasama','Nanggroe Aceh Darussalam','Kab. Bireuen',NA_real_, 96.44000,  5.21000),
  list('Pangkalan Lesung','Pos Hujan Kerjasama','Riau','Kab. Pelalawan',NA_real_,102.11920, -0.02884),
  list('Musi Landas','Pos Hujan Kerjasama','Sumatera Selatan','Kab. Banyuasin',  34,104.61340, -2.91678),
  list('Sei Balei','Pos Hujan Kerjasama','Sumatera Utara','Kab. Batu Bara',NA_real_, 99.18000,  3.11000),
  list('Pojok','Pos Hujan Kerjasama','Jawa Timur','Kab. Blitar', 259,112.20780, -8.06944),
  list('Kemantren/Jabung','Pos Hujan Kerjasama','Jawa Timur','Kab. Malang', 530,112.75190, -7.94639),
  list('Sri Gunung','Pos Hujan Kerjasama','Sumatera Selatan','Kab. Musi Banyuasin',  13,104.08600, -2.58656),
  list('Patilangio','Pos Hujan Kerjasama','Gorontalo','Kab. Pahuwato',NA_real_,121.86000,  0.51000),
  list('Gembungan/Slemanan','Pos Hujan Kerjasama','Jawa Timur','Kab. Blitar', 195,112.09970, -8.00194),
  list('Kappang','Pos Hujan Kerjasama','Sulawesi Selatan','Kab. Maros', 321,119.75220, -5.01417),
  list('Kedungkamal','Pos Hujan Kerjasama','Jawa Tengah','Kab. Purworejo',NA_real_,109.92060, -7.79101),
  list('Tangunan','Pos Hujan Kerjasama','Jawa Timur','Kab. Mojokerto',  70,112.46420, -7.54306),
  list('Sungai Apit','Pos Hujan Kerjasama','Riau','Kab. Siak',  29,102.13330,  0.75000),
  list('Bpp Ngadirejo','Pos Hujan Kerjasama','Jawa Tengah','Kab. Temanggung',NA_real_,110.06100, -7.23600),
  list('Pakuan Ratu','Pos Hujan Kerjasama','Lampung','Kab. Way Kanan',NA_real_,104.72843, -4.41624),
  list('Singkawang Timur','Pos Hujan Kerjasama','Kalimantan Barat','Kota Singkawang',   3,109.02480,  0.88422),
  list('Amahai','Pos Hujan Kerjasama','Maluku','Kab. Maluku Tengah',NA_real_,128.94000, -3.33000),
  list('Tongoa','Pos Hujan Kerjasama','Sulawesi Tengah','Kab. Sigi',NA_real_,120.17100, -1.21870),
  list('Panongan','Pos Hujan Kerjasama','Jawa Barat','Kab. Majalengka',  41,108.25000, -6.73333),
  list('Ngelak','Pos Hujan Kerjasama','Jawa Tengah','Kab. Karanganyar', 799,111.12230, -7.66414),
  list('Pos Uptd Saumlaki','Pos Hujan Kerjasama','Maluku','Kab Maluku Tenggara Barat',   8,131.31000, -7.97820),
  list('Tegorejo','Pos Hujan Kerjasama','Jawa Tengah','Kab. Kendal',  32,110.14870, -6.97747),
  list('Waru','Pos Hujan Kerjasama','Kalimantan Timur','Kab. Penajam Paser Utara',   9,116.61610, -1.38910),
  list('Sawahan','Pos Hujan Kerjasama','Jawa Timur','Kab. Nganjuk', 569,111.77000, -7.72389),
  list('Stageof Winangun','Pos Hujan Kerjasama','Sulawesi Utara','Kota Manado', 124,124.84000,  1.45000),
  list('Sidoklumpuk/Sidoarjo','Pos Hujan Kerjasama','Jawa Timur','Kab. Sidoarjo',  56,112.72060, -7.44611),
  list('Kedungjati 1','Pos Hujan Kerjasama','Jawa Tengah','Kab. Grobogan',NA_real_,110.63370, -7.15889),
  list('Bpp Tanete Rilau','Pos Hujan Kerjasama','Sulawesi Selatan','Kab. Barru',   9,119.60810, -4.48056),
  list('Gebangan Slt','Pos Hujan Kerjasama','Jawa Tengah','Kab. Kendal',NA_real_,110.05800, -7.03600),
  list('Srimadono Smpk','Pos Hujan Kerjasama','Jawa Tengah','Kab. Kebumen',  26,109.78350, -7.71906),
  list('Telaga Tawang','Pos Hujan Kerjasama','Bali','Kab. Karangasem', 291,115.46000, -8.46000),
  list('Tungkal Ilir','Pos Hujan Kerjasama','Jambi','Kab. Tanjung Jabung Barat',  14,103.46810, -1.26283),
  list('Pulau Malan','Pos Hujan Kerjasama','Kalimantan Tengah','Kab. Katingan',NA_real_,113.11289, -7.73718),
  list('Boja','Pos Hujan Kerjasama','Jawa Tengah','Kab. Kendal', 316,110.27090, -7.10379),
  list('Parung Kuda','Pos Hujan Kerjasama','Jawa Barat','Kab. Sukabumi',NA_real_,106.76670, -6.85000),
  list('Pringgarata','Pos Hujan Kerjasama','Nusa Tenggara Barat','Kab. Lombok Tengah',   5,116.25000, -8.62000),
  list('Maligas','Pos Hujan Kerjasama','Sumatera Utara','Kab. Simalungun',NA_real_, 99.36300,  3.03800),
  list('Cigeulis','Pos Hujan Kerjasama','Banten','Kab. Pandeglang',NA_real_,105.64000, -6.55000),
  list('Campalagoe/ Bonto Jaya','Pos Hujan Kerjasama','Sulawesi Selatan','Kab. Bantaeng', 256,119.87440, -5.53417),
  list('Boyolangu','Pos Hujan Kerjasama','Jawa Timur','Kab. Tulungagung', 135,111.89560, -8.11056),
  list('Turangi Lama','Pos Hujan Kerjasama','Sumatera Utara','Kab. Langkat',NA_real_, 98.28000,  3.45000),
  list('Tiworo Kepulauan3','Pos Hujan Kerjasama','Sulawesi Tenggara','Kab. Muna',NA_real_,122.43300, -4.77249),
  list('Kali Orang','Pos Hujan Kerjasama','Kalimantan Timur','Kab. Kutai Timur',NA_real_,116.98750,  0.99971),
  list('Nanga Mau','Pos Hujan Kerjasama','Kalimantan Barat','Kab. Sintang', 160,113.55000, -6.94000),
  list('Balongjeruk','Pos Hujan Kerjasama','Jawa Timur','Kab. Kediri',  95,112.19030, -7.67222),
  list('Curah Malang','Pos Hujan Kerjasama','Jawa Timur','Kab. Jember',  51,113.56720, -8.24139),
  list('Sukajadi','Pos Hujan Kerjasama','Banten','Kota Serang',NA_real_,106.22960, -6.18920),
  list('Bunga Melur','Pos Hujan Kerjasama','Jawa Barat','Kab. Cianjur',1050,106.97470, -7.04762),
  list('Bulu Cenrana/Pituriawa','Pos Hujan Kerjasama','Sulawesi Selatan','Kab. Sidenreng Rappang',NA_real_,119.95640, -3.82366),
  list('Sorkam','Pos Hujan Kerjasama','Sumatera Utara','Kab. Tapanuli Tengah',NA_real_, 98.58000,  1.90000),
  list('Gembo','Pos Hujan Kerjasama','Jawa Tengah','Kab. Banjarnegara',NA_real_,109.76340, -7.39702),
  list('Oro-Oro Pule','Pos Hujan Kerjasama','Jawa Timur','Kab. Pasuruan',  84,112.83810, -7.74417),
  list('Sangata Utara','Pos Hujan Kerjasama','Kalimantan Timur','Kab. Kutai Timur',  65,117.05190,  0.54022),
  list('Pintu Air Cideng','Pos Hujan Kerjasama','DKI Jakarta','Kota Adm. Jakarta Pusat',  21,106.80670, -6.17360),
  list('Boto/Botogerdu','Pos Hujan Kerjasama','Jawa Timur','Kab. Probolinggo', 250,113.11360, -7.84833),
  list('Petungasri/Kasri','Pos Hujan Kerjasama','Jawa Timur','Kab. Pasuruan', 212,112.68810, -7.65167),
  list('Cambai','Pos Hujan Kerjasama','Kep. Bangka Belitung','Kab. Bangka Tengah',NA_real_,106.15640, -2.24583),
  list('Sememu','Pos Hujan Kerjasama','Jawa Timur','Kab. Lumajang', 171,113.13250, -8.19750),
  list('Ktr. Kel. Kairagi Ii','Pos Hujan Kerjasama','Sulawesi Utara','Kota Manado',  67,124.90000,  1.50000),
  list('Kalipentung/ Tanjungtirto','Pos Hujan Kerjasama','DI Yogyakarta','Kab. Sleman', 137,110.45850, -7.78681),
  list('Pagar Jati','Pos Hujan Kerjasama','Bengkulu','Kab. Bengkulu Tengah', 136,102.40220, -3.63610),
  list('Rejosari 2','Pos Hujan Kerjasama','Lampung','Kab. Lampung Selatan',NA_real_,105.16167, -4.30111),
  list('Bontomatene','Pos Hujan Kerjasama','Sulawesi Selatan','Kab. Jeneponto',  94,119.77500, -5.61444),
  list('Blahkiuh','Pos Hujan Kerjasama','Bali','Kab. Badung', 193,115.21000, -8.51000),
  list('Cipero','Pos Hujan Kerjasama','Jawa Tengah','Kab. Tegal',NA_real_,109.30800, -6.98782),
  list('Lebak Wangi','Pos Hujan Kerjasama','Jawa Barat','Kab. Bogor',  23,106.78000, -6.13000),
  list('Padalarang','Pos Hujan Kerjasama','Jawa Barat','Kab. Bandung Barat',NA_real_,107.47735, -6.85347),
  list('Kd Labuhan Uki','Pos Hujan Kerjasama','Sulawesi Utara','Kab. Bolaang Mongondow',   6,123.94500,  0.86139),
  list('Sumurpanggang','Pos Hujan Kerjasama','Jawa Tengah','Kota Tegal',NA_real_,109.10430, -6.87830),
  list('Mompang','Pos Hujan Kerjasama','Sumatera Utara','Kab. Mandailing Natal',NA_real_, 99.59400,  0.84600),
  list('AWS Krui','aaws','Lampung','Kab. Lampung Barat',NA_real_,103.93021, -5.12924),
  list('ARG Panti','aaws','Jawa Timur','Kab. Jember',NA_real_,113.63710, -8.14961),
  list('AAWS Rembang','aaws','Jawa Tengah','Kab. Rembang',NA_real_,111.33484, -6.70188),
  list('ARG Pandeglang','aaws','Banten','Kab. Pandeglang',NA_real_,106.10207, -6.30380),
  list('ARG Sigaluh','aaws','Jawa Tengah','Kab. Banjarnegara',NA_real_,109.75300, -7.39200),
  list('AWS Kemayoran','aaws','DKI Jakarta','Kota Adm. Jakarta Pusat',NA_real_,106.85290, -6.13343),
  list('ARG Sampang Madura','aaws','Jawa Timur','Kab. Sampang',NA_real_,113.47667, -6.89578),
  list('ARG Kintap','aaws','Kalimantan Selatan','Kab. Tanah Laut',NA_real_,115.28959, -3.86192),
  list('ARG Kedungombo','aaws','Jawa Tengah','Kab. Grobogan',NA_real_,110.83634, -7.24037),
  list('ARG Lebak Bulus','aaws','DKI Jakarta','Kota Adm. Jakarta Selatan',NA_real_,106.77737, -6.30417),
  list('AWS Kulon Progo','aaws','DI Yogyakarta','Kab. Kulon Progo',NA_real_,110.12022, -7.89177),
  list('ARG Sedinginan','aaws','Riau','Kab. Rokan Hilir',NA_real_,101.27300,  1.34210),
  list('AWS Stamet Paloh Sambas','aaws','Kalimantan Barat','Kab. Sambas',NA_real_,109.32210,  1.73740),
  list('ARG Tomang','aaws','DKI Jakarta','Kota Adm. Jakarta Barat',NA_real_,106.78000, -6.16670),
  list('AAWS Bojongmangu','aaws','Jawa Barat','Kab. Bekasi',NA_real_,107.19200, -6.43187),
  list('ARG OKU','aaws','Sumatera Selatan','Kab. Ogan Komering Ulu',NA_real_,103.97150, -4.76706),
  list('ARG Batu','aaws','Jawa Timur','Kab. Malang',NA_real_,112.52310, -7.84123),
  list('AAWS Sruweng Kebumen','aaws','Jawa Tengah','Kab. Kebumen',NA_real_,109.60476, -7.66142),
  list('ARG Buleleng','aaws','Bali','Kab. Buleleng',NA_real_,115.10050, -8.13499),
  list('AWS Lamongan','aaws','Jawa Timur','Kab. Lamongan',NA_real_,112.32972, -7.09389),
  list('AWS Juwana Pati','aaws','Jawa Tengah','Kab. Pati',NA_real_,111.14600, -6.72100),
  list('AWS Sungai Tabuk','aaws','Kalimantan Selatan','Kab. Banjar',NA_real_,114.68182, -3.32653),
  list('AWS Jagorawi Bogor','aaws','Jawa Barat','Kab. Bogor',NA_real_,106.52164, -6.27635),
  list('AAWS Tuban','aaws','Jawa Timur','Kab. Tuban',NA_real_,112.01381, -6.88317),
  list('AWS Gorontalo','aaws','Gorontalo','Kab. Gorontalo',NA_real_,122.85260,  0.63843),
  list('AWS Maritim Panjang Lampung','aaws','Lampung','Kota Bandar Lampung',NA_real_,105.31050, -5.45553),
  list('AWS Stamet Tanjung Pinang','aaws','Kep. Riau','Kota Tanjung Pinang',NA_real_,104.52600,  0.91840),
  list('ARG Pasir Pangaraian','aaws','Riau','Kab. Rokan Hulu',NA_real_,101.01000,  0.19340),
  list('ARG Solok','aaws','Sumatera Barat','Kab. Solok',NA_real_,100.62778, -0.94472),
  list('ARG Pekutatan','aaws','Bali','Kab. Jembrana',NA_real_,114.83830, -8.40194),
  list('AAWS Pondok Salam (AWS+)','aaws','Jawa Barat','Kab. Purwakarta',NA_real_,107.49780, -6.61872),
  list('AAWS Serang Curug','aaws','Banten','Kab. Serang',NA_real_,106.19020, -6.17122),
  list('AAWS Ciamis','aaws','Jawa Barat','Kab. Ciamis',NA_real_,108.36940, -7.29875),
  list('AWS Nangapinoh','aaws','Kalimantan Barat','Kab. Melawi',NA_real_,111.47000,  0.42000),
  list('AWS Tebo','aaws','Jambi','Kab. Tebo',NA_real_,102.17555, -1.38166),
  list('AWS Kuningan','aaws','Jawa Barat','Kab. Kuningan',NA_real_,108.45680, -6.99970),
  list('AWS Cisolok Sukabumi','aaws','Jawa Barat','Kab. Sukabumi',NA_real_,106.47630, -6.95950),
  list('ARG Karang Nunggal','aaws','Jawa Barat','Kab. Tasikmalaya',NA_real_,108.14000, -7.62530),
  list('AAWS Purwodadi','aaws','Jawa Tengah','Kab. Grobogan',NA_real_,110.91430, -7.10400),
  list('AWS Staklim Sampali','aaws','Sumatera Utara','Kab. Deli Serdang',NA_real_, 98.78000,  3.61600),
  list('AWS Koba','aaws','Kep. Bangka Belitung','Kab. Bangka Tengah',NA_real_,106.41670, -2.52330),
  list('ARG SMPK Lahat','aaws','Sumatera Selatan','Kab. Lahat',NA_real_,103.43900, -3.98277),
  list('AAWS Jember','aaws','Jawa Timur','Kab. Jember',NA_real_,113.44444, -8.20750),
  list('ARG SMPK Pemalang','aaws','Jawa Tengah','Kab. Pemalang',NA_real_,109.46130, -6.89682),
  list('ARG Kandis','aaws','Riau','Kab. Siak',NA_real_,101.11519,  1.02854),
  list('AWS Staklim Pondok Betung','aaws','Banten','Kota Tangerang Selata',NA_real_,106.76000, -6.25000),
  list('AAWS Waikabubak (AWS+)','aaws','Nusa Tenggara Timur','Kab. Sumba Barat',NA_real_,119.39120, -9.62411),
  list('ARG SMPK Bunga Maya','aaws','Riau','Kab. Siak',NA_real_,102.11670,  0.78330),
  list('ARG SMPK Tanah Datar','aaws','Sumatera Barat','Kab. Tanah Datar',NA_real_,100.56778, -0.47444),
  list('AWS Selagai lingga','aaws','Lampung','Kab. Lampung Tengah',NA_real_,104.68037, -5.10120),
  list('ARG Cibaliung','aaws','Banten','Kab. Pandeglang',NA_real_,106.01540, -6.49350),
  list('ARG Pronojiwo','aaws','Jawa Timur','Kab. Lumajang',NA_real_,113.00430, -8.17200),
  list('AWS Tamiang Hulu','aaws','Nanggroe Aceh Darussalam','Kab. Aceh Tamiang',NA_real_, 97.86441,  4.27744),
  list('AWS Stamet Banyuwangi','aaws','Jawa Timur','Kab. Banyuwangi',NA_real_,114.36944, -8.21028),
  list('ARG Cariu','aaws','Jawa Barat','Kab. Bogor',NA_real_,107.14294, -6.16028),
  list('AWS Singkil','aaws','Sulawesi Utara','Kota Manado',NA_real_,124.85389,  1.50306),
  list('ARG Cisalak','aaws','Banten','Kab. Lebak',NA_real_,106.30059, -6.37080),
  list('AWS Stamet Branti','aaws','Lampung','Kab. Lampung Selatan',NA_real_,105.17571, -5.24253),
  list('ARG Aceh','aaws','Nanggroe Aceh Darussalam','Kab. Aceh Timur',NA_real_, 97.78570,  4.92426),
  list('AAWS GAW Agam','aaws','Sumatera Barat','Kab. Agam',NA_real_,100.38222,  0.25028),
  list('ARG Sekotong','aaws','Nusa Tenggara Barat','Kab. Lombok Barat',NA_real_,116.07900, -8.83769),
  list('ARG SMPK Musi Rawas','aaws','Sumatera Selatan','Kab. Musi Rawas',NA_real_,102.92247, -3.20283),
  list('ARG Tuban','aaws','Jawa Timur','Kab. Tuban',NA_real_,112.48361, -6.90039),
  list('AAWS Lape Sumbawa (AWS+)','aaws','Nusa Tenggara Barat','Kab. Sumbawa',NA_real_,117.61720, -8.62531),
  list('ARG Tanjung Jabung Barat','aaws','Jambi','Kab. Tanjung Jabung Barat',NA_real_,100.42669,  0.84200),
  list('ARG Waikabubak','aaws','Nusa Tenggara Timur','Kab. Belu',NA_real_,124.89250,  9.10833),
  list('ARG Waduk Mrica','aaws','Jawa Tengah','Kab. Banjarnegara',NA_real_,109.61300, -7.39100),
  list('AWS Ranowulu Bitung','aaws','Sulawesi Utara','Kota Bitung',NA_real_,125.12389,  1.46167),
  list('ARG SMPK Tasikmalaya','aaws','Jawa Barat','Kota Tasikmalaya',NA_real_,108.20310, -7.33082),
  list('ARG Rendang','aaws','Bali','Kab. Karangasem',NA_real_,115.42020, -8.41583),
  list('ARG Manokwari','aaws','Papua Barat','Kab. Manokwari',NA_real_,133.72890, -0.79186),
  list('AWS Borobudur','aaws','Jawa Tengah','Kab. Magelang',NA_real_,110.21397, -7.61412),
  list('AWS Mojokerto','aaws','Jawa Timur','Kab. Mojokerto',NA_real_,112.92000, -7.65944),
  list('AAWS Kandat Kediri (AWS+)','aaws','Jawa Timur','Kab. Kediri',NA_real_,112.03361, -7.93389),
  list('ARG Siak','aaws','Riau','Kab. Siak',NA_real_,102.11100,  1.11490),
  list('ARG Ratahan','aaws','Sulawesi Utara','Kab. Minahasa Tenggara',NA_real_,124.79028,  1.02250),
  list('ARG SMPK Bojong Picung','aaws','Jawa Barat','Kab. Cianjur',NA_real_,107.27380, -6.83699),
  list('ARG Cidaun','aaws','Jawa Barat','Kab. Cianjur',NA_real_,107.36083, -7.49167),
  list('ARG SMPK Jombang','aaws','Jawa Timur','Kab. Jombang',NA_real_,112.23300, -7.49361),
  list('ARG Sorolangun','aaws','Jambi','Kab. Sarolangun',NA_real_,102.92231, -1.07486),
  list('AWS Stamet Serang','aaws','Banten','Kota Serang',NA_real_,106.13180, -6.11127),
  list('AAWS Kroya Cilacap','aaws','Jawa Tengah','Kab. Cilacap',NA_real_,109.24522, -7.59953),
  list('AWS Stageof Karangkates','aaws','Jawa Timur','Kab. Malang',NA_real_,112.77194, -8.15218),
  list('ARG Alas Sumbawa','aaws','Nusa Tenggara Barat','Kab. Sumbawa',NA_real_,116.99369, -8.51358),
  list('AWS Sorong','aaws','Papua Barat','Kota Sorong',NA_real_,131.29030, -0.89167),
  list('AAWS Bantul','aaws','DI Yogyakarta','Kab. Bantul',NA_real_,110.37470, -7.72167),
  list('AAWS Singingi Riau','aaws','Riau','Kab. Kuantan Singingi',NA_real_,101.48859, -0.59824),
  list('AWS Bamaang Hulu','aaws','Kalimantan Tengah','Kab. Kotawaringin Timur',NA_real_,112.97500, -2.50182),
  list('ARG Poso','aaws','Sulawesi Tengah','Kab. Poso',NA_real_,120.90830, -1.91720),
  list('AAWS Gerokgak (AWS+)','aaws','Bali','Kab. Buleleng',NA_real_,114.68000, -8.15000),
  list('AAWS Oelbubuk','aaws','Nusa Tenggara Timur','Kab. Timor Tengah Selatan',NA_real_,124.26700, -9.75900),
  list('ARG Sumber Probolinggo','aaws','Jawa Timur','Kab. Probolinggo',NA_real_,113.37889, -7.78194),
  list('AAWS Tabalong','aaws','Kalimantan Selatan','Kab. Tabalong',NA_real_,115.39008, -2.13466),
  list('AWS Kanigoro Blitar','aaws','Jawa Timur','Kab. Blitar',NA_real_,112.22260, -8.14288),
  list('ARG Pinolosian','aaws','Sulawesi Utara','Kab. Bolaang Mongondow Se',NA_real_,124.11750,  0.39417),
  list('ARG Tapanuli','aaws','Sumatera Utara','Kab. Tapanuli Selatan',NA_real_, 99.22000,  1.40000),
  list('AWS Karanganyar','aaws','Jawa Tengah','Kab. Karanganyar',NA_real_,111.06700, -7.61500),
  list('ARG Bengkalis','aaws','Riau','Kab. Bengkalis',NA_real_,102.08400,  1.29580),
  list('ARG Tulis','aaws','Jawa Tengah','Kab. Banjarnegara',NA_real_,109.79800, -7.32500),
  list('AWS Pekalongan','aaws','Jawa Tengah','Kab. Pekalongan',NA_real_,109.58700, -7.10400),
  list('AWS Stamet Cilacap','aaws','Jawa Tengah','Kab. Cilacap',NA_real_,109.01480, -7.71860),
  list('AAWS Sukodono (AWS+)','aaws','Jawa Tengah','Kab. Sragen',NA_real_,110.96620, -7.32750),
  list('AAWS Penajam Paser Utara','aaws','Kalimantan Timur','Kab. Penajam Paser Utara',NA_real_,116.65650, -1.36163),
  list('AWS Tiris Probolinggo','aaws','Jawa Timur','Kab. Probolinggo',NA_real_,113.40111, -7.96306),
  list('ARG Nusa Penida','aaws','Bali','Kab. Klungkung',NA_real_,115.49940, -8.67455),
  list('ARG Manggarai','aaws','DKI Jakarta','Kota Adm. Jakarta Selatan',NA_real_,106.84870, -6.20748),
  list('ARG SMPK Sragen','aaws','Jawa Tengah','Kab. Sragen',NA_real_,111.04380, -7.52544),
  list('AWS Sintang','aaws','Kalimantan Barat','Kab. Sintang',NA_real_,111.47000,  0.06000),
  list('AAWS Rokan Hilir','aaws','Riau','Kab. Rokan Hilir',NA_real_,101.02467,  1.85911),
  list('AAWS Tabanan','aaws','Bali','Kab. Tabanan',NA_real_,115.13917, -8.46083),
  list('ARG Baturiti Tabanan','aaws','Bali','Kab. Tabanan',NA_real_,115.16000, -8.26150),
  list('AWS Staklim Banjarbaru','aaws','Kalimantan Selatan','Kota Banjarbaru',NA_real_,114.84084, -3.46225),
  list('AAWS Kampar Kiri Riau','aaws','Riau','Kab. Kampar',NA_real_,100.25300,  0.34170),
  list('ARG Kangean','aaws','Jawa Timur','Kab. Sumenep',NA_real_,115.47111, -6.85867),
  list('AWS Tambilahan','aaws','Riau','Kab. Indragiri Hilir',NA_real_,103.13473, -0.28358),
  list('ARG Karanggede','aaws','Jawa Tengah','Kab. Boyolali',NA_real_,110.64800, -7.36200),
  list('ARG Ciomas','aaws','Banten','Kab. Serang',NA_real_,106.25109, -6.20400),
  list('ARG Majenang','aaws','Jawa Tengah','Kab. Cilacap',NA_real_,108.75810, -7.29584),
  list('AWS Putussibau','aaws','Kalimantan Barat','Kab. Kapuas Hulu',NA_real_,112.93470,  0.83895),
  list('AAWS Rappang','aaws','Sulawesi Selatan','Kab. Sidenreng Rappang',NA_real_,119.80465, -3.88929),
  list('ARG Kolaka','aaws','Sulawesi Tenggara','Kab. Kolaka',NA_real_,121.67747, -4.11241),
  list('AWS Cepu Blora','aaws','Jawa Tengah','Kab. Blora',NA_real_,111.57599, -7.17012),
  list('AAWS Banyumas','aaws','Jawa Tengah','Kab. Banyumas',NA_real_,109.30510, -7.52870),
  list('ARG Muaro Jambi','aaws','Jambi','Kab. Muaro Jambi',NA_real_,103.47430, -1.98370),
  list('ARG Muara Bungo','aaws','Jambi','Kab. Bungo',NA_real_,101.88840, -1.60930),
  list('AAWS Banyuwangi','aaws','Jawa Timur','Kab. Banyuwangi',NA_real_,114.14445, -8.37910),
  list('AWS Stamet Citeko','aaws','Jawa Barat','Kab. Bogor',NA_real_,106.93510, -6.69796),
  list('AWS Stamet Bima','aaws','Nusa Tenggara Barat','Kab. Bima',NA_real_,118.67940, -8.51661),
  list('ARG SMPK Singamerta','aaws','Banten','Kab. Serang',NA_real_,106.24058, -6.11928),
  list('AWS Staklim Maros','aaws','Sulawesi Selatan','Kab. Maros',NA_real_,119.57879, -4.94593),
  list('AWS Maritim Merak','aaws','Banten','Kota Cilegon',NA_real_,105.98530, -5.91889),
  list('ARG Minasatene','aaws','Sulawesi Selatan','Kab. Pangkajene Kepulauan',NA_real_,119.57420, -4.83161),
  list('ARG Mojokerto','aaws','Jawa Timur','Kab. Mojokerto',NA_real_,112.42778, -7.59139),
  list('ARG Betun','aaws','Nusa Tenggara Timur','Kab. Belu',NA_real_,124.90750,  9.56500),
  list('ARG Pulosari','aaws','Jawa Tengah','Kab. Pemalang',NA_real_,109.27074, -7.18731),
  list('AAWS Sukoharjo','aaws','Jawa Tengah','Kab. Sukoharjo',NA_real_,110.71680, -7.57697),
  list('ARG Gandus','aaws','Sumatera Selatan','Kota Palembang',NA_real_,104.72190, -3.01480),
  list('AWS Galela','aaws','Maluku Utara','Kab. Halmahera Utara',NA_real_,127.83333,  1.81667),
  list('AWS Mayang Jember','aaws','Jawa Timur','Kab. Jember',NA_real_,114.33778, -8.17693),
  list('AWS Prancak Jembrana','aaws','Bali','Kab. Jembrana',NA_real_,114.64000, -8.40000),
  list('AWS Bondowoso','aaws','Jawa Timur','Kab. Bondowoso',NA_real_,114.26556, -7.85608),
  list('ARG Saptosari','aaws','DI Yogyakarta','Kab. Gunung Kidul',NA_real_,110.51010, -8.04752),
  list('ARG Randublatung','aaws','Jawa Tengah','Kab. Blora',NA_real_,111.41400, -7.19400),
  list('AWS Murung','aaws','Kalimantan Tengah','Kab. Murung Raya',NA_real_,114.75280,  0.50550),
  list('AAWS Lore Piore (AWS+)','aaws','Sulawesi Tengah','Kab. Poso',NA_real_,120.19000,  1.31000),
  list('AWS Labuha','aaws','Maluku Utara','Kab. Halmahera Selatan',NA_real_,127.50080, -0.63796),
  list('AWS Stageof Ternate','aaws','Maluku Utara','Kota Ternate',NA_real_,127.36630,  0.76580),
  list('AWS Patean Kendal','aaws','Jawa Tengah','Kab. Kendal',NA_real_,110.09000, -7.09100),
  list('AWS Tanggul Jember','aaws','Jawa Timur','Kab. Jember',NA_real_,113.75806, -8.16300),
  list('AAWS Selogiri (AWS+)','aaws','Jawa Tengah','Kab. Wonogiri',NA_real_,110.88970, -7.78156),
  list('Stasiun Meteorologi Maritim Lampung','UPT','Lampung','Kota Bandar Lampung',   1,105.33000, -5.47000),
  list('Stasiun Meteorologi Binaka','UPT','Sumatera Utara','Kab. Nias', 175, 97.71000,  1.16502),
  list('Stasiun Meteorologi Gamar Malamo','UPT','Maluku Utara','Kab. Halmahera Utara',  29,127.50000,  1.49000),
  list('Stasiun Meteorologi Hasanuddin','UPT','Sulawesi Selatan','Kab. Maros',  14,119.55000, -5.07000),
  list('Stasiun Meteorologi Rahadi Oesman','UPT','Kalimantan Barat','Kab. Ketapang',   9,109.97000, -1.80000),
  list('Stasiun Meteorologi Mararena','UPT','Papua','Kab. Sarmi',   3,138.75000, -1.85000),
  list('Stasiun Meteorologi Maritim Bitung','UPT','Sulawesi Utara','Kota Bitung',   3,125.18000,  1.44306),
  list('Stasiun Geofisika Tuntungan','UPT','Sumatera Utara','Kab. Deli Serdang',  86, 98.56000,  3.50100),
  list('Stasiun Meteorologi Sentani','UPT','Papua','Kab. Jayapura',  96,140.48000, -2.57000),
  list('Stasiun Geofisika Parapat','UPT','Sumatera Utara','Kab. Simalungun',1061, 98.92000,  2.69395),
  list('Stasiun Meteorologi Maritim Perak II','UPT','Jawa Timur','Kota Surabaya',   3,112.74000, -7.20937),
  list('Stasiun Meteorologi Mathilda Batlayeri','UPT','Maluku','Kab Maluku Tenggara Barat',  24,131.30000, -7.98000),
  list('Stasiun Meteorologi Frans Kaisiepo','UPT','Papua','Kab. Biak Numfor',   3,136.10000, -1.19069),
  list('Stasiun Meteorologi Sultan Muhammad Kaharuddin','UPT','Nusa Tenggara Barat','Kab. Sumbawa',   4,117.43000, -8.52000),
  list('Stasiun Geofisika Jakarta','UPT','DKI Jakarta','Kota Adm. Jakarta Pusat',   4,106.84000, -6.15568),
  list('Stasiun Meteorologi Yuvai Semaring','UPT','Kalimantan Timur','Kab. Nunukan',1090,115.68000,  3.86000),
  list('Stasiun Meteorologi Tjut Nyak Dien Meulaboh','UPT','Nanggroe Aceh Darussalam','Kab. Nagan Raya',   3, 96.25000,  4.04898),
  list('Stasiun Meteorologi Malikussaleh','UPT','Nanggroe Aceh Darussalam','Kab. Aceh Utara',  28, 96.93000,  5.22800),
  list('Balai Besar Meteorologi, Klimatologi, dan Geofisika Wilayah I','UPT','Sumatera Utara','Kota Medan',NA_real_, 98.64000,  3.53970),
  list('Stasiun Meteorologi Syamsudin Noor','UPT','Kalimantan Selatan','Kota Banjarmasin',  32,114.76000, -3.45000),
  list('Stasiun Meteorologi David Constantijn Saudale','UPT','Nusa Tenggara Timur','Kab. Rote Ndao', 174,123.07000,-10.73000),
  list('Stasiun Meteorologi Rendani','UPT','Papua Barat','Kab. Manokwari',   3,134.05000, -0.88000),
  list('Stasiun Geofisika Tangerang','UPT','Banten','Kota Tangerang',  14,106.38000, -6.10000),
  list('Stasiun Meteorologi Frans Sales Lega','UPT','Nusa Tenggara Timur','Kab. Manggarai',1070,120.45000, -8.63000),
  list('Stasiun Meteorologi Sanggu','UPT','Kalimantan Tengah','Kab. Barito Selatan',  37,114.90000, -1.67000),
  list('Stasiun Geofisika Mataie','UPT','Nanggroe Aceh Darussalam','Kab. Aceh Besar',  12, 95.49000,  5.46610),
  list('Stasiun Meteorologi Kasiguncu','UPT','Sulawesi Tengah','Kab. Poso',  15,120.73000, -1.38000),
  list('Stasiun Klimatologi Maros','UPT','Sulawesi Selatan','Kab. Maros',  11,119.55000, -4.90000),
  list('Stasiun Geofisika Silaing Bawah','UPT','Sumatera Barat','Kota Padang Panjang', 773,100.40000,  0.45700),
  list('Stasiun Meteorologi Maritim','UPT','Sumatera Barat','Kota Padang',   2,100.37000, -0.99600),
  list('Stasiun Geofisika Waingapu','UPT','Nusa Tenggara Timur','Kab. Sumba Timur',  13,120.30000, -9.90000),
  list('Stasiun Meteorologi H. Asan','UPT','Kalimantan Tengah','Kab. Kotawaringin Timur',   3,112.93000, -2.55000),
  list('Stasiun Meteorologi Cut Bau Maimun Saleh','UPT','Nanggroe Aceh Darussalam','Kota Sabang', 126, 95.35000,  5.86979),
  list('Stasiun Meteorologi Sam Ratulangi','UPT','Sulawesi Utara','Kota Manado',  80,124.92000,  1.53000),
  list('Stasiun Meteorologi Juwata','UPT','Kalimantan Timur','Kota Tarakan',   6,117.57000,  3.33000),
  list('Stasiun Meteorologi Sudjarwo Tjondro Negoro','UPT','Papua','Kab. Kepulauan Yapen',   8,136.23000, -1.87000),
  list('Stasiun Klimatologi Kenten','UPT','Sumatera Selatan','Kota Palembang',   8,104.70000, -3.00000),
  list('Stasiun Geofisika Karang Panjang','UPT','Maluku','Kota Ambon',  83,128.19000, -3.69000),
  list('Stasiun Meteorologi Komodo','UPT','Nusa Tenggara Timur','Kab. Manggarai Barat',  67,119.89000, -8.49000),
  list('Stasiun Meteorologi Geser','UPT','Maluku','Kab. Seram Bagian Timur',   2,130.88000, -3.88130),
  list('Stasiun Klimatologi Semarang','UPT','Jawa Tengah','Kota Semarang', 227,110.42000, -6.98000),
  list('Stasiun Meteorologi Budiarto','UPT','Banten','Kab. Tangerang',  46,106.57000, -6.28700),
  list('Stasiun Meteorologi Soekarno Hatta','UPT','Banten','Kota Tangerang',   8,106.65000, -6.12000),
  list('Stasiun Meteorologi Maritim Tanjung Priok','UPT','DKI Jakarta','Kota Adm. Jakarta Utara',   2,106.88000, -6.10781),
  list('Stasiun Meteorologi Sultan Mahmud Badaruddin II','UPT','Sumatera Selatan','Kota Palembang',  11,104.70000, -2.90000),
  list('Stasiun Meteorologi Maritim','UPT','Sulawesi Tenggara','Kota Kendari',   2,122.60000, -3.95000),
  list('Stasiun Geofisika Saumlaki','UPT','Maluku','Kab Maluku Tenggara Barat',  26,131.30000, -7.98000),
  list('Stasiun Meteorologi Fransiskus Xaverius Seda','UPT','Nusa Tenggara Timur','Kab. Sikka',   4,122.25000, -8.63000),
  list('Stasiun Klimatologi Pondok Betung','UPT','Banten','Kota Tangerang Selata',  26,106.75000, -6.26107),
  list('Stasiun Meteorologi Aek Godang','UPT','Sumatera Utara','Kab. Padang Lawas Utara', 281, 99.45000,  1.55000),
  list('Stasiun Geofisika Sanglah','UPT','Bali','Kota Denpasar',  15,115.21000, -8.67689),
  list('Stasiun Meteorologi Umbu Mehang Kunda','UPT','Nusa Tenggara Timur','Kab. Sumba Timur',  52,120.33000, -9.67000),
  list('Stasiun Meteorologi Sangkapura','UPT','Jawa Timur','Kab. Gresik',NA_real_,112.63000, -5.85000),
  list('Stasiun Meteorologi Bandara Internasional Lombok','UPT','Nusa Tenggara Barat','Kota Mataram',  16,116.11000, -8.55000),
  list('Stasiun Meteorologi Seigun','UPT','Papua Barat','Kota Sorong',NA_real_,131.26000, -0.87726),
  list('Stasiun Meteorologi Tegal','UPT','Jawa Tengah','Kota Tegal',   3,109.15000, -6.85000),
  list('Stasiun Meteorologi Supadio','UPT','Kalimantan Barat','Kota Pontianak',   3,109.45000, -0.14206),
  list('Stasiun Meteorologi Sepinggan','UPT','Kalimantan Timur','Kota Balikpapan',   3,116.90000, -1.26000),
  list('Stasiun Meteorologi Ngurah Rai','UPT','Bali','Kab. Badung',   4,115.17000, -8.74916),
  list('Stasiun Meteorologi Djalaluddin','UPT','Gorontalo','Kab. Gorontalo',  18,122.83000,  0.63300),
  list('Stasiun Meteorologi Sultan Thaha','UPT','Jambi','Kota Jambi',  26,103.64000, -1.63368),
  list('Stasiun Geofisika Bandung','UPT','Jawa Barat','Kota Bandung', 829,107.60000, -6.88000),
  list('Stasiun Klimatologi Sicincin','UPT','Sumatera Barat','Kab. Padang Pariaman', 137,100.71000, -0.57500),
  list('Stasiun Meteorologi Enarotali','UPT','Papua','Kab. Paniai',1774,136.38000, -3.92688),
  list('Stasiun Geofisika Kepahiyang','UPT','Bengkulu','Kab. Kepahiang', 517,102.59000, -3.63000),
  list('Stasiun Klimatologi Banjar Baru','UPT','Kalimantan Selatan','Kota Banjarbaru',  55,114.84000, -3.46225),
  list('Stasiun Meteorologi Sultan Muhammad Salahuddin','UPT','Nusa Tenggara Barat','Kota Bima',   3,118.69000, -8.54323),
  list('Stasiun Meteorologi Hang Nadim','UPT','Kep. Riau','Kota Batam',  28,104.12000,  1.11700),
  list('Stasiun Geofisika Tanjung Pandan','UPT','Kep. Bangka Belitung','Kab. Belitung',  22,107.65000, -2.75750),
  list('Stasiun Meteorologi Tanjung Harapan','UPT','Kalimantan Timur','Kab. Bulungan',   3,117.22000,  2.50000),
  list('Stasiun Meteorologi F.L Tobing','UPT','Sumatera Utara','Kab. Tapanuli Tengah',  10, 98.88000,  1.55000),
  list('Stasiun Meteorologi Pangsuma','UPT','Kalimantan Barat','Kab. Kapuas Hulu',  43,112.93000,  0.84000),
  list('Stasiun Geofisika Tual','UPT','Maluku','Kab. Maluku Tenggara',  22,132.75000, -5.64000),
  list('Stasiun Geofisika Balikpapan','UPT','Kalimantan Timur','Kota Balikpapan',  31,116.91463, -1.25530),
  list('Stasiun Meteorologi Moanamani','UPT','Papua','Kab. Nabire',   3,135.52000, -3.35210),
  list('Stasiun Geofisika Banjarnegara','UPT','Jawa Tengah','Kab. Banjarnegara', 608,109.71000, -7.31800),
  list('Stasiun Geofisika Angkasapura','UPT','Papua','Kota Jayapura', 400,140.41000, -2.52000),
  list('Stasiun Meteorologi Beringin','UPT','Kalimantan Tengah','Kab. Barito Utara',  42,114.53000, -0.56000),
  list('Stasiun Meteorologi Andi Jemma','UPT','Sulawesi Selatan','Kab. Luwu Utara',  50,120.32000, -2.55440),
  list('Stasiun Meteorologi Sultan Iskandar Muda','UPT','Nanggroe Aceh Darussalam','Kab. Aceh Besar',  20, 95.42000,  5.52100),
  list('Stasiun Geofisika Kendari','UPT','Sulawesi Tenggara','Kota Kendari',  10,122.62000, -3.96000),
  list('Stasiun Meteorologi Tanah Merah','UPT','Papua','Kab. Boven Digoel',  97,140.31000, -6.10000),
  list('Stasiun Klimatologi Lasiana','UPT','Nusa Tenggara Timur','Kota Kupang',  19,123.60000,-10.17000),
  list('Stasiun Geofisika Sorong','UPT','Papua Barat','Kab. Sorong',  64,131.26000, -0.86000),
  list('Stasiun Meteorologi Utarom','UPT','Papua Barat','Kab. Kaimana',   5,133.75000, -3.67000),
  list('Stasiun Klimatologi Sei Durian','UPT','Jambi','Kab. Muaro Jambi',  34,103.49000, -1.60000),
  list('Stasiun Meteorologi Sangia Ni Bandera','UPT','Sulawesi Tenggara','Kab. Kolaka',  10,121.60000, -4.17500),
  list('Stasiun Meteorologi Pongtiku','UPT','Sulawesi Selatan','Kab. Tana Toraja', 879,119.82000, -3.04250),
  list('Stasiun Meteorologi Paloh','UPT','Kalimantan Barat','Kab. Sambas',  15,109.30000,  1.74000),
  list('Balai Besar Meteorologi, Klimatologi, dan Geofisika Wilayah V','UPT','Papua','Kab. Merauke',NA_real_,140.68000, -2.56700),
  list('Stasiun Meteorologi Cilacap','UPT','Jawa Tengah','Kab. Cilacap',  20,109.02000, -7.73000),
  list('Stasiun Meteorologi Dabo','UPT','Kep. Riau','Kab. Lingga',  29,104.58000, -0.48000),
  list('Stasiun Meteorologi Citeko','UPT','Jawa Barat','Kab. Bogor', 920,106.85000, -6.70000),
  list('Stasiun Meteorologi Japura','UPT','Riau','Kab. Indragiri Hulu',  19,102.32000,  0.33000),
  list('Stasiun Meteorologi Susilo','UPT','Kalimantan Barat','Kab. Sintang',  31,111.47000,  0.06000),
  list('Stasiun Meteorologi Mozez Kilangin','UPT','Papua','Kab. Mimika',  30,136.43000, -4.72000),
  list('Stasiun Meteorologi Raja Haji Abdullah Tanjung Balai Karimun','UPT','Kep. Riau','Kab. Karimun',   1,103.38000,  1.03000),
  list('Stasiun Meteorologi Dumatubun','UPT','Maluku','Kab. Maluku Tenggara',  16,132.73000, -5.66000),
  list('Stasiun Meteorologi Oesman Sadik','UPT','Maluku Utara','Kab. Halmahera Selatan',  14,127.50000, -0.63784),
  list('Stasiun Meteorologi H. Asan Hananjoedin','UPT','Kep. Bangka Belitung','Kab. Belitung',  47,107.75000, -2.75000),
  list('Stasiun Klimatologi Darmaga','UPT','Jawa Barat','Kota Bogor', 207,106.75000, -6.50000),
  list('Stasiun Geofisika Palu','UPT','Sulawesi Tengah','Kota Palu', 149,119.84000, -0.90519),
  list('Stasiun Klimatologi Masgar','UPT','Lampung','Kab. Pesawaran',  71,105.18000, -5.17236),
  list('Stasiun Meteorologi Naha','UPT','Sulawesi Utara','Kab. Kepulauan Sangihe',   6,125.51000,  3.68000),
  list('Stasiun Meteorologi Eltari','UPT','Nusa Tenggara Timur','Kota Kupang', 102,123.40000,-10.10000),
  list('Stasiun Meteorologi Jatiwangi','UPT','Jawa Barat','Kab. Majalengka',  85,108.27000, -6.75000),
  list('Pos Meteorologi Penggung','UPT','Jawa Barat','Kab. Cirebon',  23,108.33205, -6.45190),
  list('Stasiun Meteorologi Tarempa','UPT','Kep. Riau','Kab. Kepulauan Anambas',   2,106.25000,  3.20000),
  list('Stasiun Geofisika Gowa','UPT','Sulawesi Selatan','Kab. Gowa',  28,119.50000, -5.20000),
  list('Stasiun Klimatologi Negara','UPT','Bali','Kab. Jembrana',  24,114.62000, -8.34000),
  list('Stasiun Meteorologi Kemayoran','UPT','DKI Jakarta','Kota Adm. Jakarta Pusat',   4,106.84000, -6.15559),
  list('Stasiun Klimatologi Indrapuri','UPT','Nanggroe Aceh Darussalam','Kab. Aceh Besar',  51, 95.46000,  5.40000),
  list('Stasiun Geofisika Winangun','UPT','Sulawesi Utara','Kota Manado', 128,124.84000,  1.44390),
  list('Stasiun Geofisika Kota Bumi','UPT','Lampung','Kab. Lampung Utara',  60,104.87000, -4.83631),
  list('Stasiun Meteorologi Kijang','UPT','Kep. Riau','Kota Tanjung Pinang',  18,104.53000,  0.91816),
  list('Stasiun Meteorologi Mopah','UPT','Papua','Kab. Merauke',NA_real_,140.38000, -8.47000),
  list('Stasiun Meteorologi Perak I','UPT','Jawa Timur','Kota Surabaya',   3,112.46000, -7.13000),
  list('Stasiun Meteorologi Majene','UPT','Sulawesi Barat','Kab. Majene',  29,119.20000, -3.56194),
  list('Stasiun Geofisika Gunung Sitoli','UPT','Sumatera Utara','Kab. Nias', 175, 97.58000,  1.30000),
  list('Stasiun Meteorologi Sultan Syarif Kasim II','UPT','Riau','Kota Pekanbaru',  27,101.43000,  0.46700),
  list('Stasiun Meteorologi Kualanamu','UPT','Sumatera Utara','Kab. Deli Serdang',  23, 98.89000,  3.64200),
  list('Stasiun Meteorologi Mutiara Sis-Al Jufri','UPT','Sulawesi Tengah','Kota Palu', 228,119.91000, -0.91589),
  list('Stasiun Meteorologi Amahai','UPT','Maluku','Kab. Maluku Tengah',  10,128.88000, -3.35000),
  list('Stasiun Klimatologi Genyem','UPT','Papua','Kab. Jayapura',  70,140.16000, -2.60000),
  list('Stasiun Meteorologi Kalianget','UPT','Jawa Timur','Kab. Sumenep',NA_real_,113.97000, -7.05000),
  list('Stasiun Klimatologi Siantan','UPT','Kalimantan Barat','Kab. Pontianak',   2,109.19000,  0.07500),
  list('Stasiun Geofisika Yogyakarta','UPT','DI Yogyakarta','Kab. Sleman', 153,110.30000, -7.82000),
  list('Stasiun Meteorologi Nangapinoh','UPT','Kalimantan Barat','Kab. Melawi',  40,111.47000, -0.42000),
  list('Stasiun Meteorologi Beto Ambari','UPT','Sulawesi Tenggara','Kota Bau Bau', 233,122.62000, -5.47000),
  list('Stasiun Meteorologi Juanda','UPT','Jawa Timur','Kab. Sidoarjo',NA_real_,112.77000, -7.37000),
  list('Stasiun Meteorologi Temindung','UPT','Kalimantan Timur','Kota Samarinda',  10,117.16000, -0.48000),
  list('Stasiun Meteorologi Gusti Syamsir Alam','UPT','Kalimantan Selatan','Kab. Kotabaru',   2,116.17000, -3.30000),
  list('Stasiun Meteorologi Sultan Baabullah','UPT','Maluku Utara','Kota Ternate',  33,127.38000,  0.83330),
  list('Stasiun Geofisika Kampung Baru','UPT','Nusa Tenggara Timur','Kota Kupang',   4,123.61000,-10.15000),
  list('Stasiun Meteorologi Ahmad Yani','UPT','Jawa Tengah','Kota Semarang',NA_real_,110.38000, -6.98000),
  list('Stasiun Meteorologi Raden Inten II','UPT','Lampung','Kab. Lampung Selatan',  81,105.18000, -5.27000),
  list('Stasiun Meteorologi Maritim Pontianak','UPT','Kalimantan Barat','Kota Pontianak',   4,109.34000, -0.03000),
  list('Stasiun Geofisika Sawahan','UPT','Jawa Timur','Kab. Nganjuk', 835,111.79000, -7.74000),
  list('Stasiun Meteorologi Sultan Bantilan','UPT','Sulawesi Tengah','Kab. Toli Toli',   2,120.79000,  1.12114),
  list('Stasiun Meteorologi Dok II Jayapura','UPT','Papua','Kota Jayapura',   3,140.72000, -2.53000),
  list('Stasiun Geofisika Tretes','UPT','Jawa Timur','Kota Pasuruan', 832,112.64000, -7.70000),
  list('Stasiun Pemantau Atmosfir Global Bukit Koto Tabang','UPT','Sumatera Barat','Kab. Agam', 865,100.32000, -0.20000),
  list('Stasiun Meteorologi Ranai','UPT','Kep. Riau','Kab. Natuna',   2,108.39000,  3.91206),
  list('Stasiun Geofisika Ternate','UPT','Maluku Utara','Kota Ternate',  43,127.35000,  0.78100),
  list('Stasiun Meteorologi Iskandar','UPT','Kalimantan Tengah','Kab. Kotawaringin Barat',  22,111.66000, -2.73000),
  list('Stasiun Klimatologi Sampali','UPT','Sumatera Utara','Kab. Deli Serdang',  25, 98.78000,  3.61687),
  list('Stasiun Meteorologi Tardamu','UPT','Nusa Tenggara Timur','Kab. Kupang',  15,121.83000,-10.50000),
  list('Stasiun Meteorologi Banyuwangi','UPT','Jawa Timur','Kab. Banyuwangi',  43,114.35000, -8.21000),
  list('Stasiun Meteorologi Maritim Belawan','UPT','Sumatera Utara','Kota Medan',   3, 98.71000,  3.78824),
  list('Stasiun Klimatologi Kediri','UPT','Nusa Tenggara Barat','Kab. Lombok Barat',  52,116.17000, -8.63640)
))
names(df_sta) <- c('Nama','Jenis','Provinsi','Kabupaten','Elevasi','lon','lat')
df_sta$lon <- as.numeric(df_sta$lon)
df_sta$lat <- as.numeric(df_sta$lat)
df_sta$Elevasi <- as.numeric(df_sta$Elevasi)

pal <- colorFactor(
  palette = c("#2196F3", "#E53935", "#43A047"),
  levels  = c("Pos Hujan Kerjasama", "aaws", "UPT")
)

leaflet(df_sta) %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  addCircleMarkers(
    lng         = ~lon, lat = ~lat,
    radius      = 4, stroke = FALSE, fillOpacity = 0.6,
    color       = ~pal(Jenis),
    # Label muncul otomatis — visibilitasnya diatur JS berdasarkan zoom
    label       = ~Nama,
    labelOptions = labelOptions(
      noHide    = TRUE,       # selalu tampil (JS yang kontrol show/hide)
      textOnly  = TRUE,       # hanya teks, tanpa box
      direction = "top",
      offset    = c(0, -6),
      style     = list(
        "font-size"   = "10px",
        "font-weight" = "600",
        "color"       = "#1A3A5C",
        "text-shadow" = "0 0 3px #fff, 0 0 3px #fff"
      )
    ),
    popup = ~paste0(
      "<b>", Nama, "</b><br>",
      "Jenis: ", Jenis, "<br>",
      "Provinsi: ", Provinsi, "<br>",
      "Elevasi: ", ifelse(is.na(Elevasi), "&mdash;", paste0(Elevasi, " mdpl"))
    )
  ) %>%
  addLegend("bottomright",
    pal     = pal,
    values  = ~Jenis,
    title   = "Jenis Stasiun",
    opacity = 0.8,
    labFormat = labelFormat(
      transform = function(x) c(
        "Pos Hujan Kerjasama (5.918)",
        "AAWS (391)",
        "UPT (178)"
      )[match(x, c("Pos Hujan Kerjasama","aaws","UPT"))]
    )
  ) %>%
  # Label hanya muncul saat zoom >= 9
  htmlwidgets::onRender("
    function(el, x) {
      var map = this;
      var ZOOM_LABEL = 9;

      function updateLabels() {
        var zoom = map.getZoom();
        // Semua tooltip leaflet ada di class .leaflet-tooltip
        var labels = el.querySelectorAll('.leaflet-tooltip');
        labels.forEach(function(lbl) {
          lbl.style.display = zoom >= ZOOM_LABEL ? 'block' : 'none';
        });
      }

      // Jalankan saat zoom berubah
      map.on('zoomend', updateLabels);

      // Jalankan saat peta selesai render (initial state)
      map.on('layeradd', function() {
        setTimeout(updateLabels, 100);
      });

      // Pastikan state awal benar
      setTimeout(updateLabels, 500);
    }
  ")

Gambar 1: Sebaran stasiun hujan BMKG di Indonesia (sampel 450 dari 6.487 stasiun). Biru = Pos Hujan Kerjasama, merah = AAWS, hijau = UPT.

2.2 Kerapatan Stasiun per Wilayah

Standar WMO (2008) merekomendasikan kerapatan minimum 1 pos per 250 km² untuk kawasan pegunungan tropis.

Tampilkan kode R

luas_pulau <- data.frame(
  Wilayah   = c("Jawa","Sumatera","Kalimantan","Sulawesi","Maluku","Papua"),
  Luas_km2  = c(128297, 473481, 748168, 188522, 31982, 421981),
  N_Stasiun = c(2765, 1321, 842, 1023, 213, 323),
  stringsAsFactors = FALSE
) %>%
  mutate(
    Kerapatan = round(Luas_km2 / N_Stasiun),
    Status    = case_when(
      Kerapatan <= 250  ~ "Memenuhi WMO",
      Kerapatan <= 500  ~ "Di bawah standar",
      TRUE              ~ "Jauh di bawah standar"
    )
  )

datatable(
  luas_pulau,
  colnames = c("Wilayah","Luas (km²)","Jumlah Stasiun",
               "Kerapatan (km²/pos)","Status"),
  options  = list(dom="t", pageLength=10),
  rownames = FALSE
) %>%
  formatStyle("Status",
    backgroundColor = styleEqual(
      c("Memenuhi WMO","Di bawah standar","Jauh di bawah standar"),
      c("#C8E6C9","#FFF9C4","#FFCDD2")
    )
  )

Tabel 1: Perbandingan kerapatan pos hujan per wilayah vs standar WMO

2.3 Rekaman Data yang Pendek dan Tidak Kontinu

Implikasi untuk analisis AMDAL

Untuk analisis banjir kala ulang Q25, Q50, Q100 yang dibutuhkan dalam perencanaan infrastruktur, data dengan rekaman < 15 tahun menghasilkan ketidakpastian yang tinggi. Penggunaan GPM IMERG (tersedia sejak 2000, ~25 tahun) dapat meningkatkan keandalan analisis secara signifikan.

3 GPM IMERG sebagai Alternatif

Global Precipitation Measurement (GPM) IMERG adalah produk curah hujan berbasis satelit hasil kolaborasi NASA dan JAXA. Tersedia dalam tiga varian: Early (latency 4 jam), Late (14 jam), dan Final (3.5 bulan, sudah dikoreksi dengan data pos hujan).

Parameter	Nilai
Resolusi spasial	0.1° × 0.1° (~11 km)
Resolusi temporal	30 menit
Ketersediaan data	Juni 2000 – sekarang
Cakupan	60°N – 60°S (seluruh Indonesia)
Akses	Gratis via NASA Earthdata

3.1 Download GPM IMERG dengan R

Tampilkan kode R

library(httr)
library(lubridate)

download_imerg_daily <- function(tanggal, user, pass,
                                  output_dir = "data/imerg/") {
  dir.create(output_dir, showWarnings = FALSE, recursive = TRUE)
  thn  <- format(tanggal, "%Y")
  doy  <- format(tanggal, "%j")
  tgl  <- format(tanggal, "%Y%m%d")
  nama_file <- paste0("3B-DAY.MS.MRG.3IMERG.", tgl, ".V07B.nc4")
  output    <- file.path(output_dir, nama_file)

  if (file.exists(output)) {
    message("Sudah ada: ", nama_file); return(invisible(output))
  }

  resp <- GET(
    url = paste0("https://gpm.nasa.gov/data/imerg/final/",
                 thn, "/", doy, "/", nama_file),
    authenticate(user, pass),
    write_disk(output, overwrite = FALSE),
    progress()
  )

  if (status_code(resp) == 200) message("OK: ", nama_file) else {
    message("Gagal (", status_code(resp), ")")
    file.remove(output)
  }
  invisible(output)
}

# Download 1 tahun untuk DAS target
for (tgl in seq(as.Date("2020-01-01"), as.Date("2020-12-31"), by="day")) {
  download_imerg_daily(as.Date(tgl, origin="1970-01-01"),
                       user="NASA_USER", pass="NASA_PASS")
  Sys.sleep(0.5)
}

4 Validasi GPM IMERG

4.1 Metrik Statistik

Tampilkan kode R

hitung_validasi <- function(obs, sim) {
  mean_obs <- mean(obs, na.rm=TRUE)
  mean_sim <- mean(sim, na.rm=TRUE)
  nse      <- 1 - sum((obs-sim)^2,na.rm=TRUE) /
                  sum((obs-mean_obs)^2,na.rm=TRUE)
  pbias    <- (sum(sim-obs,na.rm=TRUE)/sum(obs,na.rm=TRUE))*100
  rmse     <- sqrt(mean((obs-sim)^2,na.rm=TRUE))
  r        <- cor(obs, sim, use="complete.obs")
  alpha    <- sd(sim,na.rm=TRUE)/sd(obs,na.rm=TRUE)
  beta     <- mean_sim/mean_obs
  kge      <- 1 - sqrt((r-1)^2+(alpha-1)^2+(beta-1)^2)
  list(NSE=round(nse,3), PBIAS=round(pbias,1),
       RMSE=round(rmse,1), KGE=round(kge,3))
}

Grafik Interaktif — Koreksi Bias

Geser slider di bawah untuk mensimulasikan koreksi bias pada data GPM IMERG. Nilai metrik validasi (NSE, KGE, PBIAS, RMSE) akan berubah otomatis.

Tampilkan kode R

imerg_kor = imerg_r.map(v => Math.round(v * bias * 10) / 10)

function hitung(obs, sim) {
  const n   = obs.length
  const mo  = obs.reduce((a,b)=>a+b,0)/n
  const ms  = sim.reduce((a,b)=>a+b,0)/n
  const ssTot = obs.reduce((s,o)=>s+Math.pow(o-mo,2),0)
  const ssRes = obs.reduce((s,o,i)=>s+Math.pow(o-sim[i],2),0)
  const nse   = 1 - ssRes/ssTot
  const pbias = (sim.reduce((a,b)=>a+b,0)-obs.reduce((a,b)=>a+b,0))/
                obs.reduce((a,b)=>a+b,0)*100
  const rmse  = Math.sqrt(ssRes/n)
  const num   = obs.reduce((s,o,i)=>s+(o-mo)*(sim[i]-ms),0)
  const d1    = Math.sqrt(obs.reduce((s,o)=>s+Math.pow(o-mo,2),0))
  const d2    = Math.sqrt(sim.reduce((s,v)=>s+Math.pow(v-ms,2),0))
  const r     = num/(d1*d2)
  const alpha = d2/Math.sqrt(n) / (d1/Math.sqrt(n))
  const kge   = 1-Math.sqrt(Math.pow(r-1,2)+Math.pow(alpha-1,2)+Math.pow(ms/mo-1,2))
  return {nse,kge,pbias,rmse}
}

mtk = hitung(obs_r, imerg_kor)

function sta(v, t) {
  if (t==="nse")   return v>=0.75?"🟢 Sangat baik":v>=0.65?"🟡 Baik":v>=0.5?"🟠 Memuaskan":"🔴 Tidak memuaskan"
  if (t==="kge")   return v>=0.75?"🟢 Sangat baik":v>=0.5?"🟡 Baik":"🔴 Perlu koreksi"
  if (t==="pbias") return Math.abs(v)<=10?"🟢 Sangat baik":Math.abs(v)<=15?"🟡 Baik":Math.abs(v)<=25?"🟠 Memuaskan":"🔴 Tidak memuaskan"
  return "—"
}

// Grafik Observable Plot
Plot.plot({
  height: 300,
  marginBottom: 40,
  marginLeft: 55,
  x: {
    label: null,
    tickFormat: (_, i) => i % 4 === 0 ? bulan_r[i] : "",
    tickSize: 0
  },
  y: { label: "Curah Hujan (mm/bulan)", grid: true },
  marks: [
    Plot.areaY(
      obs_r.map((v, i) => ({x: i, y: v})),
      {x: "x", y: "y", fill: "#1565C0", fillOpacity: 0.12}
    ),
    Plot.lineY(
      obs_r.map((v, i) => ({x: i, y: v})),
      {x: "x", y: "y", stroke: "#1565C0", strokeWidth: 2,
       tip: true, title: (_, i) => `${bulan_r[i]}\nObservasi: ${obs_r[i]} mm`}
    ),
    Plot.lineY(
      imerg_kor.map((v, i) => ({x: i, y: v})),
      {x: "x", y: "y", stroke: "#C62828", strokeWidth: 2,
       strokeDasharray: "6,4",
       tip: true, title: (_, i) => `${bulan_r[i]}\nGPM IMERG: ${imerg_kor[i]} mm`}
    )
  ],
  style: {fontSize: 11}
})

Tampilkan kode R

// Tabel metrik reaktif
html`<div style="display:flex;gap:20px;flex-wrap:wrap;align-items:flex-start;margin-top:0.5rem;">

  <table style="border-collapse:collapse;font-size:0.88rem;flex:1;min-width:340px;">
    <thead>
      <tr style="background:#1A3A5C;color:white;">
        <th style="padding:7px 12px;text-align:left;">Metrik</th>
        <th style="padding:7px 12px;text-align:center;">Nilai</th>
        <th style="padding:7px 12px;text-align:left;">Kriteria Baik</th>
        <th style="padding:7px 12px;text-align:left;">Status</th>
      </tr>
    </thead>
    <tbody>
      <tr>
        <td style="padding:6px 12px;font-weight:600;">NSE</td>
        <td style="padding:6px 12px;text-align:center;font-family:monospace;font-weight:700;font-size:1rem;">
          ${mtk.nse.toFixed(3)}</td>
        <td style="padding:6px 12px;color:#555;">&gt; 0.65</td>
        <td style="padding:6px 12px;">${sta(mtk.nse,"nse")}</td>
      </tr>
      <tr style="background:#f8f8f8;">
        <td style="padding:6px 12px;font-weight:600;">KGE</td>
        <td style="padding:6px 12px;text-align:center;font-family:monospace;font-weight:700;font-size:1rem;">
          ${mtk.kge.toFixed(3)}</td>
        <td style="padding:6px 12px;color:#555;">&gt; 0.50</td>
        <td style="padding:6px 12px;">${sta(mtk.kge,"kge")}</td>
      </tr>
      <tr>
        <td style="padding:6px 12px;font-weight:600;">PBIAS (%)</td>
        <td style="padding:6px 12px;text-align:center;font-family:monospace;font-weight:700;font-size:1rem;">
          ${mtk.pbias.toFixed(1)}</td>
        <td style="padding:6px 12px;color:#555;">± 10–15%</td>
        <td style="padding:6px 12px;">${sta(mtk.pbias,"pbias")}</td>
      </tr>
      <tr style="background:#f8f8f8;">
        <td style="padding:6px 12px;font-weight:600;">RMSE (mm)</td>
        <td style="padding:6px 12px;text-align:center;font-family:monospace;font-weight:700;font-size:1rem;">
          ${mtk.rmse.toFixed(1)}</td>
        <td style="padding:6px 12px;color:#555;">Sekecil mungkin</td>
        <td style="padding:6px 12px;">—</td>
      </tr>
    </tbody>
  </table>

  <div style="background:#EFF6FF;border-left:3px solid #1565C0;padding:0.75rem 1rem;
              border-radius:0 6px 6px 0;font-size:0.83rem;min-width:200px;max-width:260px;">
    <strong style="color:#1A3A5C;">Kriteria Moriasi et al. (2007):</strong>
    <table style="margin-top:0.4rem;font-size:0.82rem;border-collapse:collapse;">
      <tr><td style="padding:2px 6px 2px 0;color:#555;">NSE &gt; 0.75</td><td>Sangat baik</td></tr>
      <tr><td style="padding:2px 6px 2px 0;color:#555;">NSE 0.65–0.75</td><td>Baik</td></tr>
      <tr><td style="padding:2px 6px 2px 0;color:#555;">NSE 0.50–0.65</td><td>Memuaskan</td></tr>
      <tr><td style="padding:2px 6px 2px 0;color:#555;">NSE &lt; 0.50</td><td>Tidak memuaskan</td></tr>
    </table>
  </div>

</div>`

Sumber data: DAS Akelamo, Pulau Obi, 2020–2022 (data simulasi untuk ilustrasi)

5 Kesimpulan

Kerapatan pos hujan di Indonesia sangat tidak merata — hanya Jawa yang mendekati standar WMO, sementara Kalimantan, Papua, dan Maluku jauh di bawah.
GPM IMERG Final Run merupakan alternatif yang andal untuk DAS tanpa pos hujan memadai, dengan catatan perlu divalidasi dan dikoreksi bias terlebih dahulu.
Metrik NSE, KGE, dan PBIAS dapat digunakan sebagai ukuran kelayakan sesuai panduan Moriasi dkk. (2007).

Unduh data & skrip

Skrip R lengkap tersedia di irpanchumaedi.com. Data stasiun: inventarisasi BMKG 2017, 6.487 titik seluruh Indonesia.

Referensi

Moriasi, D. N., J. G. Arnold, M. W. Van Liew, R. L. Bingner, R. D. Harmel, dan T. L. Veith. 2007. “Model Evaluation Guidelines for Systematic Quantification of Accuracy in Watershed Simulations.” Transactions of the ASABE 50 (3): 885–900. https://doi.org/10.13031/2013.23153.