search

nextgis-borsch / lib_gdal

GDAL is a translator library for raster and vector geospatial data formats that is released under an X/MIT style Open Source license by the Open Source Geospatial Foundation. Reorganized GDAL repository
http://gdal.org
28 stars 9 forks source link

wrapdateline reprojection problem #22

Open simgislab opened 6 years ago

simgislab commented 6 years ago 
ogr2ogr -t_srs "EPSG:4326" -datelineoffset 5 -wrapdateline output.shp input.shp

При перепроецировании две проблемы:

  1. пропадает один или несколько полигонов
  2. один полигон не режется как нужно

Пробовали разные datelineoffset

В процессе ругается на ошибку топологии:

ERROR 1: TopologyException: Input geom 0 is invalid: Self-intersection at or near point 180 66.690594954016618 at 180 66.690594954016618
ERROR 1: TopologyException: Input geom 0 is invalid: Self-intersection at or near point 180 66.690594954016618 at 180 66.690594954016618

Проверено в:

ogr2ogr --version
GDAL 2.2.1, released 2017/06/23
GEOS: 3.6.1-CAPI-1.10.1

ogr2ogr --version
GDAL 2.2.2, released 2017/09/15
GEOS: 3.5.0-CAPI-1.9.0

Исходная проекция данных:

PROJCS["Asia_North_Equidistant_Conic_Krass",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Equidistant_Conic"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",95.0],PARAMETER["Standard_Parallel_1",49.0],PARAMETER["Standard_Parallel_2",68.5],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]]

Исходные данные: input.zip. Можно так же экспериментировать с одним проблемным полигоном: bad-poly.zip Проблем с топологией предварительной проверкой выявлено не было. Успешное перепроецирование с разрезкой конкурирующим продуктом): output.zip

simgislab commented 6 years ago

добавил шейп с одним полигоном который неправильно режется