[GIS] SHP을 SQL 문으로 변환해 주는 툴 – SHP2SQL

SHP 파일의 좌표 정보와 속성 정보를 DBMS에 테이블 구조를 생성하고 값을 추가하는데 필요한 SQL문을 생성해 주는 툴입니다. 현재는 포인트 타입에 대한 SHP 파일에 대해서 적용할 수 있습니다. 먼저 실행 화면은 아래와 같습니다.

사용자 삽입 이미지
위의 UI 화면에서 Export Fields에 DBMS로 내보내고자 하는 필드를 선택하면 됩니다.. 필요없는 필드값에 대해서도 제외시켜도 됩니다. 이 기능을 이용해서 DBMS에 따라 필드명으로써 사용할 수 없는 필드는 이곳을 통해 제외시켜줄 수 도 있습니다.. 위의 화면 예시처럼 입력값들을 지정한 후 변환 버튼을 누르게 되면 다음과 같은 sql 문이 담긴 텍스트 파일이 만들어 집니다..

사용자 삽입 이미지
테이블을 생성하고 생성된 테이블에 속성값을 추가하기 위한 Ansi SQL 구문으로 구성됩니다.. 이 내용을 수정해 줄 필요(테이블명)도 있을 것입니다..  이렇게 생성된 SQL 파일을 각 DBMS에서 불러와 실행해주면 간단하게 DBMS에 SHP 파일에 대한 테이블이 생성됩니다.. MySQL의 경우 soruce 명령어를 이용해서 외부의 sql 문이 담긴 파일을 실행할 수 있습니다.

이 툴의 실행을 위해서는 최신 버전의 듀라맵을 설치하셔야 합니다. 다음 url을 통해 듀라맵을 설치하시기 바랍니다.

http://www.gisdeveloper.co.kr/notice/574

소스 코드와 컴파일된 실행 파일 모두를 제공해 드리니.. 개발자 분이라면 목적에 맞게 커스터마징해 사용하시길 바랍니다..

아래의 이미지는 위에서 만들어진 sql 구문이 담긴 text 파일을 MySQL에서 실행시켜 만들어진 테이블을 조회해 본 화면입니다.

사용자 삽입 이미지
위의 이미지를 보시면.. X, Y 필드명으로 POI의 좌표가 자동으로 추가된 것을 확인할 수 있습니다.

[GIS] SHP 파일을 TXT 파일로 변환해 주는 툴 – Shp2Txt

SHP 파일은 바이너리 파일인지라.. 간단히 읽기에는 어려움이 있습니다. 해서.. SHP 파일을 일반적인 텍스트 파일로 변환하여 좌표나 속성값을 쉽게 읽거나.. 아니면 이 텍스트 파일을 읽어 자신이 원하는 또 다른 형식으로 변환하여 활용할 수 있도록 개발자를 돕는 툴입니다. 실행은 간단합니다. 아래처럼 입력할 SHP 파일을 지정하고 저장할 TXT 파일 이름을 지정해 주고 변환 버튼을 클릭하면 됩니다.

사용자 삽입 이미지
만들어진 txt 파일을 살펴보면 아래와 같습니다. 처음 도형에 대한 좌표 리스트가 나오며.. 이 도형에 대한 속성값이 ‘|’ 문자로 구분되어 출력됩니다. 빨간색으로 표시된 부분이 바로 도형에 대한 정보이며 파랑색으로 표시된 부분이 해당 도형의 속성에 대한 정보입니다.

사용자 삽입 이미지
끝으로… 이 프로그램은 듀라맵을 기반으로 수행됩니다. 실행을 위해 먼저 듀라맵을 설치(http://www.gisdeveloper.co.kr/notice/574)하셔야 합니다. 이 툴에 대한 소스 파일과 실행 파일을 모두 올려 드립니다. 개발이 가능하신 분이라면 자신이 원하는 텍스트 형식으로 저장할 수 있도록 수정이 가능할 것입니다..

[GIS] DuraMap-Xr, 그리드 레이어 응용(IDW 보간)

그리드 레이어에 대한 간단한 소개에 이어… 그럼 어떻게 그리드 레이어를 사용하는지를 하나의 예를 통해 살펴보도록 하겠습니다. IDW 보간 기법인데요. 공간상에 분포된 기지점(이미 알고 있는 공간상의 위치와 속성값)을 이용해서 미지점에 대한 속성값을 보간해 보는 알고리즘이 바로 IDW 입니다. IDW에 대한 설명과 수학적 공식은 이 블로그에서 IDW로 검색해 살펴보시기 바랍니다. 이후부터는 IDW에 대한 알고리즘을 이해하고 있다는 전제 아래에서 설명하겠습니다만……. 이해하지 못한다고 해도 그리드 레이어를 응용하는 것에 대한 감은 충분히 익히실 수 있을거라 판단됩니다. 뭐니 뭐니해도… 코딩을 통해 살펴 보는 것이 머리도 맑고.. 정신 건강에도 좋을듯합니다.

사용자 삽입 이미지
각 버튼에 대해 간단히 설명을 드리면… Add Points 버튼은 기지점에 대한 포인트를 추가하는 버튼이며.. Add Grid Layer는 그리드 레이어를 추가하는 버튼이며.. IDW Analysis는 기지점과 추가한 그리드 레이어를 이용해 IDW 보간을 수행하는 버튼입니다. 끝으로 Color Setting은 분석 결과 그리드를 표현하는 색상을 변경하는 버튼입니다. 먼저 Add Points 코드에 대해 간단히 살펴보도록 하겠습니다.

bool bOK = axXr1.Layers.AddShapeMapLayer("points", "d:/__data__/Points.shp");
if(bOK) 
{
    axXr1.WaitForAllConnections();

    bOK = axXr1.Labels.AddLabel("points", "{Value2}");
    if(bOK)
    {
        XrMapLib.Label lbl = axXr1.Labels.GetLabel("points", "{Value2}");
        lbl.Font.Name = "맑은 고딕";
        lbl.Font.Color = (uint)ColorTranslator.ToOle(Color.White);
        lbl.Font.Size = 21;
        lbl.Effect.OutlineColor = (uint)ColorTranslator.ToOle(Color.DarkGray);
        lbl.Effect.Enable = true;
        lbl.Effect.OutlineWidth = 3;
    }

    axXr1.ZoomFullExtent();
}

네..! points라는 레이어 명으로 지정된 shp 파일을 레이어로 추가하는 코드입니다. 덧붙여.. 이 shp 파일에는 Value2라는 수치값 속성이 존재합니다. 즉, 이 shp 파일에는 공간상의 (x, y) 지점과 해당 수치값을 담고 있는 공간 데이터입니다. 어떤 값이 담겨 있는지를 사용자에게 보여주기 위해 9번~15번 코드를 통해 라벨을 표시하도록 했습니다. 결과는 아래 이미지와 같습니다.

사용자 삽입 이미지
이제 다음으로 Add Grid Layer 버튼에 대한 코드를 살펴보도록 하겠습니다. 코드는 다음과 같습니다. 이미 그리드 레이어에 대한 소개 글에서 살펴보았으므로 역시 간단하게 설명하도록 하겠습니다.

Extent MBR = axXr1.MapExtent;

int ColumnCount = 500;
double Resolution = (MBR.MaxX - MBR.MinX) / ColumnCount;
int RowCount = (int)((MBR.MaxY - MBR.MinY) / Resolution);
string ConnectionString =
        "OnMemory=1;MinX=" + MBR.MinX + ";MinY=" + MBR.MinY +
        ";Resolution=" + Resolution +
        ";ColumnsCount=" + ColumnCount + ";RowsCount=" + RowCount + 
        ";NullValue=-9999";

bool bOK = axXr1.Layers.AddGridMapLayer("grid", ConnectionString);
if (bOK)
{
    axXr1.WaitForAllConnections();
    axXr1.ZoomFullExtent();
}

위의 코드는 그리드 레이어를 추가하는 코드로써 먼저 그리드 레이어에 대한 정보를 수집하는 코드가 실행됩니다. 1번 코드는 현재 화면에 대한 MBR을 얻는 코드이고 3번 코드는 그리드 레이어의 가로 셀 개수를 500개로 지정하기 위한 변수입니다. 4번 코드는 셀의 해상도 값입니다. 세로 셀 개수는 앞서 구한 MBR과 가로 셀 개수 그리고 해상도를 통해 계산하며 5번 코드에 해당됩니다. 이렇게 지정한 값을 통해 6번 코드에서 그리드 레이어의 연결문자열을 지정합니다. 메모리 상에서만 생성되며 파일을 저장하지 않습니다. 보다 자세한 내용은 그리드 레이어 소개 글을 참조하시기 바랍니다. 아무런 데이터 값도 지정하지 않은 그리드 레이어를 추가했으므로 화면상의 결과는 달라지지 않습니다. 이제 IDW Analysis 버튼에 대한 코드를 살펴보겠습니다.

ShapeMapLayer pointsLyr = axXr1.Layers.GetLayerAsShapeMap("points");
GridMapLayer gridLyr = axXr1.Layers.GetLayerAsGridMap("grid");
if (pointsLyr != null && gridLyr != null)
{
    AttributeTable attTbl = pointsLyr.AttributeTable;
    ShapeTable shpTbl = pointsLyr.ShapeTable;

    int FieldIndex = attTbl.FieldSet.GetFieldIndex("Value2");
    int RowCount = shpTbl.RowCount;
    int nCellRowCount = gridLyr.Cells.RowCount;
    int nCellColumnCount = gridLyr.Cells.ColumnCount;
    double GridMinX = gridLyr.MBR.MinX;
    double GridMinY = gridLyr.MBR.MinY;
    double Resolution = gridLyr.Cells.Resolution;
    double power = 1.3;

    for (int iRow = 0; iRow < nCellRowCount; ++iRow)
    {
        for (int iColumn = 0; iColumn < nCellColumnCount; ++iColumn)
        {
            double CellX = GridMinX + (double)iColumn * Resolution 
                + Resolution / 2.0;
            double CellY = GridMinY + (double)iRow * Resolution
                + Resolution / 2.0;

            double A = 0.0;
            double B = 0.0;

            for (int FID = 0; FID < RowCount; ++FID)
            {
                ShapeRow sRow = shpTbl.GetRow(FID);
                AttributeRow aRow = attTbl.GetRow(FID);

                if (sRow.Load() && aRow.Load())
                {
                    Coordinate coord = sRow.GetVertex(0);
                    double value = aRow.GetValueAsFloat(FieldIndex);

                    double W = 1.0 / Math.Pow(Math.Sqrt((coord.X - CellX) 
                        * (coord.X - CellX) + (coord.Y - CellY) 
                        * (coord.Y - CellY)), power);
                    A += (W * value);
                    B += (W);
                }
            }

            gridLyr.Cells.SetValue(iRow, iColumn, A/B);
        }
    }

    gridLyr.Cells.Normalize();
    gridLyr.EnableCacheBitmap = true;
    axXr1.Update();
}

Shape 지도 레이어의 공간 데이터(x,y 좌표)와 속성 데이터에 접근하기 위해 ShapeTable과 AttributeTable 객체를 얻어옵니다. 그리고 이 객체를 통해 IDW 계산식을 수행합니다. IDW에 대한 계산식은 이 블로그에서 IDW로 검색해 보시기 바랍니다. 소스 코드는 길지만.. IDW 계산을 위한 코드가 대부분을 차지 하므로 이 부분에 대한 자세한 설명은 생략하도록 하겠습니다. IDW를 계산하기 위한 power 값을 1.3으로 지정했다는 것만 언급하겠습니다. 또한 위의 코드는 IDW 보간 연산을 빠르게 하기 위한 전처리를 하지 않고 최대한 코드를 단순하게 작성하여 속도가 제법 느립니다만.. 약간의 전처리 연산을 선행하면 훨씬 빠른 연산 속도를 얻을 수 있습니다. 결과는 아래와 같습니다.

사용자 삽입 이미지
네.. 결과를 보시면 기지점을 이용하여 그 이외의 미지점에 대한 값을 보간하여 각 셀에 보간된 값을 담고 있습니다.그리드 레이어는 기본적으로 위의 결과처럼 Greyscale로 셀값을 표현하지만.. 사용자가 원하는 색상으로도 표현이 가능하며 이 기능이 바로 Color Setting 버튼입니다. 코드를 살펴보면…

GridMapLayer gridLyr = axXr1.Layers.GetLayerAsGridMap("grid");
if (gridLyr != null)
{
    gridLyr.ColorTable.Clear();

    int cntColorStep = 64;
    float R1 = 0, G1 = 255, B1 = 0;
    float R2 = 255, G2 = 255, B2 = 0;
    float R3 = 255, G3 = 0, B3 = 0;

    float R, G, B;
    float StepR, StepG, StepB;

    R = R1; 
    G = G1; 
    B = B1;

    StepR = (R2 - R1) / cntColorStep; 
    StepG = (G2 - G1) / cntColorStep;
    StepB = (B2- B1) / cntColorStep; 

    for (int i = 0; i < cntColorStep; i++)
    {
        R += StepR;
        G += StepG;
        B += StepB;
        gridLyr.ColorTable.AddColor(255, (byte)R, (byte)G, (byte)B);
    }

    R = R2;
    G = G2;
    B = B2;

    StepR = (R3 - R2) / cntColorStep;
    StepG = (G3 - G2) / cntColorStep;
    StepB = (B3 - B2) / cntColorStep;

    for (int i = 0; i < cntColorStep; i++)
    {
        gridLyr.ColorTable.AddColor(255, (byte)R, (byte)G, (byte)B);
        R += StepR;
        G += StepG;
        B += StepB;
    }

    gridLyr.EnableCacheBitmap = false;
    gridLyr.EnableCacheBitmap = true;
    axXr1.Update();
}

색상을 구성하는 핵심 키 포인트는 먼저 기존의 Greyscale 색상들을 제거하기 위해 4번 코드에서 처럼 그리드 레이어의 ColorTable의 Clear 매서드를 호출하고 이 ColorTable에 사용자가 원하는 색상을 ColorTable의 AddColor 매서드를 통해 추가하기만 하면됩니다. 결과는 아래와 같습니다.

사용자 삽입 이미지

[GIS] DuraMap-Xr, 그리드 레이어 소개

DuraMap-Xr은 DEM나 그리드 데이터를 생성하고 읽어들이면서 그리드 데이터의 값을 수정할 수 있는 그리드 레이어 기능을 제공합니다. 그리드 레이어로써 쉽게 ArcGIS에서 제공하는 FLT 파일을 DuraMap-Xr에서 읽어 표시하는 코드는 아래와 같습니다.

bool bOK = Xr1.Layers.AddGridMapLayer("grid", 
        "FileName=d:/__data__/flt/rastert_r151.flt;OnMemory=0");
if(bOK) 
{
    Xr1.WaitForAllConnections();
    GML.EnableCacheBitmap = true;
    Xr1.ZoomFullExtent();
 }

코드를 살펴보면…. 먼저 1번 코드에서 AddGridMapLayer 매서드를 통해 그리드 레이어를 추가하며 이 매서드의 첫번째 인자는 레이어의 이름입니다. 중요한 것이.. 두번째 인자인데요. 오픈할 파일명과 이 파일에 대한 데이터를 메모리로 올릴지에 대한 여부를 지정하고 있습니다. 파일명을 지정하기 위해 FileName= 다음에 파일 경로를 지정하면 되고.. OnMemory= 뒤의 값이 1이면 메모리에 올린다는 의미이고 0이면 메모리에 올리지 않고 파일에서 바로 바로 데이터 값을 읽겠다는 의미입니다. 위의 코드에 대한 실행 결과는 아래와 같습니다.

사용자 삽입 이미지
실제로는 이렇게 이미 존재하는 파일을 읽는것뿐만이 아니라 새롭게 그리드 파일을 생성하는 경우가 빈번합니다. 아래의 코드는 새로운 그리드 데이터를 생성하는 API에 대한 예제입니다.

bool bOK = Xr1.Layers.AddGridMapLayer("grid", 
"FileName=d:/__data__/flt/test.flt;OnMemory=1;MinX=1000;MinY=1000;Resolution=10;ColumnsCount=100;RowsCount=100;NullValue=-9999");

if(bOK)
{
    Xr1.WaitForAllConnections();
}

1번 코드에서 AddGridMapLayer의 두번째 인자가 그리드 데이터를 생성할때 필요한 인자값들입니다. 각 인자값에 대해 간단히 설명을 하면 다음과 같습니다.

  • FileName 새롭게 생성할 그리드 데이터를 저장할 파일명이며, 이 FileName이 지정되지 않을 경우 메모리 상에서만 생성되고 파일로 저장하지 않겠다는 의미로 해석됩니다.
  • OnMemory 새롭게 생성할 그리드 데이터를 메모리 상에 올릴지에 대한 여부로 1이면 메모리로 올리겠다는 의미이고 0이면 올리지 않겠다는 의미입니다. 만약 값이 0일 경우 FileName을 만든시 지정해 주어야 합니다.
  • MinX 그리드 데이터에 대한 영역(MBR)에서 좌측하단에 대한 X 좌표입니다.
  • MinY 그리드 데이터에 대한 영역(MBR)에서 좌측하단에 대한 Y 좌표입니다.
  • Resolution 그리드 데이터에 대한 셀 해상도입니다.
  • ColumnCount 그리드 데이터에 대한 가로측 방향의 셀 개수입니다.
  • RowCount 그리드 데이터에 대한 세로측 방향의 셀 개수입니다.
  • NullValue 그리드 데이터에 대한 셀에 값을 지정하지 않았을 경우 할당될 값입니다.

새롭게 생성된 그리드 데이터는 값이 전혀 할당되지 않았으므로 화면상에 아무것도 표시되지 않지만 분명히 주어진 공간상에 위치하고 있습니다. 이렇게 생성된 그리드 데이터에 대한 셀에 값을 읽고 쓰는 코드는 아래와 같습니다.

GridMapLayer GML = Xr1.Layers.GetLayerAsGridMap("grid");
GridCells Cell = GML.Cells;

Cell.SetValue(100, 200, 45364.453);
double value = Cell.GetValue(100, 200);

먼저 레이어의 이름으로 그리드 레이어 객체를 얻어오고 그리드 레이어의 셀값을 담고 있는 GridCells 객체를 얻어옵니다. 이렇게 얻어온 GridCells 객체를 통해 값을 설정하고 얻어올 수 있습니다.

이상으로 간단하게 그리드 레이어에 대한 소개를 마치고 다음으로 그리드 레이어를 응용해 공간상의 포인트 데이터를 이용해 IDW 분석을 하는 예제를 통해 그리드 레이어를 활용하는 것에 대해 살펴보도록 하겠습니다.