<< VBユーザーのためのGDAL/OGR | top | ArcGISの横暴さがイヤ >>

Ubuntu(8.04 Hardy Haron)でPython環境構築

Ubuntu上での環境構築は、Windows、Mac OS X に比べてとても簡単である。
Synapticパッケージマネージャですべて導入できる。
依存するパッケージの導入もしてくれるので、次のものを指定するのが一番楽かも。

・spe
・python-gdal
・python-mapscript
・python-qt4
・qt4-designer

・今回は、GeoTiff を読み込ませようと思い、以下からGeoTiffのデータをダウンロードしました。このGeoTiffは、ヘッダの部分に座標系が埋め込まれているはずです。これを表示させてみましょう。

 ■MapTools Download: /dl/geotiff/samples/other
 http://dl.maptools.org/dl/geotiff/samples/other/
 この中の、erdas_spnad83.tif をダウンロードします。 

・SPEを起動して、

import mapscript
import gdal,gdalconst,ogr,osr

aa = mapscript.mapObj("")

print aa.name

gdal.AllRegister
dataset = gdal.Open('/Users/bugna/Projects/Python/erdas_spnad83.tif', gdal.GA_ReadOnly)


print 'Driver: ', dataset.GetDriver().ShortName,'/', dataset.GetDriver().LongName
print dataset.GetProjection()


などとして実行すると、

MS
Driver: GTiff / GeoTIFF
PROJCS["NAD83 / Georgia East",GEOGCS["NAD83",DATUM["North_American_Datum_1983",SPHEROID["GRS 1980",6378137,298.2572221010002,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","6269"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4269"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",30],PARAMETER["central_meridian",-82.16666666666667],PARAMETER["scale_factor",0.9999],PARAMETER["false_easting",656166.6666666666],PARAMETER["false_northing",0],UNIT["US survey foot",0.3048006096012192,AUTHORITY["EPSG","9003"]],AUTHORITY["EPSG","26966"]]

のような結果を得ることができます。

※GDAL のバージョンが、1.4 なので、Namespace osgeo を『つけ
ない』。import 行に注意。

bugna | Python | 07:49 | comments(0) | trackbacks(0) | - | - |
Comment









Trackback
URL:

11
--
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
--
>>
<<
--