ArcMap 라이센스없이 Python을 사용하여 Shapefile의 내용을 볼 수 있습니까?

ArcMap 라이센스없이 Python을 사용하여 shapefile의 내용을 볼 수 있는지 궁금합니다. 상황은 ESRI 소프트웨어뿐만 아니라 다양한 응용 프로그램에서 셰이프 파일을 만들 수 있다는 것입니다. 공간 참조, 기능 유형, 속성 이름 및 정의 및 셰이프 파일의 필드 내용을 확인하고 허용 가능한 값 집합과 비교하는 Python 스크립트를 만들고 싶습니다. 조직에 ESRI 라이센스가없는 경우에도이 스크립트가 작동하기를 원합니다. 이와 같은 작업을 수행하려면 ArcPy를 사용해야합니까, 아니면 ArcPy를 사용하지 않고 모양 파일을 파헤칠 수 있습니까?

벡터 및 래스터 데이터를 모두 사용하려면 Python GDAL / OGR API에 익숙해지는 것이 좋습니다 . GDAL / OGR 사용을 시작하는 가장 쉬운 방법은 python (x, y) , Anaconda 또는 OSGeo4W 와 같은 python 배포판을 사용하는 것 입니다.

특정 작업에 GDAL 사용에 대한 자세한 내용 :

• 셰이프 파일 필드 및 유형 가져 오기
• 투사하기

또한 시작을 위해 USU의 다음 자습서를 권장합니다.

• 오픈 소스 GIS를 사용하여 파이썬으로 지오 프로세싱

위 예제에서 차용 한 다음 스크립트는 FOSS 도구를 사용하여 다음 작업을 수행합니다.

1. 공간 참조 확인
2. shapefile 필드 및 유형 가져 오기
3. 사용자 정의 필드의 행에 값이 포함되어 있는지 확인

# Import the necessary modules
from  osgeo import ogr, osr

driver = ogr.GetDriverByName('ESRI Shapefile')
shp = driver.Open(r'C:\your\shapefile.shp')

# Get Projection from layer
layer = shp.GetLayer()
spatialRef = layer.GetSpatialRef()
print spatialRef

# Get Shapefile Fields and Types
layerDefinition = layer.GetLayerDefn()

print "Name  -  Type  Width  Precision"
for i in range(layerDefinition.GetFieldCount()):
    fieldName =  layerDefinition.GetFieldDefn(i).GetName()
    fieldTypeCode = layerDefinition.GetFieldDefn(i).GetType()
    fieldType = layerDefinition.GetFieldDefn(i).GetFieldTypeName(fieldTypeCode)
    fieldWidth = layerDefinition.GetFieldDefn(i).GetWidth()
    GetPrecision = layerDefinition.GetFieldDefn(i).GetPrecision()
    print fieldName + " - " + fieldType+ " " + str(fieldWidth) + " " + str(GetPrecision)

# Check if rows in attribute table meet some condition
inFeature = layer.GetNextFeature()
while inFeature:

    # get the cover attribute for the input feature
    cover = inFeature.GetField('cover')

    # check to see if cover == grass
    if cover == 'trees':
        print "Do some action..."

    # destroy the input feature and get a new one
    inFeature = None
    inFeature = inLayer.GetNextFeature()

Python에서 ArcPy보다 오래된 shapefile을 읽는 많은 모듈이 있습니다 .PyPi (Python Package Index) : shapefiles를보십시오 . GIS SE에는 많은 예제가 있습니다 (예 : [Python] Fiona 검색 ).

모두 지오메트리, 필드 및 투영을 읽을 수 있습니다.

• 나이가 를 OSGeo (GDAL / OGR) , 상기보기 파이썬 GDAL / OGR 요리 책 , 예를 들어
• 또 다른 솔루션은 GDAL / OGR을 기반으로하지만 사용하기 쉬운 (파이썬 사전 : GeoJSON 형식) Fiona 입니다.
• pyshp (shapefile) 는 순수한 파이썬 솔루션입니다
• GeoPanda 는 데이터 분석 도구 를 위해 Fiona 를 사용 하여 모양 파일과 Pandas 를 읽고 씁니다 . 그것을 사용하려면 팬더 를 알아야 합니다.

그러나 PySAL과 같은 다른 모듈 : Python Spatial Analysis Library , Cartopy ( pyshp 사용 ) 또는 Matplotlib Basemap 도 모양 파일을 읽을 수 있습니다.

가장 사용하기 쉬운이다 피오나 ,하지만 당신은 단지 ArcPy, 사용 알고있는 경우 pyshp을 하기 때문에, 를 OSGeo 와 피오나는 해야 GDAL C / C ++ 라이브러리를 설치, GeoPandas이 필요 팬더의 모듈 및 PySAL가 너무 큰 (많은, 많은 다른 치료)

모양 파일의 내용 만 읽으려면 복잡한 것을 필요로하지 않고 ArcPy ( ArcPy : AsShape )로 구현 된 Geo 인터페이스 프로토콜 (GeoJSON)을 사용하면 됩니다.

피오나 (파이썬 사전) :

import fiona
with fiona.open('a_shape.shp') as shp:
     # schema of the shapefile
     print shp.schema
     {'geometry': 'Point', 'properties': OrderedDict([(u'DIP', 'int:2'), (u'DIP_DIR', 'int:3'), (u'TYPE', 'str:10')])}
     # projection
     print shp.crs
     {u'lon_0': 4.367486666666666, u'ellps': u'intl', u'y_0': 5400088.438, u'no_defs': True, u'proj': u'lcc', u'x_0': 150000.013, u'units': u'm', u'lat_2': 49.8333339, u'lat_1': 51.16666723333333, u'lat_0': 90}
     for feature in shp:
        print feature              
{'geometry': {'type': 'Point', 'coordinates': (272070.600041, 155389.38792)}, 'type': 'Feature', 'id': '0', 'properties': OrderedDict([(u'DIP', 30), (u'DIP_DIR', 130), (u'TYPE', u'incl')])}
{'geometry': {'type': 'Point', 'coordinates': (271066.032148, 154475.631377)}, 'type': 'Feature', 'id': '1', 'properties': OrderedDict([(u'DIP', 55), (u'DIP_DIR', 145), (u'TYPE', u'incl')])}
{'geometry': {'type': 'Point', 'coordinates': (273481.498868, 153923.492988)}, 'type': 'Feature', 'id': '2', 'properties': OrderedDict([(u'DIP', 40), (u'DIP_DIR', 155), (u'TYPE', u'incl')])}

pyshp 사용 (Python 사전으로)

import shapefile
reader= shapefile.Reader("a_shape.shp")
# schema of the shapefile
print dict((d[0],d[1:]) for d in reader.fields[1:])
{'DIP_DIR': ['N', 3, 0], 'DIP': ['N', 2, 0], 'TYPE': ['C', 10, 0]}
fields = [field[0] for field in reader.fields[1:]]
for feature in reader.shapeRecords():
    geom = feature.shape.__geo_interface__
    atr = dict(zip(fields, feature.record))
    print geom, atr
{'type': 'Point', 'coordinates': (272070.600041, 155389.38792)} {'DIP_DIR': 130, 'DIP': 30, 'TYPE': 'incl'}
{'type': 'Point', 'coordinates': (271066.032148, 154475.631377)} {'DIP_DIR': 145, 'DIP': 55, 'TYPE': 'incl'}
{'type': 'Point', 'coordinates': (273481.498868, 153923.492988)} {'DIP_DIR': 155, 'DIP': 40, 'TYPE': 'incl'}

osgeo / ogr 사용 (Python 사전)

from osgeo import ogr
reader = ogr.Open("a_shape.shp")
layer = reader.GetLayer(0)
for i in range(layer.GetFeatureCount()):
    feature = layer.GetFeature(i)
    print feature.ExportToJson()
{"geometry": {"type": "Point", "coordinates": [272070.60004, 155389.38792]}, "type": "Feature", "properties": {"DIP_DIR": 130, "DIP": 30, "TYPE": "incl"}, "id": 0}
{"geometry": {"type": "Point", "coordinates": [271066.032148, 154475.631377]}, "type": "Feature", "properties": {"DIP_DIR": 145, "DIP": 55, "TYPE": "incl"}, "id": 1}
{"geometry": {"type": "Point", "coordinates": [273481.49887, 153923.492988]}, "type": "Feature", "properties": {"DIP_DIR": 155, "DIP": 40, "TYPE": "incl"}, "id": 2}

GeoPandas 사용 (Pandas 데이터 프레임)

import geopandas as gp
shp = gp.GeoDataFrame.from_file('a_shape.shp')
print shp
        DIP_DIR    DIP  TYPE                       geometry
0         130       30  incl          POINT (272070.600041 155389.38792)
1         145       55  incl          POINT (271066.032148 154475.631377)
2         155       40  incl          POINT (273481.498868 153923.492988)

* 지오 판다에 대한 참고 사항 이전 버전의 Fiona 및 GDAL을 함께 사용해야합니다. 그렇지 않으면 설치되지 않습니다. GDAL : 1.11.2 피오나 : 1.6.0 Geopandas : 0.1.0.dev-

웹과 책에 대한 많은 자습서가 있습니다 ( Python Geospatial Development , Python으로 지리 공간 분석 학습 및 Python 으로 지오 프로세싱 )

좀 더 일반적으로 ArcPy없이 Python을 사용하려면 Python을 사용 하여 shapefile의 간단한 주제 매핑을 살펴보십시오 .

ArcPy 외에 이러한 기능을 제공하는 지리 공간 파이썬 라이브러리가 있습니다. 다음은 두 가지 예입니다.

파이썬 셰이프 파일 라이브러리 (pyshp)

지오 판다

다른 라이브러리에 관심이 있다면 필수 Python Geospatial 라이브러리 에 대한 이 게시물을 참조하십시오 .