Нужна помощь в создании prj-файла

[witec]

Добрый день всем!

Есть shp-файл который был создан в ArcView(и заданы определенные параметры в виде Projection: Transverse Mercator; Spheroid: Krasovsky; Central Meridian: 27; Reference Latitude:0; Scale Factor: 1; False Easting: 5500000; False Northing: 0).
Я пробовал создавать prj в Arcview Projection Utility и ArcCatalog. В обоих случаях он создавался, но координаты не совпадают с теми которые показываются в ArcView(при установлении настроек в ввиде) и шеп-файл искажается в ArcMap, возможно я не правильно установил проекцию или систему координат?!
Помогите разобраться

[Максим Дубинин]

как именно создавались данные?

[witec]

как именно создавались данные?

Устанавливались параметры проекции в виде ArcView(Projection: Transverse Mercator; Spheroid: Krasovsky; Central Meridian: 27; Reference Latitude:0; Scale Factor: 1; False Easting: 5500000; False Northing: 0), затем добавлялись: привязанный растр(Gauss Kruger-Pulkovo1942 GK Zone 5) и создавался shp-файл(в ArcView, полигональный) и производилась вектаризация.

[Максим Дубинин]

в этом случае у вас векторные данные в Geographic/Pulkovo1942

загрузите слой в чистый проект и повозите мышкой, если градусы, то так и есть и создавать prj нужно соответствующий
http://gis-lab.info/qa/proj-sk-faq.html#4

[ggis]

может так:

[witec]

Да, так и есть при добавлении shp-файла в чистый проект(ArcMap Рис1) показываются градусы и они соответствуют реальным данным. Но вот в чем проблема, полигоны искажены(и при установки Geographic/Pulkovo1942). И координаты другие, в ArcView метры(если в чистом проекте открыть то показывает 24 градуса и 52(рис2), а если в ArcView и с параметрамиProjection: Transverse Mercator; Spheroid: Krasovsky; Central Meridian: 27; Reference Latitude:0; Scale Factor: 1; False Easting: 5500000; False Northing: 0, то показывает 5345... и 58201....
И похоже здесь проблема в растровой основе по которому создавался вектор(я точно не уверен в какой проекции он(растр) привязывался т. к. координаты он показывает 5345... и 58201....(рис 2), а если поверить в Google.maps то должны быть 24 градуса и 52 ).

Помогите плиз правильно создать prj, для того чтобы в дальнейшем можно было конвертнуть.
Прикрепил: часть shp и привязанную растровую основу.

[witec]

Растр

[ggis]

Export_Output.shp-в географической СК Pulkovo42, для него ее надо прописать в арккаталоге, потом начитать в аркмап и установить проекцию TM и они будут отображаться вместе с растровой картой (для нее также предварительно в арккаталоге прописать проекцию ТМ, ту что описали выше) нормально, можно сохранить Export_Output.shp в новый шейп уже в метрах (спроецировать), прикрепил архив.

[Александр Мурый]

IMHO, гораздо проще пойти через GDAL + базу проекций EPSG. Зачем вручную вбивать параметры, если за нас уже поработали?

-- Зададим растру "потерянную" проекцию Pulkovo 1942 / Gauss-Kruger zone 5 :
gdalwarp -s_srs "EPSG:28405" Береза.bmp Bereza.tif
Получается:

Код: Выделить всё

Driver: GTiff/GeoTIFF
Files: Bereza.tif
Size is 3510, 3814
Coordinate System is:
PROJCS["Pulkovo 1942 / Gauss-Kruger zone 5",
    GEOGCS["Pulkovo 1942",
        DATUM["Pulkovo_1942",
            SPHEROID["Krassowsky 1940",6378245,298.2999999999998,
                AUTHORITY["EPSG","7024"]],
            AUTHORITY["EPSG","6284"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4284"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",27],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",5500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","28405"]]
Origin = (5329724.974118598736823,5840590.890793216414750)
Pixel Size = (10.000000000000000,-10.000000000000000)
Corner Coordinates:
Upper Left  ( 5329724.974, 5840590.891) ( 24d28'58.63"E, 52d39'59.60"N)
Lower Left  ( 5329724.974, 5802450.891) ( 24d30' 8.79"E, 52d19'26.91"N)
Upper Right ( 5364824.974, 5840590.891) ( 25d 0' 5.11"E, 52d40'35.17"N)
Lower Right ( 5364824.974, 5802450.891) ( 25d 1' 0.85"E, 52d20' 2.05"N)
Center      ( 5347274.974, 5821520.891) ( 24d45' 3.37"E, 52d30' 1.96"N)
Band 1 Block=3510x2 Type=Byte, ColorInterp=Palette
  Color Table (RGB with 256 entries)

-- Для перепроецирования шейпа в ГК зону 5:
ogr2ogr -s_srs "EPSG:4326" -t_srs "EPSG:28405" Export_Output_GKzone5.shp Export_Output.shp

Самое интересное, что при наложении вектора на растр непонятно, что оцифровывалось (или это только у меня?)

[witec]

Все получилось! Всем БОЛЬШОЕ СПАСИБО за помощь('Максим Дубинин', 'ggis', 'amuriy')!!!

Самое интересное, что при наложении вектора на растр непонятно, что оцифровывалось (или это только у меня?)

Да, так и есть т.к. было несколько растров, но прикрепленный основной, к которому все и привязывалось.

А как обратно преобразовать shp, т.е. из преобразованного(TM) в тот, который был изначально( в географической СК Pulkovo42)?

[ggis]

2 способа:
1. Задать в аркмап для View проэкцию географическую СК Pulkovo42 и сохранить в новый шейп, задав при сохранении "сохранить в единицах активного View"
2. С помощью арктулбокса Data Management tools-Projection and Transformation-Features-Project (указать шейп файл (TM) и указать куда сохранить новый шейп-файл и задать географическую СК Pulkovo42 )

[Александр Мурый]

В географ. координаты можно ещё и так:

Код: Выделить всё

ogr2ogr -t_srs "EPSG:4326" Export_Output_wgs84.shp Export_Output_GKzone5.shp

[witec]

Спасибо за помощь!!!