Daha çox

Gdal / Ogr'da Koordinat Dəyişikliyi (EPSG: 4326 - EPSG: 3068)


Məncə gdal / ogr (/ proj4) epsg tanımır: 3068 düzdür?

Budur cəhd etdiyim şey:

public static void transformSRS () {// EPSG yaradılması: 4326 SpatialReference sourceSRS = yeni SpatialReference (); sourceSRS.ImportFromEPSG (4326); // EPSG yaratmaq: 3068 SpatialReference targetSRS = yeni SpatialReference (); sourceSRS.ImportFromEPSG (3068); // həm SRS-in System.out.println (sourceSRS) çap; // System.out.println (targetSRS) yazdırılır; // çap olunmur !? // EPSG-də bir nöqtə yaratmaq: 3068 koordinat 850, -150 Həndəsə nöqtəsi = Həndəsə.CreateFromWkt ("POINT (850 -150)"); // sistemin EPSG-də olduğunu bilməsi üçün nöqtəyə EPSG: 3068 təyin etmək: 3068 nöqtə.AssignSpatialReference (sourceSRS); // onu EPSG-yə çevirməyə çalışır: 4326, lakin istisna meydana gəlir. point.TransformTo (targetSRS); }

və bu çıxışdır (ilk System.out.println () işləyir, ancaq ikincisi sıfır və istisna olmaqla)point.TransformTo (targetSRS)):

GEOGCS ["WGS 84", DATUM ["WGS_1984", SPHEROID ["WGS 84", 6378137,298.257223563, Authority ["EPSG", "7030"]], Authority ["EPSG", "6326"]], PRIMEM [ "Greenwich", 0, YETKİ ["EPSG", "8901"]], UNIT ["dərəcə", 0.0174532925199433, YETKİ [[EPSG "," 9122 "]], YAZI [[EPSG", "4326"]] XATA 6: Boş bir SRS-dən PROJ.4 formatına heç bir tərcümə məlum deyil. "Ana" java.lang.RuntimeException işində istisna: OGR Xətası: org.gdal.ogr.ogrJNI.Geometry_TransformTo (Native Method) org.gdal.ogr.Geometry.TransformTo (Geometry.java:409) at org .geotools.mavenGdal.Test.transformSRS (Test.java:880) at org.geotools.mavenGdal.Test.main (Test.java:109)

Bunu gdalsrsinfo ilə test edə bilərsiniz http://www.gdal.org/gdalsrsinfo.html

gdalsrsinfo epsg: 3068 PROJ.4: '+ proj = cass + lat_0 = 52.41864827777778 + lon_0 = 13.62720366666667 + x_0 = 40000 + y_0 = 10000 + ellps = bessel + towgs84 = 598.1,73.7,418.2,0.202,045, -2. 6.7 + vahid = m + no_defs 'OGC WKT: PROJCS ["DHDN / Soldner Berlin", GEOGCS ["DHDN", DATUM ["Deutsches_Hauptdreiecksnetz", SPHEROID ["Bessel 1841", 6377397.155,299.152812, 7004 "]], TOWGS84 [598.1,73.7,418.2,0.202,0.045, -2.455,6.7], YETKİ [[EPSG", "6314"]], PRIMEM ["Greenwich", 0, OTORITY ["EPSG", " 8901 "]], UNIT [" derece ", 0.0174532925199433, YETKİ [[EPSG", "9122"]], YETKİ [[EPSG "," 4314 "]], PROJEKSİYA [" Cassini_Soldner "], PARAMETER [" enlem_of_origin ", 52.41864827777778], PARAMETER ["central_meridian", 13.62720366666667], PARAMETER ["false_easting", 40000], PARAMETER ["false_northing", 10000], UNIT ["meter", 1, Authority ["EPSG", "9001) Səlahiyyət ["EPSG", "3068"]]

Digər bir test, gdaltransform http://www.gdal.org/gdaltransform.html ilə koordinatları çevirməkdir

gdaltransform -s_srs epsg: 3068 -t_srs epsg: 4326 850 -150 13.0512110820125 52.3246380119762 41.6757647842169

Nəticə: Heç olmasa mənim GDAL EPSG-ni tanıyır: 3068.


Giriş

Koordinat istinad sistemlərinin (CRS) təmsilçiliyində və koordinatlar üzərində əməliyyatlarda on illərdir davam edən dəyişikliklər, məkan obyektlərinin R paketlərində işlənməsi qaydasında həyata keçirilməlidir. 1990-cı illərə qədər, əksər məkan məlumatları, sadəcə yerli xəritələşdirmə orqanı tərəfindən verilən koordinatları istifadə edirdi, Meuse bank məlumat dəsti EPSG: 28992, “Amersfoort / RD New” olduğu ortaya çıxan sayğaclarda planar bir təsvir istifadə etdi. CRS-in niyə Clifford J. Mugnier-in sütunları olduqları kimi göstərildiyini öyrənmək üçün böyük bir qaynaqdır Fotogrammetrik Mühəndislik & amp Uzaqdan Algılama, istinadlar mövcuddur rgdal Hollandiya, Fevral 2003 sütununda yer aldı:

Əksər milli xəritələşdirmə agentlikləri öz standart coğrafi və proqnozlaşdırılan CRS-lərini müəyyənləşdirsələr də, hərbi ittifaqlar və müstəmləkə idarələri kimi millətlərarası qurumlar tez-tez milli sərhədlər üzərindəki əməliyyatları asanlaşdırmaq üçün bir sıra qanunauyğunluqlar tətbiq etdilər. Bu, eyni zamanda Avropa Petrol Tədqiqat Qrupunun (EPSG) yaradılmasına gətirib çıxardı, çünki dəniz yurisdiksiyası nizamlı deyildi və eyni sahil şelfini paylaşan ölkələrin ziddiyyətli kəşfiyyat güzəştləri təyin etməyə çalışdıqları zaman əhəmiyyətli idi. Neft şirkətlərindən olan mütəxəssislər, Beynəlxalq Standartlar Təşkilatı (ISO, xüsusilə TC 211) və Açıq Yerleşim Konsorsiumu (OGC) ilə təmin olunan böyük təcrübə toplamışlar.

Digər məlumatları fərqli bir CRS ilə inteqrasiya edərkən və veb xəritə fonunda göstərmək üçün CRS-nin müəyyənləşdirilməsi zərurətə çevrildi. ESRI Shapefile formatı kimi bir çox köhnə fayl formatı, mövqe məlumatlarının CRS-in daxil edilməsini məcbur etməmişdir. Daha sonra açıq mənbəli proqramların əksəriyyəti PROJ.4 simlərindən çevik bir nümayəndəlik kimi istifadə edirdi, lakin beynəlxalq səviyyədə qəbul edilmiş standartlara, xüsusən də ISO 19111-ə çatdıqdan və zamanla təkrarlanaraq yaxşılaşdırıldığından, müasir olaraq tanınan müasir bir mətn təmsilçiliyinə keçmək lazımdır. WKT2 (2019). İndi sanki demək olar ki, bütün şirkətlər və xəritələşdirmə agentlikləri bu nümayəndəliyi qəbul edir və qəbul edilmişdir sp vasitəsilə rgdal, sf və digər paketlər.


Bu RFC-nin əsas problemi, EPSG koordinat sistemlərindən geniş istifadə edən, lakin ox istiqamətlərini ləğv edən və bunlara baxılması lazım olduğunu düşünən mövcud faylları və kodu pozmadan, coğrafi koordinat sistemləri üçün EPSG istədiyi ox sifarişini yerinə yetirmək də daxil olmaqla, balta sifarişlərinə dəstək əlavə etməkdir. EPSG-nin dediklərindən asılı olmayaraq uzun, lat.

Xüsusilə, EPSG: 4326 kimi bir coğrafi koordinat sistemindəki bir sənədin nə zaman lat / uzun və nə vaxt uzun / lat kimi qəbul ediləcəyinə qərar vermək üçün uyğun siyasət və mexanizmlər təklif edirik. Mövcud geniş təcrübə sayəsində keçmiş təcrübə ilə yanılmağımızı və EPSG ox sifarişinə hörmətlə yanaşmağı "seçməyimizi" tələb edirik.


1564 saylı coğrafi 3D ofset istifadə edərək transformasiyaların aparılması

GDAL / OGR / Proj vasitəsi ilə bir tətbiqdə iş vaxtı SRS-lərin bir alt hissəsindən WGS84-ə şaquli bir məlumat çevrilməsini həyata keçirməyə çalışıram. Bunu etmək üçün, WGS84 (epsg: 4326) -ni üfüqi sənəd kimi istifadə vəziyyətim üçün EGM96 (epsg: 5773) olduğumuzu düşünə biləcəyimiz şaquli verilənlər bazası kimi istifadə etdiyim bir kompleks CS inşa edirəm. WGS84 (epsg: 4979) ilə uyğun bir 3D koordinat sisteminə layihə hazırlamaq istəyirəm.

Transformasiyamı belə qururam:

Mənbə istinad çərçivəsindən təyinata qədər yüksəklik məlumatlarını proyeksiya etmək üçün bu çevrilmədən istifadə etməyə çalışıram:

Dönüşüm yaradıldığında bir istisna alıram:

koordinateoperasiya.cpp. Oradakı kod yollarının ətrafında gəzdikdən sonra metodum EPSGCode-un 9660 olduğunu və məni əvvəlcədən işləyən bayraqlarla çıxarılan bir kodepatdan salacağımı bir yerə yığa bildim:

Orada niyə qismən tətbiq olunmasının xüsusi bir səbəbi varmı? Bu çevrilmə növü Proj tərəfindən dəstəklənirmi? Bu çevrilməni həyata keçirmək üçün daha asan bir yol varmı?

Mətn uğurla yeniləndi, lakin bu səhvlərlə qarşılaşdı:

Hazırda tapşırığı bir problemə çevirə bilmirik. Zəhmət olmasa bir daha cəhd edin.

Məsələ uğurla yaradıldı, lakin şərhi hazırda yeniləyə bilmirik.


Koordinat Sistem obyektləri¶

Məkan İstinad ¶

Məkan istinad obyektləri aşağıdakılardan biri ola bilən srs_input-da başlanğıc edilir:

  • OGC yaxşı bilinən mətn (WKT) (bir simli)
  • EPSG kodu (tam və ya simli)
  • PROJ.4 simli
  • Tanınmış standartlar üçün stenoqrafiya sətri ('WGS84', 'WGS72', 'NAD27', 'NAD83')

Verilən simli atribut qovşağının dəyərini qaytarır, Düyün yoxdursa, Yoxdur. Uşağın WKT-dəki atributun indeksinin olduğu bir parametr olaraq (hədəf, uşaq) bir toplu da götürə bilər. Misal üçün:

Verilən hədəf nodu üçün atribut dəyəri (məs. 'PROJCS'). İndeks açar sözü qayıdacaq uşaq qovşağının indeksini göstərir.

Verilmiş sətir hədəf düyünü üçün səlahiyyət adını qaytarır.

Verilmiş simli hədəf düyününün səlahiyyət kodunu qaytarır.

Bu məkan istinad obyektinin klonunu qaytarır.

Bu metod bu SpatialReference-ın WKT-sini yoxlayır və EPSG identifikatorunun tətbiq olunduğu yerdə EPSG səlahiyyət düyünlərini əlavə edəcəkdir.

Bu Məkan Referansını ESRI & # 8217s formatından EPSG-yə dəyişdirir

Bu Məkan Referansını ESRI & # 8217s formatına çevirir.

Verilən məkan istinadının etibarlı olub olmadığını yoxlayır, istisna deyilsə.

EPSG kodundan məkan istinadını idxal edin.

PROJ.4 sətirindən məkan istinadını idxal edin.

import_user_input (user_input) [mənbə] ¶ import_wkt (wkt) [mənbə] ¶

WKT-dən məkan istinadını idxal edin.

XML-dən məkan istinadını idxal edin.

Bu Məkan İstinadının adını qaytarır.

Üst səviyyəli səlahiyyətin SRID-i qaytarır və ya təyin edilmədiyi təqdirdə Yoxdur.

Xətti vahidlərin adını qaytarır.

Xətti vahidlərin dəyərini qaytarır.

Açısal vahidlərin adını qaytarır. & # 8221

Açısal vahidlərin dəyərini qaytarır.

Vahidlər dəyərinin və vahidlərin adının 2 cütlüyünü qaytarır və xətti və ya açısal vahidlərin qaytarılmasını avtomatik olaraq təyin edəcəkdir.

Bu məkan referansı üçün elipsoid parametrlərinin toplağını qaytarır: (yarı böyük ox, yarı yarım ox və tərs düzləşmə)

Bu məkan referansı üçün elipsoidin yarı əsas oxunu qaytarır.

Bu məkan referansı üçün elipsoidin yarı kiçik oxunu qaytarır.

Bu məkan referansı üçün elipsoidin tərs düzlənməsini qaytarır.

Bu məkan referansı coğrafi olduqda True (Kök qovşağı GEOGCS) qaytarır.

Bu məkan referansı lokaldırsa (kök nodu LOCAL_CS) True qaytarır.

Bu məkan referansı proqnozlaşdırılan bir koordinat sistemi olduqda (kök nodu PROJCS) True qaytarır.

Bu məkan istinadının WKT təqdimatını qaytarır.

WKT-nin & # 8216pretty & # 8217 nümayəndəliyini qaytarır.

Bu məkan referansı üçün PROJ.4 təqdimatını qaytarır.

Bu məkan istinadının XML təqdimatını qaytarır.

CoordTransform ¶

Koordinat sistemi çevrilməsini təmsil edir. Müvafiq olaraq mənbə və hədəf koordinat sistemlərini təmsil edən iki SpatialReference ilə başlanğıc edilir. Bu obyektlər fərqli həndəsələrdə eyni koordinat çevrilməsini dəfələrlə həyata keçirərkən istifadə olunmalıdır:


Sponsorlar

Səy 20 qurumdan səxavətli sponsorluq aldı və 15 MAY 2018-ci il tarixinə qədər 144.000 ABŞ dolları toplamaq məqsədini 125.000 ABŞ dollarından keçdi.

Təşkilat Məbləğ
anonim $5,000
AGI $2,000
AppGeo $5,000
Tətbiqi Təsəvvür $5,000
Azavea $5,000
Hüdudsuz $5,000
Cizgi filmi $5,000
Dynamic Graphics, Inc $5,000
ESRI $30,000
GeoCue $5,000
İntevasiya $1,000
LizardTech $5,000
LINZ $7,000
Xəritə qutusu $5,000
Mobil Coğrafiya $500
Oslandiya $2,000
Təhlükəsiz proqram $31,000
Septima $1,500
Sparkgeo $5,000
Məkan Şəbəkələri $1,000
RIEGL $7,000
Radiant Solutions $5,000
Terranodo $500
Timoney Qrupu $500
Ümumi $144,000


Gdal / Ogr'da Koordinat Dəyişikliyi (EPSG: 4326 - EPSG: 3068) - Coğrafi İnformasiya Sistemləri

Mənbə və hədəf koordinat sistemlərini təyin etmək üçün girişlər üçün yalnız EPSG kodlarını alan bir skript üçün Python üçün GDAL API istifadə edirəm.

Cs2cs əmri Şimali Amerikadakı NAD27 (EPSG: 4267) və NAD83 (EPSG: 4326) koordinatları arasında olduqca yaxşı bir çevrilməni idarə etsə də, AGD66 (EPSG: 4202) üçün GDA94 (EPSG: 4283) ilə oxşar imkanlara sahib olmaq istərdim. və AGD84 (EPSG: 4203) ilə Avstraliyada GDA94 (EPSG: 4283) arasında.

-V seçimi istifadə olunan proyeksiya sətirinə baxmaq üçün cs2cs ilə istifadə edildikdə aşağıdakıları əldə edirəm:

Nadgrids parametrləri, cs2cs-ə konus, Alyaska, ntv2_0.gsb və ntv1_can.dat-dan ən yaxşı uyğunlaşdırma datur grid shift faylını istifadə etməsini söyləyir.

AGD66 (EPSG: 4202) üçün cs2cs qaytarır:

Nadgrids parametri yox, bunun yerinə bir towgs84 yoxdur. "A66 National (13.09.01) .gsb" "ICSM-dən bir grid shift faylı, AGD66 və GDA94 arasındakı dəyişikliyi ehtiva edir.

EPSG kodlarını istifadə edərkən GDAL / cs2cs / proj-in bu faylı necə müəyyənləşdirə biləcəyəm?

GDAL / cs2cs / projesi osgeo4w installer-dən istifadə edərək A66 National (13.09.01) .gsb ilə yenidən quraşdırmadan əvvəl digər tarix dəyişikliyi sənədlərini ehtiva edən qovluqda yenidən qurmağa çalışdım, ancaq müvəffəq olmadı.


OGR Həndəsələri¶

OGRGeometriya ¶

OGRGeometry obyektləri GEOSGeometry obyektləri ilə oxşar funksiyanı bölüşür və OGR & # 8217s daxili həndəsə nümayişi ətrafında incə sarmaldır. Beləliklə, DataSource istifadə edərkən məlumatlara daha səmərəli giriş imkanı verirlər. GEOS həmkarından fərqli olaraq, OGRGeometry məkan istinad sistemlərini və koordinat çevrilməsini dəstəkləyir:

Bu obyekt OGR Həndəsə sinfi üçün bir sarğıdır. Bu obyektlər birbaşa WKT, HEX, GeoJSON, WKB məlumatlarını ehtiva edən bir tampon və ya OGRGeomType obyekti ola bilən bir geom_input parametrindən təsirlənir. Bu obyektlər Layer-dən vektor məlumatlarını oxuduqda (bu da bir DataSource-un bir hissəsidir) Feature.geom atributundan qaytarılır.

sinif metodu qutudan_boks ( bbox ) ¶

Verilmiş məhdudlaşdırma qutusundan (4 köşeli) bir çoxbucaqlı qurur.

LineStringdəki nöqtələrin sayını, Çoxbucaqdakı üzüklərin sayını və ya GeometryCollectiondakı həndəsələrin sayını qaytarır. Digər həndəsə növləri üçün tətbiq olunmur.

LineStringdəki nöqtələr, çoxbucaqlı halqalar və ya GeometryCollectiondakı həndəsələr üzərində təkrarlanır. Digər həndəsə növləri üçün tətbiq olunmur.

LineString üçün göstərilən indeksdəki nöqtəni, bir çoxbucaqlı üçün göstərilən indeksdəki daxili halqanı və ya GeometryCollection-da göstərilən indeksdəki həndəsəni qaytarır. Digər həndəsə növləri üçün tətbiq olunmur.

Həndəsənin koordinasiyalı ölçülərinin sayını qaytarır, yəni nöqtələr üçün 0, xətlər üçün 1 və s.

Bu həndəsənin koordinat ölçüsünü qaytarır və ya təyin edir. Məsələn, iki ölçülü həndəsələr üçün dəyər 2 olardı.

Bu həndəsədəki elementlərin sayını qaytarır:

Bu həndəsi təsvir etmək üçün istifadə olunan nöqtələrin sayını qaytarır:

OGRGeomType obyekti olaraq bu həndəsənin növünü qaytarır.

Bu həndəsənin növünün adını qaytarır:

Bu həndəsənin sahəsini və ya sahəsi olmayan həndəsələr üçün 0 qaytarır:

Bu həndəsənin zərfini zərf obyekti olaraq qaytarır.

Bu həndəsənin zərfini zərf obyekti əvəzinə 4 köşeli olaraq qaytarır:

Bu xüsusiyyət, bu həndəsə üçün məkan istinadını və ya ona heç bir məkan istinad sistemi təyin edilmədiyi təqdirdə Yoxdur. Əgər təyin olunarsa, bu xassəyə daxil olmaq bir SpatialReference obyektini qaytarır. Başqa bir SpatialReference obyekti və ya SpatialReference-in qəbul etdiyi hər hansı bir giriş ilə təyin edilə bilər. Misal:

Bu həndəsənin SpatialReference-a uyğun olan məkan istinad identifikatorunu qaytarır və ya təyin edir. Bu həndəsə ilə əlaqəli bir məkan istinad məlumatı olmadıqda və ya bir SRID təyin edilə bilmədikdə Yoxdur.

Bu həndəsəyə uyğun bir GEOSGeometry obyektini qaytarır.

Bu həndəsənin GML formatında bir simli təsvirini qaytarır:

Bu geometriyanın HEX WKB formatında bir simli təsvirini qaytarır:

JSON formatında bu həndəsənin bir simli təsvirini qaytarır:

Bu həndəsənin KML formatında bir simli təsvirini qaytarır.

Bu həndəsənin WKB təsvirini saxlamaq üçün lazım olan WKB tamponunun ölçüsünü qaytarır:

Bu həndəsənin WKB təsvirini ehtiva edən bir tamponu qaytarır.

WKT formatında bu həndəsənin simli təsvirini qaytarır.

Bu həndəsənin EWKT təsvirini qaytarır.

Bu həndəsə obyektinin yeni bir OGRGeometry klonunu qaytarır.

Bu həndəsə içərisində bağlanmamış hər hansı bir üzük varsa, bu rutin sonu başlanğıc nöqtəsini əlavə edərək edəcəkdir:

Bu həndəsəni fərqli bir məkan istinad sisteminə çevirir. Bir CoordTransform obyekti, bir SpatialReference obyekti və ya SpatialReference tərəfindən qəbul edilmiş hər hansı digər giriş (məkan referansı WKT və PROJ.4 simləri və ya tam bir SRID daxil olmaqla) ala bilər. Varsayılan olaraq heç bir şey qaytarılmır və həndəsə yerində çevrilir. Lakin, klon açar sözü True olaraq təyin olunarsa, bunun əvəzinə bu həndəsənin çevrilmiş klonu qaytarılır.

Bu həndəsə digəri ilə kəsişirsə True qaytarır, əks halda False qaytarır.

Bu həndəsə digərinə bərabərdirsə, əks halda False qaytararsa True qaytarır.

Bu həndəsə digərinə (yəni kəsişməsə) məkan baxımından ayrılıbsa, əks halda False qaytararsa True qaytarır.

Bu həndəsə digərinə toxunduqda True qaytarır, əks halda False qaytarır.

Bu həndəsə digərini keçərsə, əks halda False qaytararsa True qaytarır.

Bu həndəsə digərinin içərisindədirsə True qaytarır, əks halda False qaytarır.

Bu həndəsə digərini ehtiva edərsə True qaytarır, əks halda False qaytarır.

Bu həndəsə digərinin üst-üstə düşdüyü təqdirdə True qaytarır, əks halda False qaytarır.

Bu həndəsənin sərhədi, yeni bir OGRGeometry obyekti olaraq.

Bu həndəsəni ehtiva edən ən kiçik qabarıq çoxbucaqlı, yeni bir OGRGeometry obyekti olaraq.

Bu həndəsə ilə digərinin fərqindən ibarət olan bölgəni yeni OGRGeometry obyekti olaraq qaytarır.

Bu həndəsə ilə digərinin kəsişməsindən ibarət olan bölgəni yeni OGRGeometry obyekti olaraq qaytarır.

Bu həndəsənin və digərinin simmetrik fərqindən ibarət olan bölgəni yeni OGRGeometry obyekti olaraq qaytarır.

Bu həndəsənin və digərinin birləşməsindən ibarət olan bölgəni yeni OGRGeometry obyekti olaraq qaytarır.

Nöqtə həndəsəsinin koordinatlarını bir xəndək, bir xətt həndəsəsinin koordinatlarını bir güzgü kimi və s. Qaytarır:

Bu nöqtənin X koordinatını qaytarır:

Bu nöqtənin Y koordinatını qaytarır:

Bu nöqtənin Z koordinatını qaytarır və ya nöqtənin Z koordinatına sahib olmadığı təqdirdə Yoxdur:

Bu sətirdə X koordinatlarının siyahısını qaytarır:

Bu sətirdə Y koordinatlarının siyahısını qaytarır:

Bu sətirdə Z koordinatlarının siyahısını və ya sətirdə Z koordinatları yoxdursa Yoxdur:

LinearRing həndəsi olaraq bu çoxbucağın qabığını və ya xarici halqasını qaytarır.

Bu çoxbucağın sentroidini təmsil edən bir nöqtəni qaytarır.

sinif GeometryCollection ¶ əlavə et ( geom ) ¶

Bu həndəsə kolleksiyasına həndəsə əlavə edir. Digər həndəsə növləri üçün tətbiq olunmur.

OGRGeomType ¶

Bu sinif bir OGR həndəsə növünün bir neçə yolla təmsil olunmasına imkan verir:

OGR Həndəsə tipinin qısa bir simli formasını qaytarır:

OGR həndəsə növünə uyğun ədədi qaytarır:

Bu OGR tipini saxlamaq üçün istifadə etmək üçün Django sahə növünü (GeometryField alt sinfi) və ya uyğun Django növü olmadıqda Yoxdur qaytarır:

Zərf ¶

Bir düzbucaqlı məhdudlaşdıran qutu üçün minimum və maksimum X, Y koordinatlarını ehtiva edən bir OGR zərf quruluşunu təmsil edir. Dəyişənlərin adlanması OGR Zərf C quruluşu ilə uyğundur.

Minimum X koordinatının dəyəri.

Maksimum X koordinatının dəyəri.

Minimum Y koordinatının dəyəri.

Maksimum Y koordinatının dəyəri.

Yuxarı sağ koordinat, bir topl kimi.

Sol alt alt koordinat.

Zərfi təmsil edən küp.

Bu zərfi WKT formatında çoxbucaqlı olaraq təmsil edən bir simli.

genişləndirmək_ daxil etmək ( özünü, * args ) ¶


Raster və vektor formatında dəstək¶

Yazı zamanı GDAL / OGR tərəfindən idarə olunan heç bir format koordinat dövrünü standartlaşdırma üsuluna malik deyildir. Nəticədə mövcud oxucularla mümkün qədər geriyə uyğun olmaq məqsədi ilə onu necə kodlayacağımızı seçdik. Bu konsepsiya nəzərə alınmaq üçün müvafiq rəsmi spesifikasiyalar inkişaf edərsə, bu kodlamalar dəyişə bilər. Koordinat dövrü yalnız yaradılan təbəqənin / verilənlər bazasının SRS-yə qoşulduqda yazılır.

FlatGeoBuf¶

Koordinat dövrü, FlatGeoBuf faylının Crs.wkt başlıq sahəsindəki COORDINATEMETADATA quruluşunun EPOCH alt qovluğundan istifadə edərək WKT: 2019 sətri kimi kodlanır.

Bu cür quruluş GDAL & lt 3.4 tərəfindən başa düşülməyəcək, lakin CRS ilə əlaqəli bir EPSG kodu varsa, bu köhnə GDAL versiyalarında problem yaratmayacaqdır.

GeoPackage vektoru / raster¶

Əlaqələndirilmiş bir koordinat dövrü olan hər bir vektor / raster cədvəli, CRS WKT uzantısının genişləndirilmiş bir versiyasından (https://github.com/opengeospatial/geopackage/pull/600) istifadə edərək gpkg_spatial_ref_sys cədvəlinin dövr sütununda kodlaşdırır.

GeoTIFF¶

Koordinat dövrü yeni bir GeoTIFF GeoKey, 5120 kodunun CoordinateEpochGeoKey kimi kodlanır və DOUBLE yazın.

JPEG2000¶

GeoJP2 qutuları yuxarıda göstərilən GeoTIFF kodlamasından istifadə edir.

Davamlı Köməkçi Metadata (.aux.xml) ¶

Koordinat dövrü SRS elementinin coordinateEpoch atributu kimi kodlanır.

GDAL VRT¶

Koordinat dövrü SRS elementinin coordinateEpoch atributu kimi kodlanır.


GeoJSON-a ixrac edərkən crs göstərin

Məlumat varsa, v.out.ogr, GeoJSON-a ixrac edərkən, ehtimal ki, "crs" atributunu göstərməlidir. ogr2ogr bu xüsusiyyəti dəstəkləyir.

Heç bir xərci olmayan bir GeoJSON faylı:

Yeri gəlmişkən, crs müəyyən edilmədikdə ogrinfo EPSG-yə görə default olaraq görünür: 4326, ən ağlabatan seçimdir, baxmayaraq ki, daha açıq bir şeyə üstünlük verəcəyimi düşünürəm (məsələn, bir xəbərdarlıq):

История изменений (4)

Şərh: 1 5 il əvvəl pmav99 tərəfindən dəyişdirildi

Hmm. ogr2ogr 'un gizli seçimlərinin səhvlərə yol aça biləcəyi görünür:

3857.json-dan CRS sahəsini silsək:

Ogrinfo, EPSG'nin 4326 olduğunu qəbul edir:

reprojecting koordinatlarda heç bir dəyişiklik etmir (təxmin etdiyim gözlənilir):

lakin bu xüsusi (koordinatlar, crs) birləşməsi açıq-aşkar zibildir:

Şərh: 2 5 il əvvəl rouault tərəfindən dəyişdirildi

Geojson spesifikasından sitat gətirmək (http://geojson.org/geojson-spec.html#coordinate-reference-system-objects): "Varsayılan CRS, WGS84 datumunu istifadə edərək və uzunluq və enlem ilə coğrafi koordinat istinad sistemidir. ondalık dərəcə vahidləri. " . Beləliklə, OGR davranışı spesifikasiyanı izləyir.

Və yeni geojson spesifikası (https://www.ietf.org/id/draft-ietf-geojson-00.txt) digər CRS-lərin istifadəsini ləğv edir: "Digər koordinat istinad sistemləri, o cümlədən CRS obyektləri tərəfindən təsvir edilənlər [GJ2008] -də müəyyən edilmiş Tövsiyə edilmir "

Şərh: 3 Nikos Alexandris tərəfindən 5 il əvvəl dəyişdirilib

Düzgün başa düşsəm, bu bilet bağlana bilər?

Şərh: 4 5 il əvvəl dəyişdirilib neteler

Görünür xüsusiyyətlərə hörmət olunur, burada çox şey edə bilmərik. Bağlanır.


Videoya baxın: Sentinel2 images exploration and processing with Python and Rasterio - Tutorial (Oktyabr 2021).