Tuesday, August 2, 2016

Research Note #5 - Algoritma Array Koordinat Resolusi 0.1 Deg (GSMaP)

Algoritma ini digunakan untuk membuat array koordinat (geo) dengan resolusi spasial 0.1 derajat yang biasa digunakan oleh data GSMaP. Pada data ini, terdapat 3600 komponen x dan 1200 komponen y yang masing-masing merepresentasikan bujur dan lintang.

Pada tulisan ini, bahasa pemrograman yang digunakan adalah Fortran, walaupun algoritma ini sejatinya bisa digunakan pada semua bahasa pemrograman.


      integer :: x,y
      real :: o_lat,o_lon,delta
      delta=0.1

      do y=1,1200
        do x=1,3600
          o_lon=0.05+delta*(x-1)

c---- Bila bujur melewati dateline (>180E) ----- 

          if (o_lon>180) then
            o_lon=-360+(0.05+delta*(x-1))
          end if
c-----------------------------------------------

          o_lat=59.95-delta*(y-1)
          print *,y,'',x,'',o_lat,'',o_lon
        end do
      end do

Hasilnya adalah sebagai berikut :

Gambar pertama adalah output pada awal loop, gambar kedua adalah output ketika bujur melewati dateline (180E) dan yang terakhir adalah output pada akhir loop.

1 comment: