Displaying 2 results from an estimated 2 matches for "gis_data".
Did you mean:
gid_data
2012 Apr 04
0
crosstabs and histograms with flexible binning of dates
...0 0
Below is my code to handle the arbitrary bins, but I'm guessing there are
useful libraries and more elegant approaches. Any pointers would be
appreciated.
import(foreign)
# LOAD FILE
#parcels=read.dbf() #depending on source file
parcels=read.delim("~/Projects/GIS_DATA/Parcels_NSP_BlockGroup.txt")
attach(parcels)
# DEFINE BINNING
basedate=as.Date("2011/05/11")
currentdate=basedate
interval=3 #width of interval in months. 3 = quarterly
num_intervals=5 #how many intervals to include after basedate
for (i in c(1:num_intervals)) {
startdate=current...
2012 May 31
0
Incorporate a shapefile with an package animation
...--------------------------------------------------------------------------------------------------------------------------------
### Animation of EDEN real time uncorrected data.
library(animation)
library(ncdf)
library(fields)
library(maptools)
ENP_WCA <-
readShapePoly(system.file("U:\\GIS_Data\\GISlayers\\Boundaries_Lines\\EVERareas_Dissolve.shp",
package="maptools")[1], proj4String=CRS("+proj=utm +zone17
+datum=WGS84"))
# here an error is reported.
(Error in getinfo.shape(filen) : Error opening SHP file
# change to the working directory where the N...