On Spatial

OGC SLD Converter for ArcGIS

2012-02-26T02:58:00.002+09:00

자체 프로젝트로 진행된 ArcMap에서 설정한 레이어의 렌더러 + 심볼을 OGC SLD 포맷으로 변환하는 도구 구현 결과물 ArcObjects의 Symbol 인터페이스 및 SLD에 대해 다시한번 공부할 수 있어서 의미 있는 프로젝트였음 ▣ 관련 내용 - ArcGIS to SLD(Styled Layer Descriptor) 기능정의 - SLD(Styled Layer Descriptor) Converter & Editor? - [GeoTools]Filter Encoding Quick Reference - GeoServer & SLD(Styled Layer Descriptor) ▣ OGC SLD Converter for ArcGIS - GeoServer : SLD 1.0.0 + Vendor

ArcGIS to SLD(Styled Layer Descriptor) 기능정의

2012-02-22T09:00:00.000+09:00

ArcMap의 레이어를 SLD(Styled Layer Descriptor) 파일로 내보내기 위한 도구 개발 SLD 변환도구 이전글: SLD(Styled Layer Descriptor) Converter & Editor? 참조 ▣ SLD 지원 대상 - GeoServer : SLD 1.0.0 + Vendor Option - ArcGIS 10 + .NET C# + Custom ICommand 구현 ▣ GeoServer에서의 SLD 지원 현황 - 현재 GeoSErver에서 SLD 인코딩은 ISO-8859-1을 사용하며 utf-8로 저장 후 GeoServer Administrator 페이지에서 Style을 수정하는 경우 헤더 부분 및 한글이 깨지는 현상 발생 - 그러나 저장 후 GeoServer

[GeoTools]Filter Encoding Quick Reference

2012-02-22T08:30:00.000+09:00

SLD Converter 구현 중에 Filter Encoding이 포함되어 있어 간략하게나마 Filter Encoding을 도식화해 보았다. SLD 변환도구는 이전글: SLD(Styled Layer Descriptor) Converter & Editor? 참조 현재 OGC Filter Encoding 표준은 2.0 버전이 최신이지만 GeoTools에서는 Filter 1.0, Filter 1.1을 지원하고 있다. ▣ Filter ▣ BinaryLogicOperator ▣ BinaryComparisonOperator ▣ SpatialOperator ▣ 참고 - http://docs.geotools.org/latest/userguide/library/opengis/filter.html

GeoServer & SLD(Styled Layer Descriptor)

2012-02-22T08:00:00.000+09:00

현재 GeoTools와 GeoServer는 공식적으로는 SLD version 1.0.0을 지원하고 있다고 하지만 SLD 1.0, 1.1, GeoServer옵션(Vendor 옵션)을 함께 사용하고 있다. SLD Converter 구현을 위해서 간략하게 SLD를 정리해 보았다. SLD 변환도구는 이전글: SLD(Styled Layer Descriptor) Converter & Editor? 참조 ▣ SLD의 기본 구조 예) <?xml version="1.0" encoding="utf-8"?> <sld:StyledLayerDescriptor xmlns="http://www.opengis.net/sld" xmlns:sld="http://www.opengis.net/sld" xmlns:ogc="http://

[ArcGIS]FeatureClass XY Domain..

2012-02-21T23:00:00.001+09:00

> XY Domain에 대한 자료를 찾던 중 알게 된 내용을 정리하였습니다. 보다 자세한 내용을 알고 싶었으나, 짧은 영어로 인한 독해의 압박에 그만.. 혹시 한글로 된 기술문서를 가지고 계시거나 사이트를 아시는 분은 코멘트 부탁 드립니다. ㅜㅜ << XY Domain >> > ShapeFile 의 경우에는 FeatureClass의 XY Domain이 없어 Edit 시에도 범위의 구애를 받지 않는다. 하지만 PersonalGDB(mdb), SDE의 경우에는 XY Domain의 범위를 관리한다. > 위 그림은 ArcCatalog에서 FeatureClass의 속성 창에 표시되는 Domain 정보이다. ■ PersonalGDB(mdb) : mdb파일을 MS Access로 연다. >

[ArcObjects]두 FeatureClass의 XY Domain 범위 합치기..

2012-02-21T21:42:00.000+09:00

> 두 개의 A, B FeatureClass를 Merge하여 하나의 FeatureClass로 만드는 과정에서.. > A의 Fields를 사용하여 새로운 FeatureClass를 만들고 A와 B의 FeatureClass의 Geometry를 Insert 하는 과정에서 A의 SpatialReference에 정의된 XY Domain의 범위에서 벗어난 Geometry가 B에 존재한다면 “The coordinates or measures are out of bounds.” 이라는 오류가 발생했다. > 물런 IBasicGeoprocessor 인터페이스의 Merge 메소드를 사용하면 별 문제 없이 작업이 가능하겠지만 단순한 Merge 작업이 아닌 중간에 어떠한 프로세스가 들어간다면 어쩔 수없이 수동으로

[ArcObjects]ArcSDE Connection 정보 저장…

2012-02-20T16:38:00.000+09:00

>> ArcObjects로 프로젝트를 진행하던 도중.. FeatureClass나 Table의 정보를 저장하였다가 불러와서 사용해야 하는 기능이 있었다. 로컬 ShapeFile이나 Personal GDB, File GDB의 경우에는 파일의 Path를 저장하여 나중에 해당 Path로 WorkSpace를 잡아주면 되는데… 문제는 SDE상에 올라가 있는 DataSet이다. >> SDE에 올라가있는 데이터의 경우에 필요한 것이 서버정보(IP), 인스턴스(esri_sde), USER, PASSWORD 인데.. 다른건 모두 저장하면 되는데.. 문제는.. PassWord 였다.. 사용자에게 입력받기도 그렇고.. 그렇다고 매번 데이터를 사용할때 마다 PassWord를 입력받기도 그렇고.. >> 일단 다음 코드로 SDE

SLD(Styled Layer Descriptor) Converter & Editor?

2012-02-18T18:13:00.000+09:00

Styling은 정말 손이 많이 가는 작업이다. 아래 도구들을 하나씩 사용해 보면서 부족한 점들을 나열해 보면 - 여러 레이어를 쉽게 스타일링 작업하기에는 불편함 - OpenGeo Suite의 경우 Web Browser에서 배치로 레이어를 Import하면서 Default SLD 를 생성해 주고 이를 GeoExplorer를 이용해서 하나씩 편집해 나갈 수 있어 그나마 편리함 - uDig은 Desktop 환경이기 때문에 스타일링 작업은 편리하지만 일일이 하나씩 SLD를 내보내기해야 하는 불편함 조만간 데스크탑 환경에서 레이어 스타일링 후 배치로 SLD 파일로 내보내는 기능을 만들어볼까 한다. 대상 소프트웨어는 - ArcGIS Desktop - Free 또는 Open Source 중 Map

[PostGIS]Proximity & Distance Funtion

2012-02-15T23:03:00.000+09:00

OSGeo(Open Source GeoSpatial) 한국어 지부 활동을 하면서 Open Source Winter School 프로그램 중 PostGIS 강사로 참여하게 되었습니다. 1차 강의를 마친 후 메일로 질문도 보내 주시고 합니다. 블로그를 통해 질문에 답하는 자리를 만들어볼까 합니다. 혹시 2차 교육에 관심 있으신 분들은 아래 URL에서 확인바랍니다. * 제1차 GeoSpatial Open Source Winter School : 2012년 2월 9일 ~ 10일 * 제2차 GeoSpatial Open Source Winter School : 2012년 2월 28일 ~ 29일 Geometry의 거리(Distance)는 Geometry의 공간좌표체계(spatial reference system

[PostGIS]ERROR: LWGEOM_estimated_extent

2012-01-31T23:50:00.001+09:00

GeoServer 등 GeoTools를 사용하는 프로그램의 경우 PostGIS Data Store를 사용하여 레이어의Extent를 구할 때 다음과 같은 오류가 발생하는 경우가 있다. ▣ GeoTools - PostGIS Map<String, Serializable> params = new HashMap<String, Serializable>(); params.put(JDBCDataStoreFactory.DBTYPE.key, "postgis"); params.put(JDBCDataStoreFactory.HOST.key, "localhost"); params.put(JDBCDataStoreFactory.PORT.key, 5432); params.put(JDBCDataStoreFactory.SCHEMA.key

[ArcObjects] Layer와 Feature Class/Table의 필드 별칭 설정

2011-12-25T13:31:00.000+09:00

▣ 미션 - ArcMap에서 아래 그림과 같이 테이블과 속성정보 조회시 영문필드에 한글필드로 별칭을 보여주고 싶다. - 면적 등과 같이 Numeric 필드일 경우 소숫점 2째자리까지만 보여주고 천단위 구분자를 두고싶다. - GeoDatabase(Personal, File, ArcSDE)에서 레이어를 불러올 경우 한글별칭을 기본값으로 사용하고 싶다. ▣ 설명 - Layer의 필드정보는 Feature Layer 수준, FeatureClass, Table의 필드정보는 GeoDatabase 수준에서 변경. - 따라서 FeatureClass, Table의 필드정보 변경은 물리적으로 저장되고 Layer의 필드정보는 ArcMap과 같이 Application 내에서만 임시 적용됨. -

[ArcObjects]How to identify on a map layer

2011-11-23T20:13:00.000+09:00

▣ ArcObjects Interface - IIdentify ▣ Code Snippet using ESRI.ArcGIS.Carto; using ESRI.ArcGIS.Display; using ESRI.ArcGIS.Geometry; using ESRI.ArcGIS.Geodatabase; using ESRI.ArcGIS.esriSystem; ■ Identify on a map public static void IdentifyMap(IMap focusMap, IPoint mapPoint) { if (focusMap.LayerCount == 0) return; ESRI.ArcGIS.esriSystem.UID uid = new UIDClass(); uid.Value = "{

[ArcObjects]GRID VAT(Value Attribute Table) 활용

2011-11-14T00:21:00.001+09:00

ArcGIS의 Grid 포맷을 사용할 때 Integer Grid일 경우 아래 그림과 같이 테이블을 열어 조작할 수 있습니다. 다음은 Grid의 VAT(Value Attribute Table)를 C# 코드에서 활용하는 예제입니다. - VAT에 필드 추가하기 - VAT에 셀값의 면적 계산하기 ▣ Sample Integer GRID ▣ ArcObjects Interface - IGridTableOp ▣ Code Snippet using System; using ESRI.ArcGIS.Geodatabase; using ESRI.ArcGIS.DataSourcesRaster; using ESRI.ArcGIS.GeoAnalyst; ■ VAT에 필드 추가하기 public static bool

[ArcGIS]지자체별 평균표고 미만의 셀 갯수와 면적 계산하기

2011-11-10T00:56:00.003+09:00

스크립트나 프로그램을 작성하지 않고 DEM을 이용하여 지자체별 평균표고 미만의 셀 갯수와 면적을 계산하는 과정입니다. 1. 레이어 준비 ArcMap을 실행하여 다음과 같이 시군구 및 dem 레이어를 추가한다. - 지자체 레이어 명: 시군구 - DEM 레이어 명: dem(셀사이즈 30을 가정함) 2. 필드 추가 TOC에서 시군구 레이어를 선택 후 오른쪽 마우스를 눌러 테이블을 열고 다음과 같이 필드를 추가한다. DEM_MEAN (Double) - 지자체별 평균표고 DEM_MIN (Double) - 지자체별 최소표고 DEM_MAX (Double) - 지자체별 최대표고 UND_CNT (Long Integer) - 지자체별 평균표고 미만의 셀 갯수 UND_ARA (Double) - 평균표고 미만 면적 3.

[GeoTools]GridCoverage Statistics

2011-11-06T07:25:00.002+09:00

GeoTools에서 GridCoverage(Raster)의 단순통계를 계산하는 Java 샘플 코드입니다. - JAI의 RectIter 인터페이스를 이용하여 계산 - AWT의 RenderedImage인터페이스를 이용하여 계산 - 단일 밴드에 대한 샘플입니다. ▣ Code Snippet import java.awt.image.Raster; import java.awt.image.RenderedImage; import java.text.DecimalFormat; import javax.media.jai.PlanarImage; import javax.media.jai.iterator.RectIter; import javax.media.jai.iterator.RectIterFactory; import

[ArcObjects]Raster Statistics & RasterCursor

2011-11-06T06:00:00.000+09:00

ArcObjects에서 RasterDataset의 단순통계를 계산하는 C# 샘플 코드입니다. - IRasterStatistics 인터페이스를 이용하여 계산 - IRasterCursor 인터페이스를 이용하여 계산, 이 코드에서는 Raster cursor를 이용해서 raster dataset의 pixel data에 어떻게 접근하는지도 잘 살펴 보시기 바랍니다. ▣ Code Snippet ■ Reference using System; using System.Diagnostics; using ESRI.ArcGIS.Geodatabase; using ESRI.ArcGIS.DataSourcesRaster; ■ IRasterStatistics void RasterStatistics1(IRasterDataset2

GeoTools와 GeoServer를 이용한 KOPSS Open API의 구현

2011-11-05T06:17:00.001+09:00

2011년 10월 28일 일산 킨텍스에서 열린 FOSS4G Korea 2011년 발표자료입니다. View more presentations from MinPa Lee

시군별 해안선 길이 계산해보기

2011-10-20T23:07:00.000+09:00

전국의 해안선 길이는 어떻게 계산할까? ■ 공식자료 * 국토지킴이 - 국토해양부 공식 블로그 - 해안선의 길이는 어떻게 측정할까요? * 국립해양조사원 해안선 조사측량 * 국립해양조사원 전국 해안선 연장길이 통계 자료 ■ 시군구 행정구역 경계로 시군별 해안선 길이 계산해보기 시군구 행정경계가 해안선과 일치하지는 않으나(실제보다 훨씬 짧음) 아래에서는 일반적인 방법론을 소개한다. ① 시군구 행정경계 준비하기 한번에 작업을 끝내려면 시군별로 Dissolve를 한번 해 주는 것이 좋다. ② Rectangle Polygon Shapefile 생성 ArcCatalog에서 Polygon Feature Type으로 rectangle Shapefile을 하나 생성한다. Edit... 버튼을 눌러 좌표체계는 위

Combining multiple rasters

2011-10-20T21:49:00.003+09:00

ArcGIS의 Combine (Spatial Analyst) 도구에 대한 이해 ■ 개념 - 여러 래스터의 값을 조합하여 각 셀에 중첩되는 Unique한 값을 결과 래스터에 저장하며, input raster의 값도 모두 유지된다.- 적지분석 등 래스터 중첩분석에 많이 활용된다. - 한 래스터의 각각의 래스터에 대한 정보가 포함되어 있어 데이터 조합이나 통계정보 생성이 쉽다. ■ 제약사항 - Attribute Table을 가진 Integer raster만 사용하며, 만약 Floating point raster의 경우에는 처리과정에서 자동으로 Integer로 변환한다. 따라서 Output Raster는 항상 Integer raster로 생성된다. - 한번에 처리할 수 있는 raster는 20개로 제한 -

여러개의 래스터를 하나의 래스터로 합치기

2011-10-19T19:46:00.006+09:00

■ What is a mosaic? 여러개의 Raster 파일을 하나의 Raster로 합치는 과정 [그림출처: http://help.arcgis.com/en/arcgisdesktop/10.0/help/009t/GUID-BCA5B031-B811-424B-9F54-BAB2224FBAD0-web.gif]■ 모자이크 전에 고려해야 할 사항 ① 좌표체계 on-the-fly projection을 지원하긴 하지만 미리 좌표체계를 통일하는 것이 좋다. ② 서로 다른 해상도 Mosaic 대상 Raster가 서로 다른 해상도일 경우에는 미리 Resample을 하는 것이 좋다. Continuous data type을 경우에는 Bilinear interpolation or Cubic convolution을 Discrete

격자 폴리곤을 생성하고 각 셀에 도로와의 거리를 계산하기

2011-10-17T19:07:00.001+09:00

■ 미션 격자 폴리곤을 생성하고 각 셀에 도로와의 거리를 계산하기 ■ 데이터 준비 - 도로데이터의 이름 : line ■ ArcToolbox 처리과정 1. Euclidean Distance 분석 ArcToolbox --> Spatial Analyst Tools --> Distance --> Euclidean Distance 도구를 실행한다. 셀의 크기는 500(원본 데이터 단위, 여기에서는 미터)으로 설정하였다. Direction 래스터가 필요하면 옵션으로 설정해도 된다. 또는 Spatial Analyst Toolbar --> Distance --> Straight Line 역시 Euclidean Distance 도구와 같은 분석 도구이다. 결과는 다음과 같은 도로와의 거리값이 계산된

레이어를 격자 폴리곤으로 잘라내기

2011-10-17T00:01:00.000+09:00

■ 미션 도로 데이터셋을 일정 격자 크기 단위로 잘라서 Shapefile로 저장하기 ■ ArcToolbox에서 생성하기1. 데이터셋 준비잘라낼 데이터셋(여기에서는 도로 )을 준비한다.2. Fishnet 생성 ArcToolbox --> Data Management Tools --> Feature Class --> Create Fishnet 도구를 실행한다Create Fishnet 도구의 실행 결과는 폴리곤이 아닌 Fishnet 라인(옵션으로 중심점 포인트를 생성할 수 있음)만을 생성한다. Extent는 레이어의 Extent를 참조하거나 임의로 직접 입력할 수 있다.Extent를 정의하고 나면 셀의 크기 또는 셀의 갯수를 정의한다.아래 예에서는 구경계 Extent를 기준으로 가로 4개, 세로 4개의 셀 갯수로

[ArcObjects]Define Projection...VBA

2011-10-12T15:03:00.004+09:00

ArcToolBox의 Define Projection 은 하나의 레이어를 대상으로 하기때문에.. 여러개의 레이어를 한번에 Define 할수 있는 VBA 입니다. FeatureClass를 대상으로 하지만 IGeoDataset 인터페이스를 사용함으로 코드를 조금수정하면 Raster도 같은 방법으로 할수 있습니다. Public Sub DefineProjection() Dim sPrjPath As String sPrjPath = "C:\XXXXXXX.prj" Dim pSR As ISpatialReference Dim pSRS As ISpatialReferenceFactory2 Set pSRS = New SpatialReferenceEnvironment Set pSR =

OSGeo(Open Source GeoSpatial) 한국어 지부 홈페이지

2011-10-08T15:35:00.000+09:00

OSGeo 한국어지부의 홈페이지가 개설되었습니다. 이곳에서는 OSGeo 한국어 지부와 오픈 소스 GIS 관련 정보를 제공합니다. ◎ OSGeo 한국어지부의 홈페이지 --> http://osgeo-kr.tistory.com/ ◎ OSGeo 한국어지부 메일링리스트 --> http://groups.google.com/group/osgeo-kr

스포츠정신 랜스 암스트롱

2011-09-27T21:39:00.003+09:00

[원문출처: 풋살레볼루션]