Viac

Zistite, ktoré hodnoty zemepisnej šírky a dĺžky klesajú v ktorých rastrových bunkách


Som úplný začiatočník v analýze geografických údajov s R.

Mám databázu s dĺžkou a šírkou asi milióna budov v krajine. Chcel by som zistiť priemernú hustotu týchto budov. Za to som pomocou Rasteru rozdelil celú krajinu na niekoľko buniek a použil projekciu s príslušným profilom proj4. Teraz chcem skontrolovať, ktoré budovy prichádzajú pod ktorou bunkou, aby sme našli hustotu.

> metersCoordinates <- totalMeterDatabase [, c ("longitude", "latitude")]> Coordinates (metersCoordinates) <- ~ longitude + latitude> r <- raster (ncols = 6000, nrows = 2000)> r [] <- 0 > crs (r) <- "+ proj = somerc + lat_0 = 46,95240555555556 + lon_0 = 7,439583333333333 + k_0 = 1 + x_0 = 6000 + y_0 = 2000 + ellps = bessel + towgs84 = 674,374,15,056,405,346,0,0,0,0 , 0 + jednotky = m + no_defs "> projekcia (r)

Môžete mi prosím povedať, ako to nájdem?


Ak vašej otázke rozumiem, môžete použiť funkciu extract (). Najskôr pri vytváraní rastra dajte každej rastrovej bunke jedinečnú hodnotu. Tiež by som odporučil dať rastru rozsah vašich súradníc. Takže:

metersCoordinates <- totalMeterDatabase [, c ("longitude", "latitude")] Coordinates (metersCoordinates) <- ~ longitude + latitude r <- raster (matrix (seq (1,6000 * 2000), ncol = 6000, nrow = 2000) ), xmx = max (meterCoordinates $ longitude), xmn = min (metrovCoordinates $ longitude), ymx = max (metrovCoordinates $ latitude), ymn = min (metrovCoordinates $ latitude)) totalMeterDatabase $ cellnumber <- výpis (r, metersCoordinates)

A pomocou toho urobte svoju agregáciu / analýzu.