翻譯與PIL的圖像，並保持

我有我需要翻譯2500px左右，1500px的了一個形象，我用這個代碼的所有內容：翻譯與PIL的圖像，並保持

from PIL import Image 
img = Image.open('decoupe_test4.tif') 
a = 1 
b = 0 
c = 2500 # +left/-right 
d = 0 
e = 1 
f = 1500 # +up/-down 
translate = img.transform(img.size, Image.AFFINE, (a, b, c, d, e, f)) 
translate.save('translated.tif')

圖像的邊框並沒有改變但內容一樣，所以我結束了與此圖像，右邊：

result of translation

有誰知道我怎麼可以得到所有的圖像翻譯正確，遠離黑眼球？

謝謝！

下面是完整的代碼，它可能是更清晰：

import subprocess 
gm = os.path.join('C:\\','OSGeo4W64','bin','gdal_merge.py') 
merge_command = ["python", gm, "-o", "mergedB04.tif", "S2A_tile_20170410_31TFJ_0\B04.jp2", "S2A_tile_20170410_31TGH_0\B04.jp2", "S2A_tile_20170410_31TGJ_0\B04.jp2", "S2A_tile_20170420_31TFH_0\B04.jp2", "S2A_tile_20170420_31TFJ_0\B04.jp2", "S2A_tile_20170420_31TGH_0\B04.jp2"] 
merge_command2 = ["python", gm, "-o", "mergedB08.tif", "S2A_tile_20170410_31TFJ_0\B08.jp2", "S2A_tile_20170410_31TGH_0\B08.jp2", "S2A_tile_20170410_31TGJ_0\B08.jp2", "S2A_tile_20170420_31TFH_0\B08.jp2", "S2A_tile_20170420_31TFJ_0\B08.jp2", "S2A_tile_20170420_31TGH_0\B08.jp2"] 
subprocess.call(merge_command,shell=True) 
subprocess.call(merge_command2,shell=True) 
from PIL import Image 
import rasterio 
from rasterio import Affine 
outfile = r'test4.tif' 
Image.MAX_IMAGE_PIXELS = 1000000000 
im = Image.open('mergedB08.tif') 
half_the_width = im.size[0]/2 
half_the_height = im.size[1]/2 
decoupe = im.crop(
    (
     half_the_width - 2500, 
     half_the_height - 1500, 
     half_the_width + 2500, 
     half_the_height + 1500 
    ) 
) 

decoupe.save('decoupe_test4.tif') 

################### test de translation 
img = Image.open('decoupe_test4.tif') 
a = 1 
b = 0 
c = 2500 #left/right (i.e. 5/-5) 
d = 0 
e = 1 
f = 1500 #up/down (i.e. 5/-5) 
translate = img.transform(img.size, Image.AFFINE, (a, b, c, d, e, f)) 
translate.save('translated.tif') 

# remise dans le bon scr de l'image 
profile = {'driver': 'GTiff', 
      'dtype': 'float32', 
      'nodata': None, 
      'width': decoupe.size[0], 
      'height': decoupe.size[1], 
      'count': 1, 
      'crs': {'init': 'epsg:32631'}, 
      'transform': Affine(10.0, 0.0, 704880.000,0.0, -10.0, 4795110.000) 
      } 
profile.update(driver='GTiff') 
profile.update(dtype=rasterio.float32) 

with rasterio.open('translated.tif') as decoupe_RAS: 
    RAS = decoupe_RAS.read().astype(float) 
with rasterio.open(outfile, 'w', **profile) as dst: 
    dst.write(RAS.astype(rasterio.float32))

來源

2017-06-21 hort_98663

什麼是輸入圖像是什麼樣子？輸出圖像應該是什麼樣子？ –

輸入圖像是較大圖像的切割圖像（通過合併一些衛星圖像sentinel2獲得），但剪切後，此圖像不再被引用。我設法再給它一個乾淨的，但仍然有一個抵消。輸入和輸出需要是相同的圖像，只需通過翻譯進行區分 –

我們需要的是輸入圖像，輸出圖像和預期輸出圖像以及再現問題的最少量代碼。 –

你需要移動它後，裁剪圖像：

img = Image.open('decoupe_test4.tif') 
x_shift = 2500 
y_shift = 1500 
a = 1 
b = 0 
c = x_shift #left/right (i.e. 5/-5) 
d = 0 
e = 1 
f = y_shift #up/down (i.e. 5/-5) 
translate = img.transform(img.size, Image.AFFINE, (a, b, c, d, e, f)) 
# Calculate the size after cropping 
size = (translate.size[0] - x_shift, translate.size[1] - y_shift) 
# Crop to the desired size 
translate = translate.transform(size, Image.EXTENT, (0, 0, size[0], size[1])) 
translate.save('translated.tif')

來源

2017-06-21 19:08:33

翻譯與PIL的圖像，並保持

回答

相關問題