python gdal 튜토리얼: 기하학과 투영

빈 도형 객체 생성: ogr.Geometry

다양한 도형을 정의하는 데 사용되는 방법은 다릅니다(점, 선, 다각형 등)

새 점을 생성하고 다음 방법을 사용합니다. AddPoint(, , []). z 좌표는 일반적으로 생략되며 기본값은 0

예:

point = ogr.Geometry(ogr.wkbPoint)

point.AddPoint(10, 20 )

새 줄

AddPoint(, , [])를 사용하여 포인트를 추가합니다

SetPoint(< index> , , , []) 점의 좌표를 변경합니다

예를 들어 다음 코드는 점 0의 좌표를 변경합니다.

line = ogr.Geometry(ogr.wkbLineString)

line.AddPoint(10,10)

line.AddPoint(20,20)

line.SetPoint(0 ,30, 30) #(10,10) -> (30,30)

모든 포인트 수 계산

line.GetPointCount()

읽기 point 0 x 좌표와 y 좌표

print line.GetX(0)

print line.GetY(0)

새 다각형을 만들려면 먼저 새 다각형을 만듭니다. 링을 선택한 다음 다각형 개체에 링을 추가합니다.

반지는 어떻게 만드나요? 먼저 새 링 개체를 만든 다음 점을 하나씩 추가합니다.

ring = ogr.Geometry(ogr.wkbLinearRing)

ring.AddPoint(0,0)

ring.AddPoint(100,0)

ring.AddPoint(100,100)

ring.AddPoint(0,100)

마지막으로 CloseRings를 사용하여 링을 닫거나 마지막 지점의 좌표를 첫 번째 지점과 동일하게 설정합니다. 가리키다.

ring.CloseRings()

ring.AddPoint(0,0)

다음은 상자를 생성하는 예입니다. 이것은 두 개의 링 레이어로 구성된 다각형 개체입니다.

outring = ogr.Geometry(ogr.wkbLinearRing)

outring.AddPoint(0,0)

outring.AddPoint(100,0)

outring.AddPoint(100,100)

outring.AddPoint(0,100)

outring.AddPoint(0,0)

inring = ogr.Geometry(ogr.wkbLinearRing)inring = ogr.Geometry(ogr.wkbLinearRing)

inring.AddPoint(25,25)

inring.AddPoint(75,25)

inring.AddPoint(75,75)

inring.AddPoint(25,75)

inring.CloseRings()

polygon = ogr.Geometry(ogr.wkbPolygon)

polygon.AddGeometry (outring)

polygon.AddGeometry(inring)

마지막 세 문장이 더 중요합니다. 즉, 먼저 다각형 개체를 만든 다음 외부 링과 내부 링을 추가합니다

다음 문장은 다각형이 가질 수 있는 고리 수를 계산하는 데 도움이 될 수 있습니다.

polygon.GetGeometryCount() 인쇄

폴리곤에서 고리, 인덱스 순서를 읽고 다각형을 만들 때 추가합니다. 링의 순서는 동일합니다

outring =polygon.GetGeometryRef(0)

inring =polygon.GetGeometryRef(1)

복합 기하학 다중 기하학 만들기

예를 들어 MultiPoint, MultiLineString, MultiPolygon

AddGeometry를 사용하여 일반 기하학적 모양을 복합 기하학적 모양에 추가합니다. 예:

multipoint = ogr.Geometry(ogr.wkbMultiPoint)

point = ogr.Geometry(ogr.wkbPoint)point = ogr.Geometry(ogr.wkbPoint)

point.AddPoint(10,10)

multipoint.AddGeometry(point)

point.AddPoint(20,20)

multipoint.AddGeometry(point)

MultiGeometry에서 기하학을 읽는 방법은 Polygon에서 고리를 읽는 것과 동일할 수 있습니다. Polygon은 내장된 MultiGeometry라고 말했습니다.

기존 기능의 도형을 삭제하지 마십시오. Python이 충돌할 수 있습니다.

스크립트 실행 중에 생성된 도형만 삭제할 수 있습니다. 다른 함수를 호출하여 생성되었습니다. 이 형상을 사용하여 다른 기능을 생성한 경우에도 삭제할 수 있습니다.

예: Polygon.Destroy()

투영 투영의 경우 SpatialReference 개체를 사용합니다.

다양한 투영, GDAL은 WKT, PROJ.4, ESPG, USGS를 지원합니다. , ESRI.prj

는 레이어 및 형상에서 투영을 읽을 수 있습니다. 예:

spatialRef = layer.GetSpatialRef()

spatialRef = geom.GetSpatialReference()

프로젝션 정보는 일반적으로 .prj 파일에 저장됩니다. 해당 파일이 없으면 위 함수는 None을 반환합니다.

새 프로젝션 만들기:

먼저 osr 라이브러리를 가져오고 그런 다음 osr.SpatialReference()를 사용합니다. SpatialReference 개체

를 만든 후 다음 명령문을 사용하여 투영 정보를 SpatialReference 개체

로 가져옵니다. •ImportFromWkt()

•ImportFromEPSG()

•ImportFromProj4()

•ImportFromESRI()

•ImportFromPCI(,

  )

 •ImportFromUSGS(, )

•ImportFromXML( )

프로젝션 내보내기, 사용 다음 명령문을 문자열로 내보낼 수 있습니다

•ExportToWkt()

•ExportToPrettyWkt()

•ExportToProj4()

•ExportToPCI()

•ExportToUSGS()

•ExportToXML()

기하학적 모양 Geometry에 대한 투영 변환을 수행하려면 먼저 두 개의 Projections를 초기화한 다음 CoordinateTransformation 개체를 생성하고 이를 변환에 사용해야 합니다.

sourceSR = osr.SpatialReference()

sourceSR .ImportFromEPSG(32612) #UTM 12N WGS84

targetSR = osr.SpatialReference()

targetSR.ImportFromEPSG(4326) #Geo WGS84

coordTrans = osr.CoordinateTransformation(sourceSR , targetSR)

geom.Transform(coordTrans)

하지만 이 코드는 매우 고통스럽습니다! Windows에서는 작동하지 않습니다. 외국인 포럼에서는 리눅스는 문제가 없는데 윈도우는 죽었다는 논의가 있습니다. . .

주의해야 할 몇 가지 사항이 있습니다.

적절한 시점에 형상을 편집하는 것이 투영 변환 후에는 이동하지 않는 것이 가장 좋습니다.

DataSource의 모든 Geometry에 대해 투영 변환을 수행하려면 하나씩 수행해야 합니다. 루프를 사용하십시오.

프로젝션을 .prj 파일에 작성하는 것은 실제로 매우 간단합니다. 먼저 MorphToESRI()를 문자열로 변환한 다음 텍스트 파일을 열고 그 안에 씁니다. 예:

targetSR.MorphToESRI()

file = open('test.prj', 'w')

file.write(targetSR.ExportToWkt())

ffile.close()

위 내용은 Python gdal 튜토리얼의 내용입니다: 기하학과 투영에 대한 자세한 내용은 PHP 중국어 웹사이트(www.php. CN)!