[GIS] ArcObjects, 속성값으로 쿼리(Query)하기

GIS의 강력한 기능 중에 하나는 속성 데이터에 기반한 공간적인 처리 및 가시화인데, 이 글은 ArcObjects에서 속성값 조건에 대해 쿼리(Query)하는 방법에 대해 설명하겠습니다. 일단…. 속성값에 대한 쿼리가 목적이므로 먼저 속성값에 접근할 수 있는 ITable 인터페이스에 접근할 필요가 있습니다. 다음은 d:\__data__ 폴더 안의 seoul.dbf 파일로부터 ITable 인터페이스의 구체 클래스 인스턴스를 얻는 코드입니다.

ESRI.ArcGIS.Geodatabase.IWorkspaceFactory pWF;
ESRI.ArcGIS.Geodatabase.IFeatureWorkspace pFW;

pWF = new ESRI.ArcGIS.DataSourcesFile.ShapefileWorkspaceFactory();
pFW = pWF.OpenFromFile("d:/__data__", 0) 
                as ESRI.ArcGIS.Geodatabase.IFeatureWorkspace;

ESRI.ArcGIS.Geodatabase.ITable pTable;
pTable = pFW.OpenTable("seoul");

ITable의 구체 클래스의 인스턴스인 pTable은 속성값으로 검색할 수 있는 매서드인 Search 함수를 가지고 있습니다. 이 Search 함수는 두개의 인수를 가지며 첫번째 인수의 타입은 IQueryFilter로써 SQL 문의 WHERE 절의 조건문을 지정하기 위해 사용됩니다. 그리고 두번째 인수는 검색 결과에 대한 Row에 대한 재활용 여부에 대한 bool 형입니다. 이에 대해서는 이 글의 마지막에 간단히 다시 설명하겠습니다.

속성값으로 쿼리하기 위하여 조건문을 만들기 위해서… 먼저 앞의 코드를 통해 오픈한 seoul.dbf 파일의 필드 구조를 알 필요가 있습니다. 힌트를 준다면 seoul.dbf의 속성 중의 하나는 SGG_NM으로써 서울시를 구성하는 행정구의 이름값이 저장되어 있습니다. 여기서는 행정구 이름에 ‘동’자가 들어가는 모든 Row를 검색하려고 합니다. 이런 조건을 위한 ITable의 Search의 첫번째 인자인 IQueryFilter의 CoClass(구체 클래스)의 정의는 다음과 같습니다.

ESRI.ArcGIS.Geodatabase.IQueryFilter pQF;
pQF = new ESRI.ArcGIS.Geodatabase.QueryFilter();
pQF.WhereClause = "SGG_NM LIKE '%동%'";

IQueryFilter의 WhereClause 프로퍼티가 바로 ‘그것’인데.. 일반적인 SQL의 SELECT 문의 WHERE 절의 바로 정확히 ‘그것’에 해당합니다. SQL문의 SELECT문을 알고 있는 개발자라면 아~!라는 탄식이 나올정도로 ArcObjects는 속성값으로 데이터를 쿼리하는 탁월한 방법을 사용하고 있습니다. 이렇게 만든 IQueryFilter를 이용하여 쿼리하는 코드는 다음과 같습니다.

ESRI.ArcGIS.Geodatabase.IFeatureCursor pFC;
pFC = pTable.Search(pQF, true) as ESRI.ArcGIS.Geodatabase.IFeatureCursor;

ITable의 Search 함수를 사용하며 검색 결과 Row는 IFeatureCursor 타입으로 반환됩니다. 검색 조건과 일치하는 Row를 얻기 위해 IFeatureCursor 타입의 인스턴스 변수인 pFC를 사용하는데.. 아래의 코드가 바로 그것이며, listBox1 이름의 ListBox 컨트롤에 그 결과를 채웁니다.

ESRI.ArcGIS.Geodatabase.IFeature pFeature;
int cntFields = pFC.Fields.FieldCount;

listBox1.Items.Clear();
pFeature = pFC.NextFeature();
while(pFeature != null) {
    for(int iField=0; iField

5번째 코드에서 가장 먼저 IFeatureClass의 NextFeature 매서드를 호출하여 첫번째 검색 결과의 Row에 해당되는 IFeature 클래스 타입의 인스턴스를 반환받고 첫번째 Row에 대한 값을 얻기 위하여 IFeature의 get_Value 매서드를 사용하여 값을 얻습니다. 그 다음 결과 Row를 얻기 위하여 다시 NextFeature를 호출합니다. NextFeature의 반환값이 null이면 더 이상 조건에 맞는 Row가 없다는 의미입니다. 실행 결과는 아래와 같습니다. 강동구, 성동구, 동작구, 동대문구와 같이 '동'자가 들어간 데이터만이 검색된 것을 알 수 있습니다.

사용자 삽입 이미지
끝으로 앞서 언급한... 이 글의 핵심 매서드인 ITable의 Search 매서드의 두번째 인자의 의미에 대해서 설명하면... 만약 true라면 IFeatureClass의 NextFeature의 결과에 대한 인스턴스를 재활용한다는 의미이며 false이면 재활용하지 않고 NextFeature를 호출할때 마다 새로운 인스턴스를 생성한다는 의미입니다. 이를 목적에 대해 생각해보면 true이면 Feature를 읽기 전용으로만 사용해야 하며 false일 경우 IFeature를 통해 속성값을 변경할 수 있습니다.

[GIS] ArcObjects, 속성 테이블(DBF)의 필드 정보 얻기

ArcGIS에서 Feature는 좌표에 대한 도형 데이터와 속성으로 묶어진 단위입니다. ESRI에서 정의한 Shape 파일의 경우 좌표에 대한 도형 데이터는 .SHP 파일에 저장되며 속성은 .DBF 파일에 저장됩니다. 속성은 하나의 DB 테이블 개념과 정확히 일치합니다. DB 테이블의 필드 정보와 각 필드에 대한 값들의 집합의 구성 요소인 로우(Row) 또는 레코드(Record)로 구성됩니다.

필드 정보라 함은 필드의 이름과 필드의 타입 그리고 필드 타입에 대한 데이터 저장 가능 길이 등으로 구성됩니다. 이 포스트는 Shape 파일에 대해서 .DBF에 저장된 필드 정보를 얻어내는 ArcObjects의 API에 대해서 설명합니다.

가장 먼저 해야할 작업은 .DBF 파일을 여는 것에서 시작합니다. 다음 코드는 .DBF 파일을 여는 방법입니다.

ESRI.ArcGIS.Geodatabase.IWorkspaceFactory pWF;
ESRI.ArcGIS.Geodatabase.IFeatureWorkspace pFW;

pWF = new ESRI.ArcGIS.DataSourcesFile.ShapefileWorkspaceFactory();
pFW = pWF.OpenFromFile("d:/__data__", 0) 
              as ESRI.ArcGIS.Geodatabase.IFeatureWorkspace;

ESRI.ArcGIS.Geodatabase.ITable pTable;
pTable = pFW.OpenTable("seoul");

.SHP 파일을 여는 것과 매우 동일하다는 것을 알 수 있습니다. 다른 부분은 8, 9번 라인에서 ITable 타입으로 하여 테이블을 여는 것입니다. 위의 예는 d:/__data__ 폴더 안에 seoul.dbf 파일을 열어 그 테이블을 ITable 타입인 pTable에 저장해 놓는 경우입니다.

이렇게 얻은 pTable을 통해 필드 정보를 얻기 위한 코드는 아래와 같습니다.

ESRI.ArcGIS.Geodatabase.IFields pFields;
pFields = pTable.Fields;

int cntFields = pFields.FieldCount;

테이블의 필드 정보는 IFields 타입에 저장되며 pTable의 Fields 속성을 통해 접근이 가능합니다. 4번째 줄은 필드의 전체 개수를 얻는 코드입니다. 이제 이렇게 얻은 pFields를 통해 각 필드의 이름, 별칭(Alias), 필드길이, 정밀도, 타입을 얻어 ListBox 컨트롤에 채워 보는 코드는 아래와 같습니다.

listBox1.Items.Clear();
for(int iField=0; iField<cntFields; ++iField) 
{
    ESRI.ArcGIS.Geodatabase.IField pField;
    pField = pFields.get_Field(iField);

    string AliasName = pField.AliasName;
    string FieldName = pField.Name;
    int Length = pField.Length;
    int Precision = pField.Precision;
    ESRI.ArcGIS.Geodatabase.esriFieldType FT = pField.Type;
                
    listBox1.Items.Add(FieldName + "(" + AliasName + ")");
    listBox1.Items.Add("  Length: " + Length.ToString());
    listBox1.Items.Add("  Precision: " + Precision.ToString());
    listBox1.Items.Add("  Type: " + FT.ToString());
    listBox1.Items.Add("");
}

참고로 ListBox 컨트롤의 이름은 listBox1입니다. 앞서 구한 전체 필드의 개수를 통해 반복문에서 필드 정보를 IField 타입으로 가져와 필드의 정보를 구할 수 있습니다. 이름은 IField 타입의 Name 속성으로 알 수 있고, 별칭명은 AliasName, 길이는 Length, 정밀도는 Precision 그리고 필드의 타입은 Type으로 알 수 있습니다. 참고로 아래는 실행 결과 화면입니다.

사용자 삽입 이미지

[GIS] ArcObjects, Shape 파일로부터 레이어 추가하기

ArcObjects는 GIS 어플리케이션 개발에 적용할 수 있는 매우 방대하며 매우 강력한 라이브러리입니다. 기본적으로 ArcGIS를 설치하면 개발자가 쉽게 사용할 수 있도록 되어 있습니다. 필자는 2002년 초에 처음 ArcObjects를 접했는데… 그 방대함과 강력함에 매료되었습니다. 제가 아는 GIS 개발 범위에서 필요한 기능은 모두 ArcObjects에서 구할 수 있었다고 해도 과언이 아닐 정도로 말입니다.

간단한 수치지도 도시, 네트워크 및 위상 분석, 그리드 분석, 3D에 대한 틴, 모델 등의 표현과 관련된 공간 분석 그리고 강력한 편집 기능 등… 이루 다 말할 수 없을 정도로 방대한 GIS 개발 도구를 제공합니다. 처음 ArcObjects를 접하고 수년이 흘렀지만 여전이 ArcObjects는 GIS 어플리케이션 개발에서 매우 큰 영역을 차지하고 있으며 향후에도 개발자와 GIS 개발 회사에 미치는 영향력은 더 커지거나 지금처럼 계속 유지될 것으로 기대됩니다.

ArcObjects를 이용한 GIS 시스템 개발에 대한 일을 맡으면서 ArcObjects를 활용해야할 일이 생겼습니다. 과거 ArcObjects를 사용하여 개발했던 프로젝트의 경험에 대한 기억이 가물가물해져 버렸습니다. 해서… 그때의 실수가 되풀이 되지 않도록 ArcObjects에 대한 경험을 글로 하나 하나 남겨 놔야 겠다는 생각이 들어 ArcObjects 관련 글을 포스팅하기로 하였습니다.

그 첫번째는 GIS에서 수치지도 교환에 거의 표준으로 활용되고 있는 ESRI의 SHP 파일을 추가하는 방법입니다. 개발언어는 C#를 사용할 것이며 Form 기반의 어플리케이션의 폼에 ArcObjects에서 제공하는 MapControl을 올려 놓습니다. MapControl의 정확한 클래스 타입은 AxMapControl입니다. 그리고 버턴을 하나 올려 놓습니다. 이 버턴 안의 코드는 지정된 폴더에 존재하는 SHP 파일을 MapControl에 추가하는 코드입니다.

ESRI.ArcGIS.Geodatabase.IWorkspaceFactory pWF;
ESRI.ArcGIS.Geodatabase.IFeatureWorkspace pFW;
ESRI.ArcGIS.Carto.IFeatureLayer pFL;

pWF = new ESRI.ArcGIS.DataSourcesFile.ShapefileWorkspaceFactory();
pFW = pWF.OpenFromFile("d:/__data__", 0) 
              as ESRI.ArcGIS.Geodatabase.IFeatureWorkspace;
pFL = new ESRI.ArcGIS.Carto.FeatureLayer();
pFL.FeatureClass = pFW.OpenFeatureClass("seoul");
pFL.Name = pFL.FeatureClass.AliasName;

axMapControl1.AddLayer(pFL);

SHP 파일 하나를 추가하 위한 절차가 단순하지 않다는 것을 느낄 수 있습니다. 이는 ArcObjects가 매우 다양한 데이터 포맷을 지원하는 이유이며 로컬 데이터 뿐만 아니라 네트워크로부터 데이터를 받을 수 있도록 하기 위함입니다.

중요한 부분만을 짚어 설명한다면 6번 줄의 코드에서 OpenFromFile 매서드의 첫번째 인자 값인 d:/__data__ 라는 폴더를 작업공간으로 지정하겠다는 의미입니다. 즉, 추가하고자 하는 SHP 파일이 존재하는 경로입니다. 두번째 인자는 부모 윈도우의 핸들로써 NULL(0) 값으로 주어도 관계없습니다. 실제 SHP 파일명은 9번째 줄의 OpenFeatureClass 매서드의 인자로 지정하며 주의할 점은 확장자인 .shp를 생략하고 있다는 점입니다.

이런 절차로 생성한 ILayer 타입인 IFeatureLayer를 MapControl의 AddLayer 매서드의 인자로 주게 되어 추가할 수 있습니다. 결과는 아래와 같습니다.

사용자 삽입 이미지

레이어를 추가하자 마자 MapControl은 추가된 레이어의 전체를 도시하게 됩니다. 이상으로 간단히 SHP 파일을 ArcObjects에서 제공하는 MapControl에 추가하여 도시하는 방법을 살펴보았습니다.

[GIS] TXT 파일을 SHP 파일로 변환하기(Txt2Shp)

가끔 좌표 데이터와 속성 데이터가 텍스트 파일 형태로 되어 있는 경우가 있습니다. 하나의 예로 서울시의 버스 노선에 대한 데이터인데, 서울시에서 버스 노선에 대한 정보을 엑셀로 저장하여 사용자에게 제공하고 있습니다. 이 엑셀 파일을 쉽게 텍스트 파일로 저장할 수 있으며 다시 이 텍스트 파일을 읽어서 SHP 파일로 변환해 주는 툴입니다. 참고로 서울시 버스 노선에 대한 엑셀 파일은 아래 URL을 통해 다운로드 받을 수 있습니다.

위의 엑셀 파일을 엑셀 프로그램에서 읽어 [다른 이름으로 저장하기]에서 “텍스트 (탭으로 분리)” 형식으로 저장한 파일을 살펴보면 다음과 같습니다.

사용자 삽입 이미지

첫번째 줄에 필드 명 들이 텝문자로 구분되어 나타나고.. 이후의 줄은 각 필드에 대한 값이며 역시 텝으로 분리되어 있습니다. 이 텍스트 파일을 SHP 파일로 만드는 툴의 사용은 아래와 같습니다.

사용자 삽입 이미지

살짝 설명을 하면, Fields에 내보내고자 하는 항목을 체크하고 X, Y 좌표로써 사용할 필드를 지정해 줍니다. 위의 서울시 노선 데이터와 툴을 사용해 만든 SHP 파일을 ArcMap에서 열어 보면 아래 화면과 같습니다.

사용자 삽입 이미지
툴에서 지정한 속성과 좌표가 제대로 표시되는 것을 확인할 수 있습니다. 간단히 TXT 파일을 SHP 파일로 변환할 경우가 있을 때 요긴하게 사용하시길 바라면서 만들어 올려 봅니다. 참고로 이 프로그램은 GIS 오픈소스 라이브러리인 shapelib을 사용해 만들어 었습니다. shapelib의 자세한 사용법은 cafe.naver.com에서 오픈소스 라이브러리 카테고리를 참고하시면 도움을 받으실 수 있습니다.

[GIS] 표준 GML/KML/WMS/WFS에 대한 요약

자주 솔루션에 대한 요구사항으로 표준을 준수하는가… 라는 대응 항목 중 gml, kml, wms, wfs 이 있습니다. 새로운 gis 데이터 서버 개발에 이러한 표준 항목을 준수하기 위한 준비로 해당 항목에 대해 개략적으로 파악할 목적으로 작성한 내용입니다.

GML (Geography Markup Language)
http://www.opengeospatial.org/standards/gml

GML은 지리적 특성을 표현하기 위한 XML 문법입니다. GML은 인터넷에서 지리 데이터 송수신을 위한 개방된 교환 포맷임과 동시에 지리데이터 시스템을 위한 모델링 언어로써 사용됩니다. 문법에 기반한 대부분의 XML과 같이 두 부분으로 나눠져 있습니다.

  1. 문서의 설명 부분
  2. 실제 데이터 부분

GML 문서는 GML 스키마를 사용하여 기술됩니다. 이는 사용자와 개발자에게 포인트, 라인 그리고 폴리곤을 포함하는 일반적인 지리데이터셋을 설명할 수 있도록 합니다. 그러나 GML의 개발자들은 GML의 확장을 특화시키는 커뮤니티 지향적인 어플리케이션 스키마를 정의하기 위한 커뮤니티 작업을 상상합니다. 어플리케이션 스키마를 사용하여 사용자는 일반적인 포인트, 라인, 폴리곤을 길, 고속도로 그리고 다리로써 참조할 수 있습니다. 만약 커뮤니티의 모든 사람이 동일한 스키마를 사용하는 것에 대해 동의한다면 데이터를 쉽게 교환할 수 있으며 도로에 대한 라인을 표시할때 도로로 표시할 수 있게 됩니다.

KML(Keyhole Markup Language)
http://earth.google.com/intl/ko/userguide/v4/ug_kml.html
http://en.wikipedia.org/wiki/Keyhole_Markup_Language

KML은 현재 또는 미래의 웹 기반의 2차원과 3차원 브라우저에서 지리 데이터의 주기와 가시화를 위한 XML 기반의 스키마입니다. KML은 Google 어스, Google 지도 및 기타 응용 프로그램에 표시하기 위해 점, 선, 이미지, 다각형 및 모델과 같은 지형 기능을 모델링하고 저장하기 위한 XML 문법 및 파일 형식입니다. KML을 사용하여 Google 어스 및 Google 지도의 다른 사용자와 장소 및 정보를 공유할 수 있습니다. KML은 구글어스에서 사용하기 위해 개발되어졌으며 원래 이름은 Keyhole Earth Viewer였습니다. Keyhole 사에서 만들어졌으며 이 회사는 2004년에 구글로 인수되었습니다. 현재 KML은 오픈 GIS 컨소시엄의 국제 표준입니다. KML을 표현하고 편집할 수 있는 첫번째 프로그램은 구글어스이며 KML을 지원하는 툴이 다양하게 제공되고 있습니다.

참고 : KMZ 파일은 KML 파일의 압축된 버전입니다.

WMS(Web Map Service)
http://en.wikipedia.org/wiki/Web_Map_Service

WMS는 GIS 데이터베이스의 데이터를 사용하는 지도 서버로부터 생성되어진 지오레퍼런싱된 지도 이미지를 서비스하기 위한 표준 프로토콜입니다. 스펙은 1999년에 오픈 GIS 컨소시엄에서 처음 개발어 공개되었습니다.

WMS는 다양한 요청 타입을 가지고 있으며 WMS 서버는 다음 2가지 요청 타입을 제공해야 합니다.

  1. GetCapabilities : WMS에 대한 인자값과 이용가능한 지도 레이어 목록 반환
  2. GetMap : 해당 인자에 대한 지도 이미지

다음은 WMS 서버에서 옵션으로 제공하는 요청 타입입니다.

  1. GetFeatureInfo
  2. DescribeLayer
  3. GetLegendGraphic

WMS는 클라이언트가 인터넷으로부터 접근할 수 있는 지도와 GIS 데이터 포맷을 폭넓게 지원합니다. WMS를 지원하는 주요 상업적인 GIS 및 맵핑 소프트웨어는 다음과 같습니다.

  • Bentley Systems’s GIS products
  • ESRI’s ArcGIS products
  • MapInfo Professional
  • GeoMedia
  • Manifold System
  • Google Earth

또한 WMS를 지원하는 오픈소스 제품군은 다음과 같습니다.

  • Quantum GIS
  • uDig
  • OpenJUMP
  • MapGuide Open Source
  • NASA World Wind
  • GRASS GIS
  • gvSIG
  • OpenLayers

WFS(Web Feature Service)
http://en.wikipedia.org/wiki/Web_Feature_Service

WFS는 플랫폼에 독립적인 호출을 이용하여 웹에서 지리데이터에 대한 피쳐(feature)를 요청할 수 있는 인터페이스를 제공합니다. WMS가 단순히 지도 이미지만을 서비스하는데 반해 WFS는 공간 분석 등과 같은 연산을 위해 필요한 데이터를 서비스 합니다.

WFS의 서비스 기능은 다음과 같습니다.

  1. 공간 또는 비공간 조건에 기반한 피쳐(feature)의 쿼리
  2. 새로운 피쳐 생성
  3. 피쳐 삭제
  4. 피쳐 업데이트

가장 기본적인 WFS 서버는 반드시 공간 또는 비공간 조건에 기반한 피쳐의 쿼리 기능을 지원해야 하며, 트랜젝션 WFS(WFS-T)는 피쳐의 생성과 삭제 그리고 업데이트를 제공합니다.