In [1]:
setwd('/Users/erinlarson/Dropbox/Documents/GIT/DDESSS_Stream_Dream_Team')

library(dplyr)


Attaching package: ‘dplyr’

The following objects are masked from ‘package:stats’:

    filter, lag

The following objects are masked from ‘package:base’:

    intersect, setdiff, setequal, union


In [7]:
data <- read.csv("Lat_Long_P_LandUse_StreamOrder.csv", head=T)
t(data[1,])


1
ActivityStartDate2/7/01
mo2
day7
date2001-02-07T00:00:00Z
MonitoringLocationIdentifierUSGS-09498400
HydrologicConditionStable, normal stage
HydrologicEventRoutine sample
ResultDetectionConditionTextNot Detected
CharacteristicNamePhosphorus
ResultSampleFractionTextTotal
ResultMeasureValueNA
ResultMeasure20
ResultMeasure_MeasureUnitCode
ResultStatusIdentifierAccepted
HUCEightDigitCode15060103
DrainageAreaMeasure_MeasureValue195
DrainageAreaMeasure_MeasureUnitCodesq mi
HUC815060103
PFOR273.918
PWETL20.01
PURB20.132
PAGT20.076
OrganizationIdentifierUSGS-AZ
OrganizationFormalNameUSGS Arizona Water Science Center
MonitoringLocationIdentifier.1USGS-09498400
MonitoringLocationNamePINAL CREEK AT INSPIRATION DAM, NR GLOBE, AZ.
MonitoringLocationTypeNameStream
MonitoringLocationDescriptionTextNA
HUCEightDigitCode.115060103
DrainageAreaMeasure_MeasureValue.1195
DrainageAreaMeasure_MeasureUnitCode.1sq mi
ContributingDrainageAreaMeasure_MeasureValueNA
ContributingDrainageAreaMeasure_MeasureUnitCode
LatitudeMeasure33.57311
LongitudeMeasure-110.9012
SourceMapScaleNumericNA
HorizontalAccuracyMeasure.MeasureValue10
HorizontalAccuracyMeasure.MeasureUnitCodeseconds
HorizontalCollectionMethodNameInterpolated from MAP.
HorizontalCoordinateReferenceSystemDatumNameNAD83
VerticalMeasure.MeasureValue2740
VerticalMeasure.MeasureUnitCodefeet
VerticalAccuracyMeasure.MeasureValue20
VerticalAccuracyMeasure.MeasureUnitCodefeet
VerticalCollectionMethodNameInterpolated from topographic map.
VerticalCoordinateReferenceSystemDatumNameNGVD29
CountryCodeUS
StateCode4
CountyCode7
AquiferNameNA
FormationTypeTextNA
AquiferTypeNameNA
ConstructionDateTextNA
WellDepthMeasure.MeasureValueNA
WellDepthMeasure.MeasureUnitCodeNA
WellHoleDepthMeasure.MeasureValueNA
WellHoleDepthMeasure.MeasureUnitCodeNA
ProviderNameNWIS
HUC8.115060103
STRAHLER22.4

In [9]:
group_by(data, MonitoringLocationIdentifier) %>% 
    summarise(AvgLat = mean(LatitudeMeasure, na.rm = TRUE), 
              AvgLong = mean(LongitudeMeasure, na.rm = TRUE)
             ) %>% write.csv("SiteLatLong.csv")