Zonal Statistics or Extracting Raster Data - Revision history

Ianm: Set PropertyValue: TargetDate = 2015-01-20

2015-01-21T04:23:28Z

Set PropertyValue: TargetDate = 2015-01-20

Ianm: Set PropertyValue: StartDate = 2015-01-20

2015-01-21T04:23:15Z

Set PropertyValue: StartDate = 2015-01-20

Ianm: Set PropertyValue: Progress = 100

2015-01-21T04:23:10Z

Set PropertyValue: Progress = 100

Ianm at 04:22, 21 January 2015

2015-01-21T04:22:55Z

Ianm at 04:22, 21 January 2015

2015-01-21T04:22:28Z

Ianm at 04:21, 21 January 2015

2015-01-21T04:21:47Z

Ianm at 04:20, 21 January 2015

2015-01-21T04:20:37Z

Ianm at 04:19, 21 January 2015

2015-01-21T04:19:46Z

Ianm at 04:19, 21 January 2015

2015-01-21T04:19:17Z

Ianm at 04:18, 21 January 2015

2015-01-21T04:18:31Z

@@ Line 84: / Line 84: @@
 	Progress=100|
 	StartDate=2015-01-20|
 	Type=Low}}

@@ Line 83: / Line 83: @@
 	Participants=Ian_McCullough|
 	Progress=100|
 	Type=Low}}

@@ Line 77: / Line 77: @@
 <!-- Do NOT Edit below this Line -->
 {{#set:
 	Expertise=R_programming|
 	Expertise=Raster|
 	Owner=Ian_McCullough|
 	Participants=Ian_McCullough|
 	Type=Low}}

@@ Line 73: / Line 73: @@
 write.csv(bus, "H:/Ian_GIS/Tejon/FullStudyArea/Mar_precip_SA.csv", na="NA")
-</nowiki>
 <!-- Add any wiki Text above this Line -->
 <!-- Do NOT Edit below this Line -->

@@ Line 49: / Line 49: @@
 library(maptools)
 points = readShapePoints("H:/Ian_GIS/Tejon_data/TejonPlots/TejonPoints")
 studyarea = readShapePoly("H:/Ian_GIS/Tejon_data/FullStudyArea/Tejon_study_area",IDvar="ID")

@@ Line 34: / Line 34: @@
 library(raster)
-# Make list of asci files in already specified working directory
+<nowiki>#</nowiki> Make list of asci files in already specified working directory
 Clim_list = list.files(pattern = "mar.asc$") # "mar.asc$" queries files in directory ending in this pattern
-# Check number of files in raster list
+<nowiki>#</nowiki> Check number of files in raster list
 length = length(Clim_list)
-# Crete stack of rasters in list
+<nowiki>#</nowiki> Crete stack of rasters in list
 Clim_stack = stack(Clim_list)
-# Create spatial dataset from shapefile (use points, polygons or lines)
+<nowiki>#</nowiki> Create spatial dataset from shapefile (use points, polygons or lines)
-library(maptools) # (may have to install)
 points = readShapePoints("H:/Ian_GIS/Tejon_data/TejonPlots/TejonPoints")
 studyarea = readShapePoly("H:/Ian_GIS/Tejon_data/FullStudyArea/Tejon_study_area",IDvar="ID")
-# Extract raster pixels from point dataset
+<nowiki>#</nowiki> Extract raster pixels from point dataset
 for(j in 1:length) {
    ex = extract(Clim_stack,points,method='simple',layer=1,nl=length, df=T)
    ex$ID = NULL}
-# Extract study area mean from polygon dataset
+<nowiki>#</nowiki> Extract study area mean from polygon dataset
 for(j in 1:length) {
    bus = extract(Clim_stack,studyarea, fun=mean, layer=1, nl=length, df=T, na.rm=T)
@@ Line 59: / Line 63: @@
 bus = t(bus)
-# Output csv file (point data)
+<nowiki>#</nowiki> Output csv file (point data)
 write.csv(ex, "H:/Ian_GIS/Tejon/TejonPlots/Mar_precip_pts.csv", row.names=plots$Id, na="NA")
-# Output csv file (zonal polygon data)
+<nowiki>#</nowiki> Output csv file (zonal polygon data)
 write.csv(bus, "H:/Ian_GIS/Tejon/FullStudyArea/Mar_precip_SA.csv", na="NA")

@@ Line 27: / Line 27: @@
 <nowiki>#</nowiki> Update working directory
 setwd("H:/Ian_GIS/Climate/precip/Tejon/Historic")
 <nowiki>#</nowiki> Load raster package (may have to install first)
 library(raster)

← Older revision		Revision as of 04:19, 21 January 2015
Line 24:		Line 24:


−	# Import a shapefile to R and use it to extract zonal statistics from a directory of rasters	+	<nowiki>#</nowiki> Import a shapefile to R and use it to extract zonal statistics from a directory of rasters

	# Update working directory		# Update working directory

← Older revision		Revision as of 04:18, 21 January 2015
Line 24:		Line 24:


−	~~<nowiki>~~# Import a shapefile to R and use it to extract zonal statistics from a directory of rasters	+	# Import a shapefile to R and use it to extract zonal statistics from a directory of rasters

	# Update working directory		# Update working directory