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.

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

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. 


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.


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. 


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. 


Mapping the Syrian Conflict

The violence associated with the Syrian revolution continues to escalate as recent reports have indicated that Bashar Al-Assad's forces have begun to use SCUD missiles in attempts to drive back rebel forces in the area near the Syrian Capital. SCUD missiles are infamous for their use by former Iraqi dictator Saddam Hussein during the Persian Gulf War. Equally as infamous are the U.S. produced Patriot Missile systems that are designed to eliminate surface-to-surface missile attacks. Recent reports have indicated that American, German, and Dutch patriot missile battalions have been sent to Turkey to eliminate the risk to the Turkish people from possible SCUD missile attacks from across the Syrian border. Check out the article on the Huffington Post about the recent visit from U.S. Defense Secretary Leon Panetta's visit to Turkey regarding the Syrian conflict and the U.S. Patriot Missile response. See below for a diagram from Raytheon for how the Patriot Missile system works. 

In order to add more context to situation we have created a few maps to highlight the areas where the Patriot Missile systems are being deployed. The Patriot Missile systems are being deployed to Kahramanmaras, Turkey, which is about 70 miles North of the Syrian border. We can confirm the distance between Kahramanmaras and Syria by using Cartographica's Live Maps and measurement tools. See the image below for an example. To add a Live Map, Choose Tools > Add Live Map, and to use the measurement tools click on the measure tool button and then drag a line to determine the distance. 

We can also look at the radar and missile ranges of the Patriot Missile systems by using buffers. The radar used by the Patriot Missile system has a detection range of about 180 miles in all directions around the system. The map below shows how far the systems deployed in Kahramanmaras will be able to detect incoming SCUD missiles. To create the map you first need to add a new layer by choosing Layer > Add Layer. Add a new feature to the layer by choosing Edit > Add Feature and then select to add a new Point feature. Create a buffer for the point by Choosing Tools > Create Buffer for Layer's Features. Make the buffer for 180 miles. 

You can add to the map above by also creating a buffer for the effective range of the Patriot Missiles. The range is about 70 miles. Using the same steps as above create a second buffer for 70 miles. The map below shows the results.