ClueTrust GIS Products


Macintosh GIS

Try Cartographica, the Macintosh GIS today. Explore, analyze, and present geospatial data and maps on your Macintosh.

Get CartoMobile, our field data entry and viewing tool for iOS.

Updates
Blog Index
The journal that this archive was targeting has been deleted. Please update your configuration.
Navigation
Monday
Feb042013

Mapping Shale Gas and Shale Fracturing (Fracking) Sites in the United States

Fracking has become a hot button topic among environmentalists and politicians over the past 10 years. Shale Gas is a type of natural gas found deep within shale rock formations. In 2000 shale gas provided only 1% of U.S. natural gas production; by 2010 it was over 20% and the U.S. government's Energy Information Administration predicts that by 2035 46% of the United States' natural gas supply will come from shale gas (Wikipedia, 2013). There is significant debate about whether fracking should be increased to these levels due to the environmental risks associated with the fracking process. Opponents of fracking argue that the extraction and use of shale gas can affect the environment through the leaking of extraction chemicals and waste into water supplies, the leaking of greenhouse gasses during extraction, and the pollution caused by the improper processing of natural gas. A challenge to preventing pollution is that shale gas extractions varies widely in this regard, even between different wells in the same project; the processes that reduce pollution sufficiently in one extraction may not be enough in another (Wikipedia, 2013). However, proponents of fracking argue that shale gas represents a significant improvement when compared to other fossil fuels, and may even help stem growth in greenhouse gases. Also, shale gas and the associated industry can greatly improve the United States' ability to become energy efficient and self sustaining. No matter you perspective on this issue it may be of interest to you to find out where fracking and the related industries are taking hold in the U.S. In this post we use data from U.S. Energy Information Administration. The image below shows the basic idea behind fracking. 

The first map below shows the shale basins within the U.S. These are areas where there is enough shale gas fracking to occur. To download this data click on the EIA link above and then scroll down about half way down the page until you see the "Geospatial Data in Shapefile (.shp) Format". Click on the the Shapefiles for Basin Boundaries under the Data for Shale Plays Map heading. To import the data choose File > Import Vector Data. To add the live map as a basemap choose File > Add Live Map

 

 The second map below show the areas known as Shale Plays, which are areas that have shale gas currently being harvested using fracking. The areas being harvested are dark red.

 

 

Recently, fracking made the news when NASA released satellite imagery highlighting the growth in fracking in undeveloped areas around the country. Specifically, NASA highlighted the growth in nighttime lights in the North Dakota region, which is a relatively undeveloped part of the county. Despite the low levels of development there are significant nighttime lights visible in areas where shale gas fracking operations are underway. Download a Geo.tiff image created by the Visible Infrared Imaging Radiometer Suite (VIRUS). Once you download the image you can import it into Cartographica by choosing File > Import Raster Data. Below is an image showing the geo.tiff with the Shale Plays layer. 

And a closer look at the increased nighttime lights due to fracking operations in North Dakota. 

 

Wednesday
Jan302013

Analysis with Live Maps: Mapping the Location of Important Naval Ships

Live Maps can be used for many purposes. On this blog we have highlighted the use of Live Maps for georeferencing images, identifying geological features, and providing context to local area studies. Another useful purpose of Live Maps is locating objects that are tied to specific places. What does that mean? It means we can look at places where we expect things to occur and make observations at those locations. In addition to making observations and identifying objects we can also use Cartographica to analyze what we see. 

If you have followed this blog at all you might have notice that I have a light obsession with all things military. Part of the reason is that military objects (especially naval) make interesting maps. Over time I have determined the locations of several interesting military ships and before Bing Maps gets updated I wanted to provide a look at these locations while also conducting some analysis and highlight a few of the functions of Cartographica. Below are a few descriptions of some of the most recent additions to the United States and Chinese Navies. 

The U.S.S. Gerald R. Ford (CVN 78): 

The U.S.S. Gerald R. Ford is the newest addition of the U.S. Navy Aircraft Carrier and is the first ship in the new Gerald R. Ford Class of supercarrier. The new carrier comes at a cost of $13.5 billion and it includes numerous improvement over past classes of carriers. One of the biggest improvements is the aircraft launch system which moves from steam power to electromagnets. Like other nuclear powered carriers the GRF will have an unlimited service range for a period of 25-30 years and will only need to come to port for supplies and regular maintenance. The ship is being built by Huntington Ingalls Industries in Newport News, VA. Below is an image of the port where the ship is being constructed. Like a game of ISpy, do you see the ship?

To highlight the location of the of the GRF add a new feature by choosing Layer > New Layer and then Edit > Add Feature. Select to add a polygon feature and then trace the outline of the ship. 

A closer look at the U.S.S. Gerald R. Ford under construction.

Chinese Aircraft Carrier: Liaoning

Recently the Chinese military acquired an 67,500 ton Soviet era aircraft carrier and has spent the past several years refurbishing and upgrading the ship to make it battle ready. Last November, China landed its first plane on the surface of the Liaoning. See this CNN video of the plane landing. Below is an image of Dailan, China where the Liaoning has been under construction for several years (its has completed tests in the Yellow Sea). Again, like a game of ISpy, do you see the ship?

We can again create a new feature to show the location of the ship choose Layer > New Layer and then Edit > Add Feature. Select to add a new polygon feature and then draw the outline of the ship. 

A closer look:

Based on the polygons that we have create we can use Cartographica's table tools to help enhance what we know about the ships. First, add a new area column to each of the new polygons so that we can see the size of each ship's deck space. To add the area column choose Tools > Add Area Column. Based on this analysis the Gerald R. Ford has a deck space of 32,356 square meters and the Lioaning has a deck space of 27,494 square meters. Also, we can add polygon coordinate columns which will give us the coordinates for each of the ships. To add coordinate columns choose Tools > Add Centroid Coordinate Columns. The final map below shows the general location of the ships. 

Monday
Jan282013

Creating Elevation Contour Maps from Digital Elevation Models

A new feature in Cartographica version 1.4 is the ability to create contour maps. In GIS, a contour line joins areas of equal elevation above a given level, such as mean sea level. contour map is a map illustrated with contour lines, which show valleys and hills, and the steepness of slopes. In this sense a contour map can be very useful in many contexts. One method for creating contour lines is by using Digital Elevation Models (DEM). A DEM is a digital model or 3D representation of a terrain's surface — commonly for a planet (including Earth), moon, Mars, or asteroid. Typically DEMs are created based on data that are retrieved through remote sensing technology. Remote sensing technology typically refers to and includes specialized sensors that are attached to various satellites or aerial vehicles that makes observations on the surfaces of their target object. 

In this post we emphasize the use of the new contour mapping functions in version 1.4 by using data from the Kentucky Division of Geographic Information. The data on the DGI website are free for download and include DEMs for the entire state of Kentucky.  To illustrate the contour mapping we will use data from the Middlesboro, KY area.  An interesting geological feature about Middlesboro is that experts believe that its location between Pine Mountain and the Cumberland Mountains is actually an ancient crater from an asteroid impact. This fact makes Middelsboro among the few cities in the world that is seated within an impact crater. 

In order to view the entire area for Middlesboro we actually need to download 4 separate DEMs. We need the cells U47, U48, V47, and v48 from the DGI link listed above. To download the DEMs individually control-click within each cell and then select Download Linked File as. Save the Files to your Desktop.   

To import the DEMs choose File > Import Raster Data. The following image shows what your map should look like. 

Notice the blue colored circular area near the center of the map. That is the location of Middlesboro and the fairly clear outline of an impact crater. Next, we are going to create contour lines for each of the DEMs that will show the differences in elevations throughout the map. To create contour lines select the DEM layer in the layer stack and then choose Tools > Create Contour Lines. This step has to be done individually for each of the DEM layers. A window will appear that allows you to choose the increments and the base. Select 50 as the increment and leave the base at 0. See below for the contour window. 

When the contour lines are created they will all be black in color. In order to enhance the visibility of the contours we need to adjust the color schemes. Double-click on the U47 DEM in the Layer Stack. Change the Based on selection to Elevation. Click the + button below the table 6 times and then click the Gear box and select Distribute with Natural Breaks (Jenks). Next, choose Window > Show Uber Browser, click on the Palettes tab, and while holding down the option key click and drag a color palette of your choice into the table within the Layer Styles Window. In the map presented below a blue-red scheme was used where blue = lower elevations and red = higher elevations. Repeat this process for each of the DEMs. See the images below for an example of the Layer Styles window and the map.

Layer Styles Window:

Contour Line Map:

The next image is a closer look at the Impact Basin. Also the following image has the DEM layers turned off and a Live Map added. To add a live map choose File > Add Live Map.

Thursday
Jan032013

Street Segments as Units of Analysis: Spatially Joining Points and Lines

Recently, authors David Weisburd, Elizabeth Groff, and Sue-Ming Yang published a book titled The Criminology of Place: Street Segments and Our Understanding of the Crime Problem. The crux of the book is that in Criminology and other disciples have a long history of using geographic units of analysis that are area based. That is, we have become accustomed to using polygons as units of analysis. Part of the reason for the use of area based units of analysis is that demographic and social data are available through sources like the United States Census Bureau. In the book, the authors make a compelling argument that area based studies miss a lot of spatial variation that occurs at lower levels of analysis. In an attempt to combat the loss of information through aggregation the authors show that using street segments as units of analysis may be a viable alternative to area based units of analysis. 

In the spirit of using street segments as units of analysis we thought highlighting how to use Cartographica to aggregate data points to street segments would be a fun and interesting exercise. Recent upgrades to Cartographica 1.4 include the ability to spatially join data based on their spatial association. This includes joining points to lines. 

The data being used in this example are from Washington D.C.'s GIS repository DCGIS. The data are 2011 crime incident locations and a street file containing all of the street segments in the city. To import the data into Cartographica choose File > Import Vector Data. Below is an image of the data described above. At first glance it appears as though the points are "on top of" the line segments. Observe. 

However, upon closer inspection you can see that the points are not directly on top of the lines. As a result, we will need to use a distance based spatial join to join the crime points to the line segments. 

Because we want to join the crime points to the street segments we first need to select the Streets layer in the layer stack. Once selected, choose Tools > Spatial Join. The following image shows the set up for the Spatial Join. Uncheck the Discard Unmatched Features so that the new join layer will have all of the street segments including those that have zero crimes. Note that the "Within Distance" operation is selected and a distance of 50 meters is the designated distance. Also note in the table that the 2011 crime data fields are set to ignore. The reason for this is that at this point we do not need all of the excess crime information. All we want is the counts of crimes on the various street segments. Finally, note that the last field in the table "Join Count" is selected to be copied. This will be the new field on the street segments layer that will contain the crime counts. 

 

After the points are counted and aggregated to the street segments you can reclassify the counts into groups and then vary the street segment colors by the groups. Double-click on the streets join layer in the Layer Stack to brig up the Layer Styles Window. Add three categories to the table by click on the + button. Click on the Gear Box and select Distribute with Natural Breaks (Jenks). Change the colors of each of the categories by clicking on the fill box and choosing a series of colors. See the map below for the final result. 

Friday
Dec282012

Best Open Geodata Release of 2012: Philadelphia Open Data 

With the end of the year approaching everyone is putting out their "Best of" list for 2012. Emily Badger of The Atlantic Cities recently released a list of the "Best Open Data Releases of 2012". These open data sources are available through various agency and government websites and allow users access to data sources that include geospatial data. 

We decided to check out a few of the Open Data sources and wanted to create a few maps using Cartographica. The Philadelphia Open Data website was listed as the best Open Data website of 2012. We were able to download a .csv file for all of the crime incidents in Philly since 2006 from the website. Download the Philly Crime data.  

To start, add a Live Map to use as a baseamp by choosing File > Add Live Map

To geocode the data using Cartographica choose File > Import Table Data. The Import File window will appear. Select the coordinates tab in the top right. Change the Map to selection for the Point X and Point Y fields to X (or longitude) and Y (or latitude) and then click Import. See below for an example of the set up. 

*When we geocoded the Philly Crime data 164 incidents were geocoded to the Cincinnati area. Upon further inspection the points were geocoded to the appropriate location based on the data entries. The X and Y coordinates for all of the 164 incorrect points were, -84.69369121 longitude and 39.1207947 latitude, which places the points at the corner of Salyer Ave and Goodrich Lane just West of Cincinnati. These points were considered erroneous and were discarded. To delete the points select them using the identify tool and then choose Edit > Delete. 

The map below shows the geocoded crime incident points in Philadelphia since 2006. The Philly_Crime layer contains 592,064 crime points. 

You can filter the data by using the Filter bar. Change the Filter bar selection to Text_General_Code and then type in Homicide. The data will be filtered to only show homicide incident locations. Since 2006, there were 3051 homicides in Philadelphia. Use the identify tool to select all of the points and then hold down the option button and choose Tools > Create Kernel Density Map for Selection. The result of the KDM analysis is shown below. The red areas indicate places where homicide was highly concentrated.

 

A closer look at central Philadelphia.

Another way to filter the data is by date. The Philly_Crime layer contains a field called Dispatch_2 that contains the date of the offense. With the Philly Crime layer selected in the layer stack type in 12/31/2010. A small box will appear at the top of the Map Window. You can filter based on multiple criteria to limit the dataset to any year that you want. To select only the 2011 data use the settings shown below. Notice that we are selecting crime incidents that fall between 12/31/2010 and 1/1/2012. 

Use the identify tool to select all of the 2011 crime points and then choose Layer > Create Layer from Selection. The new layer should have 82197 crime points. You can repeat these steps so that you can create a new layer for each of the years in the data set. Create another Kernel Density map by Choosing Tools > Create Kernel Density Map. The result of the KDM analysis for the year 2011 data are shown below.