def test_semivariogram_calculation(self):
     V = Variogram()
     x = np.random.random(100)
     lon = np.random.random(100)
     lat = np.random.random(100)
     lags = np.asarray([1., 2., 3., 4.])
     V.semivariogram(x, lon, lat, lags, 1.)
Exemple #2
0
 def test_semivariogram_calculation(self):
     V = Variogram()
     x = np.random.random(100)
     lon = np.random.random(100)
     lat = np.random.random(100)
     lags = np.asarray([1.,2.,3.,4.])
     V.semivariogram(x, lon, lat, lags, 1.)
    def test_variogram_semivariance(self):

        V = Variogram()

        x = np.random.random(100)
        lon = np.random.random(100) * 10. - 90.
        lat = np.random.random(100) * 10. - 90.

        h_km = 20.
        dh_km = 2.

        # just test if it works; no reference solution yet
        g = V._semivariance(x, lon, lat, np.asarray([h_km]), dh_km)
Exemple #4
0
    def test_variogram_semivariance(self):

        V = Variogram()

        x = np.random.random(100)
        lon = np.random.random(100)*10.-90.
        lat = np.random.random(100)*10.-90.

        h_km = 20.
        dh_km = 2.

        # just test if it works; no reference solution yet
        g = V._semivariance(x, lon, lat, np.asarray([h_km]), dh_km)
Exemple #5
0
    def test_variogram_orthodrome(self):
        # example distance Berlin-Tokio
        # http://de.wikipedia.org/wiki/Orthodrome

        lat_berlin = 52.517
        lon_berlin = 13.4
        lat_tokio = 35.70
        lon_tokio = 139.767

        V = Variogram()
        r1 = V._orthodrome(lon_tokio, lat_tokio, lon_berlin, lat_berlin, radius=6370.*1000.)
        self.assertTrue(abs(r1-8918000.)<1000.)
        r2 = V._orthodrome(lon_berlin, lat_berlin, lon_tokio, lat_tokio, radius=6370.*1000.)
        self.assertTrue(abs(r2-8918000.)<1000.)
    def test_variogram_orthodrome(self):
        # example distance Berlin-Tokio
        # http://de.wikipedia.org/wiki/Orthodrome

        lat_berlin = 52.517
        lon_berlin = 13.4
        lat_tokio = 35.70
        lon_tokio = 139.767

        V = Variogram()
        r1 = V._orthodrome(lon_tokio,
                           lat_tokio,
                           lon_berlin,
                           lat_berlin,
                           radius=6370. * 1000.)
        self.assertTrue(abs(r1 - 8918000.) < 1000.)
        r2 = V._orthodrome(lon_berlin,
                           lat_berlin,
                           lon_tokio,
                           lat_tokio,
                           radius=6370. * 1000.)
        self.assertTrue(abs(r2 - 8918000.) < 1000.)