Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Hands on training on developing Hands on training on developing
ground water level map using ground water level map using geo-statistical analystgeo-statistical analyst
Dr. A.K.M. Saiful IslamInstitute of Water and Flood Management (IWFM)Bangladesh University of Engineering and Technology (BUET)
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Geo-statistical Analyst of ArcGISGeo-statistical Analyst of ArcGIS
This training will be on:
1.1. Represent data Represent data 2.2. Explore dataExplore data3.3. Fit a interpolation Model Fit a interpolation Model 4.4. Diagnosis outputDiagnosis output5.5. Create ground water level mapsCreate ground water level maps
Input Data
Groundwater well data of Dinajpur district of Bangladesh
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Study Area and Data
• Study area – Seven upazillas of Dinajpur
District of Bangladesh
• Data– Data from 27 Groundwater
observation Wells as shape file “gwowell_bwdb.shp”. Weekly data from December to May for 1994 to 2003
– Upazilla shape file “upazila.shp”
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Activate Geo-statistical Analyst
• Turn on Geostatistaical Anaylst of ArcGIS
Enable toolbarEnable Extension
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Add Data
• Add both shape files
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am 1. Represent Data
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Groundwater well data
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am 2. Explore Data
a) Histogram
b) Normal Q-Q Plot
c) Trend Analysis
d) Voronoi Map
e) Semivariogram
f) Covariance cloud
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am a) Histogram
• Select attribute: any data e.g. DEC05_1994 • We can change no of bars or bin size• Distribution is normal
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Transformation
• Log- transformation doesn’t change distribution pattern
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am b) Normal Q-Q Plot
• Normal Q-Q plot is straight line which represents normal distribution
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am
c) Trend Analysis
• Shows trend in both X and Y direction since the projection lines (blue and green) are not straight.
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am d) Voronoi map
• Shows the zone of influence of known data points
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am e) Semi-variogram
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Shows search directional
• Exhibits directional influence in different angle (arrows)
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am f) Co-variance
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am 3. Fit interpolation model
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Kriging Geo-statistical method
• Select ordinary kriging
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Semivariogram modeling
• Select spherical method
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Searching neighbour
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Cross validation
• Root mean square error is 1.437
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Report of output layer
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Prediction map
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Extent of Map
• Set Extend as upazilla
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Export as Raster
• Select cell size as 100
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Prediction map as raster
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Zonal statistics
• Zonal statistics from Spatial Analyst
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Mean ground water level
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am
Change color for Mean ground water level of Dinajpur
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am
Thank you !
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Glossary
Methods of Interpolation in geostatisticsal Analysis o Inverse Distance Weighting (IDW)o Global Polynomial (GP)o Local Polynomial (LP)o Radial Basis Functions (RBF)o Krigingo Cokriging
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Inverse Distance Weighting (IDW)
• Inverse Distance Weighting (IDW) is a quick deterministic interpolator that is exact. There are very few decisions to make regarding model parameters. It can be a good way to take a first look at an interpolated surface. However, there is no assessment of prediction errors, and IDW can produce "bulls eyes" around data locations. There are no assumptions required of the data.
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Global Polynomial (GP)
• Global Polynomial (GP) is a quick deterministic interpolator that is smooth (inexact). There are very few decisions to make regarding model parameters. It is best used for surfaces that change slowly and gradually. However, there is no assessment of prediction errors and it may be too smooth. Locations at the edge of the data can have a large effect on the surface. There are no assumptions required of the data.
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Local Polynomial (LP)
• Local Polynomial (LP) is a moderately quick deterministic interpolator that is smooth (inexact). It is more flexible than the global polynomial method, but there are more parameter decisions. There is no assessment of prediction errors. The method provides prediction surfaces that are comparable to kriging with measurement errors. Local polynomial methods do not allow you to investigate the autocorrelation of the data, making it less flexible and more automatic than kriging. There are no assumptions required of the data.
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Radial Basis Functions (RBF)
• Radial Basis Functions (RBF) are moderately quick deterministic interpolators that are exact. They are much more flexible than IDW, but there are more parameter decisions. There is no assessment of prediction errors. The method provides prediction surfaces that are comparable to the exact form of kriging. Radial Basis Functions do not allow you to investigate the autocorrelation of the data, making it less flexible and more automatic than kriging. Radial Basis Functions make no assumptions about the data.
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Kriging
• Kriging is a moderately quick interpolator that can be exact or smoothed depending on the measurement error model. It is very flexible and allows you to investigate graphs of spatial autocorrelation. Kriging uses statistical models that allow a variety of map outputs including predictions, prediction standard errors, probability, etc. The flexibility of kriging can require a lot of decision-making. Kriging assumes the data come from a stationary stochastic process, and some methods assume normally-distributed data.
Re
mo
te S
en
sin
g a
nd
GIS
in W
ate
r M
an
ag
em
en
t @
Dr.
A.K
.M.
Sa
iful I
sla
mR
em
ote
Se
nsi
ng
an
d G
IS in
Wa
ter
Ma
na
ge
me
nt
@ D
r. A
.K.M
. S
aifu
l Isl
am Cokriging
• Cokriging is a moderately quick interpolator that can be exact or smoothed depending on the measurement error model. Cokriging uses multiple datasets and is very flexible, allowing you to investigate graphs of cross-correlation and autocorrelation. Cokriging uses statistical models that allow a variety of map outputs including predictions, prediction standard errors, probability, etc. The flexibility of cokriging requires the most decision-making. Cokriging assumes the data come from a stationary stochastic process, and some methods assume normally-distributed data.