tuGIS 2010 Keynote - The More GIS Changes, the More it Remains the Same

Photo by clarksworth - http://ﬂic.kr/p/61giBZ The More GIS Changes, the
More it Remains the Same James Fee WeoGeo.com

Photo by JKönig - http://ﬂic.kr/p/s9CUo Retro

Early Map Makers

Classic Cartography

Not Much Happens for 300 Years

Photo by sylvar - http://ﬂic.kr/p/o8EML

Photo by sylvar - http://ﬂic.kr/p/o8EML Photo by drmacro - http://ﬂic.kr/p/83i2aa

Computer Science

Esri in the 70’s

ArcInfo

ARC/INFO Required Prime Computer

ArcGIS Desktop

Why The History Lesson?

Photo by dullhunk - http://ﬂic.kr/p/zDLZm

A Simpler Time

Problems with Wizards • Wizards cause us to conform •
Options are either preset or confusing • Lack of documenting the process

Familiarity is a Tough Nut to Crack Most GIS Professionals
only know ArcToolbox

With the Command Line • Freedom to work with the
data as you want • Better understanding of the commands • Combined with scripting you have a recipe for future work. AKA documentation!

The “New” Way

The “New” Way • Start up ArcCatalog (wait 2-4 min
while it loads)

while it loads) • Browse for toolbox (clicking through the endless tree)

while it loads) • Browse for toolbox (clicking through the endless tree) • Open tool dialog

while it loads) • Browse for toolbox (clicking through the endless tree) • Open tool dialog • Browse for input dataset

while it loads) • Browse for toolbox (clicking through the endless tree) • Open tool dialog • Browse for input dataset • Browse for any optional data for analysis

while it loads) • Browse for toolbox (clicking through the endless tree) • Open tool dialog • Browse for input dataset • Browse for any optional data for analysis • Browse for an output location

while it loads) • Browse for toolbox (clicking through the endless tree) • Open tool dialog • Browse for input dataset • Browse for any optional data for analysis • Browse for an output location • Fill out options

while it loads) • Browse for toolbox (clicking through the endless tree) • Open tool dialog • Browse for input dataset • Browse for any optional data for analysis • Browse for an output location • Fill out options • Click OK and pray it doesn’t crash

Photo by TheNixer - http://ﬂic.kr/p/yZvyU Photo by TheNixer - http://ﬂic.kr/p/yZvyU

The “Classic” Way

The “Classic” Way • Start up Command Prompt (about a
second)

second) • Type simple command (BUILD, CLIP, etc)

second) • Type simple command (BUILD, CLIP, etc) • Hit Enter

second) • Type simple command (BUILD, CLIP, etc) • Hit Enter • Tell everyone how hard it was

Python!

Scripting With GIS

Scripting With GIS • AML

Scripting With GIS • AML • SML

Scripting With GIS • AML • SML • Avenue

Scripting With GIS • AML • SML • Avenue •
Visual Basic for Applications *

Visual Basic for Applications * • VBScript

Visual Basic for Applications * • VBScript • JavaScript

Visual Basic for Applications * • VBScript • JavaScript • Perl

Visual Basic for Applications * • VBScript • JavaScript • Perl • Tcl

Visual Basic for Applications * • VBScript • JavaScript • Perl • Tcl • Python

Where is Python Used? • GIS Applications (ArcGIS, QGIS, Safe
FME, Autodesk) • Web Frameworks (WeoGeo.com, Google) • Anywhere (Corel, Civilization, Gmail) • Become the de-facto scripting language

ArcPy, the New Command Line • Much like the Arc
prompt and AML, ArcPy and Python give us a powerful analysis and scripting method. • ArcPy builds on the commands we learned 15 years ago. • ArcPy is integrated in ArcMap • Python can be extended to other GIS platforms such as QGIS.

AML Scripting in ARC/INFO clip soils studbndy stdysoil

Python Scripting in ArcGIS arcpy.Clip_analysis(soils.shp, studbndy.shp, stdysoil)

Python Scripting With ArcPy import arcpy from arcpy import
env env.workspace = "c:/workspace" # variables  in_features = "soils.shp" clip_features = "study_boundary.shp" out_feature_class = "c:/workspace/output/study_area_soils.shp" xy_tolerance = "" # Execute Clip arcpy.Clip_analysis(in_features, clip_features,  out_feature_class, xy_tolerance)

http://www.ﬂickr.com/photos/wwarby/2989238425/ Python

Export to KML import arcpy arcpy.CheckOutExtension(“3D”) env.workspace =
"c:/workspace" # variables  in_feature = "c:/data/TIGER2009/04/ARIZONA/tl_2009_04_county.lyr" out_feature = "c:/temp/output.kmz" # Execute KML Export arcpy.LayerToKML_conversion(in_feature, out_feature,1)

Export to KML arcpy.LayerToKML_conversion(input.shp, output.kml, scale)

Python Libraries

Python Scripting With WeoGeo import WeoGeoAPI #Browse my WeoGeo
Library s = WeoGeoAPI.weoSession('http://jamesfee.weogeo.com', 'user', 'pass') s.connectToMarket() #Find Data Near Tempe, AZ data = s.getDatasets('JSON', '&data_type=VECTOR&per_page=2&page=1&north=33.5&south=33.2&west=-‐111.99&eas t=-‐111.91') #print the JSON response print data

Python Scripting With WeoGeo { 'per_page':
2, 'total_entries': 34, 'current_page': 1, 'total_pages': 17, 'items': [ { 'rating': 0.0, 'projection': 'geo', 'provider_margin': 1995.0, 'uncompressed_misc_files_size': 7683891, 'spatial_resolution': 0, 'children_count': 0, 'datum': 'WGS84', 'library': { 'name': 'Pitney Bowes -‐ Business Insight', 'id': 112 }, 'kml_file_size': 0, 'hosted': True, 'market': 'Complete', 'center_lat': 38.8051135, 'layers': [ 'all' ], 'east': -‐74.986282, 'votes': 0, 'content_license': { 'url': 'http: //licenses.weogeo.com/licenses/8/original.PDF?1273263090', 'name': 'PBBI Software and Data End User License v. April 2008' }, 'data_type': 'VECTOR', 'royalty_model': 'CREDITED', 'west': -‐79.487651, 'scales': '6;7;8;9;10;11;12', 'provider_discount_expires_at': None, 'boundaries': { 'geo': { 'proj4': '+proj=latlong +datum=wgs84', 'north': '39.723622', 'west': '-‐79.487651', 'datum': 'WGS84', 'projection_datum': 'geo-‐wgs84', 'east': '-‐74.986282', 'south': '37.886605' }, 'tiles': { 'number_of_lines': '316', 'number_of_samples': '316', 'datum': 'WGS84', 'line_pixel_size': '-‐830.510836842', 'sample_pixel_size': '1585.72818023', 'proj4': '+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs', 'projection_datum': 'spherical_mercator', 'west': '-‐8848524.83367', 'north': '4825860.68838', 'east': '-‐8347434.72872', 'south': '4563419.26394' }, 'native': { 'number_of_lines': '316', 'number_of_samples': '316', 'datum': 'WGS84', 'line_pixel_size': '-‐830.510836842', 'sample_pixel_size': '1585.72818023', 'proj4': '+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs', 'projection_datum': 'spherical_mercator', 'west': '-‐8848524.83367', 'north': '4825860.68838', 'east': '-‐8347434.72872', 'south': '4563419.26394' }, 'data': { 'proj4': '', 'datum': 'WGS84', 'projection_datum': 'geo-‐WGS84' }, 'baseimage': { 'number_of_lines': 0, 'number_of_samples': 0, 'west': '-‐8848524.833673440000000', 'line_pixel_size': 0, 'sample_pixel_size': 0, 'proj4': 'spherical_mercator', 'projection_datum': 'spherical_mercator', 'north': '4945185.028635530000000', 'east': '-‐8347434.728720820000000', 'south': '4444094.923682910000000' } }, 'tile_layer_type': 'xyz', 'provider_discount_rate': 100, 'provider_discount_expire_option': True, 'x_conv': 1, 'parents_count': 0, 'status': 'Approved', 'north': 39.723622, 'description': 'StreetPro USA \nVersion 2009.12 \nStreetPro offers users a premier street-‐level data product featuring accuracy and street display quality unparalleled in the industry. It reflects real world geographic conditions

Rainbows and Unicorns

Rails, Java, Python

Escape to the Cloud!

Photo by Ezra.Wolfe - http://ﬂic.kr/p/7UUnB9

Hosted GIS

It’s teh awesome! http://www.ﬂickr.com/photos/soundfromwayout/143822346

Analysis

Analysis Storage and ETL

Consumer Analysis

5IBOLT 2VFTUJPOT$PNNFOUT KGFF!XFPHFPDPN !DBHFZKBNFT TQBUJBMMZBEKVTUFEDPN XFPHFPDPN

tuGIS 2010 Keynote - The More GIS Changes, the ...

tuGIS 2010 Keynote - The More GIS Changes, the More it Remains the Same

More Decks by James Fee

Other Decks in Technology

Featured

Transcript