Worksheet 7 - GIS

Question 1 - Qgis

Lets begin with a brief introduction to GIS mapping using an open source application Quantus GIS. This will be a pictorial demonstration of some of the fundemental principals of mapping data in a GIS application.

Qgis

Q1-1. Lets begin by constructing a map of the Royal Botanical Gardens Cranbourne.
adding a GeoTIFF raster layer that represents an arial photograph of the Royal Botanical Gardens Cranbourne.
  1. Add a GeoTIFF raster layer that represents an arial photograph of the Royal Botanical Gardens Cranbourne (Cranbourne_HighLevel_0,7m.tif)
  2. Add a vector layer that represents the vegetation classes of the Royal Botanical Gardens Cranbourne (vegetation.shp)
  3. Adjust the properties of this layer. In particular, alter the legend type to Unique value and modify the classification Field to EVC_DESC. This will therefore color each of the vegetation types differently. You might want to alter the fill colors and styles of the individual vegetation classes. You might also want to alter the transparency of the entire layer.
  4. Add a vector layer that represents the walking and management tracks of the Royal Botanical Gardens Cranbourne (sandtracks.shp)
  5. Add a vector layer that represents the roads of the Royal Botanical Gardens Cranbourne (mainroads.shp)

Q1-2. As part of an honours project investigating habitat use and activity patterns of southern brown bandicoots, a student last year tracked a number of male and female bandicoots
  1. Add a vector layer that represents the radio tracking fixes of the six female and seven male bandicoots (coots.shp). You might want to alter the colors and or symbols so that different individuals can be highlighted.
  2. Add a new layer

Question 2 - Multidimensional scaling

The following example is designed to help you appreciate the link between distance measures and ordination space (MDS). The data set consists of distances (km) between major Australia cities (as the crow flies), and is in the form of a triangular matrix.

Format of austcities.csv data file
 CanberraSydneyMelbourne..
Canberra0NANA..
Sydney2460NA..
Melbourne4677130..
Adelaide9581160653..
Perth309032902720..
..........
Australian cities

Open the austcities data file. Note the format of the file, it is a triangular distance matrix.

While the file is a distance matrix, at this stage R is unaware of it, we must manually make it aware (a round about way of saying that we must type a command to force R to treat the data set as a distance matrix. We do this by typing the following command in the top window (Script Window) of R Commander (actually we could alternatively type it on the command prompt of RGui).
aust.cities.dis <- as.dist(aust.cities)
where aust.cities is the name you gave the Australian cities data set when you imported it, as.dist is the name of the function that forces a data set to be considered a distance matrix, and aust.cities.dis is a new name that I am defining to store the distance matrix.
If you chose to enter the statement in the R Commander command log window, once you have completed typing the above statement, press the Submit button in the top right hand corner of R Commander. This will case the statement to be evaluated and the task to be completed.

We are now ready to perform the MDS for the purpose of examining the ordination plot.

Q2-1. Perform an MDS with 2 dimensions on the city distances matrix.
  1. What was the final stress value?
  2. What does this stress value suggest about the success of the MDS?
  3. The Sheppard diagram (plot) represents the relationship between the original distances (y-axis) and the new MDS distances (x-axis). Does this and the stress value indicate that the patterns present in the original distance matrix (crow flies distances between cities) are adequately reproduced from the 2 new dimensions?
  4. The final ordination plot summarizes the relationship between the cities. Does this ordination plot approximate the true geographical arrangement of the cities?
  5. In this case, what might the two new MDS dimensions (variables) represent? (hint think of the ordination plot as a map)

Question 3 - Multidimensional scaling

A fish ecologist investigated differences in fish community composition associated with flow regulation in rivers of four major New South Wales catchements (Murray, Darling, North Coast and South Coast). A subset of the community data is presented in the file riverfish.csv.

Format of riverfish.csv data file
 H.compG.macM.flu..
Darling Unreg0062..
Darling Reg0042..
Murray Unreg000..
Murray Reg002..
North Coast Reg165700..
..........
trout

Open the riverfish data file.

Q3-1. Examine the whole data set. Before proceeding with any formal analysis, lets explore the patterns in the data visually.
  1. Based solely on Hypseleotris compressa which sites are most similar?
  2. What do the abundances of Nematalosa.erebi suggest about the sites?

Of course, we could explore the patterns amongst sites according to each separate fish species. However, in the full investigation, the fish ecologist had 51 species of freshwater fish, each of which yield slightly different patterns. Furthermore, the fish ecologist was not interested in the patterns of any one species of fish. What he was interested in was whether (and how) the fish communities differed between catchements and flow regimes. Consequently, an MDS was used to explore the patterns amongst sites based on all the fish species.

Q3-2. Calculate the Bray-Curtis(Czekanowski) dissimilarity coefficients amongst the sites using all six fish species (note we are only using a subset of the entire data set).
  1. Which sites are most similar?

Q3-3. Perform an MDS with 2 dimensions on the dissimilarity matrix.
  1. What was the final stress value?
  2. What does this stress value suggest about the success of the MDS?
  3. The Sheppard diagram (plot) represents the relationship between the original distances (y-axis) and the new MDS distances (x-axis). How would you describe the shape of this curve, and base on this is metric or non-metric MDS more appropriate?
  4. The final ordination plot summarizes the relationship between the sites. What would you conclude from this?

Question 4 - Multidimensional scaling

Ludwig and Reynolds (1988) described a data set in which a number of sites where characterized by the abundances of five species of cockroach.

Format of Ludwig.csv data file
 C.sppC.cucullatusC.delicatulus..
BCI01428..
LC040..
FORT011..
BOQ010..
MIR041..
CORG100..
cockroach

Open the Ludwig data file.

Q4-1. Examine the whole data set. Are there any patterns

Q4-2. Calculate the Bray-Curtis(Czekanowski) dissimilarity coefficients amongst the sites using all cockroach species.
  1. Which sites are most similar?

Q4-3. Perform an MDS with 2 dimensions on the dissimilarity matrix.
  1. What was the final stress value?
  2. What does this stress value suggest about the success of the MDS?
  3. The Sheppard diagram (plot) represents the relationship between the original distances (y-axis) and the new MDS distances (x-axis). How would you describe the shape of this curve, and base on this is metric or non-metric MDS more appropriate?
  4. The final ordination plot summarizes the relationship between the sites. What would you conclude from this?

Welcome to the end of Worksheet 6