[태그:] GIS

작성일자

[GIS] ArcObjects, 도형의 좌표 구하기

ArcObjects를 통해 SHP 파일을 읽어서 만든 레이어를 화면상에 표시한 후… 이 레이어를 구성하는 도형(Shape) 좌표를 구하는 방법을 정리해 보았습니다. ArcObjects에서 제공하는 Map 컨트롤을 통해 SHP 레이어가 하나 추가되었다고 가정하겠습니다. 참고로 이 포스트는 레이어가 폴리곤 도형으로 구성된 레이어를 대상으로 합니다.

가장 먼저 해야할 일은 지도 컨트롤이 가지고 있는 레이어를 가져오는 일입니다. 가져올때 도형 정보에 접근해야 하므로 IFeatureLayer 인터페이스 타입으로 가져와야 합니다.

ESRI.ArcGIS.Carto.IFeatureLayer pFL;
pFL = axMapControl1.get_Layer(0) as ESRI.ArcGIS.Carto.IFeatureLayer;

이렇게 가져온 IFeatureLayer 변수인 pFL을 통해 각각의 Feature를 얻어오기 위해 IFeatureClass와 IFeatureCursor를 이용합니다.

ESRI.ArcGIS.Geodatabase.IFeatureClass pFClass;
pFClass = pFL.FeatureClass;

ESRI.ArcGIS.Geodatabase.IFeatureCursor pFCursor;
pFCursor = pFClass.Search(null, false);

이제 레이어를 구성하는 Feature를 읽을 수 있는 준비가 되었습니다. 아래의 반복문을 통해 각각의 Feature를 순회하면서 원하는 좌표를 읽게 됩니다.

ESRI.ArcGIS.Geodatabase.IFeature pFeature;            
while((pFeature=pFCursor.NextFeature()) != null) {
    ....
}

아래의 코드가 위의 코드의 반복문 사이(while문)에 들어가는 코드입니다. 도형은 Point의 집합이므로 IPointCollection 인터페이스 타입이 필요하며 각 포인트 좌표값과 관련 정보를 얻기 위해 IEnumVertex 인터페이스가 필요합니다.

    ESRI.ArcGIS.Geometry.IPointCollection pPC;
    pPC = pFeature.Shape as ESRI.ArcGIS.Geometry.IPointCollection;

    ESRI.ArcGIS.Geometry.IEnumVertex pEV = pPC.EnumVertices;

    ESRI.ArcGIS.Geometry.IPoint pPt;
    int outPartIndex;
    int vertexIndex;
    listBox1.Items.Add("FID: " + pFeature.OID);
    for(int i=0; i<pPC.PointCount; ++i) {
        pEV.Next(out pPt, out outPartIndex, out vertexIndex);
        listBox1.Items.Add("    Part Idx :" + outPartIndex + 
            " Vertex Idx : " + vertexIndex + " (" + pPt.X + "," + pPt.Y + ")");
    }

비록 이 포스트의 글은 폴리곤 레이어에 대한 내용이지만 이 보다 간단한 포인트나 폴리라인에 대해서도 쉽게 좌표를 읽을 수 있습니다. 마지막으로 아래는 위의 코드들을 실행했을때의 결과 입니다.

사용자 삽입 이미지

작성일자

[GIS] ArcObjects, 속상값에 따라 다른 색상 지정(UniqueValueRenderer)

ArcObjects는 공간분석등을 포함한 지오프로세싱과 공간 데이터 편집에도 탁월한 기능을 제공하지만 공간 데이터의 표현에도 매우 뛰어난 기능을 제공합니다. 여기서는 공간 데이터에 대한 각 속성값에 따라 다른 그리기 심벌을 지정할 수 있는 UniqueValueRenderer에 대해 정리해 보도록 하겠습니다. UniqueValueRenderer는 특정 값에 대한 그리기 심벌을 정의할 수 있는 렌더러입니다. 이 포스트는 기본적으로 SHP 파일을 통한 레이어를 ArcObjects에서 제공하는 지도 컨트롤에 추가하고 있다고 가장합니다.

ESRI.ArcGIS.Carto.IFeatureLayer pFL;
pFL = axMapControl1.get_Layer(0) as ESRI.ArcGIS.Carto.IFeatureLayer;

ESRI.ArcGIS.Geodatabase.IFeatureClass pFClass;
pFClass = pFL.FeatureClass;

ESRI.ArcGIS.Geodatabase.IQueryFilter pQF;
pQF = new ESRI.ArcGIS.Geodatabase.QueryFilter();

ESRI.ArcGIS.Geodatabase.IFeatureCursor pFCur;
pFCur = pFClass.Search(pQF, false);

앞의 코드는 먼저 지도 컨트롤에 추가된 SHP 데이터에 대한 레어이의 IFeatureLayer 인터페이스를 구하고 있습니다. IFeatureLayer를 통해 속성값에 접근할 수 있는 IFeatureClass 인터페이스를 구할 수 있으며 이 IFeatureClass의 검색 메서드인 Search에 인자로 IQueryFilter 타입의 값을 넘겨주면 레이어의 전체 속성 레코드를 얻어오게 되고 그 결과는 IFeautreCursor 인터페이스를 통해 접근할 수 있습니다.

ESRI.ArcGIS.Display.IRandomColorRamp pRCR;
pRCR = new ESRI.ArcGIS.Display.RandomColorRamp();
pRCR.MinSaturation = 0;
pRCR.MaxSaturation = 100;
pRCR.MinValue = 0;
pRCR.MaxValue = 100;
pRCR.StartHue = 0;
pRCR.EndHue = 360;
pRCR.UseSeed = true;
pRCR.Seed = 87;

앞의 코드는 색상 심벌을 손쉽게 만들어 낼 수 있는 유틸리티 인터페이스인 IRandomColorRamp를 인스턴스화 하고 있습니다. 이 인터페이스는 사용자가 지정한 색조와 채도등의 값을 지정하고 지정된 값의 범위 내에서 임의의 난수 발생을 통해 색상 심벌을 뽑아 낼 수 있습니다. 참고로 지정할 수 있는 Hue의 최소와 최대값은 각각 0, 360이며 Saturation의 최소와 최대값은 각각 0, 100 그리고 Value의 최소, 최대값은 0, 100입니다. 또한 난수 발생에 대한 Seed값을 지정할 수 있는 속성을 제공합니다.

ESRI.ArcGIS.Carto.IUniqueValueRenderer pRender;
pRender = new ESRI.ArcGIS.Carto.UniqueValueRenderer();

pRender.FieldCount = 1;
pRender.set_Field(0, "SGG_NM");

ESRI.ArcGIS.Display.ISimpleFillSymbol pSFS;
pSFS = new ESRI.ArcGIS.Display.SimpleFillSymbol();
pSFS.Style = ESRI.ArcGIS.Display.esriSimpleFillStyle.esriSFSSolid;
pSFS.Outline.Width = 0.4;
pRender.DefaultSymbol = pSFS as ESRI.ArcGIS.Display.ISymbol;
pRender.UseDefaultSymbol = true;

앞의 코드는 이 포스트의 주인공인 UniqueValueRenderer를 생성하고 속성값을 설정하는 코드의 일부입니다. 앞서 언급했듯이 UniqueValueRenderer는 속성값에 기반하여 서로 다른 속성값에 따라 각기 다른 색상 심벌을 지정한다고 했습니다. 4번과 5번 코드는 속성값을 얻어올 필드의 개수의 필드명(여기서는 SGG_NM)을 지정합니다. 또한 지정하지 못한 속성값을 가지는 도형에 대해 기본적으로 사용할 그리기 심벌을 지정하는 코드가 7~12번의 코드입니다.

ESRI.ArcGIS.Geodatabase.IFeature pFeat;
long n = pFClass.FeatureCount(pQF);
ESRI.ArcGIS.Geodatabase.IFields pFields;
pFields = pFClass.Fields;
int iField = pFields.FindField("SGG_NM");
for(long i = 0; i

앞의 코드부분은 각 속성값에 대한 그리기 심벌(정확히 말하면 SimpleFillSymbol)을 지정합니다. 동일한 값의 중복되는 속성값에 대해 이중으로 그리기 심벌을 지정하지 않도록 valFound 변수를 통해 걸러내고 있습니다. 코드 라인 별로 살펴보면, 2번 코드는 레이어가 가지고 있는 전체 속성 레코드의 개수를 얻어와 n 변수에 저장하고 있습니다. 그리고 3번, 5번 코드는 속성값을 가져올 필드명에 대한 인덱스를 얻어와 iField에 저장합니다. 6번 코드의 for문은 모든 레코드에 대해서 앞서 지정한 필드의 값을 가져오고 그 필드 값에 대해 그리기 심벌을 생성하여 지정하고 있습니다. 15번 라인의 for문은 앞서 저장해둔 속성값들 중에 중복된 속성값에 대해서는 이중으로 심벌을 지정하지 않도록 하는 코드입니다. 여기서 눈여겨 봐야 할 점은 생성해 지정한 그리기 심벌에 대한 색상값을 지정하고 있지 않다는 점입니다. 일단 그리기 심벌을 생성해 두기만 하고 다음에 색상값을 지정합니다.


pRCR.Size = pRender.ValueCount;
bool bOK = true;
pRCR.CreateRamp(out bOK);
ESRI.ArcGIS.Display.IEnumColors pEnumClrs;
pEnumClrs = pRCR.Colors;
pEnumClrs.Reset();
for(int ny=0; ny

앞의 코드는 이전 부분에서 생성해둔 그리기 심벌을 IRandomColorRamp 인터페이스 타입으로 생성해둔 객체를 통해 색상을 얻어와 실제 지정하는 코드입니다. 앞서서 pRCR이라는 변수 명으로 IRandomColorRamp 인터페이스의 CoClass를 생성해 놓았습니다. 이 pRCR로부터 생성할 색상의 개수를 지정하고 CreateRamp 매서드를 호출하여 원하는 개수만큼 원하는 색상 계열로 색상을 생성해 두는 코드가 1번~3번 코드입니다. 이제 생성해 놓은 색상을 얻기 위해 Iterator 디자인 패턴의 방법을 통해 얻어 오는 부분이 4,5,6번 코드와 12번 코드입니다. 7번 코드의 for문은 고유한 필드값으로 저장된 그리기 심벌들을 가져와서 이 심벌의 색상을 pRCR에서 얻어온 색상으로 지정하기 위한 반복문입니다. 이 부분까지가 UniqueValueRenderer에 대한 필요한 속성값들을 설정하기 위한 모든 과정입니다. 이렇게 구성한 UniqueValueRenderer를 레이어에 지정하는 코드는 아래와 같습니다.


pRender.ColorScheme = "Custom";
pRender.set_FieldType(0, true);

ESRI.ArcGIS.Carto.IGeoFeatureLayer pGFL = 
    pFL as ESRI.ArcGIS.Carto.IGeoFeatureLayer;
pGFL.Renderer = pRender as ESRI.ArcGIS.Carto.IFeatureRenderer;

axMapControl1.ActiveView.Refresh();


레이어의 렌더러는 IGeoFeatureLayer 인터페이스를 통해 접근할 수 있으므로 4번~6번 코드와 같은 형태로 렌더러를 지정할 수 있습니다. 최종적인 실행 결과는 아래와 같습니다.


사용자 삽입 이미지
작성일자

[GIS] WKB(Well-Known Binary) Format

공간 데이터를 처리하는 DBMS는 포인트, 폴리라인, 폴리곤과 같은 공간 좌표 데이터를 효율적으로 저장하기 위해 BLOB 타입으로써 WKB라는 형태로 저장합니다. 개발하고 있는 지도 렌더러를 위한 맵엔진(코드명: Xr)에서 사용할 수치지도 데이터 포맷에 대해 고민하던 차에… GIS 데이터 표현의 범용 표준 중에 하나인 WKB에 대한 자료를 정리해 보았습니다.

WKB, 즉 Well-known Binary는 OpenGIS 스펙에 의해 정의된 지오메트리 값을 표현하기 위해 사용됩니다. KWB는 지오메트리 WKB 정보를 담고 있는 BLOB 타입의 값으로 표현되는 바이너리 스트림으로써 지오메트리 데이터를 상호간에 교환하기 위해 사용됩니다. WKB는 1바이트의 부호없는 정수와 4바이트의 부호없는 정수 그리고 8바이트의 배정밀 소수(IEEE 754)를 사용합니다. 여기서 1바이트는 8비트입니다.

예를 들면… WKT인 POINT(1 1)을 동일한 WKB로 변환해 본다면… 다음과 같은 총 21개로 구성된 연속된 바이트입니다.

0101000000000000000000F03F000000000000F03F

앞의 연속된 바이트를 해석 해보면… 다음과 같습니다.

사용자 삽입 이미지

즉, 첫번째 1바이트의 값이 1이면 바이트 정렬이 Little-Endian이며 0이면 Big-Endian이라는 의미입니다. 그리고 다음 2바이트는 지오메트리 타입으로써 현재까지 제안된 타입은 총 7개로써 Point, LineString, Polygon, MultiPoint, MultiLineString, MultiPolygon 그리고 GemetryCollection입니다. 끝으로 포인트의 좌표인 X와 Y 값이 각각 8바이트씩 오게됩니다.

앞의 예는 Point 타입에 대한 예이며 다른 지오메트리 타입에 대해서, 보다 직관적이고 효율적으로 파악하기 위해 C언어의 구조체 형식으로 보면 다음과 같습니다.

// Basic Type definitions
// byte : 1 byte
// uint32 : 32 bit unsigned integer  (4 bytes)
// double : double precision number (8 bytes)
// Building Blocks : Point, LinearRing
 
Point {
    double x;
    double y;
};

LinearRing   {
    uint32 numPoints;
    Point points[numPoints];
};

enum wkbGeometryType {
    wkbPoint = 1,
    wkbLineString = 2,
    wkbPolygon = 3,
    wkbMultiPoint = 4,
    wkbMultiLineString = 5,
    wkbMultiPolygon = 6
};

enum wkbByteOrder {
    wkbXDR = 0, // Big Endian
    wkbNDR = 1 // Little Endian
};

WKBPoint {
    byte byteOrder;
    uint32 wkbType; // 1=wkbPoint
    Point point;
};

WKBLineString {
  byte byteOrder;
  uint32 wkbType; // 2=wkbLineString
  uint32 numPoints;
  Point points[numPoints];
};
 
WKBPolygon {
  byte byteOrder;
  uint32 wkbType; // 3=wkbPolygon
  uint32 numRings;
  LinearRing rings[numRings];
};

WKBMultiPoint {
  byte byteOrder;
  uint32 wkbType; // 4=wkbMultipoint
  uint32 num_wkbPoints;
  WKBPoint WKBPoints[num_wkbPoints];
};

WKBMultiLineString {
  byte byteOrder;
  uint32 wkbType; // 5=wkbMultiLineString
  uint32 num_wkbLineStrings;
  WKBLineString WKBLineStrings[num_wkbLineStrings];
};
 
wkbMultiPolygon {
  byte byteOrder;
  uint32 wkbType; // 6=wkbMultiPolygon
  uint32 num_wkbPolygons;
  WKBPolygon wkbPolygons[num_wkbPolygons];
};
 
WKBGeometry  {
  union {
    WKBPoint point;
    WKBLineString linestring;
    WKBPolygon polygon;
    WKBMultiPoint mpoint;
    WKBMultiLineString mlinestring;
    WKBMultiPolygon mpolygon;
  }
};

바이너리는 크기면에서 데이터를 매우 최적화된 형태로 저장할 수 있는 포맷입니다. GIS 엔진을 설계하고 개발할때 자체 엔진에 특화된 데이터 구조를 설계하는 것 대신에 이러한 표준을 수용할 경우 얻을 수 있는 이점은 매우 크다고 할 수 있겠습니다.

이 글의 일부는 다음 URL에서 참고하였습니다.
http://dev.mysql.com/doc/refman/5.0/en/gis-wkb-format.html

작성일자

[GIS] ArcObjects, 지도 레이어 그리기 심벌 지정하기(SimpleRenderer)

ArcObjects를 이용하여 Map 컨트롤에 지도레이어를 추가하게 되면, 지도 레이어를 그리기 위해 사용되는 그리기 기능으로 기본적으로 SimpleRenderer라는 렌더러가 할당됩니다. 이 렌더러는 해당 레어어를 구성하는 모든 도형에 대해 동일한 그리기 심벌을 이용해 그리게 됩니다. 예를 들어서 폴리곤 지도 레이어의 경우 채움색과 외곽선색 등이 모든 도형에 대해 동일하게 적용되어 그려지게 됩니다.

ArcObjects가 지원하는 렌더러는 매우 다양한데… ClassBreaksRenderer나 DotDensityRenderer, UniqueValueRenderer 등이 있습니다. 이 포스트에서는 렌더러 중 가장 쉬우며 기본적으로 설정되어 있는 SimpleRenderer의 심벌을 지정하여 사용자가 원하는 색상으로 지도 레이어를 그리는 방법에 대해 정리해 보았습니다.

먼저 폴리곤 지도 레이어가 Map 컨트롤에 하나 추가되어 있다가 가정하겠습니다. 그리고 이 지도 레이어의 채움색을 초록색으로 표시하고 외곽선 색을 검정색으로 표시해 보도록 하겠습니다.

ESRI.ArcGIS.Carto.IGeoFeatureLayer pGFL;
pGFL = axMapControl1.get_Layer(0) as ESRI.ArcGIS.Carto.IGeoFeatureLayer;

ESRI.ArcGIS.Carto.ISimpleRenderer pSR;
pSR = pGFL.Renderer as ESRI.ArcGIS.Carto.ISimpleRenderer;

앞의 코드는 가장 먼저 추가한 폴리곤 지도 레이어를 IGeoFeatureLayer 타입으로 가져옵니다. 이 타입을 통해 렌더러 객체에 접근할 수 있기 때문입니다. SHP 파일을 통해 새로 추가한 지도 레이어의 그리기 렌더러는 기본적으로 ISimpleRenderer가 지정되므로 이 렌더러를 가져와 pSR이라는 변수에 담아 두고 있습니다.

ESRI.ArcGIS.Display.IRgbColor pFillRGB;
pFillRGB = new ESRI.ArcGIS.Display.RgbColor();
pFillRGB.Red = 0;
pFillRGB.Green = 255;
pFillRGB.Blue = 0;

앞의 코드는 채움 색상을 지정합니다. 앞서 언급했듯이 채움색으로써 녹색(RGB(0,255,0))을 지정하고 있습니다.

ESRI.ArcGIS.Display.IRgbColor pLineRGB;
pLineRGB = new ESRI.ArcGIS.Display.RgbColor();
pLineRGB.Red = 0;
pLineRGB.Green = 0;
pLineRGB.Blue = 0;

ESRI.ArcGIS.Display.ISimpleLineSymbol pSLS;
pSLS = new ESRI.ArcGIS.Display.SimpleLineSymbol();
pSLS.Color = pLineRGB;

앞의 코드는 폴리곤의 외곽선을 지정하기 위한 코드입니다. 외곽선은 ISimpleLineSymbol이며 라인의 폭(Width 프로퍼티)와 라인의 스타일(Style 프로퍼티) 그리고 색상(Color 프로퍼트)를 저장할 수 있습니다. 여기서는 검정색 색상으로 외곽선의 색을 지정하기 위고 있습니다.

ESRI.ArcGIS.Display.ISimpleFillSymbol pSFS;
pSFS = new ESRI.ArcGIS.Display.SimpleFillSymbol();

pSFS.Color = pFillRGB;
pSFS.Outline = pSLS;

앞의 코드는 앞서 만들어 놓은 채움색과 외곽선에 대한 객체를 실제 폴리곤 도형 지도 레이어에 적용하기 위한 채움 스타일에 대한 타입인 SimpleFillSymbol을 생성하고 있습니다. 이 SimpleFillSymbol 타입의 객체의 변수명은 pSFS라고 해 놓았고 앞서 구한 렌더러의 Symbol 속성에 할당해주면 우리가 원하는 결과를 얻을 수 있게 됩니다.

pSR.Symbol = pSFS as ESRI.ArcGIS.Display.ISymbol;
axMapControl1.ActiveView.Refresh();

앞의 코드가 바로 렌더러의 Symbol 속성에 심벌을 지정하는 코드입니다. 실행해 보면 원하는 색상 심벌로 지도 레이어가 그려지는 것을 확인할 수 있습니다.

사용자 삽입 이미지