Daha çox

GeoTIFF nümunəsinə əsasən biçin şəklini kəsin və yenidən işləyin


Bu xüsusiyyətlərə sahib GeoTIFF var:

Ölçü 495, 309 Koordinat Sistemi: PROJCS ["adsız", GEOGCS ["adsız ellips", DATUM ["bilinməyən", SPHEROID ["adsız", 6370000,0]], PRIMEM ["Greenwich", 0], UNIT ["dərəcə", 0.0174532925199433]], PROJEKSİYA ["Lambert_Conformal_Conic_1SP"], PARAMETER ["enlem_of_origin", 47.5], PARAMETER ["central_meridian", 4], PARAMETER ["scale_factor", 1], PARET ], PARAMETER ["false_northing", 0], UNIT ["meter", 1, Authority ["EPSG", "9001"]]] PROJ.4 string: '+ proj = lcc + lat_1 = 47.5 + lat_0 = 47.5 + lon_0 = 4 + k_0 = 1 + x_0 = 0 + y_0 = 0 + a = 6370000 + b = 6370000 + ədəd = m + no_defs 'Mənşəyi = (-2963997.870569999800000, -1848004.200800000000000) Piksel Boyutu = (11975.754948886731) Metadata: SIKIŞTIRMA = LZW INTERLEAVE = BAND Köşe Koordinatları: Yuxarı Sol (-2963997.871, -1848004.201) (24d36'23.04 "W, 26d22 '5.88" N) Aşağı Sol (-2963997.871, 1848004.097) (45d39'38.75' W, 55d 35.82 "N) Üst Sağ (2964000.829, -1848004.201) (32d36'23.13" E, 26d22 '5.85 "N) Aşağı Sağ (2964000.829) , 1848004.097) (53d39'38.88 "E, 55d16'35.76" N) Center (1.4791350, -0.0520200) (4d 0 '0.07 "E, 47d30' 0.00" N) Band 1 Blok = 495x4 Type = Float32, ColorInterp = Grey NoData Dəyər = 9.969209968386869e + 036

və Təbii Yer şəklini örtmək istəyirəm.

Buna görə qlobal Natural Earth şəklini bu GeoTIFF nümunəsinə əkib yenidən redaktə etməliyəm, amma bunu belə bir proyeksiya ilə necə edəcəyimi başa düşə bilmirəm və bir növ təlimat verə biləcəyimi düşünürdüm.ogr2ogrGeoTIFF-i biçmək və shapefile yenidən işləmək üçün istinad kimi istifadə etmək?

Həm də başqa yanaşmalarla maraqlanıramogr2ogrbu xüsusiyyəti yoxdur.


GeoTIFF etiketlərində saxlanılan Proj4 sətri düzgündürsə, bu asan bir vəziyyət olmalıdır. Qırpmadan belə görünür

ogr2ogr -f "ESRI Shapefile" -s_srs epsg: 4326 -t_srs "+ proj = lcc + lat_1 = 47.5 + lat_0 = 47.5 + lon_0 = 4 + k_0 = 1 + x_0 = 0 + y_0 = 0 + a = 6370000 + b = 6370000 + ədəd = m + no_defs "projected_out.shp natural_earth_in.shp

Hedef ölçülərə qədər əkin etmək üçün gdalinfo hesabatında künc koordinatlarını görə bilərsiniz. İn ogr2ogr parametrindən istifadə edin -clipdst xmin ymin xmax ymax. Əlavə məlumatı http://www.gdal.org/ogr2ogr.html.

Redaktə edin

Miqyasları və hədəf proyeksiyasını hədəf alan dünya ölkələrini təkzib edən və klip çəkən tam əmrdir

ogr2ogr -f "ESRI Shapefile" -s_srs epsg: 4326 -t_srs "+ proj = lcc + lat_1 = 47.5 + lat_0 = 47.5 + lon_0 = 4 + k_0 = 1 + x_0 = 0 + y_0 = 0 + a = 6370000 + b = 6370000 + ədəd = m + no_defs "projed_out.shp ne_110m_admin_0_countries.shp -clipdst -2963997 18480 04 2964000 -1848004 --fiqırlama OGR_ENABLE_PARTIAL_REPROJECTION TRUE -skipfailures

Bu halda məcburi olan GDAL konfiqurasiya seçiminin OGR_ENABLE_PARTIAL_REPROJECTION istifadəsinə diqqət yetirin, çünki bütün mənbəni yalnız məhdud bir ərazi üçün etibarlı olan proyeksiyaya yenidən proqnozlaşdırmaq mümkün deyil. Həm də mənbə məlumatlarının topoloji səhvləri olduğu üçün və ya təkrar proyeksiya prosesində öz-özünə kəsişmələr yaradıldığı üçün -skip çatışmazlıqlarına ehtiyac var idi. Nəticə hələ də mükəmməl deyildi - Rusiya itkin düşmüşdür və Antarktida bütün Avropanı əhatə edir, bu da bəzi araşdırmalara ehtiyac olacağını göstərir. Ancaq başqa bir suala dəyər. Aşağıdakı şəkildəki açıq fon, proqnozlaşdırılan Antarktidadır.


Videoya baxın: Convert GeoTiff to XYZ using QGIS (Oktyabr 2021).