Using Python, GDAL and NumPy for spatial analysis使用Python空间分析的GDAL和NumPy.pptVIP

Using Python, GDAL and NumPy for spatial analysis and modeling Overview of GDAL and NumPy Where to use these instead of ArcGIS, Imagine, etc. Description of a real-world application Step through two examples Mess around with NumPy functionality How to install the environment using Python 2.4 GDAL (Geospatial Data Abstraction Library) GDAL is a “translator library for raster geospatial data formats” Open source Used in many applications: GRASS, UMN MapServer, Google Earth, ArcGIS 9.2, etc. Can handle many image formats for read and slightly fewer for write: AI Grid, Imagine, GeoTiff, JPEG, PNG, NetCDF, etc. GDAL (Geospatial Data Abstraction Library) Presents an “abstract data model” for processing spatial data Can be used directly from C/C++ and can be “wrapped” for use with Python, Perl, VB, C#, R, Java … Early developers have chosen Python as their scripting language and documentation is relatively good for this. NumPy(Numerical Python) An array/matrix package for Python Well suited for image processing – i.e. one function can operate on the entire array Slicing by dimensions and applying functions to these slices is concise and straightforward Nearly 400 methods defined for use with NumPy arrays (e.g. type conversions, mathematical, logical, etc.) GDAL and NumPy Since GDAL 1.3(?), GDAL has implemented NG (New Generation) Python bindings which includes NumPy Process: Open GDALDataset Get raster band(s) Convert the raster band(s) to a NumPy array using ReadAsArray() Process the raster band(s) using NumPy functionality Convert the NumPy array to GDAL raster bands using WriteAsArray() Write out GDALDataset Why use GDAL/NumPy instead of canned GIS software? Not advisable if what you want to do is easily handled within ArcGIS/Imagine/etc. – there is a lot of programming overhead Well suited for process model applications where the logic at a cell based is too complex Example: Grid algebra : grid1 + grid2 (probably use GIS) Finding NN in multidimensional space (may