<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr" class="gmail_signature"><div dir="ltr"><div>Hello,</div><div><br></div><div>I'm using GDAL version 1.11.4 and I have a problem how to use gdal to translate a bounding box correctly [Lat(62:63), Long(10:11)] from source CRS EPSG:4326 to target CRS EPSG:32632.</div><div><br></div><div>I've tried so far with gdaltransform to translate coordinates: Long min, Lat min and Long max, Lat max (EPSG:4326) to E min, N min and E max, N Max (EPSG:32632) like this:</div><div><br></div><div><div>gdaltransform -s_srs EPSG:4326 -t_srs EPSG:32632</div><div> </div><div>10 61.9999999999999 0<br></div><div>552375.799656895 6874583.72713382<br></div><div><br></div><div>11 62.9999999999954 0<br></div><div>601293.020582477 6987164.64881585</div><div><br></div></div><div>Then, I use these translated coordinates to create a bounding box in EPSG:32632 E(552375.799656895:601293.020582477), N(6874583.72713382:6987164.64881585) in EPSG:32632 for subsetting a 32632 tiff file and then using gdalwarp the cropped result to EPSG:4326.</div><div><br></div><div>The problem is the results of 2 cropped tiff files (both in EPSG:4326) are significantly different and I could not understand the reason (I think bounding box in EPSG:32632 are wrong but how to find the correct bounding box in EPSG:32632 from bounding box EPSG:4326)?<br></div><div><br></div><div>32632.tiff is a tiff file with grid size: 6820, 4632 and 4326.tiff is a file created by using gdalwarp -t_srs EPSG:4326 32632.tiff 4326.tiff with grid size: 8213, 2654.</div><div><br></div><div><pre> Test 1:
 # Subset on 32632 tiff by 32632 subset which are translated point by point
 from 4326 subset
 gdalwarp -te 552375.799656895 6874583.72713382 601293.020582477
 6987164.64881585 32632.tiff subset_32632.tiff
 Creating output file that is 1631P x 3753L.

 # Warp cropped 32632 tiff file to 4326
 gdalwarp -t_srs EPSG:4326 subset_32632.tiff subset_4326.tiff
 Creating output file that is 3009P x 3070L.

 Test 2:
 # Subset 4326 on 4326 tiff file directly
 gdalwarp -te 10 61.9999999999999 11 62.9999999999954 4326.tiff
 subset_4326_org.tif
 Creating output file that is 2019P x 2019L.
</pre></div><div>Thanks,</div><div><br></div><div><br></div></div></div></div></div></div></div></div>