GDAL介绍.ppt

Ex. 3: Generate Near Table tool For each point, continued: Get the distance to the first line and its FID Loop through the rest of the lines and if the distance is less than the “remembered” distance, then remember the new distance and FID Write out the point FID and the remembered line FID and distance to a text file Close the shapefiles and the text file Needs an ArcInfo license in ArcMap (can probably write your own script, though) Ex. 5: Compute NDVI Set the georeferencing and projection for the new image Build pyramids for the new image Needs the Spatial Analyst extension, and I don’t know if you could do it with anything in the geoprocessor without SA – I’m sure you could do it with ArcObjects Ex. 6: NDVI in ArcMap Run the NDVI script in ArcMap Change the script to accept parameters Change the script to send error messages to the geoprocessor Add the script to a toolbox * Form data using Language=Python tag in asp. Def getdata() retrieves variables from form. Getdata() gets variables for use in Python * Project() transforms xy from one projection to another written with Response.Write * PLSS uses pyodbc shows ability to query Microsoft SQL database and return results * Shapefile() uses spatial queries to return data from dbf * Geocoding Service * Create shapefile and write prj file Introduction to Open Source RS/GIS programming with Python Chris Garrard RS/GIS Laboratory Utah State University Why bother with FOSS? Free and Open Source Software (FOSS) Affordable (free, as in free beer) Fits the budget of small companies or individuals Open (as in free speech) Use it for whatever you want Distribute it if you want Modify it if you want Why bother with FOSS? Helpful user community Including the developers themselves Fast bug fixes Not limited to Windows Lightweight Cons Lightweight Not the ESRI geoprocessor Smaller user community than some other packages Free != inferior Apache, MySQL, Linux, Firefox, etc. ESRI likes Python GDAL is used by popular commercial